The USGS 3D Seismic Velocity Model for Northern California
NASA Astrophysics Data System (ADS)
Brocher, T. M.; Aagaard, B.; Simpson, R. W.; Jachens, R. C.
2006-12-01
We present a new regional 3D seismic velocity model for Northern California for use in strong motion simulations of the 1906 San Francisco and other earthquakes. The model includes compressional-wave velocity (Vp), shear-wave velocity (Vs), density, and intrinsic attenuation (Qp, Qs). These properties were assigned for each rock type in a 3D geologic model derived from surface outcrops, boreholes, gravity and magnetic data, and seismic reflection, refraction, and tomography studies. A detailed description of the model, USGS Bay Area Velocity Model 05.1.0, is available online [http://www.sf06simulation.org/geology/velocitymodel]. For ground motion simulations Vs and Qs are more important parameters than Vp and Qp because the strongest ground motions are generated chiefly by shear and surface wave arrivals. Because Vp data are more common than Vs data, however, we first developed Vp versus depth relations for each rock type and then converted these to Vs versus depth relations. For the most important rock types in Northern California we compiled measurements of Vp versus depth using borehole logs, laboratory measurements on hand samples, seismic refraction profiles, and tomography models. These rock types include Salinian and Sierran granitic rocks, metagraywackes and greenstones of the Franciscan Complex, Tertiary and Mesozoic sedimentary and volcanic rocks, and Quaternary and Holocene deposits (Brocher, USGS OFR 05-1317, 2005). Vp versus depth curves were converted to Vs versus depth curves using new empirical nonlinear relations between Vs and Vp (Brocher, BSSA, 2005). These relations, showing that Poisson's ratio is a nonlinear function of Vp, were similarly based on compilations of diverse Vs and Vp measurements on a large suite of rock types, mainly from California and the Pacific Northwest. The model is distributed in a discretized form with routines to query the model using C++, C, and Fortran 77 programming languages. The geologic model was discretized at
Towards an Anisotropic Whole Mantle 3D Elastic Velocity Model
NASA Astrophysics Data System (ADS)
Panning, M. P.; Romanowicz, B.; Gung, Y.
2001-12-01
Many studies have documented the existence of anisotropy in the earth's upper mantle, concentrated in the top 200 km. This evidence comes from the study of surface waves as well as shear wave splitting. There is also evidence for shear wave splitting in D", at least in well sampled regions. There are some hints of anisotropy at the base of the transition zone. Tomographic models of the upper mantle have been developed with simplifying assumptions about the nature of the anisotropy, in order to minimize the number of free parameters in the inversions. Some assume transverse isotropy (e.g Ekström and Dziewonski, 1997), others include additional degrees of freedom with some realistic constraints on mineralogy (e.g. Montagner and Tanimoto, 1991). Our goal is to investigate anisotropy in the whole mantle, using the framework of waveform inversion, and the nonlinear asymptotic mode coupling theory (NACT), previously developed and applied to the construction of whole-mantle SH velocity models (Li and Romanowicz, 1996; Mégnin and Romanowicz, 2000). For this we require a 3 component dataset, and we have extended our automatic transverse (T) component wavepicking procedures to the vertical (Z) and longitudinal (L) component - a non-trivial task given the large number of phases present in the coupled P-SV system. A useful initial assumption, for which the theory has been readily adapted, is that of transverse isotropy. As a first step towards this, we have been investigating inversions using T component and Z,L component data separately. In particular, this allows us to explore the sampling that can be achieved with Z,L component data alone in the deepest part of the mantle. Indeed, D" is in general much better sampled in SH than in SV, owing to the availability of SHdiff at large distances, while SVdiff decays more rapidly due to mantle-core coupling. We present the results of our resolution experiments and discuss the differences between the 3D SV model obtained in well
Uncertainty assessment of 3D instantaneous velocity model from stack velocities
NASA Astrophysics Data System (ADS)
Emanuele Maesano, Francesco; D'Ambrogi, Chiara
2015-04-01
3D modelling is a powerful tool that is experiencing increasing applications in data analysis and dissemination. At the same time the need of quantitative uncertainty evaluation is strongly requested in many aspects of the geological sciences and by the stakeholders. In many cases the starting point for 3D model building is the interpretation of seismic profiles that provide indirect information about the geology of the subsurface in the domain of time. The most problematic step in the 3D modelling construction is the conversion of the horizons and faults interpreted in time domain to the depth domain. In this step the dominant variable that could lead to significantly different results is the velocity. The knowledge of the subsurface velocities is related mainly to punctual data (sonic logs) that are often sparsely distributed in the areas covered by the seismic interpretation. The extrapolation of velocity information to wide extended horizons is thus a critical step to obtain a 3D model in depth that can be used for predictive purpose. In the EU-funded GeoMol Project, the availability of a dense network of seismic lines (confidentially provided by ENI S.p.A.) in the Central Po Plain, is paired with the presence of 136 well logs, but few of them have sonic logs and in some portion of the area the wells are very widely spaced. The depth conversion of the 3D model in time domain has been performed testing different strategies for the use and the interpolation of velocity data. The final model has been obtained using a 4 layer cake 3D instantaneous velocity model that considers both the initial velocity (v0) in every reference horizon and the gradient of velocity variation with depth (k). Using this method it is possible to consider the geological constraint given by the geometries of the horizons and the geo-statistical approach to the interpolation of velocities and gradient. Here we present an experiment based on the use of set of pseudo-wells obtained from the
Congruence of 3-D Whole Mantle Models of Shear Velocity
NASA Astrophysics Data System (ADS)
Dziewonski, A. M.; Lekic, V.; Romanowicz, B. A.
2012-12-01
The range of shear velocity anomalies in published whole mantle models is considerable. This impedes drawing conclusions of importance for geodynamic modeling and for interpretation of mineral physics results. However, if one considers only the models that were built using data that are sensitive to mantle structure at all depths, these models show robust features in their power spectra as a function of depth. On this basis we propose that there are five depth intervals with distinct spectral characteristics. 1. Heterosphere (Moho - 300 km) is characterized by strong power spectrum relatively flat up to degree 6. With lateral shear wavespeed variations as large as 15%, this zone accounts for more than 50% of the entire heterogeneity in the mantle. Differences among models for different tectonic regions decrease rapidly below 300 km depth. 2. Upper mantle buffer zone (300- 500 km) has a flat spectrum and the overall power of heterogeneity drops by an order of magnitude compared to the region above. There may be still weak difference between continents and oceans, but the oceanic regions lose their age dependence. The spectral characteristics do not change across the 410 km discontinuity. 3. Transition zone (500 - 650 km) The degree 2 anomaly becomes dominant. There are long wavelength anomalies in regions of the fastest plate subduction during the last 15-20 Ma, suggesting slab ponding above the 650 km discontinuity. Several slower-than-average anomalies of unknown origin are present in this depth range. 4. Lower mantle buffer zone (650 - 2300 km) has a weak, flat spectrum without long wavelength velocity anomalies that could be interpreted as unfragmented subducted slabs. However, there are three relatively narrow and short high velocity anomalies under Peru, Tonga and Indonesia that may indicate limited slab penetration. 5 Abyssal layer (2300 - CMB) Strong spectrum dominated by degrees 2 and 3. The amplitude is the largest at the CMB and decreases rapidly up to
UCVM: Open Source Software for Understanding and Delivering 3D Velocity Models
NASA Astrophysics Data System (ADS)
Gill, D.; Small, P.; Maechling, P. J.; Jordan, T. H.; Shaw, J. H.; Plesch, A.; Chen, P.; Lee, E. J.; Taborda, R.; Olsen, K. B.; Callaghan, S.
2014-12-01
Physics-based ground motion simulations can calculate the propagation of earthquake waves through 3D velocity models of the Earth. The Southern California Earthquake Center (SCEC) has developed the Unified Community Velocity Model (UCVM) framework to help researchers build structured or unstructured velocity meshes from 3D velocity models for use in wave propagation simulations. The UCVM software framework makes it easy to extract P and S wave propagation speeds and other material properties from 3D velocity models by providing a common interface through which researchers can query earth models for a given location and depth. Currently, the platform supports multiple California models, including SCEC CVM-S4 and CVM-H 11.9.1, and has been designed to support models from any region on earth. UCVM is currently being use to generate velocity meshes for many SCEC wave propagation codes, including AWP-ODC-SGT and Hercules. In this presentation, we describe improvements to the UCVM software. The current version, UCVM 14.3.0, released in March of 2014, supports the newest Southern California velocity model, CVM-S4.26, which was derived from 26 full-3D tomographic iterations using CVM-S4 as the starting model (Lee et al., this meeting), and the Broadband 1D velocity model used in the CyberShake 14.2 study. We have ported UCVM to multiple Linux distributions and OS X. Also included in this release is the ability to add small-scale stochastic heterogeneities to extract Cartesian meshes for use in high-frequency ground motion simulations. This tool was built using the C language open-source FFT library, FFTW. The stochastic parameters (Hurst exponent, correlation length, and the horizontal/vertical aspect ratio) can be customized by the user. UCVM v14.3.0 also provides visualization scripts for constructing cross-sections, horizontal slices, basin depths, and Vs30 maps. The interface allows researchers to visually review velocity models . Also, UCVM v14.3.0 can extract
UCVM: An Open Source Software Package for Querying and Visualizing 3D Velocity Models
NASA Astrophysics Data System (ADS)
Gill, D.; Small, P.; Maechling, P. J.; Jordan, T. H.; Shaw, J. H.; Plesch, A.; Chen, P.; Lee, E. J.; Taborda, R.; Olsen, K. B.; Callaghan, S.
2015-12-01
Three-dimensional (3D) seismic velocity models provide foundational data for ground motion simulations that calculate the propagation of earthquake waves through the Earth. The Southern California Earthquake Center (SCEC) has developed the Unified Community Velocity Model (UCVM) package for both Linux and OS X. This unique framework provides a cohesive way for querying and visualizing 3D models. UCVM v14.3.0, supports many Southern California velocity models including CVM-S4, CVM-H 11.9.1, and CVM-S4.26. The last model was derived from 26 full-3D tomographic iterations on CVM-S4. Recently, UCVM has been used to deliver a prototype of a new 3D model of central California (CCA) also based on full-3D tomographic inversions. UCVM was used to provide initial plots of this model and will be used to deliver CCA to users when the model is publicly released. Visualizing models is also possible with UCVM. Integrated within the platform are plotting utilities that can generate 2D cross-sections, horizontal slices, and basin depth maps. UCVM can also export models in NetCDF format for easy import into IDV and ParaView. UCVM has also been prototyped to export models that are compatible with IRIS' new Earth Model Collaboration (EMC) visualization utility. This capability allows for user-specified horizontal slices and cross-sections to be plotted in the same 3D Earth space. UCVM was designed to help a wide variety of researchers. It is currently being use to generate velocity meshes for many SCEC wave propagation codes, including AWP-ODC-SGT and Hercules. It is also used to provide the initial input to SCEC's CyberShake platform. For those interested in specific data points, the software framework makes it easy to extract P and S wave propagation speeds and other material properties from 3D velocity models by providing a common interface through which researchers can query earth models for a given location and depth. Also included in the last release was the ability to add small
3-D crustal velocity model for Lithuania and its application to local event studies
NASA Astrophysics Data System (ADS)
Budraitis, M.; Kozlovskaya, E.; Janutyte, I.; Motuza, G.
2009-12-01
PASSEQ 2006-2008 project (PASsive Seismic Experiment in TESZ) aimed at studying the lithosphere-asthenosphere system around the TransEuropean Suture Zone (TESZ)- the transition between old Proterozoic platform of north and east Europe and younger Phanerozoic platform in central and western Europe. The experiment was a seismic array research aiming to retrieve the structure of the crust and Earth's mantle down to the mantle transition zone, including mapping of upper mantle seismic velocity variations and discontinuities (Moho, lithosphere-asthenosphere boundary, mantle transition zone) using all available techniques. During the experiment 26 seismic stations (including four broadband stations) were installed in Lithuania and operated since June, 2006 till January, 2008. One of the main reasons of PASSEQ deployment in Lithuania is identification and characterisation of the local seismic activity. During the data acquisition period a number of local seismic events was identified and preliminary event location was made using LocSat and VELEST algorithms and 1-D velocity models. These standard procedures is not enough precise for Lithuania, however, because the thickness of the crust varies significantly in the region (from 45 to 55 km). Another problem was low quality of S-wave arrivals due to thick (up to 2 km) sediments in most part of Lithuania. In order to improve event location, we compiled a 3-D seismic velocity model of the crust down to a depth of 60 km. The model, consisting of four major layers (sediments, upper crust, middle crust, lower crust and uppermost mantle) was interpolated from 2-D velocity models along previous wide-angle reflection and refraction profiles into a regular grid. The quality of the approximation was analysed using comparison of travel times of P-waves recorded by controlled source experiments and calculated travel times through the 3-D velocity model. The model was converted into a density model using a special procedure, in which
Probabilistic earthquake location and 3-D velocity models in routine earthquake location
NASA Astrophysics Data System (ADS)
Lomax, A.; Husen, S.
2003-12-01
Earthquake monitoring agencies, such as local networks or CTBTO, are faced with the dilemma of providing routine earthquake locations in near real-time with high precision and meaningful uncertainty information. Traditionally, routine earthquake locations are obtained from linearized inversion using layered seismic velocity models. This approach is fast and simple. However, uncertainties derived from a linear approximation to a set of non-linear equations can be imprecise, unreliable, or even misleading. In addition, 1-D velocity models are a poor approximation to real Earth structure in tectonically complex regions. In this paper, we discuss the routine location of earthquakes in near real-time with high precision using non-linear, probabilistic location methods and 3-D velocity models. The combination of non-linear, global search algorithms with probabilistic earthquake location provides a fast and reliable tool for earthquake location that can be used with any kind of velocity model. The probabilistic solution to the earthquake location includes a complete description of location uncertainties, which may be irregular and multimodal. We present applications of this approach to determine seismicity in Switzerland and in Yellowstone National Park, WY. Comparing our earthquake locations to earthquake locations obtained using linearized inversion and 1-D velocity models clearly demonstrates the advantages of probabilistic earthquake location and 3-D velocity models. For example, the more complete and reliable uncertainty information of non-linear, probabilistic earthquake location greatly facilitates the identification of poorly constrained hypocenters. Such events are often not identified in linearized earthquake location, since the location uncertainties are determined with a simplified, localized and approximate Gaussian statistic.
Development of a Regional Velocity Model Using 3D Broadband Waveform Sensitivity
NASA Astrophysics Data System (ADS)
Panning, M. P.; Romanowicz, B. A.; Kim, A.
2005-12-01
We are developing a new approach which relies on a cascade of increasingly accurate theoretical approximations for computation of the seismic wavefield to develop a model of regional seismic velocity structure for eastern Eurasia using full seismic waveforms. The selected area is particularly suitable for the purpose of this experiment, as it is highly heterogeneous, presenting a challenge for standard modeling techniques, but it is well surrounded by earthquake sources and a significant number of high quality broadband digital stations exist, for which data are readily accessible through IRIS (Incorporated Research Institutions for Seismology) and the FDSN (Federation of Digital Seismic Networks). The initial model is derived from a large database of teleseismic long period waveforms (surface waves and overtone wavepackets) using well-developed theoretical approximations, the Path Average Approximation (PAVA) and Nonlinear Asymptotic Coupling Theory (NACT). These approaches assume waveforms are only sensitive to the 1D (PAVA) and 2D (NACT) structure in the vertical plane between source and receiver, which is adequate for the development of a smooth initial 3D velocity model. We refine this model using a more accurate theoretical approach. We utilize an implementation of a 3D Born approximation, which takes into account the contribution to the waveform from single scattering throughout the model, giving full 3D waveform sensitivity kernels. We perform verification tests of this approach for synthetic models, and show that it can accurately represent the wavefield as predicted by numerical approaches in several situations where approximations such as PAVA and NACT are insufficient. The Born 3D waveform sensitivity kernels are used to perform a higher resolution inversion of regional waveforms for a smaller subregion between longitudes 90 and 150 degrees E, and latitudes 15 and 40 degrees N. To further increase the accuracy of this model, we intend to utilize a very
3-D crustal velocity model for Lithuania and its application to local event studies
NASA Astrophysics Data System (ADS)
Kozlovskaya, Elena; Budraitis, Mantas; Janutyte, Ilma; Motuza, Gediminas; Lazauskiene, Jurga; Passeq-Working Group
2010-05-01
PASSEQ 2006-2008 project (PASsive Seismic Experiment in TESZ) aimed at studying the lithosphere-asthenosphere system around the TransEuropean Suture Zone (TESZ)- the transition between old Proterozoic platform of north and east Europe and younger Phanerozoic platform in central and western Europe. The experiment was a seismic array research aiming to retrieve the structure of the crust and Earth's mantle down to the mantle transition zone, including mapping of upper mantle seismic velocity variations and discontinuities (Moho, lithosphere-asthenosphere boundary, mantle transition zone) using all available techniques. During the experiment 26 seismic stations (including four broadband stations) were installed in Lithuania. One of the main targets of PASSEQ deployment in Lithuania was identification and characterization of the local seismic activity. The PASSEQ stations in Lithuania were in operation since June, 2006 till January, 2008. During this period a number of local seismic events was recorded and preliminary event location was made using the LocSat algorithm and 1-D velocity model. This standard procedure is not enough precise for Lithuania, however, because the thickness of the crust varies significantly in the region (from 45 to 55 km). In order to improve event location, we separated the events into several groups and located each group separately using a VELEST algorithms and own 1-D velocity model for each group. We also compiled a 3-D seismic velocity of the crust down to a depth of 60 km. The model, consisting of four major layers (sediments, upper crust, middle crust, lower crust and uppermost mantle) was interpolated from 2-D velocity models along previous wide-angle reflection and refraction profiles into a regular grid. The quality of the approximation was analysed using comparison of travel times of P-waves recorded along previous controlled source profiles and synthetic travel times calculated using the 3-D velocity model. The model was converted
Stephenson, William J.
2007-01-01
INTRODUCTION In support of earthquake hazards and ground motion studies in the Pacific Northwest, three-dimensional P- and S-wave velocity (3D Vp and Vs) and density (3D rho) models incorporating the Cascadia subduction zone have been developed for the region encompassed from about 40.2?N to 50?N latitude, and from about -122?W to -129?W longitude. The model volume includes elevations from 0 km to 60 km (elevation is opposite of depth in model coordinates). Stephenson and Frankel (2003) presented preliminary ground motion simulations valid up to 0.1 Hz using an earlier version of these models. The version of the model volume described here includes more structural and geophysical detail, particularly in the Puget Lowland as required for scenario earthquake simulations in the development of the Seattle Urban Hazards Maps (Frankel and others, 2007). Olsen and others (in press) used the model volume discussed here to perform a Cascadia simulation up to 0.5 Hz using a Sumatra-Andaman Islands rupture history. As research from the EarthScope Program (http://www.earthscope.org) is published, a wealth of important detail can be added to these model volumes, particularly to depths of the upper-mantle. However, at the time of development for this model version, no EarthScope-specific results were incorporated. This report is intended to be a reference for colleagues and associates who have used or are planning to use this preliminary model in their research. To this end, it is intended that these models will be considered a beginning template for a community velocity model of the Cascadia region as more data and results become available.
Azimuthally Anisotropic 3D Velocity Continuation
Burnett, William; Fomel, Sergey
2011-01-01
We extend time-domain velocity continuation to the zero-offset 3D azimuthally anisotropic case. Velocity continuation describes how a seismic image changes given a change in migration velocity. This description turns out to be of a wave propagation process, in which images change along a velocity axis. In the anisotropic case, the velocity model is multiparameter. Therefore, anisotropic image propagation is multidimensional. We use a three-parameter slowness model, which is related to azimuthal variations in velocity, as well as their principal directions. This information is useful for fracture and reservoir characterization from seismic data. We provide synthetic diffraction imaging examples to illustratemore » the concept and potential applications of azimuthal velocity continuation and to analyze the impulse response of the 3D velocity continuation operator.« less
Validating 3D Seismic Velocity Models Using the Spectral Element Method
NASA Astrophysics Data System (ADS)
Maceira, M.; Rowe, C. A.; Allen, R. M.; Obrebski, M. J.
2010-12-01
As seismic instrumentation, data storage and dissemination and computational power improve, seismic velocity models attempt to resolve smaller structures and cover larger areas. However, it is unclear how accurate these velocity models are and, while the best models available are used for event determination, it is difficult to put uncertainties on seismic event parameters. Model validation is typically done using resolution tests that assume the imaging theory used is accurate and thus only considers the impact of the data coverage on resolution. We present the results of a more rigorous approach to model validation via full three-dimensional waveform propagation using Spectral Element Methods (SEM). This approach makes no assumptions about the theory used to generate the models but require substantial computational resources. We first validate 3D tomographic models for the Western USA generated using both ray-theoretical and finite-frequency methods. The Dynamic North America (DNA) Models of P- and S- velocity structure (DNA09-P and DNA09-S) use teleseismic body-wave traveltime residuals recorded at over 800 seismic stations provided by the Earthscope USArray and regional seismic networks. We performed systematic computations of synthetics for the dataset used to generate the DNA models. Direct comparison of these synthetic seismograms to the actual observations allows us to accurately assess and validate the models. Implementation of the method for a densely instrumented region such as that covered by the DNA model provides a useful testbed for the validation methods that we will subsequently apply to other, more challenging study areas.
Using Adjoint Methods to Improve 3-D Velocity Models of Southern California
NASA Astrophysics Data System (ADS)
Liu, Q.; Tape, C.; Maggi, A.; Tromp, J.
2006-12-01
We use adjoint methods popular in climate and ocean dynamics to calculate Fréchet derivatives for tomographic inversions in southern California. The Fréchet derivative of an objective function χ(m), where m denotes the Earth model, may be written in the generic form δχ=int Km(x) δln m(x) d3x, where δln m=δ m/m denotes the relative model perturbation. For illustrative purposes, we construct the 3-D finite-frequency banana-doughnut kernel Km, corresponding to the misfit of a single traveltime measurement, by simultaneously computing the 'adjoint' wave field s† forward in time and reconstructing the regular wave field s backward in time. The adjoint wave field is produced by using the time-reversed velocity at the receiver as a fictitious source, while the regular wave field is reconstructed on the fly by propagating the last frame of the wave field saved by a previous forward simulation backward in time. The approach is based upon the spectral-element method, and only two simulations are needed to produce density, shear-wave, and compressional-wave sensitivity kernels. This method is applied to the SCEC southern California velocity model. Various density, shear-wave, and compressional-wave sensitivity kernels are presented for different phases in the seismograms. We also generate 'event' kernels for Pnl, S and surface waves, which are the Fréchet kernels of misfit functions that measure the P, S or surface wave traveltime residuals at all the receivers simultaneously for one particular event. Effectively, an event kernel is a sum of weighted Fréchet kernels, with weights determined by the associated traveltime anomalies. By the nature of the 3-D simulation, every event kernel is also computed based upon just two simulations, i.e., its construction costs the same amount of computation time as an individual banana-doughnut kernel. One can think of the sum of the event kernels for all available earthquakes, called the 'misfit' kernel, as a graphical
Coronal Outflow Velocities in a 3D Coronal Model Determined from UVCS Doppler Dimming Observations
NASA Astrophysics Data System (ADS)
Strachan, L.; Panasyuk, A. V.; Dobrzycka, D.; Gibson, S.; Biesecker, D. A.; Ko, Y.-K.; Galvin, A. B.; Romoli, M.; Kohn, J. L.
1998-04-01
We constrain coronal outflow velocity solutions, resolved along the line-of-sight, by using Doppler dimming models of H I Lyman alpha and O VI 1032/1037 Angstrom emissivities obtained with data from the Ultraviolet Coronagraph Spectrometer (UVCS) on SOHO. The local emissivities, from heliocentric heights of 1.5 to 3.0 radii, were determined from 3-D reconstructions of line-of-sight intensities obtained during the Whole Sun Month Campaign (10 Aug. -- 8 Sep. 1996). The models use electron densities derived from polarized brightness measurements made with the visible light coronagraphs on UVCS and LASCO, supplemented with data from Mark III at NCAR/MLSO. Electron temperature profiles are derived from `freezing-in' temperatures obtained from an analysis of charge state data from SWICS/Ulysses. The work concentrates on O5+ outflow velocities which are determined from an analysis of the the O VI line ratios. This analysis is less sensitive to the uncertainties in the electron density and independent of the ionization balance and elemental abundance than the analyses which use individual spectral lines. This work is supported in part by NASA under grant NAG-3192 to the Smithsonian Astrophysical Observatory, by the Italian Space Agency and by Swiss funding agencies.
NASA Astrophysics Data System (ADS)
Begnaud, M. L.; Ballard, S.; Young, C. J.; Hipp, J. R.; Encarnacao, A.; Phillips, W. S.; Chael, E. P.; Rowe, C. A.
2012-12-01
We are developing a global 3D P wave velocity model of the Earth's crust and mantle using seismic tomography to assess improvement to seismic event locations obtained using high quality 3D Earth models in lieu of 1D and 2/2.5D models. We present the most recent version of SALSA3D (SAndia LoS Alamos 3D) version 1.9, and demonstrate its ability to reduce mislocations for a large set of realizations derived from a carefully chosen set of globally-distributed ground truth (GT) events. Our model is derived from the latest version of the GT catalog of P/Pn travel-time picks assembled by Los Alamos National Laboratory. For this current version, we employ more robust data quality control measures than previously used, as well as additional global GT data sources. To prevent over-weighting due to ray path redundancy and to reduce the computational burden, we cluster rays into representative rays. The model is represented using the triangular tessellation system described by Ballard et al. (2009), which incorporates variable resolution in both the geographic and radial dimensions. For our starting model, we use a simplified layer crustal model derived from the NNSA Unified model in Eurasia and Crust 2.0 model everywhere else, overlying a uniform ak135 mantle. Sufficient damping is used to reduce velocity adjustments so that ray path changes between iterations are small. We obtain proper model smoothness by using progressive grid refinement, refining the grid only in areas where the data warrant such a refinement. In previous versions, we based this refinement on velocity changes from previous model iterations. For the current version, we utilize the diagonal of the model resolution matrix to control where grid refinement occurs, resulting in more consistent and continuous areas of refinement than before. In addition to the changes in grid refinement, we also employ a more robust convergence criterion between successive grid refinements, allowing a better fit to first broader
NASA Astrophysics Data System (ADS)
Begnaud, M. L.; Ballard, S.; Young, C. J.; Hipp, J. R.; Chang, M.; Encarnacao, A.; Rowe, C. A.; Phillips, W. S.; Steck, L.
2011-12-01
To test the hypothesis that high quality 3D Earth models will produce seismic event locations that are more accurate and more precise than currently used 1D and 2/2.5D models, we are developing a global 3D P wave velocity model of the Earth's crust and mantle using seismic tomography. In this paper, we present the most recent version of our model, SALSA3D (SAndia LoS Alamos 3D) version 1.7, and demonstrate its ability to reduce mislocations for a large set of realizations derived from a carefully chosen set of globally-distributed ground truth (GT) events, compared to existing models and/or systems. Our model is derived from the latest version of the GT catalog of P and Pn travel time picks assembled by Los Alamos National Laboratory. To prevent over-weighting due to ray path redundancy and to reduce the computational burden, we cluster rays to produce representative rays. Reduction in the total number of ray paths is ~50%. The model is represented using the triangular tessellation system described by Ballard et al. (2009), which incorporates variable resolution in both the geographic and radial dimensions. For our starting model, we use a simplified layer crustal model derived from the NNSA Unified model in Eurasia and Crust 2.0 model elsewhere, over a uniform ak135 mantle. Sufficient damping is used to reduce velocity adjustments so that ray path changes between iterations are small. We obtain proper model smoothness by using progressive grid refinement, refining the grid only in areas where the data warrant it. In previous versions of SALSA3D, we based this refinement on velocity changes from previous model iterations. For version 1.7, we utilize the diagonal of the model resolution matrix to control where grid refinement occurs, resulting in more consistent and continuous areas of refinement than before. Our approach produces a smooth, multi-resolution model with node density appropriate to both ray coverage and the velocity gradients required by the data. We
A New Regional 3-D Velocity Model of the India-Pakistan Region for Improved Event Location Accuracy
NASA Astrophysics Data System (ADS)
Reiter, D.; Vincent, C.; Johnson, M.
2001-05-01
A 3-D velocity model for the crust and upper mantle (WINPAK3D) has been developed to improve regional event location in the India-Pakistan region. Results of extensive testing demonstrate that the model improves location accuracy for this region, specifically for the case of small regionally recorded events, for which teleseismic data may not be available. The model was developed by integrating the results of more than sixty previous studies related to crustal velocity structure in the region. We evaluated the validity of the 3-D model using the following methods: (1) cross validation analysis for a variety of events, (2) comparison of model determined hypocenters with known event location, and (3) comparison of model-derived and empirically-derived source-specific station corrections (SSSC) generated for the International Monitoring System (IMS) auxiliary seismic station located at Nilore. The 3-D model provides significant improvement in regional location compared to both global and regional 1-D models in this area of complex structural variability. For example, the epicenter mislocation for an event with a well known location was only 6.4 km using the 3-D model, compared with a mislocation of 13.0 km using an average regional 1-D model and 15.1 km for the IASPEI91 model. We will present these and other results to demonstrate that 3-D velocity models are essential to improving event location accuracy in regions with complicated crustal geology and structures. Such 3-D models will be a prerequisite for achieving improved location accuracies for regions of high monitoring interest.
NASA Astrophysics Data System (ADS)
Polkowski, Marcin; Grad, Marek
2015-04-01
The 3D P-wave seismic velocity model was obtained by combining data from multiple studies during past 50 years. Data sources included refraction seismology, reflection seismology, geological boreholes, vertical seismic profiling, magnetotellurics and gravimetry. Use of many data sources allowed creation of detailed 3D P-wave velocity model that reaches to depth of 60 km and includes 6-layers of sediments and 3-layers of the crust. Purpose of this study is to analyze how 3D model influences local (accuracy of location and source time estimation for local events), regional (identification of wide-angle seismic phases) and global (teleseismic tomography) seismic travel times. Additionally we compare results of forward seismic wave propagation with signals observed on short period and broadband stations. National Science Centre Poland provided financial support for this work by NCN grant DEC-2011/02/A/ST10/00284.
SALSA3D - A Global 3D P-Velocity Model of the Earth's Crust and Mantle for Improved Event Location
NASA Astrophysics Data System (ADS)
Ballard, S.; Begnaud, M. L.; Young, C. J.; Hipp, J. R.; Chang, M.; Encarnacao, A. V.; Rowe, C. A.; Phillips, W. S.; Steck, L.
2010-12-01
To test the hypothesis that high quality 3D Earth models will produce seismic event locations which are more accurate and more precise, we are developing a global 3D P wave velocity model of the Earth’s crust and mantle using seismic tomography. In this paper, we present the most recent version of our model, SALSA3D version 1.5, and demonstrate its ability to reduce mislocations for a large set of realizations derived from a carefully chosen set of globally-distributed ground truth events. Our model is derived from the latest version of the Ground Truth (GT) catalog of P and Pn travel time picks assembled by Los Alamos National Laboratory. To prevent over-weighting due to ray path redundancy and to reduce the computational burden, we cluster rays to produce representative rays. Reduction in the total number of ray paths is ~50%. The model is represented using the triangular tessellation system described by Ballard et al. (2009), which incorporates variable resolution in both the geographic and radial dimensions.. For our starting model, we use a simplified two layer crustal model derived from the Crust 2.0 model over a uniform AK135 mantle. Sufficient damping is used to reduce velocity adjustments so that ray path changes between iterations are small. We obtain proper model smoothness by using progressive grid refinement, refining the grid only around areas with significant velocity changes from the starting model. At each grid refinement level except the last one we limit the number of iterations to prevent convergence thereby preserving aspects of broad features resolved at coarser resolutions. Our approach produces a smooth, multi-resolution model with node density appropriate to both ray coverage and the velocity gradients required by the data. This scheme is computationally expensive, so we use a distributed computing framework based on the Java Parallel Processing Framework, providing us with ~400 processors. Resolution of our model is assessed using a
NASA Astrophysics Data System (ADS)
Ramachandran, K.
2011-12-01
Three dimensional velocity models constructed through seismic tomography are seldom digitally processed further for imaging structural features. A study conducted to evaluate the potential for imaging subsurface discontinuities in horizontal and vertical direction from three dimensional velocity models using image processing/computer vision techniques has provided significant results. Three-dimensional velocity models constructed through tomographic inversion of active source and/or earthquake traveltime data are generally built from an initial 1-D velocity model that varies only with depth. Regularized tomographic inversion algorithms impose constraints on the roughness of the model that help to stabilize the inversion process. Final velocity models obtained from regularized tomographic inversions have smooth three-dimensional structures that are required by the data. Final velocity models are usually analyzed and interpreted either as a perturbation velocity model or as an absolute velocity model. Compared to perturbation velocity model, absolute velocity model has an advantage of providing constraints on lithology. Both velocity models lack the ability to provide sharp constraints on subsurface faults. However, results from the analysis of the 3-D velocity model from northern Cascadia using Roberts, Prewitt, Sobel, and Canny operators show that subsurface faults that are not clearly interpretable from velocity model plots can be identified through this approach. This analysis resulted in inferring the locations of Tacoma Fault, Seattle Fault, Southern Whidbey Island Fault, and Darrington Devils Mountain fault much clearly. The Coast Range Boundary Fault, previously hypothesized on the basis of sedimentological and tectonic observations is inferred clearly from processed images. Many of the fault locations so imaged correlate with earthquake hypocenters indicating their seismogenic nature.
NASA Astrophysics Data System (ADS)
Floriane, Provost; Jean-Philippe, Malet; Cécile, Doubre; Julien, Gance; Alessia, Maggi; Agnès, Helmstetter
2015-04-01
Characterizing the micro-seismic activity of landslides is an important parameter for a better understanding of the physical processes controlling landslide behaviour. However, the location of the seismic sources on landslides is a challenging task mostly because of (a) the recording system geometry, (b) the lack of clear P-wave arrivals and clear wave differentiation, (c) the heterogeneous velocities of the ground. The objective of this work is therefore to test whether the integration of a 3D velocity model in probabilistic seismic source location codes improves the quality of the determination especially in depth. We studied the clay-rich landslide of Super-Sauze (French Alps). Most of the seismic events (rockfalls, slidequakes, tremors...) are generated in the upper part of the landslide near the main scarp. The seismic recording system is composed of two antennas with four vertical seismometers each located on the east and west sides of the seismically active part of the landslide. A refraction seismic campaign was conducted in August 2014 and a 3D P-wave model has been estimated using the Quasi-Newton tomography inversion algorithm. The shots of the seismic campaign are used as calibration shots to test the performance of the different location methods and to further update the 3D velocity model. Natural seismic events are detected with a semi-automatic technique using a frequency threshold. The first arrivals are picked using a kurtosis-based method and compared to the manual picking. Several location methods were finally tested. We compared a non-linear probabilistic method coupled with the 3D P-wave model and a beam-forming method inverted for an apparent velocity. We found that the Quasi-Newton tomography inversion algorithm provides results coherent with the original underlaying topography. The velocity ranges from 500 m.s-1 at the surface to 3000 m.s-1 in the bedrock. For the majority of the calibration shots, the use of a 3D velocity model
Development of a State-Wide 3-D Seismic Tomography Velocity Model for California
NASA Astrophysics Data System (ADS)
Thurber, C. H.; Lin, G.; Zhang, H.; Hauksson, E.; Shearer, P.; Waldhauser, F.; Hardebeck, J.; Brocher, T.
2007-12-01
We report on progress towards the development of a state-wide tomographic model of the P-wave velocity for the crust and uppermost mantle of California. The dataset combines first arrival times from earthquakes and quarry blasts recorded on regional network stations and travel times of first arrivals from explosions and airguns recorded on profile receivers and network stations. The principal active-source datasets are Geysers-San Pablo Bay, Imperial Valley, Livermore, W. Mojave, Gilroy-Coyote Lake, Shasta region, Great Valley, Morro Bay, Mono Craters-Long Valley, PACE, S. Sierras, LARSE 1 and 2, Loma Prieta, BASIX, San Francisco Peninsula and Parkfield. Our beta-version model is coarse (uniform 30 km horizontal and variable vertical gridding) but is able to image the principal features in previous separate regional models for northern and southern California, such as the high-velocity subducting Gorda Plate, upper to middle crustal velocity highs beneath the Sierra Nevada and much of the Coast Ranges, the deep low-velocity basins of the Great Valley, Ventura, and Los Angeles, and a high- velocity body in the lower crust underlying the Great Valley. The new state-wide model has improved areal coverage compared to the previous models, and extends to greater depth due to the data at large epicentral distances. We plan a series of steps to improve the model. We are enlarging and calibrating the active-source dataset as we obtain additional picks from investigators and perform quality control analyses on the existing and new picks. We will also be adding data from more quarry blasts, mainly in northern California, following an identification and calibration procedure similar to Lin et al. (2006). Composite event construction (Lin et al., in press) will be carried out for northern California for use in conventional tomography. A major contribution of the state-wide model is the identification of earthquakes yielding arrival times at both the Northern California Seismic
Hartzell, S.; Harmsen, S.; Williams, R.A.; Carver, D.; Frankel, A.; Choy, G.; Liu, P.-C.; Jachens, R.C.; Brocher, T.M.; Wentworth, C.M.
2006-01-01
A 3D seismic velocity and attenuation model is developed for Santa Clara Valley, California, and its surrounding uplands to predict ground motions from scenario earthquakes. The model is developed using a variety of geologic and geophysical data. Our starting point is a 3D geologic model developed primarily from geologic mapping and gravity and magnetic surveys. An initial velocity model is constructed by using seismic velocities from boreholes, reflection/refraction lines, and spatial autocorrelation microtremor surveys. This model is further refined and the seismic attenuation is estimated through waveform modeling of weak motions from small local events and strong-ground motion from the 1989 Loma Prieta earthquake. Waveforms are calculated to an upper frequency of 1 Hz using a parallelized finite-difference code that utilizes two regions with a factor of 3 difference in grid spacing to reduce memory requirements. Cenozoic basins trap and strongly amplify ground motions. This effect is particularly strong in the Evergreen Basin on the northeastern side of the Santa Clara Valley, where the steeply dipping Silver Creek fault forms the southwestern boundary of the basin. In comparison, the Cupertino Basin on the southwestern side of the valley has a more moderate response, which is attributed to a greater age and velocity of the Cenozoic fill. Surface waves play a major role in the ground motion of sedimentary basins, and they are seen to strongly develop along the western margins of the Santa Clara Valley for our simulation of the Loma Prieta earthquake.
NASA Astrophysics Data System (ADS)
Persaud, P.; Stock, J. M.; Fuis, G. S.; Hole, J. A.; Goldman, M.; Scheirer, D. S.
2014-12-01
We have analyzed explosive shot data from the 2011 Salton Seismic Imaging Project (SSIP) across a 2-D seismic array and 5 profiles in the Coachella Valley to produce a 3-D P-wave velocity model that will be used in calculations of strong ground shaking. Accurate maps of seismicity and active faults rely both on detailed geological field mapping and a suitable velocity model to accurately locate earthquakes. Adjoint tomography of an older version of the SCEC 3-D velocity model shows that crustal heterogeneities strongly influence seismic wave propagation from moderate earthquakes (Tape et al., 2010). These authors improve the crustal model and subsequently simulate the details of ground motion at periods of 2 s and longer for hundreds of ray paths. Even with improvements such as the above, the current SCEC velocity model for the Salton Trough does not provide a match of the timing or waveforms of the horizontal S-wave motions, which Wei et al. (2013) interpret as caused by inaccuracies in the shallow velocity structure. They effectively demonstrate that the inclusion of shallow basin structure improves the fit in both travel times and waveforms. Our velocity model benefits from the inclusion of known location and times of a subset of 126 shots detonated over a 3-week period during the SSIP. This results in an improved velocity model particularly in the shallow crust. In addition, one of the main challenges in developing 3-D velocity models is an uneven stations-source distribution. To better overcome this challenge, we also include the first arrival times of the SSIP shots at the more widely spaced Southern California Seismic Network (SCSN) in our inversion, since the layout of the SSIP is complementary to the SCSN. References: Tape, C., et al., 2010, Seismic tomography of the Southern California crust based on spectral-element and adjoint methods: Geophysical Journal International, v. 180, no. 1, p. 433-462. Wei, S., et al., 2013, Complementary slip distributions
NASA Astrophysics Data System (ADS)
Gaite, B.; Villaseñor, A.; Iglesias, A.; Herraiz, M.; Jiménez-Munt, I.
2015-02-01
We use group velocities from earthquake tomography together with group and phase velocities from ambient noise tomography (ANT) of Rayleigh waves to invert for the 3-D shear-wave velocity structure (5-70 km) of the Caribbean (CAR) and southern North American (NAM) plates. The lithospheric model proposed offers a complete image of the crust and uppermost-mantle with imprints of the tectonic evolution. One of the most striking features inferred is the main role of the Ouachita-Marathon-Sonora orogeny front on the crustal seismic structure of the NAM plate. A new imaged feature is the low crustal velocities along the USA-Mexico border. The model also shows a break of the east-west mantle velocity dichotomy of the NAM and CAR plates beneath the Isthmus of the Tehuantepec and the Yucatan Block. High upper-mantle velocities along the Mesoamerican Subduction Zone coincide with inactive volcanic areas while the lowest velocities correspond to active volcanic arcs and thin lithospheric mantle regions.
NASA Astrophysics Data System (ADS)
Coviello, Velio; Manconi, Andrea; Occhiena, Cristina; Arattano, Massimo; Scavia, Claudio
2013-04-01
Rock-falls are one of the most common and hazardous phenomena occurring in mountainous areas. The formation of cracks in rocks is often accompanied by a sudden release of energy, which propagates in form of elastic waves and can be detected by a suitable transducer array. Therefore, geophones are among the most effective monitoring devices to investigate eventual precursors of rock-fall phenomena. However, the identification of an efficient procedure to forecast rock-fall occurrence in space and time is still an open challenge. In this study, we aim at developing an efficient procedure to locate microseismic sources relevant to cracking mechanisms, and thus gather indications on eventual precursors of rock-fall phenomena. Common seismic location tools usually implement homogeneous or multilayered velocity models but, in case of high slope gradients and heavily fractured rock masses, these simplifications may lead to errors on the correct estimation of the source location. Thus, we analyzed how the consideration of 3D material properties on the propagation medium may influence the location. In the framework of the Alcotra 2007-2013 Project MASSA (Medium And Small Size rock-fall hazard Assessment), a monitoring system composed by 8 triaxial geophones was installed in 2010 at the J.A. Carrel hut (3829 m a.s.l., Matterhorn, NW Italian Alps) and during the first year of operation the network recorded more than 600 natural events that exceeded a fixed threshold [1]. Despite the harsh environmental conditions of the study area, eighteen points distributed as uniformly as possible in space were selected for hammering. The artificial source dataset of known coordinates was used to constrain a 3D heterogeneous velocity model through a Simultaneous Iterative Reconstructive Technique. In order to mitigate the intrinsic uncertainties of the inversion procedure, bootstrapping was performed to extend the dataset and a statistical analysis was issued to improve the model
NASA Astrophysics Data System (ADS)
Trovato, Claudio; Aochi, Hideo; De Martin, Florent
2014-05-01
Understanding the source mechanism of long-period (LP) seismic signals on volcanoes is an important key point in volcanology and for the hazard forecasting. In the last decades, moment tensor inversions have led to various descriptions of the kinematic source mechanism. These inversions suppose a relatively simple structure of the medium. However, the seismic wave propagation in a realistic 3-D volcano model should be taken into account for understanding the complicated physical processes of magma and gas behaviors at depth. We are studying Etna volcano, Italy, to understand the volcanic processes during different stages of activity. We adopt a spectral element method (SEM), a code EFISPEC3D (De Martin, BSSA, 2011), which shows a good accuracy and numerical stability in the simulations of seismic wave propagation. First we construct the geometrical model. We use a digital elevation model (DEM) to generate finite element meshes with a spacing of 50 m on the ground surface. We aim to calculate the ground motions until 3 Hz for the shallowest layer with Vs = ~500 m/s. The minimal size of the hexahedral elements is required to be around 100 m, with a total number of elements n = ~2 10 ^ 6 for the whole model. We compare different velocity structure configurations. We start with a homogeneous medium and add complexities taking in account the shallow low velocity structure. We also introduce a velocity gradient towards depth. Simulations performed in the homogeneous medium turn in approximately 20 hours for calculations parallelized on 16 CPUs. Complex velocity models should take approximately the same time of computation. We then try to simulate the ground motion from the LP sources (0.1-1.5 Hz) obtained by the inversion for the Etna volcano in 2008 (De Barros, GRL, 2009 and De Barros, JGR, 2011). Some vertical and horizontal structures can be added to reproduce injected dikes or sills respectively.
NASA Astrophysics Data System (ADS)
Kohler, M. D.; Bowden, D. C.; Tsai, V. C.; Weeraratne, D. S.
2014-12-01
The Pacific-North America plate boundary in Southern California extends far west of the coastline, and a 12-month ocean bottom seismometer (OBS) array spanned the western side of the plate boundary in order to image seismic velocities in the lithosphere. Velocities are modeled through stacked cross correlations of ambient noise data. The offshore data come primarily from the OBS array that collected 12 months of continuous data during 2010-2011, combined with Southern California Seismic Network (SCSN) station data. The cross correlations were stacked for noise correlation functions and examined using standard time- and frequency-domain methods to determine phase velocity and group velocity dispersion curves. Signals between the vertical-component OBS and co-located horizontal-component OBS observations associated with tilt noise, and pressure gauge observations associated with infragravity waves, were examined to further improve signals. The non-elastic noise was estimated by calculating the transfer functions between the vertical-to-horizontal and vertical-to-pressure components, and subtracting the coherent signal between the two from the vertical-component time series. We find that these effects are small in our dataset. We are simultaneously inverting all measureable dispersion curves to solve for 3D crustal velocity structure. Shear-wave velocities comprise the direct solution, and Vp/Vs ratios are constrained as much as the data allow. Calculations on data from 780 OBS-OBS, SCSN-SCSN, and OBS-SCSN pairs filtered around multiple narrow bands between 5 and 50 s show clear propagating waves traveling at group velocities between 1.2 and 3.5 km/s. The longer-term outcome of this work will comprise a 3D crustal and uppermost mantle velocity model with areal coverage not attainable before the deployment of the ocean bottom seismometers. The results define the transition in three dimensions from continental lithospheric structure in the near-shore region to oceanic
NASA Astrophysics Data System (ADS)
Kohler, M. D.; Bowden, D. C.; Tsai, V. C.; Weeraratne, D. S.
2015-12-01
The Pacific-North America plate boundary in Southern California extends far west of the coastline, and a 12-month ocean bottom seismometer (OBS) array spanned the western side of the plate boundary to image lithospheric seismic velocities. Velocities are modeled through stacked cross correlations of ambient noise data. Twelve months of continuous data were used from 22 OBS stations and ~30 coastal and island Southern California Seismic Network stations. Particular attention has been paid to improving signal-to-noise ratios in the noise correlations with OBS stations by removing the effects of instrument tilt and infragravity waves. Different applications of preprocessing techniques allow us to distinguish the fundamental and first higher order Rayleigh modes, especially in deep water OBS pairs where the water layer dominates crustal sensitivity of the fundamental mode. Standard time domain and frequency domain methods are used to examine surface wave dispersion curves for group and phase velocities between 5 and 50 second periods, and these are inverted for 3D velocity structure. The results define the transition in three dimensions from continental lithospheric structure in the near-shore region to oceanic structure west of the continental borderland. While the most prominent features of the model relate to thinning of the crust west of the Patton Escarpment, other notable anomalies are present north-to-south throughout the continental borderland and along the coast from the Los Angeles Basin to the Peninsular Ranges. The velocity model will help describe the region's tectonic history, as well as provide new constraints for determination of earthquake relocations and rupture styles.
Nugroho, Hendro; Widiyantoro, Sri; Nugraha, Andri Dian
2013-09-09
Determination of earthquake hypocenter in Indonesia conducted by the Meteorological, Climatological, and Geophysical Agency (MCGA) has still used a 1-D seismic velocity model. In this research, we have applied a Fast Grid Search (FGM) method and a 3-D velocity model resulting from tomographic imaging to relocate earthquakes in the Sumatran region. The data were taken from the MCGA data catalog from 2009 to 2011 comprising of subduction zone and on land fault earthquakes with magnitude greater than 4 Mw. Our preliminary results show some significant changes in the depths of the relocated earthquakes which are in general deeper than the depths of hypocenters from the MCGA data catalog. The residual times resulting from the relocation process are smaller than those prior to the relocation. Encouraged by these results, we will continue to conduct hypocenter relocation for all events from the MCGA data catalog periodically in order to produce a new data catalog with good quality. We hope that the new data catalog will be useful for further studies.
NASA Astrophysics Data System (ADS)
Font, Yvonne; Segovia, Monica; Vaca, Sandro; Theunissen, Thomas
2013-04-01
To improve earthquake location, we create a 3-D a priori P-wave velocity model (3-DVM) that approximates the large velocity variations of the Ecuadorian subduction system. The 3-DVM is constructed from the integration of geophysical and geological data that depend on the structural geometry and velocity properties of the crust and the upper mantle. In addition, specific station selection is carried out to compensate for the high station density on the Andean Chain. 3-D synthetic experiments are then designed to evaluate the network capacity to recover the event position using only P arrivals and the MAXI technique. Three synthetic earthquake location experiments are proposed: (1) noise-free and (2) noisy arrivals used in the 3-DVM, and (3) noise-free arrivals used in a 1-DVM. Synthetic results indicate that, under the best conditions (exact arrival data set and 3-DVM), the spatiotemporal configuration of the Ecuadorian network can accurately locate 70 per cent of events in the frontal part of the subduction zone (average azimuthal gap is 289° ± 44°). Noisy P arrivals (up to ± 0.3 s) can accurately located 50 per cent of earthquakes. Processing earthquake location within a 1-DVM almost never allows accurate hypocentre position for offshore earthquakes (15 per cent), which highlights the role of using a 3-DVM in subduction zone. For the application to real data, the seismicity distribution from the 3-D-MAXI catalogue is also compared to the determinations obtained in a 1-D-layered VM. In addition to good-quality location uncertainties, the clustering and the depth distribution confirm the 3-D-MAXI catalogue reliability. The pattern of the seismicity distribution (a 13 yr record during the inter-seismic period of the seismic cycle) is compared to the pattern of rupture zone and asperity of the Mw = 7.9 1942 and the Mw = 7.7 1958 events (the Mw = 8.8 1906 asperity patch is not defined). We observe that the nucleation of 1942, 1958 and 1906 events coincides with
NASA Astrophysics Data System (ADS)
Lin, F. C.; Schmandt, B.
2015-12-01
Imaging the crust and lithosphere structure beneath North America is one of the primary targets for the NSF-funded EarthScope project. In this study, we apply the recently developed ambient noise and surface wave tomography methods to construct a detailed 3D crustal model across the entire contiguous US using USArray data between January 2007 and May 2015. By using both Rayleigh wave phase velocity and ellipticity measurements between 8 and 100 sec period, the shear velocity structure can be well resolved within the five crustal layers we modeled: three upper crust, one middle crust, and one lower crust. Clear correlations are observed between the resolved velocity anomalies and known geological features at all depths. In the uppermost crust, slow Vs anomalies are observed within major sedimentary environments such as the Williston Basin, Denver Basin, and Mississippi embayment, and fast Vs anomalies are observed in environments with deeply exhumed bedrock outcrops at the surface including the Laurentian Highlands, Ouachita-Ozark Interior Highlands, and Appalachian Highlands. In the deeper upper crust, slow anomalies are observed in deep sedimentary basins such as the Green River Basin, Appalachian Basin, Southern Oklahoma Aulacogen, and areas surrounding the Gulf of Mexico. Fast anomalies, on the other hand, are observed in the Colorado Plateau, within the Great Plains between the Front Ranges and Midcontinental Rift, and east of the Appalachian Mountains. At this depth, the Midcontinental Rift and Grenville Front clearly correlate well with various velocity structure boundaries. In the middle crust, slow anomalies are mostly observed in the tectonically active areas in the western US, but relatively slow anomalies are also observed southeast of the Precambrian Rift Margins. At this depth, fast anomalies are observed beneath various deep sedimentary basins such as the Southern Oklahoma Aulacogen, Appalachian Basin, and Central Valley. In the lower crust, a clear
A Global 3D P-Velocity Model of the Earth's Crust and Mantle for Improved Event Location
NASA Astrophysics Data System (ADS)
Ballard, S.; Young, C. J.; Hipp, J. R.; Chang, M.; Lewis, J.; Begnaud, M. L.; Rowe, C. A.
2009-12-01
further refinement takes place around adjusted nodes to form a new model, and the process is repeated until no more improvement can be obtained. We thus produce a smooth, multi-resolution model with node density appropriate to both ray coverage and the velocity gradients required by the data. This scheme is computationally expensive, so we use a Java-based distributed computing framework developed by Sandia National Laboratories (SNL), providing us with 300+ processors having an efficiency of better than 90% for the calculations. We evaluate our model both in terms of travel time residual variance reduction and in location improvement for GT events. For the latter, we use a new multi-threaded version of the SNL-developed LocOO code modified to use 3D velocity models.
NASA Astrophysics Data System (ADS)
Ballard, S.; Begnaud, M. L.; Hipp, J. R.; Chael, E. P.; Encarnacao, A.; Maceira, M.; Yang, X.; Young, C. J.; Phillips, W.
2013-12-01
SALSA3D is a global 3D P wave velocity model of the Earth's crust and mantle developed specifically to provide seismic event locations that are more accurate and more precise than are locations from 1D and 2.5D models. In this paper, we present the most recent version of our model, for the first time jointly derived from multiple types of data: body wave travel times, surface wave group velocities, and gravity. The latter two are added to provide information in areas with poor body wave coverage, and are down-weighted in areas where body wave coverage is good. To constrain the inversions, we invoked empirical relations among the density, S velocity, and P velocity. We demonstrate the ability of the new SALSA3D model to reduce mislocations and generate statistically robust uncertainty estimates for a large set of realizations derived from a carefully chosen set of globally-distributed ground truth events. We obtain path-dependent travel time prediction uncertainties for our model by computing the full 3D model covariance matrix of our tomographic system and integrating the model slowness variance and covariance along paths of interest. This approach yields very low travel time prediction uncertainties for well-sampled paths through the Earth and higher uncertainties for paths that are poorly represented in the data set used to develop the model. While the calculation of path-dependent prediction uncertainties with this approach is computationally expensive, uncertainties can be pre-computed for a network of stations and stored in 3D lookup tables that can be quickly and efficiently interrogated using GeoTess software.
NASA Astrophysics Data System (ADS)
Ayadi, A.; Dorbath, C.; Ousadou, F.; Maouche, S.; Chikh, M.; Bounif, M. A.; Meghraoui, M.
2008-09-01
We analyze the aftershocks sequence of the Zemmouri thrust faulting earthquake (21 May 2003, Mw 6.8) located east of Algiers in the Tell Atlas. The seismic sequence located during ˜2 months following the mainshock is made of more than 1500 earthquakes and extends NE-SW along a ˜60-km fault rupture zone crossing the coastline. The earthquake relocation was performed using handpicked P and S phases located with the tomoDD in a detailed 3D velocity structure of the epicentral area. Contrasts between velocity patches seem to correlate with contacts between granitic-volcanic basement rocks and the sedimentary formation of the eastern Mitidja basin. The aftershock sequence exhibits at least three seismic clouds and a well-defined SE-dipping main fault geometry that reflects the complex rupture. The distribution of seismic events presents a clear contrast between a dense SW zone and a NE zone with scattered aftershocks. We observe that the mainshock locates between the SW and NE seismic zones; it also lies at the NNS-SSE contact that separates a basement block to the east and sedimentary formations to the west. The aftershock distribution also suggests fault bifurcation at the SW end of the fault rupture, with a 20-km-long ˜N 100° trending seismic cluster, with a vertical fault geometry parallel to the coastline juxtaposed. Another aftershock cloud may correspond to 75° SE dipping fault. The fault geometry and related SW branches may illustrate the interference between pre-existing fault structures and the SW rupture propagation. The rupture zone, related kinematics, and velocity contrasts obtained from the aftershocks distribution are in agreement with the coastal uplift and reflect the characteristics of an active zone controlled by convergent movements at a plate boundary.
3D finite element simulations of high velocity projectile impact
NASA Astrophysics Data System (ADS)
Ožbolt, Joško; İrhan, Barış; Ruta, Daniela
2015-09-01
An explicit three-dimensional (3D) finite element (FE) code is developed for the simulation of high velocity impact and fragmentation events. The rate sensitive microplane material model, which accounts for large deformations and rate effects, is used as a constitutive law. In the code large deformation frictional contact is treated by forward incremental Lagrange multiplier method. To handle highly distorted and damaged elements the approach based on the element deletion is employed. The code is then used in 3D FE simulations of high velocity projectile impact. The results of the numerical simulations are evaluated and compared with experimental results. It is shown that it realistically predicts failure mode and exit velocities for different geometries of plain concrete slab. Moreover, the importance of some relevant parameters, such as contact friction, rate sensitivity, bulk viscosity and deletion criteria are addressed.
NASA Astrophysics Data System (ADS)
Lozano, Lucía; Cantavella, Juan Vicente; Barco, Jaime; Carranza, Marta; Burforn, Elisa
2016-04-01
The Atlantic margin of the SW Iberian Peninsula and northern Morocco has been subject of study during the last 30 years. Many seismic reflection and refraction profiles have been carried out offshore, providing detailed information about the crustal structure of the main seafloor tectonic domains in the region, from the South Portuguese Zone and the Gulf of Cadiz to the Abyssal Plains and the Josephine Seamount. The interest to obtain a detailed and realistic velocity model for this area, integrating the available data from these studies, is clear, mainly to improve real-time earthquake hypocentral location and for tsunami and earthquake early warning. Since currently real-time seismic location tools allow the implementation of 3D velocity models, we aim to generate a full 3D crustal model. For this purpose we have reviewed more than 50 profiles obtained in different seismic surveys, from 1980 to 2008. Data from the most relevant and reliable 2D seismic velocity published profiles were retrieved. We first generated a Moho depth map of the studied area (latitude 32°N - 41°N and longitude 15°W - 5°W) by extracting Moho depths along each digitized profile with a 10 km spacing, and then interpolating this dataset using ordinary kriging method and generating the contour isodepth map. Then, a 3D crustal velocity model has been obtained. Selected vertical sections at different distances along each profile were considered to retrieve P-wave velocity values at each interface in order to reproduce the geometry and the velocity gradient within each layer. A double linear interpolation, both in distance and depth, with sampling rates of 10 km and 1 km respectively, was carried out to generate a (latitude, longitude, depth, velocity) matrix. This database of all the profiles was interpolated to obtain the P-wave velocity distribution map every kilometer of depth. The new 3D velocity model has been integrated in NonLinLoc location program to relocate several representative
NASA Astrophysics Data System (ADS)
Béthoux, Nicole; Theunissen, Thomas; Beslier, Marie-Odile; Font, Yvonne; Thouvenot, François; Dessa, Jean-Xavier; Simon, Soazig; Courrioux, Gabriel; Guillen, Antonio
2016-02-01
The region between the inner zones of the Alps and Corsica juxtaposes an overthickened crust to an oceanic domain, which makes difficult to ascertain the focal depth of seismic events using routine location codes and average 1D velocity models. The aim of this article is to show that, even with a rather lose monitoring network, accurate routine locations can be achieved by using realistic 3D modelling and advanced location techniques. Previous earthquake tomography studies cover the whole region with spatial resolutions of several tens of kilometres on land, but they fail to resolve the marine domain due to the absence of station coverage and sparse seismicity. To overcome these limitations, we first construct a 3D a-priori P and S velocity model integrating known geophysical and geological information. Significant progress has been achieved in the 3D numerical modelling of complex geological structures by the development of dedicated softwares (e.g. 3D GeoModeller), capable at once of elaborating a 3D structural model from geological and geophysical constraints and, possibly, of refining it by inversion processes (Calcagno et al., 2008). Then, we build an arrival-time catalogue of 1500 events recorded from 2000 to 2011. Hypocentres are then located in this model using a numerical code based on the maximum intersection method (Font et al., 2004), updated by Theunissen et al. (2012), as well as another 3D location technique, the NonLinLoc software (Lomax and Curtis, 2001). The reduction of arrival-time residuals and uncertainties (dh, dz) with respect to classical 1D locations demonstrates the improved accuracy allowed by our approach and confirms the coherence of the 3D geological model built and used in this study. Our results are also compared with previous works that benefitted from the installation of dense temporary networks surrounding the studied epicentre area. The resulting 3D location catalogue allows us to improve the regional seismic hazard assessment
NASA Astrophysics Data System (ADS)
Środa, Piotr; Dec, Monika
2016-04-01
The area of Trans-European Suture Zone in SE Poland represents a contact of major tectonic units of different consolidation age - from the Precambrian East European Craton, through Palaeozoic West European Platform to Cenozoic Carpathian orogen. The region was built by several phases of crustal accretion, which resulted in a complex collage of tectonic blocks. In 2000, this region was studied by several seismic wide-angle profiles of CELEBRATION 2000 experiment, providing a dense coverage of seismic data in SE Poland and allowing for detailed investigations of the crustal structure and properties in this area. Beneath the marginal part of the EEC, the 2-D modelling of in-line data form several CELEBRATION profiles revealed a prominent high P-wave velocity anomaly in the upper crust, with Vp of 6.7-7.1 km/s, starting at 10-16 km depth (e.g., Środa et al., 2006). Anomalously high velocities are observed in the area located approximately beneath Lublin trough, to the NE of Teisseyre-Tornquist Zone. Based on 3-D tomography of first arrivals of in- and off-line CELEBRATION 2000 recordings (Malinowski et al., 2008), elevated velocities are also reported in the same area and seem to continue to the SW, off the craton margin. Gravimetric modelling also revealed anomalously high density in the same region at similar depths. High seismic velocities and densities are interpreted as indicative for a pronounced mafic intrusion, possibly related to extensional processes at the EEC margin. Previous 3-D models of the high-velocity intrusion were based on first arrivals (crustal refractions) only. In this study, also off-line reflections (not modelled up to now) are used, in order to enlarge the data set and to better constrain the geometry and properties of the velocity anomaly. A code for 3-D joint tomographic inversion of refracted and reflected arrivals, with model parametrization allowing for velocity discontinuities was used (Rawlinson, 2007). With this approach, besides the
NASA Astrophysics Data System (ADS)
Gaite, B.; Villaseñor, A.; Iglesias, A.; Herraiz, M.; Jiménez-Munt, I.
2014-10-01
We use group velocities from earthquake tomography together with group and phase velocities from ambient noise tomography (ANT) of Rayleigh-waves to invert for the 3-D shear-wave velocity structure (5-70 km) of the Caribbean (CAR) and southern North American (NAM) plates. The lithospheric model proposed offers a complete image of the crust and uppermost-mantle with imprints of the tectonic evolution. One of the most striking features inferred is the main role of the Ouachita-Marathon-Sonora orogeny front on the crustal seismic structure of NAM plate. A new imaged feature is the low crustal velocities along USA-Mexico border. The model also shows a break of the E-W mantle velocity dichotomy of the NAM and CAR plates beneath the Isthmus of Tehuantepec and Yucatan Block. High upper-mantle velocities along the Mesoamerican Subduction Zone coincide with inactive volcanic areas while the lowest velocities correspond to active volcanic arcs and thin lithospheric mantle regions.
A high resolution 3D velocity model beneath the Tokyo Metropolitan area by MeSO-net
NASA Astrophysics Data System (ADS)
Nakagawa, S.; Sakai, S.; Honda, R.; Kimura, H.; Hirata, N.
2015-12-01
Beneath the Tokyo metropolitan area, the Philippine Sea Plate (PSP) subducts and causes devastating mega-thrust earthquakes, such as the 1703 Genroku earthquake (M8.0) and the 1923 Kanto earthquake (M7.9). An M7 or greater (M7+) earthquake in this area at present has high potential to produce devastating serious loss of life and property with even greater global economic repercussions. The Central Disaster Management Council of Japan estimates that an M7+ earthquake will cause 23,000 fatalities and 95 trillion yen (about 1 trillion US$) economic loss. We have launched the Special Project for Reducing Vulnerability for Urban Mega Earthquake Disasters in collaboration with scientists, engineers, and social-scientists in nationwide institutions since 2012. We analyze data from the dense seismic array called Metropolitan Seismic Observation network (MeSO-net), which has 296 seismic stations with spacing of 5 km (Sakai and Hirata, 2009; Kasahara et al., 2009). We applied the double-difference tomography method (Zhang and Thurber, 2003) and estimated the velocity structure and the upper boundary of PSP (Nakagawa et al., 2010). The 2011 Tohoku-oki earthquake (M9.0) has activated seismicity also in Kanto region, providing better coverage of ray paths for tomographic analysis. We obtain much higher resolution velocity models from whole dataset observed by MeSO-net between 2008 and 2015. A detailed image of tomograms shows that PSP contacts Pacific plate at a depth of 50 km beneath northern Tokyo bay. A variation of velocity along the oceanic crust suggests dehydration reaction to produce seismicity in a slab, which may related to the M7+ earthquake. Acknowledgement: This study was supported by the Special Project for Reducing Vulnerability for Urban Mega Earthquake Disasters of MEXT, Japan and the Earthquake Research Institute cooperative research program.
BEAMS3D Neutral Beam Injection Model
Lazerson, Samuel
2014-04-14
With the advent of applied 3D fi elds in Tokamaks and modern high performance stellarators, a need has arisen to address non-axisymmetric effects on neutral beam heating and fueling. We report on the development of a fully 3D neutral beam injection (NBI) model, BEAMS3D, which addresses this need by coupling 3D equilibria to a guiding center code capable of modeling neutral and charged particle trajectories across the separatrix and into the plasma core. Ionization, neutralization, charge-exchange, viscous velocity reduction, and pitch angle scattering are modeled with the ADAS atomic physics database [1]. Benchmark calculations are presented to validate the collisionless particle orbits, neutral beam injection model, frictional drag, and pitch angle scattering effects. A calculation of neutral beam heating in the NCSX device is performed, highlighting the capability of the code to handle 3D magnetic fields.
García-Ramos, F. Javier; Malón, Hugo; Aguirre, A. Javier; Boné, Antonio; Puyuelo, Javier; Vidal, Mariano
2015-01-01
A computational fluid dynamics (CFD) model of the air flow generated by an air-assisted sprayer equipped with two axial fans was developed and validated by practical experiments in the laboratory. The CFD model was developed by considering the total air flow supplied by the sprayer fan to be the main parameter, rather than the outlet air velocity. The model was developed for three air flows corresponding to three fan blade settings and assuming that the sprayer is stationary. Actual measurements of the air velocity near the sprayer were taken using 3D sonic anemometers. The workspace sprayer was divided into three sections, and the air velocity was measured in each section on both sides of the machine at a horizontal distance of 1.5, 2.5, and 3.5 m from the machine, and at heights of 1, 2, 3, and 4 m above the ground The coefficient of determination (R2) between the simulated and measured values was 0.859, which demonstrates a good correlation between the simulated and measured data. Considering the overall data, the air velocity values produced by the CFD model were not significantly different from the measured values. PMID:25621611
NASA Astrophysics Data System (ADS)
Simmons, N. A.; Myers, S. C.; Johannesson, G.; Matzel, E.
2011-12-01
LLNL-G3D is a global-scale model of P-wave velocity designed to accurately predict seismic travel times at regional and teleseismic distances simultaneously. The underlying goal of the model is to provide enhanced seismic event location capabilities. Previous versions of LLNL-G3D (versions 1 and 2) provide substantial improvements in event location accuracy via 3-D ray tracing. The latest models are based on ~2.7 million P and Pn arrivals that are re-processed using our global multi-event locator known as BayesLoc. Bayesloc is a formulation of the joint probability distribution across multiple-event location parameters, including hypocenters, travel time corrections, pick precision, and phase labels. Modeling the whole multiple-event system results in accurate locations and an internally consistent data set that is ideal for tomography. Our recently developed inversion approach (called Progressive Multi-level Tessellation Inversion or PMTI) captures regional trends and fine details where data warrant. Using PMTI, we model multiple heterogeneity scale lengths without defining parameter grids with variable densities based on some ad hoc criteria. LLNL-G3Dv3 (version 3) is produced with data generated with the BayesLoc procedure, recently modified to account for localized travel time trends via a multiple event clustering technique. We demonstrate the significance of BayesLoc processing, the impact on the resulting tomographic images, and the application of LLNL-G3D to seismic event location. This work performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344. LLNL-ABS-491805.
NASA Astrophysics Data System (ADS)
Wilde-Piorko, Monika; Polkowski, Marcin; Grad, Marek
2015-04-01
Geological and seismic structure under area of Poland is well studied by over one hundred thousand boreholes, over thirty deep seismic refraction and wide angle reflection profiles and by vertical seismic profiling, magnetic, gravity, magnetotelluric and thermal methods. Compilation of these studies allowed to create a high-resolution 3D P-wave velocity model down to 60 km depth in the area of Poland (Polkowski et al. 2014). Model also provides details about the geometry of main layers of sediments (Tertiary and Quaternary, Cretaceous, Jurassic, Triassic, Permian, old Paleozoic), consolidated/crystalline crust (upper, middle and lower) and uppermost mantle. This model gives an unique opportunity for calculation synthetic receiver function and compering it with observed receiver function calculated for permanent and temporary seismic stations. Modified ray-tracing method (Langston, 1977) can be used directly to calculate the response of the structure with dipping interfaces to the incoming plane wave with fixed slowness and back-azimuth. So, 3D P-wave velocity model has been interpolated to 2.5D P-wave velocity model beneath each seismic station and back-azimuthal sections of components of receiver function have been calculated. Vp/Vs ratio is assumed to be 1.8, 1.67, 1.73, 1.77 and 1.8 in the sediments, upper/middle/lower consolidated/crystalline crust and uppermost mantle, respectively. Densities were calculated with combined formulas of Berteussen (1977) and Gardner et al. (1974). Additionally, to test a visibility of the lithosphere-asthenosphere boundary phases at receiver function sections models have been extended to 250 km depth based on P4-mantle model (Wilde-Piórko et al., 2010). National Science Centre Poland provided financial support for this work by NCN grant DEC-2011/02/A/ST10/00284 and by NCN grant UMO-2011/01/B/ST10/06653.
NASA Astrophysics Data System (ADS)
Simpson, D. G.; Lipatov, A. S.; Sittler, E. C.; Hartle, R. E.; Cooper, J. F.
2013-12-01
Wave-particle interactions play a very important role in the plasma dynamics near Titan: mass loading, excitation of the low-frequency waves and the formation of the particle velocity distribution function, e.g. ring/shell-like distributions, etc. The kinetic approach is important for estimation of the collision processes e.g. a charge exchange. The particle velocity distribution function also plays a key role for understanding the observed particle fluxes. In this report we discuss the ion velocity distribution function dynamics from 3D hybrid modeling. The modeling is based on recent analysis of the Cassini Plasma Spectrometer (CAPS) ion measurements during the TA flyby. In our model the background ions, all pickup ions, and ionospheric ions are considered as particles, whereas the electrons are described as a fluid. Inhomogeneous photoionization, electron-impact ionization and charge exchange are included in our model. The temperatures of the background electrons and pickup electrons were also included into the generalized Ohm's law. We also take into account the collisions between the ions and neutrals. We use Chamberlain profiles for the exosphere's components and include a simple ionosphere model with M=28 ions that were generated inside the ionosphere. The moon is considered as a weakly conducting body. Our modeling shows that interaction between background plasma and pickup ions H+, H2+, CH4+ and N2+ has a more complicated structure than was observed in the T9 flyby and modeling due to the large gyroradius of the background O+ ions [1,2,3,4]. Special attention will be paid to comparing the simulated pickup ion velocity distribution with CAPS TA observations. We also compare our kinetic modeling with other hybrid and MHD modeling of Titan's environment. References [1] Sittler, E.C., et al., Energy Deposition Processes in Titan's Upper Atmosphere and Its Induced Magnetosphere. In: Titan from Cassini-Huygens, Brown, R.H., Lebreton J.P., Waite, J.H., Eds
NASA Astrophysics Data System (ADS)
Shiddiqi, Hasbi Ash; Widiyantoro, Sri; Nugraha, Andri Dian; Ramdhan, Mohamad; Wandono, Sutiyono, Handayani, Titi; Nugroho, Hendro
2015-04-01
We have relocated hypocenters of earthquakes occurring in the Molucca collision zone and surrounding region taken from the BMKG catalog using teleseismic double-difference relocation algorithm (teletomoDD). We used P-wave arrival times of local, regional, and teleseismic events recorded at 304 recording stations. Over 7,000 earthquakes were recorded by the BMKG seismographicnetworkin the study region from April, 2009 toJune, 2014. We used a 3D regional-global nested velocity modelresulting fromprevious global tomographystudy. In this study, the3D seismic velocity model was appliedto theIndonesian region, whilethe1D seismicvelocity model (ak135)wasused for regions outside of Indonesia. Our relocation results show a better improvement in travel-time RMS residuals comparedto those of the BMKG catalog.Ourresultsalso show that relocation shifts were dominated intheeast-west direction, whichmaybeinfluenced by theexistingvelocity anomaly related to the reversed V-shaped slabbeneaththestudy region. Our eventrelocation results refine the geometry of slabs beneath the Halmahera and Sangihe arcs.
Shiddiqi, Hasbi Ash E-mail: h.a.shiddiqi@gmail.com; Widiyantoro, Sri; Nugraha, Andri Dian; Ramdhan, Mohamad; Wandono,; Sutiyono,; Handayani, Titi; Nugroho, Hendro
2015-04-24
We have relocated hypocenters of earthquakes occurring in the Molucca collision zone and surrounding region taken from the BMKG catalog using teleseismic double-difference relocation algorithm (teletomoDD). We used P-wave arrival times of local, regional, and teleseismic events recorded at 304 recording stations. Over 7,000 earthquakes were recorded by the BMKG seismographicnetworkin the study region from April, 2009 toJune, 2014. We used a 3D regional-global nested velocity modelresulting fromprevious global tomographystudy. In this study, the3D seismic velocity model was appliedto theIndonesian region, whilethe1D seismicvelocity model (ak135)wasused for regions outside of Indonesia. Our relocation results show a better improvement in travel-time RMS residuals comparedto those of the BMKG catalog.Ourresultsalso show that relocation shifts were dominated intheeast-west direction, whichmaybeinfluenced by theexistingvelocity anomaly related to the reversed V-shaped slabbeneaththestudy region. Our eventrelocation results refine the geometry of slabs beneath the Halmahera and Sangihe arcs.
NASA Astrophysics Data System (ADS)
Hicks, S. P.; Rietbrock, A.; Ryder, I. M.; Miller, M.; Lee, C.
2013-12-01
Knowledge of seismic properties along a subduction megathrust can shed light on the composition and structure of rocks along the fault. By comparing seismic velocity structure with models of interseismic locking, co-seismic slip and afterslip, we can begin to understand how physical properties may affect fault dynamics throughout the subduction seismic cycle. The Maule earthquake, which hit the coast of central Chile in 2010, is the 6th largest earthquake ever recorded, rupturing a 500 x 80 km area of the Chilean megathrust. Published models demonstrate a complex bilateral rupture, with most co-seismic slip occurring to the north of the mainshock epicentre, although significant slip likely stopped short of the trench and the continental Moho. Here, we show a new high-resolution 3D velocity model (vp and vp/vs ratio) of the central Chilean margin Our velocity model is based on manually picked P- and S-wave arrival times from 670 aftershocks recorded by the International Maule Aftershock Deployment (IMAD) network. Seismic properties of the marine forearc are poorly understood in subduction zones, but by incorporating picks from two ocean-bottom seismometer (OBS) networks, we can resolve the velocity structure of the megathrust as far as the trench. In total, the catalogue used for the tomographic inversion yields a total of ~50,000 high quality P- and S-wave picks. We analyse the quality of our model by analysis of the resolution matrix and by testing characteristic models. The 3D velocity model shows the main structures associated within a subduction forearc: the marine forearc basin (vp < 6.0 km/s), continental mantle (vp > 7.5 km/s), and subducting oceanic crust (vp ~ 7.7 km/s). The plate interface is well defined by relocated aftershock seismicity. P-wave velocities along the megathrust range from 6.5 km/s beneath the marine forearc to 7.7 km/s at the intersection of the megathrust with the continental Moho. We infer several high vp anomalies within the South
NASA Astrophysics Data System (ADS)
Calo, M.; Bodin, T.; Yuan, H.; Romanowicz, B. A.; Larmat, C. S.; Maceira, M.
2013-12-01
this work we propose instead to directly tackle the non-linearity of the inverse problem by using stochastic methods to construct a 3D starting model with a good estimate of the depths of the main layering interfaces. We present preliminary results of the construction of such a starting 3D model based on: (1) Regionalizing the study area to define provinces within which lateral variations are smooth; (2) Applying trans-dimensional stochastic inversion (Bodin et al., 2012) to obtain accurate 1D models in each province as well as the corresponding error distribution, constrained by receiver function and surface wave dispersion data as well as the previously constructed 3D model (name), and (3) connecting these models laterally using data-driven smoothing operators to obtain a starting 3D model with errors. References Bodin, T.,et al. 2012, Transdimensional inversion of receiver functions and surface wave dispersion, J. Geophys. Res., 117, B02301, doi:10.1029/2011JB008560. Yuan and Romanowicz, 2013, in revison. Yuan, H., et al. 2011, 3-D shear wave radially and azimuthally anisotropic velocity model of the North American upper mantle. Geophysical Journal International, 184: 1237-1260. doi: 10.1111/j.1365-246X.2010.04901.x Yuan, H. & Romanowicz, B., 2010. Lithospheric layering in the North American Craton, Nature, 466, 1063-1068.
NASA Astrophysics Data System (ADS)
Muench, Thomas; Koch, Manfred; Schlittenhard, Jörg
2010-05-01
On December 5, 2004 a strong earthquake occurred near the city of Waldkirch, about 30 km's north of Freiburg, with a local magnitude of ML = 5.4. This seismic event was one of the strongest observed since the ML = 5.7 'Schwäbische Alb' event of September 3, 1978, 30 years before. In the aftermath of the event several institutions (Bens, BGR, LGBR, LED, SED and NEIC) have attempted to relocate this earthquake that came up with a hypocentral depth range of 9 - 12 km which. In fact, as the exact hypocentral location of the Waldkirch - and other events in the area - namely, the seismic depths, are of utmost importance for the further understanding of the seismotectonics as well as of the seismic hazard in the upper Rhinegraben area, one cannot over stress the necessity for a hypocenter relocation as best as possible. This requires a careful analysis of all factors that may impede an unbiased relocation of such an event. In the present talk we put forward the question whether the Waldkirch seismic event can be relocated with sufficient accuracy by a regional network when, additionally, improved regional 1D- and 3D seismic velocity models for the crust and upper mantle that take into consideration Pn-anisotropy of the upper mantle beneath Germany are employed in the hypocentral determination process. The seismological work starts with a comprehensive analysis of the dataset available for the relocation of the event. By means of traveltime curves a reevaluation of the observed phases is done and it is shown that some of the big observed traveltime residuals are most likely the consequence of wrongly associated phases as well as of the neglect of the anisotropic Pn traveltime correction for the region. Then hypcocenter relocations are done for 1D vertically inhomogeneous and 3D laterally inhomogeneous seismic velocity models, without and with the anisotropic Pn-traveltime correction included. The effects of the - often not well-known - Moho depth and of the VP
ERIC Educational Resources Information Center
Manos, Harry
2016-01-01
Visual aids are important to student learning, and they help make the teacher's job easier. Keeping with the "TPT" theme of "The Art, Craft, and Science of Physics Teaching," the purpose of this article is to show how teachers, lacking equipment and funds, can construct a durable 3-D model reference frame and a model gravity…
NASA Astrophysics Data System (ADS)
Taramón, Jorge M.; Morgan, Jason P.; Pérez-Gussinyé, Marta
2016-04-01
The treatment of far-field boundary conditions is one of the most poorly resolved issues for regional modeling of geodynamic processes. In viscous flow, the choice of far-field boundary conditions often strongly shapes the large-scale structure of a geosimulation. The mantle velocity field along the sidewalls and base of a modeling region is typically much more poorly known than the geometry of past global motions of the surface plates as constrained by global plate motion reconstructions. For regional rifting models it has become routine to apply highly simplified 'plate spreading' or 'uniform rifting' boundary conditions to a 3-D model that limits its ability to simulate the geodynamic evolution of a specific rifted margin. One way researchers are exploring the sensitivity of regional models to uncertain boundary conditions is to use a nested modeling approach in which a global model is used to determine a large-scale flow pattern that is imposed as a constraint along the boundaries of the region to be modeled. Here we explore the utility of a different approach that takes advantage of the ability of finite element models to use unstructured meshes than can embed much higher resolution sub-regions. Here we demonstrate the workflow and code tools that we created to generate this unstructured mesh: solver based on springs, guide-mesh and routines to improve the quality, e.g., closeness to a regular tetrahedron, of the tetrahedral elements of the mesh. Note that the same routines are used to generate a new mesh in the remeshing of a distorted Lagrangian mesh. In our initial project to validate this approach, we create a global spherical shell mesh in which a higher resolution sub-region is created around the nascent South Atlantic Rifting Margin. Global Plate motion BCs and plate boundaries are applied for the time of the onset of rifting, continuing through several 10s of Ma of rifting. Thermal, compositional, and melt-related buoyancy forces are only non
Boatwright, John; Blair, Luke; Catchings, Rufus; Goldman, Mark; Perosi, Fabio; Steedman, Clare
2004-01-01
Campbell (1983) demonstrated that site amplification correlates with depths to the 1.0, 1.5, and 2.5 km/s S-wave velocity horizons. To estimate these depths for the Bay Area stations in the PEER/NGA database, we compare the depths to the 3.2 and 4.4 km/s P-wave velocities in the Brocher and others (1997) 3D velocity model with the depths to these horizons determined from 6 refraction lines shot in the Bay Area from 1991 to 2003. These refraction lines range from two recent 20 km lines that extend from Los Gatos to downtown San Jose, and from downtown San Jose into Alum Rock Park, to two older 200 km lines than run axially from Hollister up the San Francisco Peninsula to Inverness and from Hollister up the East Bay across San Pablo Bay to Santa Rosa. Comparison of these cross-sections with the Brocher and others (1997) model indicates that the 1.5 km/s S-wave horizon, which we correlate with the 3.2 km/s P-wave horizon, is the most reliable horizon that can be extracted from the Brocher and others (1997) velocity model. We determine simple adjustments to bring the Brocher and others (1997) 3.2 and 4.4 km/s P-wave horizons into an average agreement with the refraction results. Then we apply these adjustments to estimate depths to the 1.5 and 2.5 km/s S-wave horizons beneath the strong motion stations in the PEER/NGA database.
NASA Astrophysics Data System (ADS)
Provost, Floriane; Malet, Jean-Philippe; Helmstetter, Agnès; Doubre, Cécile; Gance, Julien
2016-04-01
Microseismicity monitoring has proven to be an important tool for a better understanding of the deformation occurring in slow-sliding landslides. However locating the seismic sources generated by a landslide remains a challenging problem due to (1) the small sizes of the landslide, (b) the heterogenous and time-changing petro-physical properties of the landslide material, (c) the complexity of the recorded signals with unclear discriminations of the wave onsets, and (d) the difficulties to install and maintain a dense seismological network on-site close to the seismic sources. We studied the seismic sources generated by the deformation of the clay-rich Super-Sauze landslide (South French Alps). Previous studies show that the most active zone is the uphill part of the landslide within a zone of 300x300m2. Two seismic antennas have been installed on the sides of this zone and a seismic campaign was conducted to build a 3D velocity model of the area. Calibration shots were performed to test the performance of the location method. We show that the use of a 3D velocity model integrated in a beam forming location method slightly improves the accuracy of the shot location epicenter. However, this approach does not help to interpret with confidence the location of the natural events because the horizontal error remains larger than 50m for more than 50% of the shots. Nevertheless, adding station corrections and constraining the grid search area with additional informations based on the signal and the landslide behavior such as SNR, seismic event typology, and surface kinematics of the landslide allow obtaining reliable results. More than 70% of the calibration shots could be located with a horizontal error of less than 40m. The lack of sensor installed in depth as well as the the lack of calibration shots realized at different depths does not permit us to identify the depth of the sources.
NASA Astrophysics Data System (ADS)
Manos, Harry
2016-03-01
Visual aids are important to student learning, and they help make the teacher's job easier. Keeping with the TPT theme of "The Art, Craft, and Science of Physics Teaching," the purpose of this article is to show how teachers, lacking equipment and funds, can construct a durable 3-D model reference frame and a model gravity well tailored to specific class lessons. Most of the supplies are readily available in the home or at school: rubbing alcohol, a rag, two colors of spray paint, art brushes, and masking tape. The cost of these supplies, if you don't have them, is less than 20.
NASA Astrophysics Data System (ADS)
Lin, F. C.; Schmandt, B.; Tsai, V. C.
2014-12-01
The EarthScope USArray Transportable Array (TA) has provided a great opportunity for imaging the detailed lithospheric structure beneath the continental US. In this presentation, we will report our recent progress on constructing detailed 3D isotropic and anisotropic crustal models of the contiguous US using Rayleigh wave phase velocity and ellipticity measurements across TA. In particular, we will discuss our recent methodology development of extracting short period Rayleigh wave ellipticity, or Rayleigh-wave H/V (horizontal to vertical) amplitude ratios, using multicomponent noise cross-correlations. To retain the amplitude ratio information between vertical and horizontal components, for each station, we perform daily noise pre-processing (temporal normalization and spectrum whitening) simultaneously for all three components. For each station pair, amplitude measurements between cross-correlations of different components (radial-radial, radial-vertical, vertical-radial and vertical-vertical) are then used to determine the Rayleigh-wave H/V ratios at the two station locations. Measurements from all available station pairs are used to determine isotropic and directionally dependent Rayleigh-wave H/V ratios at each location between 8- and 24-second period. The isotropic H/V ratio maps, combined with previous longer period Rayleigh-wave H/V ratio maps from earthquakes and Rayleigh-wave phase velocity maps from both ambient noise and earthquakes, are used to invert for a new 3-D isotropic crustal and upper-mantle model in the western United States. The new model has an outstanding vertical resolution in the upper crust and tradeoffs between different parameters are mitigated. A clear 180-degree periodicity is observed in the directionally dependent H/V ratio measurements for many locations where upper crustal anisotropy is likely strong. Across the US, good correlation is observed between the inferred fast directions in the upper crust and documented maximum
NASA Astrophysics Data System (ADS)
Gerke, Kirill; Vasilyev, Roman; Khirevich, Siarhei; Karsanina, Marina; Collins, Daniel; Korost, Dmitry; Mallants, Dirk
2015-04-01
In this contribution we introduce a novel free software which solves the Stokes equation to obtain velocity fields for low Reynolds-number flows within externally generated 3D pore geometries. Provided with velocity fields, one can calculate permeability for known pressure gradient boundary conditions via Darcy's equation. Finite-difference schemes of 2nd and 4th order of accuracy are used together with an artificial compressibility method to iteratively converge to a steady-state solution of Stokes' equation. This numerical approach is much faster and less computationally demanding than the majority of open-source or commercial softwares employing other algorithms (finite elements/volumes, lattice Boltzmann, etc.) The software consists of two parts: 1) a pre and post-processing graphical interface, and 2) a solver. The latter is efficiently parallelized to use any number of available cores (the speedup on 16 threads was up to 10-12 depending on hardware). Due to parallelization and memory optimization our software can be used to obtain solutions for 300x300x300 voxels geometries on modern desktop PCs. The software was successfully verified by testing it against lattice Boltzmann simulations and analytical solutions. To illustrate the software's applicability for numerous problems in Earth Sciences, a number of case studies have been developed: 1) identifying the representative elementary volume for permeability determination within a sandstone sample, 2) derivation of permeability/hydraulic conductivity values for rock and soil samples and comparing those with experimentally obtained values, 3) revealing the influence of the amount of fine-textured material such as clay on filtration properties of sandy soil. This work was partially supported by RSF grant 14-17-00658 (pore-scale modelling) and RFBR grants 13-04-00409-a and 13-05-01176-a.
Measuring the Stellar Halo Velocity Anisotropy With 3D Kinematics
NASA Astrophysics Data System (ADS)
Cunningham, Emily C.; Deason, Alis J.; Guhathakurta, Puragra; Rockosi, Constance M.; van der Marel, Roeland P.; Sohn, S. Tony
2016-08-01
We present the first measurement of the anisotropy parameter β using 3D kinematic information outside of the solar neighborhood. Our sample consists of 13 Milky Way halo stars with measured proper motions and radial velocities in the line of sight of M31. Proper motions were measured using deep, multi-epoch HST imaging, and radial velocities were measured from Keck II/DEIMOS spectra. We measure β = -0.3-0.9 +0.4, which is consistent with isotropy, and inconsistent with measurements in the solar neighborhood. We suggest that this may be the kinematic signature of a relatively early, massive accretion event, or perhaps several such events.
NASA Astrophysics Data System (ADS)
Geenen, T.; Heister, T.; Van Den Berg, A. P.; Jacobs, M.; Bangerth, W.
2011-12-01
We present high resolution 3D results of the complex mineral phase distribution in the transition zone obtained by numerical modelling of mantle convection. We extend the work by [Jacobs and van den Berg, 2011] to 3D and illustrate the efficiency of adaptive mesh refinement for capturing the complex spatial distribution and sharp phase transitions as predicted by their model. The underlying thermodynamical model is based on lattice dynamics which allows to predict thermophysical properties and seismic wave speeds for the applied magnesium-endmember olivine-pyroxene mineralogical model. The use of 3D geometry allows more realistic prediction of phase distribution and seismic wave speeds resulting from 3D flow processes involving the Earth's transition zone and more significant comparisons with interpretations from seismic tomography and seismic reflectivity studies aimed at the transition zone. Model results are generated with a recently developed geodynamics modeling application based on dealII (www.dealii.org). We extended this model to incorporate both a general thermodynamic model, represented by P,T space tabulated thermophysical properties, and a solution strategy that allows for compressible flow. When modeling compressible flow in the so called truncated anelastic approximation framework we have to adapt the solver strategy that has been proven by several authors to be highly efficient for incompressible flow to incorporate an extra term in the continuity equation. We present several possible solution strategies and discuss their implication in terms of robustness and computational efficiency.
MT3D was first developed by Chunmiao Zheng in 1990 at S.S. Papadopulos & Associates, Inc. with partial support from the U.S. Environmental Protection Agency (USEPA). Starting in 1990, MT3D was released as a pubic domain code from the USEPA. Commercial versions with enhanced capab...
Crowdsourcing Based 3d Modeling
NASA Astrophysics Data System (ADS)
Somogyi, A.; Barsi, A.; Molnar, B.; Lovas, T.
2016-06-01
Web-based photo albums that support organizing and viewing the users' images are widely used. These services provide a convenient solution for storing, editing and sharing images. In many cases, the users attach geotags to the images in order to enable using them e.g. in location based applications on social networks. Our paper discusses a procedure that collects open access images from a site frequently visited by tourists. Geotagged pictures showing the image of a sight or tourist attraction are selected and processed in photogrammetric processing software that produces the 3D model of the captured object. For the particular investigation we selected three attractions in Budapest. To assess the geometrical accuracy, we used laser scanner and DSLR as well as smart phone photography to derive reference values to enable verifying the spatial model obtained from the web-album images. The investigation shows how detailed and accurate models could be derived applying photogrammetric processing software, simply by using images of the community, without visiting the site.
Effects of 3D random correlated velocity perturbations on predicted ground motions
Hartzell, S.; Harmsen, S.; Frankel, A.
2010-01-01
Three-dimensional, finite-difference simulations of a realistic finite-fault rupture on the southern Hayward fault are used to evaluate the effects of random, correlated velocity perturbations on predicted ground motions. Velocity perturbations are added to a three-dimensional (3D) regional seismic velocity model of the San Francisco Bay Area using a 3D von Karman random medium. Velocity correlation lengths of 5 and 10 km and standard deviations in the velocity of 5% and 10% are considered. The results show that significant deviations in predicted ground velocities are seen in the calculated frequency range (≤1 Hz) for standard deviations in velocity of 5% to 10%. These results have implications for the practical limits on the accuracy of scenario ground-motion calculations and on retrieval of source parameters using higher-frequency, strong-motion data.
3-D seismic velocity and attenuation structures in the geothermal field
Nugraha, Andri Dian; Syahputra, Ahmad; Fatkhan,; Sule, Rachmat
2013-09-09
We conducted delay time tomography to determine 3-D seismic velocity structures (Vp, Vs, and Vp/Vs ratio) using micro-seismic events in the geothermal field. The P-and S-wave arrival times of these micro-seismic events have been used as input for the tomographic inversion. Our preliminary seismic velocity results show that the subsurface condition of geothermal field can be fairly delineated the characteristic of reservoir. We then extended our understanding of the subsurface physical properties through determining of attenuation structures (Qp, Qs, and Qs/Qp ratio) using micro-seismic waveform. We combined seismic velocities and attenuation structures to get much better interpretation of the reservoir characteristic. Our preliminary attanuation structures results show reservoir characterization can be more clearly by using the 3-D attenuation model of Qp, Qs, and Qs/Qp ratio combined with 3-D seismic velocity model of Vp, Vs, and Vp/Vs ratio.
NASA Astrophysics Data System (ADS)
Dambrink, Roula; Gunnink, Jan; Stafleu, Jan; de Lange, Ger; Kruiver, Pauline
2016-04-01
The Groningen gas field in the Netherlands is one of the largest gas fields of Europe and has been in production since the 1960's. Due to the progressive depletion of the reservoir, induced seismic activity has increased in recent years. In 2012, an earthquake of magnitude 3.6 initiated further research in prediction and management of risks related to man-induced earthquakes. Last year the government decided to reduce the gas extraction for this reason. One of the topics of concern is the large difference in earthquake-related damage to buildings which, in addition to the distance to the epicenter, appears to be also related to the composition of the shallow subsurface. To improve the spatial distribution of Shear Wave Velocities (Vs) in the shallow subsurface, used for hazard prediction, the Geological Survey of the Netherlands and Deltares constructed a Vs30 map of the upper 30 m of the gas field. In this map a high-resolution geological model (GeoTOP) is combined with seismic cone penetration tests (SCPT) from the area. The GeoTOP model is a 3D voxel model of the upper 50 m, in which each voxel (100x100x0.5 m) is attributed with lithostratigraphy and the most likely lithological class (peat, clay, fine sand, etc.). To obtain statistical distributions (with mean and standard deviation) of Vs for each combination of lithostratigraphical unit and lithoclass, 60 SCPTs were analyzed. In this way, it was possible to assign a specific Vs to each voxel in the model. For each voxel in the stack of voxels that covers the upper 30 m (i.e. 60 voxels), a Vs value was randomly drawn from the statistical distribution of the lithostratigraphical - lithoclass combination it belongs to. The Vs30 for each voxelstack is then calculated using the harmonic mean of the Vs of the 60 voxels. By repeating this procedure 100 times, an (average) Vs30 map and the uncertainty in Vs30 has been constructed. Using the procedure described above we were able to delineate zones with distinct Vs30
Pore-scale intermittent velocity structure underpinning anomalous transport through 3-D porous media
NASA Astrophysics Data System (ADS)
Kang, Peter K.; Anna, Pietro; Nunes, Joao P.; Bijeljic, Branko; Blunt, Martin J.; Juanes, Ruben
2014-09-01
We study the nature of non-Fickian particle transport in 3-D porous media by simulating fluid flow in the intricate pore space of real rock. We solve the full Navier-Stokes equations at the same resolution as the 3-D micro-CT (computed tomography) image of the rock sample and simulate particle transport along the streamlines of the velocity field. We find that transport at the pore scale is markedly anomalous: longitudinal spreading is superdiffusive, while transverse spreading is subdiffusive. We demonstrate that this anomalous behavior originates from the intermittent structure of the velocity field at the pore scale, which in turn emanates from the interplay between velocity heterogeneity and velocity correlation. Finally, we propose a continuous time random walk model that honors this intermittent structure at the pore scale and captures the anomalous 3-D transport behavior at the macroscale.
Visualizing 3D velocity fields near contour surfaces
Max, N.; Crawfis, R.; Grant, C.
1994-03-01
Vector field rendering is difficult in 3D because the vector icons overlap and hide each other. We propose four different techniques for visualizing vector fields only near surfaces. The first uses motion blurred particles in a thickened region around the surface. The second uses a voxel grid to contain integral curves of the vector field. The third uses many antialiased lines through the surface, and the fourth uses hairs sprouting from the surface and then bending in the direction of the vector field. All the methods use the graphite pipeline, allowing real time rotation and interaction, and the first two methods can animate the texture to move in the flow determined by the velocity field.
The Esri 3D city information model
NASA Astrophysics Data System (ADS)
Reitz, T.; Schubiger-Banz, S.
2014-02-01
With residential and commercial space becoming increasingly scarce, cities are going vertical. Managing the urban environments in 3D is an increasingly important and complex undertaking. To help solving this problem, Esri has released the ArcGIS for 3D Cities solution. The ArcGIS for 3D Cities solution provides the information model, tools and apps for creating, analyzing and maintaining a 3D city using the ArcGIS platform. This paper presents an overview of the 3D City Information Model and some sample use cases.
Bayesian 3D velocity field reconstruction with VIRBIUS
NASA Astrophysics Data System (ADS)
Lavaux, Guilhem
2016-03-01
I describe a new Bayesian-based algorithm to infer the full three dimensional velocity field from observed distances and spectroscopic galaxy catalogues. In addition to the velocity field itself, the algorithm reconstructs true distances, some cosmological parameters and specific non-linearities in the velocity field. The algorithm takes care of selection effects, miscalibration issues and can be easily extended to handle direct fitting of e.g. the inverse Tully-Fisher relation. I first describe the algorithm in details alongside its performances. This algorithm is implemented in the VIRBIUS (VelocIty Reconstruction using Bayesian Inference Software) software package. I then test it on different mock distance catalogues with a varying complexity of observational issues. The model proved to give robust measurement of velocities for mock catalogues of 3000 galaxies. I expect the core of the algorithm to scale to tens of thousands galaxies. It holds the promises of giving a better handle on future large and deep distance surveys for which individual errors on distance would impede velocity field inference.
3-D P Wave Velocity Structure of Marmara Region Using Local Earthquake Tomography
NASA Astrophysics Data System (ADS)
Işık, S. E.; Gurbuz, C.
2014-12-01
The 3D P wave velocity model of upper and lower crust of the Marmara Region between 40.200- 41.200N and 26.500- 30.500E is obtained by tomographic inversion (Simulps) of 47034 P wave arrivals of local earthquakes recorded at 90 land stations between October 2009 and December 2012 and 30 OBO stations and 14162 shot arrivals recorded at 35 OBO stations (Seismarmara Survey, 2001). We first obtained a 1D minimum model with Velest code in order to obtain an initial model for 3D inversion with 648 well located earthquakes located within the study area. After several 3D inversion trials we decided to create a more adequate initial model for 3D inversion. Choosing the initial model we estimated the 3D P wave velocity model representing the whole region both for land and sea. The results are tested by making Checkerboard , Restoring Resolution and Characteristic Tests, and the reliable areas of the resulting model is defined in terms of RDE, DWS, SF and Hit count distributions. By taking cross sections from the resulting model we observed the vertical velocity change along profiles crossing both land and sea. All the profiles crossing the basins showed that the high velocities of lower crust make extensions towards the basin area which looks like the force that gives a shape to the basins. These extensions of lower crust towards the basins appeared with an average velocity of 6.3 km/s which might be the result of the deformation due the shearing in the region. It is also interpreted that the development of these high velocities coincide with the development of the basins. Thus, both the basins and the high velocity zones around them might be resulted from the entrance of the NAF into the Marmara Sea and at the same time a shear regime was dominated due to the resistance of the northern Marmara Region (Yılmaz, 2010). The seismicity is observed between 5 km and 15 km after the 3D location of the earthquakes. The locations of the earthquakes improved and the seismogenic zone
Modeling Cellular Processes in 3-D
Mogilner, Alex; Odde, David
2011-01-01
Summary Recent advances in photonic imaging and fluorescent protein technology offer unprecedented views of molecular space-time dynamics in living cells. At the same time, advances in computing hardware and software enable modeling of ever more complex systems, from global climate to cell division. As modeling and experiment become more closely integrated, we must address the issue of modeling cellular processes in 3-D. Here, we highlight recent advances related to 3-D modeling in cell biology. While some processes require full 3-D analysis, we suggest that others are more naturally described in 2-D or 1-D. Keeping the dimensionality as low as possible reduces computational time and makes models more intuitively comprehensible; however, the ability to test full 3-D models will build greater confidence in models generally and remains an important emerging area of cell biological modeling. PMID:22036197
Effects of 3D Velocity and Attenuation in the Tonga-Fiji Subduction Zone
NASA Astrophysics Data System (ADS)
Savage, B.; Wiens, D. A.; Tromp, J.
2005-12-01
The current understanding of a subduction zone's temperature and composition is limited. Much of our recent knowledge of subduction zones comes from earthquake locations, geochemical measurements, and lab based experiments. Recently, two studies of the Tonga-Fiji subduction zone have presented tomographic images of velocity and attenuation (Roth et al., 1999; Zhao et al., 1997). Roth et al. (2000) then combined these two tomographic models of the Tonga-Fiji subduction zone to derive an empirical relationship between changes in velocity and attenuation. This relationship agrees well with two independent, experimental data sets (Jackson et al., 1992; Sato et al., 1989). Using the tomographic velocity model and the empirical relationship between velocity and attenuation we create synthetic seismograms for the Tonga-Fiji subduction zone to test whether a simple increase in velocity accurately depicts this subduction zone. To construct the model we use the tomographic model of Zhao et al. (1997) to create a shear velocity model using a simple Vs/Vp ratio. Following Roth et al. (2000) these tomographic models are combined with the empirical relation between velocity and attenuation to create an attenuation model. The resulting synthetics are compared to recorded data to validate the tomographic velocity model and the empirical relation between velocity and attenuation. Any mismatch in this comparison will provide a basis for further refinement of the tomographic models and the velocity-attenuation relation. The synthetics are created using the SPECFEM3D global code (Komatitsch et al., 2002) with the new addition of a three-dimensional attenuation operator. Attenuation is simulated by a set of standard linear solids over the desired frequency range as described in Liu et al. (1976). Our initial results at a minimum period of 3.3 seconds suggest that the attenuation structure plays a minor role for the present source-receiver geometry. The addition of the 3D attenuation
Radiosity diffusion model in 3D
NASA Astrophysics Data System (ADS)
Riley, Jason D.; Arridge, Simon R.; Chrysanthou, Yiorgos; Dehghani, Hamid; Hillman, Elizabeth M. C.; Schweiger, Martin
2001-11-01
We present the Radiosity-Diffusion model in three dimensions(3D), as an extension to previous work in 2D. It is a method for handling non-scattering spaces in optically participating media. We present the extension of the model to 3D including an extension to the model to cope with increased complexity of the 3D domain. We show that in 3D more careful consideration must be given to the issues of meshing and visibility to model the transport of light within reasonable computational bounds. We demonstrate the model to be comparable to Monte-Carlo simulations for selected geometries, and show preliminary results of comparisons to measured time-resolved data acquired on resin phantoms.
NASA Astrophysics Data System (ADS)
Delph, Jonathan R.; Zandt, George; Beck, Susan L.
2015-12-01
We present a new approach to the joint inversion of surface wave dispersion data and receiver functions by utilizing Common Conversion Point (CCP) stacking to reconcile the different sampling domains of the two datasets. Utilizing CCP stacking allows us to suppress noise in the data by waveform stacking, and correct for backazimuthal variations and complex crustal structure by mapping receiver functions back to their theoretical location. When applied to eastern Turkey, this approach leads to a higher resolution image of the subsurface and clearly delineates different tectonic features in eastern Turkey that were not apparent using other approaches. We observe that the slow seismic velocities near the Karliova Triple Junction correlate to moderate strain rates and high heat flow, which leads to a rheologically weak crust that has allowed for the upward propagation of Miocene and younger volcanics near the triple junction. We find seismically fast, presumably rigid blocks located in the southeastern Anatolian Plate and Arabian Plate are separated by a band of low shear wave velocities that correspond to the East Anatolian Fault Zone, which is consistent with the presence of fluids in the fault zone. We observe that the Arabian Plate has underthrust the Eurasian Plate as far as the northern boundary of the Bitlis Massif, which can explain the high exhumation rates in the Bitlis Massif as a result of slab break-off of the Arabian oceanic lithosphere. We also find a shallow (~ 33 km) anomaly beneath eastern Turkey that we interpret as a localized wedge of mantle that was underthrust by a crustal fragment during the collision of Arabia and Eurasia. These observations are possible because of the high-resolution images obtained by combining common conversion point receiver function stacks with ambient noise dispersion data to create a data-driven three-dimensional shear wave velocity model.
3D model reconstruction of underground goaf
NASA Astrophysics Data System (ADS)
Fang, Yuanmin; Zuo, Xiaoqing; Jin, Baoxuan
2005-10-01
Constructing 3D model of underground goaf, we can control the process of mining better and arrange mining work reasonably. However, the shape of goaf and the laneway among goafs are very irregular, which produce great difficulties in data-acquiring and 3D model reconstruction. In this paper, we research on the method of data-acquiring and 3D model construction of underground goaf, building topological relation among goafs. The main contents are as follows: a) The paper proposed an efficient encoding rule employed to structure the field measurement data. b) A 3D model construction method of goaf is put forward, which by means of combining several TIN (triangulated irregular network) pieces, and an efficient automatic processing algorithm of boundary of TIN is proposed. c) Topological relation of goaf models is established. TIN object is the basic modeling element of goaf 3D model, and the topological relation among goaf is created and maintained by building the topological relation among TIN objects. Based on this, various 3D spatial analysis functions can be performed including transect and volume calculation of goaf. A prototype is developed, which can realized the model and algorithm proposed in this paper.
Rodgers, A; Xie, X
2008-01-09
This project seeks to compute ground motions for large (M>6.5) scenario earthquakes on the Hayward Fault using realistic pseudodynamic ruptures, the USGS three-dimensional (3D) velocity model and anelastic finite difference simulations on parallel computers. We will attempt to bound ground motions by performing simulations with suites of stochastic rupture models for a given scenario on a given fault segment. The outcome of this effort will provide the average, spread and range of ground motions that can be expected from likely large earthquake scenarios. The resulting ground motions will be based on first-principles calculations and include the effects of slip heterogeneity, fault geometry and directivity, however, they will be band-limited to relatively low-frequency (< 1 Hz).
NASA Astrophysics Data System (ADS)
Frankel, A. D.; Stephenson, W. J.; Carver, D.; Odum, J.; Williams, R. A.; Rhea, S.
2010-12-01
We have produced probabilistic seismic hazard maps of Seattle for 1 Hz spectral acceleration, using over five hundred 3D finite-difference simulations of earthquakes on the Seattle fault, Southern Whidbey Island fault, and Cascadia subduction zone, as well as for random deep and shallow earthquakes at various locations. The 3D velocity model was validated by modeling the observed waveforms for the 2001 M6.8 Nisqually earthquake and several smaller events in the region. At these longer periods (≥ 1 sec) that are especially important to the response of buildings of ten stories or higher, seismic waves are strongly influenced by sedimentary basins and rupture directivity. We are investigating how random spatial variations in the 3D velocity model affect the simulated ground motions for M6.7 earthquakes on the Seattle fault. A fractal random variation of shear-wave velocity with a Von Karman correlation function produces spatial variations of peak ground velocity with multiple scale lengths. We find that a 3D velocity model with a 10% standard deviation in shear-wave velocity in the top 1.5 km and 5% standard deviation from 1.5-10 km depth produces variations in peak ground velocities of as much as a factor of two, relative to the case with no random variations. The model with random variations generally reduces the peak ground velocity of the forward rupture directivity pulse for sites near the fault where basin-edge focusing of S-waves occurs. It also tends to reduce the peak velocity of localized areas where basin surface waves are focused. However, the medium with random variations also causes small-scale amplification of ground motions over distances of a few kilometers. We are also evaluating alternative methods of characterizing the aleatory uncertainty in the probabilistic hazard calculations.
3D Simulation of Velocity Profile of Turbulent Flow in Open Channel with Complex Geometry
NASA Astrophysics Data System (ADS)
Kamel, Benoumessad; Ilhem, Kriba; Ali, Fourar; Abdelbaki, Djebaili
Simulation of open channel flow or river flow presents unique challenge to numerical simulators, which is widely used in the applications of computational fluid dynamics. The prediction is extremely difficult because the flow in open channel is usually transient and turbulent, the geometry is irregular and curved, and the free-surface elevation is varying with time. The results from a 3D non-linear k- ɛ turbulence model are presented to investigate the flow structure, the velocity distribution and mass transport process in a meandering compound open channel and a straight open channel. The 3D numerical model for calculating flow is set up in cylinder coordinates in order to calculate the complex boundary channel. The finite volume method is used to disperse the governing equations and the SIMPLE algorithm is applied to acquire the coupling of velocity and pressure. The non-linear k- ɛ turbulent model has good useful value because of taking into account the anisotropy and not increasing the computational time. The main contributions of this study are developing a numerical method that can be applied to predict the flow in river bends with various bend curvatures and different width-depth ratios. This study demonstrates that the 3D non-linear k- ɛ turbulence model can be used for analyzing flow structures, the velocity distribution and pollutant transport in the complex boundary open channel, this model is applicable for real river and wetland problem.
3D Modeling Engine Representation Summary Report
Steven Prescott; Ramprasad Sampath; Curtis Smith; Timothy Yang
2014-09-01
Computers have been used for 3D modeling and simulation, but only recently have computational resources been able to give realistic results in a reasonable time frame for large complex models. This summary report addressed the methods, techniques, and resources used to develop a 3D modeling engine to represent risk analysis simulation for advanced small modular reactor structures and components. The simulations done for this evaluation were focused on external events, specifically tsunami floods, for a hypothetical nuclear power facility on a coastline.
NASA Astrophysics Data System (ADS)
Spica, Zack; Perton, Mathieu; Calò, Marco; Legrand, Denis; Córdoba Montiel, Francisco; Iglesias, Arturo
2016-07-01
This work presents an innovative strategy to enhance the resolution of surface wave tomography obtained from ambient noise cross-correlation (C1) by bridging asynchronous seismic networks through the correlation of coda of correlations (C3). Rayleigh wave group dispersion curves show consistent results between synchronous and asynchronous stations. Rayleigh wave group travel times are inverted to construct velocity-period maps with unprecedented resolution for a region covering Mexico and the southern United States. The resulting period maps are then used to regionalize dispersion curves in order to obtain local 1-D shear velocity models (VS) of the crust and uppermost mantle in every cell of a grid of 0.4°. The 1-D structures are obtained by iteratively adding layers until reaching a given misfit, and a global tomography model is considered as an input for depths below 150 km. Finally, a high-resolution 3-D VS model is obtained from these inversions. The major structures observed in the 3-D model are in agreement with the tectonic-geodynamic features and with previous regional and local studies. It also offers new insights to understand the present and past tectonic evolution of the region.
NASA Astrophysics Data System (ADS)
Spica, Zack; Perton, Mathieu; Calò, Marco; Legrand, Denis; Córdoba-Montiel, Francisco; Iglesias, Arturo
2016-09-01
This work presents an innovative strategy to enhance the resolution of surface wave tomography obtained from ambient noise cross-correlation (C1) by bridging asynchronous seismic networks through the correlation of coda of correlations (C3). Rayleigh wave group dispersion curves show consistent results between synchronous and asynchronous stations. Rayleigh wave group traveltimes are inverted to construct velocity-period maps with unprecedented resolution for a region covering Mexico and the southern United States. The resulting period maps are then used to regionalize dispersion curves in order to obtain local 1-D shear velocity models (VS) of the crust and uppermost mantle in every cell of a grid of 0.4°. The 1-D structures are obtained by iteratively adding layers until reaching a given misfit, and a global tomography model is considered as an input for depths below 150 km. Finally, a high-resolution 3-D VS model is obtained from these inversions. The major structures observed in the 3-D model are in agreement with the tectonic-geodynamic features and with previous regional and local studies. It also offers new insights to understand the present and past tectonic evolution of the region.
BEAMS3D Neutral Beam Injection Model
NASA Astrophysics Data System (ADS)
McMillan, Matthew; Lazerson, Samuel A.
2014-09-01
With the advent of applied 3D fields in Tokamaks and modern high performance stellarators, a need has arisen to address non-axisymmetric effects on neutral beam heating and fueling. We report on the development of a fully 3D neutral beam injection (NBI) model, BEAMS3D, which addresses this need by coupling 3D equilibria to a guiding center code capable of modeling neutral and charged particle trajectories across the separatrix and into the plasma core. Ionization, neutralization, charge-exchange, viscous slowing down, and pitch angle scattering are modeled with the ADAS atomic physics database. Elementary benchmark calculations are presented to verify the collisionless particle orbits, NBI model, frictional drag, and pitch angle scattering effects. A calculation of neutral beam heating in the NCSX device is performed, highlighting the capability of the code to handle 3D magnetic fields. Notice: this manuscript has been authored by Princeton University under Contract Number DE-AC02-09CH11466 with the US Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes.
Analysis of the rupture process of the 1995 Kobe earthquake using a 3D velocity structure
NASA Astrophysics Data System (ADS)
Guo, Yujia; Koketsu, Kazuki; Ohno, Taichi
2013-12-01
A notable feature of the 1995 Kobe (Hyogo-ken Nanbu) earthquake is that violent ground motions occurred in a narrow zone. Previous studies have shown that the origin of such motions can be explained by the 3D velocity structure in this zone. This indicates not only that the 3D velocity structure significantly affects strong ground motions, but also that we should consider its effects in order to determine accurately the rupture process of the earthquake. Therefore, we have performed a joint source inversion of strong-motion, geodetic, and teleseismic data, where 3D Green's functions were calculated for strong-motion and geodetic data in the Osaka basin. Our source model estimates the total seismic moment to be about 2.1 × 1019 N m and the maximum slip reaches 2.9 m near the hypocenter. Although the locations of large slips are similar to those reported by Yoshida et al. (1996), there are quantitative differences between our results and their results due to the differences between the 3D and 1D Green's functions. We have also confirmed that our source model realized a better fit to the strong motion observations, and a similar fit as Yoshida et al. (1996) to the observed static displacements.
Solar abundances and 3D model atmospheres
NASA Astrophysics Data System (ADS)
Ludwig, Hans-Günter; Caffau, Elisabetta; Steffen, Matthias; Bonifacio, Piercarlo; Freytag, Bernd; Cayrel, Roger
2010-03-01
We present solar photospheric abundances for 12 elements from optical and near-infrared spectroscopy. The abundance analysis was conducted employing 3D hydrodynamical (CO5BOLD) as well as standard 1D hydrostatic model atmospheres. We compare our results to others with emphasis on discrepancies and still lingering problems, in particular exemplified by the pivotal abundance of oxygen. We argue that the thermal structure of the lower solar photosphere is very well represented by our 3D model. We obtain an excellent match of the observed center-to-limb variation of the line-blanketed continuum intensity, also at wavelengths shortward of the Balmer jump.
The 3D rocket combustor acoustics model
NASA Technical Reports Server (NTRS)
Priem, Richard J.; Breisacher, Kevin J.
1992-01-01
The theory and procedures for determining the characteristics of pressure oscillations in rocket engines with prescribed burning rate oscillations are presented. Analyses including radial and hub baffles and absorbers can be performed in one, two, and three dimensions. Pressure and velocity oscillations calculated using this procedure are presented for the SSME to show the influence of baffles and absorbers on the burning rate oscillations required to achieve neutral stability. Comparisons are made between the results obtained utilizing 1-D, 2-D, and 3-D assumptions with regards to capturing the physical phenomena of interest and computational requirements.
NASA Astrophysics Data System (ADS)
Pujol, J. M.; Chiu, J. M.
2014-12-01
When applying seismic tomography to local arrival times from an area with a low-velocity sedimentary cover, the effect of the sediments on travel times should be taken into account. If that is not done, the resulting velocity model(s) cannot be assumed to be correct. This fairly obvious statement has been challenged recently by Powell et al. (JGR, 2010), who claimed that the sediments that cover the New Madrid seismic zone (NMSZ, central United States) can be ignored. This claim is examined here and shown to be incorrect. The NMSZ is covered by low-velocity, poorly consolidated sediments (Vp=1.8 km/s, Vs=3), which are underlain by Paleozoic rocks of much higher velocity. In the central NMSZ the sediment thickness varies between about 0.1 and 0.7 km. The JHD analysis of the data collected in that area by a portable network (PANDA) showed that the P- and S-wave station corrections spanned large ranges (0.35 and 0.63 s, respectively, Pujol et al., Eng. Geol., 1997). This study also showed that a Vp/Vs of 3 for the sediments would be too high if the lateral velocity variations were confined to the sedimentary cover. Here we generate synthetic traveltimes for a model with a sedimentary cover having variable depth (as determined from boreholes) underlain by the high-velocity layers in the 1-D model used for the JHD analysis. The synthetic data were generated for the station and event distributions corresponding to the Panda data. The tomographic inversion of the synthetic times produces spurious anomalies in Vp, Vs, and Vp/Vs, from the surface to a depth of 10 km. In addition, the events are mislocated in depth, with errors between 0 and 1 km for most of them. These results should dispel the notion that the effect of the unconsolidated sediments can be ignored. On the other hand, the inversion of the actual Panda data results in velocity anomalies similar to the synthetic anomalies, although larger, which is consistent with the conclusions of Pujol et al. (1997
ERIC Educational Resources Information Center
Bradley, Joan; Farland-Smith, Donna
2010-01-01
Allowing a student to "see" through touch what other students see through a microscope can be a challenging task. Therefore, author Joan Bradley created three-dimensional (3-D) models with one student's visual impairment in mind. They are meant to benefit all students and can be used to teach common high school biology topics, including the…
Debris Dispersion Model Using Java 3D
NASA Technical Reports Server (NTRS)
Thirumalainambi, Rajkumar; Bardina, Jorge
2004-01-01
This paper describes web based simulation of Shuttle launch operations and debris dispersion. Java 3D graphics provides geometric and visual content with suitable mathematical model and behaviors of Shuttle launch. Because the model is so heterogeneous and interrelated with various factors, 3D graphics combined with physical models provides mechanisms to understand the complexity of launch and range operations. The main focus in the modeling and simulation covers orbital dynamics and range safety. Range safety areas include destruct limit lines, telemetry and tracking and population risk near range. If there is an explosion of Shuttle during launch, debris dispersion is explained. The shuttle launch and range operations in this paper are discussed based on the operations from Kennedy Space Center, Florida, USA.
Illustrative visualization of 3D city models
NASA Astrophysics Data System (ADS)
Doellner, Juergen; Buchholz, Henrik; Nienhaus, Marc; Kirsch, Florian
2005-03-01
This paper presents an illustrative visualization technique that provides expressive representations of large-scale 3D city models, inspired by the tradition of artistic and cartographic visualizations typically found in bird"s-eye view and panoramic maps. We define a collection of city model components and a real-time multi-pass rendering algorithm that achieves comprehensible, abstract 3D city model depictions based on edge enhancement, color-based and shadow-based depth cues, and procedural facade texturing. Illustrative visualization provides an effective visual interface to urban spatial information and associated thematic information complementing visual interfaces based on the Virtual Reality paradigm, offering a huge potential for graphics design. Primary application areas include city and landscape planning, cartoon worlds in computer games, and tourist information systems.
Measurements of 3D slip velocities and plasma column lengths of a gliding arc discharge
Zhu, Jiajian; Gao, Jinlong; Ehn, Andreas; Aldén, Marcus; Li, Zhongshan E-mail: alpers@ma.tum.de; Moseev, Dmitry; Kusano, Yukihiro; Salewski, Mirko; Alpers, Andreas E-mail: alpers@ma.tum.de; Gritzmann, Peter; Schwenk, Martin
2015-01-26
A non-thermal gliding arc discharge was generated at atmospheric pressure in an air flow. The dynamics of the plasma column and tracer particles were recorded using two synchronized high-speed cameras. Whereas the data analysis for such systems has previously been performed in 2D (analyzing the single camera image), we provide here a 3D data analysis that includes 3D reconstructions of the plasma column and 3D particle tracking velocimetry based on discrete tomography methods. The 3D analysis, in particular, the determination of the 3D slip velocity between the plasma column and the gas flow, gives more realistic insight into the convection cooling process. Additionally, with the determination of the 3D slip velocity and the 3D length of the plasma column, we give more accurate estimates for the drag force, the electric field strength, the power per unit length, and the radius of the conducting zone of the plasma column.
A 3-D crustal velocity structure across the southeastern Carpathians of Romania
NASA Astrophysics Data System (ADS)
Landes, M.; Hauser, F.; Ritter, J. R. R.; Fielitz, W.; Popa, M.
2003-04-01
The Vrancea zone in the southeastern Carpathians is one of the most active seismic regions in Europe. In order to study the crustal and upper-mantle structure in this region, two seismic refraction experiments were carried out in 1999 and 2001. The 1999 campaign comprised a 320 km long N-S profile and a 80 km long transverse profile. All shots were recorded simultaneously on both profiles. The profile conducted in 2001 extended in E-W direction from the Hungarian border across the Vrancea zone to the Black Sea. A first ever 3-D crustal velocity model of the south-eastern Carpathians within a 115 x 235 km wide region around the Vrancea zone is presented. This model was generated by application of a 3-D refraction and reflection tomography algorithm (Hole 1992, 1995). In order to enhance the model resolution, first arrival data from local earthquakes were also included. The results indicate a high-velocity structure above the Vrancea zone extending from shallow levels to depths of about 11 km. A possible relation to the Trotus and Capidava-Ovidiu faults, which converge to the north of it, is deemed unlikely. However,the existence of the outstanding high velocities may be explained by crystalline basement thrust onto the sub-Carpathian nappes. The high-velocity region is surrounded by the lower velocity Focsani and Brasov basins. The sedimentary succession beneath the southern part of the model area extends to 18 km depth, while in the north sediment thickness varies between 10 and 15 km. Further results of the interface modelling of prominent reflections show that the mid-crustal and Moho interfaces shallow northwards from 30 km to 22 km and from 42 km to 38 km, respectively.
3D velocity structure of upper crust beneath NW Bohemia/Vogtland
NASA Astrophysics Data System (ADS)
Javad Fallahi, Mohammad; Mousavi, Sima; Korn, Michael; Sens-Schönfelder, Christoph; Bauer, Klaus; Rößler, Dirk
2013-04-01
The 3D structure of the upper crust beneath west Bohemia/Vogtland region, analyzed with travel time tomography and ambient noise surface wave tomography using existing data. This region is characterized by a series of phenomena like occurrence of repeated earthquake swarms, surface exhalation, CO2 enriched fluids, mofettes, mineral springs and enhanced heat flow, and has been proposed as an excellent location for an ICDP drilling project targeted to a better understanding of the crust in an active magmatic environment. We performed a 3D tomography using P-and S-wave travel times of local earthquakes and explosions. The data set were taken from permanent and temporary seismic networks in Germany and Czech Republic from 2000 to 2010, as well as active seismic experiments like Celebration 2000 and quarry blasts. After picking P and S wave arrival times, 399 events which were recorded by 9 or more stations and azimuthal gap<160° were selected for inversion. A simultaneous inversion of P and S wave 1D velocity models together with relocations of hypocenters and station corrections was performed. The obtained minimum 1D velocity model was used as starting model for the 3D Vp and Vp/Vs velocity models. P and S wave travel time tomography employs damped least-square method and ray tracing by pseudo-bending algorithm. For model parametrization different cell node spacings have been tested to evaluate the resolution in each node. Synthetic checkerboard tests have been done to check the structural resolution. Then Vp and Vp/Vs in the preferred 3D grid model have been determined. Earthquakes locations in iteration process change till the hypocenter adjustments and travel time residuals become smaller than the defined threshold criteria. Finally the analysis of the resolution depicts the well resolved features for interpretation. We observed lower Vp/Vs ratio in depth of 5-10 km close to the foci of earthquake swarms and higher Vp/Vs ratio is observed in Saxoturingian zone and
Integrated Biogeomorphological Modeling Using Delft3D
NASA Astrophysics Data System (ADS)
Ye, Q.; Jagers, B.
2011-12-01
The skill of numerical morphological models has improved significantly from the early 2D uniform, total load sediment models (with steady state or infrequent wave updates) to recent 3D hydrodynamic models with multiple suspended and bed load sediment fractions and bed stratigraphy (online coupled with waves). Although there remain many open questions within this combined field of hydro- and morphodynamics, we observe an increasing need to include biological processes in the overall dynamics. In riverine and inter-tidal environments, there is often an important influence by riparian vegetation and macrobenthos. Over the past decade more and more researchers have started to extend the simulation environment with wrapper scripts and other quick code hacks to estimate their influence on morphological development in coastal, estuarine and riverine environments. Although one can in this way quickly analyze different approaches, these research tools have generally not been designed with reuse, performance and portability in mind. We have now implemented a reusable, flexible, and efficient two-way link between the Delft3D open source framework for hydrodynamics, waves and morphology, and the water quality and ecology modules. The same link will be used for 1D, 2D and 3D modeling on networks and both structured and unstructured grids. We will describe the concepts of the overall system, and illustrate it with some first results.
Sensing and compressing 3-D models
Krumm, J.
1998-02-01
The goal of this research project was to create a passive and robust computer vision system for producing 3-D computer models of arbitrary scenes. Although the authors were unsuccessful in achieving the overall goal, several components of this research have shown significant potential. Of particular interest is the application of parametric eigenspace methods for planar pose measurement of partially occluded objects in gray-level images. The techniques presented provide a simple, accurate, and robust solution to the planar pose measurement problem. In addition, the representational efficiency of eigenspace methods used with gray-level features were successfully extended to binary features, which are less sensitive to illumination changes. The results of this research are presented in two papers that were written during the course of this project. The papers are included in sections 2 and 3. The first section of this report summarizes the 3-D modeling efforts.
NASA Astrophysics Data System (ADS)
Berchtold, Waldemar; Schäfer, Marcel; Rettig, Michael; Steinebach, Martin
2014-02-01
3D models and applications are of utmost interest in both science and industry. With the increment of their usage, their number and thereby the challenge to correctly identify them increases. Content identification is commonly done by cryptographic hashes. However, they fail as a solution in application scenarios such as computer aided design (CAD), scientific visualization or video games, because even the smallest alteration of the 3D model, e.g. conversion or compression operations, massively changes the cryptographic hash as well. Therefore, this work presents a robust hashing algorithm for 3D mesh data. The algorithm applies several different bit extraction methods. They are built to resist desired alterations of the model as well as malicious attacks intending to prevent correct allocation. The different bit extraction methods are tested against each other and, as far as possible, the hashing algorithm is compared to the state of the art. The parameters tested are robustness, security and runtime performance as well as False Acceptance Rate (FAR) and False Rejection Rate (FRR), also the probability calculation of hash collision is included. The introduced hashing algorithm is kept adaptive e.g. in hash length, to serve as a proper tool for all applications in practice.
3D modeling of optically challenging objects.
Park, Johnny; Kak, Avinash
2008-01-01
We present a system for constructing 3D models of real-world objects with optically challenging surfaces. The system utilizes a new range imaging concept called multi-peak range imaging, which stores multiple candidates of range measurements for each point on the object surface. The multiple measurements include the erroneous range data caused by various surface properties that are not ideal for structured-light range sensing. False measurements generated by spurious reflections are eliminated by applying a series of constraint tests. The constraint tests based on local surface and local sensor visibility are applied first to individual range images. The constraint tests based on global consistency of coordinates and visibility are then applied to all range images acquired from different viewpoints. We show the effectiveness of our method by constructing 3D models of five different optically challenging objects. To evaluate the performance of the constraint tests and to examine the effects of the parameters used in the constraint tests, we acquired the ground truth data by painting those objects to suppress the surface-related properties that cause difficulties in range sensing. Experimental results indicate that our method significantly improves upon the traditional methods for constructing reliable 3D models of optically challenging objects. PMID:18192707
Fallon FORGE 3D Geologic Model
Doug Blankenship
2016-03-01
An x,y,z scattered data file for the 3D geologic model of the Fallon FORGE site. Model created in Earthvision by Dynamic Graphic Inc. The model was constructed with a grid spacing of 100 m. Geologic surfaces were extrapolated from the input data using a minimum tension gridding algorithm. The data file is tabular data in a text file, with lithology data associated with X,Y,Z grid points. All the relevant information is in the file header (the spatial reference, the projection etc.) In addition all the fields in the data file are identified in the header.
3D Models of Symbiotic Binaries
NASA Astrophysics Data System (ADS)
Mohamed, S.; Booth, R.; Podsiadlowski, Ph.; Ramstedt, S.; Vlemmings, W.; Maercker, M.
2015-12-01
Symbiotic binaries consist of a cool, mass-losing giant and an accreting, compact companion. We present 3D Smoothed Particle Hydrodynamics (SPH) models of two such interacting binaries, RS Oph and Mira AB. RS Oph is also a recurrent nova system, thus we model multiple quiescent mass transfer-nova outburst cycles. The resulting circumstellar structures of both systems are highly complex with the formation of spirals, arcs, shells, equatorial and bipolar outflows. We compare the models to recent observations and discuss the implications of our results for related systems, e.g., bipolar nebulae and jets, chemically peculiar stars, and the progenitors of Type Ia supernovae.
Inferential modeling of 3D chromatin structure
Wang, Siyu; Xu, Jinbo; Zeng, Jianyang
2015-01-01
For eukaryotic cells, the biological processes involving regulatory DNA elements play an important role in cell cycle. Understanding 3D spatial arrangements of chromosomes and revealing long-range chromatin interactions are critical to decipher these biological processes. In recent years, chromosome conformation capture (3C) related techniques have been developed to measure the interaction frequencies between long-range genome loci, which have provided a great opportunity to decode the 3D organization of the genome. In this paper, we develop a new Bayesian framework to derive the 3D architecture of a chromosome from 3C-based data. By modeling each chromosome as a polymer chain, we define the conformational energy based on our current knowledge on polymer physics and use it as prior information in the Bayesian framework. We also propose an expectation-maximization (EM) based algorithm to estimate the unknown parameters of the Bayesian model and infer an ensemble of chromatin structures based on interaction frequency data. We have validated our Bayesian inference approach through cross-validation and verified the computed chromatin conformations using the geometric constraints derived from fluorescence in situ hybridization (FISH) experiments. We have further confirmed the inferred chromatin structures using the known genetic interactions derived from other studies in the literature. Our test results have indicated that our Bayesian framework can compute an accurate ensemble of 3D chromatin conformations that best interpret the distance constraints derived from 3C-based data and also agree with other sources of geometric constraints derived from experimental evidence in the previous studies. The source code of our approach can be found in https://github.com/wangsy11/InfMod3DGen. PMID:25690896
NASA Astrophysics Data System (ADS)
Yao, H.; Fang, H.; Li, C.; Liu, Y.; Zhang, H.; van der Hilst, R. D.; Huang, Y. C.
2014-12-01
Ambient noise tomography has provided essential constraints on crustal and uppermost mantle shear velocity structure in global seismology. Recent studies demonstrate that high frequency (e.g., ~ 1 Hz) surface waves between receivers at short distances can be successfully retrieved from ambient noise cross-correlation and then be used for imaging near surface or shallow crustal shear velocity structures. This approach provides important information for strong ground motion prediction in seismically active area and overburden structure characterization in oil and gas fields. Here we propose a new tomographic method to invert all surface wave dispersion data for 3-D variations of shear wavespeed without the intermediate step of phase or group velocity maps.The method uses frequency-dependent propagation paths and a wavelet-based sparsity-constrained tomographic inversion. A fast marching method is used to compute, at each period, surface wave traveltimes and ray paths between sources and receivers. This avoids the assumption of great-circle propagation that is used in most surface wave tomographic studies, but which is not appropriate in complex media. The wavelet coefficients of the velocity model are estimated with an iteratively reweighted least squares (IRLS) algorithm, and upon iterations the surface wave ray paths and the data sensitivity matrix are updated from the newly obtained velocity model. We apply this new method to determine the 3-D near surface wavespeed variations in the Taipei basin of Taiwan, Hefei urban area and a shale and gas production field in China using the high-frequency interstation Rayleigh wave dispersion data extracted from ambient noisecross-correlation. The results reveal strong effects of off-great-circle propagation of high-frequency surface waves in these regions with above 30% shear wavespeed variations. The proposed approach is more efficient and robust than the traditional two-step surface wave tomography for imaging complex
3D MHD Models of Active Region Loops
NASA Technical Reports Server (NTRS)
Ofman, Leon
2004-01-01
Present imaging and spectroscopic observations of active region loops allow to determine many physical parameters of the coronal loops, such as the density, temperature, velocity of flows in loops, and the magnetic field. However, due to projection effects many of these parameters remain ambiguous. Three dimensional imaging in EUV by the STEREO spacecraft will help to resolve the projection ambiguities, and the observations could be used to setup 3D MHD models of active region loops to study the dynamics and stability of active regions. Here the results of 3D MHD models of active region loops are presented, and the progress towards more realistic 3D MHD models of active regions. In particular the effects of impulsive events on the excitation of active region loop oscillations, and the generation, propagations and reflection of EIT waves are shown. It is shown how 3D MHD models together with 3D EUV observations can be used as a diagnostic tool for active region loop physical parameters, and to advance the science of the sources of solar coronal activity.
Park, J.; Morgan, J.K.; Zelt, C.A.; Okubo, P.G.; Peters, L.; Benesh, N.
2007-01-01
We present a 3-D P-wave velocity model of the combined subaerial and submarine portions of the southeastern part of the Island of Hawaii, based on first-arrival seismic tomography of marine airgun shots recorded by the onland seismic network. Our model shows that high-velocity materials (6.5-7.0??km/s) lie beneath Kilauea's summit, Koae fault zone, and the upper Southwest Rift Zone (SWRZ) and upper and middle East Rift Zone (ERZ), indicative of magma cumulates within the volcanic edifice. A separate high-velocity body of 6.5-6.9??km/s within Kilauea's lower ERZ and upper Puna Ridge suggests a distinct body of magma cumulates, possibly connected to the summit magma cumulates at depth. The two cumulate bodies within Kilauea's ERZ may have undergone separate ductile flow seaward, influencing the submarine morphology of Kilauea's south flank. Low velocities (5.0-6.3??km/s) seaward of Kilauea's Hilina fault zone, and along Mauna Loa's seaward facing Kao'iki fault zone, are attributed to thick piles of volcaniclastic sediments deposited on the submarine flanks. Loihi seamount shows high-velocity anomalies beneath the summit and along the rift zones, similar to the interpreted magma cumulates below Mauna Loa and Kilauea volcanoes, and a low-velocity anomaly beneath the oceanic crust, probably indicative of melt within the upper mantle. Around Kilauea's submarine flank, a high-velocity anomaly beneath the outer bench suggests the presence of an ancient seamount that may obstruct outward spreading of the flank. Mauna Loa's southeast flank is also marked by a large, anomalously high-velocity feature (7.0-7.4??km/s), interpreted to define an inactive, buried volcanic rift zone, which might provide a new explanation for the westward migration of Mauna Loa's current SWRZ and the growth of Kilauea's SWRZ. ?? 2007 Elsevier B.V. All rights reserved.
Imaging 3D seismic velocity along the seismogenic zone of Algarve region (southern Portugal)
NASA Astrophysics Data System (ADS)
Rocha, João.; Bezzeghoud, Mourad; Caldeira, Bento; Dias, Nuno; Borges, José; Matias, Luís.; Dorbath, Catherine; Carrilho, Fernando
2010-05-01
The present seismic tomographic study is focused around Algarve region, in South of Portugal. To locate the seismic events and find the local velocity structure of epicentral area, the P and S arrival times at 38 stations are used. The data used in this study were obtained during the Algarve campaign which worked from January/2006 to July/2007. The preliminary estimate of origin times and hypocentral coordinates are determined by the Hypoinverse program. Linearized inversion procedure was applied to comprise the following two steps: 1) finding the minimum 1D velocity model using Velest and 2) simultaneous relocation of hypocenters and determination of local velocity structure. The velocity model we have reached is a 10 layer model which gave the lowest RMS, after several runnings of eight different velocity models that we used "a priori". The model parameterization assumes a continuous velocity field between 4.5 km/s and 7.0 km/s until 30 km depth. The earth structure is represented in 3D by velocity at discrete points, and velocity at any intervening point is determined by linear interpolation among the surrounding eight grid points. A preliminary analysis of the resolution capabilities of the dataset, based on the Derivative Weight Sum (DWS) distribution, shows that the velocity structure is better resolved in the West part of the region between the surface to15 km. The resulting tomographic image has a prominent low-velocity anomaly that shows a maximum decrease in P-wave velocity in the first 12 kms in the studied region. We also identified the occurrence of local seismic events of reduced magnitude not catalogued, in the neighbourhood of Almodôvar (low Alentejo). The spatial distribution of epicentres defines a NE-SW direction that coincides with the strike of the mapped geological faults of the region and issued from photo-interpretation. Is still expectable to refine the seismicity of the region of Almodôvar and establish more rigorously its role in the
3D Model of Surfactant Replacement Therapy
NASA Astrophysics Data System (ADS)
Grotberg, James; Tai, Cheng-Feng; Filoche, Marcel
2015-11-01
Surfactant Replacement Therapy (SRT) involves instillation of a liquid-surfactant mixture directly into the lung airway tree. Though successful in neonatal applications, its use in adults had early success followed by failure. We present the first mathematical model of 3D SRT where a liquid plug propagates through the tree from forced inspiration. In two separate modeling steps, the plug first deposits a coating film on the airway wall which subtracts from its volume, a ``coating cost''. Then the plug splits unevenly at the airway bifurcation due to gravity. The steps are repeated until a plug ruptures or reaches the tree endpoint alveoli/acinus. The model generates 3D images of the resulting acinar distribution and calculates two global indexes, efficiency and homogeneity. Simulating published literature, the earlier successful adult SRT studies show comparatively good index values, while the later failed studies do not. Those unsuccessful studies used smaller dose volumes with higher concentration mixtures, apparently assuming a well mixed compartment. The model shows that adult lungs are not well mixed in SRT due to the coating cost and gravity effects. Returning to the higher dose volume protocols could save many thousands of lives annually in the US. Supported by NIH Grants HL85156, HL84370 and Agence Nationale de la Recherche, ANR no. 2010-BLAN-1119-05.
MOSSFRAC: An anisotropic 3D fracture model
Moss, W C; Levatin, J L
2006-08-14
Despite the intense effort for nearly half a century to construct detailed numerical models of plastic flow and plastic damage accumulation, models for describing fracture, an equally important damage mechanism still cannot describe basic fracture phenomena. Typical fracture models set the stress tensor to zero for tensile fracture and set the deviatoric stress tensor to zero for compressive fracture. One consequence is that the simple case of the tensile fracture of a cylinder under combined compressive radial and tensile axial loads is not modeled correctly. The experimental result is a cylinder that can support compressive radial loads, but no axial load, whereas, the typical numerical result is a cylinder with all stresses equal to zero. This incorrect modeling of fracture locally also has a global effect, because material that is fracturing produces stress release waves, which propagate from the fracture and influence the surrounding material. Consequently, it would be useful to have a model that can describe the stress relief and the resulting anisotropy due to fracture. MOSSFRAC is a material model that simulates three-dimensional tensile and shear fracture in initially isotropic elastic-plastic materials, although its framework is also amenable to initially anisotropic materials. It differs from other models by accounting for the effects of cracks on the constitutive response of the material, so that the previously described experiment, as well as complicated fracture scenarios are simulated more accurately. The model is implemented currently in the LLNL hydrocodes DYNA3D, PARADYN, and ALE3D. The purpose of this technical note is to present a complete qualitative description of the model and quantitative descriptions of salient features.
A 3-D crustal velocity structure across the Variscides of southwest Ireland
NASA Astrophysics Data System (ADS)
Landes, M.; Readman, P. W.; O'Reilly, B. M.; Shannon, P. M.
2003-04-01
In the VARNET-96 experiment three seismic refraction profiles were acquired to examine the crustal structure in the south-west of Ireland. The shotpoint geometry allowed for both in-line and off-line fan shot recordings on the three profiles. Results of 3-D inversion modelling illustrate that there is pervasive lateral heterogeneity of the sedimentary and crustal velocity structure south of the Shannon Estuary. Palaeozoic strata at the south coast are about 5-6 km thick associated with the sedimentary infill of the Munster and South Munster Basins. To the north, shallow upper crust in the vicinity of the Killarney-Mallow Fault Zone is followed by a 3-4 km thick sedimentary succession in the Dingle-Shannon Basin. A zone of high-velocity upper crust (6.4-6.6 km/s) beneath the South Munster Basin correlates with a gravity high between the Kenmare-Killarney and the Leinster Granite gravity lows. Other high-velocity zones beneath Dingle Bay and the Kenmare River region may be associated with the deep traces of the Killarney-Mallow Fault Zone and the Cork-Kenmare Line. The 3-D velocity model was taken as a basis for the computation of PmP reflected arrivals from the crust-mantle boundary. The Moho depth varies from 28-29 km at the south coast to 32-33 km beneath the Dingle-Shannon Basin. Pervasive Variscan deformation appears to be confined to the sedimentary and upper crustal structure thus supporting a thin-skinned tectonic model for Variscan deformation. Deep-crustal variations only occur where they can be correlated with major tectonic features such as the Caledonian Iapetus Suture near the Shannon Estuary. The shallowing of the Moho towards the coast may result from Mesozoic crustal extension in the adjacent offshore sedimentary basins.
3D Stratigraphic Modeling of Central Aachen
NASA Astrophysics Data System (ADS)
Dong, M.; Neukum, C.; Azzam, R.; Hu, H.
2010-05-01
Since 1980s, advanced computer hardware and software technologies, as well as multidisciplinary research have provided possibilities to develop advanced three dimensional (3D) simulation software for geosciences application. Some countries, such as USA1) and Canada2) 3), have built up regional 3D geological models based on archival geological data. Such models have played huge roles in engineering geology2), hydrogeology2) 3), geothermal industry1) and so on. In cooperating with the Municipality of Aachen, the Department of Engineering Geology of RWTH Aachen University have built up a computer-based 3D stratigraphic model of 50 meter' depth for the center of Aachen, which is a 5 km by 7 km geologically complex area. The uncorrelated data from multi-resources, discontinuous nature and unconformable connection of the units are main challenges for geological modeling in this area. The reliability of 3D geological models largely depends on the quality and quantity of data. Existing 1D and 2D geological data were collected, including 1) approximately 6970 borehole data of different depth compiled in Microsoft Access database and MapInfo database; 2) a Digital Elevation Model (DEM); 3) geological cross sections; and 4) stratigraphic maps in 1m, 2m and 5m depth. Since acquired data are of variable origins, they were managed step by step. The main processes are described below: 1) Typing errors of borehole data were identified and the corrected data were exported to Variowin2.2 to distinguish duplicate points; 2) The surface elevation of borehole data was compared to the DEM, and differences larger than 3m were eliminated. Moreover, where elevation data missed, it was read from the DEM; 3) Considerable data were collected from municipal constructions, such as residential buildings, factories, and roads. Therefore, many boreholes are spatially clustered, and only one or two representative points were picked out in such areas; After above procedures, 5839 boreholes with -x
Spatial parallelism of a 3D finite difference, velocity-stress elastic wave propagation code
Minkoff, S.E.
1999-12-01
Finite difference methods for solving the wave equation more accurately capture the physics of waves propagating through the earth than asymptotic solution methods. Unfortunately, finite difference simulations for 3D elastic wave propagation are expensive. The authors model waves in a 3D isotropic elastic earth. The wave equation solution consists of three velocity components and six stresses. The partial derivatives are discretized using 2nd-order in time and 4th-order in space staggered finite difference operators. Staggered schemes allow one to obtain additional accuracy (via centered finite differences) without requiring additional storage. The serial code is most unique in its ability to model a number of different types of seismic sources. The parallel implementation uses the MPI library, thus allowing for portability between platforms. Spatial parallelism provides a highly efficient strategy for parallelizing finite difference simulations. In this implementation, one can decompose the global problem domain into one-, two-, and three-dimensional processor decompositions with 3D decompositions generally producing the best parallel speedup. Because I/O is handled largely outside of the time-step loop (the most expensive part of the simulation) the authors have opted for straight-forward broadcast and reduce operations to handle I/O. The majority of the communication in the code consists of passing subdomain face information to neighboring processors for use as ghost cells. When this communication is balanced against computation by allocating subdomains of reasonable size, they observe excellent scaled speedup. Allocating subdomains of size 25 x 25 x 25 on each node, they achieve efficiencies of 94% on 128 processors. Numerical examples for both a layered earth model and a homogeneous medium with a high-velocity blocky inclusion illustrate the accuracy of the parallel code.
Spatial Parallelism of a 3D Finite Difference, Velocity-Stress Elastic Wave Propagation Code
MINKOFF,SUSAN E.
1999-12-09
Finite difference methods for solving the wave equation more accurately capture the physics of waves propagating through the earth than asymptotic solution methods. Unfortunately. finite difference simulations for 3D elastic wave propagation are expensive. We model waves in a 3D isotropic elastic earth. The wave equation solution consists of three velocity components and six stresses. The partial derivatives are discretized using 2nd-order in time and 4th-order in space staggered finite difference operators. Staggered schemes allow one to obtain additional accuracy (via centered finite differences) without requiring additional storage. The serial code is most unique in its ability to model a number of different types of seismic sources. The parallel implementation uses the MP1 library, thus allowing for portability between platforms. Spatial parallelism provides a highly efficient strategy for parallelizing finite difference simulations. In this implementation, one can decompose the global problem domain into one-, two-, and three-dimensional processor decompositions with 3D decompositions generally producing the best parallel speed up. Because i/o is handled largely outside of the time-step loop (the most expensive part of the simulation) we have opted for straight-forward broadcast and reduce operations to handle i/o. The majority of the communication in the code consists of passing subdomain face information to neighboring processors for use as ''ghost cells''. When this communication is balanced against computation by allocating subdomains of reasonable size, we observe excellent scaled speed up. Allocating subdomains of size 25 x 25 x 25 on each node, we achieve efficiencies of 94% on 128 processors. Numerical examples for both a layered earth model and a homogeneous medium with a high-velocity blocky inclusion illustrate the accuracy of the parallel code.
Subduction zone guided waves: 3D modelling and attenuation effects
NASA Astrophysics Data System (ADS)
Garth, T.; Rietbrock, A.
2013-12-01
Waveform modelling is an important tool for understanding complex seismic structures such as subduction zone waveguides. These structures are often simplified to 2D structures for modelling purposes to reduce computational costs. In the case of subduction zone waveguide affects, 2D models have shown that dispersed arrivals are caused by a low velocity waveguide, inferred to be subducted oceanic crust and/or hydrated outer rise normal faults. However, due to the 2D modelling limitations the inferred seismic properties such as velocity contrast and waveguide thickness are still debated. Here we test these limitations with full 3D waveform modelling. For waveguide effects to be observable the waveform must be accurately modelled to relatively high frequencies (> 2 Hz). This requires a small grid spacing due to the high seismic velocities present in subduction zones. A large area must be modelled as well due to the long propagation distances (400 - 600 km) of waves interacting with subduction zone waveguides. The combination of the large model area and small grid spacing required means that these simulations require a large amount of computational resources, only available at high performance computational centres like the UK National super computer HECTOR (used in this study). To minimize the cost of modelling for such a large area, the width of the model area perpendicular to the subduction trench (the y-direction) is made as small as possible. This reduces the overall volume of the 3D model domain. Therefore the wave field is simulated in a model ';corridor' of the subduction zone velocity structure. This introduces new potential sources of error particularly from grazing wave side reflections in the y-direction. Various dampening methods are explored to reduce these grazing side reflections, including perfectly matched layers (PML) and more traditional exponential dampening layers. Defining a corridor model allows waveguide affects to be modelled up to at least 2
3D Modeling of Equatorial Plasma Bubbles
NASA Astrophysics Data System (ADS)
Huba, Joseph; Joyce, Glenn; Krall, Jonathan
2011-10-01
Post-sunset ionospheric irregularities in the equatorial F region were first observed by Booker and Wells (1938) using ionosondes. This phenomenon has become known as equatorial spread F (ESF). During ESF the equatorial ionosphere becomes unstable because of a Rayleigh-Taylor-like instability: large scale (10s km) electron density ``bubbles'' can develop and rise to high altitudes (1000 km or greater at times). Understanding and modeling ESF is important because of its impact on space weather: it causes radio wave scintillation that degrades communication and navigation systems. In fact, it is the focus of of the Air Force Communications/Navigation Outage Forecast Satellite (C/NOFS) mission. We will describe 3D simulation results from the NRL ionosphere models SAMI3 and SAMI3/ESF of this phenomenon. In particular, we will examine the causes of the day-to-day ariability of ESF which is an unresolved problem at this time. Post-sunset ionospheric irregularities in the equatorial F region were first observed by Booker and Wells (1938) using ionosondes. This phenomenon has become known as equatorial spread F (ESF). During ESF the equatorial ionosphere becomes unstable because of a Rayleigh-Taylor-like instability: large scale (10s km) electron density ``bubbles'' can develop and rise to high altitudes (1000 km or greater at times). Understanding and modeling ESF is important because of its impact on space weather: it causes radio wave scintillation that degrades communication and navigation systems. In fact, it is the focus of of the Air Force Communications/Navigation Outage Forecast Satellite (C/NOFS) mission. We will describe 3D simulation results from the NRL ionosphere models SAMI3 and SAMI3/ESF of this phenomenon. In particular, we will examine the causes of the day-to-day ariability of ESF which is an unresolved problem at this time. Research supported by ONR.
NASA Astrophysics Data System (ADS)
Asmerom, B. B.; Chiu, J.; Pujol, J.; Smalley, R.
2010-12-01
Lateral crustal velocity variations across the Andean Foreland in San Juan Argentina are explored by joint hypocentral determination (JHD) analysis and 3D velocity inversion. JHD results show consistent positive station corrections beneath Precordillera and negative station corrections beneath Pie de Palo, corresponding to regions of low and high velocity, respectively. These observations are supported by the results from the 3D velocity inversion. A 20% increase in velocity is observed from the Precordilleras in the west to Pie de Palo in the east. The tomography result also reveals a narrow east dipping and NNE trending high velocity anomalous zone bisecting the southern half of Pie de Palo. This anomalous zone was previously identified by a magnetic study and was interpreted to represent the structure corresponding to the Grenvillian Precordillera-Pie de Palo tectonic boundary zone. Finally, P and S station corrections are calculated from the synthetic travel time obtained by using the resultant 3D P- and S- wave velocity model. The observed pattern and magnitude of the P- and S-wave station corrections are recovered successfully from the synthetic calculation, indicating that the resultant 3D velocity model is close to the real earth structure in the Andean Foreland region. Relocation of all intermediate events from the flat subducting slab using this newly acquired 3D velocity model shows a significant change in the slab geometry. The relocated hypocenter distribution is more clustered than previous studies obtained using a 1D model. The slab is simply flat and it resumes a normal subduction angle towards the east of the study area.
Reservoir geology using 3D modelling tools
Dubrule, O.; Samson, P.; Segonds, D.
1996-12-31
The last decade has seen tremendous developments in the area of quantitative geological modelling. These developments have a significant impact on the current practice of constructing reservoir models. A structural model can first be constructed on the basis of depth-converted structural interpretations produced on a seismic interpretation workstation. Surfaces and faults can be represented as geological objects, and interactively modified. Once the tectonic framework has been obtained, intermediate stratigraphic surfaces can be constructed between the main structural surfaces. Within each layer, reservoir attributes can be represented using various techniques. Examples show how the distribution of different facies (i.e. from fine to coarse grain) can be represented, or how various depositional units (for instance channels, crevasses and lobes in a turbidite setting) can be modelled as geological {open_quotes}objects{close_quotes} with complex geometries. Elf Aquitaine, in close co-operation with the GOCAD project in Nancy (France) is investigating how geological models can be made more realistic by developing interactive functionalities. Examples show that, contrary to standard deterministic or geostatistical modelling techniques (which tend to be difficult to control) the use of new 3D tools allows the geologist to interactively modify geological surfaces (including faults) or volumetric properties. Thus, the sensitivity of various economic parameters (oil in place, connected volumes, reserves) to major geological uncertainties can be evaluated. It is argued that future breakthroughs in geological modelling techniques are likely to happen in the development of interactive approaches rather than in the research of new mathematical algorithms.
Reservoir geology using 3D modelling tools
Dubrule, O. ); Samson, P. ); Segonds, D. )
1996-01-01
The last decade has seen tremendous developments in the area of quantitative geological modelling. These developments have a significant impact on the current practice of constructing reservoir models. A structural model can first be constructed on the basis of depth-converted structural interpretations produced on a seismic interpretation workstation. Surfaces and faults can be represented as geological objects, and interactively modified. Once the tectonic framework has been obtained, intermediate stratigraphic surfaces can be constructed between the main structural surfaces. Within each layer, reservoir attributes can be represented using various techniques. Examples show how the distribution of different facies (i.e. from fine to coarse grain) can be represented, or how various depositional units (for instance channels, crevasses and lobes in a turbidite setting) can be modelled as geological [open quotes]objects[close quotes] with complex geometries. Elf Aquitaine, in close co-operation with the GOCAD project in Nancy (France) is investigating how geological models can be made more realistic by developing interactive functionalities. Examples show that, contrary to standard deterministic or geostatistical modelling techniques (which tend to be difficult to control) the use of new 3D tools allows the geologist to interactively modify geological surfaces (including faults) or volumetric properties. Thus, the sensitivity of various economic parameters (oil in place, connected volumes, reserves) to major geological uncertainties can be evaluated. It is argued that future breakthroughs in geological modelling techniques are likely to happen in the development of interactive approaches rather than in the research of new mathematical algorithms.
3D Babcock-Leighton Solar Dynamo Models
NASA Astrophysics Data System (ADS)
Miesch, Mark S.; Hazra, Gopal; Karak, Bidya Binay; Teweldebirhan, Kinfe; Upton, Lisa
2016-05-01
We present results from the new STABLE (Surface flux Transport and Babcock Leighton) Dynamo Model. STABLE is a 3D Babcock-Leighton/Flux Transport dynamo model in which the source of poloidal field is the explicit emergence, distortion, and dispersal of bipolar magnetic regions (BMRs). In this talk I will discuss initial results with axisymmetric flow fields, focusing on the operation of the model, the general features of the cyclic solutions, and the challenge of achieving supercritical dynamo solutions using only the Babcock-Leighton source term. Then I will present dynamo simulations that include 3D convective flow fields based on the observed velocity power spectrum inferred from photospheric Dopplergrams. I'll use these simulations to assess how the explicit transport and amplification of fields by surface convection influences the operation of the dynamo. I will also discuss the role of surface magnetic fields in regulating the subsurface toroidal flux budget.
NASA Astrophysics Data System (ADS)
Brocher, T. M.; Frankel, A. D.; Oppenheimer, D. H.; Fletcher, J. B.; Luetgert, J. H.
2010-12-01
USGS 3D seismic velocity model. For the past year, long period motions have also been recorded by the NetQuakes project which is installing a large number (hundreds) of small, relatively inexpensive seismographs in 1 to 2 story homes and businesses. The instruments have an 18-bit resolution recorder and ± 3 g internal tri-axial MEMS accelerometers. NetQuakes instruments have successfully recorded long period motions from the January 9, 2010 M6.5 quake in the Gorda Plate near the Mendocino Triple Junction, the largest earthquake in northern California since June 17, 2005.
Regional geothermal 3D modelling in Denmark
NASA Astrophysics Data System (ADS)
Poulsen, S. E.; Balling, N.; Bording, T. S.; Nielsen, S. B.
2012-04-01
In the pursuit of sustainable and low carbon emission energy sources, increased global attention has been given to the exploration and exploitation of geothermal resources within recent decades. In 2009 a national multi-disciplinary geothermal research project was established. As a significant part of this project, 3D temperature modelling is to be carried out, with special emphasis on temperatures of potential geothermal reservoirs in the Danish area. The Danish subsurface encompasses low enthalpy geothermal reservoirs of mainly Triassic and Jurassic age. Geothermal plants at Amager (Copenhagen) and Thisted (Northern Jutland) have the capacity of supplying the district heating network with up to 14 MW and 7 MW, respectively, by withdrawing warm pore water from the Gassum (Lower Jurassic/Upper Triassic) and Bunter (Lower Triassic) sandstone reservoirs, respectively. Explorative studies of the subsurface temperature regime typically are based on a combination of observations and modelling. In this study, the open-source groundwater modelling code MODFLOW is modified to simulate the subsurface temperature distribution in three dimensions by taking advantage of the mathematical similarity between saturated groundwater flow (Darcy flow) and heat conduction. A numerical model of the subsurface geology in Denmark is built and parameterized from lithological information derived from joint interpretation of seismic surveys and borehole information. Boundary conditions are constructed from knowledge about the heat flow from the Earth's interior and the shallow ground temperature. Matrix thermal conductivities have been estimated from analysis of high-resolution temperature logs measured in deep wells and porosity-depth relations are included using interpreted main lithologies. The model takes into account the dependency of temperature and pressure on thermal conductivity. Moreover, a transient model based correction of the paleoclimatic thermal disturbance caused by the
NASA Astrophysics Data System (ADS)
Lestari, Titik; Nugraha, Andri Dian
2015-04-01
Southern Sumatra region has a high level of seismicity due to the influence of the subduction system, Sumatra fault, Mentawai fault and stretching zone activities. The seismic activities of Southern Sumatra region are recorded by Meteorological Climatological and Geophysical Agency (MCGA's) Seismograph network. In this study, we used earthquake data catalog compiled by MCGA for 3013 events from 10 seismic stations around Southern Sumatra region for time periods of April 2009 - April 2014 in order to invert for the 3-D seismic velocities structure (Vp, Vs, and Vp/Vs ratio). We applied double-difference seismic tomography method (tomoDD) to determine Vp, Vs and Vp/Vs ratio with hypocenter adjustment. For the inversion procedure, we started from the initial 1-D seismic velocity model of AK135 and constant Vp/Vs of 1.73. The synthetic travel time from source to receiver was calculated using ray pseudo-bending technique, while the main tomographic inversion was applied using LSQR method. The resolution model was evaluated using checkerboard test and Derivative Weigh Sum (DWS). Our preliminary results show low Vp and Vs anomalies region along Bukit Barisan which is may be associated with weak zone of Sumatran fault and migration of partial melted material. Low velocity anomalies at 30-50 km depth in the fore arc region may indicated the hydrous material circulation because the slab dehydration. We detected low seismic seismicity in the fore arc region that may be indicated as seismic gap. It is coincides contact zone of high and low velocity anomalies. And two large earthquakes (Jambi and Mentawai) also occurred at the contact of contrast velocity.
Lestari, Titik; Nugraha, Andri Dian
2015-04-24
Southern Sumatra region has a high level of seismicity due to the influence of the subduction system, Sumatra fault, Mentawai fault and stretching zone activities. The seismic activities of Southern Sumatra region are recorded by Meteorological Climatological and Geophysical Agency (MCGA’s) Seismograph network. In this study, we used earthquake data catalog compiled by MCGA for 3013 events from 10 seismic stations around Southern Sumatra region for time periods of April 2009 – April 2014 in order to invert for the 3-D seismic velocities structure (Vp, Vs, and Vp/Vs ratio). We applied double-difference seismic tomography method (tomoDD) to determine Vp, Vs and Vp/Vs ratio with hypocenter adjustment. For the inversion procedure, we started from the initial 1-D seismic velocity model of AK135 and constant Vp/Vs of 1.73. The synthetic travel time from source to receiver was calculated using ray pseudo-bending technique, while the main tomographic inversion was applied using LSQR method. The resolution model was evaluated using checkerboard test and Derivative Weigh Sum (DWS). Our preliminary results show low Vp and Vs anomalies region along Bukit Barisan which is may be associated with weak zone of Sumatran fault and migration of partial melted material. Low velocity anomalies at 30-50 km depth in the fore arc region may indicated the hydrous material circulation because the slab dehydration. We detected low seismic seismicity in the fore arc region that may be indicated as seismic gap. It is coincides contact zone of high and low velocity anomalies. And two large earthquakes (Jambi and Mentawai) also occurred at the contact of contrast velocity.
3-D Isotropic and Anisotropic S-velocity Structure in the North American Upper Mantle
NASA Astrophysics Data System (ADS)
Yuan, H.; Marone, F.; Romanowicz, B.; Abt, D.; Fischer, K.
2008-12-01
The tectonic diversity of the North American continent has led to a number of geological, tectonic and geodynamical models, many of which can be better tested with high resolution 3-d tomographic models of the isotropic and anisotropic mantle structure of the continent. In the framework of non-linear asymptotic coupling theory (NACT), we recently developed tools to invert long period seismic waveforms combined with SKS splitting data, for both isotropic and radial and azimuthal anisotropic S-wave velocity structure in the upper mantle at the continental scale (Marone et al., 2007; Marone and Romanowicz, 2007). Striking differences in both isotropic and anisotropic velocity structure were observed: beneath the high velocity stable cratonic region a distinct two-layer anisotropic domain is present, with the bottom layer fast axis direction aligned with the absolute plate motion, and a shallower lithospheric layer with north pointing fast axis most likely showing records of past tectonic history; under the active western US the direction of tomographically inferred anisotropy is stable with depth and compatible with the absolute plate motion direction. Here we present an updated model which includes nearly five more years of data, including data from newly operative USArray stations, and a somewhat more extended frequency band. Our new model confirms our previous results, and reveals greater yet complex details of the anisotropic velocity structure beneath the western U.S.. We also show initial results of incorporating constraints on the depth to the lithosphere-asthenosphere boundary (LAB) using teleseismic receiver functions. We discuss the different anisotropic domains resolved both laterally and in depth, in the context of tectonic history of the north American continent.
The crustal and mantle velocity structure in central Asia from 3D traveltime tomography
NASA Astrophysics Data System (ADS)
Sun, Y.; Martin, R. V.; Toksoz, M. N.; Pei, S.
2010-12-01
The lithospheric structure in central Asia features large blocks such as the Indian plate, the Afghan block, the Turan plate, and the Tarim block. This geologically and tectonically complicated area is also one of the most seismically active regions in the world. We developed P- and S- wave velocity structures of the central Asia in the crust using the traveltime data from Kyrgyzstan, Tajikistan, Kazakhstan, and Uzbek. We chose the events and stations between 32N65E and 45N85E and focused on the areas of Pamir and western Tianshan. In this data set, there are more than 6000 P and S arrivals received at 80 stations from about 300 events. The double difference tomography is applied to relocate events and to invert for seismic structures simultaneously. Our results provide accurate locations of earthquakes and high resolution crustal structure in this region. To extend the model deeper into the mantle through the upper mantle transition zone, ISC/EHB data for P and PP phases are combined with the ABCE data. To counteract the “smearing effect,” the crust and upper mantle velocity structure, derived from regional travel-times, is used. An adaptive grid method based on ray density is used in the inversion. A P-wave velocity model extending down to a depth of 2000 km is obtained. regional-teleseismic tomography provides a high-resolution, 3-D P-wave velocity model for the crust, upper mantle, and the transition zone. The crustal models correlate well with geologic and tectonic features. The upper mantle tomograms show the images of Tian Shan. The slab geometry is quite complex, reflecting the history of the changes in the plate motions and collision processes. Vp/Vs tomography was also determined in the study region, and an attenuation tomography was obtained as well.
Simulation of 3D Global Wave Propagation Through Geodynamic Models
NASA Astrophysics Data System (ADS)
Schuberth, B.; Piazzoni, A.; Bunge, H.; Igel, H.; Steinle-Neumann, G.
2005-12-01
This project aims at a better understanding of the forward problem of global 3D wave propagation. We use the spectral element program "SPECFEM3D" (Komatitsch and Tromp, 2002a,b) with varying input models of seismic velocities derived from mantle convection simulations (Bunge et al., 2002). The purpose of this approach is to obtain seismic velocity models independently from seismological studies. In this way one can test the effects of varying parameters of the mantle convection models on the seismic wave field. In order to obtain the seismic velocities from the temperature field of the geodynamical simulations we follow a mineral physics approach. Assuming a certain mantle composition (e.g. pyrolite with CMASF composition) we compute the stable phases for each depth (i.e. pressure) and temperature by system Gibbs free energy minimization. Elastic moduli and density are calculated from the equations of state of the stable mineral phases. For this we use a mineral physics database derived from calorimetric experiments (enthalphy and entropy of formation, heat capacity) and EOS parameters.
3-D Crustal Velocity Structure Across the Vrancea Zone in Romania, Derived From Seismic Data
NASA Astrophysics Data System (ADS)
Landes, M.; Hauser, F.; Popa, M.
2002-12-01
The Vrancea zone in the south-eastern Carpathians is one of the most active seismic zones in Europe. In order to study the crustal and upper-mantle structure in this region, two seismic refraction/wide-angle reflection experiments were carried out in 1999 and 2001. The 1999 campaign comprised a 320 km long N-S profile and a 80 km long transverse profile (E-W). All shots were recorded simultaneously on both profiles. The profile conducted in 2001 extended in E-W direction from the Hungarian border across the Vrancea zone to the Black Sea. We present an application of a 3-D refraction and reflection tomography algorithm (Hole 1992, 1995), elaborating the crustal velocity and interface structure within a 115 x 235 km wide region around the Vrancea zone. In order to enhance the model resolution, first arrival data from local earthquakes of the CALIXTO-99 teleseismic project were also included. The results indicate a high-velocity structure beneath the northern part of the Vrancea zone extending from shallow levels to depths of about 11 km. This structure may be related to the Trotus and Capidava-Ovidiu faults, which converge to the north of it. The high-velocity region is surrounded by the lower velocity Focsani and Brasov basins. The sedimentary succession beneath the southern part of the model extends to 18 km depth, while in the north sediment thickness varies between 10 and 15 km. Further results of the interface modelling of prominent reflections show that the mid-crustal and Moho interfaces shallow northwards from 30 km to 22 km and from 42 km to 38 km, respectively. This correlates well with previous results of Hauser et al. (2001).
A hybrid experimental-numerical technique for determining 3D velocity fields from planar 2D PIV data
NASA Astrophysics Data System (ADS)
Eden, A.; Sigurdson, M.; Mezić, I.; Meinhart, C. D.
2016-09-01
Knowledge of 3D, three component velocity fields is central to the understanding and development of effective microfluidic devices for lab-on-chip mixing applications. In this paper we present a hybrid experimental-numerical method for the generation of 3D flow information from 2D particle image velocimetry (PIV) experimental data and finite element simulations of an alternating current electrothermal (ACET) micromixer. A numerical least-squares optimization algorithm is applied to a theory-based 3D multiphysics simulation in conjunction with 2D PIV data to generate an improved estimation of the steady state velocity field. This 3D velocity field can be used to assess mixing phenomena more accurately than would be possible through simulation alone. Our technique can also be used to estimate uncertain quantities in experimental situations by fitting the gathered field data to a simulated physical model. The optimization algorithm reduced the root-mean-squared difference between the experimental and simulated velocity fields in the target region by more than a factor of 4, resulting in an average error less than 12% of the average velocity magnitude.
3-D physical models of amitosis (cytokinesis).
Cheng, Kang; Zou, Changhua
2005-01-01
Based on Newton's laws, extended Coulomb's law and published biological data, we develop our 3-D physical models of natural and normal amitosis (cytokinesis), for prokaryotes (bacterial cells) in M phase. We propose following hypotheses: Chromosome rings exclusion: No normally and naturally replicated chromosome rings (RCR) can occupy the same prokaryote, a bacterial cell. The RCR produce spontaneous and strong electromagnetic fields (EMF), that can be alternated environmentally, in protoplasm and cortex. The EMF is approximately a repulsive quasi-static electric (slowly variant and mostly electric) field (EF). The EF forces between the RCR are strong enough, and orderly accumulate contractile proteins that divide the procaryotes in the cell cortex of division plane or directly split the cell compartment envelope longitudinally. The radial component of the EF forces could also make furrows or cleavages of procaryotes. The EF distribution controls the protoplasm partition and completes the amitosis (cytokinesis). After the cytokinesis, the spontaneous and strong EF disappear because the net charge accumulation becomes weak, in the protoplasm. The exclusion is because the two sets of informative objects (RCR) have identical DNA codes information and they are electro magnetically identical, therefore they repulse from each other. We also compare divisions among eukaryotes, prokaryotes, mitochondria and chloroplasts and propose our hypothesis: The principles of our models are applied to divisions of mitochondria and chloroplasts of eucaryotes too because these division mechanisms are closer than others in a view of physics. Though we develop our model using 1 division plane (i.e., 1 cell is divided into 2 cells) as an example, the principle of our model is applied to the cases with multiple division planes (i.e., 1 cell is divided into multiple cells) too. PMID:15533619
3D P and S Wave Velocity Structure and Tremor Locations in the Parkfield Region
NASA Astrophysics Data System (ADS)
Zeng, X.; Thurber, C. H.; Shelly, D. R.; Bennington, N. L.; Cochran, E. S.; Harrington, R. M.
2014-12-01
We have assembled a new dataset to refine the 3D seismic velocity model in the Parkfield region. The S arrivals from 184 earthquakes recorded by the Parkfield Experiment to Record MIcroseismicity and Tremor array (PERMIT) during 2010-2011 were picked by a new S wave picker, which is based on machine learning. 74 blasts have been assigned to four quarries, whose locations were identified with Google Earth. About 1000 P and S wave arrivals from these blasts at permanent seismic network were also incorporated. Low frequency earthquakes (LFEs) occurring within non-volcanic tremor (NVT) are valuable for improving the precision of NVT location and the seismic velocity model at greater depths. Based on previous work (Shelley and Hardebeck, 2010), waveforms of hundreds of LFEs in same family were stacked to improve signal qualify. In a previous study (McClement et al., 2013), stacked traces of more than 30 LFE families at the Parkfileld Array Seismic Observatory (PASO) have been picked. We expanded our work to include LFEs recorded by the PERMIT array. The time-frequency Phase Weight Stacking (tf-PWS) method was introduced to improve the stack quality, as direct stacking does not produce clear S-wave arrivals on the PERMIT stations. This technique uses the coherence of the instantaneous phase among the stacked signals to enhance the signal-to-noise ratio (SNR) of the stack. We found that it is extremely effective for picking LFE arrivals (Thurber et al., 2014). More than 500 P and about 1000 S arrivals from 58 LFE families were picked at the PERMIT and PASO arrays. Since the depths of LFEs are much deeper than earthquakes, we are able to extend model resolution to lower crustal depths. Both P and S wave velocity structure have been obtained with the tomoDD method. The result suggests that there is a low velocity zone (LVZ) in the lower crust and the location of the LVZ is consistent with the high conductivity zone beneath the southern segment of the Rinconada fault that
3D Models of Stellar Interactions
NASA Astrophysics Data System (ADS)
Mohamed, S.; Podsiadlowski, Ph.; Booth, R.; Maercker, M.; Ramstedt, S.; Vlemmings, W.; Harries, T.; Mackey, J.; Langer, N.; Corradi, R.
2014-04-01
Symbiotic binaries consist of a cool, evolved mass-losing giant and an accreting compact companion. As symbiotic nebulae show similar morphologies to those in planetary nebulae (so much so that it is often difficult to distinguish between the two), they are ideal laboratories for understanding the role a binary companion plays in shaping the circumstellar envelopes in these evolved systems. We will present 3D Smoothed Particle Hydrodynamics (SPH) models of interacting binaries, e.g. R Aquarii and Mira, and discuss the formation of spiral outflows, arcs, shells and equatorial density enhancements.We will also discuss the implications of the former for planetary nebulae, e.g. the Egg Nebula and Cat's Eye, and the latter for the formation of bipolar geometries, e.g. M2-9. We also investigate accretion and angular momentum evolution in symbiotic binaries which may be important to understand the formation of jets and more episodic mass-loss features we see in circumstellar envelopes and the orbital characteristics of binary central stars of planetary nebulae.
NASA Astrophysics Data System (ADS)
Chong, J.; Yuan, H.; French, S. W.; Romanowicz, B. A.; Ni, S.
2011-12-01
Southeast Asia as a special region in the world which is seismically active and is surrounded by active tectonic belts, such as the Himalaya collision zone, western Pacific subduction zones and the Tianshan- Baikal tectonic belt. Seismic anisotropic tomography can shade light on the complex crust and upper mantle dynamics of this region, which is the subject of much debate. In this study, we applied full waveform time domain tomography to image 3D isotropic and anisotropic upper mantle shear velocity structure of Southeast Asia. Three component waveforms of teleseismic and far regional events (15 degree ≤ Δ≤ 165 degree) with magnitude ranges from Mw6.0 to Mw7.0 are collected from 91 permanent and 438 temporary broadband seismic stations in SE Asia. Wavepackets of both fundamental and overtone modes, filtered between 60 and 400 sec, are selected automatically according to the similarity between data and synthetic waveforms (Panning & Romanowicz, 2006). Wavepackets corresponding to event-station paths that sample the region considered are weighted according to path redundancy and signal to noise ratio. Higher modes and fundamental mode wavepackets are weighted separately in order to enhance the contribution of higher modes which are more sensitive to deeper structure compared to the fundamental mode. Synthetic waveforms and broadband sensitivity kernels are computed using normal mode asymptotic coupling theory (NACT, Li & Romanowicz, 1995). As a starting model, we consider a global anisotropic upper mantle shear velocity model based on waveform inversion using the Spectral Element Method (Lekic & Romanowicz, 2011), updated for more realistic crustal thickness (French et al., 2011) as our starting model, we correct waveforms for the effects of 3D structure outside of the region, and invert them for perturbations in the 3D structure of the target region only. We start with waveform inversion down to 60sec and after several iterations, we include shorter period
Multi-view and 3D deformable part models.
Pepik, Bojan; Stark, Michael; Gehler, Peter; Schiele, Bernt
2015-11-01
As objects are inherently 3D, they have been modeled in 3D in the early days of computer vision. Due to the ambiguities arising from mapping 2D features to 3D models, 3D object representations have been neglected and 2D feature-based models are the predominant paradigm in object detection nowadays. While such models have achieved outstanding bounding box detection performance, they come with limited expressiveness, as they are clearly limited in their capability of reasoning about 3D shape or viewpoints. In this work, we bring the worlds of 3D and 2D object representations closer, by building an object detector which leverages the expressive power of 3D object representations while at the same time can be robustly matched to image evidence. To that end, we gradually extend the successful deformable part model [1] to include viewpoint information and part-level 3D geometry information, resulting in several different models with different level of expressiveness. We end up with a 3D object model, consisting of multiple object parts represented in 3D and a continuous appearance model. We experimentally verify that our models, while providing richer object hypotheses than the 2D object models, provide consistently better joint object localization and viewpoint estimation than the state-of-the-art multi-view and 3D object detectors on various benchmarks (KITTI [2] , 3D object classes [3] , Pascal3D+ [4] , Pascal VOC 2007 [5] , EPFL multi-view cars[6] ). PMID:26440264
3D-GNOME: an integrated web service for structural modeling of the 3D genome
Szalaj, Przemyslaw; Michalski, Paul J.; Wróblewski, Przemysław; Tang, Zhonghui; Kadlof, Michal; Mazzocco, Giovanni; Ruan, Yijun; Plewczynski, Dariusz
2016-01-01
Recent advances in high-throughput chromosome conformation capture (3C) technology, such as Hi-C and ChIA-PET, have demonstrated the importance of 3D genome organization in development, cell differentiation and transcriptional regulation. There is now a widespread need for computational tools to generate and analyze 3D structural models from 3C data. Here we introduce our 3D GeNOme Modeling Engine (3D-GNOME), a web service which generates 3D structures from 3C data and provides tools to visually inspect and annotate the resulting structures, in addition to a variety of statistical plots and heatmaps which characterize the selected genomic region. Users submit a bedpe (paired-end BED format) file containing the locations and strengths of long range contact points, and 3D-GNOME simulates the structure and provides a convenient user interface for further analysis. Alternatively, a user may generate structures using published ChIA-PET data for the GM12878 cell line by simply specifying a genomic region of interest. 3D-GNOME is freely available at http://3dgnome.cent.uw.edu.pl/. PMID:27185892
3D-GNOME: an integrated web service for structural modeling of the 3D genome.
Szalaj, Przemyslaw; Michalski, Paul J; Wróblewski, Przemysław; Tang, Zhonghui; Kadlof, Michal; Mazzocco, Giovanni; Ruan, Yijun; Plewczynski, Dariusz
2016-07-01
Recent advances in high-throughput chromosome conformation capture (3C) technology, such as Hi-C and ChIA-PET, have demonstrated the importance of 3D genome organization in development, cell differentiation and transcriptional regulation. There is now a widespread need for computational tools to generate and analyze 3D structural models from 3C data. Here we introduce our 3D GeNOme Modeling Engine (3D-GNOME), a web service which generates 3D structures from 3C data and provides tools to visually inspect and annotate the resulting structures, in addition to a variety of statistical plots and heatmaps which characterize the selected genomic region. Users submit a bedpe (paired-end BED format) file containing the locations and strengths of long range contact points, and 3D-GNOME simulates the structure and provides a convenient user interface for further analysis. Alternatively, a user may generate structures using published ChIA-PET data for the GM12878 cell line by simply specifying a genomic region of interest. 3D-GNOME is freely available at http://3dgnome.cent.uw.edu.pl/. PMID:27185892
James E. Fisher; Cliff B. Davis; Walter L. Weaver
2005-06-01
A compressor model has been implemented in the RELAP5-3D© code. The model is similar to that of the existing pump model, and performs the same function on a gas as the pump performs on a single-phase or two-phase fluid. The compressor component consists of an inlet junction and a control volume, and optionally, an outlet junction. This feature permits cascading compressor components in series. The equations describing the physics of the compressor are derived from first principles. These equations are used to obtain the head, the torque, and the energy dissipation. Compressor performance is specified using a map, specific to the design of the machine, in terms of the ratio of outlet-to-inlet total (or stagnation) pressure and adiabatic efficiency as functions of rotational velocity and flow rate. The input quantities are specified in terms of dimensionless variables, which are corrected to stagnation density and stagnation sound speed. A small correction was formulated for the input of efficiency to account for the error introduced by assumption of constant density when integrating the momentum equation. Comparison of the results of steady-state operation of the compressor model to those of the MIT design calculation showed excellent agreement for both pressure ratio and power.
3D P-Wave Velocity Structure of the Deep Galicia Rifted Margin
NASA Astrophysics Data System (ADS)
Bayrakci, Gaye; Minshull, Timothy; Davy, Richard; Sawyer, Dale; Klaeschen, Dirk; Papenberg, Cord; Reston, Timothy; Shillington, Donna; Ranero, Cesar
2015-04-01
The combined wide-angle reflection-refraction and multi-channel seismic (MCS) experiment, Galicia 3D, was carried out in 2013 at the Galicia rifted margin in the northeast Atlantic Ocean, west of Spain. The main geological features within the 64 by 20 km (1280 km²) 3D box investigated by the survey are the peridotite ridge (PR), the fault bounded, rotated basement blocks and the S reflector, which has been interpreted to be a low angle detachment fault. 44 short period four-component ocean bottom seismometers and 28 ocean bottom hydrophones were deployed in the 3D box. 3D MCS profiles sampling the whole box were acquired with two airgun arrays of 3300 cu.in. fired alternately every 37.5 m. We present the results from 3D first-arrival time tomography that constrains the P-wave velocity in the 3D box, for the entire depth sampled by reflection data. Results are validated by synthetic tests and by the comparison with Galicia 3D MCS lines. The main outcomes are as follows: 1- The 3.5 km/s iso-velocity contour mimics the top of the acoustic basement observed on MCS profiles. Block bounding faults are imaged as velocity contrasts and basement blocks exhibit 3D topographic variations. 2- On the southern profiles, the top of the PR rises up to 5.5 km depth whereas, 20 km northward, its basement expression (at 6.5 km depth) nearly disappears. 3- The 6.5 km/s iso-velocity contour matches the topography of the S reflector where the latter is visible on MCS profiles. Within a depth interval of 0.6 km (in average), velocities beneath the S reflector increase from 6.5 km/s to 7 km/s, which would correspond to a decrease in the degree of serpentinization from ~45 % to ~30 % if these velocity variations are caused solely by variations in hydration. At the intersections between the block bounding normal faults and the S reflector, this decrease happens over a larger depth interval (> 1 km), suggesting that faults act as conduit for the water flow in the upper mantle.
3D fast wavelet network model-assisted 3D face recognition
NASA Astrophysics Data System (ADS)
Said, Salwa; Jemai, Olfa; Zaied, Mourad; Ben Amar, Chokri
2015-12-01
In last years, the emergence of 3D shape in face recognition is due to its robustness to pose and illumination changes. These attractive benefits are not all the challenges to achieve satisfactory recognition rate. Other challenges such as facial expressions and computing time of matching algorithms remain to be explored. In this context, we propose our 3D face recognition approach using 3D wavelet networks. Our approach contains two stages: learning stage and recognition stage. For the training we propose a novel algorithm based on 3D fast wavelet transform. From 3D coordinates of the face (x,y,z), we proceed to voxelization to get a 3D volume which will be decomposed by 3D fast wavelet transform and modeled after that with a wavelet network, then their associated weights are considered as vector features to represent each training face . For the recognition stage, an unknown identity face is projected on all the training WN to obtain a new vector features after every projection. A similarity score is computed between the old and the obtained vector features. To show the efficiency of our approach, experimental results were performed on all the FRGC v.2 benchmark.
A 3-D shape model of Interamnia
NASA Astrophysics Data System (ADS)
Sato, Isao
2015-08-01
A 3-D shape model of the sixth largest of the main belt asteroids, (704) Interamnia, is presented. The model is reproduced from its two stellar occultation observations and six lightcurves between 1969 and 2011. The first stellar occultation was the occultation of TYC 234500183 on 1996 December 17 observed from 13 sites in the USA. An elliptical cross section of (344.6±9.6km)×(306.2±9.1km), for position angle P=73.4±12.5 was fitted. The lightcurve around the occultation shows that the peak-to-peak amplitude was 0.04 mag. and the occultation phase was just before the minimum. The second stellar occultation was the occultation of HIP 036189 on 2003 March 23 observed from 39 sites in Japan and Hawaii. An elliptical cross section of (349.8±0.9km)×(303.7±1.7km), for position angle P=86.0±1.1 was fitted. A companion of 8.5 mag. of the occulted star was discovered whose separation is 12±2 mas (milli-arcseconds), P=148±11 . A combined analysis of rotational lightcurves and occultation chords can return more information than can be obtained with either technique alone. From follow-up photometric observations of the asteroid between 2003 and 2011, its rotation period is determined to be 8.728967167±0.00000007 hours, which is accurate enough to fix the rotation phases at other occultation events. The derived north pole is λ2000=259±8, β2000=-50±5 (retrograde rotation); the lengths of the three principal axes are 2a=361.8±2.8km, 2b=324.4±5.0km, 2c=297.3±3.5km, and the mean diameter is D=326.8±3.0km. Supposing the mass of Interamnia as (3.5±0.9)×10-11 solar masses, the density is then ρ=3.8±1.0 g cm-3.
Northern California Seismic Attenuation: 3-D Qp and Qs models
NASA Astrophysics Data System (ADS)
Eberhart-Phillips, D. M.
2015-12-01
The northern California crust exhibits a wide range of rock types and deformation processes which produce pronounced heterogeneity in regional attenuation. Using local earthquakes, 3-D Qp and Qs crustal models have been obtained for this region which includes the San Andreas fault system, the Central Valley, the Sierra Nevada batholith, and the Mendocino subduction volcanic system. Path attenuation t* values were determined from P and S spectra of 959 spatially distributed earthquakes, magnitude 2.5-6.0 from 2005-2014, using 1254 stations from NCEDC networks and IRIS Mendocino and Sierra Nevada temporary arrays. The t* data were used in Q inversions, using existing hypocenters and 3-D velocity models, with basic 10-km node spacing. The uneven data coverage was accounted for with linking of nodes into larger areas in order to provide useful Q images across the 3-D volume. The results at shallow depth (< 2 km) show very low Q in the Sacramento Delta, the Eureka area, and parts of the Bay Area. In the brittle crust, fault zones that have high seismicity exhibit low Q. In the lower crust, low Q is observed along fault zones that have large cumulative displacement and have experienced grain size reduction. Underlying active volcanic areas, low Q features are apparent below 20-km depth. Moderately high Q is associated with igneous rocks of the Sierra Nevada and Salinian block, while the Franciscan subduction complex shows moderately low Q. The most prominent high Q feature is related to the Great Valley Ophiolite.
Effect of postural changes on 3D joint angular velocity during starting block phase.
Slawinski, Jean; Dumas, Raphaël; Cheze, Laurence; Ontanon, Guy; Miller, Christian; Mazure-Bonnefoy, Alice
2013-01-01
Few studies have focused on the effect of posture during sprint start. The aim of this study was to measure the effect of the modification of horizontal distance between the blocks during sprint start on three dimensional (3D) joint angular velocity. Nine trained sprinters started using three different starting positions (bunched, medium and elongated). They were equipped with 63 passive reflective markers, and an opto-electronic Motion Analysis system was used to collect the 3D marker trajectories. During the pushing phase on the blocks, norm of the joint angular velocity (NJAV), 3D Euler angular velocity (EAV) and pushing time on the blocks were calculated. The results demonstrated that the decrease of the block spacing induces an opposite effect on the angular velocity of joints of the lower and the upper limbs. The NJAV of the upper limbs is greater in the bunched start, whereas the NJAV of the lower limbs is smaller. The modifications of NJAV were due to a combination of the movement of the joints in the different degrees of freedom. The medium start seems to be the best compromise because it leads, in a short pushing time, to a combination of optimal joint velocities for upper and lower segments. PMID:23062070
NASA Astrophysics Data System (ADS)
Przyborska, Anna; Kosecki, Szymon; Jakacki, Jaromir
2014-05-01
Kongsfjorden is a West Svalbard fjord with a surface area of about 210 km2. It is obvious that the depths of the outer and central basins are influenced by the open sea, under influence of West Spitsbergen Current (WSC), which curry out warm Atlantic water and cold East Spitsbergen Current, while the shallower, inner basin has a large glacial outflow and its maximum depths do not exceed 100 m. Freshwater stored in Spitsbergen glaciers have strong influence on local hydrology and physical fjord conditions. Both, local and shelf conditions have impact on state of the fjord. External forces like tides, velocities at the boundary and atmospheric forces together with sources of cold and dens fresh water in the fjords will give reliable representation of physical conditions in Kongsfjorden. Modeling could help to solve this problem and we have hope that we find answer which one is the most important for local conditions in fjord. Calculations of balances between cold fresh water and warm and salt will provide additional information that could help to answer the main question of the GAME (Growing of the Arctic Marine Ecosystem) project - what is the reaction of physically controlled Arctic marine ecosystem to temperature rise.
NASA Astrophysics Data System (ADS)
Fadel, I.; van der Meijde, M.; Kerle, N.; Lauritsen, N.
2015-03-01
Non-uniqueness of satellite gravity interpretation has traditionally been reduced by using a priori information from seismic tomography models. This reduction in the non-uniqueness has been based on velocity-density conversion formulas or user interpretation of the 3D subsurface structures (objects) based on the seismic tomography models and then forward modelling these objects. However, this form of object-based approach has been done without a standardized methodology on how to extract the subsurface structures from the 3D models. In this research, a 3D object-oriented image analysis (3D OOA) approach was implemented to extract the 3D subsurface structures from geophysical data. The approach was applied on a 3D shear wave seismic tomography model of the central part of the East African Rift System. Subsequently, the extracted 3D objects from the tomography model were reconstructed in the 3D interactive modelling environment IGMAS+, and their density contrast values were calculated using an object-based inversion technique to calculate the forward signal of the objects and compare it with the measured satellite gravity. Thus, a new object-based approach was implemented to interpret and extract the 3D subsurface objects from 3D geophysical data. We also introduce a new approach to constrain the interpretation of the satellite gravity measurements that can be applied using any 3D geophysical model.
3D Modeling Techniques for Print and Digital Media
NASA Astrophysics Data System (ADS)
Stephens, Megan Ashley
In developing my thesis, I looked to gain skills using ZBrush to create 3D models, 3D scanning, and 3D printing. The models created compared the hearts of several vertebrates and were intended for students attending Comparative Vertebrate Anatomy. I used several resources to create a model of the human heart and was able to work from life while creating heart models from other vertebrates. I successfully learned ZBrush and 3D scanning, and successfully printed 3D heart models. ZBrush allowed me to create several intricate models for use in both animation and print media. The 3D scanning technique did not fit my needs for the project, but may be of use for later projects. I was able to 3D print using two different techniques as well.
Numerical modeling of Tibetan Plateau formation: Thin-sheet versus fully 3D models
NASA Astrophysics Data System (ADS)
Lechmann, S. M.; Schmalholz, S. M.; Kaus, B. J. P.
2009-04-01
Knowledge about the tectonic evolution of the Tibetan Plateau is still incomplete and many open questions remain concerning the deformation style of the crustal thickening, causing the abnormally high elevation of the Tibetan Plateau. Different models have been suggested explaining the crustal thickening by (1) homogeneous, continuous deformation using thin-sheet models, (2) discrete movement along thrusts developing crustal wedges and (3) lateral crustal flow due to pressure gradients resulting from topography. Most existing models are not fully three-dimensional (3D) models (e.g. thin-sheet models) and assume a certain deformation style a priori, which makes it difficult to judge the applicability of such constrained models to the formation of the Tibetan Plateau. We present a comparison of deformation styles during continent indentation resulting from a fully 3D numerical model and a thin-sheet model. The rheology for both models is power-law. The 3D model consists of four layers representing a simplified lithosphere: strong upper crust, weak lower crust, strong upper mantle and weak lower mantle. From the effective viscosity distribution of the 3D model a vertically averaged effective viscosity is calculated and used for the thin-sheet model to make direct comparisons between the two models. Simulating indentation is achieved by assigning free slip at one lateral side of the model, and fixing two other sides. The boundary at which indentation is taking place, exhibits a tripartite velocity profile: Next to the free slip side a section with constant horizontal velocity is applied. The velocity then gradually decreases towards zero, applying a cosine-function. The last section of the indenting boundary next to the fixed side is also fixed. The 3D model additionally exhibits a free surface and a bottom boundary allowing free slip. The 3D code employs the finite element method with a mixed velocity-pressure formulation to simulate incompressible flow. A Lagrangian
3D visualization of endocardial peak velocities during systole and diastole
NASA Astrophysics Data System (ADS)
Eusemann, Christian D.; Ritman, Erik L.; Robb, Richard A.
2002-04-01
Quantitative assessment of regional heart motion has the potential to provide diagnostic data for assessment of cardiac malfunction. Local heart motion may be obtained with various medical imaging scanners, so the goal is to provide an imaging modality-independent display/analysis technique. In this study, 3D reconstructions of a canine heart before and after infarction were obtained from the Dynamic Spatial Reconstructor (DSR) at 15 time points throughout one cardiac cycle. Deformable models of each time point were created. Through this process regional excursions and velocities in the mesh can be assigned to represent a piece of endocardium, which can be calculated for each time-point interval. These calculations are based on the distance change between a single vertex of the mesh and the model centerline from LV apex to aortic/mitral valve separation. This allows computation of color maps corresponding to regional values of contraction or dilation motion of the endocardium relative to the LV long axis (centerline) during systole and/or diastole. These color maps can be illustrated through model animations and multi view static images. Using functional parametric mappings of disturbances in regional contractility and relaxation facilitates appreciation of the effect of altered structure-to-function relationships in the myocardium.
3D modeling of metallic grain growth
George, D.; Carlson, N.; Gammel, J.T.; Kuprat, A.
1999-06-01
This paper will describe simulating metallic grain growth using the Gradient Weighted Moving Finite Elements code, GRAIN3D. The authors also describe the set of mesh topology change operations developed to respond to changes in the physical topology such as the collapse of grains and to maintain uniform calculational mesh quality. Validation of the method is demonstrated by comparison to analytic calculations. The authors present results of multigrain simulations where grain boundaries evolve by mean curvature motion and include results which incorporate grain boundary orientation dependence.
NASA Astrophysics Data System (ADS)
Bayrakci, Gaye; Minshull, Timothy A.; Davy, Richard G.; Karplus, Marianne S.; Kaeschen, Dirk; Papenberg, Cord; Krabbenhoeft, Anne; Sawyer, Dale; Reston, Timothy J.; Shillington, Donna J.; Ranero, César R.
2014-05-01
Galicia 3D, a reflection-refraction and long offset seismic experiment was carried out from May through September 2013, at the Galicia rifted margin (in the northeast Atlantic Ocean, west of Spain) as a collaboration between US, UK, German and Spanish groups. The 3D multichannel seismic acquisition conducted by R/V Marcus Langseth covered a 64 km by 20 km (1280 km2) zone where the main geological features are the Peridotite Ridge (PR), composed of serpentinized peridotite and thought be upper mantle exhumed to the seafloor during rifting, and the S reflector which has been interpreted to be a low angle detachment fault overlain by fault bounded, rotated, continental crustal blocks. In the 3D box, two airgun arrays of 3300 cu.in. were fired alternately (in flip-flop configuration) every 37.5 m. All shots are recorded by 44 short period four component ocean bottom seismometers (OBS) and 26 ocean bottom hydrophones (OBH) deployed and recovered by R/V Poseidon, as well as four 6 km hydrophone streamers with 12.5 m channel spacing towed by R/V Marcus Langseth. We present the preliminary results of the first arrival time tomography study which is carried out with a subset of the wide-angle dataset, in order to generate a 3D P-wave velocity volume for the entire depth sampled by the reflection data. After the relocation of OBSs and OBHs, an automatic first-arrival time picking approach is applied to a subset of the dataset, which comprises more than 5.5 million source-receiver pairs. Then, the first-arrival times are checked visually, in 3-dimensions. The a priori model used for the first-arrival time tomography is built up using information from previous seismic surveys carried out at the Galicia margin (e.g. ISE, 1997). The FAST algorithm of Zelt and Barton (1998) is used for the first-arrival time inversion. The 3D P-wave velocity volume can be used in interpreting the reflection dataset, as a starting point for migration, to quantify the thinning of the crustal layers
Towards Forward Modeling of 3D Heterogeneity in D" region
NASA Astrophysics Data System (ADS)
To, A.; Capdeville, Y.; Romanowicz, B.
2002-12-01
The presence of strong lateral heterogeneity in D" is now well documented. While tomographic modeling provides constraints on the large scale patterns, strong variations on shorter scales are best addressed by forward modeling. Appropriate tools are needed for forward modeling that will handle strong 3D heterogeneity, at relatively short periods and including diffracted waves. We use a coupled mode/SEM (Spectral Element Method) to compute synthetic seismograms in 3D models of the D" layer down to 1/12s. This coupled method (Capdeville, 2001) affords faster computations than SEM in cases where heterogeneity can be restricted to a specific layer. We compare them with observed waveforms for several events in the Western Pacific. Observed and synthetic travel time trends are very consistent, although in most cases the observed residuals are significantly larger. Waveform amplitudes are less consistent. In order to understand the origin of the amplitude difference, we test the effect of 3D heterogeneity on Sdiff phase. In particular, the results show opposite trends in the amplitude of Sdiff due to heterogeneity located near the CMB or well above it. This provides constraints on the location of the causative velocity heterogeneity. Because the forward modeling approach requires many iterations, the coupled mode/SEM approach is still computationally intensive. It is more efficient to use a less accurate traditional approach to first get closer to a final model, and only then use coupled mode/SEM to refine the model. Ray theory is the most expedient way to calculate travel times. However, it is an infinite frequency approximation and not appropriate to handle diffracting waves. We show that ray theory predicts larger travel time anomaly for Sdiff phase than the one obtained by coupled mode/SEM. Although it is based on a weak heterogeneity assumption, Non-linear Asymptotic Coupling Theory(NACT) (Li and Romanowicz, 1995) helps to overcome this difficulty. It can handle
3D Face Modeling Using the Multi-Deformable Method
Hwang, Jinkyu; Yu, Sunjin; Kim, Joongrock; Lee, Sangyoun
2012-01-01
In this paper, we focus on the problem of the accuracy performance of 3D face modeling techniques using corresponding features in multiple views, which is quite sensitive to feature extraction errors. To solve the problem, we adopt a statistical model-based 3D face modeling approach in a mirror system consisting of two mirrors and a camera. The overall procedure of our 3D facial modeling method has two primary steps: 3D facial shape estimation using a multiple 3D face deformable model and texture mapping using seamless cloning that is a type of gradient-domain blending. To evaluate our method's performance, we generate 3D faces of 30 individuals and then carry out two tests: accuracy test and robustness test. Our method shows not only highly accurate 3D face shape results when compared with the ground truth, but also robustness to feature extraction errors. Moreover, 3D face rendering results intuitively show that our method is more robust to feature extraction errors than other 3D face modeling methods. An additional contribution of our method is that a wide range of face textures can be acquired by the mirror system. By using this texture map, we generate realistic 3D face for individuals at the end of the paper. PMID:23201976
Life in 3D is never flat: 3D models to optimise drug delivery.
Fitzgerald, Kathleen A; Malhotra, Meenakshi; Curtin, Caroline M; O' Brien, Fergal J; O' Driscoll, Caitriona M
2015-10-10
The development of safe, effective and patient-acceptable drug products is an expensive and lengthy process and the risk of failure at different stages of the development life-cycle is high. Improved biopharmaceutical tools which are robust, easy to use and accurately predict the in vivo response are urgently required to help address these issues. In this review the advantages and challenges of in vitro 3D versus 2D cell culture models will be discussed in terms of evaluating new drug products at the pre-clinical development stage. Examples of models with a 3D architecture including scaffolds, cell-derived matrices, multicellular spheroids and biochips will be described. The ability to simulate the microenvironment of tumours and vital organs including the liver, kidney, heart and intestine which have major impact on drug absorption, distribution, metabolism and toxicity will be evaluated. Examples of the application of 3D models including a role in formulation development, pharmacokinetic profiling and toxicity testing will be critically assessed. Although utilisation of 3D cell culture models in the field of drug delivery is still in its infancy, the area is attracting high levels of interest and is likely to become a significant in vitro tool to assist in drug product development thus reducing the requirement for unnecessary animal studies. PMID:26220617
Assessing the RELAPS-3D Heat Conduction Enclosure Model
McCann, Larry D.
2008-09-30
Three heat conduction problems that have exact solutions are modeled with RELAP5-3D using the conduction enclosure model. These comparisons are designed to be used in the RELAP5-3D development assessment scheduled to be completed in 2009. It is shown that with proper input choices and adequate model detail the exact solutions can be matched. In addition, this analysis identified an error and the required correction in the cylindrical and spherical heat conductor models in RELAP5-3D which will be corrected in a future version of RELAP5-3D.
3D Seismic Velocity Structure in the Rupture Area of the 2010 Maule Mw=8.8 Earthquake
NASA Astrophysics Data System (ADS)
Hicks, S. P.; Rietbrock, A.; Ryder, I. M.; Nippress, S.; Haberland, C. A.
2011-12-01
The 2010 Mw=8.8 Maule, Chile earthquake is one of the largest subduction zone earthquakes ever recorded. Up to now numerous co-seismic and some post-seismic slip models have been published based entirely on seismological, geodetic, or tsunami run-up heights, or combinations of these data. Most of these models use a simplified megathrust geometry derived mainly from global earthquake catalogues, and also simplified models of seismic parameters (e.g. shear modulus). By using arrival times for a vast number of aftershocks that have been recorded on a temporary seismic array, we present a new model for the slab geometry based on earthquake locations together with a new 3D seismic velocity model of the region, for both vp and vp/vs. We analyzed 3552 aftershocks that occurred between 18 March and 24 May 2011, recorded by the International Maule Aftershock Dataset (IMAD) seismic network. Event selection from a catalogue of automatically-determined events was based on 20 or more arrival times, from which at least 10 are S-wave observations. In total over 170,000 arrival times (~125,000 and 45,000 P and S wave arrival times respectively) are used for the tomographic reconstructions. Initially, events were relocated in a 2D velocity model based on a previously published model for the southern end of the rupture area (Haberland et al., 2009). Afterwards a staggered inversion scheme is implemented, starting with a 2D inversion followed by a coarse 3D and a subsequent fine 3D inversion. Based on our preliminary inversions we conclude that aftershock seismicity is mainly concentrated between 20 and 35 km depth along the subduction interface. A second band of seismicity between 40 and 50 km depth is also observed. Low seismic velocities and an increased vp/vs ratio characterize the marine forearc. The obtained velocity model will be discussed.
3D Seismic Velocity Structure Around Philippine Sea Slab Subducting Beneath Kii Peninsula, Japan
NASA Astrophysics Data System (ADS)
Shibutani, T.; Imai, M.; Hirahara, K.; Nakao, S.
2013-12-01
Kii Peninsula is a part of the source area of Nankai Trough megaquakes and the region through which the strong seismic waves propagate to big cities in Kansai such as Osaka, Kyoto, Nara, Kobe, and so on. Moreover, the rupture starting point is thought to be possibly at off the peninsula. Therefore, it is important for simulations of the megaquakes and the strong motions to estimate accurately the configuration of the Philippine Sea slab and the seismic velocity structure around the slab and to investigate properties and conditions of the plate boundary surface. Deep low frequency events (DLFEs) are widely distributed from western Shikoku to central Tokai at 30 - 40 km depths on the plate boundary (Obara, 2002). Results from seismic tomography and receiver function analyses revealed that the oceanic crust of the Philippine Sea plate had a low velocity and a high Vp/Vs ratio (Hirose et al., 2007; Ueno et al., 2008). Hot springs with high 3He/4He ratios are found in an area between central Kinki and Kii Peninsula despite in the forearc region (Sano and Wakita, 1985). These phenomena suggest the process that H2O subducting with the oceanic crust dehydrates at the depths, causes the DLFEs, and moves to shallower depths. We carried out linear array seismic observations in the Kii Peninsula since 2004 in order to estimate the structure of the Philippine Sea slab and the surrounding area. We have performed receiver function analyses for four profile lines in the dipping direction of the slab and two lines in the perpendicular direction so far. We estimated three dimensional shapes of seismic velocity discontinuities such as the continental Moho, the upper surface of the oceanic crust and the oceanic Moho (Imai et al., 2013, this session). In addition, we performed seismic tomography with a velocity model embedded the discontinuities and observed travel times at stations in the linear arrays, and successfully estimated 3D seismic velocity structure around the Philippine Sea
3D scene modeling from multiple range views
NASA Astrophysics Data System (ADS)
Sequeira, Vitor; Goncalves, Joao G. M.; Ribeiro, M. Isabel
1995-09-01
This paper presents a new 3D scene analysis system that automatically reconstructs the 3D geometric model of real-world scenes from multiple range images acquired by a laser range finder on board of a mobile robot. The reconstruction is achieved through an integrated procedure including range data acquisition, geometrical feature extraction, registration, and integration of multiple views. Different descriptions of the final 3D scene model are obtained: a polygonal triangular mesh, a surface description in terms of planar and biquadratics surfaces, and a 3D boundary representation. Relevant experimental results from the complete 3D scene modeling are presented. Direct applications of this technique include 3D reconstruction and/or update of architectual or industrial plans into a CAD model, design verification of buildings, navigation of autonomous robots, and input to virtual reality systems.
Visualization of 3D Geological Models on Google Earth
NASA Astrophysics Data System (ADS)
Choi, Y.; Um, J.; Park, M.
2013-05-01
Google Earth combines satellite imagery, aerial photography, thematic maps and various data sets to make a three-dimensional (3D) interactive image of the world. Currently, Google Earth is a popular visualization tool in a variety of fields and plays an increasingly important role not only for private users in daily life, but also for scientists, practitioners, policymakers and stakeholders in research and application. In this study, a method to visualize 3D geological models on Google Earth is presented. COLLAborative Design Activity (COLLADA, an open standard XML schema for establishing interactive 3D applications) was used to represent different 3D geological models such as borehole, fence section, surface-based 3D volume and 3D grid by triangle meshes (a set of triangles connected by their common edges or corners). In addition, we designed Keyhole Markup Language (KML, the XML-based scripting language of Google Earth) codes to import the COLLADA files into the 3D render window of Google Earth. The method was applied to the Grosmont formation in Alberta, Canada. The application showed that the combination of COLLADA and KML enables Google Earth to effectively visualize 3D geological structures and properties.; Visualization of the (a) boreholes, (b) fence sections, (c) 3D volume model and (d) 3D grid model of Grossmont formation on Google Earth
Comparing a quasi-3D to a full 3D nearshore circulation model: SHORECIRC and ROMS
NASA Astrophysics Data System (ADS)
Haas, Kevin A.; Warner, John C.
Predictions of nearshore and surf zone processes are important for determining coastal circulation, impacts of storms, navigation, and recreational safety. Numerical modeling of these systems facilitates advancements in our understanding of coastal changes and can provide predictive capabilities for resource managers. There exists many nearshore coastal circulation models, however they are mostly limited or typically only applied as depth integrated models. SHORECIRC is an established surf zone circulation model that is quasi-3D to allow the effect of the variability in the vertical structure of the currents while maintaining the computational advantage of a 2DH model. Here we compare SHORECIRC to ROMS, a fully 3D ocean circulation model which now includes a three dimensional formulation for the wave-driven flows. We compare the models with three different test applications for: (i) spectral waves approaching a plane beach with an oblique angle of incidence; (ii) monochromatic waves driving longshore currents in a laboratory basin; and (iii) monochromatic waves on a barred beach with rip channels in a laboratory basin. Results identify that the models are very similar for the depth integrated flows and qualitatively consistent for the vertically varying components. The differences are primarily the result of the vertically varying radiation stress utilized by ROMS and the utilization of long wave theory for the radiation stress formulation in vertical varying momentum balance by SHORECIRC. The quasi-3D model is faster, however the applicability of the fully 3D model allows it to extend over a broader range of processes, temporal, and spatial scales.
Comparing a quasi-3D to a full 3D nearshore circulation model: SHORECIRC and ROMS
Haas, K.A.; Warner, J.C.
2009-01-01
Predictions of nearshore and surf zone processes are important for determining coastal circulation, impacts of storms, navigation, and recreational safety. Numerical modeling of these systems facilitates advancements in our understanding of coastal changes and can provide predictive capabilities for resource managers. There exists many nearshore coastal circulation models, however they are mostly limited or typically only applied as depth integrated models. SHORECIRC is an established surf zone circulation model that is quasi-3D to allow the effect of the variability in the vertical structure of the currents while maintaining the computational advantage of a 2DH model. Here we compare SHORECIRC to ROMS, a fully 3D ocean circulation model which now includes a three dimensional formulation for the wave-driven flows. We compare the models with three different test applications for: (i) spectral waves approaching a plane beach with an oblique angle of incidence; (ii) monochromatic waves driving longshore currents in a laboratory basin; and (iii) monochromatic waves on a barred beach with rip channels in a laboratory basin. Results identify that the models are very similar for the depth integrated flows and qualitatively consistent for the vertically varying components. The differences are primarily the result of the vertically varying radiation stress utilized by ROMS and the utilization of long wave theory for the radiation stress formulation in vertical varying momentum balance by SHORECIRC. The quasi-3D model is faster, however the applicability of the fully 3D model allows it to extend over a broader range of processes, temporal, and spatial scales. ?? 2008 Elsevier Ltd.
a Fast Method for Measuring the Similarity Between 3d Model and 3d Point Cloud
NASA Astrophysics Data System (ADS)
Zhang, Zongliang; Li, Jonathan; Li, Xin; Lin, Yangbin; Zhang, Shanxin; Wang, Cheng
2016-06-01
This paper proposes a fast method for measuring the partial Similarity between 3D Model and 3D point Cloud (SimMC). It is crucial to measure SimMC for many point cloud-related applications such as 3D object retrieval and inverse procedural modelling. In our proposed method, the surface area of model and the Distance from Model to point Cloud (DistMC) are exploited as measurements to calculate SimMC. Here, DistMC is defined as the weighted distance of the distances between points sampled from model and point cloud. Similarly, Distance from point Cloud to Model (DistCM) is defined as the average distance of the distances between points in point cloud and model. In order to reduce huge computational burdens brought by calculation of DistCM in some traditional methods, we define SimMC as the ratio of weighted surface area of model to DistMC. Compared to those traditional SimMC measuring methods that are only able to measure global similarity, our method is capable of measuring partial similarity by employing distance-weighted strategy. Moreover, our method is able to be faster than other partial similarity assessment methods. We demonstrate the superiority of our method both on synthetic data and laser scanning data.
A 3D Geometry Model Search Engine to Support Learning
ERIC Educational Resources Information Center
Tam, Gary K. L.; Lau, Rynson W. H.; Zhao, Jianmin
2009-01-01
Due to the popularity of 3D graphics in animation and games, usage of 3D geometry deformable models increases dramatically. Despite their growing importance, these models are difficult and time consuming to build. A distance learning system for the construction of these models could greatly facilitate students to learn and practice at different…
3-D Modeling of a Nearshore Dye Release
NASA Astrophysics Data System (ADS)
Maxwell, A. R.; Hibler, L. F.; Miller, L. M.
2006-12-01
The usage of computer modeling software in predicting the behavior of a plume discharged into deep water is well established. Nearfield plume spreading in coastal areas with complex bathymetry is less commonly studied; in addition to geometry, some of the difficulties of this environment include: tidal exchange, temperature, and salinity gradients. Although some researchers have applied complex hydrodynamic models to this problem, nearfield regions are typically modeled by calibration of an empirical or expert system model. In the present study, the 3D hydrodynamic model Delft3D-FLOW was used to predict the advective transport from a point release in Sequim Bay, Washington. A nested model approach was used, wherein a coarse model using a mesh extending to nearby tide gages (cell sizes up to 1 km) was run over several tidal cycles in order to provide boundary conditions to a smaller area. The nested mesh (cell sizes up to 30 m) was forced on two open boundaries using the water surface elevation derived from the coarse model. Initial experiments with the uncalibrated model were conducted in order to predict plume propagation based on the best available field data. Field experiments were subsequently carried out by releasing rhodamine dye into the bay at near-peak flood tidal current and near high slack tidal conditions. Surface and submerged releases were carried out from an anchored vessel. Concurrently collected data from the experiment include temperature, salinity, dye concentration, and hyperspectral imagery, collected from boats and aircraft. A REMUS autonomous underwater vehicle was used to measure current velocity and dye concentration at varying depths, as well as to acquire additional bathymetric information. Preliminary results indicate that the 3D hydrodynamic model offers a reasonable prediction of plume propagation speed and shape. A sensitivity analysis is underway to determine the significant factors in effectively using the model as a predictive tool
Energy flow in passive and active 3D cochlear model
Wang, Yanli; Steele, Charles; Puria, Sunil
2015-12-31
Energy flow in the cochlea is an important characteristic of the cochlear traveling wave, and many investigators, such as von Békésy and Lighthill, have discussed this phenomenon. Particularly after the discovery of the motility of the outer hair cells (OHCs), the nature of the power gain of the cochlea has been a fundamental research question. In the present work, direct three-dimensional (3D) calculations of the power on cross sections of the cochlea and on the basilar membrane are performed based on a box model of the mouse cochlea. The distributions of the fluid pressure and fluid velocity in the scala vestibuli are presented. The power output from the OHCs and the power loss due to fluid viscous damping are calculated along the length of the cochlea. This work provides a basis for theoretical calculations of the power gain of the OHCs from mechanical considerations.
Energy flow in passive and active 3D cochlear model
NASA Astrophysics Data System (ADS)
Wang, Yanli; Puria, Sunil; Steele, Charles
2015-12-01
Energy flow in the cochlea is an important characteristic of the cochlear traveling wave, and many investigators, such as von Békésy and Lighthill, have discussed this phenomenon. Particularly after the discovery of the motility of the outer hair cells (OHCs), the nature of the power gain of the cochlea has been a fundamental research question. In the present work, direct three-dimensional (3D) calculations of the power on cross sections of the cochlea and on the basilar membrane are performed based on a box model of the mouse cochlea. The distributions of the fluid pressure and fluid velocity in the scala vestibuli are presented. The power output from the OHCs and the power loss due to fluid viscous damping are calculated along the length of the cochlea. This work provides a basis for theoretical calculations of the power gain of the OHCs from mechanical considerations.
Evaluation of 3D-Jury on CASP7 models
Kaján, László; Rychlewski, Leszek
2007-01-01
Background 3D-Jury, the structure prediction consensus method publicly available in the Meta Server , was evaluated using models gathered in the 7th round of the Critical Assessment of Techniques for Protein Structure Prediction (CASP7). 3D-Jury is an automated expert process that generates protein structure meta-predictions from sets of models obtained from partner servers. Results The performance of 3D-Jury was analysed for three aspects. First, we examined the correlation between the 3D-Jury score and a model quality measure: the number of correctly predicted residues. The 3D-Jury score was shown to correlate significantly with the number of correctly predicted residues, the correlation is good enough to be used for prediction. 3D-Jury was also found to improve upon the competing servers' choice of the best structure model in most cases. The value of the 3D-Jury score as a generic reliability measure was also examined. We found that the 3D-Jury score separates bad models from good models better than the reliability score of the original server in 27 cases and falls short of it in only 5 cases out of a total of 38. We report the release of a new Meta Server feature: instant 3D-Jury scoring of uploaded user models. Conclusion The 3D-Jury score continues to be a good indicator of structural model quality. It also provides a generic reliability score, especially important for models that were not assigned such by the original server. Individual structure modellers can also benefit from the 3D-Jury scoring system by testing their models in the new instant scoring feature available in the Meta Server. PMID:17711571
3D velocity measurements in a premixed flame by tomographic PIV
NASA Astrophysics Data System (ADS)
Tokarev, M. P.; Sharaborin, D. K.; Lobasov, A. S.; Chikishev, L. M.; Dulin, V. M.; Markovich, D. M.
2015-06-01
Tomographic particle image velocimetry (PIV) has become a standard tool for 3D velocity measurements in non-reacting flows. However, the majority of the measurements in flows with combustion are limited to small resolved depth compared to the size of the field of view (typically 1 : 10). The limitations are associated with inhomogeneity of the volume illumination and the non-uniform flow seeding, the optical distortions and errors in the 3D calibration, and the unwanted flame luminosity. In the present work, the above constraints were overcome for the tomographic PIV experiment in a laminar axisymmetric premixed flame. The measurements were conducted for a 1 : 1 depth-to-size ratio using a system of eight CCD cameras and a 200 mJ pulsed laser. The results show that camera calibration based on the triangulation of the tracer particles in the non-reacting conditions provided reliable accuracy for the 3D image reconstruction in the flame. The modification of the tomographic reconstruction allowed a posteriori removal of unwanted bright objects, which were located outside of the region of interest but affected the reconstruction quality. This study reports on a novel experience for the instantaneous 3D velocimetry in laboratory-scale flames by using tomographic PIV.
3D density model of the Central Andes
NASA Astrophysics Data System (ADS)
Prezzi, Claudia B.; Götze, Hans-Jürgen; Schmidt, Sabine
2009-12-01
We developed a 3D density model of the continental crust, the subducted plate and the upper mantle of the Central Andes between 20-29°S and 74-61°W through the forward modelling of Bouguer anomaly. The goal of this contribution is to gain insight on the lithospheric structure integrating the available information (geophysical, geologic, petrologic, and geochemical) in a single model. The geometry of our model is defined and constrained by hypocentre location, reflection and refraction on and offshore seismic lines, travel time and attenuation tomography, receiver function analysis, magnetotelluric studies, thermal models and balanced structural cross-sections. The densities allocated to the different bodies are calculated considering petrologic and geochemical data and pressure and temperature conditions. The model consists of 31 parallel E-W vertical planes, where the continental crust comprises distinct bodies, which represent the different morphotectonic units of the Central Andes. We include a partial melting zone at midcrustal depths under the Altiplano-Puna (low-velocity zone) and consider the presence of a rheologically strong block beneath the Salar de Atacama basin, according to recent seismic studies. Contour maps of the depth of the continental Moho, the thickness of the lower crust and the depth to the bottom of the lithosphere below South America are produced. The possible percentage of partial melt in the Central Andes low-velocity zone is estimated. The residual anomaly is calculated by subtracting from the Bouguer anomaly the gravimetric effect of the modelled subducted slab and of the modelled Moho. Isostatic anomalies are calculated from regional and local isostatic Mohos calculated with and without internal loads, derived from our gravity model, which are then compared to the modelled continental Moho. This study contributes to a more detailed knowledge of the lithospheric structure of this region of the Andes and provides an integrated 3D
Computational modeling of RNA 3D structures and interactions.
Dawson, Wayne K; Bujnicki, Janusz M
2016-04-01
RNA molecules have key functions in cellular processes beyond being carriers of protein-coding information. These functions are often dependent on the ability to form complex three-dimensional (3D) structures. However, experimental determination of RNA 3D structures is difficult, which has prompted the development of computational methods for structure prediction from sequence. Recent progress in 3D structure modeling of RNA and emerging approaches for predicting RNA interactions with ions, ligands and proteins have been stimulated by successes in protein 3D structure modeling. PMID:26689764
3D P-Wave Velocity Structure of the Crust and Relocation of Earthquakes in 21 the Lushan Source Area
NASA Astrophysics Data System (ADS)
Yu, X.; Wang, X.; Zhang, W.
2014-12-01
The double difference seismic tomography method is applied to the absolute first arrival P wave arrival times and high quality relative P arrival times of the Lushan seismic sequence to determine the detailed crustal 3D P wave velocity structure and the hypocenter parameters in the Lushan seismic area. The results show that the Lushan mainshock locates at 30.28 N, 103.98 E, with the depth of 16.38 km. The leading edge of aftershock in the northeast of mainshock present a spade with a steep dip angle, the aftershocks' extended length is about 12 km. In the southwest of the Lushan mainshock, the leading edge of aftershock in low velocity zone slope gently, the aftershocks' extended length is about 23 km. The P wave velocity structure of the Lushan seismic area shows obviously lateral heterogeneity. The P wave velocity anomalies represent close relationship with topographic relief and geological structure. In Baoxing area the complex rocks correspond obvious high-velocity anomalies extending down to 15 km depth，while the Cenozoic rocks are correlated with low-velocity anomalies. Our high-resolution tomographic model not only displays the general features contained in the previous models, but also reveals some new features. An obvious high-velocity anomaly is visible in Daxing area. The high-velocity anomalies beneath Baoxing and Daxing connect each other in 10 km depth, which makes the contrast between high and low velocity anomalies more sharp. Above 20 km depth the velocity structure in southwest and northeast segment of the mainshock shows a big difference: low-velocity anomalies are dominated the southwest segment, while high-velocity anomalies rule the northeast segment. The Lushan mainshock locates at the leading edge of a low-velocity anomaly surrounded by the Baoxing and Daxing high-velocity anomalies. The Lushan aftershocks in southwest are distributed in low-velocity anomalies or the transition belt: the footwall represents low-velocity anomalies, while
An Automated 3d Indoor Topological Navigation Network Modelling
NASA Astrophysics Data System (ADS)
Jamali, A.; Rahman, A. A.; Boguslawski, P.; Gold, C. M.
2015-10-01
Indoor navigation is important for various applications such as disaster management and safety analysis. In the last decade, indoor environment has been a focus of wide research; that includes developing techniques for acquiring indoor data (e.g. Terrestrial laser scanning), 3D indoor modelling and 3D indoor navigation models. In this paper, an automated 3D topological indoor network generated from inaccurate 3D building models is proposed. In a normal scenario, 3D indoor navigation network derivation needs accurate 3D models with no errors (e.g. gap, intersect) and two cells (e.g. rooms, corridors) should touch each other to build their connections. The presented 3D modeling of indoor navigation network is based on surveying control points and it is less dependent on the 3D geometrical building model. For reducing time and cost of indoor building data acquisition process, Trimble LaserAce 1000 as surveying instrument is used. The modelling results were validated against an accurate geometry of indoor building environment which was acquired using Trimble M3 total station.
Highway 3D model from image and lidar data
NASA Astrophysics Data System (ADS)
Chen, Jinfeng; Chu, Henry; Sun, Xiaoduan
2014-05-01
We present a new method of highway 3-D model construction developed based on feature extraction in highway images and LIDAR data. We describe the processing road coordinate data that connect the image frames to the coordinates of the elevation data. Image processing methods are used to extract sky, road, and ground regions as well as significant objects (such as signs and building fronts) in the roadside for the 3D model. LIDAR data are interpolated and processed to extract the road lanes as well as other features such as trees, ditches, and elevated objects to form the 3D model. 3D geometry reasoning is used to match the image features to the 3D model. Results from successive frames are integrated to improve the final model.
Simulation of 3D infrared scenes using random fields model
NASA Astrophysics Data System (ADS)
Shao, Xiaopeng; Zhang, Jianqi
2001-09-01
Analysis and simulation of smart munitions requires imagery for the munition's sensor to view. The traditional infrared background simulations are always limited in the plane scene studies. A new method is described to synthesize the images in 3D view and with various terrains texture. We develop the random fields model and temperature fields to simulate 3D infrared scenes. Generalized long-correlation (GLC) model, one of random field models, will generate both the 3D terrains skeleton data and the terrains texture in this work. To build the terrain mesh with the random fields, digital elevation models (DEM) are introduced in the paper. And texture mapping technology will perform the task of pasting the texture in the concavo-convex surfaces of the 3D scene. The simulation using random fields model is a very available method to produce 3D infrared scene with great randomicity and reality.
An Automatic Registration Algorithm for 3D Maxillofacial Model
NASA Astrophysics Data System (ADS)
Qiu, Luwen; Zhou, Zhongwei; Guo, Jixiang; Lv, Jiancheng
2016-09-01
3D image registration aims at aligning two 3D data sets in a common coordinate system, which has been widely used in computer vision, pattern recognition and computer assisted surgery. One challenging problem in 3D registration is that point-wise correspondences between two point sets are often unknown apriori. In this work, we develop an automatic algorithm for 3D maxillofacial models registration including facial surface model and skull model. Our proposed registration algorithm can achieve a good alignment result between partial and whole maxillofacial model in spite of ambiguous matching, which has a potential application in the oral and maxillofacial reparative and reconstructive surgery. The proposed algorithm includes three steps: (1) 3D-SIFT features extraction and FPFH descriptors construction; (2) feature matching using SAC-IA; (3) coarse rigid alignment and refinement by ICP. Experiments on facial surfaces and mandible skull models demonstrate the efficiency and robustness of our algorithm.
Extending 3D city models with legal information
NASA Astrophysics Data System (ADS)
Frank, A. U.; Fuhrmann, T.; Navratil, G.
2012-10-01
3D city models represent existing physical objects and their topological and functional relations. In everyday life the rights and responsibilities connected to these objects, primarily legally defined rights and obligations but also other socially and culturally established rights, are of importance. The rights and obligations are defined in various laws and it is often difficult to identify the rules applicable for a certain case. The existing 2D cadastres show civil law rights and obligations and plans to extend them to provide information about public law restrictions for land use are in several countries under way. It is tempting to design extensions to the 3D city models to provide information about legal rights in 3D. The paper analyses the different types of information that are needed to reduce conflicts and to facilitate decisions about land use. We identify the role 3D city models augmented with planning information in 3D can play, but do not advocate a general conversion from 2D to 3D for the legal cadastre. Space is not anisotropic and the up/down dimension is practically very different from the two dimensional plane - this difference must be respected when designing spatial information systems. The conclusions are: (1) continue the current regime for ownership of apartments, which is not ownership of a 3D volume, but co-ownership of a building with exclusive use of some rooms; such exclusive use rights could be shown in a 3D city model; (2) ownership of 3D volumes for complex and unusual building situations can be reported in a 3D city model, but are not required everywhere; (3) indicate restrictions for land use and building in 3D city models, with links to the legal sources.
3d Velocity Tomography of The Kos - Nisyros Volcanic Area - East Aegean Sea
NASA Astrophysics Data System (ADS)
Nikolova, S.; Ilinski, D.; Makris, J.; Chonia, T.; Stavrakakis, J.
Since June 2000, active and passive seismic observations have been carried out by IfG, GeoPro GmbH, Hamburg and Institute of Geodynamics, Athens within the frame of the project GEOWARN (Geo-Spacial Warning Systems Nisyros Volcano, Greece: An Emergency Case Study of the Volcanic Area of Nisyros) supported by the European Community. In the active experiment 48 recording seismic units were deployed and recorded more than 7000 shots in 3D array. The Nisyros volcano has been identified as an apophytic intrusion of much larger volcanic structure with a caldera of 35 km diameter, extending between the southern coasts of the islands of Kos and Nisyros. To obtain 3-D velocity structure of the area a tomographic inversion was made using 6800 rays which probed the area with a very high ray density. The method applied and the high accuracy of active tomographic data allowed to resolve the high velocity bodies in the caldera. The complex volcanic structure is identified by high velocity rocks in- truding through the upper crust and penetrating the volcanic cone to depth of approx. 1.0 km to 1.8 km below the surface. Particularly high velocity bodies were identified below the islands of Yali and the central caldera of Nisyros. The high velocity bodies at shallow depth were interpreted as high-density cumulates of solidified magma intru- sion in the caldera. These intrusions explain very high temperature of 300C observed in the lower aquifer in the caldera at 1.5 km depth as confirmed by drilling. The vol- canic edifices of Kos, Yali, Nisyros and Strongily are part of a major volcanic caldera nearly 35 km in diameter. This size of the volcanic caldera explains the large volume of ignimbrites erupted 160 000 years ago. By combining geodetic, geophysical, geo- chemical and geological observations it is intended to correlate magma movements and associated changes of physical and chemical parameters of the recent volcanism.
3D scanning modeling method application in ancient city reconstruction
NASA Astrophysics Data System (ADS)
Ren, Pu; Zhou, Mingquan; Du, Guoguang; Shui, Wuyang; Zhou, Pengbo
2015-07-01
With the development of optical engineering technology, the precision of 3D scanning equipment becomes higher, and its role in 3D modeling is getting more distinctive. This paper proposed a 3D scanning modeling method that has been successfully applied in Chinese ancient city reconstruction. On one hand, for the existing architectures, an improved algorithm based on multiple scanning is adopted. Firstly, two pieces of scanning data were rough rigid registered using spherical displacers and vertex clustering method. Secondly, a global weighted ICP (iterative closest points) method is used to achieve a fine rigid registration. On the other hand, for the buildings which have already disappeared, an exemplar-driven algorithm for rapid modeling was proposed. Based on the 3D scanning technology and the historical data, a system approach was proposed for 3D modeling and virtual display of ancient city.
3-D model-based Bayesian classification
Soenneland, L.; Tenneboe, P.; Gehrmann, T.; Yrke, O.
1994-12-31
The challenging task of the interpreter is to integrate different pieces of information and combine them into an earth model. The sophistication level of this earth model might vary from the simplest geometrical description to the most complex set of reservoir parameters related to the geometrical description. Obviously the sophistication level also depend on the completeness of the available information. The authors describe the interpreter`s task as a mapping between the observation space and the model space. The information available to the interpreter exists in observation space and the task is to infer a model in model-space. It is well-known that this inversion problem is non-unique. Therefore any attempt to find a solution depend son constraints being added in some manner. The solution will obviously depend on which constraints are introduced and it would be desirable to allow the interpreter to modify the constraints in a problem-dependent manner. They will present a probabilistic framework that gives the interpreter the tools to integrate the different types of information and produce constrained solutions. The constraints can be adapted to the problem at hand.
Opportunity Landing Spot Panorama (3-D Model)
NASA Technical Reports Server (NTRS)
2004-01-01
The rocky outcrop traversed by the Mars Exploration Rover Opportunity is visible in this three-dimensional model of the rover's landing site. Opportunity has acquired close-up images along the way, and scientists are using the rover's instruments to closely examine portions of interest. The white fragments that look crumpled near the center of the image are portions of the airbags. Distant scenery is displayed on a spherical backdrop or 'billboard' for context. Artifacts near the top rim of the crater are a result of the transition between the three-dimensional model and the billboard. Portions of the terrain model lacking sufficient data appear as blank spaces or gaps, colored reddish-brown for better viewing. This image was generated using special software from NASA's Ames Research Center and a mosaic of images taken by the rover's panoramic camera.
[figure removed for brevity, see original site] Click on image for larger view
The rocky outcrop traversed by the Mars Exploration Rover Opportunity is visible in this zoomed-in portion of a three-dimensional model of the rover's landing site. Opportunity has acquired close-up images along the way, and scientists are using the rover's instruments to closely examine portions of interest. The white fragments that look crumpled near the center of the image are portions of the airbags. Distant scenery is displayed on a spherical backdrop or 'billboard' for context. Artifacts near the top rim of the crater are a result of the transition between the three-dimensional model and the billboard. Portions of the terrain model lacking sufficient data appear as blank spaces or gaps, colored reddish-brown for better viewing. This image was generated using special software from NASA's Ames Research Center and a mosaic of images taken by the rover's panoramic camera.
Virtual 3d City Modeling: Techniques and Applications
NASA Astrophysics Data System (ADS)
Singh, S. P.; Jain, K.; Mandla, V. R.
2013-08-01
3D city model is a digital representation of the Earth's surface and it's related objects such as Building, Tree, Vegetation, and some manmade feature belonging to urban area. There are various terms used for 3D city models such as "Cybertown", "Cybercity", "Virtual City", or "Digital City". 3D city models are basically a computerized or digital model of a city contains the graphic representation of buildings and other objects in 2.5 or 3D. Generally three main Geomatics approach are using for Virtual 3-D City models generation, in first approach, researcher are using Conventional techniques such as Vector Map data, DEM, Aerial images, second approach are based on High resolution satellite images with LASER scanning, In third method, many researcher are using Terrestrial images by using Close Range Photogrammetry with DSM & Texture mapping. We start this paper from the introduction of various Geomatics techniques for 3D City modeling. These techniques divided in to two main categories: one is based on Automation (Automatic, Semi-automatic and Manual methods), and another is Based on Data input techniques (one is Photogrammetry, another is Laser Techniques). After details study of this, finally in short, we are trying to give the conclusions of this study. In the last, we are trying to give the conclusions of this research paper and also giving a short view for justification and analysis, and present trend for 3D City modeling. This paper gives an overview about the Techniques related with "Generation of Virtual 3-D City models using Geomatics Techniques" and the Applications of Virtual 3D City models. Photogrammetry, (Close range, Aerial, Satellite), Lasergrammetry, GPS, or combination of these modern Geomatics techniques play a major role to create a virtual 3-D City model. Each and every techniques and method has some advantages and some drawbacks. Point cloud model is a modern trend for virtual 3-D city model. Photo-realistic, Scalable, Geo-referenced virtual 3
Graves, Robert W.; Aagaard, Brad T.
2011-01-01
Using a suite of five hypothetical finite-fault rupture models, we test the ability of long-period (T>2.0 s) ground-motion simulations of scenario earthquakes to produce waveforms throughout southern California consistent with those recorded during the 4 April 2010 Mw 7.2 El Mayor-Cucapah earthquake. The hypothetical ruptures are generated using the methodology proposed by Graves and Pitarka (2010) and require, as inputs, only a general description of the fault location and geometry, event magnitude, and hypocenter, as would be done for a scenario event. For each rupture model, two Southern California Earthquake Center three-dimensional community seismic velocity models (CVM-4m and CVM-H62) are used, resulting in a total of 10 ground-motion simulations, which we compare with recorded ground motions. While the details of the motions vary across the simulations, the median levels match the observed peak ground velocities reasonably well, with the standard deviation of the residuals generally within 50% of the median. Simulations with the CVM-4m model yield somewhat lower variance than those with the CVM-H62 model. Both models tend to overpredict motions in the San Diego region and underpredict motions in the Mojave desert. Within the greater Los Angeles basin, the CVM-4m model generally matches the level of observed motions, whereas the CVM-H62 model tends to overpredict the motions, particularly in the southern portion of the basin. The variance in the peak velocity residuals is lowest for a rupture that has significant shallow slip (<5 km depth), whereas the variance in the residuals is greatest for ruptures with large asperities below 10 km depth. Overall, these results are encouraging and provide confidence in the predictive capabilities of the simulation methodology, while also suggesting some regions in which the seismic velocity models may need improvement.
Venusian Applications of 3D Convection Modeling
NASA Technical Reports Server (NTRS)
Bonaccorso, Timary Annie
2011-01-01
This study models mantle convection on Venus using the 'cubed sphere' code OEDIPUS, which models one-sixth of the planet in spherical geometry. We are attempting to balance internal heating, bottom mantle viscosity, and temperature difference across Venus' mantle, in order to create a realistic model that matches with current planetary observations. We also have begun to run both lower and upper mantle simulations to determine whether layered (as opposed to whole-mantle) convection might produce more efficient heat transfer, as well as to model coronae formation in the upper mantle. Upper mantle simulations are completed using OEDIPUS' Cartesian counterpart, JOCASTA. This summer's central question has been how to define a mantle plume. Traditionally, we have defined a hot plume the region with temperature at or above 40% of the difference between the maximum and horizontally averaged temperature, and a cold plume as the region with 40% of the difference between the minimum and average temperature. For less viscous cases (1020 Pa?s), the plumes generated by that definition lacked vigor, displaying buoyancies 1/100th of those found in previous, higher viscosity simulations (1021 Pa?s). As the mantle plumes with large buoyancy flux are most likely to produce topographic uplift and volcanism, the low viscosity cases' plumes may not produce observable deformation. In an effort to eliminate the smallest plumes, we experimented with different lower bound parameters and temperature percentages.
Image based 3D city modeling : Comparative study
NASA Astrophysics Data System (ADS)
Singh, S. P.; Jain, K.; Mandla, V. R.
2014-06-01
3D city model is a digital representation of the Earth's surface and it's related objects such as building, tree, vegetation, and some manmade feature belonging to urban area. The demand of 3D city modeling is increasing rapidly for various engineering and non-engineering applications. Generally four main image based approaches were used for virtual 3D city models generation. In first approach, researchers were used Sketch based modeling, second method is Procedural grammar based modeling, third approach is Close range photogrammetry based modeling and fourth approach is mainly based on Computer Vision techniques. SketchUp, CityEngine, Photomodeler and Agisoft Photoscan are the main softwares to represent these approaches respectively. These softwares have different approaches & methods suitable for image based 3D city modeling. Literature study shows that till date, there is no complete such type of comparative study available to create complete 3D city model by using images. This paper gives a comparative assessment of these four image based 3D modeling approaches. This comparative study is mainly based on data acquisition methods, data processing techniques and output 3D model products. For this research work, study area is the campus of civil engineering department, Indian Institute of Technology, Roorkee (India). This 3D campus acts as a prototype for city. This study also explains various governing parameters, factors and work experiences. This research work also gives a brief introduction, strengths and weakness of these four image based techniques. Some personal comment is also given as what can do or what can't do from these softwares. At the last, this study shows; it concluded that, each and every software has some advantages and limitations. Choice of software depends on user requirements of 3D project. For normal visualization project, SketchUp software is a good option. For 3D documentation record, Photomodeler gives good result. For Large city
NoSQL Based 3D City Model Management System
NASA Astrophysics Data System (ADS)
Mao, B.; Harrie, L.; Cao, J.; Wu, Z.; Shen, J.
2014-04-01
To manage increasingly complicated 3D city models, a framework based on NoSQL database is proposed in this paper. The framework supports import and export of 3D city model according to international standards such as CityGML, KML/COLLADA and X3D. We also suggest and implement 3D model analysis and visualization in the framework. For city model analysis, 3D geometry data and semantic information (such as name, height, area, price and so on) are stored and processed separately. We use a Map-Reduce method to deal with the 3D geometry data since it is more complex, while the semantic analysis is mainly based on database query operation. For visualization, a multiple 3D city representation structure CityTree is implemented within the framework to support dynamic LODs based on user viewpoint. Also, the proposed framework is easily extensible and supports geoindexes to speed up the querying. Our experimental results show that the proposed 3D city management system can efficiently fulfil the analysis and visualization requirements.
Modelling Polymer Deformation during 3D Printing
NASA Astrophysics Data System (ADS)
McIlroy, Claire; Olmsted, Peter
Three-dimensional printing has the potential to transform manufacturing processes, yet improving the strength of printed parts, to equal that of traditionally-manufactured parts, remains an underlying issue. The fused deposition modelling technique involves melting a thermoplastic, followed by layer-by-layer extrusion to fabricate an object. The key to ensuring strength at the weld between layers is successful inter-diffusion. However, prior to welding, both the extrusion process and the cooling temperature profile can significantly deform the polymer micro-structure and, consequently, how well the polymers are able to ``re-entangle'' across the weld. In particular, polymer alignment in the flow can cause de-bonding of the layers and create defects. We have developed a simple model of the non-isothermal extrusion process to explore the effects that typical printing conditions and material rheology have on the conformation of a polymer melt. In particular, we incorporate both stretch and orientation using the Rolie-Poly constitutive equation to examine the melt structure as it flows through the nozzle, the subsequent alignment with the build plate and the resulting deformation due to the fixed nozzle height, which is typically less than the nozzle radius.
The Prevalence of Gas Outflows in Type 2 AGNs. II. 3D Biconical Outflow Models
NASA Astrophysics Data System (ADS)
Bae, Hyun-Jin; Woo, Jong-Hak
2016-09-01
We present 3D models of biconical outflows combined with a thin dust plane for investigating the physical properties of the ionized gas outflows and their effect on the observed gas kinematics in type 2 active galactic nuclei (AGNs). Using a set of input parameters, we construct a number of models in 3D and calculate the spatially integrated velocity and velocity dispersion for each model. We find that three primary parameters, i.e., intrinsic velocity, bicone inclination, and the amount of dust extinction, mainly determine the simulated velocity and velocity dispersion. Velocity dispersion increases as the intrinsic velocity or the bicone inclination increases, while velocity (i.e., velocity shifts with respect to systemic velocity) increases as the amount of dust extinction increases. Simulated emission-line profiles well reproduce the observed [O iii] line profiles, e.g., narrow core and broad wing components. By comparing model grids and Monte Carlo simulations with the observed [O iii] velocity–velocity dispersion distribution of ∼39,000 type 2 AGNs, we constrain the intrinsic velocity of gas outflows ranging from ∼500 to ∼1000 km s‑1 for the majority of AGNs, and up to ∼1500–2000 km s‑1 for extreme cases. The Monte Carlo simulations show that the number ratio of AGNs with negative [O iii] velocity to AGNs with positive [O iii] velocity correlates with the outflow opening angle, suggesting that outflows with higher intrinsic velocity tend to have wider opening angles. These results demonstrate the potential of our 3D models for studying the physical properties of gas outflows, applicable to various observations, including spatially integrated and resolved gas kinematics.
3D PIC Modeling of Microcavity Discharge
NASA Astrophysics Data System (ADS)
Hopkins, Matthew; Manginell, Ronald; Moore, Christopher; Yee, Benjamin; Moorman, Matthew
2015-09-01
We present a number of techniques and challenges in simulating the transient behavior of a microcavity discharge. Our microcavities are typically cylindrical with diameters approximately 50 - 100 μm, heights of 50 - 200 μm, pressure near atmospheric, and operate at a few hundred volts. We employ a fully kinetic simulation methodology, the Particle-in-Cell (PIC) method, with interparticle collisions handled via methods based on direct simulation Monte Carlo (DSMC). In particular, we explicitly include kinetic electrons. Some of the challenges we encounter include variations in number densities, external circuit coupling, and time step resolution constraints. By employing dynamic particle weighting (particle weights vary over time by species and location) we can mitigate some of the challenges modeling systems with 107 variations in number densities. Smoothing mechanisms have been used to attempt to mitigate external circuit response. We perform our simulations on hundreds or thousands of processing cores to accommodate the computational work inherent in using relatively small time step sizes (e.g., 50 fs for a 100 ns calculation). In addition, particle weighting issues inherent to three-dimensional low temperature plasma systems will be mentioned. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's NNSA under Contract DE-AC04-94AL85000.
3-D Eutrophication Modeling for Lake Simcoe, Canada
NASA Astrophysics Data System (ADS)
Lu, Q.; Duckett, F.; Nairn, R.; Brunton, A.
2006-12-01
The Lake Simcoe Region Conservation Authority (LSRCA) and the Province of Ontario are undertaking a series of studies to facilitate management of the pressures of population growth in the Lake Simcoe watershed. With rapid population growth and urban development comes additional land clearing, storm water runoff and the discharge of treated sewage, all of which are sources of increased phosphorus loading to Lake Simcoe. Depressed oxygen levels were linked to phosphorous enrichment of the lake, with the resultant stimulation of algal growth in the sunlit upper waters of the lake, and its subsequent senescence and settling into the hypolimnion where bacterial decomposition consumes oxygen from the stratified waters. This poster describes a 3-D hydrodynamic, thermal and water quality model of Lake Simcoe developed using the Danish Hydraulics Institute (DHI) MIKE3 model. The hydrodynamic module includes wind-driven circulation, temperature variation, development of the thermocline and thermal stratification, and hydraulic forcing from inflowing tributaries. This is linked to the water quality module which simulates the eutrophication processes in the response of the lake to loadings of phosphorus, such as algal growth, the growth of aquatic plants and subsequent oxygen consumption. The model has been calibrated against Acoustic Doppler Current Profiler velocity data, plus measured temperature and water quality data at MOE stations in the lake and water intakes. The model is an important assessment tool for the management of the lake and its watersheds, allowing assessment of the impacts of the urban growth and land use change on the water quality in Lake Simcoe.
The 3D model: explaining densification and deformation mechanisms by using 3D parameter plots.
Picker, Katharina M
2004-04-01
The aim of the study was to analyze very differently deforming materials using 3D parameter plots and consequently to gain deeper insights into the densification and deformation process described with the 3D model in order to define an ideal tableting excipient. The excipients used were dicalcium phosphate dihydrate (DCPD), sodium chloride (NaCl), microcrystalline cellulose (MCC), xylitol, mannitol, alpha-lactose monohydrate, maltose, hydroxypropyl methylcellulose (HPMC), sodium carboxymethylcellulose (NaCMC), cellulose acetate (CAC), maize starch, potato starch, pregelatinized starch, and maltodextrine. All of the materials were tableted to graded maximum relative densities (rhorel, max) using an eccentric tableting machine. The data which resulted, namely force, displacement, and time, were analyzed by the application of 3D modeling. Different particle size fractions of DCPD, CAC, and MCC were analyzed in addition. Brittle deforming materials such as DCPD exhibited a completely different 3D parameter plot, with low time plasticity, d, and low pressure plasticity, e, and a strong decrease in omega values when densification increased, in contrast to the plastically deforming MCC, which had much higher d, e, and omega values. e and omega values changed only slightly when densification increased for MCC. NaCl showed less of a decrease in omega values than DCPD did, and the d and e values were between those of MCC and DCPD. The sugar alcohols, xylitol and mannitol, behaved in a similar fashion to sodium chloride. This is also valid for the crystalline sugars, alpha-lactose monohydrate, and maltose. However, the sugars are more brittle than the sugar alcohols. The cellulose derivatives, HPMC, NaCMC, and CAC, are as plastic as MCC, however, their elasticity depends on substitution indicated by lower (more elastic) or higher (less elastic) omega values. The native starches, maize starch and potato starch, are very elastic, and pregelatinized starch and maltodextrine are
Verification and Validation of the k-kL Turbulence Model in FUN3D and CFL3D Codes
NASA Technical Reports Server (NTRS)
Abdol-Hamid, Khaled S.; Carlson, Jan-Renee; Rumsey, Christopher L.
2015-01-01
The implementation of the k-kL turbulence model using multiple computational uid dy- namics (CFD) codes is reported herein. The k-kL model is a two-equation turbulence model based on Abdol-Hamid's closure and Menter's modi cation to Rotta's two-equation model. Rotta shows that a reliable transport equation can be formed from the turbulent length scale L, and the turbulent kinetic energy k. Rotta's equation is well suited for term-by-term mod- eling and displays useful features compared to other two-equation models. An important di erence is that this formulation leads to the inclusion of higher-order velocity derivatives in the source terms of the scale equations. This can enhance the ability of the Reynolds- averaged Navier-Stokes (RANS) solvers to simulate unsteady ows. The present report documents the formulation of the model as implemented in the CFD codes Fun3D and CFL3D. Methodology, veri cation and validation examples are shown. Attached and sepa- rated ow cases are documented and compared with experimental data. The results show generally very good comparisons with canonical and experimental data, as well as matching results code-to-code. The results from this formulation are similar or better than results using the SST turbulence model.
3D Droplet velocities and sizes in the Ranque-Hilsch vortex tube
NASA Astrophysics Data System (ADS)
Liew, R.; Zeegers, J. C. H.; Kuerten, J. G. M.; Michalek, W. R.
2012-11-01
The Ranque-Hilsch vortex tube is a known device that is used to generate spot cooling. In this study, we experimentally investigate the behavior of small water droplets in the vortex tube by means of Phase Doppler Particle Analysis. In an experimental vortex tube, droplets were injected together with a carrier gas to form a fast rotating (up to 80.000 rpm) droplet-gas mixture. Droplet sizes, 3D velocity components, and turbulent properties were measured, showing high intensity isotropic turbulence in the core region. To investigate the cause of the high intensity turbulence, a frequency analysis was applied on the measured velocity. The frequency spectrum of the velocity is presented and indicates that wobbling of the vortex axis is the cause of the high turbulence intensity. It was expected that larger droplets have a higher radial velocity because of the larger centrifugal force. Results show, however, that small and lager droplets behave similar. This research is supported by the Dutch Technology Foundation STW, which is the applied science division of NWO, and the Technology Programme of the Ministry of Economic Affairs.
Numerical Results of 3-D Modeling of Moon Accumulation
NASA Astrophysics Data System (ADS)
Khachay, Yurie; Anfilogov, Vsevolod; Antipin, Alexandr
2014-05-01
For the last time for the model of the Moon usually had been used the model of mega impact in which the forming of the Earth and its sputnik had been the consequence of the Earth's collision with the body of Mercurial mass. But all dynamical models of the Earth's accumulation and the estimations after the Pb-Pb system, lead to the conclusion that the duration of the planet accumulation was about 1 milliard years. But isotopic results after the W-Hf system testify about a very early (5-10) million years, dividing of the geochemical reservoirs of the core and mantle. In [1,2] it is shown, that the account of energy dissipating by the decay of short living radioactive elements and first of all Al26,it is sufficient for heating even small bodies with dimensions about (50-100) km up to the iron melting temperature and can be realized a principal new differentiation mechanism. The inner parts of the melted preplanets can join and they are mainly of iron content, but the cold silicate fragments return to the supply zone and additionally change the content of Moon forming to silicates. Only after the increasing of the gravitational radius of the Earth, the growing area of the future Earth's core can save also the silicate envelope fragments [3]. For understanding the further system Earth-Moon evolution it is significant to trace the origin and evolution of heterogeneities, which occur on its accumulation stage.In that paper we are modeling the changing of temperature,pressure,velocity of matter flowing in a block of 3d spherical body with a growing radius. The boundary problem is solved by the finite-difference method for the system of equations, which include equations which describe the process of accumulation, the Safronov equation, the equation of impulse balance, equation Navier-Stocks, equation for above litho static pressure and heat conductivity in velocity-pressure variables using the Businesque approach.The numerical algorithm of the problem solution in velocity
On the Quality of Velocity Interpolation Schemes for Marker-In-Cell Methods on 3-D Staggered Grids
NASA Astrophysics Data System (ADS)
Kaus, B.; Pusok, A. E.; Popov, A.
2015-12-01
The marker-in-cell method is generally considered to be a flexible and robust method to model advection of heterogenous non-diffusive properties (i.e. rock type or composition) in geodynamic problems or incompressible Stokes problems. In this method, Lagrangian points carrying compositional information are advected with the ambient velocity field on an immobile, Eulerian grid. However, velocity interpolation from grid points to marker locations is often performed without preserving the zero divergence of the velocity field at the interpolated locations (i.e. non-conservative). Such interpolation schemes can induce non-physical clustering of markers when strong velocity gradients are present (Jenny et al., 2001) and this may, eventually, result in empty grid cells, a serious numerical violation of the marker-in-cell method. Solutions to this problem include: using larger mesh resolutions and/or marker densities, or repeatedly controlling the marker distribution (i.e. inject/delete), but which does not have an established physical background. To remedy this at low computational costs, Jenny et al. (2001) and Meyer and Jenny (2004) proposed a simple, conservative velocity interpolation (CVI) scheme for 2-D staggered grid, while Wang et al. (2015) extended the formulation to 3-D finite element methods. Here, we follow up with these studies and report on the quality of velocity interpolation methods for 2-D and 3-D staggered grids. We adapt the formulations from both Jenny et al. (2001) and Wang et al. (2015) for use on 3-D staggered grids, where the velocity components have different node locations as compared to finite element, where they share the same node location. We test the different interpolation schemes (CVI and non-CVI) in combination with different advection schemes (Euler, RK2 and RK4) and with/out marker control on Stokes problems with strong velocity gradients, which are discretized using a finite difference method. We show that a conservative formulation
NASA Astrophysics Data System (ADS)
Obermann, Anne; Planès, Thomas; Hadziioannou, Céline; Campillo, Michel
2016-07-01
In the context of seismic monitoring, recent studies made successful use of seismic coda waves to locate medium changes on the horizontal plane. Locating the depth of the changes, however, remains a challenge. In this paper, we use 3-D wavefield simulations to address two problems: firstly, we evaluate the contribution of surface and body wave sensitivity to a change at depth. We introduce a thin layer with a perturbed velocity at different depths and measure the apparent relative velocity changes due to this layer at different times in the coda and for different degrees of heterogeneity of the model. We show that the depth sensitivity can be modelled as a linear combination of body- and surface-wave sensitivity. The lapse-time dependent sensitivity ratio of body waves and surface waves can be used to build 3-D sensitivity kernels for imaging purposes. Secondly, we compare the lapse-time behavior in the presence of a perturbation in horizontal and vertical slabs to address, for instance, the origin of the velocity changes detected after large earthquakes.
NASA Astrophysics Data System (ADS)
Yu, Xiangwei; Wang, Xiaona; Zhang, Wenbo
2016-04-01
Many researchers have investigated the Lushan source area with geological and geophysical approaches since the 2013 Lushan, China, earthquake happened. Compared with the previous tomographic studies, we have used a much large data set and an updated tomographic method to determine a small scale three-dimensional P wave velocity structure with spatial resolution less than 5km, which plays the important role for understanding the deep structure and the genetic mechanism beneath the Lushan area. The double difference seismic tomography method is applied to 50,711 absolute first arrival P wave arrival times and 7,294,691 high quality relative P arrival times of 5,285 events of Lushan seismic sequence to simultaneously determine the detailed crustal 3D P wave velocity structure and the hypocenter parameters in the Lushan seismic area. This method takes account of the path anomaly biases explicitly by making full use of valuable information of seismic wave propagation jointly with absolute and relative arrival time data. Our results show that the Lushan mainshock locates at 30.28N, 103.98E, with the depth of 16.38km. The front edge of aftershock in the northeast of mainshock present a spade with a steep dip angle, the aftershocks' extended length is about 12km. In the southwest of Lushan mainshock, the front edge of aftershock in low velocity zone slope gently, the aftershocks' extended length is about 23km. Our high-resolution tomographic model not only displays the general features contained in the previous models, but also reveals some new features. The Tianquan, Shuangshi and Daguan line lies in the transition zone between high velocity anomalies to the southeast and low velocity anomalies to the northwest at the ground surface. An obvious high-velocity anomaly is visible in Daxing area. With the depth increasing, Baoxing high velocity anomaly extends to Lingguan, while the southeast of the Tianquan, Shuangshi and Daguan line still shows low velocity. The high-velocity
Collision of continental corner from 3-D numerical modeling
NASA Astrophysics Data System (ADS)
Li, Zhong-Hai; Xu, Zhiqin; Gerya, Taras; Burg, Jean-Pierre
2013-10-01
Continental collision has been extensively investigated with 2-D numerical models assuming infinitely wide plates or insignificant along-strike deformation in the third dimension. However, the corners of natural collision zones normally have structural characteristics that differ from linear parts of mountain belt. We conducted 3-D high-resolution numerical simulations to study the dynamics of a continental corner (lateral continental/oceanic transition zone) during subduction/collision. The results demonstrate different modes between the oceanic subduction side (continuous subduction and retreating trench) and the continental collision side (slab break-off and topography uplift). Slab break-off occurs at a depth (⩽100 km to ˜300 km) that depends on the convergence velocity. The numerical models produce lateral extrusion of the overriding crust from the collisional side to the subduction side, which is also a phenomenon recognized around natural collision of continental corners, for instance around the western corner of the Arabia-Asia collision zone and around the eastern corner of the India-Asia collision zone. Modeling results also indicate that extrusion tectonics may be driven both from above by the topography and gravitational potentials and from below by the trench retreat and asthenospheric mantle return flow, which supports the link between deep mantle dynamics and shallower crustal deformation.
3D-model building of the jaw impression
NASA Astrophysics Data System (ADS)
Ahmed, Moumen T.; Yamany, Sameh M.; Hemayed, Elsayed E.; Farag, Aly A.
1997-03-01
A novel approach is proposed to obtain a record of the patient's occlusion using computer vision. Data acquisition is obtained using intra-oral video cameras. The technique utilizes shape from shading to extract 3D information from 2D views of the jaw, and a novel technique for 3D data registration using genetic algorithms. The resulting 3D model can be used for diagnosis, treatment planning, and implant purposes. The overall purpose of this research is to develop a model-based vision system for orthodontics to replace traditional approaches. This system will be flexible, accurate, and will reduce the cost of orthodontic treatments.
3D model-based still image object categorization
NASA Astrophysics Data System (ADS)
Petre, Raluca-Diana; Zaharia, Titus
2011-09-01
This paper proposes a novel recognition scheme algorithm for semantic labeling of 2D object present in still images. The principle consists of matching unknown 2D objects with categorized 3D models in order to infer the semantics of the 3D object to the image. We tested our new recognition framework by using the MPEG-7 and Princeton 3D model databases in order to label unknown images randomly selected from the web. Results obtained show promising performances, with recognition rate up to 84%, which opens interesting perspectives in terms of semantic metadata extraction from still images/videos.
A 3D numerical model for Kepler's supernova remnant
NASA Astrophysics Data System (ADS)
Toledo-Roy, J. C.; Esquivel, A.; Velázquez, P. F.; Reynoso, E. M.
2014-07-01
We present new 3D numerical simulations for Kepler's supernova remnant. In this work we revisit the possibility that the asymmetric shape of the remnant in X-rays is the product of a Type Ia supernova explosion which occurs inside the wind bubble previously created by an AGB companion star. Due to the large peculiar velocity of the system, the interaction of the strong AGB wind with the interstellar medium results in a bow shock structure. In this new model we propose that the AGB wind is anisotropic, with properties such as mass-loss rate and density having a latitude dependence, and that the orientation of the polar axis of the AGB star is not aligned with the direction of motion. The ejecta from the Type Ia supernova explosion is modelled using a power-law density profile, and we let the remnant evolve for 400 yr. We computed synthetic X-ray maps from the numerical results. We find that the estimated size and peculiar X-ray morphology of Kepler's supernova remnant are well reproduced by considering an AGB mass-loss rate of 10-5 M⊙ yr-1, a wind terminal velocity of 10 km s-1, an ambient medium density of 10-3 cm-3 and an explosion energy of 7 × 1050 erg. The obtained total X-ray luminosity of the remnant in this model reaches 6 × 1050 erg, which is within a factor of 2 of the observed value, and the time evolution of the luminosity shows a rate of decrease in recent decades of ˜2.4 per cent yr-1 that is consistent with the observations.
Summary on Several Key Techniques in 3D Geological Modeling
2014-01-01
Several key techniques in 3D geological modeling including planar mesh generation, spatial interpolation, and surface intersection are summarized in this paper. Note that these techniques are generic and widely used in various applications but play a key role in 3D geological modeling. There are two essential procedures in 3D geological modeling: the first is the simulation of geological interfaces using geometric surfaces and the second is the building of geological objects by means of various geometric computations such as the intersection of surfaces. Discrete geometric surfaces that represent geological interfaces can be generated by creating planar meshes first and then spatially interpolating; those surfaces intersect and then form volumes that represent three-dimensional geological objects such as rock bodies. In this paper, the most commonly used algorithms of the key techniques in 3D geological modeling are summarized. PMID:24772029
Formal representation of 3D structural geological models
NASA Astrophysics Data System (ADS)
Wang, Zhangang; Qu, Honggang; Wu, Zixing; Yang, Hongjun; Du, Qunle
2016-05-01
The development and widespread application of geological modeling methods has increased demands for the integration and sharing services of three dimensional (3D) geological data. However, theoretical research in the field of geological information sciences is limited despite the widespread use of Geographic Information Systems (GIS) in geology. In particular, fundamental research on the formal representations and standardized spatial descriptions of 3D structural models is required. This is necessary for accurate understanding and further applications of geological data in 3D space. In this paper, we propose a formal representation method for 3D structural models using the theory of point set topology, which produces a mathematical definition for the major types of geological objects. The spatial relationships between geologic boundaries, structures, and units are explained in detail using the 9-intersection model. Reasonable conditions for describing the topological space of 3D structural models are also provided. The results from this study can be used as potential support for the standardized representation and spatial quality evaluation of 3D structural models, as well as for specific needs related to model-based management, query, and analysis.
Metrological validation for 3D modeling of dental plaster casts.
Brusco, Nicola; Andreetto, Marco; Lucchese, Luca; Carmignato, Simone; Cortelazzo, Guido M
2007-11-01
The contribution of this paper is twofold: (1) it presents an automatic 3D modeling technique and (2) it advances a procedure for its metrological evaluation in the context of a medical application, the 3D modeling of dental plaster casts. The motivation for this work is the creation of a "virtual gypsotheque" where cumbersome dental plaster casts can be replaced by numerical 3D models, thereby alleviating storage and access problems and allowing dentists and orthodontists the use of novel and unprecedented software tools for their medical evaluations. Modeling free-form surfaces of anatomical interest is an intriguing mixture of open issues concerning 3D modeling, geometrical metrology, and medicine. Of general interest is both the fact that a widespread use of 3D modeling in non-engineering applications requires automatic procedures of the kind presented in this work and the adopted validation paradigm for free-form surfaces, rather useful for practical purposes. In this latter respect, the metrological analysis we advance is the first seminal attempt in the field of 3D modeling and can be readily extended to contexts other than the medical one discussed in this paper. PMID:17126062
3D Modeling from Photos Given Topological Information.
Kim, Young Min; Cho, Junghyun; Ahn, Sang Chul
2016-09-01
Reconstructing 3D models given a single-view 2D information is inherently an ill-posed problem and requires additional information such as shape prior or user input.We introduce a method to generate multiple 3D models of a particular category given corresponding photographs when the topological information is known. While there is a wide range of shapes for an object of a particular category, the basic topology usually remains constant.In consequence, the topological prior needs to be provided only once for each category and can be easily acquired by consulting an existing database of 3D models or by user input. The input of topological description is only connectivity information between parts; this is in contrast to previous approaches that have required users to interactively mark individual parts. Given the silhouette of an object and the topology, our system automatically finds a skeleton and generates a textured 3D model by jointly fitting multiple parts. The proposed method, therefore, opens the possibility of generating a large number of 3D models by consulting a massive number of photographs. We demonstrate examples of the topological prior and reconstructed 3D models using photos. PMID:26661474
Performance Evaluation of 3d Modeling Software for Uav Photogrammetry
NASA Astrophysics Data System (ADS)
Yanagi, H.; Chikatsu, H.
2016-06-01
UAV (Unmanned Aerial Vehicle) photogrammetry, which combines UAV and freely available internet-based 3D modeling software, is widely used as a low-cost and user-friendly photogrammetry technique in the fields such as remote sensing and geosciences. In UAV photogrammetry, only the platform used in conventional aerial photogrammetry is changed. Consequently, 3D modeling software contributes significantly to its expansion. However, the algorithms of the 3D modelling software are black box algorithms. As a result, only a few studies have been able to evaluate their accuracy using 3D coordinate check points. With this motive, Smart3DCapture and Pix4Dmapper were downloaded from the Internet and commercial software PhotoScan was also employed; investigations were performed in this paper using check points and images obtained from UAV.
Automatic Texture Mapping of Architectural and Archaeological 3d Models
NASA Astrophysics Data System (ADS)
Kersten, T. P.; Stallmann, D.
2012-07-01
Today, detailed, complete and exact 3D models with photo-realistic textures are increasingly demanded for numerous applications in architecture and archaeology. Manual texture mapping of 3D models by digital photographs with software packages, such as Maxon Cinema 4D, Autodesk 3Ds Max or Maya, still requires a complex and time-consuming workflow. So, procedures for automatic texture mapping of 3D models are in demand. In this paper two automatic procedures are presented. The first procedure generates 3D surface models with textures by web services, while the second procedure textures already existing 3D models with the software tmapper. The program tmapper is based on the Multi Layer 3D image (ML3DImage) algorithm and developed in the programming language C++. The studies showing that the visibility analysis using the ML3DImage algorithm is not sufficient to obtain acceptable results of automatic texture mapping. To overcome the visibility problem the Point Cloud Painter algorithm in combination with the Z-buffer-procedure will be applied in the future.
Gis-Based Smart Cartography Using 3d Modeling
NASA Astrophysics Data System (ADS)
Malinverni, E. S.; Tassetti, A. N.
2013-08-01
3D City Models have evolved to be important tools for urban decision processes and information systems, especially in planning, simulation, analysis, documentation and heritage management. On the other hand existing and in use numerical cartography is often not suitable to be used in GIS because not geometrically and topologically correctly structured. The research aim is to 3D structure and organize a numeric cartography for GIS and turn it into CityGML standardized features. The work is framed around a first phase of methodological analysis aimed to underline which existing standard (like ISO and OGC rules) can be used to improve the quality requirement of a cartographic structure. Subsequently, from this technical specifics, it has been investigated the translation in formal contents, using an owner interchange software (SketchUp), to support some guide lines implementations to generate a GIS3D structured in GML3. It has been therefore predisposed a test three-dimensional numerical cartography (scale 1:500, generated from range data captured by 3D laser scanner), tested on its quality according to the previous standard and edited when and where necessary. Cad files and shapefiles are converted into a final 3D model (Google SketchUp model) and then exported into a 3D city model (CityGML LoD1/LoD2). The GIS3D structure has been managed in a GIS environment to run further spatial analysis and energy performance estimate, not achievable in a 2D environment. In particular geometrical building parameters (footprint, volume etc.) are computed and building envelop thermal characteristics are derived from. Lastly, a simulation is carried out to deal with asbestos and home renovating charges and show how the built 3D city model can support municipal managers with risk diagnosis of the present situation and development of strategies for a sustainable redevelop.
The 3D Space and Spin Velocities of a Gamma-ray Pulsar
NASA Astrophysics Data System (ADS)
Romani, Roger W.
2016-04-01
PSR J2030+4415 is a LAT-discovered 0.5My-old gamma-ray pulsar with an X-ray synchrotron trail and a rare Halpha bowshock. We have obtained GMOS IFU spectroscopic imaging of this shell, and show a sweep through the remarkable Halpha structure, comparing with the high energy emission. These data provide a unique 3D map of the momentum distribution of the relativistic pulsar wind. This shows that the pulsar is moving nearly in the plane of the sky and that the pulsar wind has a polar component misaligned with the space velocity. The spin axis is shown to be inclined some 95degrees to the Earth line of sight, explaining why this is a radio-quiet, gamma-only pulsar. Intriguingly, the shell also shows multiple bubbles that suggest that the pulsar wind power has varied substantially over the past 500 years.
Combined registration of 3D tibia and femur implant models in 3D magnetic resonance images
NASA Astrophysics Data System (ADS)
Englmeier, Karl-Hans; Siebert, Markus; von Eisenhart-Rothe, Ruediger; Graichen, Heiko
2008-03-01
The most frequent reasons for revision of total knee arthroplasty are loosening and abnormal axial alignment leading to an unphysiological kinematic of the knee implant. To get an idea about the postoperative kinematic of the implant, it is essential to determine the position and orientation of the tibial and femoral prosthesis. Therefore we developed a registration method for fitting 3D CAD-models of knee joint prostheses into an 3D MR image. This rigid registration is the basis for a quantitative analysis of the kinematics of knee implants. Firstly the surface data of the prostheses models are converted into a voxel representation; a recursive algorithm determines all boundary voxels of the original triangular surface data. Secondly an initial preconfiguration of the implants by the user is still necessary for the following step: The user has to perform a rough preconfiguration of both remaining prostheses models, so that the fine matching process gets a reasonable starting point. After that an automated gradient-based fine matching process determines the best absolute position and orientation: This iterative process changes all 6 parameters (3 rotational- and 3 translational parameters) of a model by a minimal amount until a maximum value of the matching function is reached. To examine the spread of the final solutions of the registration, the interobserver variability was measured in a group of testers. This variability, calculated by the relative standard deviation, improved from about 50% (pure manual registration) to 0.5% (rough manual preconfiguration and subsequent fine registration with the automatic fine matching process).
3D Bioprinting of Tissue/Organ Models.
Pati, Falguni; Gantelius, Jesper; Svahn, Helene Andersson
2016-04-01
In vitro tissue/organ models are useful platforms that can facilitate systematic, repetitive, and quantitative investigations of drugs/chemicals. The primary objective when developing tissue/organ models is to reproduce physiologically relevant functions that typically require complex culture systems. Bioprinting offers exciting prospects for constructing 3D tissue/organ models, as it enables the reproducible, automated production of complex living tissues. Bioprinted tissues/organs may prove useful for screening novel compounds or predicting toxicity, as the spatial and chemical complexity inherent to native tissues/organs can be recreated. In this Review, we highlight the importance of developing 3D in vitro tissue/organ models by 3D bioprinting techniques, characterization of these models for evaluating their resemblance to native tissue, and their application in the prioritization of lead candidates, toxicity testing, and as disease/tumor models. PMID:26895542
3D web visualization of huge CityGML models
NASA Astrophysics Data System (ADS)
Prandi, F.; Devigili, F.; Soave, M.; Di Staso, U.; De Amicis, R.
2015-08-01
Nowadays, rapid technological development into acquiring geo-spatial information; joined to the capabilities to process these data in a relative short period of time, allows the generation of detailed 3D textured city models that will become an essential part of the modern city information infrastructure (Spatial Data Infrastructure) and, can be used to integrate various data from different sources for public accessible visualisation and many other applications. One of the main bottlenecks, which at the moment limit the use of these datasets to few experts, is a lack on efficient visualization systems through the web and interoperable frameworks that allow standardising the access to the city models. The work presented in this paper tries to satisfy these two requirements developing a 3D web-based visualization system based on OGC standards and effective visualization concepts. The architectural framework, based on Services Oriented Architecture (SOA) concepts, provides the 3D city data to a web client designed to support the view process in a very effective way. The first part of the work is to design a framework compliant to the 3D Portrayal Service drafted by the of the Open Geospatial Consortium (OGC) 3D standardization working group. The latter is related to the development of an effective web client able to render in an efficient way the 3D city models.
NASA Astrophysics Data System (ADS)
Kanaujia, Jyotima; Kumar, Ashwani; Gupta, S. C.
2016-02-01
We investigate the upper crustal velocity structure beneath the Tehri region of the Garhwal Himalaya. The investigated region is situated within the 700-km-long central seismic gap of the Himalaya that has experienced three gap-filling earthquakes since 1991 including the recent 2015 Nepal earthquake (Mw 7.8). The local tomographic inversion is based on a dataset of 1365 events collected from January 2008 to December 2012 by a 12-station local network that covers an area of about 100 × 80 km around Tehri Dam. We perform a simultaneous inversion for P- and S-wave velocity anomalies. Tomograms are interpreted in the backdrop of the regional geological and tectonic framework of the region. The spatial distribution of relocated events from the 3- D velocity model has shed new light on the pattern of seismicity in the vicinity of the Main Central thrust (MCT), and has elucidated the structure of the underthrusting Indian plate. Our model exhibits a significant negative velocity anomaly up to ˜5 per cent beneath the central part of the Garhwal Inner Lesser Himalaya, and a P-wave low velocity anomaly near the Chamoli region. The seismicity zone around the Chamoli region may be attributed to the presence of fluid filled rocks. Furthermore, an area with˜3-4 per cent positive velocity anomaly is delineated to the northwest of the Uttarkashi thrust in the vicinity of the MCT. Significant findings of the study include: a flat-ramp-flat type sub-surface geometry of the underthrusting Indian plate below the Garhwal Himalaya, high velocity images representing the trend and configuration of Delhi-Haridwar-ridge below the Sub Himalaya and Lesser Himalaya, and a seismically active zone representing geometrical asperity on the basement thrust in the vicinity of the MCT.
3-D electromagnetic modeling of wakefields in accelerator components
Poole, B.R.; Caporaso, G.J.; Ng, Wang C.; Shang, C.C.; Steich, D.
1996-09-18
We discuss the use of 3-D finite-difference time-domain (FDTD) electromagnetic codes for modeling accelerator components. Computational modeling of cylindrically symmetric structures such as induction accelerator cells has been very successful in predicting the wake potential and wake impedances of these structures, but full 3-D modeling of complex structures has been limited due to substantial computer resources required for a full 3-D model. New massively parallel 3-D time domain electromagnetic codes now under development using conforming unstructured meshes allow a substantial increase in the geometric fidelity of the structures being modeled. Development of these new codes are discussed in context of applicability to accelerator problems. Various 3-D structures are tested with an existing cubical cell FDTD code and wake impedances compared with simple analytic models for the structures; results will be used as benchmarks for testing the new time time domain codes. Structures under consideration include a stripline beam position monitor as well as circular and elliptical apertures in circular waveguides. Excellent agreement for monopole and dipole impedances with models were found for these structures below the cutoff frequency of the beam line.
3-D electromagnetic modeling of wakefields in accelerator components
Poole, Brian R.; Caporaso, George J.; Ng, Wang C.; Shang, Clifford C.; Steich, David
1997-02-01
We discuss the use of 3-D finite-difference time-domain (FDTD) electromagnetic codes for the modeling of accelerator components. Computational modeling of cylindrically symmetric structures such as induction accelerator cells has been extremely successful in predicting the wake potential and wake impedances of these structures, but fully 3-D modeling of complex structures has been limited due to the substantial computer resources required for a fully 3-D model. New massively parallel 3-D time domain electromagnetic codes now under development using conforming unstructured meshes allow a substantial increase in the geometric fidelity of the structures being modeled. Development of these new codes will be discussed in the context of their applicability to accelerator problems. A variety of 3-D structures are tested with an existing cubical cell FDTD code and the wake impedances are compared with simple analytic models for the structures. These results will provide a set of benchmarks for testing the new time domain codes. Structures under consideration include a stripline beam position monitor as well as circular and elliptical apertures in circular waveguides. Excellent agreement for the monopole and dipole impedances with the models are found for these structures below the cutoff frequency of the beam line.
Multivariate 3D modelling of Scottish soil properties
NASA Astrophysics Data System (ADS)
Poggio, Laura; Gimona, Alessandro
2015-04-01
Information regarding soil properties across landscapes at national or continental scales is critical for better soil and environmental management and for climate regulation and adaptation policy. The prediction of soil properties variation in space and time and their uncertainty is an important part of environmental modelling. Soil properties, and in particular the 3 fractions of soil texture, exhibit strong co-variation among themselves and therefore taking into account this correlation leads to spatially more accurate results. In this study the continuous vertical and lateral distributions of relevant soil properties in Scottish soils were modelled with a multivariate 3D-GAM+GS approach. The approach used involves 1) modelling the multivariate trend with full 3D spatial correlation, i.e., exploiting the values of the neighbouring pixels in 3D-space, and 2) 3D kriging to interpolate the residuals. The values at each cell for each of the considered depth layers were defined using a hybrid GAM-geostatistical 3D model, combining the fitting of a GAM (generalised Additive Models) to estimate multivariate trend of the variables, using a 3D smoother with related covariates. Gaussian simulations of the model residuals were used as spatial component to account for local details. A dataset of about 26,000 horizons (7,800 profiles) was used for this study. A validation set was randomly selected as 25% of the full dataset. Numerous covariates derived from globally available data, such as MODIS and SRTM, are considered. The results of the 3D-GAM+kriging showed low RMSE values, good R squared and an accurate reproduction of the spatial structure of the data for a range of soil properties. The results have an out-of-sample RMSE between 10 to 15% of the observed range when taking into account the whole profile. The approach followed allows the assessment of the uncertainty of both the trend and the residuals.
Perception-based shape retrieval for 3D building models
NASA Astrophysics Data System (ADS)
Zhang, Man; Zhang, Liqiang; Takis Mathiopoulos, P.; Ding, Yusi; Wang, Hao
2013-01-01
With the help of 3D search engines, a large number of 3D building models can be retrieved freely online. A serious disadvantage of most rotation-insensitive shape descriptors is their inability to distinguish between two 3D building models which are different at their main axes, but appear similar when one of them is rotated. To resolve this problem, we present a novel upright-based normalization method which not only correctly rotates such building models, but also greatly simplifies and accelerates the abstraction and the matching of building models' shape descriptors. Moreover, the abundance of architectural styles significantly hinders the effective shape retrieval of building models. Our research has shown that buildings with different designs are not well distinguished by the widely recognized shape descriptors for general 3D models. Motivated by this observation and to further improve the shape retrieval quality, a new building matching method is introduced and analyzed based on concepts found in the field of perception theory and the well-known Light Field descriptor. The resulting normalized building models are first classified using the qualitative shape descriptors of Shell and Unevenness which outline integral geometrical and topological information. These models are then put in on orderly fashion with the help of an improved quantitative shape descriptor which we will term as Horizontal Light Field Descriptor, since it assembles detailed shape characteristics. To accurately evaluate the proposed methodology, an enlarged building shape database which extends previous well-known shape benchmarks was implemented as well as a model retrieval system supporting inputs from 2D sketches and 3D models. Various experimental performance evaluation results have shown that, as compared to previous methods, retrievals employing the proposed matching methodology are faster and more consistent with human recognition of spatial objects. In addition these performance
Evaluating Biomaterial- and Microfluidic-Based 3D Tumor Models.
Carvalho, Mariana R; Lima, Daniela; Reis, Rui L; Correlo, Vitor M; Oliveira, Joaquim M
2015-11-01
Cancer is a major cause of morbidity and mortality worldwide, with a disease burden estimated to increase over the coming decades. Disease heterogeneity and limited information on cancer biology and disease mechanisms are aspects that 2D cell cultures fail to address. Here, we review the current ‘state-of-the-art’ in 3D tissue-engineering (TE) models developed for, and used in, cancer research. We assess the potential for scaffold-based TE models and microfluidics to fill the gap between 2D models and clinical application. We also discuss recent advances in combining the principles of 3D TE models and microfluidics, with a special focus on biomaterials and the most promising chip-based 3D models. PMID:26603572
Microfluidic 3D cell culture: from tools to tissue models.
van Duinen, Vincent; Trietsch, Sebastiaan J; Joore, Jos; Vulto, Paul; Hankemeier, Thomas
2015-12-01
The transition from 2D to 3D cell culture techniques is an important step in a trend towards better biomimetic tissue models. Microfluidics allows spatial control over fluids in micrometer-sized channels has become a valuable tool to further increase the physiological relevance of 3D cell culture by enabling spatially controlled co-cultures, perfusion flow and spatial control over of signaling gradients. This paper reviews most important developments in microfluidic 3D culture since 2012. Most efforts were exerted in the field of vasculature, both as a tissue on its own and as part of cancer models. We observe that the focus is shifting from tool building to implementation of specific tissue models. The next big challenge for the field is the full validation of these models and subsequently the implementation of these models in drug development pipelines of the pharmaceutical industry and ultimately in personalized medicine applications. PMID:26094109
Eck, Simon; Wörz, Stefan; Müller-Ott, Katharina; Hahn, Matthias; Biesdorf, Andreas; Schotta, Gunnar; Rippe, Karsten; Rohr, Karl
2016-08-01
The genome is partitioned into regions of euchromatin and heterochromatin. The organization of heterochromatin is important for the regulation of cellular processes such as chromosome segregation and gene silencing, and their misregulation is linked to cancer and other diseases. We present a model-based approach for automatic 3D segmentation and 3D shape analysis of heterochromatin foci from 3D confocal light microscopy images. Our approach employs a novel 3D intensity model based on spherical harmonics, which analytically describes the shape and intensities of the foci. The model parameters are determined by fitting the model to the image intensities using least-squares minimization. To characterize the 3D shape of the foci, we exploit the computed spherical harmonics coefficients and determine a shape descriptor. We applied our approach to 3D synthetic image data as well as real 3D static and real 3D time-lapse microscopy images, and compared the performance with that of previous approaches. It turned out that our approach yields accurate 3D segmentation results and performs better than previous approaches. We also show that our approach can be used for quantifying 3D shape differences of heterochromatin foci. PMID:27037463
Adhikari, Deepak; Gemmell, Brad J; Hallberg, Michael P; Longmire, Ellen K; Buskey, Edward J
2015-11-01
We describe an automated, volumetric particle image velocimetry (PIV) and tracking method that measures time-resolved, 3D zooplankton trajectories and surrounding volumetric fluid velocity fields simultaneously and non-intrusively. The method is demonstrated for groups of copepods flowing past a wall-mounted cylinder. We show that copepods execute escape responses when subjected to a strain rate threshold upstream of a cylinder, but the same threshold range elicits no escape responses in the turbulent wake downstream. The method was also used to document the instantaneous slip velocity of zooplankton and the resulting differences in trajectory between zooplankton and non-inertial fluid particles in the unsteady wake flow, showing the method's capability to quantify drift for both passive and motile organisms in turbulent environments. Applications of the method extend to any group of organisms interacting with the surrounding fluid environment, where organism location, larger-scale eddies and smaller-scale fluid deformation rates can all be tracked and analyzed. PMID:26486364
NASA Astrophysics Data System (ADS)
Harrington, R. M.; Kwiatek, G.; Moran, S. C.
2014-12-01
Shallow low frequency seismic events are common features associated with restless and erupting volcanoes. The physical mechanisms generating their characteristic low frequency, and often extended duration signals remain poorly understood. Here we present new double-difference relocations and spectral scaling of a group of ~400 shallow low-frequency seismic events occurring within the Mount St. Helens edifice during its 2004-2008 dome-building eruption, as recorded by a temporary seismic array for a month within the crater in 2006. Relocation results suggest that the majority of earthquakes occurred in the center of the crater close to the vent at depths < 500 m, with some events potentially locating under the new dome but ~200-300m southwest of the vent. Low-frequency events exhibit moment-corner frequency scaling roughly consistent with a constant static stress-drop, similar to tectonic earthquakes occurring elsewhere in shallow crustal faults. The scaling suggests that the ~400 events result from stick-slip behavior, and that the low frequency character of the waveforms may result from a combination of path effects and slow rupture speeds. For relocation, we divide the 400 events into eight families based on waveform similarity, and use a subset of nearly 40 earthquakes with clear first arrivals ranging in moment magnitude from 0.4 - 1.8 to calculate hypocenters. We then relocate these events using a double-difference method with a three-dimensional velocity model of the edifice from Waite and Moran (2009). The relocated events are then used to estimate source parameters of the remaining earthquakes via a spectral ratio technique. Spectral corner frequency estimations based on this spectral ratio approach produce stress-drop values of ~1 MPa assuming a shear-wave velocity of 1500 m/s. The estimations also indicate a constant stress-drop scaling for all events, with two event families having lower estimated stress drops of ~0.1 MPa. While localized lithological
NASA Astrophysics Data System (ADS)
Mousavi, Sima; Bauer, Klaus; Korn, Michael
2014-05-01
3-D Vp and Vp/Vs structure of the geodynamically active NW Bohemia/Vogtland area, located at the border region between Germany and Czech republic, has been determined from local earthquake tomography using 543 earthquakes which have been recorded during 2000 to 2010. This region is known for the occurrence of earthquake swarms that are supposed to be triggered by fluid upwelling in the crust, although fluid behaviour and migration paths in the subsurface of NW Bohemia is still poorly known. The events used in this study were selected based on a minimum 12 P and S phase observations and an azimuthal gap less than 160º. This data set is employed to derive a minimum 1-D velocity model and to relocate the hypocenters. The minimum 1-D velocity model is then used as an initial model in non-linear inversion to derive 3-D P-velocity and Vp/Vs ratio. Using synthetic tests, it can be shown that a high resolution is obtained in the central part of the studied region with the given source and receiver configuration. Two branches of high Vp/Vs ratio anomalies have been detected above the swarm quakes' focal zone. These anomalies support the existence of two main fluid passages toward Bad Brambach and Bublak moffette. Another interesting result is a high Vp/Vs line-like anomaly along Mariánské Lázně fault where most of the swarm quakes occur, which could be due to a fluid saturated area around the cracked zone of the fault plain. Hypocenters in the swarm region are located in a low Vp and Vp/Vs anomaly. The correlation between the detected Vp and Vp/Vs anomalies and the location of earthquake swarm suggests a model in which CO2 as part of magmatic fluids exist in a vast area beneath NW Bohemia and frequently migrate up to the surface.
NASA Astrophysics Data System (ADS)
Toloui, Mostafa; Brajkovic, David; Hong, Jiarong
2014-11-01
Digital In-line Holography is employed to image 3D flow structures in the vicinity of a transparent rough surface consisting of closely packed biomimetic shark skin denticles as roughness elements. The 3D printed surface replicates the morphological features of real shark skin, and the denticles have a geometrical scale of 2 mm, i.e. 10 times of the real ones. In order to minimize optical aberrations near the fluid-roughness interface and enable flow measurements around denticles, the optical refractive index of the fluid medium is maintained the same as that of the denticle model in an index-matched flow facility using NaI solution as the working fluid. The experiment is conducted in a 1.2 m long test section with 50 mm × 50 mm cross section. The sampling volume is located in the downstream region of a shark skin replica of 12'' stretch where the turbulent flow is fully-developed and the transitional effect from smooth to the rough surface becomes negligible. Several instantaneous realizations of the 3D velocity field are obtained and are used to illustrate turbulent coherent structures induced by shark-skin denticles. This information will provide insights on the hydrodynamic function of shark's unique surface ornamentation.
NASA Astrophysics Data System (ADS)
Ghofrani Tabari, M.; Goodfellow, S. D.; Nasseri, M. B.; Young, R.
2013-12-01
A cubic specimen of water saturated Fontainebleau Sandstone is tested in the laboratory under true-triaxial loading where three different principal stresses are applied under drained conditions. Due to the loading arrangement, closure and opening of the pre-existing cracks in the rock, as well as creation and growth of the aligned cracks cause elliptical anisotropy and distributed heterogeneities. A Geophysical Imaging Cell equipped with an Acoustic Emission monitoring system is employed to image velocity structure of the sample during the experiment through repeated transducer to transducer non-destructive ultrasonic surveys. Apparent P-wave velocities along the rock body are calculated in different directions and shown in stereonet plots which demonstrate an overall anisotropy of the sample. The apparent velocities in the main three orthogonal cubic directions are used as raw data for building a mean spatial distribution model of anisotropy ratios. This approach is based on the concept of semi-principal axes in an elliptical anisotropic model and appointing two ratios between the three orthogonal velocities in each of the cubic grid cells. The spatial distribution model of anisotropy ratios are used to calculate the anisotropic ray-path segment matrix elements (Gij). These contain segment lengths of the ith ray in the jth cell in three dimensions where, length of each ray in each cell is computed for one principal direction based on the dip and strike of the ray and these lengths differ from the ones in an isotropic G Matrix. 3D strain of the squeezed rock and the consequent geometrical deformation is also included in the ray-path segment matrix. A Singular Value Decomposition (SVD) method is used for inversion from the data space of apparent velocities to the model space of P-wave propagation velocities in the three principal directions. Finally, spatial variation and temporal evolution of induced damages in the rock, representing uniformly distributed or
NASA Astrophysics Data System (ADS)
Uznir, U.; Anton, F.; Suhaibah, A.; Rahman, A. A.; Mioc, D.
2013-09-01
The advantages of three dimensional (3D) city models can be seen in various applications including photogrammetry, urban and regional planning, computer games, etc.. They expand the visualization and analysis capabilities of Geographic Information Systems on cities, and they can be developed using web standards. However, these 3D city models consume much more storage compared to two dimensional (2D) spatial data. They involve extra geometrical and topological information together with semantic data. Without a proper spatial data clustering method and its corresponding spatial data access method, retrieving portions of and especially searching these 3D city models, will not be done optimally. Even though current developments are based on an open data model allotted by the Open Geospatial Consortium (OGC) called CityGML, its XML-based structure makes it challenging to cluster the 3D urban objects. In this research, we propose an opponent data constellation technique of space-filling curves (3D Hilbert curves) for 3D city model data representation. Unlike previous methods, that try to project 3D or n-dimensional data down to 2D or 3D using Principal Component Analysis (PCA) or Hilbert mappings, in this research, we extend the Hilbert space-filling curve to one higher dimension for 3D city model data implementations. The query performance was tested using a CityGML dataset of 1,000 building blocks and the results are presented in this paper. The advantages of implementing space-filling curves in 3D city modeling will improve data retrieval time by means of optimized 3D adjacency, nearest neighbor information and 3D indexing. The Hilbert mapping, which maps a subinterval of the [0, 1] interval to the corresponding portion of the d-dimensional Hilbert's curve, preserves the Lebesgue measure and is Lipschitz continuous. Depending on the applications, several alternatives are possible in order to cluster spatial data together in the third dimension compared to its
3D modeling of carbonates petro-acoustic heterogeneities
NASA Astrophysics Data System (ADS)
Baden, Dawin; Guglielmi, Yves; Saracco, Ginette; Marié, Lionel; Viseur, Sophie
2015-04-01
Characterizing carbonate reservoirs heterogeneity is a challenging issue for Oil & Gas Industry, CO2 sequestration and all kinds of fluid manipulations in natural reservoirs, due to the significant impact of heterogeneities on fluid flow and storage within the reservoir. Although large scale (> meter) heterogeneities such as layers petrophysical contrasts are well addressed by computing facies-based models, low scale (< meter) heterogeneities are often poorly constrained because of the complexity in predicting their spatial arrangement. In this study, we conducted petro-acoustic measurements on cores of different size and diameter (Ø = 1", 1.5" and 5") in order to evaluate anisotropy or heterogeneity in carbonates at different laboratory scales. Different types of heterogeneities which generally occur in carbonate reservoir units (e.g. petrographic, diagenetic, and tectonic related) were sampled. Dry / wet samples were investigated with different ultrasonic apparatus and using different sensors allowing acoustic characterization through a bandwidth varying from 50 to 500 kHz. Comprehensive measurements realized on each samples allowed statistical analyses of petro-acoustic properties such as attenuation, shear and longitudinal wave velocity. The cores properties (geological and acoustic facies) were modeled in 3D using photogrammetry and GOCAD geo-modeler. This method successfully allowed detecting and imaging in three dimensions differential diagenesis effects characterized by the occurrence of decimeter-scale diagenetic horizons in samples assumed to be homogeneous and/or different diagenetic sequences between shells filling and the packing matrix. We then discuss how small interfaces such as cracks, stylolithes and laminations which are also imaged may have guided these differential effects, considering that understanding the processes may be taken as an analogue to actual fluid drainage complexity in deep carbonate reservoir.
MR image denoising method for brain surface 3D modeling
NASA Astrophysics Data System (ADS)
Zhao, De-xin; Liu, Peng-jie; Zhang, De-gan
2014-11-01
Three-dimensional (3D) modeling of medical images is a critical part of surgical simulation. In this paper, we focus on the magnetic resonance (MR) images denoising for brain modeling reconstruction, and exploit a practical solution. We attempt to remove the noise existing in the MR imaging signal and preserve the image characteristics. A wavelet-based adaptive curve shrinkage function is presented in spherical coordinates system. The comparative experiments show that the denoising method can preserve better image details and enhance the coefficients of contours. Using these denoised images, the brain 3D visualization is given through surface triangle mesh model, which demonstrates the effectiveness of the proposed method.
STELLOPT Modeling of the 3D Diagnostic Response in ITER
Lazerson, Samuel A
2013-05-07
The ITER three dimensional diagnostic response to an n=3 resonant magnetic perturbation is modeled using the STELLOPT code. The in-vessel coils apply a resonant magnetic perturbation (RMP) fi eld which generates a 4 cm edge displacement from axisymmetry as modeled by the VMEC 3D equilibrium code. Forward modeling of flux loop and magnetic probe response with the DIAGNO code indicates up to 20 % changes in measured plasma signals. Simulated LIDAR measurements of electron temperature indicate 2 cm shifts on the low field side of the plasma. This suggests that the ITER diagnostic will be able to diagnose the 3D structure of the equilibria.
Potential of 3D City Models to assess flood vulnerability
NASA Astrophysics Data System (ADS)
Schröter, Kai; Bochow, Mathias; Schüttig, Martin; Nagel, Claus; Ross, Lutz; Kreibich, Heidi
2016-04-01
Vulnerability, as the product of exposure and susceptibility, is a key factor of the flood risk equation. Furthermore, the estimation of flood loss is very sensitive to the choice of the vulnerability model. Still, in contrast to elaborate hazard simulations, vulnerability is often considered in a simplified manner concerning the spatial resolution and geo-location of exposed objects as well as the susceptibility of these objects at risk. Usually, area specific potential flood loss is quantified on the level of aggregated land-use classes, and both hazard intensity and resistance characteristics of affected objects are represented in highly simplified terms. We investigate the potential of 3D City Models and spatial features derived from remote sensing data to improve the differentiation of vulnerability in flood risk assessment. 3D City Models are based on CityGML, an application scheme of the Geography Markup Language (GML), which represents the 3D geometry, 3D topology, semantics and appearance of objects on different levels of detail. As such, 3D City Models offer detailed spatial information which is useful to describe the exposure and to characterize the susceptibility of residential buildings at risk. This information is further consolidated with spatial features of the building stock derived from remote sensing data. Using this database a spatially detailed flood vulnerability model is developed by means of data-mining. Empirical flood damage data are used to derive and to validate flood susceptibility models for individual objects. We present first results from a prototype application in the city of Dresden, Germany. The vulnerability modeling based on 3D City Models and remote sensing data is compared i) to the generally accepted good engineering practice based on area specific loss potential and ii) to a highly detailed representation of flood vulnerability based on a building typology using urban structure types. Comparisons are drawn in terms of
3D head model classification using optimized EGI
NASA Astrophysics Data System (ADS)
Tong, Xin; Wong, Hau-san; Ma, Bo
2006-02-01
With the general availability of 3D digitizers and scanners, 3D graphical models have been used widely in a variety of applications. This has led to the development of search engines for 3D models. Especially, 3D head model classification and retrieval have received more and more attention in view of their many potential applications in criminal identifications, computer animation, movie industry and medical industry. This paper addresses the 3D head model classification problem using 2D subspace analysis methods such as 2D principal component analysis (2D PCA[3]) and 2D fisher discriminant analysis (2DLDA[5]). It takes advantage of the fact that the histogram is a 2D image, and we can extract the most useful information from these 2D images to get a good result accordingingly. As a result, there are two main advantages: First, we can perform less calculation to obtain the same rate of classification; second, we can reduce the dimensionality more than PCA to obtain a higher efficiency.
Numerical Results of Earth's Core Accumulation 3-D Modelling
NASA Astrophysics Data System (ADS)
Khachay, Yurie; Anfilogov, Vsevolod
2013-04-01
For a long time as a most convenient had been the model of mega impact in which the early forming of the Earth's core and mantle had been the consequence of formed protoplanet collision with the body of Mercurial mass. But all dynamical models of the Earth's accumulation and the estimations after the Pb-Pb system, lead to the conclusion that the duration of the planet accumulation was about 1 milliard years. But isotopic results after the W-Hf system testify about a very early (5-10) million years, dividing of the geochemical reservoirs of the core and mantle. In [1,3] it is shown, that the account of energy dissipating by the decay of short living radioactive elements and first of all Al,it is sufficient for heating even small bodies with dimensions about (50-100) km up to the iron melting temperature and can be realized a principal new differentiation mechanism. The inner parts of the melted preplanets can join and they are mainly of iron content, but the cold silicate fragments return to the supply zone. Only after the increasing of the gravitational radius, the growing area of the future core can save also the silicate envelope fragments. All existing dynamical accumulation models are constructed by using a spherical-symmetrical model. Hence for understanding the further planet evolution it is significant to trace the origin and evolution of heterogeneities, which occur on the planet accumulation stage. In that paper we are modeling distributions of temperature, pressure, velocity of matter flowing in a block of 3D- spherical body with a growing radius. The boundary problem is solved by the finite-difference method for the system of equations, which include equations which describe the process of accumulation, the Safronov equation, the equation of impulse balance, equation Navier-Stocks, equation for above litho static pressure and heat conductivity in velocity-pressure variables using the Businesque approach. The numerical algorithm of the problem solution in
Vhrs Stereo Images for 3d Modelling of Buildings
NASA Astrophysics Data System (ADS)
Bujakiewicz, A.; Holc, M.
2012-07-01
The paper presents the project which was carried out in the Photogrammetric Laboratory of Warsaw University of Technology. The experiment is concerned with the extraction of 3D vector data for buildings creation from 3D photogrammetric model based on the Ikonos stereo images. The model was reconstructed with photogrammetric workstation - Summit Evolution combined with ArcGIS 3D platform. Accuracy of 3D model was significantly improved by use for orientation of pair of satellite images the stereo measured tie points distributed uniformly around the model area in addition to 5 control points. The RMS for model reconstructed on base of the RPC coefficients only were 16,6 m, 2,7 m and 47,4 m, for X, Y and Z coordinates, respectively. By addition of 5 control points the RMS were improved to 0,7 m, 0,7 m 1,0 m, where the best results were achieved when RMS were estimated from deviations in 17 check points (with 5 control points)and amounted to 0,4 m, 0,5 m and 0,6 m, for X, Y, and Z respectively. The extracted 3D vector data for buildings were integrated with 2D data of the ground footprints and afterwards they were used for 3D modelling of buildings in Google SketchUp software. The final results were compared with the reference data obtained from other sources. It was found that the shape of buildings (in concern to the number of details) had been reconstructed on level of LoD1, when the accuracy of these models corresponded to the level of LoD2.
3-D Velocity Structure of Southwestern British Columbia and Northern Washington
NASA Astrophysics Data System (ADS)
Ramachandran, K.; Ramachandran, K.; Spence, G. D.; Dosso, S.; Hyndman, R. D.; Hyndman, R. D.; Brocher, T. M.; Fisher, M. M.
2001-12-01
A seismic tomography analysis in S.W. British Columbia and N. Washington has been used to define the velocity structure of the forearc crust and underlying subducting Juan de Fuca plate, and to obtain precise earthquake locations. First arrival travel-times from earthquakes and from the large airgun array used in the `Seismic Hazards Investigation of Puget Sound' (SHIPS) 1998 experiment, were simultaneously inverted for hypocentral parameters and velocity structure. Approximately 16,000 picks from 1,400 earthquakes recorded at 46 permanent stations, and 35,000 picks from the SHIPS experiment were used in the inversion. The velocity model was parameterized in the forward/inverse step by a node/cell spacing of 3 X 3 X 3 km over a volume of 360 X 450 X 93 km depth. The starting and final RMS travel time misfits were 479 ms and 120 ms respectively. Checkerboard tests conducted on the final velocity model imply good lateral resolution ranging from 30 to 50 km. The SHIPS airgun data mainly constrained the upper ~12 km and the earthquake data the deeper structure. The high velocity mafic Crescent Terrane that dips beneath the margin is well mapped in the velocity model on a regional scale. Its thickness beneath southern Vancouver Island is interpreted to reach ~20 km. Three high velocity structures above the subducting Juan de Fuca plate, having mafic to ultramafic velocities of 7.25-7.5 km/s, occur beneath southern Vancouver Island and Puget Sound at a depth of ~25 km. They may be associated with deeper parts of the Crescent Terrane, or with structures such as seamounts on the subducting Juan de Fuca plate. At the southern tip of Vancouver Islands, the Leech River Fault, Southern Whidbey Island Fault, and the Devils Mountain Fault appear to correlate with mapped seismicity. The subducting Juan de Fuca plate is well mapped beneath southern Vancouver Island, Olympic Peninsula, Strait of Georgia, and Puget Sound. The velocity model identifies the steepening dip in the
3D model of amphioxus steroid receptor complexed with estradiol
Baker, Michael E.; Chang, David J.
2009-08-28
The origins of signaling by vertebrate steroids are not fully understood. An important advance was the report that an estrogen-binding steroid receptor [SR] is present in amphioxus, a basal chordate with a similar body plan as vertebrates. To investigate the evolution of estrogen-binding to steroid receptors, we constructed a 3D model of amphioxus SR complexed with estradiol. This 3D model indicates that although the SR is activated by estradiol, some interactions between estradiol and human ER{alpha} are not conserved in the SR, which can explain the low affinity of estradiol for the SR. These differences between the SR and ER{alpha} in the steroid-binding domain are sufficient to suggest that another steroid is the physiological regulator of the SR. The 3D model predicts that mutation of Glu-346 to Gln will increase the affinity of testosterone for amphioxus SR and elucidate the evolution of steroid-binding to nuclear receptors.
Quality of 3D Models Generated by SFM Technology
NASA Astrophysics Data System (ADS)
Marčiš, Marián
2013-12-01
Using various types of automation in digital photogrammetry is associated with questions such as the accuracy of a 3D model generated on various types of surfaces and textures, the financial costs of the equipment needed, and also the time costs of the processing. This topic deals with the actual technology of computer vision, which allows the automated exterior orientation of images, camera calibration, and the generation of 3D models directly from images of the object itself, based on the automatic detection of significant points. Detailed testing is done using the Agisoft PhotoScan system, and the camera configuration is solved with respect to the accuracy of the 3D model generated and the time consumption of the calculations for the different types of textures and the different settings for the processing.
Diffusion approximation for modeling of 3-D radiation distributions
Zardecki, A.; Gerstl, S.A.W.; De Kinder, R.E. Jr.
1985-01-01
A three-dimensional transport code DIF3D, based on the diffusion approximation, is used to model the spatial distribution of radiation energy arising from volumetric isotropic sources. Future work will be concerned with the determination of irradiances and modeling of realistic scenarios, relevant to the battlefield conditions. 8 refs., 4 figs.
3D surface digitizing and modeling development at ITRI
NASA Astrophysics Data System (ADS)
Hsueh, Wen-Jean
2000-06-01
This paper gives an overview of the research and development activities in 3D surface digitizing and modeling conducted at the Industrial Technology Research Institute (ITRI) of Taiwan in the past decade. As a major technology and consulting service provider of the area, ITRI has developed 3D laser scanning digitizers ranging from low-cost compacts, industrial CAD/CAM digitizing, to large human body scanner, with in-house 3D surface modeling software to provide total solution in reverse engineering that requires processing capabilities of large number of 3D data. Based on both hardware and software technologies in scanning, merging, registration, surface fitting, reconstruction, and compression, ITRI is now exploring innovative methodologies that provide higher performances, including hardware-based correlation algorithms with advanced camera designs, animation surface model reconstruction, and optical tracking for motion capture. It is expected that the need for easy and fast high-quality 3D information in the near future will grow exponentially, at the same amazing rate as the internet and the human desire for realistic and natural images.
Modelling Gaia CCD pixels with Silvaco 3D engineering software
NASA Astrophysics Data System (ADS)
Seabroke, G. M.; Prod'Homme, T.; Hopkinson, G.; Burt, D.; Robbins, M.; Holland, A.
2011-02-01
Gaia will only achieve its unprecedented measurement accuracy requirements with detailed calibration and correction for radiation damage. We present our Silvaco 3D engineering software model of the Gaia CCD pixel and two of its applications for Gaia: (1) physically interpreting supplementary buried channel (SBC) capacity measurements (pocket-pumping and first pixel response) in terms of e2v manufacturing doping alignment tolerances; and (2) deriving electron densities within a charge packet as a function of the number of constituent electrons and 3D position within the charge packet as input to microscopic models being developed to simulate radiation damage.
Parallel Optimization of 3D Cardiac Electrophysiological Model Using GPU
Xia, Yong; Wang, Kuanquan; Zhang, Henggui
2015-01-01
Large-scale 3D virtual heart model simulations are highly demanding in computational resources. This imposes a big challenge to the traditional computation resources based on CPU environment, which already cannot meet the requirement of the whole computation demands or are not easily available due to expensive costs. GPU as a parallel computing environment therefore provides an alternative to solve the large-scale computational problems of whole heart modeling. In this study, using a 3D sheep atrial model as a test bed, we developed a GPU-based simulation algorithm to simulate the conduction of electrical excitation waves in the 3D atria. In the GPU algorithm, a multicellular tissue model was split into two components: one is the single cell model (ordinary differential equation) and the other is the diffusion term of the monodomain model (partial differential equation). Such a decoupling enabled realization of the GPU parallel algorithm. Furthermore, several optimization strategies were proposed based on the features of the virtual heart model, which enabled a 200-fold speedup as compared to a CPU implementation. In conclusion, an optimized GPU algorithm has been developed that provides an economic and powerful platform for 3D whole heart simulations. PMID:26581957
Geospatial Modelling Approach for 3d Urban Densification Developments
NASA Astrophysics Data System (ADS)
Koziatek, O.; Dragićević, S.; Li, S.
2016-06-01
With growing populations, economic pressures, and the need for sustainable practices, many urban regions are rapidly densifying developments in the vertical built dimension with mid- and high-rise buildings. The location of these buildings can be projected based on key factors that are attractive to urban planners, developers, and potential buyers. Current research in this area includes various modelling approaches, such as cellular automata and agent-based modelling, but the results are mostly linked to raster grids as the smallest spatial units that operate in two spatial dimensions. Therefore, the objective of this research is to develop a geospatial model that operates on irregular spatial tessellations to model mid- and high-rise buildings in three spatial dimensions (3D). The proposed model is based on the integration of GIS, fuzzy multi-criteria evaluation (MCE), and 3D GIS-based procedural modelling. Part of the City of Surrey, within the Metro Vancouver Region, Canada, has been used to present the simulations of the generated 3D building objects. The proposed 3D modelling approach was developed using ESRI's CityEngine software and the Computer Generated Architecture (CGA) language.
3D Model Generation From the Engineering Drawing
NASA Astrophysics Data System (ADS)
Vaský, Jozef; Eliáš, Michal; Bezák, Pavol; Červeňanská, Zuzana; Izakovič, Ladislav
2010-01-01
The contribution deals with the transformation of engineering drawings in a paper form into a 3D computer representation. A 3D computer model can be further processed in CAD/CAM system, it can be modified, archived, and a technical drawing can be then generated from it as well. The transformation process from paper form to the data one is a complex and difficult one, particularly owing to the different types of drawings, forms of displayed objects and encountered errors and deviations from technical standards. The algorithm for 3D model generating from an orthogonal vector input representing a simplified technical drawing of the rotational part is described in this contribution. The algorithm was experimentally implemented as ObjectARX application in the AutoCAD system and the test sample as the representation of the rotational part was used for verificaton.
Robust model-based 3d/3D fusion using sparse matching for minimally invasive surgery.
Neumann, Dominik; Grbic, Sasa; John, Matthias; Navab, Nassir; Hornegger, Joachim; Ionasec, Razvan
2013-01-01
Classical surgery is being disrupted by minimally invasive and transcatheter procedures. As there is no direct view or access to the affected anatomy, advanced imaging techniques such as 3D C-arm CT and C-arm fluoroscopy are routinely used for intra-operative guidance. However, intra-operative modalities have limited image quality of the soft tissue and a reliable assessment of the cardiac anatomy can only be made by injecting contrast agent, which is harmful to the patient and requires complex acquisition protocols. We propose a novel sparse matching approach for fusing high quality pre-operative CT and non-contrasted, non-gated intra-operative C-arm CT by utilizing robust machine learning and numerical optimization techniques. Thus, high-quality patient-specific models can be extracted from the pre-operative CT and mapped to the intra-operative imaging environment to guide minimally invasive procedures. Extensive quantitative experiments demonstrate that our model-based fusion approach has an average execution time of 2.9 s, while the accuracy lies within expert user confidence intervals. PMID:24505663
Test Problems for Reactive Flow HE Model in the ALE3D Code and Limited Sensitivity Study
Gerassimenko, M.
2000-03-01
We document quick running test problems for a reactive flow model of HE initiation incorporated into ALE3D. A quarter percent change in projectile velocity changes the outcome from detonation to HE burn that dies down. We study the sensitivity of calculated HE behavior to several parameters of practical interest where modeling HE initiation with ALE3D.
Geodiversity: Exploration of 3D geological model space
NASA Astrophysics Data System (ADS)
Lindsay, M. D.; Jessell, M. W.; Ailleres, L.; Perrouty, S.; de Kemp, E.; Betts, P. G.
2013-05-01
The process of building a 3D model necessitates the reconciliation of field observations, geophysical interpretation, geological data uncertainty and the prevailing tectonic evolution hypotheses and interpretations. Uncertainty is compounded when clustered data points collected at local scales are statistically upscaled to one or two points for use in regional models. Interpretation is required to interpolate between sparse field data points using ambiguous geophysical data in covered terranes. It becomes clear that multiple interpretations are possible during model construction. The various interpretations are considered as potential natural representatives, but pragmatism typically dictates that just a single interpretation is offered by the modelling process. Uncertainties are introduced into the 3D model during construction from a variety of sources and through data set optimisation that produces a single model. Practices such as these are likely to result in a model that does not adequately represent the target geology. A set of geometrical ‘geodiversity’ metrics are used to analyse a 3D model of the Gippsland Basin, southeastern Australia after perturbing geological input data via uncertainty simulation. The resulting sets of perturbed geological observations are used to calculate a suite of geological 3D models that display a range of geological architectures. The concept of biodiversity has been adapted for the geosciences to quantify geometric variability, or geodiversity, between models in order to understand the effect uncertainty has models geometry. Various geometrical relationships (depth, volume, contact surface area, curvature and geological complexity) are used to describe the range of possibilities exhibited throughout the model suite. End-member models geodiversity metrics are classified in a similar manner to taxonomic descriptions. Further analysis of the model suite is performed using principal component analysis (PCA) to determine
3D shape decomposition and comparison for gallbladder modeling
NASA Astrophysics Data System (ADS)
Huang, Weimin; Zhou, Jiayin; Liu, Jiang; Zhang, Jing; Yang, Tao; Su, Yi; Law, Gim Han; Chui, Chee Kong; Chang, Stephen
2011-03-01
This paper presents an approach to gallbladder shape comparison by using 3D shape modeling and decomposition. The gallbladder models can be used for shape anomaly analysis and model comparison and selection in image guided robotic surgical training, especially for laparoscopic cholecystectomy simulation. The 3D shape of a gallbladder is first represented as a surface model, reconstructed from the contours segmented in CT data by a scheme of propagation based voxel learning and classification. To better extract the shape feature, the surface mesh is further down-sampled by a decimation filter and smoothed by a Taubin algorithm, followed by applying an advancing front algorithm to further enhance the regularity of the mesh. Multi-scale curvatures are then computed on the regularized mesh for the robust saliency landmark localization on the surface. The shape decomposition is proposed based on the saliency landmarks and the concavity, measured by the distance from the surface point to the convex hull. With a given tolerance the 3D shape can be decomposed and represented as 3D ellipsoids, which reveal the shape topology and anomaly of a gallbladder. The features based on the decomposed shape model are proposed for gallbladder shape comparison, which can be used for new model selection. We have collected 19 sets of abdominal CT scan data with gallbladders, some shown in normal shape and some in abnormal shapes. The experiments have shown that the decomposed shapes reveal important topology features.
NASA Astrophysics Data System (ADS)
Sharkawi, K.-H.; Abdul-Rahman, A.
2013-09-01
Cities and urban areas entities such as building structures are becoming more complex as the modern human civilizations continue to evolve. The ability to plan and manage every territory especially the urban areas is very important to every government in the world. Planning and managing cities and urban areas based on printed maps and 2D data are getting insufficient and inefficient to cope with the complexity of the new developments in big cities. The emergence of 3D city models have boosted the efficiency in analysing and managing urban areas as the 3D data are proven to represent the real world object more accurately. It has since been adopted as the new trend in buildings and urban management and planning applications. Nowadays, many countries around the world have been generating virtual 3D representation of their major cities. The growing interest in improving the usability of 3D city models has resulted in the development of various tools for analysis based on the 3D city models. Today, 3D city models are generated for various purposes such as for tourism, location-based services, disaster management and urban planning. Meanwhile, modelling 3D objects are getting easier with the emergence of the user-friendly tools for 3D modelling available in the market. Generating 3D buildings with high accuracy also has become easier with the availability of airborne Lidar and terrestrial laser scanning equipments. The availability and accessibility to this technology makes it more sensible to analyse buildings in urban areas using 3D data as it accurately represent the real world objects. The Open Geospatial Consortium (OGC) has accepted CityGML specifications as one of the international standards for representing and exchanging spatial data, making it easier to visualize, store and manage 3D city models data efficiently. CityGML able to represents the semantics, geometry, topology and appearance of 3D city models in five well-defined Level-of-Details (LoD), namely LoD0
NASA Astrophysics Data System (ADS)
Yu, H.-S.; Jackson, B. V.; Hick, P. P.; Buffington, A.; Odstrcil, D.; Wu, C.-C.; Davies, J. A.; Bisi, M. M.; Tokumaru, M.
2015-09-01
The University of California, San Diego, time-dependent analyses of the heliosphere provide three-dimensional (3D) reconstructions of solar wind velocities and densities from observations of interplanetary scintillation (IPS). Using data from the Solar-Terrestrial Environment Laboratory, Japan, these reconstructions provide a real-time prediction of the global solar-wind density and velocity throughout the whole heliosphere with a temporal cadence of about one day (ips.ucsd.edu). Updates to this modeling effort continue: in the present article, near-Sun results extracted from the time-dependent 3D reconstruction are used as inner boundary conditions to drive 3D-MHD models ( e.g. ENLIL and H3D-MHD). This allows us to explore the differences between the IPS kinematic-model data-fitting procedure and current 3D-MHD modeling techniques. The differences in these techniques provide interesting insights into the physical principles governing the expulsion of coronal mass ejections (CMEs). Here we detail for the first time several specific CMEs and an induced shock that occurred in September 2011 that demonstrate some of the issues resulting from these analyses.
Creating Physical 3D Stereolithograph Models of Brain and Skull
Kelley, Daniel J.; Farhoud, Mohammed; Meyerand, M. Elizabeth; Nelson, David L.; Ramirez, Lincoln F.; Dempsey, Robert J.; Wolf, Alan J.; Alexander, Andrew L.; Davidson, Richard J.
2007-01-01
The human brain and skull are three dimensional (3D) anatomical structures with complex surfaces. However, medical images are often two dimensional (2D) and provide incomplete visualization of structural morphology. To overcome this loss in dimension, we developed and validated a freely available, semi-automated pathway to build 3D virtual reality (VR) and hand-held, stereolithograph models. To evaluate whether surface visualization in 3D was more informative than in 2D, undergraduate students (n = 50) used the Gillespie scale to rate 3D VR and physical models of both a living patient-volunteer's brain and the skull of Phineas Gage, a historically famous railroad worker whose misfortune with a projectile tamping iron provided the first evidence of a structure-function relationship in brain. Using our processing pathway, we successfully fabricated human brain and skull replicas and validated that the stereolithograph model preserved the scale of the VR model. Based on the Gillespie ratings, students indicated that the biological utility and quality of visual information at the surface of VR and stereolithograph models were greater than the 2D images from which they were derived. The method we developed is useful to create VR and stereolithograph 3D models from medical images and can be used to model hard or soft tissue in living or preserved specimens. Compared to 2D images, VR and stereolithograph models provide an extra dimension that enhances both the quality of visual information and utility of surface visualization in neuroscience and medicine. PMID:17971879
Energy cost and body centre of mass' 3D intracycle velocity variation in swimming.
Figueiredo, Pedro; Barbosa, Tiago M; Vilas-Boas, João Paulo; Fernandes, Ricardo J
2012-09-01
The purpose of this study was to examine the relationship between the energy cost (C) and the 3D intracycle velocity variation (IVV; swimming direction--x, vertical--y and lateral--z axes) throughout the 200 m front crawl event. Ten international level swimmers performed a maximal 200 m front crawl swim followed by 50, 100 and 150 m bouts at the same pace as in the 200 m splits. Oxygen consumption was measured during the bouts and blood samples were collected before and after each one. The C was calculated for each 50 m lap as the ratio of the total energy expenditure (three energy pathways) to the distance. A respiratory snorkel and valve system with low hydrodynamic resistance was used to measure pulmonary ventilation and to collect breathing air samples. Two above water and four underwater cameras videotaped the swim bouts and thereafter APAS was used to assess the centre of mass IVV (x, y and z components). The increase in the C was significantly associated with the increase in the IVV in x for the first 50 m lap (R = -0.83, P < 0.01). It is concluded that the IVV relationship with C in a competitive event does not present the direct relationship found in the literature, revealing a great specificity, which suggests that the relation between these two parameters could not be used as a performance predictor in competitive events. PMID:22262010
Teaching the geological subsurface with 3D models
NASA Astrophysics Data System (ADS)
Thorpe, Steve; Ward, Emma
2014-05-01
3D geological models have great potential as a resource when teaching geological concepts as it allows the student to visualise and interrogate UK geology. They are especially useful when dealing with the conversion of 2D field, map and GIS outputs into three dimensional geological units, which is a common problem for many students. Today's earth science students use a variety of skills and processes during their learning experience including spatial thinking, image construction, detecting patterns, making predictions and deducing the orientation of themselves. 3D geological models can reinforce spatial thinking strategies and encourage students to think about processes and properties, in turn helping the student to recognise pre-learnt geological principles in the field and to convert what they see at the surface into a picture of what is going on at depth. The British Geological Survey (BGS) has been producing digital 3D geological models for over 10 years. The models produced are revolutionising the working practices, data standards and products of the BGS. Sharing our geoscience information with academia is highlighted throughout the BGS strategy as is instilling practical skills in future geoscience professionals, such as model building and interpretation. In 2009 a project was launched to investigate the potential of the models as a teaching resource. The study included justifying if and how the models help students to learn, how models have been used historically, and how other forms of modelling are being used today. BGS now produce 3D geological models for use by anyone teaching or learning geoscience. They incorporate educational strategies that will develop geospatial skills and alleviate potential problems that some students experience. They are contained within contemporary case studies and show standard geological concepts, structures, sedimentary rocks, cross sections and field techniques. 3D geological models of the Isle of Wight and Ingleborough
3D Geological Model for "LUSI" - a Deep Geothermal System
NASA Astrophysics Data System (ADS)
Sohrabi, Reza; Jansen, Gunnar; Mazzini, Adriano; Galvan, Boris; Miller, Stephen A.
2016-04-01
Geothermal applications require the correct simulation of flow and heat transport processes in porous media, and many of these media, like deep volcanic hydrothermal systems, host a certain degree of fracturing. This work aims to understand the heat and fluid transport within a new-born sedimentary hosted geothermal system, termed Lusi, that began erupting in 2006 in East Java, Indonesia. Our goal is to develop conceptual and numerical models capable of simulating multiphase flow within large-scale fractured reservoirs such as the Lusi region, with fractures of arbitrary size, orientation and shape. Additionally, these models can also address a number of other applications, including Enhanced Geothermal Systems (EGS), CO2 sequestration (Carbon Capture and Storage CCS), and nuclear waste isolation. Fractured systems are ubiquitous, with a wide-range of lengths and scales, making difficult the development of a general model that can easily handle this complexity. We are developing a flexible continuum approach with an efficient, accurate numerical simulator based on an appropriate 3D geological model representing the structure of the deep geothermal reservoir. Using previous studies, borehole information and seismic data obtained in the framework of the Lusi Lab project (ERC grant n°308126), we present here the first 3D geological model of Lusi. This model is calculated using implicit 3D potential field or multi-potential fields, depending on the geological context and complexity. This method is based on geological pile containing the geological history of the area and relationship between geological bodies allowing automatic computation of intersections and volume reconstruction. Based on the 3D geological model, we developed a new mesh algorithm to create hexahedral octree meshes to transfer the structural geological information for 3D numerical simulations to quantify Thermal-Hydraulic-Mechanical-Chemical (THMC) physical processes.
Performance and Cognitive Assessment in 3-D Modeling
ERIC Educational Resources Information Center
Fahrer, Nolan E.; Ernst, Jeremy V.; Branoff, Theodore J.; Clark, Aaron C.
2011-01-01
The purpose of this study was to investigate identifiable differences between performance and cognitive assessment scores in a 3-D modeling unit of an engineering drafting course curriculum. The study aimed to provide further investigation of the need of skill-based assessments in engineering/technical graphics courses to potentially increase…
Tracking people and cars using 3D modeling and CCTV.
Edelman, Gerda; Bijhold, Jurrien
2010-10-10
The aim of this study was to find a method for the reconstruction of movements of people and cars using CCTV footage and a 3D model of the environment. A procedure is proposed, in which video streams are synchronized and displayed in a 3D model, by using virtual cameras. People and cars are represented by cylinders and boxes, which are moved in the 3D model, according to their movements as shown in the video streams. The procedure was developed and tested in an experimental setup with test persons who logged their GPS coordinates as a recording of the ground truth. Results showed that it is possible to implement this procedure and to reconstruct movements of people and cars from video recordings. The procedure was also applied to a forensic case. In this work we experienced that more situational awareness was created by the 3D model, which made it easier to track people on multiple video streams. Based on all experiences from the experimental set up and the case, recommendations are formulated for use in practice. PMID:20439141
A Sketching Interface for Freeform 3D Modeling
NASA Astrophysics Data System (ADS)
Igarashi, Takeo
This chapter introduces Teddy, a sketch-based modeling system to quickly and easily design freeform models such as stuffed animals and other rotund objects. The user draws several 2D freeform strokes interactively on the screen and the system automatically constructs plausible 3D polygonal surfaces. Our system supports several modeling operations, including the operation to construct a 3D polygonal surface from a 2D silhouette drawn by the user: it inflates the region surrounded by the silhouette making a wide area fat, and a narrow area thin. Teddy, our prototype system, is implemented as a Java program, and the mesh construction is done in real-time on a standard PC. Our informal user study showed that a first-time user masters the operations within 10 minutes, and can construct interesting 3D models within minutes. We also report the result of a case study where a high school teacher taught various 3D concepts in geography using the system.
Assessment of 3D Models Used in Contours Studies
ERIC Educational Resources Information Center
Alvarez, F. J. Ayala; Parra, E. B. Blazquez; Tubio, F. Montes
2015-01-01
This paper presents an experimental research focusing on the view of first year students. The aim is to check the quality of implementing 3D models integrated in the curriculum. We search to determine students' preference between the various means facilitated in order to understand the given subject. Students have been respondents to prove the…
3D Geologic Model of the Southern Great Basin
NASA Astrophysics Data System (ADS)
Wagoner, J. L.; Myers, S. C.
2006-12-01
We have constructed a regional 3D geologic model of the southern Great Basin, in support of a seismic wave propagation investigation of the 1993 Nonproliferation Experiment (NPE) at the Nevada Test Site (NTS). The model is centered on the NPE and spans longitude -119.5° to -112.6°, latitude 34.5° to 39.8°, and a depth from the surface to 150 km below sea level. Hence, the model includes the southern half of Nevada, as well as parts of eastern California, western Utah, and a portion of northwestern Arizona. The upper crust is constrained by geologic and geophysical studies, and the lower crust and upper mantle are constrained by geophysical studies. The upper crustal geologic units are Quaternary basin fill, Tertiary deposits, pre-Tertiary deposits, intrusive rocks, and calderas. The lower crust and upper mantle are parameterized with 8 layers, including the Moho. Detailed geologic data, including surface maps, borehole data, and geophysical surveys, were used to define the geology at the NTS. Digital geologic outcrop data were available for both Nevada and Arizona, whereas we scanned and hand digitized geologic maps for California and Utah. Published gravity data (2km spacing) were used to determine the thickness of the Cenozoic deposits and constrain the depth of the basins. The free surface is based on a 10m lateral resolution DEM at the NTS and a 90m resolution DEM elsewhere. The gross geophysical structure of the crust and upper mantle is taken from regional surface-wave studies. Variations in crustal thickness are based on receiver function analysis and a compilation of reflection/refraction studies. We used the Earthvision (Dynamic Graphics, Inc.) software to integrate the geologic and geophysical information into a model of x,y,z,p nodes, where p is an integer index representing the geologic unit. For regional seismic simulations we convert this realistic geologic model into elastic parameters. Upper crustal units are treated as seismically homogeneous