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
Vector quantization of 3-D point clouds
NASA Astrophysics Data System (ADS)
Sim, Jae-Young; Kim, Chang-Su; Lee, Sang-Uk
2005-10-01
A geometry compression algorithm for 3-D QSplat data using vector quantization (VQ) is proposed in this work. The positions of child spheres are transformed to the local coordinate system, which is determined by the parent children relationship. The coordinate transform makes child positions more compactly distributed in 3-D space, facilitating effective quantization. Moreover, we develop a constrained encoding method for sphere radii, which guarantees hole-free surface rendering at the decoder side. Simulation results show that the proposed algorithm provides a faithful rendering quality even at low bitrates.
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.
NASA Technical Reports Server (NTRS)
Papanyan, Valeri; Oshle, Edward; Adamo, Daniel
2008-01-01
Measurement of the jettisoned object departure trajectory and velocity vector in the International Space Station (ISS) reference frame is vitally important for prompt evaluation of the object s imminent orbit. We report on the first successful application of photogrammetric analysis of the ISS imagery for the prompt computation of the jettisoned object s position and velocity vectors. As post-EVA analyses examples, we present the Floating Potential Probe (FPP) and the Russian "Orlan" Space Suit jettisons, as well as the near-real-time (provided in several hours after the separation) computations of the Video Stanchion Support Assembly Flight Support Assembly (VSSA-FSA) and Early Ammonia Servicer (EAS) jettisons during the US astronauts space-walk. Standard close-range photogrammetry analysis was used during this EVA to analyze two on-board camera image sequences down-linked from the ISS. In this approach the ISS camera orientations were computed from known coordinates of several reference points on the ISS hardware. Then the position of the jettisoned object for each time-frame was computed from its image in each frame of the video-clips. In another, "quick-look" approach used in near-real time, orientation of the cameras was computed from their position (from the ISS CAD model) and operational data (pan and tilt) then location of the jettisoned object was calculated only for several frames of the two synchronized movies. Keywords: Photogrammetry, International Space Station, jettisons, image analysis.
3D reconstruction of tensors and vectors
Defrise, Michel; Gullberg, Grant T.
2005-02-17
Here we have developed formulations for the reconstruction of 3D tensor fields from planar (Radon) and line-integral (X-ray) projections of 3D vector and tensor fields. Much of the motivation for this work is the potential application of MRI to perform diffusion tensor tomography. The goal is to develop a theory for the reconstruction of both Radon planar and X-ray or line-integral projections because of the flexibility of MRI to obtain both of these type of projections in 3D. The development presented here for the linear tensor tomography problem provides insight into the structure of the nonlinear MRI diffusion tensor inverse problem. A particular application of tensor imaging in MRI is the potential application of cardiac diffusion tensor tomography for determining in vivo cardiac fiber structure. One difficulty in the cardiac application is the motion of the heart. This presents a need for developing future theory for tensor tomography in a motion field. This means developing a better understanding of the MRI signal for diffusion processes in a deforming media. The techniques developed may allow the application of MRI tensor tomography for the study of structure of fiber tracts in the brain, atherosclerotic plaque, and spine in addition to fiber structure in the heart. However, the relations presented are also applicable to other fields in medical imaging such as diffraction tomography using ultrasound. The mathematics presented can also be extended to exponential Radon transform of tensor fields and to other geometric acquisitions such as cone beam tomography of tensor fields.
Vector Acoustics, Vector Sensors, and 3D Underwater Imaging
NASA Astrophysics Data System (ADS)
Lindwall, D.
2007-12-01
Vector acoustic data has two more dimensions of information than pressure data and may allow for 3D underwater imaging with much less data than with hydrophone data. The vector acoustic sensors measures the particle motions due to passing sound waves and, in conjunction with a collocated hydrophone, the direction of travel of the sound waves. When using a controlled source with known source and sensor locations, the reflection points of the sound field can be determined with a simple trigonometric calculation. I demonstrate this concept with an experiment that used an accelerometer based vector acoustic sensor in a water tank with a short-pulse source and passive scattering targets. The sensor consists of a three-axis accelerometer and a matched hydrophone. The sound source was a standard transducer driven by a short 7 kHz pulse. The sensor was suspended in a fixed location and the hydrophone was moved about the tank by a robotic arm to insonify the tank from many locations. Several floats were placed in the tank as acoustic targets at diagonal ranges of approximately one meter. The accelerometer data show the direct source wave as well as the target scattered waves and reflections from the nearby water surface, tank bottom and sides. Without resorting to the usual methods of seismic imaging, which in this case is only two dimensional and relied entirely on the use of a synthetic source aperture, the two targets, the tank walls, the tank bottom, and the water surface were imaged. A directional ambiguity inherent to vector sensors is removed by using collocated hydrophone data. Although this experiment was in a very simple environment, it suggests that 3-D seismic surveys may be achieved with vector sensors using the same logistics as a 2-D survey that uses conventional hydrophones. This work was supported by the Office of Naval Research, program element 61153N.
Vector spherical harmonics application to 3-D tomography problem
NASA Astrophysics Data System (ADS)
Balandin, A. L.
2007-04-01
A method of series expansion with the aid of vector spherical harmonics intended for inverting line integrated data is proposed to investigate 3-D vector fields in the spherical plasmas. A set of numerical computations demonstrating the 3-D reconstruction of the model vector fields has been performed to assess the inversion method proposed.
Texture splats for 3D vector and scalar field visualization
Crawfis, R.A.; Max, N.
1993-04-06
Volume Visualization is becoming an important tool for understanding large 3D datasets. A popular technique for volume rendering is known as splatting. With new hardware architectures offering substantial improvements in the performance of rendering texture mapped objects, we present textured splats. An ideal reconstruction function for 3D signals is developed which can be used as a texture map for a splat. Extensions to the basic splatting technique are then developed to additionally represent vector fields.
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.
Coevrage Estimation of Geosensor in 3d Vector Environments
NASA Astrophysics Data System (ADS)
Afghantoloee, A.; Doodman, S.; Karimipour, F.; Mostafavi, M. A.
2014-10-01
Sensor deployment optimization to achieve the maximum spatial coverage is one of the main issues in Wireless geoSensor Networks (WSN). The model of the environment is an imperative parameter that influences the accuracy of geosensor coverage. In most of recent studies, the environment has been modeled by Digital Surface Model (DSM). However, the advances in technology to collect 3D vector data at different levels, especially in urban models can enhance the quality of geosensor deployment in order to achieve more accurate coverage estimations. This paper proposes an approach to calculate the geosensor coverage in 3D vector environments. The approach is applied on some case studies and compared with DSM based methods.
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.
3D vector flow using a row-column addressed CMUT array
NASA Astrophysics Data System (ADS)
Holbek, Simon; Christiansen, Thomas Lehrmann; Engholm, Mathias; Lei, Anders; Stuart, Mathias Bo; Beers, Christopher; Moesner, Lars Nordahl; Bagge, Jan Peter; Thomsen, Erik Vilain; Jensen, Jørgen Arendt
2016-04-01
This paper presents an in-house developed 2-D capacitive micromachined ultrasonic transducer (CMUT) applied for 3-D blood flow estimation. The probe breaks with conventional transducers in two ways; first, the ultrasonic pressure field is generated from thousands of small vibrating micromachined cells, and second, elements are accessed by row and/or column indices. The 62+62 2-D row-column addressed prototype CMUT probe was used for vector flow estimation by transmitting focused ultrasound into a flow-rig with a fully developed parabolic flow. The beam-to-flow angle was 90°. The received data was beamformed and processed offline. A transverse oscillation (TO) velocity estimator was used to estimate the 3-D vector flow along a line originating from the center of the transducer. The estimated velocities in the lateral and axial direction were close to zero as expected. In the transverse direction a characteristic parabolic velocity profile was estimated with a peak velocity of 0.48 m/s +/- 0.02 m/s in reference to the expected 0.54 m/s. The results presented are the first 3-D vector flow estimates obtained with a row-column CMUT probe, which demonstrates that the CMUT technology is feasible for 3-D flow estimation.
3D ultrasound image segmentation using wavelet support vector machines
Akbari, Hamed; Fei, Baowei
2012-01-01
Purpose: Transrectal ultrasound (TRUS) imaging is clinically used in prostate biopsy and therapy. Segmentation of the prostate on TRUS images has many applications. In this study, a three-dimensional (3D) segmentation method for TRUS images of the prostate is presented for 3D ultrasound-guided biopsy. Methods: This segmentation method utilizes a statistical shape, texture information, and intensity profiles. A set of wavelet support vector machines (W-SVMs) is applied to the images at various subregions of the prostate. The W-SVMs are trained to adaptively capture the features of the ultrasound images in order to differentiate the prostate and nonprostate tissue. This method consists of a set of wavelet transforms for extraction of prostate texture features and a kernel-based support vector machine to classify the textures. The voxels around the surface of the prostate are labeled in sagittal, coronal, and transverse planes. The weight functions are defined for each labeled voxel on each plane and on the model at each region. In the 3D segmentation procedure, the intensity profiles around the boundary between the tentatively labeled prostate and nonprostate tissue are compared to the prostate model. Consequently, the surfaces are modified based on the model intensity profiles. The segmented prostate is updated and compared to the shape model. These two steps are repeated until they converge. Manual segmentation of the prostate serves as the gold standard and a variety of methods are used to evaluate the performance of the segmentation method. Results: The results from 40 TRUS image volumes of 20 patients show that the Dice overlap ratio is 90.3% ± 2.3% and that the sensitivity is 87.7% ± 4.9%. Conclusions: The proposed method provides a useful tool in our 3D ultrasound image-guided prostate biopsy and can also be applied to other applications in the prostate. PMID:22755682
Systematic error analysis for 3D nanoprofiler tracing normal vector
NASA Astrophysics Data System (ADS)
Kudo, Ryota; Tokuta, Yusuke; Nakano, Motohiro; Yamamura, Kazuya; Endo, Katsuyoshi
2015-10-01
In recent years, demand for an optical element having a high degree of freedom shape is increased. High-precision aspherical shape is required for the X-ray focusing mirror etc. For the head-mounted display etc., optical element of the free-form surface is used. For such an optical device fabrication, measurement technology is essential. We have developed a high- precision 3D nanoprofiler. By nanoprofiler, the normal vector information of the sample surface is obtained on the basis of the linearity of light. Normal vector information is differential value of the shape, it is possible to determine the shape by integrating. Repeatability of sub-nanometer has been achieved by nanoprofiler. To pursue the accuracy of shapes, systematic error is analyzed. The systematic errors are figure error of sample and assembly errors of the device. This method utilizes the information of the ideal shape of the sample, and the measurement point coordinates and normal vectors are calculated. However, measured figure is not the ideal shape by the effect of systematic errors. Therefore, the measurement point coordinate and the normal vector is calculated again by feeding back the measured figure. Correction of errors have been attempted by figure re-derivation. It was confirmed theoretically effectiveness by simulation. This approach also applies to the experiment, it was confirmed the possibility of about 4 nm PV figure correction in the employed sample.
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
New techniques in 3D scalar and vector field visualization
Max, N.; Crawfis, R.; Becker, B.
1993-05-05
At Lawrence Livermore National Laboratory (LLNL) we have recently developed several techniques for volume visualization of scalar and vector fields, all of which use back-to-front compositing. The first renders volume density clouds by compositing polyhedral volume cells or their faces. The second is a ``splatting`` scheme which composites textures used to reconstruct the scalar or vector fields. One version calculates the necessary texture values in software, and another takes advantage of hardware texture mapping. The next technique renders contour surface polygons using semi-transparent textures, which adjust appropriately when the surfaces deform in a flow, or change topology. The final one renders the ``flow volume`` of smoke or dye tracer swept out by a fluid flowing through a small generating polygon. All of these techniques are applied to a climate model data set, to visualize cloud density and wind velocity.
View-dependent streamlines for 3D vector fields.
Marchesin, Stéphane; Chen, Cheng-Kai; Ho, Chris; Ma, Kwan-Liu
2010-01-01
This paper introduces a new streamline placement and selection algorithm for 3D vector fields. Instead of considering the problem as a simple feature search in data space, we base our work on the observation that most streamline fields generate a lot of self-occlusion which prevents proper visualization. In order to avoid this issue, we approach the problem in a view-dependent fashion and dynamically determine a set of streamlines which contributes to data understanding without cluttering the view. Since our technique couples flow characteristic criteria and view-dependent streamline selection we are able achieve the best of both worlds: relevant flow description and intelligible, uncluttered pictures. We detail an efficient GPU implementation of our algorithm, show comprehensive visual results on multiple datasets and compare our method with existing flow depiction techniques. Our results show that our technique greatly improves the readability of streamline visualizations on different datasets without requiring user intervention. PMID:20975200
A Reconstruction Approach for Imaging in 3D Cone Beam Vector Field Tomography
Schuster, T.; Theis, D.; Louis, A. K.
2008-01-01
3D cone beam vector field tomography (VFT) aims for reconstructing and visualizing the velocity field of a moving fluid by measuring line integrals of projections of the vector field. The data are obtained by ultrasound measurements along a scanning curve which surrounds the object. From a mathematical point of view, we have to deal with the inversion of the vectorial cone beam transform. Since the vectorial cone beam transform of any gradient vector field with compact support is identically equal to zero, we can only hope to reconstruct the solenoidal part of an arbitrary vector field. In this paper we will at first summarize important properties of the cone beam transform for three-dimensional solenoidal vector fields and then propose a solution approach based on the method of approximate inverse. In this context, we intensively make use of results from scalar 3D computerized tomography. The findings presented in the paper will continuously be illustrated by pictures from first numerical experiments done with exact, simulated data. PMID:19197391
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
Application of Vectors to Relative Velocity
ERIC Educational Resources Information Center
Tin-Lam, Toh
2004-01-01
The topic 'relative velocity' has recently been introduced into the Cambridge Ordinary Level Additional Mathematics syllabus under the application of Vectors. In this note, the results of relative velocity and the 'reduction to rest' technique of teaching relative velocity are derived mathematically from vector algebra, in the hope of providing…
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
Generation of 3D characterization databases in vector format
NASA Astrophysics Data System (ADS)
Wilkosz, Aaron; Williams, Bryan L.; Motz, Steve
2001-09-01
We discuss the methodology and techniques employed in transforming our 3D characterization databases and 3D target models from our internal 3D format to a more universal 3D format. Currently our 3D characterization databases and target models are encoded in an internal custom file format that targets specific simulators set up to receive out data. In order to make our databases available to a wider audience within the modeling and simulation community, we have developed techniques to transform our databases into the more common Open Flight file format. We outline the steps taken to accomplish this. We discuss the methodology and show examples of backgrounds, object discretes, and target models. The developed characterization databases are used in digital simulations by various customers within the US Army Aviation and Missile Command (AMCOM). These databases are used in closed loop dynamic simulations to evaluate the performance of various missile systems.
Vector algorithms for geometrically nonlinear 3D finite element analysis
NASA Technical Reports Server (NTRS)
Whitcomb, John D.
1989-01-01
Algorithms for geometrically nonlinear finite element analysis are presented which exploit the vector processing capability of the VPS-32, which is closely related to the CYBER 205. By manipulating vectors (which are long lists of numbers) rather than individual numbers, very high processing speeds are obtained. Long vector lengths are obtained without extensive replication or reordering by storage of intermediate results in strategic patterns at all stages of the computations. Comparisons of execution times with those from programs using either scalar or other vector programming techniques indicate that the algorithms presented are quite efficient.
Vectors in Use in a 3D Juggling Game Simulation
ERIC Educational Resources Information Center
Kynigos, Chronis; Latsi, Maria
2006-01-01
The new representations enabled by the educational computer game the "Juggler" can place vectors in a central role both for controlling and measuring the behaviours of objects in a virtual environment simulating motion in three-dimensional spaces. The mathematical meanings constructed by 13 year-old students in relation to vectors as objects, as a…
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.
Use of Colour and Interactive Animation in Learning 3D Vectors
ERIC Educational Resources Information Center
Iskander, Wejdan; Curtis, Sharon
2005-01-01
This study investigated the effects of two computer-implemented techniques (colour and interactive animation) on learning 3D vectors. The participants were 43 female Saudi Arabian high school students. They were pre-tested on 3D vectors using a paper questionnaire that consisted of calculation and visualization types of questions. The students…
A new video codec based on 3D-DTCWT and vector SPIHT
NASA Astrophysics Data System (ADS)
Xu, Ruiping; Li, Huifang; Xie, Sunyun
2011-10-01
In this paper, a new video coding system combining 3-D complex dual-tree discrete wavelet transform with vector SPIHT and arithmetic coding is proposed, and tested on standard video sequences. First the 3-D DTCWT of each color component is performed for video sequences. Then the wavelet coefficients are grouped to form vector, and successive refinement vector quantization techniques is used to quantize the groups. Finally experimental results are given. It shows that the proposed video codec provides better performance than the 3D-DTCWT and 3D-SPIHT codec, and the superior performance for the proposed sheme lies in not performing motion compensation.
Stars of the Big Dipper: A 3-D Vector Activity
NASA Astrophysics Data System (ADS)
Kuo, Vince H.; Beichner, Robert J.
2006-03-01
Most teachers of introductory physics will agree that many students have difficulty with vectors, so much so that we frequently spend a week at the beginning of the semester presenting material that students should know from previous mathematics courses. This review is often quite abstract, with little or no connection to familiar contexts, and seldom includes any motivation for students to "see it again." In this paper we present a vector activity that attempts to address both these issues using the stars of the Big Dipper, in the constellation Ursa Major, as a memorable context.
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.
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
Estimating 2-D vector velocities using multidimensional spectrum analysis.
Oddershede, Niels; Løvstakken, Lasse; Torp, Hans; Jensen, Jørgen Arendt
2008-08-01
Wilson (1991) presented an ultrasonic wideband estimator for axial blood flow velocity estimation through the use of the 2-D Fourier transform. It was shown how a single velocity component was concentrated along a line in the 2-D Fourier space, where the slope was given by the axial velocity. Later, it was shown that this approach could also be used for finding the lateral velocity component by also including a lateral sampling. A single velocity component would then be concentrated along a plane in the 3-D Fourier space, tilted according to the 2 velocity components. This paper presents 2 new velocity estimators for finding both the axial and lateral velocity components. The estimators essentially search for the plane in the 3- D Fourier space, where the integrated power spectrum is largest. The first uses the 3-D Fourier transform to find the power spectrum, while the second uses a minimum variance approach. Based on this plane, the axial and lateral velocity components are estimated. Several phantom measurements, for flow-to-depth angles of 60, 75, and 90 degrees, were performed. Multiple parallel lines were beamformed simultaneously, and 2 different receive apodization schemes were tried. The 2 estimators were then applied to the data. The axial velocity component was estimated with an average standard deviation below 2.8% of the peak velocity, while the average standard deviation of the lateral velocity estimates was between 2.0% and 16.4%. The 2 estimators were also tested on in vivo data from a transverse scan of the common carotid artery, showing the potential of the vector velocity estimation method under in vivo conditions. PMID:18986918
3D Linear Transformations in the Form of Matrix and Vector
NASA Astrophysics Data System (ADS)
Zhang, Hua
2008-11-01
In this article, results of four 3D linear transformations (translation, rotation, scale and shear) in the form of matrix and vector are simplified into a same 3D physical coordinates system. Comparing the simplified results of those four linear transformations, the results obtained from matrix form are exactly the same as what obtained from vector algebra in final expressions, although they are different from original expressions.
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.
Neighboring block based disparity vector derivation for multiview compatible 3D-AVC
NASA Astrophysics Data System (ADS)
Kang, Jewon; Chen, Ying; Zhang, Li; Zhao, Xin; Karczewicz, Marta
2013-09-01
3D-AVC being developed under Joint Collaborative Team on 3D Video Coding (JCT-3V) significantly outperforms the Multiview Video Coding plus Depth (MVC+D) which simultaneously encodes texture views and depth views with the multiview extension of H.264/AVC (MVC). However, when the 3D-AVC is configured to support multiview compatibility in which texture views are decoded without depth information, the coding performance becomes significantly degraded. The reason is that advanced coding tools incorporated into the 3D-AVC do not perform well due to the lack of a disparity vector converted from the depth information. In this paper, we propose a disparity vector derivation method utilizing only the information of texture views. Motion information of neighboring blocks is used to determine a disparity vector for a macroblock, so that the derived disparity vector is efficiently used for the coding tools in 3D-AVC. The proposed method significantly improves a coding gain of the 3D-AVC in the multiview compatible mode about 20% BD-rate saving in the coded views and 26% BD-rate saving in the synthesized views on average.
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.
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.
Numerical solution of 3-D magnetotelluric using vector finite element method
NASA Astrophysics Data System (ADS)
Prihantoro, Rudy; Sutarno, Doddy; Nurhasan
2015-09-01
Magnetotelluric (MT) is a passive electromagnetic (EM) method which measure natural variations of electric and magnetic vector fields at the Earth surface to map subsurface electrical conductivity/resistivity structure. In this study, we obtained numerical solution of three-dimensional (3-D) MT using vector finite element method by solving second order Maxwell differential equation describing diffusion of plane wave through the conductive earth. Rather than the nodes of the element, the edges of the element is used as a vector basis to overcome the occurrence of nonphysical solutions that usually faced by scalar (node based) finite element method. Electric vector fields formulation was used and the resulting system of equation was solved using direct solution method to obtain the electric vector field distribution throughout the earth resistivity model structure. The resulting MT response functions was verified with 1-D layered Earth and 3-D2 COMMEMI outcropping structure. Good agreement is achieved for both structure models.
Development and Calibration of New 3-D Vector VSP Imaging Technology: Vinton Salt Dome, LA
Kurt J. Marfurt; Hua-Wei Zhou; E. Charlotte Sullivan
2004-09-01
Vinton salt dome is located in Southwestern Louisiana, in Calcasieu Parish. Tectonically, the piercement dome is within the salt dome minibasin province. The field has been in production since 1901, with most of the production coming from Miocene and Oligocene sands. The goal of our project was to develop and calibrate new processing and interpretation technology to fully exploit the information available from a simultaneous 3-D surface seismic survey and 3-C, 3-D vertical seismic profile (VSP) survey over the dome. More specifically the goal was to better image salt dome flanks and small, reservoir-compartmentalizing faults. This new technology has application to mature salt-related fields across the Gulf Coast. The primary focus of our effort was to develop, apply, and assess the limitations of new 3-C, 3-D wavefield separation and imaging technology that could be used to image aliased, limited-aperture, vector VSP data. Through 2-D and 3-D full elastic modeling, we verified that salt flank reflections exist in the horizontally-traveling portion of the wavefield rather than up- and down-going portions of the wavefield, thereby explaining why many commercial VSP processing flow failed. Since the P-wave reflections from the salt flank are measured primarily on the horizontal components while P-wave reflections from deeper sedimentary horizons are measured primarily on the vertical component, a true vector VSP analysis was needed. We developed an antialiased discrete Radon transform filter to accurately model P- and S-wave data components measured by the vector VSP. On-the-fly polarization filtering embedded in our Kirchhoff imaging algorithm was effective in separating PP from PS wave images. By the novel application of semblance-weighted filters, we were able to suppress many of the migration artifacts associated with low fold, sparse VSP acquisition geometries. To provide a better velocity/depth model, we applied 3-D prestack depth migration to the surface data
Slawinski, J; Bonnefoy, A; Ontanon, G; Leveque, J M; Miller, C; Riquet, A; Chèze, L; Dumas, R
2010-05-28
The aim of the present study was to measure during a sprint start the joint angular velocity and the kinetic energy of the different segments in elite sprinters. This was performed using a 3D kinematic analysis of the whole body. Eight elite sprinters (10.30+/-0.14s 100 m time), equipped with 63 passive reflective markers, realised four maximal 10 m sprints start on an indoor track. An opto-electronic Motion Analysis system consisting of 12 digital cameras (250 Hz) was used to collect the 3D marker trajectories. During the pushing phase on the blocks, the 3D angular velocity vector and its norm were calculated for each joint. The kinetic energy of 16 segments of the lower and upper limbs and of the total body was calculated. The 3D kinematic analysis of the whole body demonstrated that joints such as shoulders, thoracic or hips did not reach their maximal angular velocity with a movement of flexion-extension, but with a combination of flexion-extension, abduction-adduction and internal-external rotation. The maximal kinetic energy of the total body was reached before clearing block (respectively, 537+/-59.3 J vs. 514.9+/-66.0 J; p< or =0.01). These results suggested that a better synchronization between the upper and lower limbs could increase the efficiency of pushing phase on the blocks. Besides, to understand low interindividual variances in the sprint start performance in elite athletes, a 3D complete body kinematic analysis shall be used. PMID:20226465
The design and implementation of stereoscopic 3D scalable vector graphics based on WebKit
NASA Astrophysics Data System (ADS)
Liu, Zhongxin; Wang, Wenmin; Wang, Ronggang
2014-03-01
Scalable Vector Graphics (SVG), which is a language designed based on eXtensible Markup Language (XML), is used to describe basic shapes embedded in webpages, such as circles and rectangles. However, it can only depict 2D shapes. As a consequence, web pages using classical SVG can only display 2D shapes on a screen. With the increasing development of stereoscopic 3D (S3D) technology, binocular 3D devices have been widely used. Under this circumstance, we intend to extend the widely used web rendering engine WebKit to support the description and display of S3D webpages. Therefore, the extension of SVG is of necessity. In this paper, we will describe how to design and implement SVG shapes with stereoscopic 3D mode. Two attributes representing the depth and thickness are added to support S3D shapes. The elimination of hidden lines and hidden surfaces, which is an important process in this project, is described as well. The modification of WebKit is also discussed, which is made to support the generation of both left view and right view at the same time. As is shown in the result, in contrast to the 2D shapes generated by the Google Chrome web browser, the shapes got from our modified browser are in S3D mode. With the feeling of depth and thickness, the shapes seem to be real 3D objects away from the screen, rather than simple curves and lines as before.
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.
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
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
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
Novel high speed method using gray level vector modulation for 3D shape measurement
NASA Astrophysics Data System (ADS)
Lin, Gui-Wen; Li, Dong; Tian, Jin-Dong
2014-11-01
Binocular Vision Technique is widely used in three-dimensional (3-D) measurement. Matching of pictures captured from two cameras is the most critical and difficult step in 3-D shape reconstruction. The method combines codedstructured light and spatial phase is usually adopted. However, being time consuming in matching, this method could not meet the requirements of real-time 3-D vision. In order to satisfy the high speed characteristic of real-time measurement, a novel method using gray level vector modulation is introduced. Combining binary code with gray coding principle, new coding patterns using gray level vector method is designed and projected onto the object surface. Each pixel corresponds to the designed sequence of gray values as a feature vector. The unique gray level vector is then dimensionally reduced to a resulting value which could be used as characteristic information for binocular matching. Experimental results further demonstrated the correctness and feasibility of the proposed method with fewer component patterns and less computational time.
Automatic classification of 3D segmented CT data using data fusion and support vector machine
NASA Astrophysics Data System (ADS)
Osman, Ahmad; Kaftandjian, Valérie; Hassler, Ulf
2011-07-01
The three dimensional X-ray computed tomography (3D-CT) has proved its successful usage as inspection method in non destructive testing. The generated 3D volume using high efficiency reconstruction algorithms contains all the inner structures of the inspected part. Segmentation of this volume reveals suspicious regions which need to be classified into defects or false alarms. This paper deals with the classification step using data fusion theory and support vector machine. Results achieved are very promising and prove the effectiveness of the data fusion theory as a method to build stronger classifier.
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
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
NASA Astrophysics Data System (ADS)
Hinrichs, Brant E.
2010-10-01
An intermediate E&M course (i.e. based on Griffiths [1]) involves the extensive integration of vector calculus concepts and notation with abstract physics concepts like field and potential. We hope that students take what they have learned in their math courses and apply it to help represent and make sense of the physics. To assess how well students are able to do this integration and application I have developed several simple concept tests on position and unit vectors in non-Cartesian coordinate systems as they are used in intermediate E&M. In this paper I describe one of these concept tests and present results that show both undergraduate physics majors and physics graduate students have difficulty using spherical unit vectors to write position vectors in 3-d space.
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.
Quantitative 3D Tracing of Gene-delivery Viral Vectors in Human Cells and Animal Tissues
Xiao, Ping-Jie; Li, Chengwen; Neumann, Aaron; Samulski, R Jude
2012-01-01
Trafficking through a variety of cellular structures and organelles is essential for the interaction between gene-delivery vectors (i.e., adeno-associated virus (AAV) and liposomes) and host cells/tissues. Here, we present a method of computer-assisted quantitative 3D biodistribution microscopy that samples the whole population of fluorescently-labeled vectors and document their trafficking routes. Using AAV as a working model, we first experimentally defined numerical parameters for the singularity of Cy5-labeled particles by combining confocal microscopy and atomic force microscopy (AFM). We then developed a robust approach that integrates single-particle fluorescence imaging with 3D deconvolution and isosurface rendering to quantitate viral distribution and trafficking in human cells as well as animal tissues at the single-particle level. Using this quantitative method, we uncovered an as yet uncharacterized rate-limiting step during viral cell entry, while delineating nuclear accumulation of virions during the first 8 hours postinfection. Further, our studies revealed for the first time that following intramuscular injection, AAV spread progressively across muscle tissues through endomysium between myofibers instead of traversing through target cells. Such 3D resolution and quantitative dissection of vector–host interactions at the subcellular level should significantly improve our ability to resolve trafficking mechanisms of gene-delivery particles and facilitate the development of enhanced viral vectors. PMID:22108857
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.
Sarrami-Foroushani, Ali; Nasr Esfahany, Mohsen; Nasiraei Moghaddam, Abbas; Saligheh Rad, Hamidreza; Firouznia, Kavous; Shakiba, Madjid; Ghanaati, Hossein; Wilkinson, Iain David; Frangi, Alejandro Federico
2015-01-01
Background: Understanding hemodynamic environment in vessels is important for realizing the mechanisms leading to vascular pathologies. Objectives: Three-dimensional velocity vector field in carotid bifurcation is visualized using TR 3D phase-contrast magnetic resonance imaging (TR 3D PC MRI) and computational fluid dynamics (CFD). This study aimed to present a qualitative and quantitative comparison of the velocity vector field obtained by each technique. Subjects and Methods: MR imaging was performed on a 30-year old male normal subject. TR 3D PC MRI was performed on a 3 T scanner to measure velocity in carotid bifurcation. 3D anatomical model for CFD was created using images obtained from time-of-flight MR angiography. Velocity vector field in carotid bifurcation was predicted using CFD and PC MRI techniques. A statistical analysis was performed to assess the agreement between the two methods. Results: Although the main flow patterns were the same for the both techniques, CFD showed a greater resolution in mapping the secondary and circulating flows. Overall root mean square (RMS) errors for all the corresponding data points in PC MRI and CFD were 14.27% in peak systole and 12.91% in end diastole relative to maximum velocity measured at each cardiac phase. Bland-Altman plots showed a very good agreement between the two techniques. However, this study was not aimed to validate any of methods, instead, the consistency was assessed to accentuate the similarities and differences between Time-resolved PC MRI and CFD. Conclusion: Both techniques provided quantitatively consistent results of in vivo velocity vector fields in right internal carotid artery (RCA). PC MRI represented a good estimation of main flow patterns inside the vasculature, which seems to be acceptable for clinical use. However, limitations of each technique should be considered while interpreting results. PMID:26793288
A 3D Vector/Scalar Visualization and Particle Tracking Package
Energy Science and Technology Software Center (ESTSC)
1999-08-19
BOILERMAKER is an interactive visualization system consisting of three components: a visualization component, a particle tracking component, and a communication layer. The software, to date, has been used primarily in the visualization of vector and scalar fields associated with computational fluid dynamics (CFD) models of flue gas flows in industrial boilers and incinerators. Users can interactively request and toggle static vector fields, dynamic streamlines, and flowing vector fields. In addition, the user can interactively placemore » injector nozzles on boiler walls and visualize massed, evaporating sprays emanating from them. Some characteristics of the spray can be adjusted from within the visualization environment including spray shape and particle size. Also included with this release is software that supports 3D menu capabilities, scrollbars, communication and navigation.« less
A 3D Vector/Scalar Visualization and Particle Tracking Package
Freitag, Lori; Disz, Terry; Papka, Mike; Heath, Daniel; Diachin, Darin; Herzog, Jim; Ryan, and Bob
1999-08-19
BOILERMAKER is an interactive visualization system consisting of three components: a visualization component, a particle tracking component, and a communication layer. The software, to date, has been used primarily in the visualization of vector and scalar fields associated with computational fluid dynamics (CFD) models of flue gas flows in industrial boilers and incinerators. Users can interactively request and toggle static vector fields, dynamic streamlines, and flowing vector fields. In addition, the user can interactively place injector nozzles on boiler walls and visualize massed, evaporating sprays emanating from them. Some characteristics of the spray can be adjusted from within the visualization environment including spray shape and particle size. Also included with this release is software that supports 3D menu capabilities, scrollbars, communication and navigation.
Validation of INSAT-3D atmospheric motion vectors for monsoon 2015
NASA Astrophysics Data System (ADS)
Sharma, Priti; Rani, S. Indira; Das Gupta, M.
2016-05-01
Atmospheric Motion Vector (AMV) over Indian Ocean and surrounding region is one of the most important sources of tropospheric wind information assimilated in numerical weather prediction (NWP) system. Earlier studies showed that the quality of Indian geo-stationary satellite Kalpana-1 AMVs was not comparable to that of other geostationary satellites over this region and hence not used in NWP system. Indian satellite INSAT-3D was successfully launched on July 26, 2013 with upgraded imaging system as compared to that of previous Indian satellite Kalpana-1. INSAT-3D has middle infrared band (3.80 - 4.00 μm) which is capable of night time pictures of low clouds and fog. Three consecutive images of 30-minutes interval are used to derive the AMVs. New height assignment scheme (using NWP first guess and replacing old empirical GA method) along with modified quality control scheme were implemented for deriving INSAT-3D AMVs. In this paper an attempt has been made to validate these AMVs against in-situ observations as well as against NCMRWF's NWP first guess for monsoon 2015. AMVs are subdivided into three different pressure levels in the vertical viz. low (1000 - 700 hPa), middle (700 - 400 hPa) and high (400 - 100 hPa) for validation purpose. Several statistics viz. normalized root mean square vector difference; biases etc. have been computed over different latitudinal belt. Result shows that the general mean monsoon circulations along with all the transient monsoon systems are well captured by INSAT-3D AMVs, as well as the error statistics viz., RMSE etc of INSAT-3D AMVs is now comparable to other geostationary satellites.
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.
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.
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 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
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
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
NASA Astrophysics Data System (ADS)
Chen, Yu; Xie, Xilin
2016-05-01
E and Liu [J. Comput. Phys. 138 (1997) 57-82] put forward a finite difference method for 3D viscous incompressible flows in the vorticity-vector potential formulation on non-staggered grids. In this paper, we will extend this method to the case of flows in the presence of a deformable surface. By use of two kinds of surface differential operators, the implementation of boundary conditions on a plane is generalized to a curved smooth surface with given velocity distribution, whether this be an inflow/outflow interface or a curved wall. To deal with the irregular and varying physical domain, time-dependent curvilinear coordinates are constructed and the corresponding tensor analysis is adopted in deriving the component form of the governing equations. Therefore, the equations can be discretized and solved in a regular and fixed parametric domain. Numerical results are presented for a 3D lid-driven cavity with a deforming surface and a 3D duct flow with a deforming boundary. A new way to validate numerical simulations is proposed based on an expression for the rate-of-strain tensor on a deformable surface.
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.
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.
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.
NASA Technical Reports Server (NTRS)
Mangalgiri, P. D.; Prabhakaran, R.
1986-01-01
An algorithm for vectorized computation of stiffness matrices of an 8 noded isoparametric hexahedron element for geometric nonlinear analysis was developed. This was used in conjunction with the earlier 2-D program GAMNAS to develop the new program NAS3D for geometric nonlinear analysis. A conventional, modified Newton-Raphson process is used for the nonlinear analysis. New schemes for the computation of stiffness and strain energy release rates is presented. The organization the program is explained and some results on four sample problems are given. The study of CPU times showed that savings by a factor of 11 to 13 were achieved when vectorized computation was used for the stiffness instead of the conventional scalar one. Finally, the scheme of inputting data is explained.
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
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.
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
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.
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.
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
NASA Astrophysics Data System (ADS)
Barnhart, Donald H.; Chan, Victor S. S.; Halliwell, Neil A.; Coupland, Jeremy M.
1999-10-01
This paper introduces a new approach to 3D displacement and velocity measurements that unifies the disciplines of holographic interferometry and holographic particle image velocimetry (HPIV). Equally applicable to fluid and solid mechanics, the overall system enables quantitative displacement measurements between two holographically recorded events from either particle or surface scattering sites, working with both pulsed and continuous-wave laser systems. The resulting measurements exhibit an accuracy corresponding to interferometric system, but with a dynamic range found with PIV systems. Most importantly, this paper introduces the novel use of an optical fiber to specify the measurement points, remove optical aberrations of windows, and eliminate directional ambiguity. An optical fiber is used to probe the recorded holographic image space at each 3D measurement point in order to extract the 3D displacement vectors. This fiber system also employs a novel optical image shifting method to eliminate the problem of directional ambiguity. In addition, the reported system uses 3D complex optical correlation rather than 2D real digital correlation. It is therefore a simple matter to directly obtain 3D displacement and velocity measurements at precisely known 3D locations in the object space. By correlating both the amplitude and phase information in the holographic image, this system can measure spatial distributions of displacements even when the presence of severe aberrations preclude the detection of sharp images.
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.
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.
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.
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.
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
PAB3D Simulations of a Nozzle with Fluidic Injection for Yaw Thrust-Vector Control
NASA Technical Reports Server (NTRS)
Deere, Karen A.
1998-01-01
An experimental and computational study was conducted on an exhaust nozzle with fluidic injection for yaw thrust-vector control. The nozzle concept was tested experimentally in the NASA Langley Jet Exit Test Facility (JETF) at nozzle pressure ratios up to 4 and secondary fluidic injection flow rates up to 15 percent of the primary flow rate. Although many injection-port geometries and two nozzle planforms (symmetric and asymmetric) were tested experimentally, this paper focuses on the computational results of the more successful asymmetric planform with a slot injection port. This nozzle concept was simulated with the Navier-Stokes flow solver, PAB3D, invoking the Shih, Zhu, and Lumley algebraic Reynolds stress turbulence model (ASM) at nozzle pressure ratios (NPRs) of 2,3, and 4 with secondary to primary injection flow rates (w(sub s)/w(sub p)) of 0, 2, 7 and 10 percent.
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.
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).
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.
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
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
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.
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
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.
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
Measuring the 3-D wind vector with a weight-shift microlight aircraft
NASA Astrophysics Data System (ADS)
Metzger, S.; Junkermann, W.; Butterbach-Bahl, K.; Schmid, H. P.; Foken, T.
2011-02-01
This study investigates whether the 3-D wind vector can be measured reliably from a highly transportable and low-cost weight-shift microlight aircraft. Therefore we draw up a transferable procedure to accommodate flow distortion originating from the aircraft body and -wing. This procedure consists of the analysis of aircraft dynamics and seven successive calibration steps. For our aircraft the horizontal wind components receive their greatest single amendment (14%, relative to the initial uncertainty) from the correction of flow distortion magnitude in the dynamic pressure computation. Conversely the vertical wind component is most of all improved (31%) by subsequent steps considering the 3-D flow distortion distribution in the flow angle computations. Therein the influences of the aircraft's aeroelastic wing (53%), as well as sudden changes in wing loading (16%) are considered by using the measured lift coefficient as explanatory variable. Three independent lines of analysis are used to evaluate the quality of the wind measurement: (a) A wind tunnel study in combination with the propagation of sensor uncertainties defines the systems input uncertainty to ≈0.6 m s-1 at the extremes of a 95% confidence interval. (b) During severe vertical flight manoeuvres the deviation range of the vertical wind component does not exceed 0.3 m s-1. (c) The comparison with ground based wind measurements yields an overall operational uncertainty (root mean square deviation) of ≈0.4 m s-1 for the horizontal and ≍0.3 m s-1 for the vertical wind components. No conclusive dependence of the uncertainty on the wind magnitude (<8 m s-1) or true airspeed (ranging from 23-30 m s-1) is found. Hence our analysis provides the necessary basis to study the wind measurement precision and spectral quality, which is prerequisite for reliable eddy-covariance flux measurements.
Measuring the 3-D wind vector with a weight-shift microlight aircraft
NASA Astrophysics Data System (ADS)
Metzger, S.; Junkermann, W.; Butterbach-Bahl, K.; Schmid, H. P.; Foken, T.
2011-07-01
This study investigates whether the 3-D wind vector can be measured reliably from a highly transportable and low-cost weight-shift microlight aircraft. Therefore we draw up a transferable procedure to accommodate flow distortion originating from the aircraft body and -wing. This procedure consists of the analysis of aircraft dynamics and seven successive calibration steps. For our aircraft the horizontal wind components receive their greatest single amendment (14 %, relative to the initial uncertainty) from the correction of flow distortion magnitude in the dynamic pressure computation. Conversely the vertical wind component is most of all improved (31 %) by subsequent steps considering the 3-D flow distortion distribution in the flow angle computations. Therein the influences of the aircraft's trim (53 %), as well as changes in the aircraft lift (16 %) are considered by using the measured lift coefficient as explanatory variable. Three independent lines of analysis are used to evaluate the quality of the wind measurement: (a) A wind tunnel study in combination with the propagation of sensor uncertainties defines the systems input uncertainty to ≈0.6 m s-1 at the extremes of a 95 % confidence interval. (b) During severe vertical flight manoeuvres the deviation range of the vertical wind component does not exceed 0.3 m s-1. (c) The comparison with ground based wind measurements yields an overall operational uncertainty (root mean square error) of ≈0.4 m s-1 for the horizontal and ≈0.3 m s-1 for the vertical wind components. No conclusive dependence of the uncertainty on the wind magnitude (<8 m s-1) or true airspeed (ranging from 23-30 m s-1) is found. Hence our analysis provides the necessary basis to study the wind measurement precision and spectral quality, which is prerequisite for reliable Eddy-Covariance flux measurements.
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
Jing, Zhang; Sheng, Kang Bao
2016-01-01
To assist physicians to quickly find the required 3D model from the mass medical model, we propose a novel retrieval method, called DRFVT, which combines the characteristics of dimensionality reduction (DR) and feature vector transformation (FVT) method. The DR method reduces the dimensionality of feature vector; only the top M low frequency Discrete Fourier Transform coefficients are retained. The FVT method does the transformation of the original feature vector and generates a new feature vector to solve the problem of noise sensitivity. The experiment results demonstrate that the DRFVT method achieves more effective and efficient retrieval results than other proposed methods. PMID:27293478
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
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
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.
Liu, Yu-Ying; Ishikawa, Hiroshi; Chen, Mei; Wollstein, Gadi; Schuman, Joel S.; Rehg, James M.
2013-01-01
We develop an automated method to determine the foveola location in macular 3D-OCT images in either healthy or pathological conditions. Structural Support Vector Machine (S-SVM) is trained to directly predict the location of the foveola, such that the score at the ground truth position is higher than that at any other position by a margin scaling with the associated localization loss. This S-SVM formulation directly minimizes the empirical risk of localization error, and makes efficient use of all available training data. It deals with the localization problem in a more principled way compared to the conventional binary classifier learning that uses zero-one loss and random sampling of negative examples. A total of 170 scans were collected for the experiment. Our method localized 95.1% of testing scans within the anatomical area of the foveola. Our experimental results show that the proposed method can effectively identify the location of the foveola, facilitating diagnosis around this important landmark. PMID:23285565
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
VECTUM. Irregular 2D Velocity Vector Field Plotting Package
McClurg, F.R.; Mousseau, V.A.
1992-05-04
VECTUM is a NCAR Graphics based package, for generating a plot of an irregular 2D velocity vector field. The program reads an ASCII database of x, y, u, v, data pairs and produces a plot in Computer Graphics Metafile (CGM) format. The program also uses an ASCII parameter file for controlling annotation details such as the plot title, arrowhead style, scale of vectors, windowing, etc. Simple geometry (i.e. lines, arcs, splines) can be defined to be included with the velocity vectors. NCAR Graphics drivers can be used to display the CGM file into PostScript, HPGL, HDF, etc, output.
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.
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.
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)
Young, D. L.; Tsai, C. H.; Wu, C. S.
2015-11-01
An alternative vector potential formulation is used to solve the Navier-Stokes (N-S) equations in 3D incompressible viscous flow problems with and without through-flow boundaries. Difficulties of the vector potential formulation include the implementation of boundary conditions for through-flow boundaries and the numerical treatment of fourth-order partial differential equations. The advantages on the other hand are the automatic satisfaction of the continuity equation; and pressure is decoupled from the velocity. The objective of this paper is to introduce the appropriate gauge and boundary conditions on the vector potential formulation by a localized meshless method. To handle the divergence-free property, a Coulomb gauge condition is enforced on the vector potential to ensure its existence and uniqueness mathematically. We further improve the algorithm to through-flow problems for the boundary conditions of vector potential by introducing the concept of Stokes' theorem. Based on this innovation, there is no need to include an additional variable to tackle the through-flow fields. This process will greatly simplify the imposition of boundary conditions by the vector potential approach. Under certain conditions, the coupled fourth-order partial differential equations can be easily solved by using this meshless local differential quadrature (LDQ) method. Due to the LDQ capability to deal with the high order differential equations, this algorithm is very attractive to solve this fourth-order vector potential formulation for the N-S equations as comparing to the conventional numerical schemes such as finite element or finite difference methods. The proposed vector potential formulation is simpler and has improved accuracy and efficiency compared to other pressure-free or pressure-coupled algorithms. This investigation can be regarded as the first complete study to obtain the N-S solutions by vector potential formulation through a LDQ method. Two classic 3D benchmark
NASA Astrophysics Data System (ADS)
Thurow, Brian; Johnson, Kyle; Kim, Taehoon; Blois, Gianluca; Best, Jim; Christensen, Ken
2014-11-01
The application of Plenoptic PIV in a Refractive Index Matched (RIM) facility housed at Illinois is presented. Plenoptic PIV is an emerging 3D diagnostic that exploits the light-field imaging capabilities of a plenoptic camera. Plenoptic cameras utilize a microlens array to measure the position and angle of light rays captured by the camera. 3D/3C velocity fields are determined through application of the MART algorithm for volume reconstruction and a conventional 3D cross-correlation PIV algorithm. The RIM facility is a recirculating tunnel with a 62.5% aqueous solution of sodium iodide used as the working fluid. Its resulting index of 1.49 is equal to that of acrylic. Plenoptic PIV was used to measure the 3D velocity field of a turbulent boundary layer flow over a smooth wall, a single wall-mounted hemisphere and a full array of hemispheres (i.e. a rough wall) with a k/ δ ~ 4.6. Preliminary time averaged and instantaneous 3D velocity fields will be presented. This material is based upon work supported by the National Science Foundation under Grant No. 1235726.
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.
NASA Astrophysics Data System (ADS)
Zhang, Mian; Huang, Cheng-li
2012-08-01
Generalized spherical harmonics (GSH) are usually applied on the problems where the Earth model is elliptical and elastic stress tensor is involved in, as stress tensor can’t be represented in vector spherical harmonics. However, the divergence of the te ns or and a vector dot - product with the tensor are only needed on computation rotation modes of the Earth which can be written in the vector spherical harmonics. We extend the equations on the spherical Earth to asymmetric 3D model by means of linear operator method. This method doesn’t use the complicated generalized spherical harmonics nor Wigner 3 - j symbol. As a validation of this method, the practical calculation of rotational modes of 3D Earth will be made and discussed.
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.
Retrieving 3D Velocity Fields of Glaciers from X-band SAR Data and Comparison with GPS Observations
NASA Astrophysics Data System (ADS)
Magnússon, E.; Nagler, T.; Hetzenecker, M.; Palsson, F.; Scharrer, K.; Floricioiu, D.; Berthier, E.; Gudmundsson, S.; Rott, H.
2013-12-01
We present 3D velocity fields obtained from time series of TerraSAR-X and TanDEM-X images acquired over the ablation area of the Breidamerkurjökull outlet glacier of Vatnjökull Ice Cap (Iceland) in 2008-2012. Coherent and incoherent offset tracking is applied to repeat pass X-Band data to obtain ice displacement in cross and along track direction. Three methods are tested and compared to extract fields of the 3D ice velocity. First, the conventional surface parallel approach, which we consider as an approximation for deriving the horizontal motion rate, but does not reveal a realistic vertical motion. Second, the combination of offset tracking results from almost simultaneous observations from ascending and descending orbits measuring the glacier motion in four different directions, allowing calculation of the 3D velocity fields without any additional approximations. Third, deriving full 3D velocity fields by using the horizontal flow direction, derived from the ascending-descending combination, as constrain on offset tracking results from a single pair of SAR images. The latter two methods reveal a measurement of the vertical ice motion plus ablation, hence equivalent to the vertical motion component measured by GPS station fixed on a platform laying on the ice surface. The results from all methods are compared with such GPS measurements recorded by permanent stations on the glacier in 2008-2012 and the errors of the different methods are calculated. Additionally, we approximate the contribution of these 3D flow fields to elevation changes (emergence/submergence velocity plus net balance) and compare it with elevation changes from surface DEMs obtained in 2008 (SPIRIT), 2010 (airborne LIDAR) and 2012 (TanDEM-X).
Creating Realistic 3D Graphics with Excel at High School--Vector Algebra in Practice
ERIC Educational Resources Information Center
Benacka, Jan
2015-01-01
The article presents the results of an experiment in which Excel applications that depict rotatable and sizable orthographic projection of simple 3D figures with face overlapping were developed with thirty gymnasium (high school) students of age 17-19 as an introduction to 3D computer graphics. A questionnaire survey was conducted to find out…
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
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.
Wiimote Experiments: 3-D Inclined Plane Problem for Reinforcing the Vector Concept
ERIC Educational Resources Information Center
Kawam, Alae; Kouh, Minjoon
2011-01-01
In an introductory physics course where students first learn about vectors, they oftentimes struggle with the concept of vector addition and decomposition. For example, the classic physics problem involving a mass on an inclined plane requires the decomposition of the force of gravity into two directions that are parallel and perpendicular to the…
On the critical one-component velocity regularity criteria to 3-D incompressible MHD system
NASA Astrophysics Data System (ADS)
Liu, Yanlin
2016-05-01
Let (u , b) be a smooth enough solution of 3-D incompressible MHD system. We prove that if (u , b) blows up at a finite time T*, then for any p ∈ ] 4 , ∞ [, there holds ∫0T* (‖u3(t‧) ‖ H ˙ 1/2 +2/p p + ‖b(t‧) ‖ H ˙ 1/2 +2/p p) dt‧ = ∞. We remark that all these quantities are in the critical regularity of the MHD system.
Krams, R; Wentzel, J J; Cespedes, I; Vinke, R; Carlier, S; van der Steen, A F; Lancee, C T; Slager, C J
1999-06-01
Novel measurement techniques based on intravenous ultrasound (IVUS) technology ('IVUS-Flowmetry') require the location of a catheter inside the coronary bed. The present study quantifies disturbances in the 3-D velocity profile induced by catheter placement inside a tube, applying computational fluid dynamics. Two curved, circular meshes (radius K = 0.025 m and K = 0.035 m) with and without a catheter inside the lumen were applied. The catheter was located at the inner curve, the outer curve and at the top position. Boundary conditions were: no slip on the wall, zero stress at the outlet, uniform inflow with entrance velocities of 0.1, 0.2 and 0.4 m/s. Curvature-associated centrifugal forces shifted the maximal velocity to the outer curve and introduced two symmetrical vortices. Additional catheter placement redistributed the 3-D axial velocity field away from the catheter, which was accompanied by the appearance of multiple low-strength vortices. In addition, peak axial velocity increased, peak secondary velocities decreased, axial pressure drop increased and shear stress increased. Flow calculations simulated to resemble IVUS-based flowmetry changed by only 1% after considering secondary velocity. In conclusion, placement of a catheter inside a curved tube resembling the human coronary system changes the velocity field and reduces secondary patterns. The present study supports the usefulness of catheter-based flowmetry during resting flow conditions. During hyperemic flow conditions, flow measurements might be accompanied by large axial pressure drops because the catheter, itself, might act as a significant stenosis. PMID:10414897
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
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
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.
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.
The effect of sliding velocity on chondrocytes activity in 3D scaffolds.
Wimmer, Markus A; Alini, Mauro; Grad, Sibylle
2009-03-11
Sliding motion and shear are important mediators for the synthesis of cartilage matrix and surface molecules. This study investigated the effects of velocity magnitude and motion path on the response of bovine chondrocytes cultured in polyurethane scaffolds and subjected to oscillation against a ceramic ball. In order to vary velocity magnitude, the ball oscillated +/-25 degrees at 0.01, 0.1, and 1Hz to generate 0.28, 2.8, and 28mm/s, respectively. The median velocity of these 'open' motion trajectories was tested against 'closed' motion trajectories in that the scaffold oscillated +/-20 degrees against the ball at 1Hz, reaching 2.8mm/s. Constructs were loaded twice a day for 1h over 5 days. Gene expression of cartilage oligomeric matrix protein (COMP), proteoglycan 4 (PRG4, lubricin), and hyaluronan synthase 1 (HAS1) and release of COMP, PRG4, and hyaluronan (HA) were analyzed. Velocity magnitude determined both gene expression and release of target molecules. Using regression analysis, there was a positive and significant relationship with all outcome variables. However, only COMP reacted significantly at 0.28mm/s, while all other measured variables were considerably up-regulated at 28mm/s. Motion path characteristics affected COMP, but not PRG4 and HAS1/HA. To conclude, velocity magnitude is a critical determinant for cellular responses in tissue engineered cartilage constructs. The motion type also plays a role. However, different molecules are affected in different ways. A molecule specific velocity threshold appears necessary to induce a significant response. This should be considered in further studies investigating the effects of continuous or intermittent motion. PMID:19152917
Kinematics and aerodynamics of the velocity vector roll
NASA Technical Reports Server (NTRS)
Durham, Wayne C.; Lutze, Frederick H.; Mason, W.
1993-01-01
The velocity vector roll is an angular rotation of an aircraft about its instantaneous velocity vector, constrained to be performed at constant angle-of-attack (AOA), no sideslip, and constant velocity. Consideration of the aerodynamic force equations leads to requirements for body-axis yawing and pitching rotations that satisfy these constraints. Here, the body axis rotations, and the constraints, are used in the moment equations to determine the aerodynamic moments required to perform the velocity vector roll. For representative tactical aircraft, the conditions for maximum pitching moment are a function of orientation, occurring at about 90 deg of bank in a level trajectory. Maximum required pitching moment occurs at peak roll rate, and is achieved at AOA above 45 deg. The conditions for maximum rolling moment depend on the value of the roll mode time constant. For a small time constant (fast response) the maximum rolling moment occurs at maximum roll acceleration and zero AOA, largely independent of aircraft orientation; for a large time constant, maximum required rolling moment occurs at maximum roll rate, at maximum AOA, and at 180 deg of bank in level flight. Maximum yawing moment occurs at maximum roll acceleration, maximum AOA, and is largely independent of airplane orientation.
The Local Stellar Velocity Field via Vector Spherical Harmonics
NASA Technical Reports Server (NTRS)
Makarov, V. V.; Murphy, D. W.
2007-01-01
We analyze the local field of stellar tangential velocities for a sample of 42,339 nonbinary Hipparcos stars with accurate parallaxes, using a vector spherical harmonic formalism.We derive simple relations between the parameters of the classical linear model (Ogorodnikov-Milne) of the local systemic field and low-degree terms of the general vector harmonic decomposition. Taking advantage of these relationships, we determine the solar velocity with respect to the local stars of (V(sub X), V(sub Y), V(sub Z)) = (10.5, 18.5, 7.3) +/- 0.1 km s(exp -1) not for the asymmetric drift with respect to the local standard of rest. If only stars more distant than 100 pc are considered, the peculiar solar motion is (V(sub X), V(sub Y), V(sub Z)) = (9.9, 15.6, 6.9) +/- 0.2 km s(exp -1). The adverse effects of harmonic leakage, which occurs between the reflex solar motion represented by the three electric vector harmonics in the velocity space and higher degree harmonics in the proper-motion space, are eliminated in our analysis by direct subtraction of the reflex solar velocity in its tangential components for each star...
Rule-based fuzzy vector median filters for 3D phase contrast MRI segmentation
NASA Astrophysics Data System (ADS)
Sundareswaran, Kartik S.; Frakes, David H.; Yoganathan, Ajit P.
2008-02-01
Recent technological advances have contributed to the advent of phase contrast magnetic resonance imaging (PCMRI) as standard practice in clinical environments. In particular, decreased scan times have made using the modality more feasible. PCMRI is now a common tool for flow quantification, and for more complex vector field analyses that target the early detection of problematic flow conditions. Segmentation is one component of this type of application that can impact the accuracy of the final product dramatically. Vascular segmentation, in general, is a long-standing problem that has received significant attention. Segmentation in the context of PCMRI data, however, has been explored less and can benefit from object-based image processing techniques that incorporate fluids specific information. Here we present a fuzzy rule-based adaptive vector median filtering (FAVMF) algorithm that in combination with active contour modeling facilitates high-quality PCMRI segmentation while mitigating the effects of noise. The FAVMF technique was tested on 111 synthetically generated PC MRI slices and on 15 patients with congenital heart disease. The results were compared to other multi-dimensional filters namely the adaptive vector median filter, the adaptive vector directional filter, and the scalar low pass filter commonly used in PC MRI applications. FAVMF significantly outperformed the standard filtering methods (p < 0.0001). Two conclusions can be drawn from these results: a) Filtering should be performed after vessel segmentation of PC MRI; b) Vector based filtering methods should be used instead of scalar techniques.
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
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
Application of edge-based finite elements and vector ABCs in 3D scattering
NASA Technical Reports Server (NTRS)
Chatterjee, A.; Jin, J. M.; Volakis, John L.
1992-01-01
A finite element absorbing boundary condition (FE-ABC) solution of the scattering by arbitrary 3-D structures is considered. The computational domain is discretized using edge-based tetrahedral elements. In contrast to the node-based elements, edge elements can treat geometries with sharp edges, are divergence-less, and easily satisfy the field continuity condition across dielectric interfaces. They do, however, lead to a higher unknown count but this is balanced by the greater sparsity of the resulting finite element matrix. Thus, the computation time required to solve such a system iteratively with a given degree of accuracy is less than the traditional node-based approach. The purpose is to examine the derivation and performance of the ABC's when applied to 2-D and 3-D problems and to discuss the specifics of our FE-ABC implementation.
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.
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
NASA Astrophysics Data System (ADS)
Deb, S. K.; Kishtawal, C. M.; Kumar, Prashant; Kiran Kumar, A. S.; Pal, P. K.; Kaushik, Nitesh; Sangar, Ghansham
2016-03-01
The advanced Indian meteorological geostationary satellite INSAT-3D was launched on 26 July 2013 with an improved imager and an infrared sounder and is placed at 82°E over the Indian Ocean region. With the advancement in retrieval techniques of different atmospheric parameters and with improved imager data have enhanced the scope for better understanding of the different tropical atmospheric processes over this region. The retrieval techniques and accuracy of one such parameter, Atmospheric Motion Vectors (AMV) has improved significantly with the availability of improved spatial resolution data along with more options of spectral channels in the INSAT-3D imager. The present work is mainly focused on providing brief descriptions of INSAT-3D data and AMV derivation processes using these data. It also discussed the initial quality assessment of INSAT-3D AMVs for a period of six months starting from 01 February 2014 to 31 July 2014 with other independent observations: i) Meteosat-7 AMVs available over this region, ii) in-situ radiosonde wind measurements, iii) cloud tracked winds from Multi-angle Imaging Spectro-Radiometer (MISR) and iv) numerical model analysis. It is observed from this study that the qualities of newly derived INSAT-3D AMVs are comparable with existing two versions of Meteosat-7 AMVs over this region. To demonstrate its initial application, INSAT-3D AMVs are assimilated in the Weather Research and Forecasting (WRF) model and it is found that the assimilation of newly derived AMVs has helped in reduction of track forecast errors of the recent cyclonic storm NANAUK over the Arabian Sea. Though, the present study is limited to its application to one case study, however, it will provide some guidance to the operational agencies for implementation of this new AMV dataset for future applications in the Numerical Weather Prediction (NWP) over the south Asia region.
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
NASA Astrophysics Data System (ADS)
Qian, Long; Lu, Yong; Zhong, Wenqi; Chen, Xi; Ren, Bing; Jin, Baosheng
2013-07-01
A novel fiber high speed photography method has been developed to measure particle phase velocity in a dense gas-solid flow. The measurement system mainly includes a fiber-optic endoscope, a high speed video camera, a metal halide light source and a powerful computer with large memory. The endoscope which could be inserted into the reactors is used to form motion images of particles within the measurement window illuminated by the metal halide lamp. These images are captured by the high speed video camera and processed through a series of digital image processing algorithms, such as calibration, denoising, enhancement and binarization in order to improve the image quality. Then particles' instantaneous velocity is figured out by tracking each particle in consecutive frames. Particle phase velocity is statistically calculated according to the probability of particle velocity in each frame within a time period. This system has been applied to the investigation of particles fluidization characteristics in a 3D spouted bed. The experimental results indicate that the particle fluidization feature in the region investigated could be roughly classified into three sections by particle phase vertical velocity and the boundary between the first section and the second is the surface where particle phase velocity tends to be 0, which is in good agreement with the results published in other literature.
Asymptotic stability and bifurcations of 3D piecewise smooth vector fields
NASA Astrophysics Data System (ADS)
Carvalho, Tiago; Teixeira, Marco Antônio; Tonon, Durval José
2016-04-01
The paper deals with the analysis of the behavior of a nonsmooth three-dimensional vector field around a typal singularity. We focus on a class of generic one-parameter families {Z_{λ}} of Filippov systems and address the persistence problem for the asymptotic stability when the parameter varies near the bifurcation value {λ = 0}.
Edge-based finite elements and vector ABCs applied to 3D scattering
NASA Technical Reports Server (NTRS)
Chatterjee, A.; Jin, J. M.; Volakis, John L.
1992-01-01
An edge based finite element formulation with vector absorbing boundary conditions is presented for scattering by composite structures having boundaries satisfying impedance and/or transition conditions. Remarkably accurate results are obtained by placing the mesh a small fraction of a wavelength away from the scatterer.
Analysis of non linear partially standing waves from 3D velocity measurements
NASA Astrophysics Data System (ADS)
Drevard, D.; Rey, V.; Svendsen, Ib; Fraunie, P.
2003-04-01
Surface gravity waves in the ocean exhibit an energy spectrum distributed in both frequency and direction of propagation. Wave data collection is of great importance in coastal zones for engineering and scientific studies. In particular, partially standing waves measurements near coastal structures and steep or barred beaches may be a requirement, for instance for morphodynamic studies. The aim of the present study is the analysis of partially standing surface waves icluding non-linear effects. According to 1st order Stokes theory, synchronous measurements of horizontal and vertical velocity components allow calculation of rate of standing waves (Drevard et al, 2003). In the present study, it is demonstrated that for deep water conditions, partially standing 2nd order Stokes waves induced velocity field is still represented by the 1st order solution for the velocity potential contrary to the surface elevation which exhibits harmonic components. For intermediate water depth, harmonic components appear not only in the surface elevation but also in the velocity fields, but their weight remains much smaller, because of the vertical decreasing wave induced motion. For irregular waves, the influence of the spectrum width on the non-linear effects in the analysis is discussed. Keywords: Wave measurements ; reflection ; non-linear effects Acknowledgements: This work was initiated during the stay of Prof. Ib Svendsen, as invited Professor, at LSEET in autumn 2002. This study is carried out in the framework of the Scientific French National Programmes PNEC ART7 and PATOM. Their financial supports are acknowledged References: Drevard, D., Meuret, A., Rey, V. Piazzola, J. And Dolle, A.. (2002). "Partially reflected waves measurements using Acoustic Doppler Velocimeter (ADV)", Submitted to ISOPE 03, Honolulu, Hawaii, May 2003.
Ultrafast vectorized multispin coding algorithm for the Monte Carlo simulation of the 3D Ising model
NASA Astrophysics Data System (ADS)
Wansleben, Stephan
1987-02-01
A new Monte Carlo algorithm for the 3D Ising model and its implementation on a CDC CYBER 205 is presented. This approach is applicable to lattices with sizes between 3·3·3 and 192·192·192 with periodic boundary conditions, and is adjustable to various kinetic models. It simulates a canonical ensemble at given temperature generating a new random number for each spin flip. For the Metropolis transition probability the speed is 27 ns per updates on a two-pipe CDC Cyber 205 with 2 million words physical memory, i.e. 1.35 times the cycle time per update or 38 million updates per second.
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)
Bonazzola, S.; Frieben, J.; Gourgoulhon, E.
1998-03-01
An analytical scheme and a numerical method in order to study the effects of general relativity on the viscosity driven secular bar mode instability of rapidly rotating stars are presented. The approach consists in perturbing an axisymmetric and stationary configuration and studying its evolution by constructing a series of triaxial quasi-equilibrium configurations. These are obtained by solution of an approximate set of field equations where only the dominant non-axisymmetric terms are taken into account. The progress with respect to our former investigation consists in a higher relativistic order of the non-axisymmetric terms included into the computation, namely the fully three-dimensional treatment of the vector part of the space-time metric tensor as opposed to the scalar part, solely, in the former case. The scheme is applied to rotating stars built on a polytropic equation of state and compared to our previous results. The 3D-vector part turns out to inhibit the symmetry breaking efficiently. Nevertheless, the bar mode instability is still possible for an astrophysically relevant mass of M_ns=3D1.4M_⊙ when a stiff polytropic equation of state with an adiabatic index of gamma =3D2.5 is employed. Triaxial neutron stars may be efficient emitters of gravitational waves and are thus potentially interesting sources for the forthcoming laser interferometric gravitational wave detectors such as LIGO, VIRGO and GEO600. >From a numerical point of view, the solution of the three-dimensional minimal-distortion shift vector equation in spherical coordinates is an important achievement of our code.
NASA Technical Reports Server (NTRS)
Stremel, Paul M.
1995-01-01
A method has been developed to accurately compute the viscous flow in three-dimensional (3-D) enclosures. This method is the 3-D extension of a two-dimensional (2-D) method developed for the calculation of flow over airfoils. The 2-D method has been tested extensively and has been shown to accurately reproduce experimental results. As in the 2-D method, the 3-D method provides for the non-iterative solution of the incompressible Navier-Stokes equations by means of a fully coupled implicit technique. The solution is calculated on a body fitted computational mesh incorporating a staggered grid methodology. In the staggered grid method, the three components of vorticity are defined at the centers of the computational cell sides, while the velocity components are defined as normal vectors at the centers of the computational cell faces. The staggered grid orientation provides for the accurate definition of the vorticity components at the vorticity locations, the divergence of vorticity at the mesh cell nodes and the conservation of mass at the mesh cell centers. The solution is obtained by utilizing a fractional step solution technique in the three coordinate directions. The boundary conditions for the vorticity and velocity are calculated implicitly as part of the solution. The method provides for the non-iterative solution of the flow field and satisfies the conservation of mass and divergence of vorticity to machine zero at each time step. To test the method, the calculation of simple driven cavity flows have been computed. The driven cavity flow is defined as the flow in an enclosure driven by a moving upper plate at the top of the enclosure. To demonstrate the ability of the method to predict the flow in arbitrary cavities, results will he shown for both cubic and curved cavities.
Xue, Zhong; Shen, Dinggang; Davatzikos, Christos
2004-10-01
Finding point correspondence in anatomical images is a key step in shape analysis and deformable registration. This paper proposes an automatic correspondence detection algorithm for intramodality MR brain images of different subjects using wavelet-based attribute vectors (WAVs) defined on every image voxel. The attribute vector (AV) is extracted from the wavelet subimages and reflects the image structure in a large neighborhood around the respective voxel in a multiscale fashion. It plays the role of a morphological signature for each voxel, and our goal is, therefore, to make it distinctive of the respective voxel. Correspondence is then determined from similarities of AVs. By incorporating the prior knowledge of the spatial relationship among voxels, the ability of the proposed algorithm to find anatomical correspondence is further improved. Experiments with MR images of human brains show that the algorithm performs similarly to experts, even for complex cortical structures. PMID:15493695
NASA Astrophysics Data System (ADS)
Kapitza, H.; Eppel, D.
1987-02-01
A conjugate gradient method for solving a 3-D Poisson equation in Cartesian unequally spaced coordinates is tested in concurrence to standard iterative methods. It is found that the tested algorithm is far superior to Red-Black-SOR with optimal parameter. In the conjugate gradient method no relaxation parameter is needed, and there are no restrictions on the number of gridpoints in the three directions. The iteration routine is vectorizable to a large extent by the compiler of a CYBER 205 without any special preparations. Utilizing some special features of vector computers it is completely vectorizable with only minor changes in the code.
Study on 3-D velocity structure of crust and upper mantle in Sichuan-yunnan region, China
Wang, C.; Mooney, W.D.; Wang, X.; Wu, J.; Lou, H.; Wang, F.
2002-01-01
Based on the first arrival P and S data of 4 625 regional earthquakes recorded at 174 stations dispersed in the Yunnan and Sichuan Provinces, the 3-D velocity structure of crust and upper mantle in the region is determined, incorporating with previous deep geophysical data. In the upper crust, a positive anomaly velocity zone exists in the Sichuan basin, whereas a negative anomaly velocity zone exists in the western Sichuan plateau. The boundary between the positive and negative anomaly zones is the Longmenshan fault zone. The images of lower crust and upper mantle in the Longmenshan fault, Xianshuihe fault, Honghe fault and others appear the characteristic of tectonic boundary, indicating that the faults litely penetrate the Moho discontinuity. The negative velocity anomalies at the depth of 50 km in the Tengchong volcanic area and the Panxi tectonic zone appear to be associated with the temperature and composition variations in the upper mantle. The overall features of the crustal and the upper mantle structures in the Sichuan-Yunnan region are the lower average velocity in both crust and uppermost mantle, the large crustal thickness variations, and the existence of high conductivity layer in the crust or/and upper mantle, and higher geothermal value. All these features are closely related to the collision between the Indian and the Asian plates. The crustal velocity in the Sichuan-Yunnan rhombic block generally shows normal.value or positive anomaly, while the negative anomaly exists in the area along the large strike-slip faults as the block boundary. It is conducive to the crustal block side-pressing out along the faults. In the major seismic zones, the seismicity is relative to the negative anomaly velocity. Most strong earthquakes occurred in the upper-mid crust with positive anomaly or normal velocity, where the negative anomaly zone generally exists below.
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)
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
3D Velocity and Hypocentre Distribution About a Cone-Volcano: Mt Taranaki, New Zealand
NASA Astrophysics Data System (ADS)
Sherburn, S.; White, R.
2003-12-01
Mt Taranaki is a 2518 m andesite cone-volcano (last eruption AD1755) within an oil-bearing sedimentary basin approximately 50 km west of the deepest part of the Benioff zone beneath the North Island of New Zealand. It is the most recent of a series of volcanoes that have erupted in the Taranaki region in the last 1.7 million years. Although a permanent six-station seismic network monitors Mt Taranaki for signs of unrest, little is known of the structure at the depths earthquakes occur and magma maybe stored. This information is vital for interpreting precursors to any future eruption. For nine months in 2001-2002, a temporary network of 75 three-component, broadband (0.03 - 50 Hz) seismographs (area c. 100 km by 100 km) was used to collect data to image crustal structure and accurately locate earthquakes in the Taranaki region. Three hundred and eighty-nine earthquakes were located using more than 15,000 phase picks (55% P and 45% S). A joint inversion for 1D Vp, Vs and hypocentres was undertaken using Velest followed by a 3D inversion for Vp, Vp/Vs ratio and hypocentres using Simul2000. The base of the seismogenic zone increases gradually from a depth of 20 km immediately west of Mt Taranaki to 35 km deep 100 km to the east, corresponding to a previously observed increase in crustal thickness. The area close to Mt Taranaki is anomalous in that there are few earthquakes and all are shallower than 10 km. Within the upper 5-10 km of the crust Vp is closely related to surface geology, being high beneath Mt Taranaki, low beneath the surrounding sedimentary basin, and very high to the east of the basin. We present the Vp and Vp/Vs structure and hypocentre distribution of the Taranaki region and discuss features that can be attributed to volcanism at Mt Taranaki and older volcanic centres.
Present-Day 3D Velocity Field of Eastern North America Based on Continuous GPS Observations
NASA Astrophysics Data System (ADS)
Goudarzi, Mohammad Ali; Cocard, Marc; Santerre, Rock
2016-07-01
The Saint Lawrence River valley in eastern Canada was studied using observations of continuously operating GPS (CGPS) stations. The area is one of the most seismically active regions in eastern North America characterized by many earthquakes, which is also subject to an ongoing glacial isostatic adjustment. We present the current three-dimensional velocity field of eastern North America obtained from more than 14 years (9 years on average) of data at 112 CGPS stations. Bernese GNSS and GITSA software were used for CGPS data processing and position time series analysis, respectively. The results show the counterclockwise rotation of the North American plate in the No-Net-Rotation model with the average of 16.8 ± 0.7 mm/year constrained to ITRF 2008. We also present an ongoing uplift model for the study region based on the present-day CGPS observations. The model shows uplift all over eastern Canada with the maximum rate of 13.7 ± 1.2 mm/year and subsidence to the south mainly over northern USA with a typical rate of -1 to -2 mm/year and the minimum value of -2.7 ± 1.4 mm/year. We compared our model with the rate of radial displacements from the ICE-5G model. Both models agree within 0.02 mm/year at the best stations; however, our model shows a systematic spatial tilt compared to ICE-5G. The misfits between two models amount to the maximum relative subsidence of -6.1 ± 1.1 mm/year to the east and maximum relative uplift of 5.9 ± 2.7 mm/year to the west. The intraplate horizontal velocities are radially outward from the centers of maximum uplift and are inward to the centers of maximum subsidence with the typical velocity of 1-1.6 ± 0.4 mm/year that is in agreement with the ICE-5G model to the first order.
Present-Day 3D Velocity Field of Eastern North America Based on Continuous GPS Observations
NASA Astrophysics Data System (ADS)
Goudarzi, Mohammad Ali; Cocard, Marc; Santerre, Rock
2016-03-01
The Saint Lawrence River valley in eastern Canada was studied using observations of continuously operating GPS (CGPS) stations. The area is one of the most seismically active regions in eastern North America characterized by many earthquakes, which is also subject to an ongoing glacial isostatic adjustment. We present the current three-dimensional velocity field of eastern North America obtained from more than 14 years (9 years on average) of data at 112 CGPS stations. Bernese GNSS and GITSA software were used for CGPS data processing and position time series analysis, respectively. The results show the counterclockwise rotation of the North American plate in the No-Net-Rotation model with the average of 16.8 ± 0.7 mm/year constrained to ITRF 2008. We also present an ongoing uplift model for the study region based on the present-day CGPS observations. The model shows uplift all over eastern Canada with the maximum rate of 13.7 ± 1.2 mm/year and subsidence to the south mainly over northern USA with a typical rate of -1 to -2 mm/year and the minimum value of -2.7 ± 1.4 mm/year. We compared our model with the rate of radial displacements from the ICE-5G model. Both models agree within 0.02 mm/year at the best stations; however, our model shows a systematic spatial tilt compared to ICE-5G. The misfits between two models amount to the maximum relative subsidence of -6.1 ± 1.1 mm/year to the east and maximum relative uplift of 5.9 ± 2.7 mm/year to the west. The intraplate horizontal velocities are radially outward from the centers of maximum uplift and are inward to the centers of maximum subsidence with the typical velocity of 1-1.6 ± 0.4 mm/year that is in agreement with the ICE-5G model to the first order.
Two high resolution velocity vector analyzers for cosmic dust particles
NASA Technical Reports Server (NTRS)
Auer, S.
1975-01-01
Two new methods are described to measure velocities and angles of incidence of charged cosmic dust particles with precisions of about 1% and 1 deg, respectively. Both methods employ four one-dimensional position-sensitive detectors in series. The first method utilizes a charge-dividing technique while the second utilizes a time-of-flight technique for determining the position of a particle inside the instrument. The velocity vectors are measured although mechanical interaction between the particle and the instrument is completely avoided. Applications to cosmic dust composition and collection experiments are discussed. The range of radii of measurable particles is from about 0.01 to 100 microns at velocities from 1 to 80 km/sec.
On the Velocity Field and the 3D Structure of the Galactic Soccer Ball Abell 43
NASA Astrophysics Data System (ADS)
Rauch, Thomas; Werner, Klaus; Ercolano, Barbara; Köppen, Joachim
2005-11-01
Planetary nebulae (PNe) and their central stars (CSs) are ideal tools to test evolutionary theory: photospheric properties of their exciting stars give stringent constraints for theoretical predictions of stellar evolution. The nebular abundances display the star's photosphere chemical composition at the time of the nebula's ejection which allows to look back into the history of stellar evolution. More importantly, they even provide a possibility to investigate on the chemical evolution of our Galaxy because most of the nuclear processed material goes back into the interstellar medium via PNe. The recent developments in observation techniques and the new three-dimensional photoionization code MOCASSIN (Ercolano et al. 2003) enable us to analyze PNe properties accurately by the construction of consistent models of PNe and CSs. In addition to PNe imaging and spectroscopy, detailed information about the velocity field within the PNe is a pre-requisite to employ de-projection techniques in modeling the physical structure of the PNe.
A Robust Method to Detect Zero Velocity for Improved 3D Personal Navigation Using Inertial Sensors
Xu, Zhengyi; Wei, Jianming; Zhang, Bo; Yang, Weijun
2015-01-01
This paper proposes a robust zero velocity (ZV) detector algorithm to accurately calculate stationary periods in a gait cycle. The proposed algorithm adopts an effective gait cycle segmentation method and introduces a Bayesian network (BN) model based on the measurements of inertial sensors and kinesiology knowledge to infer the ZV period. During the detected ZV period, an Extended Kalman Filter (EKF) is used to estimate the error states and calibrate the position error. The experiments reveal that the removal rate of ZV false detections by the proposed method increases 80% compared with traditional method at high walking speed. Furthermore, based on the detected ZV, the Personal Inertial Navigation System (PINS) algorithm aided by EKF performs better, especially in the altitude aspect. PMID:25831086
A robust method to detect zero velocity for improved 3D personal navigation using inertial sensors.
Xu, Zhengyi; Wei, Jianming; Zhang, Bo; Yang, Weijun
2015-01-01
This paper proposes a robust zero velocity (ZV) detector algorithm to accurately calculate stationary periods in a gait cycle. The proposed algorithm adopts an effective gait cycle segmentation method and introduces a Bayesian network (BN) model based on the measurements of inertial sensors and kinesiology knowledge to infer the ZV period. During the detected ZV period, an Extended Kalman Filter (EKF) is used to estimate the error states and calibrate the position error. The experiments reveal that the removal rate of ZV false detections by the proposed method increases 80% compared with traditional method at high walking speed. Furthermore, based on the detected ZV, the Personal Inertial Navigation System (PINS) algorithm aided by EKF performs better, especially in the altitude aspect. PMID:25831086
NASA Astrophysics Data System (ADS)
Moustafa, Salli; Févotte, François; Lathuilière, Bruno; Plagne, Laurent
2014-06-01
The past few years have been marked by a noticeable increase in the interest in 3D whole-core heterogeneous deterministic neutron transport solvers for reference calculations. Due to the extremely large problem sizes tackled by such solvers, they need to use adapted numerical methods and need to be efficiently implemented to take advantage of the full computing power of modern systems. As for numerical methods, one possible approach consists in iterating over resolutions of 2D and 1D MOC problems by taking advantage of prismatic geometries. The MICADO solver, developed at EDF R&D, is a parallel implementation of such a method in distributed and shared memory systems. However it is currently unable to use SIMD vectorization to leverage the full computing power of modern CPUs. In this paper, we describe our first effort to support vectorization in MICADO, typically targeting Intel© SSE CPUs. Both the 2D and 1D algorithms are vectorized, allowing for high expected speedups for the whole spatial solver. We present benchmark computations, which show nearly optimal speedups for our vectorized implementation on the TAKEDA case.
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.
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.
Optimum instantaneous impulsive orbital injection to attain a specified asymptotic velocity vector.
NASA Technical Reports Server (NTRS)
Bean, W. C.
1971-01-01
A nalysis of the necessary conditions of Battin for instantaneous orbital injection, with consideration of the uniqueness of his solution, and of the further problem which arises in the degenerate case when radius vector and asymptotic vector are separated by 180 deg. It is shown that when the angular separation between radius vector and asymptotic velocity vector satisfies theta not equal to 180 deg, there are precisely two insertion-velocity vectors which permit attainment of the target asymptotic velocity vector, one yielding posigrade, the other retrograde motion. When theta equals to 180 deg, there is a family of insertion-velocity vectors which permit attainment of a specified asymptotic velocity vector with a unique insertion-velocity vector for every arbitrary orientation of a target unit angular momentum vector.
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)
Bisi, M. M.; Jackson, B. V.; Breen, A. R.; Dorrian, G. D.; Fallows, R. A.; Clover, J. M.; Hick, P. P.
2010-08-01
The European Incoherent SCATter (EISCAT) radar has been used for remote-sensing observations of interplanetary scintillation (IPS) for a quarter of a century. During the April/May 2007 observing campaign, a large number of observations of IPS using EISCAT took place to give a reasonable spatial and temporal coverage of solar wind velocity structure throughout this time during the declining phase of Solar Cycle 23. Many co-rotating and transient features were observed during this period. Using the University of California, San Diego three-dimensional (3-D) time-dependent computer assisted tomography (C.A.T.) solar-wind reconstruction analysis, we show the velocity structure of the inner heliosphere in three dimensions throughout the time interval of 20 April through 20 May 2007. We also compare to white-light remote-sensing observations of an interplanetary coronal mass ejection (ICME) seen by the STEREO Ahead spacecraft inner Heliospheric Imager on 16 May 2007, as well as to in-situ solar-wind measurements taken with near-Earth spacebourne instrumentation throughout this interval. The reconstructions show clear co-rotating regions during this period, and the time-series extraction at spacecraft locations compares well with measurements made by the STEREO, Wind, and ACE spacecraft. This is the first time such clear structures have been revealed using this 3-D technique with EISCAT IPS data as input.
A vectorized algorithm for 3D dynamics of a tethered satellite
NASA Technical Reports Server (NTRS)
Wilson, Howard B.
1989-01-01
Equations of motion characterizing the three dimensional motion of a tethered satellite during the retrieval phase are studied. The mathematical model involves an arbitrary number of point masses connected by weightless cords. Motion occurs in a gravity gradient field. The formulation presented accounts for general functions describing support point motion, rate of tether retrieval, and arbitrary forces applied to the point masses. The matrix oriented program language MATLAB is used to produce an efficient vectorized formulation for computing natural frequencies and mode shapes for small oscillations about the static equilibrium configuration; and for integrating the nonlinear differential equations governing large amplitude motions. An example of time response pertaining to the skip rope effect is investigated.
The Local Stellar Velocity Field via Vector Spherical Harmonics
NASA Astrophysics Data System (ADS)
Makarov, V. V.; Murphy, D. W.
2007-07-01
We analyze the local field of stellar tangential velocities for a sample of 42,339 nonbinary Hipparcos stars with accurate parallaxes, using a vector spherical harmonic formalism. We derive simple relations between the parameters of the classical linear model (Ogorodnikov-Milne) of the local systemic field and low-degree terms of the general vector harmonic decomposition. Taking advantage of these relationships, we determine the solar velocity with respect to the local stars of (VX,VY,VZ)=(10.5,18.5,7.3)+/-0.1 km s-1 not corrected for the asymmetric drift with respect to the local standard of rest. If only stars more distant than 100 pc are considered, the peculiar solar motion is (VX,VY,VZ)=(9.9,15.6,6.9)+/-0.2 km s-1. The adverse effects of harmonic leakage, which occurs between the reflex solar motion represented by the three electric vector harmonics in the velocity space and higher degree harmonics in the proper-motion space, are eliminated in our analysis by direct subtraction of the reflex solar velocity in its tangential components for each star. The Oort parameters determined by a straightforward least-squares adjustment in vector spherical harmonics are A=14.0+/-1.4, B=-13.1+/-1.2, K=1.1+/-1.8, and C=-2.9+/-1.4 km s-1 kpc-1. The physical meaning and the implications of these parameters are discussed in the framework of a general linear model of the velocity field. We find a few statistically significant higher degree harmonic terms that do not correspond to any parameters in the classical linear model. One of them, a third-degree electric harmonic, is tentatively explained as the response to a negative linear gradient of rotation velocity with distance from the Galactic plane, which we estimate at ~-20 km s-1 kpc-1. A similar vertical gradient of rotation velocity has been detected for more distant stars representing the thick disk (z>1 kpc), but here we surmise its existence in the thin disk at z<200 pc. The most unexpected and unexplained term within
The Local Stellar Velocity Field via Vector Spherical Harmonics
NASA Technical Reports Server (NTRS)
Markarov, V. V.; Murphy, D. W.
2007-01-01
We analyze the local field of stellar tangential velocities for a sample of 42,339 nonbinary Hipparcos stars with accurate parallaxes, using a vector spherical harmonic formalism. We derive simple relations between the parameters of the classical linear model (Ogorodnikov-Milne) of the local systemic field and low-degree terms of the general vector harmonic decomposition. Taking advantage of these relationships, we determine the solar velocity with respect to the local stars of (V(sub X), V(sub Y), V(sub Z)) (10.5, 18.5, 7.3) +/- 0.1 km s(exp -1) not corrected for the asymmetric drift with respect to the local standard of rest. If only stars more distant than 100 pc are considered, the peculiar solar motion is (V(sub X), V(sub Y), V(sub Z)) (9.9, 15.6, 6.9) +/- 0.2 km s(exp -1). The adverse effects of harmonic leakage, which occurs between the reflex solar motion represented by the three electric vector harmonics in the velocity space and higher degree harmonics in the proper-motion space, are eliminated in our analysis by direct subtraction of the reflex solar velocity in its tangential components for each star. The Oort parameters determined by a straightforward least-squares adjustment in vector spherical harmonics are A=14.0 +/- 1.4, B=13.1 +/- 1.2, K=1.1 +/- 1.8, and C=2.9 +/- 1.4 km s(exp -1) kpc(exp -1). The physical meaning and the implications of these parameters are discussed in the framework of a general linear model of the velocity field. We find a few statistically significant higher degree harmonic terms that do not correspond to any parameters in the classical linear model. One of them, a third-degree electric harmonic, is tentatively explained as the response to a negative linear gradient of rotation velocity with distance from the Galactic plane, which we estimate at approximately -20 km s(exp -1) kpc(exp -1). A similar vertical gradient of rotation velocity has been detected for more distant stars representing the thick disk (z greater than 1 kpc
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
Volume flow in arteriovenous fistulas using vector velocity ultrasound.
Hansen, Peter Møller; Olesen, Jacob Bjerring; Pihl, Michael Johannes; Lange, Theis; Heerwagen, Søren; Pedersen, Mads Møller; Rix, Marianne; Lönn, Lars; Jensen, Jørgen Arendt; Nielsen, Michael Bachmann
2014-11-01
Volume flow in arteriovenous fistulas for hemodialysis was measured using the angle-independent ultrasound technique Vector Flow Imaging and compared with flow measurements using the ultrasound dilution technique during dialysis. Using an UltraView 800 ultrasound scanner (BK Medical, Herlev, Denmark) with a linear transducer, 20 arteriovenous fistulas were scanned directly on the most superficial part of the fistula just before dialysis. Vector Flow Imaging volume flow was estimated with two different approaches, using the maximum and the average flow velocities detected in the fistula. Flow was estimated to be 242 mL/min and 404 mL/min lower than the ultrasound dilution technique estimate, depending on the approach. The standard deviations of the two Vector Flow Imaging approaches were 175.9 mL/min and 164.8 mL/min compared with a standard deviation of 136.9 mL/min using the ultrasound dilution technique. The study supports that Vector Flow Imaging is applicable for volume flow measurements. PMID:25282482
NASA Astrophysics Data System (ADS)
Tsvetkov, A.; Vityazev, . V.; , Kumkova I. I.
2009-09-01
The paper describes the application of the 3-D vector spherical harmonics (henceforth VSH) to the investigation of stellar kinematics. The VSH technique is suitable for present and future catalogues which contain all three components of velocity vector: proper motions and radial velocities. In general, the VSH allows to detect all the systematic components in the stellar velocity field and does not depend on any model. If some physical model is used, the VSH not only determines the parameters of the model, but detects the systematic components which are beyond the model. The application of the VSH to the Hipparcos data complimented with radial velocities discovers the systematic components which are beyond the linear Ogorodnikov-milne model.
Geometric relations of eye position and velocity vectors during saccades.
Tweed, D; Vilis, T
1990-01-01
Measurements of angular position and velocity vectors of the eye in three human and three monkey subjects showed that: (1) position vectors lie roughly in a single plane, in accordance with Listing's law, between and during saccades; (2) primary position of the eye is often far from the centre of the oculomotor range. (3) saccades have nearly-fixed rotation axes, which tilt out of Listing's plane in a systematic way depending on current eye position. Findings 1 and 3 show that saccadic control signals accurately reflect the properties of three-dimensional rotations, as predicted by a new quaternion model of the saccadic system; models that approximate rotational kinematics using vectorial addition and integration do not predict these findings. PMID:2321357
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.
Kadukova, Maria; Grudinin, Sergei
2016-08-22
Here we address the problem of the assignment of atom types and bond orders in low molecular weight compounds. For this purpose, we have developed a prediction model based on nonlinear Support Vector Machines (SVM), implemented in a KNOwledge-Driven Ligand Extractor called Knodle, a software library for the recognition of atomic types, hybridization states, and bond orders in the structures of small molecules. We trained the model using an excessive amount of structural data collected from the PDBbindCN database. Accuracy of the results and the running time of our method is comparable with other popular methods, such as NAOMI, fconv, and I-interpret. On the popular Labute's benchmark set consisting of 179 protein-ligand complexes, Knodle makes five to six perception errors, NAOMI makes seven errors, I-interpret makes nine errors, and fconv makes 13 errors. On a larger set of 3,000 protein-ligand structures collected from the PDBBindCN general data set (v2014), Knodle and NAOMI have a comparable accuracy of approximately 3.9% and 4.7% of errors, I-interpret made 6.0% of errors, while fconv produced approximately 12.8% of errors. On a more general set of 332,974 entries collected from the Ligand Expo database, Knodle made 4.5% of errors. Overall, our study demonstrates the efficiency and robustness of nonlinear SVM in structure perception tasks. Knodle is available at https://team.inria.fr/nano-d/software/Knodle . PMID:27405533
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.
3D tomographic reconstruction of the internal velocity field of an immiscible drop in a shear flow
NASA Astrophysics Data System (ADS)
Kerdraon, Paul; Dalziel, Stuart B.; Goldstein, Raymond E.; Landel, Julien R.; Peaudecerf, Francois J.
2015-11-01
We study experimentally the internal flow of a drop attached to a flat substrate and immersed in an immiscible shear flow. Transport inside the drop can play a crucial role in cleaning applications. Internal advection can enhance the mass transfer across the drop surface, thus increasing the cleaning rate. We used microlitre water-glycerol drops on a hydrophobic substrate. The drops were spherical and did not deform significantly under the shear flow. An oil phase of relative viscosity 0.01 to 1 was flowed over the drop. Typical Reynolds numbers inside the drops were of the order of 0.1 to 10. Using confocal microscopy, we performed 3D tomographic reconstruction of the flow field in the drop. The in-plane velocity field was measured using micro-PIV, and the third velocity component was computed from incompressibility. To our knowledge, this study gives the first experimental measurement of the three-dimensional internal velocity field of a drop in a shear flow. Numerical simulations and theoretical models published in the past 30 years predict a toroidal internal recirculation flow, for which the entire surface flows streamwise. However, our measurements reveal a qualitatively different picture with a two-lobed recirculation, featuring two stagnation points at the surface and a reverse surface flow closer to the substrate. This finding appears to be independent of Reynolds number and viscosity ratio in the ranges studied; we conjecture that the observed flow is due to the effect of surfactants at the drop surface.
The 3-D strain patterns in Turkey using geodetic velocity fields from the RTK-CORS (TR) network
NASA Astrophysics Data System (ADS)
Kutoglu, Hakan Senol; Toker, Mustafa; Mekik, Cetin
2016-03-01
This study presents our use of GPS data to obtain and quantify the full continuous strain tensor using a 3-D velocity field in Turkey. In this study, GPS velocities improve the estimation of short-term strain tensor fields for determining the seismic hazard of Turkey. The tensorial analysis presents different aspects of deformation, such as the normal and shear strains, including their directions, the compressional and extensional strains. This analysis is appropriate for the characterizing the state of the current seismic deformation. GPS velocity data from continuous measurements (2009-2012) to estimate deformations were processed using the GAMIT/GLOBK software. Using high-rate GPS data from permanent 146 GNSS stations (RTK-CORS-TR network), the strain distribution was determined and interpolated using a biharmonic spline technique. We show the strain field patterns within axial and plane form at several critical locations, and discuss these results within the context of the seismic and tectonic deformation of Turkey. We conclude that the knowledge of the crustal strain patterns provides important information on the location of the main faults and strain accumulation for the hazard assessment. The results show an agreement between the seismic and tectonic strains confirming that there are active crustal deformations in Turkey.
NASA Astrophysics Data System (ADS)
Tronchin, Thibaut; David, Laurent; Farcy, Alain
2015-01-01
The flow around a flapping wing is characterized by an unsteady evolution of three-dimensional vortices, which are one of the main sources of loads. The difficulty in directly measuring such low forces by means of sensors and the need of the characterization of the evolution of the flow have lead to the evaluation of loads using the integral form of the momentum equation. This paper describes methods for evaluating instantaneous loads and three-dimensional pressure fields using 3D3C velocity fields only. An evaluation of the accuracy of these methods using DNS velocity fields is presented. Loads and pressure fields are then calculated using scanning tomography PIV velocity fields, around a NACA 0012 airfoil for a flapping motion in a water tank at a Reynolds number of 1,000. The results suggest a sufficient accuracy of calculated pressure fields for a global analysis of the topology of the flow and for the evaluation of loads by integrating the calculated pressure field over the surface of the wing.
Lu, Louise; Sick, Volker; Frank, Jonathan H.
2013-09-01
The University of Michigan and Sandia National Laboratories collaborated on the initial development of a compact single-camera approach for simultaneously measuring 3-D gasphase velocity and temperature fields at high frame rates. A compact diagnostic tool is desired to enable investigations of flows with limited optical access, such as near-wall flows in an internal combustion engine. These in-cylinder flows play a crucial role in improving engine performance. Thermographic phosphors were proposed as flow and temperature tracers to extend the capabilities of a novel, compact 3D velocimetry diagnostic to include high-speed thermometry. Ratiometric measurements were performed using two spectral bands of laser-induced phosphorescence emission from BaMg2Al10O17:Eu (BAM) phosphors in a heated air flow to determine the optimal optical configuration for accurate temperature measurements. The originally planned multi-year research project ended prematurely after the first year due to the Sandia-sponsored student leaving the research group at the University of Michigan.
A GRASS GIS based Spatio-Temporal Algebra for Raster-, 3D Raster- and Vector Time Series Data
NASA Astrophysics Data System (ADS)
Leppelt, Thomas; Gebbert, Sören
2015-04-01
Enhancing the well known and widely used map algebra proposed by Dr. Charles Dana Tomlin [1] with the time dimension is an ongoing research topic. The efficient processing of large time series of raster, 3D raster and vector datasets, e. g. raster datasets for temperature or precipitations on continental scale, requires a sophisticated spatio-temporal algebra that is capable of handling datasets with different temporal granularities and spatio-temporal extents. With the temporal enabled GRASS GIS [2] and the GRASS GIS Temporal Framework new spatio-temporal data types are available in GRASS GIS 7, called space time datasets. These space time datasets represent time series of raster, 3D raster and vector map layers. Furthermore the temporal framework provides a wide range of functionalities to support the implementation of a temporal algebra. While spatial capabilities of GRASS GIS are used to perform the spatial processing of the time stamped map layers that are registered in a space time dataset, the temporal processing is provided by the GRASS GIS temporal framework that supports time intervals and time instances. Mixing time instance and time intervals as well as gaps, overlapping or inclusion of intervals and instances is possible. Hence this framework allows an arbitrary layout of the time dimension. We implemented two ways to process space time datasets with arbitrary temporal layout, the temporal topology and the granularity based spatio-temporal algebra. The algebra provides the functionality to define complex spatio-temporal topological operators that process time and space in a single expression. The algebra includes methods to select map layers from space time datasets based on their temporal relations, to temporally shift time stamped map layers, to create temporal buffer and to snap time instances of time stamped map layers to create a valid temporal topology. In addition spatio-temporal operations can be evaluated within conditional statements. These
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)
Zeng, Xiangfang; Thurber, Clifford H.; Shelly, David R.; Harrington, Rebecca M.; Cochran, Elizabeth S.; Bennington, Ninfa L.; Peterson, Dana; Guo, Bin; McClement, Kara
2016-09-01
To refine the 3-D seismic velocity model in the greater Parkfield, California region, a new data set including regular earthquakes, shots, quarry blasts and low-frequency earthquakes (LFEs) was assembled. Hundreds of traces of each LFE family at two temporary arrays were stacked with time-frequency domain phase weighted stacking method to improve signal-to-noise ratio. We extend our model resolution to lower crustal depth with LFE data. Our result images not only previously identified features but also low velocity zones (LVZs) in the area around the LFEs and the lower crust beneath the southern Rinconada Fault. The former LVZ is consistent with high fluid pressure that can account for several aspects of LFE behaviour. The latter LVZ is consistent with a high conductivity zone in magnetotelluric studies. A new Vs model was developed with S picks that were obtained with a new autopicker. At shallow depth, the low Vs areas underlie the strongest shaking areas in the 2004 Parkfield earthquake. We relocate LFE families and analyse the location uncertainties with the NonLinLoc and tomoDD codes. The two methods yield similar results.
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)
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.
NASA Technical Reports Server (NTRS)
Owen, Albert K.
1993-01-01
The mathematical relations between the measured velocity fields for the same compressor rotor flow field resolved by two fringe type laser anemometers at different observational locations are developed in this report. The relations allow the two sets of velocity measurements to be combined to produce a total velocity vector field for the compressor rotor. This report presents the derivation of the mathematical relations, beginning with the specification of the coordinate systems and the velocity projections in those coordinate systems. The vector projections are then transformed into a common coordinate system. The transformed vector coordinates are then combined to determine the total velocity vector. A numerical example showing the solution procedure is included.
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
NASA Astrophysics Data System (ADS)
José Ramón, Mª; Pueyo, Emilio L.; Briz, José Luis; Caumon, Guillaume; Fernández, Óscar; Ciria, José Carlos; Pocovi, Andrés; Ros, Luis H.
2013-04-01
Three-dimensional reconstructions of the underground involve the integration of discrete and heterogeneous datasets and have significant socio- economic implications. The problem arises when there are limited data to build 3D models or when deformation processes are complex; these reasons inspired the development of restoration methods to validate subsurface reconstructions. The restoration is based on the application of simple geometric (or mechanic) laws that help reduce the uncertainty and increase geomodel accuracy. Apart from mechanical approaches, geometric methods are based on the initial assumption of global conservation of volume during deformation in addition to the paleo-horizontality of the stratigraphic horizons in the undeformed stage. The problem is that the bedding plane cannot be used as a three-dimensional reference system, because a single vector defines it and additional constraints are required. This is particularly important when dealing with complex structures, such as non-cylindrical structures and the superposition of non-coaxial geometries. In this context, paleomagnetism (known in both the deformed and undeformed stages) can contribute to building a more complete reference system and to reducing the uncertainty in restoration processes. The use of paleomagnetism in restoration tools was suggested in the early 1990's and only a few quantitative map-view applications have been developed since then. In this contribution, we introduce the two first surface restoration methods that use paleomagnetic vectors as a primary reference. The first one is a simple geometric approach based on the piecewise restoration of a triangulated surface into which paleomagnetic variables can be easily incorporated. It is valid for complexly folded structures. The surface is modelled by a mesh and the method starts from a pin-element. Triangles are laid flat and then fitted together to minimize distances between common vertices and paleomagnetic error. However
NASA Astrophysics Data System (ADS)
Zhai, Cuili; Zhang, Ting
2015-09-01
In this article, we consider the global well-posedness to the 3-D incompressible inhomogeneous Navier-Stokes equations with a class of large velocity. More precisely, assuming a 0 ∈ B˙ q , 1 /3 q ( R 3 ) and u 0 = ( u0 h , u0 3 ) ∈ B˙ p , 1 - 1 + /3 p ( R 3 ) for p, q ∈ (1, 6) with sup ( /1 p , /1 q ) ≤ /1 3 + inf ( /1 p , /1 q ) , we prove that if C a↑0↑ B˙q1/3 q α (↑u0 3↑ B˙ p , 1 - 1 + /3 p/μ + 1 ) ≤ 1 , /C μ (↑u0 h↑ B˙ p , 1 - 1 + /3 p + ↑u03↑ B˙ p , 1 - 1 + /3 p 1 - α ↑u0h↑ B˙ p , 1 - 1 + /3 p α) ≤ 1 , then the system has a unique global solution a ∈ C ˜ ( [ 0 , ∞ ) ; B˙ q , 1 /3 q ( R 3 ) ) , u ∈ C ˜ ( [ 0 , ∞ ) ; B˙ p , 1 - 1 + /3 p ( R 3 ) ) ∩ L 1 ( R + ; B˙ p , 1 1 + /3 p ( R 3 ) ) . It improves the recent result of M. Paicu and P. Zhang [J. Funct. Anal. 262, 3556-3584 (2012)], where the exponent form of the initial smallness condition is replaced by a polynomial form.
Wu, Yingchun; Wu, Xuecheng; Yao, Longchao; Gréhan, Gérard; Cen, Kefa
2015-03-20
The 3D measurement of the particles in a gas-solid pipe flow is of great interest, but remains challenging due to curved pipe walls in various engineering applications. Because of the astigmatism induced by the pipe, concentric ellipse fringes in the hologram of spherical particles are observed in the experiments. With a theoretical analysis of the particle holography by an ABCD matrix, the in-focus particle image can be reconstructed by the modified convolution method and fractional Fourier transform. Thereafter, the particle size, 3D position, and velocity are simultaneously measured by digital holographic particle tracking velocimetry (DHPTV). The successful application of DHPTV to the particle size and 3D velocity measurement in a glass pipe's flow can facilitate its 3D diagnostics. PMID:25968543
NASA Astrophysics Data System (ADS)
Asmerom, Biniam Beyene
Three-dimensional velocity structures and seismicity patterns have been studied across the Andean Foreland of San Juan Argentina using data acquired by PANDA deployment. Distinct velocity variations are revealed between Precordillera in the west and Pie de Palo in the east. The low velocity anomaly beneath Precordillera is associated with the presence of thick sedimentary rocks and thick sediment cover of Matagusanos valley. Similarly, the high velocity anomaly east of Eastern Precordillera is correlated with the presence of basement rocks. These anomalies are observed from the station corrections of Joint Hypocentral Determination (JHD) analysis. A northeast trending west dipping high velocity anomaly is imaged beneath the southern half of Pie de Palo. This anomaly represents a Grenvillian suture zone formed when Pie de Palo collided with the Precordillera. Relocated seismicity using 3-D Vp and Vs models obtained in this study revealed crustal scale buried faults beneath the Eastern Precordillera and Sierra Pie de Palo. The fault defined by the seismicity extend down to a depth of ˜ 40 km and ~35 km beneath Precordillera and Pie de Palo, respectively, defining the lower bound of the brittle to ductile transition of the crust. These results confirm that present day active crustal thickening involves the entire crust in the tectonic process and results in thick-skinned deformation beneath both the Eastern Precordillera and Pie de Palo. Based on the seismicity pattern, geomorphology, and velocity structures, Sierra Pie de Palo, a basement uplift block, can be divided into two separate semi-blocks separated by a northeast trending fracture zone. The northern block is characterized by a well-defined west dipping fault and low Vp/Vs ratio particularly at a depth of 12 to 16 km, while the southern block shows a poorly-defined east dipping fault with high Vp/Vs ratio at a depth of 20 to 26 km. Spatial distribution of the well-relocated crustal earthquakes along these
NASA Astrophysics Data System (ADS)
Perkovic, D.; Toporkov, J. V.; Sletten, M. A.; Farquharson, G.; Frasier, S. J.; Marmorino, G. O.; Judd, K. P.
2004-12-01
The Dual Beam Interferometer (DBI) developed by University of Massachusetts (UMass) consists of two C-band along-track interferometric synthetic aperture radars (ATI-SAR). The beams of this airborne system are squinted 20 degrees forward and aft of broadside allowing surface vector velocity estimation in a single aircraft pass. The instrument has been deployed several times over the period of last two years off coastal areas of Florida on a National Oceanic and Atmospheric Administration's (NOAA) WP-3D plane in collaboration with the Naval Research Laboratory (NRL). During 2002-2003 the instrument has undergone a series of tests and engineering flights with the August 2003 data producing the first interferogram. March 2004 flights were mainly focused on the western boundary of the Gulf Stream off Cape Canaveral, Florida. Simultaneous imagery of the sea-surface temperature field were obtained using NRL's Infrared (IR) camera, which was mounted in belly of the aircraft. Multiple passes over the Gulf Stream were made under a range of environmental conditions and viewing geometries relative to the Gulf Stream current and the wind. Additional flights were made over the barrier islands west of Ft. Meyers, Florida at times of near maximum ebb tidal flow. This paper will present initial estimates of the surface vector velocities for each area. The Gulf Stream velocity estimates show a current maximum of 1.5 m/s across the edge of the Stream which is consistent with estimates of the current from IR imagery using feature tracking. Estimates of the flow between the barrier islands are of the order of 1.5 to 2 m/s, which agrees well with the predicted tidal flow.
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)
Kanellopoulos, V. N.; Webb, J. P.
1993-03-01
A 3D vector analysis of plane wave scattering by a metallic sphere using finite elements and Absorbing Boundary Conditions (ABCs) is presented. The ABCs are applied on the outer surface that truncates the infinitely extending domain. Mixed order curvilinear covariantprojection elements are used to avoid spurious corruptions. The second order ABC is superior to the first at no extra computational cost. The errors due to incomplete absorption decrease as the outer surface is moved further away from the scatterer. An error of about 1% in near-field values was obtained with the second order ABC, when the outer surface was less than half a wavelength from the scatterer. Une analyse tridimensionnelle vectorielle de la diffusion d'onde plane sur une sphère métallique utilisant des éléments finis et des Conditions aux Limites Absorbantes (CLA) est présentée. Les CLA sont appliquées sur la surface exteme tronquant le domaine s'étendant à l'infini. Des éléments curvilignes mixtes utilisant des projections covariantes sont utilisés pour éviter des solutions parasites. La CLA de second ordre est supérieure à celle de premier ordre sans effort de calcul additionnel. Les erreurs dues à l'absorption incomplète décroissent à mesure que l'on déplace la surface externe à une distance croissante du diffuseur. Un taux d'erreur d'environ 1 % dans les valeurs du champ proche a été obtenu avec les CLA de second ordre lorsque la surface externe était placée à une distance inférieure à une demi-longueur de la source de diffusion.
NASA Astrophysics Data System (ADS)
Akuhara, T.; Mochizuki, K.; Nakahigashi, K.; Yamada, T.; Shinohara, M.; Sakai, S.; Kanazawa, T.; Uehira, K.; Shimizu, H.
2013-12-01
The Philippine Sea Plate subducts beneath the Eurasian Plate at a rate of ~4 cm/year along the Nankai Trough, southwest of Japan. Around the Kii Peninsula, the rupture boundary of the historical Tonankai and Nankai large earthquakes is located, and previous researches have revealed along-strike segmentation of hypocenters [Mochizuki et al., 2010], P-wave anisotropy [Ishise et al., 2009], low frequency earthquake (LFE) distribution [e.g., Obara, 2010] and subduction depth of the Philippine Sea (PHS) Plate, or there may exist a split in the PHS Plate [Ide et al., 2010]. To investigate such segmentation, in our previous work we determined 3-D velocity structure and hypocenters using P- and S-wave arrival times of earthquakes recorded by both ocean bottom seismometers (OBSs) that were deployed from 2003 to 2007 and on-land stations [Akuhara et al., 2013]. As a result, it was discovered that Vp/Vs ratio is also segmented within the oceanic crust and at the bottom of the overriding plate, which coincides with the LFE distribution: segment A is located along the Kii Channel, segment B around the western Kii Peninsula, and segment C around the eastern Kii Peninsula. In segment B, Vp/Vs ratio is low within the oceanic crust and LFE cluster characterized by an anomalously small amount of cumulative slip, compared to the other LFE clusters around the Kii Peninsula, is located [Obara, 2010]. The difference of Vp/Vs ratio and LFE activity among segments were interpreted as difference of pore fluid pressure. In fact, similar segmentation can be seen in hypocenters: Segment A with concentrated seismicity in the oceanic mantle, segment B with that in the oceanic crust, and segment C with little seismicity. To derive characteristic patterns of the hypocenters, we conducted a cluster analysis of earthquakes based on waveform similarity represented by cross-correlation coefficients (CCs) [e.g., Cattaneo, 1999], in which we took varying structural site effects among the OBS stations
NASA Astrophysics Data System (ADS)
GarcíA-Yeguas, Araceli; Koulakov, Ivan; IbáñEz, Jesús M.; Rietbrock, A.
2012-09-01
We present a high resolution 3 dimensional (3D) P wave velocity model for Tenerife Island, Canaries, covering the top of Teide volcano (3,718 m a.s.l.) down to around 8 km below sea level (b.s.l). The tomographic inversion is based on a large data set of travel times obtained from a 3D active seismic experiment using offshore shots (air guns) recorded at more than 100 onshore seismic stations. The obtained seismic velocity structure is strongly heterogeneous with significant (up to 40%) lateral variations. The main volcanic structure of the Las Cañadas-Teide-Pico Viejo Complex (CTPVC) is characterized by a high P wave velocity body, similar to many other stratovolcanoes. The presence of different high P wave velocity regions inside the CTPVC may be related to the geological and volcanological evolution of the system. The presence of high P wave velocities at the center of the island is interpreted as evidence for a single central volcanic source for the formation of Tenerife. Furthermore, reduced P wave velocities are found in a small confined region in CTPVC and are more likely related to hydrothermal alteration, as indicated by the existence of fumaroles, than to the presence of a magma chamber beneath the system. In the external regions, surrounding CTPVC a few lower P wave velocity regions can be interpreted as fractured zones, hydrothermal alterations, porous materials and thick volcaniclastic deposits.
PLOT3D/AMES, APOLLO UNIX VERSION USING GMR3D (WITH TURB3D)
NASA Technical Reports Server (NTRS)
Buning, P.
1994-01-01
PLOT3D is an interactive graphics program designed to help scientists visualize computational fluid dynamics (CFD) grids and solutions. Today, supercomputers and CFD algorithms can provide scientists with simulations of such highly complex phenomena that obtaining an understanding of the simulations has become a major problem. Tools which help the scientist visualize the simulations can be of tremendous aid. PLOT3D/AMES offers more functions and features, and has been adapted for more types of computers than any other CFD graphics program. Version 3.6b+ is supported for five computers and graphic libraries. Using PLOT3D, CFD physicists can view their computational models from any angle, observing the physics of problems and the quality of solutions. As an aid in designing aircraft, for example, PLOT3D's interactive computer graphics can show vortices, temperature, reverse flow, pressure, and dozens of other characteristics of air flow during flight. As critical areas become obvious, they can easily be studied more closely using a finer grid. PLOT3D is part of a computational fluid dynamics software cycle. First, a program such as 3DGRAPE (ARC-12620) helps the scientist generate computational grids to model an object and its surrounding space. Once the grids have been designed and parameters such as the angle of attack, Mach number, and Reynolds number have been specified, a "flow-solver" program such as INS3D (ARC-11794 or COS-10019) solves the system of equations governing fluid flow, usually on a supercomputer. Grids sometimes have as many as two million points, and the "flow-solver" produces a solution file which contains density, x- y- and z-momentum, and stagnation energy for each grid point. With such a solution file and a grid file containing up to 50 grids as input, PLOT3D can calculate and graphically display any one of 74 functions, including shock waves, surface pressure, velocity vectors, and particle traces. PLOT3D's 74 functions are organized into
PLOT3D/AMES, APOLLO UNIX VERSION USING GMR3D (WITHOUT TURB3D)
NASA Technical Reports Server (NTRS)
Buning, P.
1994-01-01
PLOT3D is an interactive graphics program designed to help scientists visualize computational fluid dynamics (CFD) grids and solutions. Today, supercomputers and CFD algorithms can provide scientists with simulations of such highly complex phenomena that obtaining an understanding of the simulations has become a major problem. Tools which help the scientist visualize the simulations can be of tremendous aid. PLOT3D/AMES offers more functions and features, and has been adapted for more types of computers than any other CFD graphics program. Version 3.6b+ is supported for five computers and graphic libraries. Using PLOT3D, CFD physicists can view their computational models from any angle, observing the physics of problems and the quality of solutions. As an aid in designing aircraft, for example, PLOT3D's interactive computer graphics can show vortices, temperature, reverse flow, pressure, and dozens of other characteristics of air flow during flight. As critical areas become obvious, they can easily be studied more closely using a finer grid. PLOT3D is part of a computational fluid dynamics software cycle. First, a program such as 3DGRAPE (ARC-12620) helps the scientist generate computational grids to model an object and its surrounding space. Once the grids have been designed and parameters such as the angle of attack, Mach number, and Reynolds number have been specified, a "flow-solver" program such as INS3D (ARC-11794 or COS-10019) solves the system of equations governing fluid flow, usually on a supercomputer. Grids sometimes have as many as two million points, and the "flow-solver" produces a solution file which contains density, x- y- and z-momentum, and stagnation energy for each grid point. With such a solution file and a grid file containing up to 50 grids as input, PLOT3D can calculate and graphically display any one of 74 functions, including shock waves, surface pressure, velocity vectors, and particle traces. PLOT3D's 74 functions are organized into
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)
Lee, Myungeun; Kim, Jong Hyo; Park, Moon Ho; Kim, Ye-Hoon; Seong, Yeong Kyeong; Cho, Baek Hwan; Woo, Kyoung-Gu
2014-03-01
In this study, we propose a computer-aided classification scheme of liver tumor in 3D ultrasound by using a combination of deformable model segmentation and support vector machine. For segmentation of tumors in 3D ultrasound images, a novel segmentation model was used which combined edge, region, and contour smoothness energies. Then four features were extracted from the segmented tumor including tumor edge, roundness, contrast, and internal texture. We used a support vector machine for the classification of features. The performance of the developed method was evaluated with a dataset of 79 cases including 20 cysts, 20 hemangiomas, and 39 hepatocellular carcinomas, as determined by the radiologist's visual scoring. Evaluation of the results showed that our proposed method produced tumor boundaries that were equal to or better than acceptable in 89.8% of cases, and achieved 93.7% accuracy in classification of cyst and hemangioma.
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)
Kumar, Prashant; Deb, Sanjib K.; Kishtawal, C. M.; Pal, P. K.
2016-01-01
The Weather Research and Forecasting (WRF) model and its three-dimensional variational data assimilation system are used in this study to assimilate the INSAT-3D, a recently launched Indian geostationary meteorological satellite derived from atmospheric motion vectors (AMVs) over the South Asian region during peak Indian summer monsoon month (i.e., July 2014). A total of four experiments were performed daily with and without assimilation of INSAT-3D-derived AMVs and the other AMVs available through Global Telecommunication System (GTS) for the entire month of July 2014. Before assimilating these newly derived INSAT-3D AMVs in the numerical model, a preliminary evaluation of these AMVs is performed with National Centers for Environmental Prediction (NCEP) final model analyses. The preliminary validation results show that root-mean-square vector difference (RMSVD) for INSAT-3D AMVs is ˜3.95, 6.66, and 5.65 ms-1 at low, mid, and high levels, respectively, and slightly more RMSVDs are noticed in GTS AMVs (˜4.0, 8.01, and 6.43 ms-1 at low, mid, and high levels, respectively). The assimilation of AMVs has improved the WRF model of produced wind speed, temperature, and moisture analyses as well as subsequent model forecasts over the Indian Ocean, Arabian Sea, Australia, and South Africa. Slightly more improvements are noticed in the experiment where only the INSAT-3D AMVs are assimilated compared to the experiment where only GTS AMVs are assimilated. The results also show improvement in rainfall predictions over the Indian region after AMV assimilation. Overall, the assimilation of INSAT-3D AMVs improved the WRF model short-range predictions over the South Asian region as compared to control experiments.
Multifractal Analysis of Velocity Vector Fields and a Continuous In-Scale Cascade Model
NASA Astrophysics Data System (ADS)
Fitton, G.; Tchiguirinskaia, I.; Schertzer, D.; Lovejoy, S.
2012-04-01
In this study we have compared the multifractal analyses of small-scale surface-layer wind velocities from two different datasets. The first dataset consists of six-months of wind velocity and temperature measurements at the heights 22, 23 and 43m. The measurements came from 3D sonic anemometers with a 10Hz data output rate positioned on a mast in a wind farm test site subject to wake turbulence effects. The location of the test site (Corsica, France) meant the large scale structures were subject to topography effects that therefore possibly caused buoyancy effects. The second dataset (Germany) consists of 300 twenty minute samples of horizontal wind velocity magnitudes simultaneously recorded at several positions on two masts. There are eight propeller anemometers on each mast, recording velocity magnitude data at 2.5Hz. The positioning of the anemometers is such that there are effectively two grids. One grid of 3 rows by 4 columns and a second of 5 rows by 2 columns. The ranges of temporal scale over which the analyses were done were from 1 to 103 seconds for both datasets. Thus, under the universal multifractal framework we found both datasets exhibit parameters α ≈ 1.5 and C1 ≈ 0.1. The parameters α and C1, measure respectively the multifractality and mean intermittency of the scaling field. A third parameter, H, quantifies the divergence from conservation of the field (e.g. H = 0 for the turbulent energy flux density). To estimate the parameters we used the ratio of the scaling moment function of the energy flux and of the velocity increments. This method was particularly useful when estimating the parameter α over larger scales. In fact it was not possible to obtain a reasonable estimate of alpha using the usual double trace moment method. For each case the scaling behaviour of the wind was almost isotropic when the scale ranges remained close to the sphero-scale. For the Corsica dataset this could be seen by the agreement of the spectral exponents of
NASA Astrophysics Data System (ADS)
Zhang, M. Q.
1989-09-01
A new Monte Carlo algorithm for 3D Kawasaki spin-exchange simulations and its implementation on a CDC CYBER 205 is presented. This approach is applicable to lattices with sizes between 4×4×4 and 256×L2×L3 ((L2+2)(L3+4)/4⩽65535) and periodic boundary conditions. It is adjustable to various kinetic models in which the total magnetization is conserved. Maximum speed on 10 million steps per second can be reached for 3-D Ising model with Metropolis rate.
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)
Louie, J. N.; Pancha, A.; Pullammanappallil, S. K.
2014-12-01
Refraction microtermor routinely assesses 1D and 2D velocity-depth profiles to shallow depths of approximately 100 m primarily for engineering applications. Estimation of both shallow and deep (>100 m) shear-velocity structure are key elements in the assessment of urban areas for potential earthquake ground shaking, damage, and the calibration of recorded ground motions. Three independent studies investigated the ability of the refraction microtremor technology to image deep velocity structure, to depths exceeding 1 km (Deep ReMi). In the first study, we were able to delineate basin thicknesses of up to 900 m and the deep-basin velocity structure beneath the Reno-area basin. Constraints on lateral velocity changes in 3D as well as on velocity profiles extended down to 1500 m, and show a possible fault offset. This deployment used 30 stand-alone wireless instruments mated to 4.5 Hz geophones, along each of five arrays 2.9 to 5.8 km long. Rayleigh-wave dispersion was clear at frequencies as low as 0.5 Hz using 120 sec ambient urban noise records. The results allowed construction of a 3D velocity model, vetted by agreement with gravity studies. In a second test, a 5.8 km array delimited the southern edge of the Tahoe Basin, with constraints from gravity. There, bedrock depth increased by 250 m in thickness over a distance of 1600 m, with definition of the velocity of the deeper basin sediments. The third study delineated the collapse region of an underground nuclear explosion within a thick sequence of volcanic extrusives, using a shear-wave minivibe in a radial direction, and horizontal geophones. Analysis showed the cavity extends to 620 m depth, with a width of 180 m and a height of 220 m. Our results demonstrate that deep velocity structure can be recovered using ambient noise. This technique offers the ability to define 2D and 3D structural representations essential for seismic hazard evaluation.
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.
NASA Astrophysics Data System (ADS)
Bai, Chao-ying; Huang, Guo-jiao; Li, Xing-wang; Greenhalgh, Stewart
2015-10-01
Traditionally, traveltime tomography entails inversion of either the velocity field and the reflector geometry sequentially, or the velocity field and the hypocenter locations simultaneously or in a cascaded fashion, but seldom are all three types (velocities, geometry of reflectors, and source locations) updated simultaneously because of the compromise between the different classes of model variable and the lack of different seismic phases to constrain these variables. By using a state-of-the-art ray-tracing algorithm for the first and later arrivals combined with a popular linearized inversion solver, it is possible to simultaneously recover the three classes of model variables. In the work discussed in this paper we combined the multistage irregular shortest-path ray-tracing algorithm with a subspace inversion solver to achieve simultaneous inversion of multi-class variables, using arrival times for different phases to concurrently obtain the velocity field, the reflector shapes, and the hypocenter locations. Simulation and comparison tests for two sets of source-receiver arrangements (one the ideal case and the other an approximated real case) indicate that the combined triple-class inversion algorithm is capable of obtaining nearly the same results as the double-class affect inversion scheme (velocity and reflector geometry, or velocity and source locations) even if a lower ray density and irregular source-receiver geometry are used to simulate the real situation. In addition, the new simultaneous inversion method is not sensitive to a modest amount of picking error in the traveltime data and reasonable uncertainty in earthquake hypocenter locations, which shows it to be a feasible and promising approach in real applications.
NASA Astrophysics Data System (ADS)
Bremner, P. M.; Panning, M. P.; Russo, R. M.; Mocanu, V. I.; Stanciu, A. C.; Torpey, M. E.; Hongsresawat, S.; VanDecar, J. C.
2014-12-01
We present the latest 3-D isotropic crustal velocity model beneath central Idaho and eastern Oregon. We produced the velocity model from vertical component Rayleigh wave group and phase velocity measurements on data from the IDaho/ORegon (IDOR) Passive seismic network, 86 broadband seismic stations, using ambient noise tomography and the methods of Gallego et. al (2010) and Lin et. al (2008). We calculated inter-station group/phase velocities in narrow frequency bands from travel-time measurements of the stacked cross-correlations (bandpass filtered between 2 and 30 seconds), which we used to invert for velocity structure beneath the network. Goals of our work include refining models of crustal structure in the accreted Blue Mountain terranes in the western study area; determining the depth extent of the Salmon River Suture/West Idaho Shear Zone (WISZ), which crosses north-south through the middle of the network; determining the architecture of the Idaho batholith, an extensive largely crustal-derived pluton; and examining the nature of the autochthonous (?) North American crust and lithosphere beneath and east of the batholith. We derived Rayleigh wave group and phase velocity maps for each frequency band using the damped least-squares inversion method of Tarantola (2005), and then jointly inverted for velocity with depth. Moho depths are prescribed in the joint inversions based on receiver functions, also from the IDOR seismic data, and provides a starting crustal velocity model. 3-D checkerboard resolution tests indicate lateral resolution of better than 40 km. Preliminary results show higher S wave velocities in the western study area, and lower velocities in the lower crust on the east side of the network, consistent with Basin-and-Range style extension there. A tabular velocity anomaly juxtaposing higher above lower seismic velocities dips shallow west in the midcrust on the west side of the network.
Group velocity locked vector dissipative solitons in a high repetition rate fiber laser.
Luo, Yiyang; Li, Lei; Liu, Deming; Sun, Qizhen; Wu, Zhichao; Xu, Zhilin; Tang, Dingyuan; Fu, Songnian; Zhao, Luming
2016-08-01
Vectorial nature of dissipative solitons (DSs) with high repetition rate is studied for the first time in a normal-dispersion fiber laser. Despite the fact that the formed DSs are strongly chirped and the repetition rate is greater than 100 MHz, polarization locked and polarization rotating group velocity locked vector DSs can be formed under 129.3 MHz fundamental mode-locking and 258.6 MHz harmonic mode-locking of the fiber laser, respectively. The two orthogonally polarized components of these vector DSs possess distinctly different central wavelengths and travel together at the same group velocity in the laser cavity, resulting in a gradual spectral edge and small steps on the optical spectrum, which can be considered as an auxiliary indicator of the group velocity locked vector DSs. Moreover, numerical simulations well confirm the experimental observations and further reveal the impact of the net cavity birefringence on the properties of the formed vector DSs. PMID:27505834
NASA Technical Reports Server (NTRS)
Wang, Ren H.
1991-01-01
A method of combined use of magnetic vector potential (MVP) based finite element (FE) formulations and magnetic scalar potential (MSP) based FE formulations for computation of three-dimensional (3D) magnetostatic fields is developed. This combined MVP-MSP 3D-FE method leads to considerable reduction by nearly a factor of 3 in the number of unknowns in comparison to the number of unknowns which must be computed in global MVP based FE solutions. This method allows one to incorporate portions of iron cores sandwiched in between coils (conductors) in current-carrying regions. Thus, it greatly simplifies the geometries of current carrying regions (in comparison with the exclusive MSP based methods) in electric machinery applications. A unique feature of this approach is that the global MSP solution is single valued in nature, that is, no branch cut is needed. This is again a superiority over the exclusive MSP based methods. A Newton-Raphson procedure with a concept of an adaptive relaxation factor was developed and successfully used in solving the 3D-FE problem with magnetic material anisotropy and nonlinearity. Accordingly, this combined MVP-MSP 3D-FE method is most suited for solution of large scale global type magnetic field computations in rotating electric machinery with very complex magnetic circuit geometries, as well as nonlinear and anisotropic material properties.
Analysis of the three-dimensional stellar velocity field using vector spherical functions
NASA Astrophysics Data System (ADS)
Vityazev, V. V.; Tsvetkov, A. S.
2009-02-01
We apply vector spherical functions to problems of stellar kinematics. Using these functions allows all of the systematic components in the stellar velocity field to be revealed without being attached to a specific physical model. Comparison of the theoretical decomposition coefficients of the equations for a particular kinematical model with observational data can provide precise information about whether the model is compatible with the observations and can reveal systematic components that are not described by this model. The formalism of vector spherical functions is particularly well suited for analyzing the present and future (e.g., GAIA) catalogs containing all three velocity vector components: the propermotions in both coordinates and the radial velocity. We show that there are systematic components in the proper motions of Hipparcos stars that cannot be interpreted in terms of the linear Ogorodnikov-Milne model. The same result is also confirmed by an analysis of the radial velocities for these stars.
NASA Astrophysics Data System (ADS)
Judenherc, S.; Zollo, A.; Auger, E.; Boschi, L.; Satriano, C.; Serapis Working Group
2003-04-01
In September 2001, the SERAPIS project was carried out as an extended active seismic survey in the gulfs of Naples and Pozzuoli. A dense array of 60 three-component on-land stations and 72 sea bottom seismographs (OBS) have been deployed to record more than 5000 air gun shots at a spacing of about 125~m. As a preliminary analysis, the first P-arrival times of a the small offset data in the central part of the region has been inverted using the codes of H.M. Benz. The linearized iterative inversion of 38000 arrival times provided a >80% variance reduction with a node spacing of 250m. At the first order, our model shows a 2-layer structure : low velocity volcanic sediments (2.5-3.5km/s) lying on an inclined high velocity limestone platform (>6km/s). The caldera itself is very well identified, the rim is characterized by a 500-1000m upward shift of the velocity isolines. The whole dataset is expected to provide a wider image with the same resolution (250m). It includes the detailed shape of the refractor beneath the caldera as well as its irregularities out of the bay which have been observered in the seismic sections.
NASA Astrophysics Data System (ADS)
Tarancioglu, Adil; Kocaoglu, Argun H.; Ozalaybey, Serdar
2014-05-01
The objective of this study is to investigate the local seismicity and obtain a detailed three-dimensional crustal velocity structure beneath the Marmara Sea in an area surrounding the North Anatolian Fault Zone (NAFZ) by tomographic inversion using both controlled-source (air-gun) and earthquake data. The tomographic inversion is carried out by using the local earthquake tomography code SIMUL2000. Two sets of seismological data, collected in 2006 (EOSMARMARA experiment) and 2001 (SEISMARMARA experiment), are re-processed and used in this study. A total of 441 high quality earthquakes and 452 air-gun shots recorded by a total of 53 Ocean Bottom Seismometers (OBS) are selected for the simultaneous inversion for velocity and hypocentral parameters. The OBS location and time-drift errors are identified from air-gun shot records by a grid search method and required corrections are made on the travel time data. The initial (reference) velocity model and earthquake locations required for the three dimensional tomographic inversion are derived from the one-dimensional velocity model obtained by using the VELEST algorithm in which a subset of earthquakes are selected such that phase readings were made by at least five stations and maximum azimuthal gap was 180o. The inversion results are checked for initial model dependence and the effect of damping factor. The reliability of the results is also evaluated in terms of derivative-weighted-sum, resolution-diagonal-elements values and checkerboard tests. The hypocenter locations of the local earthquakes have been remarkably improved by the three-dimensional velocity model obtained from the tomographic inversion. The three-dimensional velocity model shows that the Tekirdag, Central and Cinarcik Basins are characterized generally by lower Vp (3.0 - 3.5 km/s) values and most of the earthquakes across these regions are located at the depths of 10 to 17 km, about 5 km deeper than those obtained from the one-dimensional reference
NASA Astrophysics Data System (ADS)
Woollam, Jack; Fuenzallida, Amaya; Garth, Tom; Rietbrock, Andreas; Ruiz, Sergio; Tavera, Hernando
2016-04-01
Seismic velocity tomography is one of the key tools in Earth sciences to image the physical properties of the subsurface. In recent years significant advances have been made to image the Chilean subductions zone, especially in the area of the 2010 M8.8 Maule earthquake (e.g. Hicks et al., 2014), providing much needed physical constraints for earthquakes source inversions and rupture models. In 2014 the M8.2 Iquique earthquake struck the northern part of the Chilean subduction zone in close proximity to the Peruvian boarder. The pre- and aftershock sequence of this major earthquake was recorded by a densified seismological network in Northern Chile and Southern Peru, which provides an excellent data set to study in depth the 3D velocity structure along the subduction megathrust. Based on an automatic event catalogue of nearly 10,000 events spanning the time period March to May 2014 we selected approximately 450 events for a staggered 3D inversion approach. Events are selected to guarantee an even ray coverage through the inversion volume. We only select events with a minimum GAP of 200 to improve depth estimates and therefore increase resolution in the marine forearc. Additionally, we investigate secondary arrivals between the P- and S-wave arrival to improve depth location. Up to now we have processed about 450 events, from which about 150 with at least 30 P- and S-wave observations have been selected for the subsequent 3D tomography. Overall the data quality is very high, which allows arrival time estimates better than 0.05s on average. We will show results from the 1D, 2D, and preliminary 3D inversions and discuss the results together with the obtained seismicity distribution.
S. Ethier; Z. Lin
2003-09-15
Several years of optimization on the super-scalar architecture has made it more difficult to port the current version of the 3D particle-in-cell code GTC to the CRAY/NEC SX-6 vector architecture. This paper explains the initial work that has been done to port this code to the SX-6 computer and to optimize the most time consuming parts. Early performance results are shown and compared to the same test done on the IBM SP Power 3 and Power 4 machines.
Kamon, M.; Phillips, J.R.
1994-12-31
In this paper techniques are presented for preconditioning equations generated by discretizing constrained vector integral equations associated with magnetoquasistatic analysis. Standard preconditioning approaches often fail on these problems. The authors present a specialized preconditioning technique and prove convergence bounds independent of the constraint equations and electromagnetic excitation frequency. Computational results from analyzing several electronic packaging examples are given to demonstrate that the new preconditioning approach can sometimes reduce the number of GMRES iterations by more than an order of magnitude.
Trasischker, Wolfgang; Werkmeister, René M; Zotter, Stefan; Baumann, Bernhard; Torzicky, Teresa; Pircher, Michael; Hitzenberger, Christoph K
2013-11-01
We developed a three-beam Doppler optical coherence tomography (OCT) system that enables measurement of the velocity vector of moving particles in three-dimensions (3-D). The spatial orientation as well as the magnitude of motion can be determined without prior knowledge of the geometry of motion. The system combines three spectral-domain OCT interferometers whose sample beams are focused at the sample by a common focusing lens at three different angles. This provides three spatially independent velocity components simultaneously from which the velocity vector can be reconstructed. We demonstrate the system in a simple test object (rotating disc), a flow phantom, and for blood flow measurements in the retina of a healthy human subject. Measurements of blood flow at a venous bifurcation achieve a good agreement between in- and outflow and demonstrate the reliability of the method. PMID:24247747
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.
Angular velocity estimation from measurement vectors of star tracker.
Liu, Hai-bo; Yang, Jun-cai; Yi, Wen-jun; Wang, Jiong-qi; Yang, Jian-kun; Li, Xiu-jian; Tan, Ji-chun
2012-06-01
In most spacecraft, there is a need to know the craft's angular rate. Approaches with least squares and an adaptive Kalman filter are proposed for estimating the angular rate directly from the star tracker measurements. In these approaches, only knowledge of the vector measurements and sampling interval is required. The designed adaptive Kalman filter can filter out noise without information of the dynamic model and inertia dyadic. To verify the proposed estimation approaches, simulations based on the orbit data of the challenging minisatellite payload (CHAMP) satellite and experimental tests with night-sky observation are performed. Both the simulations and experimental testing results have demonstrated that the proposed approach performs well in terms of accuracy, robustness, and performance. PMID:22695598
Closed formulae to determine the angular velocity of a body-segment based on 3D measurements.
Kocsis, L; Béda, G
2001-01-01
This paper suggests a simple method to determine the global coordinates of the angular velocity and the angular acceleration of a body segment determined by the coordinates of minimum three markers. There are commonly used calculations for the angular quantities basing on the "hypothesis" of planar motion. The usage of approximate methods can result in quantitative and qualitative errors that may completely disort the reality. The method mentioned here is theoretically absolutely correct and can be well used for smoothing noisy data. PMID:11811842
NASA Astrophysics Data System (ADS)
Lokajíček, T.; Kern, H.; Svitek, T.; Ivankina, T.
2014-06-01
Ultrasonic measurements of the 3D velocity distribution of P- and S-waves were performed on a spherical sample of a biotite gneiss from the Outokumpu scientific drill hole. Measurements were done at room temperature and pressures up to 400 and 70 MPa, respectively, in a pressure vessel with oil as a pressure medium. A modified transducer/sample assembly and the installation of a new mechanical system allowed simultaneous measurements of P- and S-wave velocities in 132 independent directions of the sphere on a net in steps of 15°. Proper signals for P- and S-waves could be recorded by coating the sample surface with a high-viscosity shear wave gel and by temporal point contacting of the transmitter and receiver transducers with the sample surface during the measurements. The 3D seismic measurements revealed a strong foliation-related directional dependence (anisotropy) of P- and S-wave velocities, which is confirmed by measurements in a multi-anvil apparatus on a cube-shaped specimen of the same rock. Both experimental approaches show a marked pressure sensitivity of P- and S-wave velocities and velocity anisotropies. With increasing pressure, P- and S-wave velocities increase non-linearly due to progressive closure of micro-cracks. The reverse is true for velocity anisotropy. 3D velocity calculations based on neutron diffraction measurements of crystallographic preferred orientation (CPO) of major minerals show that the intrinsic bulk anisotropy is basically caused by the CPO of biotite constituting about 23 vol.% of the rock. Including the shape of biotite grains and oriented low-aspect ratio microcracks into the modelling increases bulk anisotropy. An important finding from this study is that the measurements on the sample sphere and on the sample cube displayed distinct differences, particularly in shear wave velocities. It is assumed that the differences are due to the different geometries of the samples and the configuration of the transducer-sample assembly
Santosa, Claudia M.; Rose, David D.; Fleming, Neal W.
2015-01-01
Left ventricular end-diastolic pressure (LVEDP) is the foundation of cardiac function assessment. Because of difficulties and risks associated with its direct measurement, correlates of LVEDP derived by pulmonary artery (PA) catheterization or transesophageal echocardiography (TEE) are commonly adopted. TEE has the advantage of being less invasive; however TEE-based estimation of LVEDP using correlates such as left ventricular end-diastolic volume (LVEDV) has technical difficulties that limit its clinical usefulness. Using intraoperative acute normovolemic hemodilution (ANH) as a controlled hemorrhagic model, we examined various mitral flow parameters and three-dimensional reconstructions of left atrial volume as surrogates of LVEDP. Our results demonstrate that peak E wave velocity and left atrial end-diastolic volume (LAEDV) correlated with known changes in intravascular volume associated with ANH. Although left atrial volumetric analysis was done offline in our study, recent advances in echocardiographic software may allow for continuous display and real-time calculation of LAEDV. Along with the ease and reproducibility of acquiring Doppler images of flow across the mitral valve, these two correlates of LVEDP may justify a more widespread use of TEE to optimize intraoperative fluid management. The clinical applicability of peak E wave velocity and LAEDV still needs to be validated during uncontrolled resuscitation. PMID:26236733
NASA Astrophysics Data System (ADS)
Jheng, Y.; Hung, S.; Zhou, Y.; Chang, Y.
2012-12-01
Surface wave travel-time tomography has been widely used as a powerful strategy to image shear wave velocity structure of the Earth's crust and upper mantle, providing comparable information other than body wave tomography. Traditionally, lateral variations of dispersive phase velocities are first obtained at multiple frequencies and then used to invert for shear wave velocity with 1-D depth-dependent sensitivity kernels. However, this approach runs short on considering the directional- and depth-dependence of scattering while surface wave propagating through laterally heterogeneous Earth. To refrain from these shortcomings, we here provide a fully 3-D finite-frequency method based on the Born scattering theory formulated with surface wave mode summation, and apply it to regional fundamental Rayleigh wave travel-time tomography in central Tibet. Our data were collected from Project Hi-CLIMB, which deployed an N-S trending linear array of over 100 broadband seismic stations with a large aperture of 800 km and very dense spacing of ~3-8 km across the Lhasa and Qiangtang terranes during 2004-2005. We follow a standard procedure of ambient noise cross correlation to extract empirical Green's functions of fundamental Rayleigh waves at 10-33 s between station pairs. A multi-taper method is employed to measure the phase differences as a function of period between observed and synthetic Rayleigh waves as well as the corresponding sensitivity kernels for the measured phase delays to 3-D shear wave velocity perturbations in a spherically-symmetric model suitable for central Tibet. A wavelet-based, multi-scale parameterization is invoked in the tomographic inversion to deal with the intrinsically multi-scale nature of unevenly distributed data and resolve the structure with data-adaptive spectral and spatial resolutions. The preliminary result shows that to the north of the Banggong-Nujiang suture (BNS), the crustal shear wave velocity beneath the Qiangtang terrane is
Illés, Tamás
2011-03-01
The EOS system is a new medical imaging device based on low-dose X-rays, gaseous detectors and dedicated software for 3D reconstruction. It was developed by Nobel prizewinner Georges Charpak. A new concept--the vertebral vector--is used to facilitate the interpretation of EOS data, especially in the horizontal plane. We studied 95 cases of idiopathic scoliosis before and after surgery by means of classical methods and using vertebral vectors, in order to compare the accuracy of the two approaches. The vertebral vector permits simultaneous analysis of the scoliotic curvature in the frontal, sagittal and horizontal planes, as precisely as classical methods. The use of the vertebral vector simplifies and facilitates the interpretation of the mass of information provided by EOS. After analyzing the horizontal data, the first goal of corrective intervention would be to reduce the lateral vertebral deviation. The reduction in vertebral rotation seems less important. This is a new element in the therapeutic management of spinal deformations. PMID:22292310
NASA Astrophysics Data System (ADS)
Bremner, P. M.; Panning, M. P.; Russo, R.; Mocanu, V. I.; Stanciu, A. C.; Torpey, M. E.; Hongsresawat, S.; VanDecar, J. C.
2015-12-01
We present new 3-D radially anisotropic and isotropic crustal velocity models beneath central Idaho and eastern Oregon. We produced the velocity models from Love and horizontal component Rayleigh wave group and phase velocity measurements on the IDaho/ORegon (IDOR) Passive seismic network, 86 broadband seismic stations, dataset using ambient noise tomography and the methods of Gallego et. al (2010) and Lin et. al (2008). We calculated inter-station group/phase velocities in narrow frequency bands from travel-time measurements of the rotated stacked horizontal component cross-correlations (bandpass filtered between 2 and 30 seconds), which we used to invert for velocity structure beneath the network. We derived group and phase velocity maps for each frequency band using the damped least-squares inversion method of Tarantola (2005), and then jointly inverted for velocity with depth. Moho depths are prescribed in the joint inversions based on receiver functions, also from the IDOR seismic data, and provides a starting crustal velocity model. Goals of our work include refining models of crustal structure in the accreted Blue Mountain terranes in the western study area; determining the depth extent of the Salmon River Suture/West Idaho Shear Zone (WISZ), which crosses north-south through the middle of the network; determining the architecture of the Idaho batholith, an extensive largely crustal-derived pluton; and examining the nature of the autochthonous (?) North American crust and lithosphere beneath and east of the batholith.
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.
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)
Mack, Tobias; Cierpka, Christian; Kähler, Christian J.
2012-11-01
Astigmatism-PTV is a method that allows to measure with a single camera the fully three-dimensional, three-component velocity field. The technique is ideally suited for microfluidic velocity measurements without errors due to in-plane and out-of-plane averaging (Cierpka et al. Meas Scie Tech 21, 2010). Recently it was shown, that the interface between two fluids or the surrounding fluid and droplets or bubbles can be estimated as well with the technique (Rossi et al., Meas Scie Tech 22, 2010). In this contribution the advantages of both techniques are combined to measure the shape of a droplet inside a micro channel along with the internal 3D flow field of the droplet induced by the surrounding fluid. For the current investigation, particles were only distributed within oil-droplets. Therefore the shape of the droplet could be later reconstructed by the volumetric particle positions and the velocity can be estimated tracking the same particles in consecutive frames of the same dataset. The procedure allows the simultaneous determination of the shape and the droplet velocity as well as the inner flow field and offers a great potential for current research.
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)
Takahashi, Tsutomu; Obana, Koichiro; Yamamoto, Yojiro; Nakanishi, Ayako; Kodaira, Shuichi; Kaneda, Yoshiyuki
2013-05-01
waves at high frequencies (>1 Hz) show collapsed and broadened wave trains caused by multiple scattering in the lithosphere. This study analyzed the envelopes of direct S waves in southwestern Japan and on the western side of the Nankai trough and estimated the spatial distribution of random inhomogeneities by assuming a von Kármán type power spectral density function (PSDF). Strongly inhomogeneous media have been mostly imaged at shallow depth (0-20 km depth) in the onshore area of southwestern Japan, and their PSDF is represented as P(m) ≈ 0.05m-3.7 km3, with m being the spatial wave number, whereas most of the other area shows weak inhomogeneities of which PSDF is P(m) ≈ 0.005m-4.5 km3. At Hyuga-nada in Nankai trough, there is an anomaly of inhomogeneity of which PSDF is estimated as P(m) ≈ 0.01m-4.5 km3. This PSDF has the similar spectral gradient with the weakly inhomogeneous media, but has larger power spectral density than other offshore areas. This anomalous region is broadly located in the subducted Kyushu Palau ridge, which was identified by using velocity structures and bathymetry, and it shows no clear correlation with the fault zones of large earthquakes in past decades. These spatial correlations suggest that possible origins of inhomogeneities at Hyuga-nada are ancient volcanic activity in the oceanic plate or deformed structures due to the subduction of the Kyushu Palau ridge.
34/45-Mbps 3D HDTV digital coding scheme using modified motion compensation with disparity vectors
NASA Astrophysics Data System (ADS)
Naito, Sei; Matsumoto, Shuichi
1998-12-01
This paper describes a digital compression coding scheme for transmitting three dimensional stereo HDTV signals with full resolution at bit-rates around 30 to 40 Mbps to be adapted for PDH networks of the CCITT 3rd digital hierarchy, 34 Mbps and 45 Mbps, SDH networks of 52 Mbps and ATM networks. In order to achieve a satisfactory quality for stereo HDTV pictures, three advanced key technologies are introduced into the MPEG-2 Multi-View Profile, i.e., a modified motion compensation using disparity vectors estimated between the left and right pictures, an adaptive rate control using a common buffer memory for left and right pictures encoding, and a discriminatory bit allocation which results in the improvement of left pictures quality without any degradation of right pictures. From the results of coding experiment conducted to evaluate the coding picture achieved by this coding scheme, it is confirmed that our coding scheme gives satisfactory picture quality even at 34 Mbps including audio and FEC data.
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)
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
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)
Arroyo, I.; Husen, S.; Flueh, E.; Alvarado, G. E.
2008-12-01
The Central Pacific sector of the erosional margin in Costa Rica shows a high seismicity rate, coincident with the subduction of rough-relief ocean floor, and generates earthquakes up to Mw 7. Precise earthquake locations and detailed knowledge of the 3-D velocity structure provide key insights into the dynamics of subduction zones. To this end, we performed a 3-D Local Earthquake Tomography using P-wave traveltimes from 595 selected events recorded by a seismological network of off- and onshore stations, deployed for 6 months in the area. The results reflect the complexity associated to subduction of bathymetric highs and the transition from normal to thickened oceanic crust (Cocos Ridge). The slab is imaged as a high-velocity anomaly with a band of low velocities (LVB) on top enclosing the intraslab events deeper than ~30 km. Below the margin slope, the LVB is locally thickened by at least two seamounts. We observe an abrupt, eastward widening of the LVB, preceded by a low-velocity anomaly under the continental shelf, which we interpret as a big seamount. The thickening coincides with an inverted basin at the inner forearc and a low-velocity anomaly under it. The latter appears in a sector where blocks of inner forearc are uplifted, possibly by underplating of eroded material against the base of the crust. The anomaly promotes seismicity by high-friction with the upper plate, and could be linked to a Mw 6.4 earthquake in 2004. In the west part of the area, the interplate seismicity forms a cluster beneath the continental shelf. Its updip limit coincides with the 150° C isotherm and an increase in Vp along the plate boundary. This further supports a proposed model in which the seismicity onset along the plate interface is mainly due to a decrease in the abundance of the fluids released by subducted sediments. Higher seismicity rates locally concur with seamounts present at the seismogenic zone, while seamounts under the margin slope may shallow the onset of
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)
Lipatov, A. S.; Sittler, E. C.; Hartle, R. E.; Cooper, J. F.; Simpson, D. G.
2011-09-01
In this report we discuss the ion velocity distribution dynamics from the 3D hybrid simulation. In our model the background, pickup, and ionospheric ions are considered as a particles, whereas the electrons are described as a fluid. Inhomogeneous photoionization, electron-impact ionization and charge exchange are included in our model. We also take into account the collisions between the ions and neutrals. The current simulation shows that mass loading by pickup ions H,H2+, CH4+ and N2+ is stronger than in the previous simulations when O + ions are introduced into the background plasma. In our hybrid simulations we use Chamberlain profiles for the atmospheric components. We also include a simple ionosphere model with average mass M = 28 amu ions that were generated inside the ionosphere. The moon is considered as a weakly conducting body. Special attention will be paid to comparing the simulated pickup ion velocity distribution with CAPS T9 observations. Our simulation shows an asymmetry of the ion density distribution and the magnetic field, including the formation of the Alfvén wing-like structures. The simulation also shows that the ring-like velocity distribution for pickup ions relaxes to a Maxwellian core and a shell-like halo.
NASA Technical Reports Server (NTRS)
Lipatov, A. S.; Sittler, E. C., Jr.; Hartle, R. E.; Cooper, J. F.; Simpson, D. G.
2011-01-01
In this report we discuss the ion velocity distribution dynamics from the 3D hybrid simulation. In our model the background, pickup, and ionospheric ions are considered as a particles, whereas the electrons are described as a fluid. Inhomogeneous photoionization, electron-impact ionization and charge exchange are included in our model. We also take into account the collisions between the ions and neutrals. The current simulation shows that mass loading by pickup ions H(+); H2(+), CH4(+) and N2(+) is stronger than in the previous simulations when O+ ions are introduced into the background plasma. In our hybrid simulations we use Chamberlain profiles for the atmospheric components. We also include a simple ionosphere model with average mass M = 28 amu ions that were generated inside the ionosphere. The moon is considered as a weakly conducting body. Special attention will be paid to comparing the simulated pickup ion velocity distribution with CAPS T9 observations. Our simulation shows an asymmetry of the ion density distribution and the magnetic field, including the formation of the Alfve n wing-like structures. The simulation also shows that the ring-like velocity distribution for pickup ions relaxes to a Maxwellian core and a shell-like halo.
NASA Astrophysics Data System (ADS)
Aucejo, M.; Totaro, N.; Guyader, J.-L.
2010-08-01
In noise control, identification of the source velocity field remains a major problem open to investigation. Consequently, methods such as nearfield acoustical holography (NAH), principal source projection, the inverse frequency response function and hybrid NAH have been developed. However, these methods require free field conditions that are often difficult to achieve in practice. This article presents an alternative method known as inverse patch transfer functions, designed to identify source velocities and developed in the framework of the European SILENCE project. This method is based on the definition of a virtual cavity, the double measurement of the pressure and particle velocity fields on the aperture surfaces of this volume, divided into elementary areas called patches and the inversion of impedances matrices, numerically computed from a modal basis obtained by FEM. Theoretically, the method is applicable to sources with complex 3D geometries and measurements can be carried out in a non-anechoic environment even in the presence of other stationary sources outside the virtual cavity. In the present paper, the theoretical background of the iPTF method is described and the results (numerical and experimental) for a source with simple geometry (two baffled pistons driven in antiphase) are presented and discussed.
NASA Astrophysics Data System (ADS)
Leitner, Christoph; Hofmann, Peter; Marschallinger, Robert
2014-05-01
Halite hopper crystals are thought to develop by displacive growth in unconsolidated mud (Gornitz & Schreiber, 1984). The Alpine Haselgebirge, but also e.g. the salt deposits of the Rhine graben (mined at the beginning of the 20th century), comprise hopper crystals with shapes of cuboids, parallelepipeds and rhombohedrons (Görgey, 1912). Obviously, they deformed under oriented stress, which had been tried to reconstruct with respect to the sedimentary layering (Leitner et al., 2013). In the present work, deformed halite hopper crystals embedded in mudrock were automated reconstructed. Object based image analysis (OBIA) has been used successfully in remote sensing for 2D images before. The present study represents the first time that the method was used for reconstruction of three dimensional geological objects. First, manually a reference (gold standard) was created by redrawing contours of the halite crystals on each HRXCT scanning slice. Then, for OBIA, the computer program eCognition was used. For the automated reconstruction a rule set was developed. Thereby, the strength of OBIA was to recognize all objects similar to halite hopper crystals and in particular to eliminate cracks. In a second step, all the objects unsuitable for a structural deformation analysis were dismissed using a support vector machine (SVM) (clusters, polyhalite-coated crystals and spherical halites) The SVM simultaneously drastically reduced the number of halites. From 184 OBIA-objects 67 well shaped remained, which comes close to the number of pre-selected 52 objects. To assess the accuracy of the automated reconstruction, the result before and after SVM was compared to the reference, i.e. the gold standard. State-of the art per-scene statistics were extended to a per-object statistics. Görgey R (1912) Zur Kenntnis der Kalisalzlager von Wittelsheim im Ober-Elsaß. Tschermaks Mineral Petrogr Mitt 31:339-468 Gornitz VM, Schreiber BC (1981) Displacive halite hoppers from the dead sea
NASA Astrophysics Data System (ADS)
Wawerzinek, B.; Ritter, J. R. R.; Roy, C.
2013-08-01
We analyse travel times of shear waves, which were recorded at the MAGNUS network, to determine the 3D shear wave velocity (vS) structure underneath Southern Scandinavia. The travel time residuals are corrected for the known crustal structure of Southern Norway and weighted to account for data quality and pick uncertainties. The resulting residual pattern of subvertically incident waves is very uniform and simple. It shows delayed arrivals underneath Southern Norway compared to fast arrivals underneath the Oslo Graben and the Baltic Shield. The 3D upper mantle vS structure underneath the station network is determined by performing non-linear travel time tomography. As expected from the residual pattern the resulting tomographic model shows a simple and continuous vS perturbation pattern: a negative vS anomaly is visible underneath Southern Norway relative to the Baltic Shield in the east with a contrast of up to 4% vS and a sharp W-E dipping transition zone. Reconstruction tests reveal besides vertical smearing a good lateral reconstruction of the dipping vS transition zone and suggest that a deep-seated anomaly at 330-410 km depth is real and not an inversion artefact. The upper part of the reduced vS anomaly underneath Southern Norway (down to 250 km depth) might be due to an increase in lithospheric thickness from the Caledonian Southern Scandes in the west towards the Proterozoic Baltic Shield in Sweden in the east. The deeper-seated negative vS anomaly (330-410 km depth) could be caused by a temperature anomaly possibly combined with effects due to fluids or hydrous minerals. The determined simple 3D vS structure underneath Southern Scandinavia indicates that mantle processes might influence and contribute to a Neogene uplift of Southern Norway.
NASA Technical Reports Server (NTRS)
Figueroa-Vinas, Adolfo; Gurgiolo, Chris A.; Nieves-Chinchilla, Teresa; Goldstein, Melvyn L.
2010-01-01
It has been suggested by a number of authors that the solar wind electron halo can be formed by the scattering of the strahl. On frequent occasions we have observed in electron angular skymaps (Phi/Theta-plots) of the electron 3D velocity distribution functions) a bursty-filament of particles connecting the strahl to the solar wind core-halo. These are seen over a very limited energy range. When the magnetic field is well off the nominal solar wind flow direction such filaments are inconsistent with any local forces and are probably the result of strong scattering. Furthermore, observations indicates that the strahl component is frequently and significantly anisotropic (Tper/Tpal approx.2). This provides a possible free energy source for the excitation of whistler waves as a possible scattering mechanism. The empirical observational evidence between the halo and the strahl suggests that the strahl population may be, at least in part, the source of the halo component.
NASA Astrophysics Data System (ADS)
Kawamura, T.; Hirata, N.; Sato, H.; Onishi, M.; Noda, K.; Saito, H.
2004-12-01
A deep seismic profiling around the Metropolitan Tokyo region, the Kanto district, started in 2002 under the project titled as the Regional Characterization of the Crust in Metropolitan Areas for Prediction of Strong Ground Motion. The deep seismic profiling, Tokyo Bay 2003, was performed along the major axis of the Tokyo Bay. Because the seismic line in the Miura peninsula passes through a densely populated area, we have a low signal-to-noise ratio data due to the cultural noise. Thus, in addition to the conventional reflection profiling, we deployed 51 off-line recorders with a 3-compornent geophone of 4.5 Hz at carefully selected, quiet receiver points. During 90 days, we had continuous records including many shot signals produced by vibrators on land and air-guns at the bay area. These data provided far-offset first arrival signals and wide angle reflections. We focus on the common receiver gather records of the Tokyo Bay 2003 off-line stations data to identify first arrival and wide angle phases. We applied the first arrival tomography method using a finite difference travel time solver (Hole, 1992) to those data to obtain a 3-D P-wave velocity structure of the uppermost crust along the profile. We obtained a velocity model in and around the Miura peninsula as follows: Across the Tokyo Bay, near surface is a layer with velocities of 2.0-2.5 km/s. A low velocity area corresponds to the fore-arc basin sediment (post Early Miocene) which extends to a depth of approximately 4 km. High velocity patches are located at a depth of approximately 6 km under the Miura peninsula, which we interpreted as Pre-Neogene basement rocks. Finally, the velocity structure obtained by the tomography analysis is used to improve the processing of the reflection profiling data to clarify the deeper structure in the peninsula, including a good velocity constraint for a pre-stack migration of the reflection profiling data.
NASA Astrophysics Data System (ADS)
Voznyuk, I.; Litman, A.; Tortel, H.
2015-08-01
A Quasi-Newton method for reconstructing the constitutive parameters of three-dimensional (3D) penetrable scatterers from scattered field measurements is presented. This method is adapted for handling large-scale electromagnetic problems while keeping the memory requirement and the time flexibility as low as possible. The forward scattering problem is solved by applying the finite-element tearing and interconnecting full-dual-primal (FETI-FDP2) method which shares the same spirit as the domain decomposition methods for finite element methods. The idea is to split the computational domain into smaller non-overlapping sub-domains in order to simultaneously solve local sub-problems. Various strategies are proposed in order to efficiently couple the inversion algorithm with the FETI-FDP2 method: a separation into permanent and non-permanent subdomains is performed, iterative solvers are favorized for resolving the interface problem and a marching-on-in-anything initial guess selection further accelerates the process. The computational burden is also reduced by applying the adjoint state vector methodology. Finally, the inversion algorithm is confronted to measurements extracted from the 3D Fresnel database.
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)
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)
Yin, Y.; Hung, S.
2007-12-01
Seismic tomographic imaging has played a key component to unravel the deep processes that caused the surface morphology and rift magmatism in the southwest United States. Several studies used teleseismic body- wave arrivals recorded by the La Ristra experiment, a dense broadband array of 950-km in length deployed during 1999-2001 and run through the Great Plains, the Rio Grande Rift, and the Colorado Plateau, to construct a 2-D tomographic image of the upper mantle structure beneath this linear array (e.g., Gao et al., 2004). However, because of the inevitable smoothing and damping imposed in the tomographic model, the resulting velocity contrast is too weak to explain distinct P and S waveform changes across the array (Song and Helmberger, 2007). In this study, we include all the data from the La Ristra and available nearby arrays and reexamine finite- frequency travel time delays measured by inter-station cross correlation of waveforms at both high- (0.3-2 Hz for P and 0.1-0.5 Hz for S) and low-frequencies (0.03-0.125 Hz for P and 0.03-0.1 Hz for S). Differing from the previous models that rely on classical ray theory and simple grid parameterization, our inversion considers more realistic 3-D sensitivity kernels for relative travel-time delays and a wavelet-based, multi-scale parameterization that enables to yield robust features with spatially-varying resolutions. Our preliminary P-wave model reveals a prominent low-velocity zone extending from near surface to the depth of 300 km beneath the Rio Grande Rift, while the upper mantle which underlies the Great Plains and the Colorado Plateau is seismically fast. We will demonstrate the difference and improvement of 3-D tomographic models through the use of finite-frequency kernels and multi-scale parameterization.
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.
3D Ion Temperature Reconstruction
NASA Astrophysics Data System (ADS)
Tanabe, Hiroshi; You, Setthivoine; Balandin, Alexander; Inomoto, Michiaki; Ono, Yasushi
2009-11-01
The TS-4 experiment at the University of Tokyo collides two spheromaks to form a single high-beta compact toroid. Magnetic reconnection during the merging process heats and accelerates the plasma in toroidal and poloidal directions. The reconnection region has a complex 3D topology determined by the pitch of the spheromak magnetic fields at the merging plane. A pair of multichord passive spectroscopic diagnostics have been established to measure the ion temperature and velocity in the reconnection volume. One setup measures spectral lines across a poloidal plane, retrieving velocity and temperature from Abel inversion. The other, novel setup records spectral lines across another section of the plasma and reconstructs velocity and temperature from 3D vector and 2D scalar tomography techniques. The magnetic field linking both measurement planes is determined from in situ magnetic probe arrays. The ion temperature is then estimated within the volume between the two measurement planes and at the reconnection region. The measurement is followed over several repeatable discharges to follow the heating and acceleration process during the merging reconnection.
NASA Astrophysics Data System (ADS)
Galve, A.; Charvis, P.; Garcia Cano, L.; Marcaillou, B.
2013-12-01
In 2005, we conducted an onshore-offshore 3D refraction and wide-angle reflection seismic experiment over the rupture zone of the 1958 subduction earthquake that occurred near the border between Colombia and Ecuador. This earthquake was part of a sequence of 3 large ruptures (1942, Mw=7.8; 1958, Mw=7.7; 1979, Mw=8.2), which successively broke from south to north the segments of the megathrust that had been ruptured in 1906 by a single, very large magnitude (8.8) earthquake. Using first arrival traveltime inversion, we constructed a well-defined Vp velocity model of the plate boundary and of the upper and lower plates, down to 25 km depth. The model reveals a 5-km thick, low velocity zone in the upper plate, located immediately above the interplate contact. Because similar low-velocity zones are commonly observed along margins made of oceanic or island-arc accreted terranes, we suggest that the low-velocity zone might result from the alteration and hydration of mafic and ultramafic rocks in the upper plate basement, rather than from hydrofracturing alone. Sediments underplated beneath the inner wedge might contribute to the low-velocity zone but it is unlikely that they are several kilometers thick. Nevertheless, fluids expelled by the compaction and dehydration of those underplated sediments possibly favor the alteration of the overlying rocks. The low-velocity zone is spatially coincident with the 1958 rupture area. Near the toe of the margin, the model shows a low velocity gradient in the outer wedge that we interpret as a zone of highly faulted and fractured rocks or of poorly consolidated sediments. This low velocity/low gradient region forms the oceanward limit of the rupture zones of both the 1958 and the 1979 earthquakes. We suggest that the two earthquake ruptures were arrested by the low velocity zone because its low rigidity contributed to dissipate most of the seismic energy and of the coseismic strain/stress. This might be the reason why the 1958
NASA Astrophysics Data System (ADS)
Pilia, S.; Rawlinson, N.; Direen, N. G.
2013-12-01
Although the notion of Rodinia is quite well accepted in the geoscience community, the location and nature of the eastern continental margin of the Gondwana fragment in Australia is still vague and remains one of the most hotly debated topics in Australian geology. Moreover, most post-Rodinian reconstructions models choose not to tackle the ';Tasmanian challenge', and focus only on the tectonic evolution of mainland southeast Australia, thereby conveniently ignoring the wider tectonic implications of Tasmania's complex geological history. One of the chief limitations of the tectonic reconstructions in this region is a lack of information on Paleozoic (possibly Proterozoic) basement structures. Vast Mesozoic-Cainozoic sedimentary and volcanic cover sequences obscure older outcrops and limit the power of direct observational techniques. In response to these challenges, our effort is focused on ambient seismic noise for imaging 3D crustal shear velocity structure using surface waves, which is capable of illuminating basement structure beneath younger cover. The data used in this study is sourced from the WOMBAT transportable seismic array, which is compounded by around 650 stations spanning the majority of southeastern Australia, including Tasmania and several islands in Bass Strait. To produce the highest quality Green's functions, careful processing of the data has been performed, after which group velocity dispersion measurements have been carried out using a frequency-time analysis method on the symmetric component of the empirical Green's functions (EGFs). Group dispersion measurements from the EGFs have been inverted using a novel hierarchical, transdimensional, Bayesian algorithm to obtain Rayleigh-wave group velocity maps at different periods from 2 to 30 s. The new approach has several advantages in that the number and distribution of model parameters are implicitly controlled by the data, in which the noise is treated as unknown in the inversion. This
NASA Technical Reports Server (NTRS)
Joughin, I.; Kwok, R.; Fahnestock, M.
1996-01-01
Satellite radar interferometry provides an importatn new tool for determining ice-flow velocity. Interferometric measurements made from a single track direction are sensitive only to a single component of the three-dimensional velocity vector. Observations from along thre different track directions would allow the full velocity vector to be determined. A north/south-looking SAR could provide these observations over large portions of the globe but not over large areas of the polar ice sheets. We develop and demonstrate a technique that allows the full three-component velocity vector to be determined from data acquired along two track directions (ascending and descending) under a surface-parallel flow assumption.
Velocity vector analysis of a juncture flow using a three component laser velocimeter
NASA Technical Reports Server (NTRS)
Meyers, J. F.; Hepner, T. E.
1984-01-01
A specialized single-axis, five-beam three-component laser velocimeter was constructed and used to study the flow field in a juncture. The juncture was defined by a blunt leading edged vertical splitter plate and a sharp leading edged horizontal plate. The investigations were conducted in the Low Turbulence Pressure Tunnel at a Mach number of 0.1 and a Reynolds number of 2.2 x 10 to the 6th per meter over the model. The three-component velocity flow field in the juncture was measured, Reynolds stresses calculated, and the velocity vector analysis performed.
Velocity Vector Analysis of a Juncture Flow Using a Three-Component Laser Velocimeter
NASA Technical Reports Server (NTRS)
Meyers, James F.; Hepner, Timothy E.
1984-01-01
A specialized single-axis, five-beam three-component laser velocimeter was constructed and used to study the flow field in a juncture. The juncture was defined by a blunt leading, edged vertical splitter plate and a sharp leading edged horizontal plate. The investigations were conducted in the Low Turbulence Pressure Tunnel at a Mach number of 0.1 and a Reynolds number of 2.2 x 10(exp 6) per meter over the model. The three-component velocity flow field in the juncture was measured, Reynolds stresses calculated and the velocity vector analysis performed.
NASA Astrophysics Data System (ADS)
Tsuji, Y.; Nakajima, J.; Kita, S.; Okada, T.; Matsuzawa, T.; Hasegawa, A.
2007-12-01
Three-dimensional heterogeneous structure beneath northeastern (NE) Japan has been investigated by previous studies and an inclined seismic low-velocity zone is imaged in the mantle wedge sub-parallel to the down-dip direction of the subducting slab (Zhao et al., 1992, Nakajima et al., 2001). However, the heterogeneous structure within the slab has not been well studied even though it is very important to understand the whole process of water transportation from the slab to the surface. Here we show a detailed 3D seismic velocity structure within the subducted Pacific slab around Japan and propose a water-transportation path from the slab to the mantle wedge. In this study, we estimated 3D velocity structure within the Pacific slab by the double-difference tomography (Zhang and Thurber, 2003). We divided the study area, from Hokkaido to Kanto, into 6 areas due to the limitation of memory and computation time. In each area, arrival-time data of 7,500-17,000 events recorded at 70-170 stations were used in the analysis. The total number of absolute travel-time data was about 140,000-312,000 for P wave and 123,000-268,000 for S wave, and differential data were about 736,000-1,920,000 for P wave and 644,000-1,488,000 for S wave. Horizontal and vertical grid separations are 10-25 km and 6.5 km, respectively. RMS residuals of travel times for P wave decreased from 0.23s to 0.09s and for S wave from 0.35s to 0.13s. The obtained results are as follows: (1) a remarkable low-Vs zone exists in the uppermost part of the subducting slab, (2) it extends down to a depth of about 80 km, (3) the termination of this low-Vs zone almost corresponds to the "seismic belt" recently detected in the upper plane of the double seismic zone (Kita et al.,2006; Hasegawa et al., 2007), (4) at depths deeper than 80 km, a low-Vs and high-Vp/Vs zone is apparently distributed in the mantle wedge, immediately above the slab crust. We consider that these features reflect water-transportation processes
NASA Technical Reports Server (NTRS)
Wang, A. H.; Wu, S. T.; Liu, Yang; Hathaway, D.
2008-01-01
We introduce a numerical simulation method for recovering the photospheric velocity field from the vector magnetograms. The traditional method is local correlation tracking (LCT) which is based on measuring the relative displacements of features in blocks of pixels between successive white-light images or magnetograms. Within this method, there are a variety of implementations. One of recently developed implementations is induction local correlation tracking (ILCT) as described by Welsch et al. (2004). They employ the normal component of magnetic induction equation as a constraint to assure consistent solutions. Our numerical method uses the fully three-dimensional MHD equations to recover the photospheric velocity field with individual vector magnetograms. We compare our method to the ILCT method using NOAA AR8210 as an example. The differences and similarities are discussed in detail.
Papadaniil, Chrysa D; Kosmidou, Vasiliki E; Tsolaki, Anthoula; Tsolaki, Magda; Kompatsiaris, Ioannis Yiannis; Hadjileontiadis, Leontios J
2016-10-01
Precise preclinical detection of dementia for effective treatment and stage monitoring is of great importance. Miscellaneous types of biomarkers, e.g., biochemical, genetic, neuroimaging, and physiological, have been proposed to diagnose Alzheimer's disease (AD), the usual suspect behind manifested cognitive decline, and mild cognitive impairment (MCI), a neuropathology prior to AD that does not affect cognitive functions. Event related potential (ERP) methods constitute a non-invasive, inexpensive means of analysis and have been proposed as sensitive biomarkers of cognitive impairment; besides, various ERP components are strongly linked with working memory, attention, sensory processing and motor responses. In this study, an auditory oddball task is employed, to acquire high density electroencephalograhy recordings from healthy elderly controls, MCI and AD patients. The mismatch negativity (MMN) and P300 ERP components are then extracted and their relationship with neurodegeneration is examined. Then, the neural activation at these components is reconstructed using the 3D vector field tomography (3D-VFT) inverse solution. The results reveal a decline of both ERPs amplitude, and a statistically significant prolongation of their latency as cognitive impairment advances. For the MMN, higher brain activation is usually localized in the inferior frontal and superior temporal gyri in the controls. However, in AD, parietal sites exhibit strong activity. Stronger P300 generators are mostly found in the frontal lobe for the controls, but in AD they often shift to the temporal lobe. Reduction in inferior frontal source strength and the switch of the maximum intensity area to parietal and superior temporal sites suggest that these areas, especially the former, are of particular significance when neurodegenerative disorders are investigated. The modulation of MMN and P300 can serve to produce biomarkers of dementia and its progression, and brain imaging can further contribute