The novel high-performance 3-D MT inverse solver
NASA Astrophysics Data System (ADS)
Kruglyakov, Mikhail; Geraskin, Alexey; Kuvshinov, Alexey
2016-04-01
We present novel, robust, scalable, and fast 3-D magnetotelluric (MT) inverse solver. The solver is written in multi-language paradigm to make it as efficient, readable and maintainable as possible. Separation of concerns and single responsibility concepts go through implementation of the solver. As a forward modelling engine a modern scalable solver extrEMe, based on contracting integral equation approach, is used. Iterative gradient-type (quasi-Newton) optimization scheme is invoked to search for (regularized) inverse problem solution, and adjoint source approach is used to calculate efficiently the gradient of the misfit. The inverse solver is able to deal with highly detailed and contrasting models, allows for working (separately or jointly) with any type of MT responses, and supports massive parallelization. Moreover, different parallelization strategies implemented in the code allow optimal usage of available computational resources for a given problem statement. To parameterize an inverse domain the so-called mask parameterization is implemented, which means that one can merge any subset of forward modelling cells in order to account for (usually) irregular distribution of observation sites. We report results of 3-D numerical experiments aimed at analysing the robustness, performance and scalability of the code. In particular, our computational experiments carried out at different platforms ranging from modern laptops to HPC Piz Daint (6th supercomputer in the world) demonstrate practically linear scalability of the code up to thousands of nodes.
3-D Inversion of MT Data for Imaging Deformation Fronts in NW Poland
NASA Astrophysics Data System (ADS)
Ślęzak, Katarzyna; Jóźwiak, Waldemar; Nowożyński, Krzysztof; Brasse, Heinrich
2016-07-01
The Pomerania region (northwest part of Poland) occupies a significant position, where the largest European tectonic boundary is situated. This is the area of the contact between the East European Craton (EEC) and the Paleozoic Platform (PP) and it is known as the Trans-European Suture Zone (TESZ). The TESZ was formed during Paleozoic time as a consequence of the collision of several crustal units and it extends from the Black Sea in the southeast to the British Isles in the northwest. It is a region of key importance for our understanding of the tectonic history of Europe. Previous magnetotelluric (MT) results, based on 2-D inverse modeling, show that the contact zone is of lithospheric discontinuity character and there are distinct differences in geoelectric structures between the Precambrian EEC, transitional zone (TESZ), and the younger PP. The presence of a significant conductor at mid and lower crustal depths was also shown. Thus, the main aim of the research presented here was to obtain detailed, 3-D images of electrical conductivity in the crust and upper mantle and its regional distribution below the TESZ in the northwest part of Poland. To accomplish this task we applied the latest 3-D inversion codes, which allowed us to get more realistic model geometries. Additionally, to confirm and complement the study, the Horizontal Magnetic Tensor (HMT) analysis was realized. This method gives us an opportunity to efficiently locate the position of well-conducting structures. As a result we obtain a clearer, three-dimensional model of conductivity distribution, where highly conductive rock complexes appear which we tentatively connected to deformation fronts.
Compartmentalization of the Coso East Flank geothermal field imaged by 3-D full-tensor MT inversion
NASA Astrophysics Data System (ADS)
Lindsey, Nathaniel J.; Kaven, Joern Ole; Davatzes, Nicholas; Newman, Gregory A.
2017-02-01
Previous magnetotelluric (MT) studies of the high-temperature Coso geothermal system in California identified a subvertical feature of low resistivity (2-5 Ohm m) and appreciable lateral extent (>1 km) in the producing zone of the East Flank field. However, these models could not reproduce gross 3-D effects in the recorded data. We perform 3-D full-tensor inversion and retrieve a resistivity model that out-performs previous 2-D and 3-D off-diagonal models in terms of its fit to the complete 3-D MT data set as well as the degree of modelling bias. Inclusion of secondary Zxx and Zyy data components leads to a robust east-dip (60†) to the previously identified conductive East Flank reservoir feature, which correlates strongly with recently mapped surface faults, downhole well temperatures, 3-D seismic reflection data, and local microseismicity. We perform synthetic forward modelling to test the best-fit dip of this conductor using the response at a nearby MT station. We interpret the dipping conductor as a fractured and fluidized compartment, which is structurally controlled by an unmapped blind East Flank fault zone.
Compartmentalization of the Coso East Flank Geothermal Field Imaged by 3-D Full-tensor MT Inversion
NASA Astrophysics Data System (ADS)
Lindsey, Nathaniel J.; Kaven, Joern Ole; Davatzes, Nicholas; Newman, Gregory A.
2016-11-01
Previous magnetotelluric (MT) studies of the high-temperature Coso geothermal system in California identified a subvertical feature of low resistivity (2 - 5 Ohm-m) and appreciable lateral extent (>1 km) in the producing zone of the East Flank field. However, these models could not reproduce gross 3-D effects in the recorded data. We perform 3-D full-tensor inversion and retrieve a resistivity model that out-performs previous 2-D and 3-D off-diagonal models in terms of its fit to the complete 3-D MT dataset as well as the degree of modeling bias. Inclusion of secondary Zxx and Zyy data components leads to a robust east-dip (60o) to the previously identified conductive East Flank reservoir feature, which correlates strongly with recently mapped surface faults, downhole well temperatures, 3-D seismic reflection data, and local microseismicity. We perform synthetic forward modeling to test the best fit dip of this conductor using the response at a nearby MT station. We interpret the dipping conductor as a fractured and fluidized compartment, which is structurally-controlled by an unmapped blind East Flank fault zone.
NASA Astrophysics Data System (ADS)
Spichak, V. V.
2011-01-01
Possibilities for three-dimensional (3D) magnetotelluric (MT) sounding of local objects contained in the Earth's crust are estimated in a case study of the magma chamber of the Vesuvius volcano. Stochastic inversion of the model MT data by the Markov Chain Monte Carlo (MCMC) method has shown that the most efficient approach is not simultaneous but successive estimation of the geometry and the depth of the anomaly and the assessment of the conductivity distribution within the anomalous region. A zone of equivalence is revealed between the a priori estimate of the depth of the anomalous zone and the a posteriori distribution of electric conductivity within it. Based on the present estimation and previous results, an algorithm for determination of the parameters of local crustal anomaly is proposed.
3D Gravity Inversion using Tikhonov Regularization
NASA Astrophysics Data System (ADS)
Toushmalani, Reza; Saibi, Hakim
2015-08-01
Subsalt exploration for oil and gas is attractive in regions where 3D seismic depth-migration to recover the geometry of a salt base is difficult. Additional information to reduce the ambiguity in seismic images would be beneficial. Gravity data often serve these purposes in the petroleum industry. In this paper, the authors present an algorithm for a gravity inversion based on Tikhonov regularization and an automatically regularized solution process. They examined the 3D Euler deconvolution to extract the best anomaly source depth as a priori information to invert the gravity data and provided a synthetic example. Finally, they applied the gravity inversion to recently obtained gravity data from the Bandar Charak (Hormozgan, Iran) to identify its subsurface density structure. Their model showed the 3D shape of salt dome in this region.
3D Electromagnetic inversion using conjugate gradients
Newman, G.A.; Alumbaugh, D.L.
1997-06-01
In large scale 3D EM inverse problems it may not be possible to directly invert a full least-squares system matrix involving model sensitivity elements. Thus iterative methods must be employed. For the inverse problem, we favor either a linear or non-linear (NL) CG scheme, depending on the application. In a NL CG scheme, the gradient of the objective function is required at each relaxation step along with a univariate line search needed to determine the optimum model update. Solution examples based on both approaches will be presented.
3-D radial gravity gradient inversion
NASA Astrophysics Data System (ADS)
Oliveira, Vanderlei C.; Barbosa, Valéria C. F.
2013-11-01
We have presented a joint inversion of all gravity-gradient tensor components to estimate the shape of an isolated 3-D geological body located in subsurface. The method assumes the knowledge about the depth to the top and density contrast of the source. The geological body is approximated by an interpretation model formed by an ensemble of vertically juxtaposed 3-D right prisms, each one with known thickness and density contrast. All prisms forming the interpretation model have a polygonal horizontal cross-section that approximates a depth slice of the body. Each polygon defining a horizontal cross-section has the same fixed number of vertices, which are equally spaced from 0° to 360° and have their horizontal locations described in polar coordinates referred to an arbitrary origin inside the polygon. Although the number of vertices forming each polygon is known, the horizontal coordinates of these vertices are unknown. To retrieve a set of juxtaposed depth slices of the body, and consequently, its shape, our method estimates the radii of all vertices and the horizontal Cartesian coordinates of all arbitrary origins defining the geometry of all polygons describing the horizontal cross-sections of the prisms forming the interpretation model. To obtain a stable estimate that fits the observed data, we impose constraints on the shape of the estimated body. These constraints are imposed through the well-known zeroth- and first-order Tikhonov regularizations allowing, for example, the estimate of vertical or dipping bodies. If the data do not have enough in-depth resolution, the proposed inverse method can obtain a set of stable estimates fitting the observed data with different maximum depths. To analyse the data resolution and deal with this possible ambiguity, we plot the ℓ2-norm of the residuals (s) against the estimated volume (vp) produced by a set of estimated sources having different maximum depths. If this s × vp curve (s as a function of vp) shows a well
3-D joint inversion of the magnetotelluric phase tensor and vertical magnetic transfer functions
NASA Astrophysics Data System (ADS)
Tietze, Kristina; Ritter, Oliver; Egbert, Gary D.
2015-11-01
With advancing computational resources, 3-D inversion techniques have become feasible in recent years and are now a more widely used tool for magnetotelluric (MT) data interpretation. Galvanic distortion caused by small-scale near-surface inhomogeneities remains an obstacle for 3-D MT inversion which so far has experienced little attention. If not considered properly, the effect on 3-D inversion can be immense and result in erroneous subsurface models and interpretations. To tackle the problem we implemented inversion of the distortion-free phase tensor into the ModEM inversion package. The dimensionless phase tensor components describe only variations of the conductivity structure. When inverting these data, particular care has to be taken of the conductivity structure in the a priori model, which provides the reference frame when transferring the information from phase tensors into absolute conductivity values. Our results obtained with synthetic data show that phase tensor inversion can recover the regional conductivity structure in presence of galvanic distortion if the a priori model provides a reasonable assumption for the regional resistivity average. Joint inversion of phase tensor data and vertical magnetic transfer functions improves recovery of the absolute resistivity structure and is less dependent on the prior model. We also used phase tensor inversion for a data set of more than 250 MT sites from the central San Andreas fault, California, where a number of sites showed significant galvanic distortion. We find the regional structure of the phase tensor inversion results compatible with previously obtained models from impedance inversion. In the vicinity of distorted sites, phase tensor inversion models exhibit more homogeneous/smoother conductivity structures.
Joint inversions of two VTEM surveys using quasi-3D TDEM and 3D magnetic inversion algorithms
NASA Astrophysics Data System (ADS)
Kaminski, Vlad; Di Massa, Domenico; Viezzoli, Andrea
2016-05-01
In the current paper, we present results of a joint quasi-three-dimensional (quasi-3D) inversion of two versatile time domain electromagnetic (VTEM) datasets, as well as a joint 3D inversion of associated aeromagnetic datasets, from two surveys flown six years apart from one another (2007 and 2013) over a volcanogenic massive sulphide gold (VMS-Au) prospect in northern Ontario, Canada. The time domain electromagnetic (TDEM) data were inverted jointly using the spatially constrained inversion (SCI) approach. In order to increase the coherency in the model space, a calibration parameter was added. This was followed by a joint inversion of the total magnetic intensity (TMI) data extracted from the two surveys. The results of the inversions have been studied and matched with the known geology, adding some new valuable information to the ongoing mineral exploration initiative.
NASA Astrophysics Data System (ADS)
Beka, Thomas I.; Smirnov, Maxim; Birkelund, Yngve; Senger, Kim; Bergh, Steffen G.
2016-08-01
Broadband (0.001-1000 s) magnetotelluric (MT) data along a crooked profile collected to investigate the geothermal potential on Spitsbergen could not be fully explained by two-dimensional (2D) models; hence we interpret the data with three-dimensional (3D) inversion herein. To better accommodate 3D features and nearby off profile resistivity structures, the full MT impedance tensor data together with the tipper were inverted. As a model control, a detailed bathymetry is systematically incorporated in the inversion. Our results from testing different inversion settings emphasised that appropriately choosing and tuning the starting model, data error floor and the model regularization together are crucial to obtain optimum benefit from MT field data. Through the 3D inversion, we reproduced out of quadrant impedance components and obtained an overall satisfactory data fit (RMS = 1.05). The final 3D resistivity model displays a complex geology of the near surface region (< 1.5 km), which suggests fractures, localized and regional fault systems and igneous intrusions in the Mesozoic platform cover deposits. The Billefjorden fault zone is revealed as a consistent and deep rooted (> 2 km) conductive anomaly, confirming the regional nature of the fault. The fault zone is positioned between two uplifted basement blocks (> 1000 Ωm) of presumably pre-Devonian (Caledonian) metamorphic rocks, and the fault may have been responsible for deformation in the overlying Paleozoic-Mesozoic unit. Upper crustal conductive anomalies (< 10 Ωm) below the Paleozoic-Mesozoic succession in the western part of the 3D model are interpreted as part of a Devonian basin fill. These conductors are laterally and vertically bounded by resistive rocks, suggesting a conducive environment for deep geothermal heat storage. Having this scenario in an area of a known high heat-flow, deep faults and a thinned lithosphere makes the hypothesis on finding a technologically exploitable geothermal resource
3D Magnetic inversion and remanence: solving the problem
NASA Astrophysics Data System (ADS)
Thomson, V.; Morris, W.
2003-04-01
3D inversion of surface magnetic data is a common processing technique when used in mineral exploration. The major drawback of most 3D inversion algorithms is that they assume that the surface magnetic anomaly is produced by induced magnetization and that there are no remanent magnetization or demagnetization effects present. This has a significant impact when modeling magnetic data that has remanent magnetization. The magnetic anomaly produced by a dipping subsurface body will be identical for a consistent relationship between the dip of the body and the dip of the magnetic vector, regardless of the actual dip of the magnetic body. For example, in the case where a subsurface body is dipping, such as a dipping dike, the dip estimated by the inversion routine will be correct only if induced magnetization is present. This has serious implications for mineral exploration. A solution to the remanence problem is to model the surface magnetic anomaly using a constrained 2D approach rather than 3D. Using a priori information on dip and strike length of a source body, it is possible to approximate the remanence direction and intensity. The 2D solutions can then be rendered into a 3D imaging package to create a model in 3D. A case study was performed on a mafic-ultramafic layered igneous intrusion located in Big Trout Lake, northwestern Ontario, Canada. Large layered igneous intrusions are known to have significant remanence. Like many other layered igneous intrusions such as the Bushveld Complex in South Africa, the Big Trout Lake Complex is highly prospective for Platinum Group Elements (PGEs). Intruded during Archean time, the Big Trout Lake Complex has been subsequently folded and faulted to near vertical. As a consequence of limited surface exposures, knowledge of layering within the pluton and the extent of deformation of the pluton is very limited. Newly acquired high-resolution aeromagnetic data shows a strongly mineralized horizon within the intrusion that
Image Appraisal for 2D and 3D Electromagnetic Inversion
Alumbaugh, D.L.; Newman, G.A.
1999-01-28
Linearized methods are presented for appraising image resolution and parameter accuracy in images generated with two and three dimensional non-linear electromagnetic inversion schemes. When direct matrix inversion is employed, the model resolution and posterior model covariance matrices can be directly calculated. A method to examine how the horizontal and vertical resolution varies spatially within the electromagnetic property image is developed by examining the columns of the model resolution matrix. Plotting the square root of the diagonal of the model covariance matrix yields an estimate of how errors in the inversion process such as data noise and incorrect a priori assumptions about the imaged model map into parameter error. This type of image is shown to be useful in analyzing spatial variations in the image sensitivity to the data. A method is analyzed for statistically estimating the model covariance matrix when the conjugate gradient method is employed rather than a direct inversion technique (for example in 3D inversion). A method for calculating individual columns of the model resolution matrix using the conjugate gradient method is also developed. Examples of the image analysis techniques are provided on 2D and 3D synthetic cross well EM data sets, as well as a field data set collected at the Lost Hills Oil Field in Central California.
Towards Automated Seismic Moment Tensor Inversion in Australia Using 3D Structural Model
NASA Astrophysics Data System (ADS)
Hingee, M.; Tkalcic, H.; Fichtner, A.; Sambridge, M.; Kennett, B. L.; Gorbatov, A.
2009-12-01
There is significant seismic activity in the region around Australia, largely due to the plate boundaries to the north and to the east of the mainland. This seismicity poses serious seismic and tsunamigenic hazard in a wider region, and risk to coastal areas of Australia, and is monitored by Geoscience Australia (GA) using a network of permanent broadband seismometers within Australia. Earthquake and tsunami warning systems were established by the Australian Government and have been using the waveforms from the GA seismological network. The permanent instruments are augmented by non-GA seismic stations based both within and outside of Australia. In particular, seismic moment tensor (MT) solutions for events around Australia as well as local distances are useful for both warning systems and geophysical studies in general. These monitoring systems, however, currently use only one dimensional, spherically-symmetric models of the Earth for source parameter determination. Recently, a novel 3D model of Australia and the surrounding area has been developed from spectral element simulations [1], taking into account not only velocity heterogeneities, but also radial anisotropy and seismic attenuation. This development, inter alia, introduces the potential of providing significant improvements in MT solution accuracy. Allowing reliable MT solutions with reduced dependence on non-GA stations is a secondary advantage. We studied the feasibility of using 1D versus 3D structural models. The accuracy of the 3D model has been investigated, confirming that these models are in most cases superior to the 1D models. A full MT inversion method using a point source approximation was developed as the first step, keeping in mind that for more complex source time functions, a finite source inversion will be needed. Synthetic experiments have been performed with random noise added to the signal to test the code in the both 1D and 3D setting, using a precomputed library of structural Greens
3D-spectral CDIs: a fast alternative to 3D inversion?
NASA Astrophysics Data System (ADS)
Macnae, James
2015-09-01
Virtually all airborne electromagnetic (AEM) data is interpreted using stitched 1D conductivity sections, derived from constrained inversion or fast but fairly accurate approximations. A small subset of this AEM data recently has been inverted using either block 3D models or thin plates, which processes have limitations in terms of cost and accuracy, and the results are in general strongly biased by the choice of starting models. Recent developments in spectral modelling have allowed fast 3D approximations of the EM response of both vortex induction and current gathering for simple geological target geometries. Fitting these spectral responses to AEM data should be sufficient to accurately locate current systems within the ground, and the behaviour of these local current systems can in theory approximately define a conductivity structure in 3D. This paper describes the results of initial testing of the algorithm in fitting vortex induction in a small target at the Forrestania test range, Western Australia, using results from a versatile time-domain electromagnetic (VTEM)-Max survey.
Computational and methodological developments towards 3D full waveform inversion
NASA Astrophysics Data System (ADS)
Etienne, V.; Virieux, J.; Hu, G.; Jia, Y.; Operto, S.
2010-12-01
Full waveform inversion (FWI) is one of the most promising techniques for seismic imaging. It relies on a formalism taking into account every piece of information contained in the seismic data as opposed to more classical techniques such as travel time tomography. As a result, FWI is a high resolution imaging process able to reach a spatial accuracy equal to half a wavelength. FWI is based on a local optimization scheme and therefore the main limitation concerns the starting model which has to be closed enough to the real one in order to converge to the global minimum. Another counterpart of FWI is the required computational resources when considering models and frequencies of interest. The task becomes even more tremendous when one tends to perform the inversion using the elastic equation instead of using the acoustic approximation. This is the reason why until recently most studies were limited to 2D cases. In the last few years, due to the increase of the available computational power, FWI has focused a lot of interests and continuous efforts towards inversion of 3D models, leading to remarkable applications up to the continental scale. We investigate the computational burden induced by FWI in 3D elastic media and propose some strategic features leading to the reduction of the numerical cost while providing a great flexibility in the inversion parametrization. First, in order to release the memory requirements, we developed our FWI algorithm in the frequency domain and take benefit of the wave-number redundancy in the seismic data to process a quite reduced number of frequencies. To do so, we extract frequency solutions from time marching techniques which are efficient for 3D structures. Moreover, this frequency approach permits a multi-resolution strategy by proceeding from low to high frequencies: the final model at one frequency is used as the starting model for the next frequency. This procedure overcomes partially the non-linear behavior of the inversion
Inverse Tomo-Lithography for Making Microscopic 3D Parts
NASA Technical Reports Server (NTRS)
White, Victor; Wiberg, Dean
2003-01-01
According to a proposal, basic x-ray lithography would be extended to incorporate a technique, called inverse tomography, that would enable the fabrication of microscopic three-dimensional (3D) objects. The proposed inverse tomo-lithographic process would make it possible to produce complex shaped, submillimeter-sized parts that would be difficult or impossible to make in any other way. Examples of such shapes or parts include tapered helices, paraboloids with axes of different lengths, and even Archimedean screws that could serve as rotors in microturbines. The proposed inverse tomo-lithographic process would be based partly on a prior microfabrication process known by the German acronym LIGA (lithographie, galvanoformung, abformung, which means lithography, electroforming, molding). In LIGA, one generates a precise, high-aspect ratio pattern by exposing a thick, x-ray-sensitive resist material to an x-ray beam through a mask that contains the pattern. One can electrodeposit metal into the developed resist pattern to form a precise metal part, then dissolve the resist to free the metal. Aspect ratios of 100:1 and patterns into resist thicknesses of several millimeters are possible.
3-D inversion of synthetic marine magnetotelluric data: resolution and sensitivity
NASA Astrophysics Data System (ADS)
Tada, N.; Baba, K.; Siripunvaraporn, W.; Uyeshima, M.; Utada, H.
2010-12-01
In recent years, seafloor magnetotelluric (MT) observation is carried out by using an increasing number of ocean bottom electromagnetometers (OBEMs) not only along a line but also in 2-D array. Thus, imaging electrical conductivity structures under the seafloor in 3-D is now feasible. A 3-D approach is indispensable especially for marine MT data, because the electric and magnetic fields observed at the seafloor are heavily distorted by the rugged seafloor topography and the distribution of land and sea which are generally 3-D. It is very important to incorporate the topography in a 3-D model for an accurate estimation of the conductivity structure beneath seafloor that is generally more resistive than seawater by several orders of magnitude. WSINV3DMT (Siripunvaraporn et al., 2005) is one of 3-D inversion codes that are now of practical use, but the original WSINV3DMT is not applicable to marine MT data because of two reasons. 1) MT responses are calculated only at the boundary corresponding to the Earth surface. 2) We have to use fine mesh design because an observation site must locate exactly at the center of the top surface of a block, which needs large memory that even a highest performance computer can not handle. We propose an extended version of the WSINV3DMT by solving the two problems shown above so that it can be applied to the marine MT data. The extended version of the WSINV3DMT is tested using synthetic models including a 3-D anomaly, seawater and topographic variation. Here shown is an example of a checkerboard test by using a model in which 10 ohm-m and 100 ohm-m blocks are put alternately in both horizontal and vertical directions. The model is composed of 5 blocks in horizontal directions and of 4 blocks in vertical direction with a background of a 31.6 ohm-m half-space below actual topography. The calculation area in the inversion is 7440 × 7440 × 1008 km, and is discretized at 35 blocks in the x and y directions, and 69 blocks in the z
Representation of 3-D surface orientation by velocity and disparity gradient cues in area MT.
Sanada, Takahisa M; Nguyenkim, Jerry D; Deangelis, Gregory C
2012-04-01
Neural coding of the three-dimensional (3-D) orientation of planar surface patches may be an important intermediate step in constructing representations of complex 3-D surface structure. Spatial gradients of binocular disparity, image velocity, and texture provide potent cues to the 3-D orientation (tilt and slant) of planar surfaces. Previous studies have described neurons in both dorsal and ventral stream areas that are selective for surface tilt based on one or more of these gradient cues. However, relatively little is known about whether single neurons provide consistent information about surface orientation from multiple gradient cues. Moreover, it is unclear how neural responses to combinations of surface orientation cues are related to responses to the individual cues. We measured responses of middle temporal (MT) neurons to random dot stimuli that simulated planar surfaces at a variety of tilts and slants. Four cue conditions were tested: disparity, velocity, and texture gradients alone, as well as all three gradient cues combined. Many neurons showed robust tuning for surface tilt based on disparity and velocity gradients, with relatively little selectivity for texture gradients. Some neurons showed consistent tilt preferences for disparity and velocity cues, whereas others showed large discrepancies. Responses to the combined stimulus were generally well described as a weighted linear sum of responses to the individual cues, even when disparity and velocity preferences were discrepant. These findings suggest that area MT contains a rudimentary representation of 3-D surface orientation based on multiple cues, with single neurons implementing a simple cue integration rule.
NASA Astrophysics Data System (ADS)
Ullmann, A.; Scheunert, M.; Afanasjew, M.; Börner, R.-U.; Siemon, B.; Spitzer, K.
2016-07-01
As a standard procedure, multi-frequency helicopter-borne electromagnetic (HEM) data are inverted to conductivity-depth models using 1-D inversion methods, which may, however, fail in areas of strong lateral conductivity contrasts (so-called induction anomalies). Such areas require more realistic multi-dimensional modelling. Since the full 3-D inversion of an entire HEM data set is still extremely time consuming, our idea is to combine fast 1-D and accurate but numerically expensive 3-D inversion of HEM data in such a way that the full 3-D inversion is only carried out for those parts of a HEM survey which are affected by induction anomalies. For all other parts, a 1-D inversion method is sufficient. We present a newly developed algorithm for identification, selection, and extraction of induction anomalies in HEM data sets and show how the 3-D inversion model of the anomalous area is re-integrated into the quasi-1-D background. Our proposed method is demonstrated to work properly on a synthetic and a field HEM data set from the Cuxhaven tunnel valley in Germany. We show that our 1-D/3-D approach yields better results compared to 1-D inversions in areas where 3-D effects occur.
3D inversion based on multi-grid approach of magnetotelluric data from Northern Scandinavia
NASA Astrophysics Data System (ADS)
Cherevatova, M.; Smirnov, M.; Korja, T. J.; Egbert, G. D.
2012-12-01
In this work we investigate the geoelectrical structure of the cratonic margin of Fennoscandian Shield by means of magnetotelluric (MT) measurements carried out in Northern Norway and Sweden during summer 2011-2012. The project Magnetotellurics in the Scandes (MaSca) focuses on the investigation of the crust, upper mantle and lithospheric structure in a transition zone from a stable Precambrian cratonic interior to a passive continental margin beneath the Caledonian Orogen and the Scandes Mountains in western Fennoscandia. Recent MT profiles in the central and southern Scandes indicated a large contrast in resistivity between Caledonides and Precambrian basement. The alum shales as a highly conductive layers between the resistive Precambrian basement and the overlying Caledonian nappes are revealed from this profiles. Additional measurements in the Northern Scandes were required. All together data from 60 synchronous long period (LMT) and about 200 broad band (BMT) sites were acquired. The array stretches from Lofoten and Bodo (Norway) in the west to Kiruna and Skeleftea (Sweden) in the east covering an area of 500x500 square kilometers. LMT sites were occupied for about two months, while most of the BMT sites were measured during one day. We have used new multi-grid approach for 3D electromagnetic (EM) inversion and modelling. Our approach is based on the OcTree discretization where the spatial domain is represented by rectangular cells, each of which might be subdivided (recursively) into eight sub-cells. In this simplified implementation the grid is refined only in the horizontal direction, uniformly in each vertical layer. Using multi-grid we manage to have a high grid resolution near the surface (for instance, to tackle with galvanic distortions) and lower resolution at greater depth as the EM fields decay in the Earth according to the diffusion equation. We also have a benefit in computational costs as number of unknowns decrease. The multi-grid forward
3D Strucutural Geological Model of the Alpi Mt. Area (Southern Italy)
NASA Astrophysics Data System (ADS)
La Bruna, Vincenzo; Lamarche, Juliette; Viseur, Sophie; Agosta, Fabrizio; Prosser, Giacomo
2016-04-01
The study area is located in the inner portion of the southern Apennines fold-and-thrust belt. The Alpi Mt. is the only portion of the Apulian domain cropping in this sector. In fact, it is considered as a structural analogue of the Val d'Agri and Tempa Rossa reservoirs (Basilicata). The Alpi Mt. tectonic unit is composed of two main cronostratigraphic intervals, represented by a 2000m-thick Mesozoic carbonate succession and a Messinian mixed carbonate-terrigenous succession. The Messinian interval is made up of a Lower Messinian sedimentary cycle, wich form a paraconformity with the underlying Mesozoic carbonates, and an Upper Messinian cycle characterized by a marked unconformity at the bottom. This study aims to better understand the role exerted by the precontractional tectonic structures during the Messinian interval, wich are responsible for the development of the sedimentary angular unconformity. To reach this goal, a 3D structural geological model was build up by using the Gocad(R) software. The construction of the 3D model was gained through the integration of several results related to geological field mapping, well log analysis and seismic reflection data. Focusing on the Upper Messinian sedimentary horizon, in order to achieve the true geometry and kinematics of the high-angle extensional faults that bound the sedimentary depocenters, the model was restored through vertical line methodology. This process allows to obtain more information about location, geometry, and sedimentary depocenter orientations. Furthermore, the 3D structural model brings some important results from the 3D fault analysis that are represented by attitude, geometry and dimensional parameters of the fault network that affect the study area.
Flow-Through Stream Modeling with MODFLOW and MT3D: Certainties and Limitations.
Ben Simon, Rose; Bernard, Stéphane; Meurville, Charles; Rebour, Vincent
2015-01-01
This paper aims to assess MODFLOW and MT3D capabilities for simulating the spread of contaminants from a river exhibiting an unusual relationship with an alluvial aquifer, with the groundwater head higher than the river head on one side and lower on the other (flow-through stream). A series of simulation tests is conducted using a simple hypothetical model so as to characterize and quantify these limitations. Simulation results show that the expected contaminant spread could be achieved with a specific configuration composed of two sets of parameters: (1) modeled object parameters (hydraulic groundwater gradient, hydraulic conductivity values of aquifer and streambed), and (2) modeling parameters (vertical discretization of aquifer, horizontal refinement of stream modeled with River [RIV] package). The influence of these various parameters on simulation results is investigated, and potential complications and errors are identified. Contaminant spread from stream to aquifer is not always reproduced by MT3D due to the RIV package's inability to simulate lateral exchange fluxes between stream and aquifer. This paper identifies the need for a MODFLOW streamflow package allowing lateral stream-aquifer interactions and streamflow routine calculations. Such developments could be of particular interest for modeling contaminated flow-through streams.
3D stochastic joint inversion of gravity and magnetic data
NASA Astrophysics Data System (ADS)
Shamsipour, Pejman; Marcotte, Denis; Chouteau, Michel
2012-04-01
A novel stochastic joint inversion method based on cokriging is applied to estimate density and magnetic susceptibility distributions from gravity and total magnetic field data. The method fully integrates the physical relations between density-gravity, on one hand, and magnetic susceptibility-total magnetic field, on the other hand. As a consequence, when the data are considered noise-free, the responses from the inverted density and susceptibility data exactly reproduce the observed data. The required density and magnetic susceptibility auto- and cross covariance are assumed to follow a linear model of coregionalization (LCM). The parameters of the LCM are estimated from v-v plot fitting of the gravity and total magnetic experimental covariances. The model is tested on two synthetic cases and one real data set, the Perseverance mine (Quebec, Canada). Joint inversions are compared to separate inversions. The joint inversions better recover the known models in the synthetic cases. With the real data set, better definition and location of the mineralized lenses are achieved by joint inversion.
NASA Astrophysics Data System (ADS)
Miller, Craig A.; Williams-Jones, Glyn
2016-06-01
A new 3D geophysical model of the Mt Tongariro Volcanic Massif (TgVM), New Zealand, provides a high resolution view of the volcano's internal structure and hydrothermal system, from which we derive implications for volcanic hazards. Geologically constrained 3D inversions of potential field data provides a greater level of insight into the volcanic structure than is possible from unconstrained models. A complex region of gravity highs and lows (± 6 mGal) is set within a broader, ~ 20 mGal gravity low. A magnetic high (1300 nT) is associated with Mt Ngauruhoe, while a substantial, thick, demagnetised area occurs to the north, coincident with a gravity low and interpreted as representing the hydrothermal system. The hydrothermal system is constrained to the west by major faults, interpreted as an impermeable barrier to fluid migration and extends to basement depth. These faults are considered low probability areas for future eruption sites, as there is little to indicate they have acted as magmatic pathways. Where the hydrothermal system coincides with steep topographic slopes, an increased likelihood of landslides is present and the newly delineated hydrothermal system maps the area most likely to have phreatic eruptions. Such eruptions, while small on a global scale, are important hazards at the TgVM as it is a popular hiking area with hundreds of visitors per day in close proximity to eruption sites. The model shows that the volume of volcanic material erupted over the lifespan of the TgVM is five to six times greater than previous estimates, suggesting a higher rate of magma supply, in line with global rates of andesite production. We suggest that our model of physical property distribution can be used to provide constraints for other models of dynamic geophysical processes occurring at the TgVM.
NASA Astrophysics Data System (ADS)
Campanyà, Joan; Ogaya, Xènia; Jones, Alan G.; Rath, Volker; Vozar, Jan; Meqbel, Naser
2016-12-01
As a consequence of measuring time variations of the electric and the magnetic field, which are related to current flow and charge distribution, magnetotelluric (MT) data in 2-D and 3-D environments are not only sensitive to the geoelectrical structures below the measuring points but also to any lateral anomalies surrounding the acquisition site. This behaviour complicates the characterization of the electrical resistivity distribution of the subsurface, particularly in complex areas. In this manuscript we assess the main advantages of complementing the standard MT impedance tensor (Z) data with interstation horizontal magnetic tensor (H) and geomagnetic transfer function (T) data in constraining the subsurface in a 3-D environment beneath a MT profile. Our analysis was performed using synthetic responses with added normally distributed and scattered random noise. The sensitivity of each type of data to different resistivity anomalies was evaluated, showing that the degree to which each site and each period is affected by the same anomaly depends on the type of data. A dimensionality analysis, using Z, H and T data, identified the presence of the 3-D anomalies close to the profile, suggesting a 3-D approach for recovering the electrical resistivity values of the subsurface. Finally, the capacity for recovering the geoelectrical structures of the subsurface was evaluated by performing joint inversion using different data combinations, quantifying the differences between the true synthetic model and the models from inversion process. Four main improvements were observed when performing joint inversion of Z, H and T data: (1) superior precision and accuracy at characterizing the electrical resistivity values of the anomalies below and outside the profile; (2) the potential to recover high electrical resistivity anomalies that are poorly recovered using Z data alone; (3) improvement in the characterization of the bottom and lateral boundaries of the anomalies with low
NASA Astrophysics Data System (ADS)
Shi, Yan-Ling; Hu, Zu-Zhi; Huang, Wen-Hui; Wei, Qiang; Zhang, Sheng; Meng, Cui-Xian; Ji, Lian-Sheng
2016-09-01
The coal-bearing strata of the deep Upper Paleozoic in the GS Sag have high hydrocarbon potential. Because of the absence of seismic data, we use electromagnetic (MT) and gravity data jointly to delineate the distribution of deep targets based on well logging and geological data. First, a preliminary geological model is established by using three-dimensional (3D) MT inversion results. Second, using the formation density and gravity anomalies, the preliminary geological model is modified by interactive inversion of the gravity data. Then, we conduct MT-constrained inversion based on the modified model to obtain an optimal geological model until the deviations at all stations are minimized. Finally, the geological model and a seismic profile in the middle of the sag is analysed. We determine that the deep reflections of the seismic profile correspond to the Upper Paleozoic that reaches thickness up to 800 m. The processing of field data suggests that the joint MT-gravity modeling and constrained inversion can reduce the multiple solutions for single geophysical data and thus improve the recognition of deep formations. The MT-constrained inversion is consistent with the geological features in the seismic section. This suggests that the joint MT and gravity modeling and constrained inversion can be used to delineate deep targets in similar basins.
Direct inversion of digital 3D Fraunhofer holography maps.
Podorov, Sergei G; Förster, Eckhart
2016-01-20
Differential Fourier holography (DFH) gives an exact mathematical solution of the inverse problem of diffraction in the Fraunhofer regime. After the first publication [Opt. Express15, 9954 (2007)], DFH was successfully applied in many experiments to obtain amplitude and phase information about two-dimensional images. In this paper, we demonstrate numerically the possibility to apply DFH also for investigation of unknown three-dimensional objects. The first simulation is made for a double-spiral structure plus a line as a reference object.
3D parallel inversion of time-domain airborne EM data
NASA Astrophysics Data System (ADS)
Liu, Yun-He; Yin, Chang-Chun; Ren, Xiu-Yan; Qiu, Chang-Kai
2016-12-01
To improve the inversion accuracy of time-domain airborne electromagnetic data, we propose a parallel 3D inversion algorithm for airborne EM data based on the direct Gauss-Newton optimization. Forward modeling is performed in the frequency domain based on the scattered secondary electrical field. Then, the inverse Fourier transform and convolution of the transmitting waveform are used to calculate the EM responses and the sensitivity matrix in the time domain for arbitrary transmitting waves. To optimize the computational time and memory requirements, we use the EM "footprint" concept to reduce the model size and obtain the sparse sensitivity matrix. To improve the 3D inversion, we use the OpenMP library and parallel computing. We test the proposed 3D parallel inversion code using two synthetic datasets and a field dataset. The time-domain airborne EM inversion results suggest that the proposed algorithm is effective, efficient, and practical.
NASA Astrophysics Data System (ADS)
Portal, Angélie; Fargier, Yannick; Lénat, Jean-François; Labazuy, Philippe
2016-04-01
The electrical resistivity tomography (ERT) method, initially developed for environmental and engineering exploration, is now commonly used for geological structures imaging. Such structures can present complex characteristics that conventional 2D inversion processes cannot perfectly integrate. Here we present a new 3D inversion algorithm named EResI, firstly developed for levee investigation, and presently applied to the study of a complex lava dome (the Puy de Dôme volcano, France). EResI algorithm is based on a conventional regularized Gauss-Newton inversion scheme and a 3D non-structured discretization of the model (double grid method based on tetrahedrons). This discretization allows to accurately model the topography of investigated structure (without a mesh deformation procedure) and also permits a precise location of the electrodes. Moreover, we demonstrate that a complete 3D unstructured discretization limits the number of inversion cells and is better adapted to the resolution capacity of tomography than a structured discretization. This study shows that a 3D inversion with a non-structured parametrization has some advantages compared to classical 2D inversions. The first advantage comes from the fact that a 2D inversion leads to artefacts due to 3D effects (3D topography, 3D internal resistivity). The second advantage comes from the fact that the capacity to experimentally align electrodes along an axis (for 2D surveys) depends on the constrains on the field (topography...). In this case, a 2D assumption induced by 2.5D inversion software prevents its capacity to model electrodes outside this axis leading to artefacts in the inversion result. The last limitation comes from the use of mesh deformation techniques used to accurately model the topography in 2D softwares. This technique used for structured discretization (Res2dinv) is prohibed for strong topography (>60 %) and leads to a small computational errors. A wide geophysical survey was carried out
Rodriguez, Brian D.; Sweetkind, Donald S.
2015-01-01
The 3-D inversion was generally able to reproduce the gross resistivity structure of the “known” model, but the simulated conductive volcanic composite unit horizons were often too shallow when compared to the “known” model. Additionally, the chosen computation parameters such as station spacing appear to have resulted in computational artifacts that are difficult to interpret but could potentially be removed with further refinements of the 3-D resistivity inversion modeling technique.
A Preprocessor for Modeling Nonpoint Sources in Fractured Media using MODFLOW and MT3D
NASA Astrophysics Data System (ADS)
Mun, Y.; Uchrin, C. G.
2002-05-01
There are a multitude of fractures in the geological structure of fractured media which act as conduits for subsurface fluid flow. The hydraulic properties of this flow are very heterogeneous even within a single unit and this heterogeneity is very localized. As a result, modeling flow in fractured media is difficult due to this heterogeneity. There are two major approaches to simulate the flow and transport of fluid flow in fractured media: the discrete fracture approach and the continuum approach. Precise characteristics such as geometry are required to use the discrete fracture approach. It, however, is difficult to determine the fluid flow through the fractures because of inaccessibility. In the continuum approach, although head distributions can match to well data, chemical concentration distributions are hard to match well sample concentration observations, because some aquifers are dominated by advective transport and others are likely to serve as reservoirs for immobile solutes. The MODFLOW preprocessor described in this paper has been developed and applied to the Cranberry Lake system in Northwestern New Jersey. Cranberry Lake has exhibited eutrophic characteristics for some time by nonpoint sources including surface water runoff, leaching from local septic systems and direct deposition. It has been estimated that 70% of the nutrient loading to the lake flows through fractured media from septic systems. The preprocessor presented in this paper utilizes percolation theory, which is concerned with the existence of ropen paths_. The percolation threshold of a body-centered cubic lattice (3D), a square lattice (2D) and several other percolation numbers are applied to make the model system represent the fractured media. The distribution of hydraulic head within groundwater is simulated by MODFLOW and the advection-dispersion equation of nitrate transport is solved by MT3D. This study also simulates boron transport as an indicator.
NASA Astrophysics Data System (ADS)
Fadel, I.; van der Meijde, M.; Kerle, N.
2013-12-01
Non-uniqueness of satellite gravity interpretation has been usually reduced by using a priori information from various sources, e.g. seismic tomography models. The reduction in non-uniqueness has been based on velocity-density conversion formulas or user interpretation for 3D subsurface structures (objects) in seismic tomography models. However, these processes introduce additional uncertainty through the conversion relations due to the dependency on the other physical parameters such as temperature and pressure, or through the bias in the interpretation due to user choices and experience. In this research, a new methodology is introduced to extract the 3D subsurface structures from 3D geophysical data using a state-of-art 3D Object Oriented Image Analysis (OOA) technique. 3D OOA is tested using a set of synthetic models that simulate the real situation in the study area of this research. Then, 3D OOA is used to extract 3D subsurface objects from a real 3D seismic tomography model. The extracted 3D objects are used to reconstruct a forward model and its response is compared with the measured satellite gravity. Finally, the result of the forward modelling, based on the extracted 3D objects, is used to constrain the inversion process of satellite gravity data. Through this work, a new object-based approach is introduced to interpret and extract the 3D subsurface objects from 3D geophysical data. This can be used to constrain modelling and inversion of potential field data using the extracted 3D subsurface structures from other methods. In summary, a new approach is introduced to constrain inversion of satellite gravity measurements and enhance interpretation capabilities.
NASA Astrophysics Data System (ADS)
Rezaie, Mohammad; Moradzadeh, Ali; Kalate, Ali Nejati; Aghajani, Hamid
2017-01-01
Inversion of gravity data is one of the important steps in the interpretation of practical data. One of the most interesting geological frameworks for gravity data inversion is the detection of sharp boundaries between orebody and host rocks. The focusing inversion is able to reconstruct a sharp image of the geological target. This technique can be efficiently applied for the quantitative interpretation of gravity data. In this study, a new reweighted regularized method for the 3D focusing inversion technique based on Lanczos bidiagonalization method is developed. The inversion results of synthetic data show that the new method is faster than common reweighted regularized conjugate gradient method to produce an acceptable solution for focusing inverse problem. The new developed inversion scheme is also applied for inversion of the gravity data collected over the San Nicolas Cu-Zn orebody in Zacatecas State, Mexico. The inversion results indicate a remarkable correlation with the true structure of the orebody that is achieved from drilling data.
Regional conductivity structure of Cascadia from 3D inversion of USArray magnetotelluric data
NASA Astrophysics Data System (ADS)
Egbert, G. D.; Patro, P. K.
2008-12-01
Magnetotelluric (MT) data are being acquired in a series of temporary arrays deployed across the continental US through EMScope, a part of the USArray component of EarthScope. Initial deployments in 2006 and 2007 acquired data at 110 sites covering the US Pacific Northwest. The MT sites, distributed with the same nominal spacing as the USArray seismic transportable array (~75 km), produced data in the period range 10- 20,000s of very good to excellent quality. The most striking and robust feature revealed by 3D inversion of this dataset is an extensive lower crustal conductor covering most of the study area southeast of a line running from the California border at the coast to the Blue Mountains of Northeastern Oregon. The conductance of this layer, which is about 15 km thick with a top at roughly 20 km depth, exceeds 3000 S beneath the he Northwest Basin and Range (BR) province of southeastern Oregon. The high conductivity in this region is inferred to result from fluids - including possibly partial melt at depth - associated with magmatic underplating and BR extension. The lower crust is much more resistive beneath the Coast Range, Willamette Valley and Puget Lowlands of Western Washington and Oregon, and beneath the Columbia Plateau. This area of resistive crust, which was derived from a large fragment of thickened oceanic lithosphere that was accreted to North America at approximately 48 Ma ("Siletzia"), is revealed by geological and geodetic studies to be strong, accommodating tectonic stresses through rigid block rotations. In contrast, the area to the southeast characterized by high conductivity in the lower crust is actively deforming, consistent with an important role for fluids in weakening of continental crust. The resistive Siletzia crust is broken by an elongated N-S zone of high conductivity beneath the Cascade volcanoes. High conductivities beneath the volcanoes also most likely reflect the presence of interconnected fluids, in this case released
A Spatially Periodic Solute Boundary for MT3DMS and PHT3D.
Laattoe, Tariq; Post, Vincent E A; Werner, Adrian D
2016-12-14
The assumption of spatial repetition is commonly made when producing bedform scale models of the hyporheic zone. Two popular solute transport codes, MT3DMS and PHT3D, do not currently provide the necessary boundary condition required to simulate spatial periodicity in hyporheic zone transport problems. In this study, we develop a spatially periodic boundary (SPB) for solutes that is compatible with a SPB that was previously developed for MODFLOW to simulate the flow component of spatially periodic problems. The approach is ideal for simulating groundwater flow and transport patterns under repeating surface features, such as ripples or dunes on the bottom of a lake or stream. The appropriate block-centered finite-difference approach to implement the boundary is presented and the necessary source code modifications are discussed. The performance of the solute SPB, operating in conjunction with the groundwater flow SPB, is explored through comparison of a multi-bedform hyporheic-zone model with a single bedform variant. The new boundary conditions perform well in situations where both dispersive effects and lateral seepage flux in the underflow regime beneath the hyporheic zone are minimal.
3D Propagation and Geoacoustic Inversion Studies in the Mid-Atlantic Bight
2016-06-07
3D Propagation and Geoacoustic Inversion Studies in the Mid-Atlantic Bight Kevin B. Smith Code PH/Sk, Department of Physics Naval Postgraduate...properties and measured transmission loss. Results from this analysis will be considered in the context of geoacoustic inversions . OBJECTIVES To...bathymetric features and ocean fronts near the shelf break of the mid-Atlantic Bight, and use of various data for geoacoutic inversion studies. The results
Rodriguez, Brian D.
2017-03-31
This report summarizes the results of three-dimensional (3-D) resistivity inversion simulations that were performed to account for local 3-D distortion of the electric field in the presence of 3-D regional structure, without any a priori information on the actual 3-D distribution of the known subsurface geology. The methodology used a 3-D geologic model to create a 3-D resistivity forward (“known”) model that depicted the subsurface resistivity structure expected for the input geologic configuration. The calculated magnetotelluric response of the modeled resistivity structure was assumed to represent observed magnetotelluric data and was subsequently used as input into a 3-D resistivity inverse model that used an iterative 3-D algorithm to estimate 3-D distortions without any a priori geologic information. A publicly available inversion code, WSINV3DMT, was used for all of the simulated inversions, initially using the default parameters, and subsequently using adjusted inversion parameters. A semiautomatic approach of accounting for the static shift using various selections of the highest frequencies and initial models was also tested. The resulting 3-D resistivity inversion simulation was compared to the “known” model and the results evaluated. The inversion approach that produced the lowest misfit to the various local 3-D distortions was an inversion that employed an initial model volume resistivity that was nearest to the maximum resistivities in the near-surface layer.
NASA Astrophysics Data System (ADS)
Simutė, S.; Fichtner, A.
2015-12-01
We present a feasibility study for seismic source inversions using a 3-D velocity model for the Japanese Islands. The approach involves numerically calculating 3-D Green's tensors, which is made efficient by exploiting Green's reciprocity. The rationale for 3-D seismic source inversion has several aspects. For structurally complex regions, such as the Japan area, it is necessary to account for 3-D Earth heterogeneities to prevent unknown structure polluting source solutions. In addition, earthquake source characterisation can serve as a means to delineate existing faults. Source parameters obtained for more realistic Earth models can then facilitate improvements in seismic tomography and early warning systems, which are particularly important for seismically active areas, such as Japan. We have created a database of numerically computed 3-D Green's reciprocals for a 40°× 40°× 600 km size area around the Japanese Archipelago for >150 broadband stations. For this we used a regional 3-D velocity model, recently obtained from full waveform inversion. The model includes attenuation and radial anisotropy and explains seismic waveform data for periods between 10 - 80 s generally well. The aim is to perform source inversions using the database of 3-D Green's tensors. As preliminary steps, we present initial concepts to address issues that are at the basis of our approach. We first investigate to which extent Green's reciprocity works in a discrete domain. Considering substantial amounts of computed Green's tensors we address storage requirements and file formatting. We discuss the importance of the initial source model, as an intelligent choice can substantially reduce the search volume. Possibilities to perform a Bayesian inversion and ways to move to finite source inversion are also explored.
Spontaneous mirror-symmetry breaking induces inverse energy cascade in 3D active fluids
Słomka, Jonasz; Dunkel, Jörn
2017-01-01
Classical turbulence theory assumes that energy transport in a 3D turbulent flow proceeds through a Richardson cascade whereby larger vortices successively decay into smaller ones. By contrast, an additional inverse cascade characterized by vortex growth exists in 2D fluids and gases, with profound implications for meteorological flows and fluid mixing. The possibility of a helicity-driven inverse cascade in 3D fluids had been rejected in the 1970s based on equilibrium-thermodynamic arguments. Recently, however, it was proposed that certain symmetry-breaking processes could potentially trigger a 3D inverse cascade, but no physical system exhibiting this phenomenon has been identified to date. Here, we present analytical and numerical evidence for the existence of an inverse energy cascade in an experimentally validated 3D active fluid model, describing microbial suspension flows that spontaneously break mirror symmetry. We show analytically that self-organized scale selection, a generic feature of many biological and engineered nonequilibrium fluids, can generate parity-violating Beltrami flows. Our simulations further demonstrate how active scale selection controls mirror-symmetry breaking and the emergence of a 3D inverse cascade. PMID:28193853
3D Airborne Electromagnetic Inversion: A case study from the Musgrave Region, South Australia
NASA Astrophysics Data System (ADS)
Cox, L. H.; Wilson, G. A.; Zhdanov, M. S.; Sunwall, D. A.
2012-12-01
Geophysicists know and accept that geology is inherently 3D, and is resultant from complex, overlapping processes related to genesis, metamorphism, deformation, alteration, weathering, and/or hydrogeology. Yet, the geophysics community has long relied on qualitative analysis, conductivity depth imaging (CDIs), 1D inversion, and/or plate modeling. There are many reasons for this deficiency, not the least of which has been the lack of capacity for historic 3D AEM inversion algorithms to invert entire surveys so as to practically affect exploration decisions. Our recent introduction of a moving sensitivity domain (footprint) methodology has been a paradigm shift in AEM interpretation. The basis of this method is that one needs only to calculate the responses and sensitivities for that part of the 3D earth model that is within the AEM system's sensitivity domain (footprint), and then superimpose all sensitivity domains into a single, sparse sensitivity matrix for the entire 3D earth model which is then updated in a regularized inversion scheme. This has made it practical to rigorously invert entire surveys with thousands of line kilometers of AEM data to mega-cell 3D models in hours using multi-processor workstations. Since 2010, over eighty individual projects have been completed for Aerodat, AEROTEM, DIGHEM, GEOTEM, HELITEM, HoisTEM, MEGATEM, RepTEM, RESOLVE, SkyTEM, SPECTREM, TEMPEST, and VTEM data from Australia, Brazil, Canada, Finland, Ghana, Peru, Tanzania, the US, and Zambia. Examples of 3D AEM inversion have been published for a variety of applications, including mineral exploration, oil sands exploration, salinity, permafrost, and bathymetry mapping. In this paper, we present a comparison of 3D inversions for SkyTEM, SPECTREM, TEMPET and VTEM data acquired over the same area in the Musgrave region of South Australia for exploration under cover.
Lithologic identification & mapping test based on 3D inversion of magnetic and gravity
NASA Astrophysics Data System (ADS)
Yan, Jiayong; Lv, Qingtian; Qi, Guang; Zhao, Jinhua; Zhang, Yongqian
2016-04-01
Though lithologic identification & mapping to achieve ore concentration district transparent within 5km depth is the main way to realize deep fine structures study, to explore deep mineral resources and to reveal metallogenic regularity of large-scale ore district . Owing to the wide covered area, high sampling density and mature three-dimensional inversion algorithm of gravity and magnetic data, so gravity and magnetic inversion become the most likely way to achieve three-dimensional lithologic mapping at the present stage. In this paper, we take Lu-zong(Lujiang county to Zongyang county in Anhui province ,east China) ore district as a case, we proposed lithologic mapping flow based 3D inversion of gravity magnetic and then carry out the lithologic mapping test. Lithologic identification & mapping flow is as follows: 1. Analysis relations between lithology and density and magnetic susceptibility by cross plot. 2.Extracting appropriate residual anomalies from high-precision Bourger gravity and aeromagnetic. 3.Use same mesh, do 3D magnetic and gravity inversion respectively under prior information constrained, and then invert susceptibility and density 3D model. 4. According setp1, construct logical topology operations between density 3D model and susceptibility. 5.Use the logical operations, identify lithogies cell by cell in 3D mesh, and then get 3D lithological model. According this flow, we obtained three-dimensional distribution of five main type lithologies in the Lu-Zong ore district within 5km depth. The result of lithologic mapping not only showed that the shallow characteristics and surface geological mapping are basically Coincide,more importantly ,it reveals the deeper lithologic changes.The lithlogical model make up the insufficient of surface geological mapping. The lithologic mapping test results in Lu-Zong ore concentration district showed that lithological mapping using 3D inversion of gravity and magnetic is a effective method to reveal the
2010-09-01
lithosphere elude us. We have been able to surmise that geologic variations here are substantial, and we know that they frustrate attempts to use robust...concepts are summarized conceptually in Figure 2, which shows the regions of the lithosphere most sensitive to the different data that we employ. To...construct an approximate 3D model of the lithosphere , we use a hybrid 1D-3D inversion. In many tomography analyses, dispersion variations are
Inverse cascades sustained by the transfer rate of angular momentum in a 3D turbulent flow.
López-Caballero, Miguel; Burguete, Javier
2013-03-22
The existence of energy cascades as signatures of conserved magnitudes is one of the universal characteristics of turbulent flows. In homogeneous 3D turbulence, the energy conservation produces a direct cascade from large to small scales, although in 2D, it produces an inverse cascade pointing towards small wave numbers. In this Letter, we present the first evidence of an inverse cascade in a fully developed 3D experimental turbulent flow where the conserved magnitude is the angular momentum. Two counterrotating flows collide in a central region where very large fluctuations are produced, generating a turbulent drag that transfers the external torque between different fluid layers.
Inverse Cascades Sustained by the Transfer Rate of Angular Momentum in a 3D Turbulent Flow
NASA Astrophysics Data System (ADS)
López-Caballero, Miguel; Burguete, Javier
2013-03-01
The existence of energy cascades as signatures of conserved magnitudes is one of the universal characteristics of turbulent flows. In homogeneous 3D turbulence, the energy conservation produces a direct cascade from large to small scales, although in 2D, it produces an inverse cascade pointing towards small wave numbers. In this Letter, we present the first evidence of an inverse cascade in a fully developed 3D experimental turbulent flow where the conserved magnitude is the angular momentum. Two counterrotating flows collide in a central region where very large fluctuations are produced, generating a turbulent drag that transfers the external torque between different fluid layers.
Laplace-domain waveform modeling and inversion for the 3D acoustic-elastic coupled media
NASA Astrophysics Data System (ADS)
Shin, Jungkyun; Shin, Changsoo; Calandra, Henri
2016-06-01
Laplace-domain waveform inversion reconstructs long-wavelength subsurface models by using the zero-frequency component of damped seismic signals. Despite the computational advantages of Laplace-domain waveform inversion over conventional frequency-domain waveform inversion, an acoustic assumption and an iterative matrix solver have been used to invert 3D marine datasets to mitigate the intensive computing cost. In this study, we develop a Laplace-domain waveform modeling and inversion algorithm for 3D acoustic-elastic coupled media by using a parallel sparse direct solver library (MUltifrontal Massively Parallel Solver, MUMPS). We precisely simulate a real marine environment by coupling the 3D acoustic and elastic wave equations with the proper boundary condition at the fluid-solid interface. In addition, we can extract the elastic properties of the Earth below the sea bottom from the recorded acoustic pressure datasets. As a matrix solver, the parallel sparse direct solver is used to factorize the non-symmetric impedance matrix in a distributed memory architecture and rapidly solve the wave field for a number of shots by using the lower and upper matrix factors. Using both synthetic datasets and real datasets obtained by a 3D wide azimuth survey, the long-wavelength component of the P-wave and S-wave velocity models is reconstructed and the proposed modeling and inversion algorithm are verified. A cluster of 80 CPU cores is used for this study.
NASA Astrophysics Data System (ADS)
Gumilar, Irwan; Abidin, Hasanuddin Z.; Putra, Andreas D.; Haerani, Nia
2015-04-01
Indonesia is a country with many volcanoes. Each volcano in Indonesia typically has its own crater characteristics. One of them is the Mt.Talaga Bodas, located in Garut, West Java. Researches regarding the crater characteristics are necessary for volcanic disaster mitigation process. One of them is the modelling of the shape of the crater. One of the methods that can be used to model the volcanic crater is using Terrestrial Laser Scanner (TLS). This research aims to create a 3 dimensional (3D) model of the crater of the Mt. Talaga Bodas, that hopefully can be utilized for volcanic disaster mitigation. The methodology used in this research is by obtaining the scanning data using TLS and GPS measurements to obtain the coordinates of the reference points. The data processing methods consist of several steps, namely target to target registration, filterization, georeference, meshing point cloud, surface making, drawing, and 3D modelling. These steps were done using the Cyclone 7 software, and also using 3DS MAX for 3D modelling. The result of this data processing is a 3D model of the crater of the Mt. Talaga Bodas which is similar with the real shape. The calculation result shows that the height of the crater is 62.522 m, the diameter of the crater is 467.231 m, and the total area is 2961054.652 m2. The main obstacle in this research is the dense vegetation which becomes the noise and affects the crater model.
Development of direct-inverse 3-D methods for applied aerodynamic design and analysis
NASA Technical Reports Server (NTRS)
Carlson, Leland A.
1988-01-01
Several inverse methods have been compared and initial results indicate that differences in results are primarily due to coordinate systems and fuselage representations and not to design procedures. Further, results from a direct-inverse method that includes 3-D wing boundary layer effects, wake curvature, and wake displacement are presented. These results show that boundary layer displacements must be included in the design process for accurate results.
Research on Joint Parameter Inversion for an Integrated Underground Displacement 3D Measuring Sensor
Shentu, Nanying; Qiu, Guohua; Li, Qing; Tong, Renyuan; Shentu, Nankai; Wang, Yanjie
2015-01-01
Underground displacement monitoring is a key means to monitor and evaluate geological disasters and geotechnical projects. There exist few practical instruments able to monitor subsurface horizontal and vertical displacements simultaneously due to monitoring invisibility and complexity. A novel underground displacement 3D measuring sensor had been proposed in our previous studies, and great efforts have been taken in the basic theoretical research of underground displacement sensing and measuring characteristics by virtue of modeling, simulation and experiments. This paper presents an innovative underground displacement joint inversion method by mixing a specific forward modeling approach with an approximate optimization inversion procedure. It can realize a joint inversion of underground horizontal displacement and vertical displacement for the proposed 3D sensor. Comparative studies have been conducted between the measured and inversed parameters of underground horizontal and vertical displacements under a variety of experimental and inverse conditions. The results showed that when experimentally measured horizontal displacements and vertical displacements are both varied within 0 ~ 30 mm, horizontal displacement and vertical displacement inversion discrepancies are generally less than 3 mm and 1 mm, respectively, under three kinds of simulated underground displacement monitoring circumstances. This implies that our proposed underground displacement joint inversion method is robust and efficient to predict the measuring values of underground horizontal and vertical displacements for the proposed sensor. PMID:25871714
3-D Inverse Teleseismic Scattered Wave Imaging using the Kirchhoff Approximation
NASA Astrophysics Data System (ADS)
Liu, K.; Levander, A.
2012-04-01
We have developed a 3-D teleseismic imaging technique for scattered elastic wavefields using the Kirchhoff approximation. Kirchhoff migration/inversion have been well developed in exploration seismology within the inverse scattering framework (e.g. Miller et al., 1987; Beylkin and Burridge, 1990) to image subsurface structure that generates secondary wavefields caused by localized heterogeneities. Application of this method in global seismology has been largely limited to 2-D images made with 1-D reference models due to high computational cost and the lack of adequately dense receiver arrays (Bostock, 2002, Poppeliers and Pavlis, 2003; Frederiksen and Revenaugh, 2004; Cao et al., 2010). The deployment of the USArray Transportable and Flexible arrays in the United States and dense array recordings in other countries motivate developing teleseismic scattered wavefield imaging with the Kirchhoff approximation for 3-D velocity models for both scalar and vector wavefields to improve upper mantle imaging. Following Bostock's development of the 2-D problem (2002), we derive the 3-D P-to-S scattering inversion formula by phrasing the inverse problem in terms of the generalized Radon transform (GRT) and singular functions of discontinuity surfaces. In the forward scattering modeling, we extend the method to utilize a 3-D migration velocity model by calculating 3-D finite-difference traveltimes, backprojected from the receivers using an eikonal solver. To demonstrate the relative accuracy of the inversion, we examine several synthetic cases with a variety of discontinuity surfaces (sinuous, dipping, dome- and crater-shaped discontinuity interfaces, point scatterers, etc.). The Kirchhoff GRT imaging can successfully recover the shapes of these structures very well. We compare our Kirchhoff approximation imaging with the Born-approximate results, as well as the common-conversion point (CCP) stacked receiver function imaging for the various synthetic cases, and show a field
3D LBFGS inversion of controlled source extremely low frequency electromagnetic data
NASA Astrophysics Data System (ADS)
Cao, Meng; Tan, Han-Dong; Wang, Kun-Peng
2016-12-01
The controlled source extremely low frequency (CSELF) electromagnetic method is characterized by extremely long and powerful sources and a huge measurement range. Its electromagnetic field can therefore be affected by the ionosphere and displacement current. Research on 3D forward modeling and inversion of CSELF electromagnetic data is currently in its infancy. This paper makes exploratory attempts to firstly calculate the 1D extremely low frequency electromagnetic field under ionosphere-air-earth coupling circumstances, and secondly analyze the propagation characteristics of the background electromagnetic field. The 3D staggered-grid finite difference scheme for solving for the secondary electric field is adopted and incorporated with the 1D modeling algorithm to complete 3D forward modeling. Considering that surveys can be carried out in the near field and transition zone for lower frequencies, the 3D Limited-memory Broyden-Fletcher-Goldfarb-Shanno (LBFGS) inversion of CSELF electromagnetic data is presented (in which the sources, or primary fields, are included), with the aim of directly inverting the impedance data, regardless of where it is acquired. Derivation of the objective functional gradient is the core component in the inversion. Synthetic tests indicate that the well-chosen approximation to the Hessian can significantly speed up the inversion. The model responses corresponding to the coexistence of conductive and resistive blocks show that the off-diagonal components of tensor impedance are much more sensitive to the resistivity variation than the diagonal components. In comparison with conventional scalar inversion, tensor inversion is superior in the recoveries of electric anomalies and background resistivity.
3D inversion of aeromagnetic Data on Las Tablas District, Panama
NASA Astrophysics Data System (ADS)
Batista-Rodríguez, José A.; Caballero, Alberto; Pérez-Flores, Marco A.; Almaguer-Carmenates, Yuri
2017-03-01
We present a 3D model of Las Tablas District, Panama, obtained from the 3D inversion of aeromagnetic data, and constrained with information from surface geology, water wells and topography. The 3D model suggests the location, boundary, shape and depths of the sedimentary basin where the Mensabé and Salados rivers hydrogeological sub-basin is located. The model shows the connections between tectonics and the sedimentary basin, suggesting the probable areas for aquifers, the relations between them, their zone of recharge and discharge, and the probable zone of pollution. The inferred faults in the model may be the main recharge and discharge conduits for the groundwater and anthropogenic pollution. The geological and geometric characteristics shown in the 3D model are fundamental data for further hydrogeological and geophysical studies such as the location for future drinking water wells.
3-D wavelet compression and progressive inverse wavelet synthesis rendering of concentric mosaic.
Luo, Lin; Wu, Yunnan; Li, Jin; Zhang, Ya-Qin
2002-01-01
Using an array of photo shots, the concentric mosaic offers a quick way to capture and model a realistic three-dimensional (3-D) environment. We compress the concentric mosaic image array with a 3-D wavelet transform and coding scheme. Our compression algorithm and bitstream syntax are designed to ensure that a local view rendering of the environment requires only a partial bitstream, thereby eliminating the need to decompress the entire compressed bitstream before rendering. By exploiting the ladder-like structure of the wavelet lifting scheme, the progressive inverse wavelet synthesis (PIWS) algorithm is proposed to maximally reduce the computational cost of selective data accesses on such wavelet compressed datasets. Experimental results show that the 3-D wavelet coder achieves high-compression performance. With the PIWS algorithm, a 3-D environment can be rendered in real time from a compressed dataset.
NASA Astrophysics Data System (ADS)
Rosas-Carbajal, M.; Linde, N.; Peacock, J.; Zyserman, F. I.; Kalscheuer, T.; Thiel, S.
2015-12-01
Surface-based monitoring of mass transfer caused by injections and extractions in deep boreholes is crucial to maximize oil, gas and geothermal production. Inductive electromagnetic methods, such as magnetotellurics, are appealing for these applications due to their large penetration depths and sensitivity to changes in fluid conductivity and fracture connectivity. In this work, we propose a 3-D Markov chain Monte Carlo inversion of time-lapse magnetotelluric data to image mass transfer following a saline fluid injection. The inversion estimates the posterior probability density function of the resulting plume, and thereby quantifies model uncertainty. To decrease computation times, we base the parametrization on a reduced Legendre moment decomposition of the plume. A synthetic test shows that our methodology is effective when the electrical resistivity structure prior to the injection is well known. The centre of mass and spread of the plume are well retrieved. We then apply our inversion strategy to an injection experiment in an enhanced geothermal system at Paralana, South Australia, and compare it to a 3-D deterministic time-lapse inversion. The latter retrieves resistivity changes that are more shallow than the actual injection interval, whereas the probabilistic inversion retrieves plumes that are located at the correct depths and oriented in a preferential north-south direction. To explain the time-lapse data, the inversion requires unrealistically large resistivity changes with respect to the base model. We suggest that this is partly explained by unaccounted subsurface heterogeneities in the base model from which time-lapse changes are inferred.
Anisotropic 3D inversion of towed streamer EM data from the Troll West oil province (Invited)
NASA Astrophysics Data System (ADS)
Mattsson, J.; Midgley, J.; Zhdanov, M. S.; ENDO, M.
2013-12-01
Obviating the need for ocean bottom receivers, the towed streamer EM system enables CSEM data to be acquired simultaneously with seismic over very large areas in frontier and mature basins for higher production rates and more cost effective than conventional marine CSEM. The towed streamer EM data are currently processed and delivered as a spectrum of frequency-domain responses. We apply a 3D anisotropic inversion methodology for towed streamer EM data that includes a moving sensitivity domain. Our implementation is based on the 3D contraction integral equation method for computing the EM responses and Fréchet derivatives, and uses the re-weighted regularized conjugate gradient method for minimizing the objective functional with focusing regularization. We present an actual case study for the 3D anisotropic inversion of towed streamer EM data from the Troll West oil province in the North Sea, and demonstrate our ability to image the Troll West Oil and Gas Provinces. We conclude that 3D anisotropic inversion of the data from the current generation of towed streamer EM system can adequately recover both the vertical and horizontal resistivities in anisotropic hydrocarbon-bearing formations.
3-D inversion of gravity data in spherical coordinates with application to the GRAIL data
NASA Astrophysics Data System (ADS)
Liang, Qing; Chen, Chao; Li, Yaoguo
2014-06-01
Three-dimensional (3-D) inversion of gravity data has been widely used to reconstruct the density distributions of ore bodies, basins, crust, lithosphere, and upper mantle. For global model of 3-D density structures of planetary interior, such as the Earth, the Moon, or Mars, it is necessary to use an inversion algorithm that operates in the spherical coordinates. We develop a 3-D inversion algorithm formulated with specially designed model objective function and radial weighting function in the spherical coordinates. We present regional and global synthetic examples to illustrate the capability of the algorithm. The inverted results show density distribution features consistent with the true models. We also apply the algorithm to a set of lunar Bouguer gravity anomaly derived from the Gravity Recovery and Interior Laboratory (GRAIL) gravity field and obtain a lunar 3-D density distribution. High-density anomalies are clearly identified underlying lunar basins, a wide region of the lateral density heterogeneities that exist beneath the South Pole-Aitken basin are found, and low-density anomalies are distributed beneath the Feldspathic Highlands Terrane on the lunar far-side. The consistency of these results with those obtained independently from other existing methods verifies the newly developed algorithm.
3D gravity inversion and uncertainty assessment of basement relief via Particle Swarm Optimization
NASA Astrophysics Data System (ADS)
Pallero, J. L. G.; Fernández-Martínez, J. L.; Bonvalot, S.; Fudym, O.
2017-04-01
Nonlinear gravity inversion in sedimentary basins is a classical problem in applied geophysics. Although a 2D approximation is widely used, 3D models have been also proposed to better take into account the basin geometry. A common nonlinear approach to this 3D problem consists in modeling the basin as a set of right rectangular prisms with prescribed density contrast, whose depths are the unknowns. Then, the problem is iteratively solved via local optimization techniques from an initial model computed using some simplifications or being estimated using prior geophysical models. Nevertheless, this kind of approach is highly dependent on the prior information that is used, and lacks from a correct solution appraisal (nonlinear uncertainty analysis). In this paper, we use the family of global Particle Swarm Optimization (PSO) optimizers for the 3D gravity inversion and model appraisal of the solution that is adopted for basement relief estimation in sedimentary basins. Synthetic and real cases are illustrated, showing that robust results are obtained. Therefore, PSO seems to be a very good alternative for 3D gravity inversion and uncertainty assessment of basement relief when used in a sampling while optimizing approach. That way important geological questions can be answered probabilistically in order to perform risk assessment in the decisions that are made.
Fast 3D inversion of airborne gravity-gradiometry data using Lanczos bidiagonalization method
NASA Astrophysics Data System (ADS)
Meng, Zhaohai; Li, Fengting; Zhang, Dailei; Xu, Xuechun; Huang, Danian
2016-09-01
We developed a new fast inversion method for to process and interpret airborne gravity gradiometry data, which was based on Lanczos bidiagonalization algorithm. Here, we describe the application of this new 3D gravity gradiometry inversion method to recover a subsurface density distribution model from the airborne measured gravity gradiometry anomalies. For this purpose, the survey area is divided into a large number of rectangular cells with each cell possessing a constant unknown density. It is well known that the solution of large linear gravity gradiometry is an ill-posed problem since using the smoothest inversion method is considerably time consuming. We demonstrate that the Lanczos bidiagonalization method can be an appropriate algorithm to solve a Tikhonov solver time cost function for resolving the large equations within a short time. Lanczos bidiagonalization is designed to make the very large gravity gradiometry forward modeling matrices to become low-rank, which will considerably reduce the running time of the inversion method. We also use a weighted generalized cross validation method to choose the appropriate Tikhonov parameter to improve inversion results. The inversion incorporates a model norm that allows us to attain the smoothing and depth of the solution; in addition, the model norm counteracts the natural decay of the kernels, which concentrate at shallow depths. The method is applied on noise-contaminated synthetic gravity gradiometry data to demonstrate its suitability for large 3D gravity gradiometry data inversion. The airborne gravity gradiometry data from the Vinton Salt Dome, USE, were considered as a case study. The validity of the new method on real data is discussed with reference to the Vinton Dome inversion result. The intermediate density values in the constructed model coincide well with previous results and geological information. This demonstrates the validity of the gravity gradiometry inversion method.
Fast 3D inversion of gravity data using solution space priorconditioned lanczos bidiagonalization
NASA Astrophysics Data System (ADS)
Rezaie, Mohammad; Moradzadeh, Ali; Kalateh, Ali Nejati
2017-01-01
Inversion of gravity data is one of the most important steps in the quantitative interpretation of practical data. Inversion is a mathematical technique that automatically constructs a subsurface geophysical model from measured data, incorporating some priori information. Inversion of gravity data is time consuming because of increase in data and model parameters. Some efforts have been made to deal with this problem, one of them is using fast algorithms for solving system of equations in inverse problem. Lanczos bidiagonalization method is a fast algorithm that works based on Krylov subspace iterations and projection method, but cannot always provide a good basis for a projection method. So in this study, we combined the Krylov method with a regularization method applied to the low-dimensional projected problem. To achieve the goal, the orthonormal basis vectors of the discrete cosine transform (DCT) were used to build the low-dimensional subspace. The forward operator matrix replaced with a matrix of lower dimension, thus, the required memory and running time of the inverse modeling is decreased by using the proposed algorithm. It is shown that this algorithm can be appropriate to solve a Tikhonov cost function for inversion of gravity data. The proposed method has been applied on a noise-corrupted synthetic data and field gravity data (Mobrun gravity data) to demonstrate its reliability for three dimensional (3D) gravity inversion. The obtained results of 3D inversion both synthetic and field gravity data (Mobrun gravity data) indicate the proposed inversion algorithm could produce density models consistent with true structures.
NASA Astrophysics Data System (ADS)
Ogaya, X.; Campanya i Llovet, J.; Jones, A. G.; Rath, V.; Vozar, J.
2015-12-01
In this work we assess the main advantages of complementing the standard magnetotelluric (MT) impedance tensor (Z) response with inter-station horizontal magnetic tensor (HMT) and geomagnetic transfer function (GTF) responses. Results were obtained from constraining the subsurface beneath a MT profile in a 3-D environment using distorted synthetic responses and real data. From the synthetic case study four main improvements were observed when performing joint inversion of Z, HMT and GTF data simultaneously, instead of only inverting the commonly used Zdata: (1) superior precision in the characterization of the electrical resistivity values of the anomalies below and outside the profile; (2) the potential to recover high electrical resistivity anomalies that are poorly recovered by using MT responses alone; (3) improvement in the characterization of the bottom and lateral boundaries of the anomalies with low electrical resistivity; and (4) superior imaging of the continuity of structures with low electrical resistivity. A corresponding inversion comparison was performed using real responses from Clare Basin (Co. Clare, Ireland) acquired during summer of 2014 within the IRECCSEM project (www.ireccsem.ie). A total of 140 sites were acquired, including audiomagnetotellurics (AMT) and broadband magnetotelluric (BBMT) time series data. The results obtained undertaking a joint inversion of Z, GTF and HMT data corroborate the results observed in the synthetic study. In addition, when confronted the results of the Clare Basin area with independent geological and geophysical data, the geoelectrical models show higher resolution when HMT and GTF data are included in the inversion process. The observed advantages offer new opportunities for the MT method by making the results from a MT profile in a 3-D environment more believable, supporting the possibility of high-resolution studies without necessarily expending a large amount of economical and computational resources.
Large scale 3-D modeling by integration of resistivity models and borehole data through inversion
NASA Astrophysics Data System (ADS)
Foged, N.; Marker, P. A.; Christansen, A. V.; Bauer-Gottwein, P.; Jørgensen, F.; Høyer, A.-S.; Auken, E.
2014-02-01
We present an automatic method for parameterization of a 3-D model of the subsurface, integrating lithological information from boreholes with resistivity models through an inverse optimization, with the objective of further detailing for geological models or as direct input to groundwater models. The parameter of interest is the clay fraction, expressed as the relative length of clay-units in a depth interval. The clay fraction is obtained from lithological logs and the clay fraction from the resistivity is obtained by establishing a simple petrophysical relationship, a translator function, between resistivity and the clay fraction. Through inversion we use the lithological data and the resistivity data to determine the optimum spatially distributed translator function. Applying the translator function we get a 3-D clay fraction model, which holds information from the resistivity dataset and the borehole dataset in one variable. Finally, we use k means clustering to generate a 3-D model of the subsurface structures. We apply the concept to the Norsminde survey in Denmark integrating approximately 700 boreholes and more than 100 000 resistivity models from an airborne survey in the parameterization of the 3-D model covering 156 km2. The final five-cluster 3-D model differentiates between clay materials and different high resistive materials from information held in resistivity model and borehole observations respectively.
Large-scale 3-D modeling by integration of resistivity models and borehole data through inversion
NASA Astrophysics Data System (ADS)
Foged, N.; Marker, P. A.; Christansen, A. V.; Bauer-Gottwein, P.; Jørgensen, F.; Høyer, A.-S.; Auken, E.
2014-11-01
We present an automatic method for parameterization of a 3-D model of the subsurface, integrating lithological information from boreholes with resistivity models through an inverse optimization, with the objective of further detailing of geological models, or as direct input into groundwater models. The parameter of interest is the clay fraction, expressed as the relative length of clay units in a depth interval. The clay fraction is obtained from lithological logs and the clay fraction from the resistivity is obtained by establishing a simple petrophysical relationship, a translator function, between resistivity and the clay fraction. Through inversion we use the lithological data and the resistivity data to determine the optimum spatially distributed translator function. Applying the translator function we get a 3-D clay fraction model, which holds information from the resistivity data set and the borehole data set in one variable. Finally, we use k-means clustering to generate a 3-D model of the subsurface structures. We apply the procedure to the Norsminde survey in Denmark, integrating approximately 700 boreholes and more than 100 000 resistivity models from an airborne survey in the parameterization of the 3-D model covering 156 km2. The final five-cluster 3-D model differentiates between clay materials and different high-resistivity materials from information held in the resistivity model and borehole observations, respectively.
NASA Astrophysics Data System (ADS)
Klibanov, Michael V.; Romanov, Vladimir G.
2016-01-01
The 3D inverse scattering problem of the reconstruction of the unknown dielectric permittivity in the generalized Helmholtz equation is considered. Applications are in imaging of nanostructures and biological cells. The main difference with the conventional inverse scattering problems is that only the modulus of the scattering wave field is measured. The phase is not measured. The initializing wave field is the incident plane wave. On the other hand, in the previous recent works of the authors about the ‘phaseless topic’ the case of the point source was considered (Klibanov and Romanov 2015 J. Inverse Ill-Posed Problem 23 415-28 J. Inverse Ill-Posed Problem 23 187-93). Two reconstruction procedures are developed.
3D Motion Planning Algorithms for Steerable Needles Using Inverse Kinematics
Duindam, Vincent; Xu, Jijie; Alterovitz, Ron; Sastry, Shankar; Goldberg, Ken
2010-01-01
Steerable needles can be used in medical applications to reach targets behind sensitive or impenetrable areas. The kinematics of a steerable needle are nonholonomic and, in 2D, equivalent to a Dubins car with constant radius of curvature. In 3D, the needle can be interpreted as an airplane with constant speed and pitch rate, zero yaw, and controllable roll angle. We present a constant-time motion planning algorithm for steerable needles based on explicit geometric inverse kinematics similar to the classic Paden-Kahan subproblems. Reachability and path competitivity are analyzed using analytic comparisons with shortest path solutions for the Dubins car (for 2D) and numerical simulations (for 3D). We also present an algorithm for local path adaptation using null-space results from redundant manipulator theory. Finally, we discuss several ways to use and extend the inverse kinematics solution to generate needle paths that avoid obstacles. PMID:21359051
NASA Astrophysics Data System (ADS)
Bell, R. E.; Morgan, J. V.; Warner, M.
2013-12-01
Our understanding of subduction margin seismogenesis has been revolutionised in the last couple of decades with the discovery that the size of the seismogenic zone may not be controlled simply by temperature and a broad spectrum of seismic behaviour exists from stick-slip to stable sliding. Laboratory and numerical experiments suggest that physical properties, particularly fluid pressure may play an important role in controlling the seismic behaviour of subduction margins. Although drilling can provide information on physical properties along subduction thrust faults at point locations at relatively shallow depths, correlations between physical properties and seismic velocity using rock physics relationships are required to resolve physical properties along the margin and down-dip. Therefore, high resolution seismic velocity models are key to recovering physical property information at subduction plate boundaries away from drill sites. 3D Full waveform inversion (FWI) is a technique pioneered by the oil industry to obtain high-resolution high-fidelity models of physical properties in the sub-surface. 3D FWI involves the inversion of low-frequency (>2 to <7 Hz), early arriving (principally transmitted) seismic data, to recover the macro (intermediate to long-wavelength) velocity structure. Although 2D FWI has been used to improve velocity models of subduction plate boundaries before, 3D FWI has not yet been attempted. 3D inversions have superior convergence and accuracy, as they sample the subsurface with multi-azimuth multiply-crossing wavefields. In this contribution we perform a suite of synthetic tests to investigate if 3D FWI could be used to better resolve physical property information along subduction margin plate boundaries using conventionally collected 3D seismic data. We base our analysis on the Muroto Basin area of the Nankai margin and investigate if the acquisition parameters and geometry of the subduction margin render 3D seismic data collected across
3-D Sound Propagation and Acoustic Inversions in Shallow Water Oceans
2011-09-01
method is used to study canonical environmental models of shelfbreak front systems and nonlinear internal wave ducts. The WHOI 3D Parabolic-Equation...localization methods with normal mode theory have been established for localizing low frequency, broadband signals in a shallow water environment. Gauss ...approach for low-frequency broadband sound source localization in a shallow-water ocean is established. Gauss -Markov inverse theory is used in both
Sensitivity study of 3-D modeling for multi-D inversion of surface NMR
NASA Astrophysics Data System (ADS)
Warsa, Grandis, Hendra
2012-06-01
Geophysical field method of surface nuclear magnetic resonance (SNMR) allows a direct determination of hydrogeological parameters of the subsurface. The amplitude of the SNMR signal is directly linked to the amount of mobile water. The relaxation behaviour of the signal correlates with pore sizes and hydraulic conductivities of an aquifer. For improving capability and reliability of SNMR method we have presented a forward modeling scheme of 3-D water content and decay time structures that can be used for multi-D interpretation. Currently SNMR is carried out mainly with a 1-D working scheme using coinciding loops. For each sounding point using a coincident circular loop antenna, the amplitudes and decay times of the SNMR signal are the product of a three dimensional distribution of the water content and decay time in the subsurface and their sensitivity to the receiver. The antenna is moved at the surface and the SNMR relaxation signal are plotted as a function of the pulse moment and sounding point. The errors might be very large by neglecting the 2-D or even 3-D geometry of the structures which have to be considered in the analysis and inversion in the future. The results show that the 3-D modeling is reliable and flexible to be integrated into the 2-D/3-D inversion scheme for inverting surface NMR data to recover a multi-D distribution of water content and decay time of an aquifer.
NASA Astrophysics Data System (ADS)
Liu, B.; Li, S. C.; Nie, L. C.; Wang, J.; L, X.; Zhang, Q. S.
2012-12-01
Traditional inversion method is the most commonly used procedure for three-dimensional (3D) resistivity inversion, which usually takes the linearization of the problem and accomplish it by iterations. However, its accuracy is often dependent on the initial model, which can make the inversion trapped in local optima, even cause a bad result. Non-linear method is a feasible way to eliminate the dependence on the initial model. However, for large problems such as 3D resistivity inversion with inversion parameters exceeding a thousand, main challenges of non-linear method are premature and quite low search efficiency. To deal with these problems, we present an improved Genetic Algorithm (GA) method. In the improved GA method, smooth constraint and inequality constraint are both applied on the object function, by which the degree of non-uniqueness and ill-conditioning is decreased. Some measures are adopted from others by reference to maintain the diversity and stability of GA, e.g. real-coded method, and the adaptive adjustment of crossover and mutation probabilities. Then a generation method of approximately uniform initial population is proposed in this paper, with which uniformly distributed initial generation can be produced and the dependence on initial model can be eliminated. Further, a mutation direction control method is presented based on the joint algorithm, in which the linearization method is embedded in GA. The update vector produced by linearization method is used as mutation increment to maintain a better search direction compared with the traditional GA with non-controlled mutation operation. By this method, the mutation direction is optimized and the search efficiency is improved greatly. The performance of improved GA is evaluated by comparing with traditional inversion results in synthetic example or with drilling columnar sections in practical example. The synthetic and practical examples illustrate that with the improved GA method we can eliminate
Kılıç, Emre Eibert, Thomas F.
2015-05-01
An approach combining boundary integral and finite element methods is introduced for the solution of three-dimensional inverse electromagnetic medium scattering problems. Based on the equivalence principle, unknown equivalent electric and magnetic surface current densities on a closed surface are utilized to decompose the inverse medium problem into two parts: a linear radiation problem and a nonlinear cavity problem. The first problem is formulated by a boundary integral equation, the computational burden of which is reduced by employing the multilevel fast multipole method (MLFMM). Reconstructed Cauchy data on the surface allows the utilization of the Lorentz reciprocity and the Poynting's theorems. Exploiting these theorems, the noise level and an initial guess are estimated for the cavity problem. Moreover, it is possible to determine whether the material is lossy or not. In the second problem, the estimated surface currents form inhomogeneous boundary conditions of the cavity problem. The cavity problem is formulated by the finite element technique and solved iteratively by the Gauss–Newton method to reconstruct the properties of the object. Regularization for both the first and the second problems is achieved by a Krylov subspace method. The proposed method is tested against both synthetic and experimental data and promising reconstruction results are obtained.
NASA Astrophysics Data System (ADS)
Cubuk-Sabuncu, Yesim; Taymaz, Tuncay; Fichtner, Andreas
2016-04-01
We present a 3D radially anisotropic velocity model of the crust and uppermost mantle structure beneath the Sea of Marmara and surroundings based on the full waveform inversion method. The intense seismic activity and crustal deformation are observed in the Northwest Turkey due to transition tectonics between the strike-slip North Anatolian Fault (NAF) and the extensional Aegean region. We have selected and simulated complete waveforms of 62 earthquakes (Mw > 4.0) occurred during 2007-2015, and recorded at (Δ < 10°) distances. Three component earthquake data is obtained from broadband seismic stations of Kandilli Observatory and Earthquake Research Center (KOERI, Turkey), Hellenic Unified Seismic Network (HUSN, Greece) and Earthquake Research Center of Turkey (AFAD-DAD). The spectral-element solver of the wave equation, SES3D algorithm, is used to simulate seismic wave propagation in 3D spherical coordinates (Fichtner, 2009). The Large Scale Seismic Inversion Framework (LASIF) workflow tool is also used to perform full seismic waveform inversion (Krischer et al., 2015). The initial 3D Earth model is implemented from the multi-scale seismic tomography study of Fichtner et al. (2013). Discrepancies between the observed and simulated synthetic waveforms are determined using the time-frequency misfits which allows a separation between phase and amplitude information (Fichtner et al., 2008). The conjugate gradient optimization method is used to iteratively update the initial Earth model when minimizing the misfit. The inversion is terminated after 19 iterations since no further advances are observed in updated models. Our analysis revealed shear wave velocity variations of the shallow and deeper crustal structure beneath western Turkey down to depths of ~35-40 km. Low shear wave velocity anomalies are observed in the upper and mid crustal depths beneath major fault zones located in the study region. Low velocity zones also tend to mark the outline of young volcanic
Simultaneous elastic parameter inversion in 2-D/3-D TTI medium combined later arrival times
NASA Astrophysics Data System (ADS)
Bai, Chao-ying; Wang, Tao; Yang, Shang-bei; Li, Xing-wang; Huang, Guo-jiao
2016-04-01
Traditional traveltime inversion for anisotropic medium is, in general, based on a "weak" assumption in the anisotropic property, which simplifies both the forward part (ray tracing is performed once only) and the inversion part (a linear inversion solver is possible). But for some real applications, a general (both "weak" and "strong") anisotropic medium should be considered. In such cases, one has to develop a ray tracing algorithm to handle with the general (including "strong") anisotropic medium and also to design a non-linear inversion solver for later tomography. Meanwhile, it is constructive to investigate how much the tomographic resolution can be improved by introducing the later arrivals. For this motivation, we incorporated our newly developed ray tracing algorithm (multistage irregular shortest-path method) for general anisotropic media with a non-linear inversion solver (a damped minimum norm, constrained least squares problem with a conjugate gradient approach) to formulate a non-linear inversion solver for anisotropic medium. This anisotropic traveltime inversion procedure is able to combine the later (reflected) arrival times. Both 2-D/3-D synthetic inversion experiments and comparison tests show that (1) the proposed anisotropic traveltime inversion scheme is able to recover the high contrast anomalies and (2) it is possible to improve the tomographic resolution by introducing the later (reflected) arrivals, but not as expected in the isotropic medium, because the different velocity (qP, qSV and qSH) sensitivities (or derivatives) respective to the different elastic parameters are not the same but are also dependent on the inclination angle.
3D inversion of land-based CSEM data from the Ketzin CO2 storage formation
NASA Astrophysics Data System (ADS)
Grayver, Alexander; Streich, Rita; Ritter, Oliver
2013-04-01
We present 3D inversion of land controlled-source electromagnetic (CSEM) data collected across the CO2 storage test site at Ketzin, Germany. The CSEM data were generated by injecting currents into the earth at eight locations using a newly developed transmitter equipped with three grounded electrodes. Electric and magnetic field responses were recorded by 39 receivers along a line approximately perpendicular to the main geological trend. The survey aimed at imaging large-scale resistivity structure beyond the near-well region monitored by higher-resolution electrical techniques. Infrastructure present in the survey area, such as pipelines with impressed-current cathodic protection systems, power lines, and wind power plants cause strong noise in the data. The noise is effectively suppressed by adopting statistically robust processing techniques known from passive magnetotellurics. A newly developed Gauss-Newton type parallel distributed inversion scheme, which is based on a direct forward solver and explicitly calculates the full sensitivity matrix, is applied to recover subsurface conductivity images. As 3D inversion is demanding on computer time and memory, we run inversions on parallel distributed machines. We achieve good scalability by distributing computations and memory uniformly among the processes involved. We carry out cumulative sensitivity and resolution analyses for the sparse CSEM acquisition geometry. These studies indicate reasonable spatial coverage along the main survey line. Synthetic studies calculated for the real survey layout and representative conductivity models indicate that the magnetic field components are practically insensitive to resistive structures, whereas the electric field components resolve resistors and conductors similarly well. Because the magnetic field contributes little subsurface information, we concentrate on inverting the electric field, which is also more computer-efficient than inverting all components. We test
The multi-scale 3D-1D compatibility scoring for inverse protein folding problem
Oniuka, Kentaro; Asai, Kiyoshi
1994-12-31
The applicability of the Multi-Scale Structure Description (MSSD) scheme to the inverse-folding problems was investigated. An MSSD represents a 3D protein structure with multiple symbolic sequences, where fine structures are represented with the sequence at low levels, the middle scale structural motifs at middle levels, and global topology at high levels. Each symbol in the symbolic sequence denotes a type of local structure of the level scale. The structure fragments are classified at each scale level respectively according to the shape and the environment around the fragments: how the structure is exposed to the solvent or buried in the molecule. I modeled the propensity of an amino-acid sequence to the structure fragment type (i.e., primary constraint) at each scale level. The local propensity is, therefore, modeled at small scale (low) levels, while the global propensity modeled at large scale (high) levels. Thus, superposing all the primary constraints, a 3D protein structure yields an amino-acid sequence profile. Evaluating the fit of an amino acid sequence to the profile derived from the known 3D protein structure, we can identify which 3D structure the given amino-acid sequence would fold into. I checked whether a sequence identifies its own structure over two hundred protein sequences. In many cases, an amino acid sequence identified its own 3D protein structure.
NASA Astrophysics Data System (ADS)
Janssen, G.; Del Val Alonso, L.; Groenendijk, P.; Griffioen, J.
2012-12-01
We developed an on-line coupling between the 1D/quasi-2D nutrient transport model ANIMO and the 3D groundwater transport model code MT3DMS. ANIMO is a detailed, process-oriented model code for the simulation of nitrate leaching to groundwater, N- and P-loads on surface waters and emissions of greenhouse gasses. It is the leading nutrient fate and transport code in the Netherlands where it is used primarily for the evaluation of fertilization related legislation. In addition, the code is applied frequently in international research projects. MT3DMS is probably the most commonly used groundwater solute transport package worldwide. The on-line model coupling ANIMO-MT3DMS combines the state-of-the-art descriptions of the biogeochemical cycles in ANIMO with the advantages of using a 3D approach for the transport through the saturated domain. These advantages include accounting for regional lateral transport, considering groundwater-surface water interactions more explicitly, and the possibility of using MODFLOW to obtain the flow fields. An additional merit of the on-line coupling concept is that it preserves feedbacks between the saturated and unsaturated zone. We tested ANIMO-MT3DMS by simulating nutrient transport for the period 1970-2007 in a Dutch agricultural polder catchment covering an area of 118 km2. The transient groundwater flow field had a temporal resolution of one day and was calculated with MODFLOW-MetaSWAP. The horizontal resolution of the model grid was 100x100m and consisted of 25 layers of varying thickness. To keep computation times manageable, we prepared MT3DMS for parallel computing, which in itself is a relevant development for a large community of groundwater transport modelers. For the parameterization of the soil, we applied a standard classification approach, representing the area by 60 units with unique combinations of soil type, land use and geohydrological setting. For the geochemical parameterization of the deeper subsurface, however, we
NASA Astrophysics Data System (ADS)
Kiyan, D.; Jones, A. G.; Fullea, J.; Ledo, J.; Siniscalchi, A.; Romano, G.
2013-12-01
The overarching objectives of the second phase of the PICASSO (Program to Investigate Convective Alboran Sea System Overturn) project and the concomitant TopoMed (Plate re-organization in the western Mediterranean: Lithospheric causes and topographic consequences - an ESF EUROSCORES TOPO-EUROPE project) project are (i) to provide new electrical conductivity constraints on the crustal and lithospheric structures of the Atlas Mountains, and (ii) to test the hypotheses for explaining the observation of a 'missing' mantle root inferred from surface heat flow, gravity and geoid anomalies, elevation and seismic data modeling (i.e. Zeyen et al., 2005; Teixell et al., 2005; Fullea et al., 2010). We present the results from three-dimensional (3-D) MT inversion of data from two MT profiles employing the parallel version of Modular system for Electromagnetic inversion (ModEM; Egbert & Kelbert, 2012) code. For the profile in eastern Morocco, passing through Midelt, a distinct conductivity difference between the Middle-High Atlas (conductive) and Anti Atlas (resistive) correlates with the South Atlas Front fault, the depth extent of which appears to be limited to the uppermost mantle (approximately 55 km). In all inverse solutions, the crust and the upper mantle show a resistive signature (750 Ωm - 1,000 Ωm) beneath the Anti Atlas to a depth of 100 km, which is the part of stable West African Craton. Our results are at variance with the proposed thin lithosphere beneath the Middle-High Atlas as we see no evidence for a shallow asthenosphere. Our second profile lies in western Morocco traversing through Marrakech. For the first time, the electrical resistivity distribution in the crust and in the upper mantle of Western High Atlas has been studied. Our 3-D resistivity model shows that conductive (1-20 Ωm) western High Atlas is confined by two resistive basins (>1,000 Ωm), Souss basin to the south and Houz basin to the north. At the southern boundary of the western High Atlas
Integrated gravity and gravity gradient 3D inversion using the non-linear conjugate gradient
NASA Astrophysics Data System (ADS)
Qin, Pengbo; Huang, Danian; Yuan, Yuan; Geng, Meixia; Liu, Jie
2016-03-01
Gravity data, which are critical in mineral, oil, and gas exploration, are obtained from the vertical component of the gravity field, while gravity gradient data are measured from changes in the gravity field in three directions. However, few studies have sought to improve exploration techniques by integrating gravity and gravity gradient data using inversion methods. In this study, we developed a new method to integrate gravity and gravity gradient data in a 3D density inversion using the non-linear conjugate gradient (NLCG) method and the minimum gradient support (MGS) functional to regularize the 3D inverse problem and to obtain a clear and accurate image of the anomalous body. The NLCG algorithm, which is suitable for solving large-scale nonlinear optimization problems and requires no memory storage, was compared to the Broyden-Fletcher-Goldfarb-Shanno (BFGS) quasi-Newton algorithm and the results indicated that the convergence rate of NLCG is slower, but that the storage requirement and computation time is lower. To counteract the decay in kernel function, we introduced a depth weighting function for anomalous bodies at the same depth, with information about anomalous body depth obtained from well log and seismic exploration data. For anomalous bodies at different depths, we introduced a spatial gradient weighting function to incorporate additional information obtained in the inversion. We concluded that the spatial gradient weighting function enhanced the spatial resolution of the recovered model. Furthermore, our results showed that including multiple components for inversion increased the resolution of the recovered model. We validated our model by applying our inversion method to survey data from Vinton salt dome, Louisiana, USA. The results showed good agreement with known geologic information; thus confirming the accuracy of this approach.
Rosas-Carbajal, Marina; Linde, Nicolas; Peacock, Jared R.; Zyserman, F. I.; Kalscheuer, Thomas; Thiel, Stephan
2015-01-01
Surface-based monitoring of mass transfer caused by injections and extractions in deep boreholes is crucial to maximize oil, gas and geothermal production. Inductive electromagnetic methods, such as magnetotellurics, are appealing for these applications due to their large penetration depths and sensitivity to changes in fluid conductivity and fracture connectivity. In this work, we propose a 3-D Markov chain Monte Carlo inversion of time-lapse magnetotelluric data to image mass transfer following a saline fluid injection. The inversion estimates the posterior probability density function of the resulting plume, and thereby quantifies model uncertainty. To decrease computation times, we base the parametrization on a reduced Legendre moment decomposition of the plume. A synthetic test shows that our methodology is effective when the electrical resistivity structure prior to the injection is well known. The centre of mass and spread of the plume are well retrieved.We then apply our inversion strategy to an injection experiment in an enhanced geothermal system at Paralana, South Australia, and compare it to a 3-D deterministic time-lapse inversion. The latter retrieves resistivity changes that are more shallow than the actual injection interval, whereas the probabilistic inversion retrieves plumes that are located at the correct depths and oriented in a preferential north-south direction. To explain the time-lapse data, the inversion requires unrealistically large resistivity changes with respect to the base model. We suggest that this is partly explained by unaccounted subsurface heterogeneities in the base model from which time-lapse changes are inferred.
3D Inversion of complex resistivity data: Case study on Mineral Exploration Site.
NASA Astrophysics Data System (ADS)
Son, Jeong-Sul; Kim, Jung-ho; Park, Sam-gyu; Park, My-Kyung
2016-04-01
Complex resistivity (CR) method is a frequency domain induced polarization (IP) method. It is also known as Spectral IP (SIP) method, if wider frequencies are used in data acquisition and interpretation. Although it takes more times than conventional time domain IP method, its data quality is more stable because its data acquisition which measures amplitude and phase is done when the source current is being injected. Our research group has been studying the modeling and inversion algorithms of complex resistivity (CR) method since several years ago and recently applied developed algorithms to various real field application. Due to tough terrain in our country, Profile survey and 2D interpretation were generally used. But to get more precise interpretation, three dimensional modeling and inversion algorithm is required. We developed three dimensional inversion algorithm for this purpose. In the inversion, we adopt the method of adaptive lagraingian multiplier which is automatically set based on the size of error misfit and model regularization norm. It was applied on the real data acquired for mineral exploration sites. CR data was acquired with the Zeta system, manufactured by Zonge Co. In the inversion, only the lower frequency data is used considering its quality and developed 3D inversion algorithm was applied to the acquired data set. Its results were compared to those of time domain IP data conducted at the same site. Resistivity image sections of CR and conventional resistivity method were almost identical. Phase anomalies were well matched with chargeability anomalies and the mining history of the test site. Each anomalies were well discriminated in 3D interpretation than those of 2D. From those experiments, we know that CR method was very effective for the mineral exploration.
NASA Astrophysics Data System (ADS)
Corall, Silke; Haraszti, Tamas; Bartoschik, Tanja; Spatz, Joachim Pius; Ludwig, Thomas; Cavalcanti-Adam, Elisabetta Ada
2014-03-01
Cell migration is a crucial event for physiological processes, such as embryonic development and wound healing, as well as for pathological processes, such as cancer dissemination and metastasis formation. Cancer cell migration is a result of the concerted action of matrix metalloproteinases (MMPs), expressed by cancer cells to degrade the surrounding matrix, and integrins, the transmembrane receptors responsible for cell binding to matrix proteins. While it is known that cell-microenvironment interactions are essential for migration, the role of the physical state of such interactions remains still unclear. In this study we investigated human fibrosarcoma cell migration in two-dimensional (2D) and three-dimensional (3D) fibronectin (FN) microenvironments. By using antibody blocking approach and cell-binding site mutation, we determined that -integrin is the main mediator of fibrosarcoma cell migration in 2D FN, whereas in 3D fibrillar FN, the binding of - and -integrins is not necessary for cell movement in the fibrillar network. Furthermore, while the general inhibition of MMPs with GM6001 has no effect on cell migration in both 2D and 3D FN matrices, we observed opposing effect after targeted silencing of a membrane-bound MMP, namely MT1-MMP. In 2D fibronectin, silencing of MT1-MMP results in decreased migration speed and loss of directionality, whereas in 3D FN matrices, cell migration speed is increased and integrin-mediated signaling for actin dynamics is promoted. Our results suggest that the fibrillar nature of the matrix governs the migratory behavior of fibrosarcoma cells. Therefore, to hinder migration and dissemination of diseased cells, matrix molecules should be directly targeted, rather than specific subtypes of receptors at the cell membrane.
Earthquake source tensor inversion with the gCAP method and 3D Green's functions
NASA Astrophysics Data System (ADS)
Zheng, J.; Ben-Zion, Y.; Zhu, L.; Ross, Z.
2013-12-01
We develop and apply a method to invert earthquake seismograms for source properties using a general tensor representation and 3D Green's functions. The method employs (i) a general representation of earthquake potency/moment tensors with double couple (DC), compensated linear vector dipole (CLVD), and isotropic (ISO) components, and (ii) a corresponding generalized CAP (gCap) scheme where the continuous wave trains are broken into Pnl and surface waves (Zhu & Ben-Zion, 2013). For comparison, we also use the waveform inversion method of Zheng & Chen (2012) and Ammon et al. (1998). Sets of 3D Green's functions are calculated on a grid of 1 km3 using the 3-D community velocity model CVM-4 (Kohler et al. 2003). A bootstrap technique is adopted to establish robustness of the inversion results using the gCap method (Ross & Ben-Zion, 2013). Synthetic tests with 1-D and 3-D waveform calculations show that the source tensor inversion procedure is reasonably reliable and robust. As initial application, the method is used to investigate source properties of the March 11, 2013, Mw=4.7 earthquake on the San Jacinto fault using recordings of ~45 stations up to ~0.2Hz. Both the best fitting and most probable solutions include ISO component of ~1% and CLVD component of ~0%. The obtained ISO component, while small, is found to be a non-negligible positive value that can have significant implications for the physics of the failure process. Work on using higher frequency data for this and other earthquakes is in progress.
Joint earthquake source inversions using seismo-geodesy and 3-D earth models
NASA Astrophysics Data System (ADS)
Weston, J.; Ferreira, A. M. G.; Funning, G. J.
2014-08-01
A joint earthquake source inversion technique is presented that uses InSAR and long-period teleseismic data, and, for the first time, takes 3-D Earth structure into account when modelling seismic surface and body waves. Ten average source parameters (Moment, latitude, longitude, depth, strike, dip, rake, length, width and slip) are estimated; hence, the technique is potentially useful for rapid source inversions of moderate magnitude earthquakes using multiple data sets. Unwrapped interferograms and long-period seismic data are jointly inverted for the location, fault geometry and seismic moment, using a hybrid downhill Powell-Monte Carlo algorithm. While the InSAR data are modelled assuming a rectangular dislocation in a homogeneous half-space, seismic data are modelled using the spectral element method for a 3-D earth model. The effect of noise and lateral heterogeneity on the inversions is investigated by carrying out realistic synthetic tests for various earthquakes with different faulting mechanisms and magnitude (Mw 6.0-6.6). Synthetic tests highlight the improvement in the constraint of fault geometry (strike, dip and rake) and moment when InSAR and seismic data are combined. Tests comparing the effect of using a 1-D or 3-D earth model show that long-period surface waves are more sensitive than long-period body waves to the change in earth model. Incorrect source parameters, particularly incorrect fault dip angles, can compensate for systematic errors in the assumed Earth structure, leading to an acceptable data fit despite large discrepancies in source parameters. Three real earthquakes are also investigated: Eureka Valley, California (1993 May 17, Mw 6.0), Aiquile, Bolivia (1998 February 22, Mw 6.6) and Zarand, Iran (2005 May 22, Mw 6.5). These events are located in different tectonic environments and show large discrepancies between InSAR and seismically determined source models. Despite the 40-50 km discrepancies in location between previous geodetic and
Advanced Multivariate Inversion Techniques for High Resolution 3D Geophysical Modeling (Invited)
NASA Astrophysics Data System (ADS)
Maceira, M.; Zhang, H.; Rowe, C. A.
2009-12-01
We focus on the development and application of advanced multivariate inversion techniques to generate a realistic, comprehensive, and high-resolution 3D model of the seismic structure of the crust and upper mantle that satisfies several independent geophysical datasets. Building on previous efforts of joint invesion using surface wave dispersion measurements, gravity data, and receiver functions, we have added a fourth dataset, seismic body wave P and S travel times, to the simultaneous joint inversion method. We present a 3D seismic velocity model of the crust and upper mantle of northwest China resulting from the simultaneous, joint inversion of these four data types. Surface wave dispersion measurements are primarily sensitive to seismic shear-wave velocities, but at shallow depths it is difficult to obtain high-resolution velocities and to constrain the structure due to the depth-averaging of the more easily-modeled, longer-period surface waves. Gravity inversions have the greatest resolving power at shallow depths, and they provide constraints on rock density variations. Moreover, while surface wave dispersion measurements are primarily sensitive to vertical shear-wave velocity averages, body wave receiver functions are sensitive to shear-wave velocity contrasts and vertical travel-times. Addition of the fourth dataset, consisting of seismic travel-time data, helps to constrain the shear wave velocities both vertically and horizontally in the model cells crossed by the ray paths. Incorporation of both P and S body wave travel times allows us to invert for both P and S velocity structure, capitalizing on empirical relationships between both wave types’ seismic velocities with rock densities, thus eliminating the need for ad hoc assumptions regarding the Poisson ratios. Our new tomography algorithm is a modification of the Maceira and Ammon joint inversion code, in combination with the Zhang and Thurber TomoDD (double-difference tomography) program.
Finite Element Based Anisotropic 3D Inversion of Marine CSEM Data
NASA Astrophysics Data System (ADS)
Chung, Y.; Byun, J.
2015-12-01
In order to interpret three-dimensional (3D) marine controlled-source electromagnetic (MCSEM) data, it is critical to accurately determine electrical anisotropy because ignoring anisotropy can produce misleading artifacts. In this study, we present an inversion method for 3D subsurface imaging in media with an inhomogeneous and anisotropic conductivity distribution. Direct solvers are incorporated both in the forward and inverse problems, For the forward problem, the vector Helmholtz equation for the secondary electric field is discretized on a hexahedral mesh using edge finite elements, then a direct sparse-matrix solver is chosen to effectively reuse its factorization both in the survey simulation and Jacobian computation. The inversion method is formulated as a functional optimization with an objective functional containing terms measuring data misfit and model structure by means of smoothness and anisotropy. These measures are efficiently incorporated through the use of an iteratively reweighted least-squares scheme. The objective functional is minimized by a Gauss-Newton approach using a direct dense-matrix solver. We demonstrate the accuracy and applicability of the algorithm by testing it on synthetic data sets.
3D inversion of time-lapse CSEM data for reservoir monitoring
NASA Astrophysics Data System (ADS)
Black, N.; Wilson, G. A.; Zhdanov, M. S.
2010-12-01
Effective reservoir monitoring requires time-lapse reservoir information throughout the interwell volume. The ability to understand and control reservoir behavior over the course of production allows for optimization of reservoir performance and production strategies. Good monitoring information makes it possible to improve the timing and location of new drilling (for both production and injection wells), to recognize flow paths, and to map oil that has been bypassed. Recent studies have inferred the feasibility of time-lapse marine controlled-source electromagnetic (CSEM) methods for the monitoring of offshore oil and gas fields. However, quantitative interpretations to ascertain what reservoir information may be recovered have not been performed. The time-lapse CSEM inverse problem can be highly constrained since the geometry of the reservoir is established prior from high resolution seismic surveys, rock and fluid properties are measured from well logs, and multiple history matched production scenarios are contained in dynamic reservoir models. We present a 3D inversion study of synthetic time-lapse CSEM data modeled from dynamic reservoir simulations. We demonstrate that even with few constraints on the model, the hydrocarbon-water front can be recovered from 3D inversion.
A hybrid method for inversion of 3D DC resistivity logging measurements.
Gajda-Zagórska, Ewa; Schaefer, Robert; Smołka, Maciej; Paszyński, Maciej; Pardo, David
This paper focuses on the application of hp hierarchic genetic strategy (hp-HGS) for solution of a challenging problem, the inversion of 3D direct current (DC) resistivity logging measurements. The problem under consideration has been formulated as the global optimization one, for which the objective function (misfit between computed and reference data) exhibits multiple minima. In this paper, we consider the extension of the hp-HGS strategy, namely we couple the hp-HGS algorithm with a gradient based optimization method for a local search. Forward simulations are performed with a self-adaptive hp finite element method, hp-FEM. The computational cost of misfit evaluation by hp-FEM depends strongly on the assumed accuracy. This accuracy is adapted to the tree of populations generated by the hp-HGS algorithm, which makes the global phase significantly cheaper. Moreover, tree structure of demes as well as branch reduction and conditional sprouting mechanism reduces the number of expensive local searches up to the number of minima to be recognized. The common (direct and inverse) accuracy control, crucial for the hp-HGS efficiency, has been motivated by precise mathematical considerations. Numerical results demonstrate the suitability of the proposed method for the inversion of 3D DC resistivity logging measurements.
Solving Dirac equations on a 3D lattice with inverse Hamiltonian and spectral methods
NASA Astrophysics Data System (ADS)
Ren, Z. X.; Zhang, S. Q.; Meng, J.
2017-02-01
A new method to solve the Dirac equation on a 3D lattice is proposed, in which the variational collapse problem is avoided by the inverse Hamiltonian method and the fermion doubling problem is avoided by performing spatial derivatives in momentum space with the help of the discrete Fourier transform, i.e., the spectral method. This method is demonstrated in solving the Dirac equation for a given spherical potential in a 3D lattice space. In comparison with the results obtained by the shooting method, the differences in single-particle energy are smaller than 10-4 MeV, and the densities are almost identical, which demonstrates the high accuracy of the present method. The results obtained by applying this method without any modification to solve the Dirac equations for an axial-deformed, nonaxial-deformed, and octupole-deformed potential are provided and discussed.
High-resolution imaging of crustal melts using 3D full-waveform seismic inversion
NASA Astrophysics Data System (ADS)
Warner, M.; Morgan, J. V.
2013-12-01
A newly practical seismic imaging technique, 3D full-waveform inversion (FWI), now has the ability to image zones of melt and melt pathways throughout the crust with a better resolution than any other geophysical method. 3D FWI has recently changed practice within the petroleum industry where it is used to obtain high-resolution high-fidelity models of physical properties in the sub-surface that are both interpreted directly and used to improve the migration of deeper reflections. This technology has been spectacularly successful in improving the imaging of reservoirs beneath shallow heterogeneities produced by, for example, gas clouds, buried fluvial channels, carbonate reefs and salt bodies. During FWI, the sub-surface model is recovered principally by using the low-frequency transmitted, refracted portion of the wavefield which is most sensitive to the macro-velocity structure. In the petroleum industry, these inversions are now routinely performed using long-offset surface-streamer and ocean-bottom data to maximum source-receiver offsets of about 15 km, leading to a maximum penetration depth of around 5 km. Using longer offsets, it is possible to extend this technology to image deeper crustal targets. Localised zones of partial melt produce large changes in p-wave and s-wave properties that are restricted in their spatial extent, and that therefore form ideal targets for 3D FWI. We have performed a suite of tests to explore the use of 3D FWI in imaging melt distribution beneath the active volcano of Montserrat. We built a model of the subsurface using a 3D travel-time tomographic model obtained from the SEA CALIPSO experiment. We added two magma chambers in accordance with a model obtained using surface-elevation changes and geochemical data. We used a wide-angle, wide-azimuth acquisition geometry to generate a fully-elastic synthetic seismic dataset, added noise, and inverted the windowed transmitted arrivals only. We used an elastic code for the forward
NASA Astrophysics Data System (ADS)
Yang, Dikun; Oldenburg, Douglas W.; Haber, Eldad
2014-03-01
Airborne electromagnetic (AEM) methods are highly efficient tools for assessing the Earth's conductivity structures in a large area at low cost. However, the configuration of AEM measurements, which typically have widely distributed transmitter-receiver pairs, makes the rigorous modelling and interpretation extremely time-consuming in 3-D. Excessive overcomputing can occur when working on a large mesh covering the entire survey area and inverting all soundings in the data set. We propose two improvements. The first is to use a locally optimized mesh for each AEM sounding for the forward modelling and calculation of sensitivity. This dedicated local mesh is small with fine cells near the sounding location and coarse cells far away in accordance with EM diffusion and the geometric decay of the signals. Once the forward problem is solved on the local meshes, the sensitivity for the inversion on the global mesh is available through quick interpolation. Using local meshes for AEM forward modelling avoids unnecessary computing on fine cells on a global mesh that are far away from the sounding location. Since local meshes are highly independent, the forward modelling can be efficiently parallelized over an array of processors. The second improvement is random and dynamic down-sampling of the soundings. Each inversion iteration only uses a random subset of the soundings, and the subset is reselected for every iteration. The number of soundings in the random subset, determined by an adaptive algorithm, is tied to the degree of model regularization. This minimizes the overcomputing caused by working with redundant soundings. Our methods are compared against conventional methods and tested with a synthetic example. We also invert a field data set that was previously considered to be too large to be practically inverted in 3-D. These examples show that our methodology can dramatically reduce the processing time of 3-D inversion to a practical level without losing resolution
3D linear inversion of magnetic susceptibility data acquired by frequency domain EMI
NASA Astrophysics Data System (ADS)
Thiesson, J.; Tabbagh, A.; Simon, F.-X.; Dabas, M.
2017-01-01
Low induction number EMI instruments are able to simultaneously measure a soil's apparent magnetic susceptibility and electrical conductivity. This family of dual measurement instruments is highly useful for the analysis of soils and archeological sites. However, the electromagnetic properties of soils are found to vary over considerably different ranges: whereas their electrical conductivity varies from ≤ 0.1 to ≥ 100 mS/m, their relative magnetic permeability remains within a very small range, between 1.0001 and 1.01 SI. Consequently, although apparent conductivity measurements need to be inverted using non-linear processes, the variations of the apparent magnetic susceptibility can be approximated through the use of linear processes, as in the case of the magnetic prospection technique. Our proposed 3D inversion algorithm starts from apparent susceptibility data sets, acquired using different instruments over a given area. A reference vertical profile is defined by considering the mode of the vertical distributions of both the electrical resistivity and of the magnetic susceptibility. At each point of the mapped area, the reference vertical profile response is subtracted to obtain the apparent susceptibility variation dataset. A 2D horizontal Fourier transform is applied to these variation datasets and to the dipole (impulse) response of each instrument, a (vertical) 1D inversion is performed at each point in the spectral domain, and finally the resulting dataset is inverse transformed to restore the apparent 3D susceptibility variations. It has been shown that when applied to synthetic results, this method is able to correct the apparent deformations of a buried object resulting from the geometry of the instrument, and to restore reliable quantitative susceptibility contrasts. It also allows the thin layer solution, similar to that used in magnetic prospection, to be implemented. When applied to field data it initially delivers a level of contrast
3D Self-Potential Inversion for Monitoring DNAPL Contaminant Distributions
NASA Astrophysics Data System (ADS)
Minsley, B. J.; Sogade, J.; Vichabian, Y.; Morgan, F. D.
2005-05-01
Self-potential (SP) data are collected over an area known to be contaminated with Dense Non-Aqueous Phase Liquids (DNAPLs) at the Savannah River Site in South Carolina. The field experiment consists of approximately 100 SP measurements on a surface grid and in four boreholes, and is repeated after one year. DNAPLs are known to undergo redox reactions during their degradation in the environment, which is often biologically mediated. Self-potential geophysics is employed in this study because of its sensitivity to the in-situ biochemical processes that degrade the contaminants. These reactions provide an electrochemical source that is manifested as an SP signature at the measurement locations remote from the contaminated areas. 3D inversion of the SP data is therefore needed to spatially locate the distribution of sources, which is related to contaminant presence. The inversion incorporates the 3D resistivity structure collected at the same site, and is better constrained in depth by using borehole data and regularization. Ground truth information taken after the first field experiment provides concentration data with depth for several DNAPL species in five boreholes. There is a good correlation between the ground truth data and SP source inversion, though this comparison is limited by several factors: the difference in resolution of the ground truth and inverted data, and the dependence of the redox processes on other constituents that were not measured during the ground truthing, such as oxygen content or microbial presence. Inversion of the second year's dataset provides information on the changes in the contaminant distribution, either due to natural degradation or ongoing remediation.
NASA Astrophysics Data System (ADS)
Heincke, B.; Moorkamp, M.; Jegen, M.; Hobbs, R. W.
2012-12-01
Imaging of sub-basalt sediments with reflection seismic techniques is limited due to absorption, scattering and transmission effects and the presence of peg-leg multiples. Although many of the difficulties facing conventional seismic profiles can be overcome by recording long offset data resolution of sub-basalt sediments in seismic sections is typically still largely restricted. Therefore multi-parametric approaches in general and joint inversion strategies in particular (e.g. Colombo et al., 2008, Jordan et al., 2012) are considered as alternative to gain additional information from sub-basalt structures. Here, we combine in a 3-D joint inversion first-arrival time tomography, FTG gravity and MT data to identify the base basalt and resolve potential sediments underneath. For sub-basalt exploration the three methods complement each other such that the null space is reduced and significantly better resolved models can be obtained than would be possible by the individual methods: The seismic data gives a robust model for the supra-basalt sediments whilst the gravity field is dominated by the high density basalt and basement features. The MT on the other hand is sensitive to the conductivity in both the supra- and sub-basalt sediments. We will present preliminary individual and joint inversion result for a FTG, seismic and MT data set located in the Faroe-Shetland basin. Because the investigated area is rather large (~75 x 40 km) and the individual data sets are relatively huge, we use a joint inversion framework (see Moorkamp et al., 2011) which is designed to handle large amount of data/model parameters. This program has moreover the options to link the individual parameter models either petrophysically using fixed parameter relationships or structurally using the cross-gradient approach. The seismic data set consists of a pattern of 8 intersecting wide-angle seismic profiles with maximum offsets of up to ~24 km. The 3-D gravity data set (size :~ 30 x 30 km) is
NASA Astrophysics Data System (ADS)
Kordy, M.; Wannamaker, P.; Maris, V.; Cherkaev, E.; Hill, G.
2016-01-01
We have developed an algorithm, which we call HexMT, for 3-D simulation and inversion of magnetotelluric (MT) responses using deformable hexahedral finite elements that permit incorporation of topography. Direct solvers parallelized on symmetric multiprocessor (SMP), single-chassis workstations with large RAM are used throughout, including the forward solution, parameter Jacobians and model parameter update. In Part I, the forward simulator and Jacobian calculations are presented. We use first-order edge elements to represent the secondary electric field (E), yielding accuracy O(h) for E and its curl (magnetic field). For very low frequencies or small material admittivities, the E-field requires divergence correction. With the help of Hodge decomposition, the correction may be applied in one step after the forward solution is calculated. This allows accurate E-field solutions in dielectric air. The system matrix factorization and source vector solutions are computed using the MKL PARDISO library, which shows good scalability through 24 processor cores. The factorized matrix is used to calculate the forward response as well as the Jacobians of electromagnetic (EM) field and MT responses using the reciprocity theorem. Comparison with other codes demonstrates accuracy of our forward calculations. We consider a popular conductive/resistive double brick structure, several synthetic topographic models and the natural topography of Mount Erebus in Antarctica. In particular, the ability of finite elements to represent smooth topographic slopes permits accurate simulation of refraction of EM waves normal to the slopes at high frequencies. Run-time tests of the parallelized algorithm indicate that for meshes as large as 176 × 176 × 70 elements, MT forward responses and Jacobians can be calculated in ˜1.5 hr per frequency. Together with an efficient inversion parameter step described in Part II, MT inversion problems of 200-300 stations are computable with total run times
3D CSEM data inversion using Newton and Halley class methods
NASA Astrophysics Data System (ADS)
Amaya, M.; Hansen, K. R.; Morten, J. P.
2016-05-01
For the first time in 3D controlled source electromagnetic data inversion, we explore the use of the Newton and the Halley optimization methods, which may show their potential when the cost function has a complex topology. The inversion is formulated as a constrained nonlinear least-squares problem which is solved by iterative optimization. These methods require the derivatives up to second order of the residuals with respect to model parameters. We show how Green's functions determine the high-order derivatives, and develop a diagrammatical representation of the residual derivatives. The Green's functions are efficiently calculated on-the-fly, making use of a finite-difference frequency-domain forward modelling code based on a multi-frontal sparse direct solver. This allow us to build the second-order derivatives of the residuals keeping the memory cost in the same order as in a Gauss-Newton (GN) scheme. Model updates are computed with a trust-region based conjugate-gradient solver which does not require the computation of a stabilizer. We present inversion results for a synthetic survey and compare the GN, Newton, and super-Halley optimization schemes, and consider two different approaches to set the initial trust-region radius. Our analysis shows that the Newton and super-Halley schemes, using the same regularization configuration, add significant information to the inversion so that the convergence is reached by different paths. In our simple resistivity model examples, the convergence speed of the Newton and the super-Halley schemes are either similar or slightly superior with respect to the convergence speed of the GN scheme, close to the minimum of the cost function. Due to the current noise levels and other measurement inaccuracies in geophysical investigations, this advantageous behaviour is at present of low consequence, but may, with the further improvement of geophysical data acquisition, be an argument for more accurate higher-order methods like those
NASA Astrophysics Data System (ADS)
Meqbel, N. M.; Egbert, G. D.; Kelbert, A.
2011-12-01
Long period (10-20,000 s) magnetotelluric (MT) data are being acquired in a series of temporary arrays deployed across the continental United States through the EMScope component of EarthScope. MT deployments in 2006-2011 have acquired data at 325 sites on an approximately regular grid, with the same nominal spacing as the USArray broadband seismic transportable array (~70 km). The MT sites span a rectangular area from NW Washington to NW Colorado. Here we present results of a 3-D inversion of the full data set. A number of conductive and resistive features appear consistently in the crust and upper mantle in essentially all of a large suite of 3-D inverse solutions. Extensive areas of high conductivity are found in the lower crust (up to a depth of ~ 40 km) beneath the Basin & Range in southeastern Oregon, as imaged by Patro and Egbert (2008). In our new model, this feature extends further to the south and to the east, where it merges with somewhat deeper (uppermost mantle) conductivities beneath the Yellowstone-Snake River Plain. This deeper feature, which extends from Yellowstone to the SW into northeastern Nevada, coincides with the track of the Yellowstone hotspot discussed e.g., in Smith et. al. (2008). The lower crust and the uppermost mantle in the northeastern part of the domain, covering the area from eastern Washington to Montana and continuing south to Wyoming, is generally resistive, with a few localized exceptions. This resistive zone coincides with high velocities discussed and interpreted, e.g., by Yang et. al. (2008) as thick, stable Proterozoic lithosphere. A number of large-scale anomalous features also appear consistently in the upper mantle, at depths of ~ 50 km to 300 km. Most striking is a zone of high resistivity on the western edge of the domain, beneath western Oregon, Washington and northern California in the area occupied by oceanic lithosphere of the Juan de Fuca Plate, which has subducted beneath the relatively more conductive
Levander, Alan Richard; Zelt, Colin A.
2015-03-17
The work plan for this project was to develop and apply advanced seismic reflection and wide-angle processing and inversion techniques to high resolution seismic data for the shallow subsurface to seismically characterize the shallow subsurface at hazardous waste sites as an aid to containment and cleanup activities. We proposed to continue work on seismic data that we had already acquired under a previous DoE grant, as well as to acquire additional new datasets for analysis. The project successfully developed and/or implemented the use of 3D reflection seismology algorithms, waveform tomography and finite-frequency tomography using compressional and shear waves for high resolution characterization of the shallow subsurface at two waste sites. These two sites have markedly different near-surface structures, groundwater flow patterns, and hazardous waste problems. This is documented in the list of refereed documents, conference proceedings, and Rice graduate theses, listed below.
The 3D inversion of airborne gamma-ray spectrometric data
NASA Astrophysics Data System (ADS)
Minty, Brian; Brodie, Ross
2016-07-01
We present a new method for the inversion of airborne gamma-ray spectrometric line data to a regular grid of radioelement concentration estimates on the ground. The method incorporates the height of the aircraft, the 3D terrain within the field of view of the spectrometer, the directional sensitivity of rectangular detectors, and a source model comprising vertical rectangular prisms with the same horizontal dimensions as the required grid cell size. The top of each prism is a plane surface derived from a best-fit plane to the digital elevation model of the earth's surface within each grid cell area. The method is a significant improvement on current methods, and gives superior interpolation between flight lines. It also eliminates terrain effects that would normally remain in the data after the conventional processing of these data assuming a flat-earth model.
NASA Astrophysics Data System (ADS)
Yang, Pengliang; Brossier, Romain; Métivier, Ludovic; Virieux, Jean
2016-10-01
In this paper, we study 3-D multiparameter full waveform inversion (FWI) in viscoelastic media based on the generalized Maxwell/Zener body including arbitrary number of attenuation mechanisms. We present a frequency-domain energy analysis to establish the stability condition of a full anisotropic viscoelastic system, according to zero-valued boundary condition and the elastic-viscoelastic correspondence principle: the real-valued stiffness matrix becomes a complex-valued one in Fourier domain when seismic attenuation is taken into account. We develop a least-squares optimization approach to linearly relate the quality factor with the anelastic coefficients by estimating a set of constants which are independent of the spatial coordinates, which supplies an explicit incorporation of the parameter Q in the general viscoelastic wave equation. By introducing the Lagrangian multipliers into the matrix expression of the wave equation with implicit time integration, we build a systematic formulation of multiparameter FWI for full anisotropic viscoelastic wave equation, while the equivalent form of the state and adjoint equation with explicit time integration is available to be resolved efficiently. In particular, this formulation lays the foundation for the inversion of the parameter Q in the time domain with full anisotropic viscoelastic properties. In the 3-D isotropic viscoelastic settings, the anelastic coefficients and the quality factors using bulk and shear moduli parametrization can be related to the counterparts using P and S velocity. Gradients with respect to any other parameter of interest can be found by chain rule. Pioneering numerical validations as well as the real applications of this most generic framework will be carried out to disclose the potential of viscoelastic FWI when adequate high-performance computing resources and the field data are available.
Inverse modeling of InSAR and ground leveling data for 3D volumetric strain distribution
NASA Astrophysics Data System (ADS)
Gallardo, L. A.; Glowacka, E.; Sarychikhina, O.
2015-12-01
Wide availability of modern Interferometric Synthetic aperture Radar (InSAR) data have made possible the extensive observation of differential surface displacements and are becoming an efficient tool for the detailed monitoring of terrain subsidence associated to reservoir dynamics, volcanic deformation and active tectonism. Unfortunately, this increasing popularity has not been matched by the availability of automated codes to estimate underground deformation, since many of them still rely on trial-error subsurface model building strategies. We posit that an efficient algorithm for the volumetric modeling of differential surface displacements should match the availability of current leveling and InSAR data and have developed an algorithm for the joint inversion of ground leveling and dInSAR data in 3D. We assume the ground displacements are originated by a stress free-volume strain distribution in a homogeneous elastic media and determined the displacement field associated to an ensemble of rectangular prisms. This formulation is then used to develop a 3D conjugate gradient inversion code that searches for the three-dimensional distribution of the volumetric strains that predict InSAR and leveling surface displacements simultaneously. The algorithm is regularized applying discontinuos first and zero order Thikonov constraints. For efficiency, the resulting computational code takes advantage of the resulting convolution integral associated to the deformation field and some basic tools for multithreading parallelization. We extensively test our algorithm on leveling and InSAR test and field data of the Northwest of Mexico and compare to some feasible geological scenarios of underground deformation.
Review on applications of 3D inverse design method for pump
NASA Astrophysics Data System (ADS)
Yin, Junlian; Wang, Dezhong
2014-05-01
The 3D inverse design method, which methodology is far superior to the conventional design method that based on geometrical description, is gradually applied in pump blade design. However, no complete description about the method is outlined. Also, there are no general rules available to set the two important input parameters, blade loading distribution and stacking condition. In this sense, the basic theory and the mechanism why the design method can suppress the formation of secondary flow are summarized. And also, several typical pump design cases with different specific speeds ranging from centrifugal pump to axial pump are surveyed. The results indicates that, for centrifugal pump and mixed pump or turbine, the ratio of blade loading on the hub to that on the shroud is more than unit in the fore part of the blade, whereas in the aft part, the ratio is decreased to satisfy the same wrap angle for hub and shroud. And the choice of blade loading type depends on the balancing of efficiency and cavitation. If the cavitation is more weighted, the better choice is aft-loaded, otherwise, the fore-loaded or mid-loaded is preferable to improve the efficiency. The stacking condition, which is an auxiliary to suppress the secondary flow, can have great effect on the jet-wake outflow and the operation range for pump. Ultimately, how to link the design method to modern optimization techniques is illustrated. With the know-how design methodology and the know-how systematic optimization approach, the application of optimization design is promising for engineering. This paper summarizes the 3D inverse design method systematically.
3D elastic full waveform inversion: case study from a land seismic survey
NASA Astrophysics Data System (ADS)
Kormann, Jean; Marti, David; Rodriguez, Juan-Esteban; Marzan, Ignacio; Ferrer, Miguel; Gutierrez, Natalia; Farres, Albert; Hanzich, Mauricio; de la Puente, Josep; Carbonell, Ramon
2016-04-01
Full Waveform Inversion (FWI) is one of the most advanced processing methods that is recently reaching a mature state after years of solving theoretical and technical issues such as the non-uniqueness of the solution and harnessing the huge computational power required by realistic scenarios. BSIT (Barcelona Subsurface Imaging Tools, www.bsc.es/bsit) includes a FWI algorithm that can tackle with very complex problems involving large datasets. We present here the application of this system to a 3D dataset acquired to constrain the shallow subsurface. This is where the wavefield is the most complicated, because most of the wavefield conversions takes place in the shallow region and also because the media is much more laterally heterogeneous. With this in mind, at least isotropic elastic approximation would be suitable as kernel engine for FWI. The current study explores the possibilities to apply elastic isotropic FWI using only the vertical component of the recorded seismograms. The survey covers an area of 500×500 m2, and consists in a receivers grid of 10 m×20 m combined with a 250 kg accelerated weight-drop as source on a displaced grid of 20 m×20 m. One of the main challenges in this case study is the costly 3D modeling that includes topography and substantial free surface effects. FWI is applied to a data subset (shooting lines 4 to 12), and is performed for 3 frequencies ranging from 15 to 25 Hz. The starting models are obtained from travel-time tomography and the all computation is run on 75 nodes of Mare Nostrum supercomputer during 3 days. The resulting models provide a higher resolution of the subsurface structures, and show a good correlation with the available borehole measurements. FWI allows to extend in a reliable way this 1D knowledge (borehole) to 3D.
Inverse cascades sustained by the transfer rate of angular momentum in a 3D turbulent flow
NASA Astrophysics Data System (ADS)
Burguete, Javier; Lopez-Caballero, Miguel
2013-11-01
The existence of energy cascades as signatures of conserved magnitudes is one of the universal characteristics of turbulent flows. In this work we present the evidence of an inverse cascade in a fully developed 3D experimental turbulent flow where the conserved magnitude is the angular momentum. We analyze the behavior of a fluid in a closed cavity where two inhomogeneous and strongly turbulent flows collide in a thin region. The experimental volume is a closed cylinder (diameter of 20 cm) where two impellers rotate in opposite directions. A key characteristic of this setup the high stability of the propellers (the instantaneous fluctuations are below 0 . 1 %). We have performed PIV and LDA measurements of the velocity fields. Typical characteristics of the turbulent flow in this setup are: turbulence intensity 50 % , the Reλ = 900 , the Taylor microscale λT = 1 . 8 mm and the integral scale LI = 15 mm. The analysis of the data series reveal that below the injection scales an inverse cascade can be identified (-1/3 in time, -7/3 in space) that can be explained as the transfer of angular momentum between the diferent fluid layers. A. de la Torre, J. Burguete, Phys Rev Lett 99 (2007) 054101. M. Lopez-Caballero, J. Burguete, Phys Rev Lett 110 (2013) 124501.
KOALA: 3-D shape of asteroids from multi-data inversion
NASA Astrophysics Data System (ADS)
Carry, B.; Kaasalainen, M.; Merline, W. J.; Drummond, J. D.; Durech, J.; Berthier, J.; Conrad, A.
2011-10-01
We describe our on-going observing program to determine the physical properties of asteroids from groundbased facilities. We combine disk-resolved images from adaptive optics, optical lightcurves, and stellar occultations to put tighter constraints on the spin, 3-D shape, and size of asteroids. We will discuss the relevance of the determination of physical properties to help understand the asteroid population (e.g., density, composition, and non-gravitational forces). We will then briefly describe our multi-data inversion algorithm KOALA (Carry et al. 2010a, Kaasalainen 2011, see also Kaasalainen et al., same meeting), which allows the determination of certain physical properties of an asteroid from the combination of different techniques of observation. A comparison of results obtained with KOALA on asteroid (21) Lutetia, prior to the ESA Rosetta flyby, with the high spatial resolution images returned from that flyby, will then be presented, showing the high accuracy of KOALA inversion. Finally, we will describe our current development of the algorithm, and focus on examples of other asteroids currently being studied with KOALA.
NASA Astrophysics Data System (ADS)
Spicer, B.; Morris, B.; Ugalde, H.
2011-09-01
Hosted within the Pacquet Harbour Group (PHG) on the Baie Verte Peninsula of north-central Newfoundland, the Rambler rhyolite is a 487 Ma unit of felsic tuffs, flows and subvolcanic intrusive rocks. The PHG has been affected by multiple phases of deformation with the youngest D4 deformation event producing broad northeast plunging upright cross folds in the Rambler rhyolite. Fold culminations on the upper bounding surface of the rhyolite host Cu +/- Au volcanogenic massive sulfide deposits (e.g. Rambler and Ming mines). Geophysical inversions of recently acquired high resolution gravity and magnetic data have been implemented to determine the extent of the fold axis (dome) at depth. To direct the outcome of the inversion process towards a more geologically reasonable solution this study outlines a procedure which permits the inclusion of known geological and geophysical constraints into the input (reference) model for inversion using the MAG3D and GRAV3D algorithms provided by the University of British Columbia Geophysical Inversion Facility. Reference model constraints included surficial geological contacts as defined by aeromagnetic data, and subsurface distribution of physical property variations from a series of drill-hole logs. The output (computed) model images the surface of the rhyolite dome as dipping roughly 40° to the northeast as a series of voxels with density values ranging from 2.71 to 2.75 g/cm3. While previously published ore deposit models parallel this structure in the near surface, results from these inversions suggest deeper exploration may be favorable. Magnetic inversion modeling has not provided any insight into dome morphology however it outlines the distribution of gabbroic dykes surrounding the dome.
An optimal transport approach for seismic tomography: application to 3D full waveform inversion
NASA Astrophysics Data System (ADS)
Métivier, L.; Brossier, R.; Mérigot, Q.; Oudet, E.; Virieux, J.
2016-11-01
the L 2 distance, in 2D and 3D contexts.
Centroid Moment Tensor Inversion in a 3D heterogeneous Earth: Application to the Australasian region
NASA Astrophysics Data System (ADS)
Hejrani, B.; Tkalcic, H.; Fichtner, A.
2015-12-01
radially anisotropic structure: new insights into present and past states of the Australasian upper mantle. Earth Planet. Sci. Lett. 290, 270-280. Hingee, M., Tkalčić, H., Fichtner A., Sambridge, M., 2011. Moment tensor inversion using a 3-D structural model: Applications for the Australian region, Geophys. J. Int., 184(2), 949-964.
NASA Astrophysics Data System (ADS)
Rosenkjaer, G. K.; Heagy, L. J.; Cockett, R.; Kang, S.
2015-12-01
The MagnetoTelluric (MT) method is widely used in geophysical exploration, especially in geothermal exploration. The extensive frequency range allows mapping the resistivity structure to greater depths than typically achieved using controlled sources (on the order of 10 km's). This is advantageous in identifying structures of interest, such as those related to hydrothermal alterations in a geothermal system and the heat source of the system. To obtain information about subsurface resistivity structures from MT data, an inversion modelling is commonly performed. The goal of an inversion is to recover a resistivity model that explains the data (within a degree of accuracy), given sets of constraints on the model. An attention to detail, often through an iterative process, is required during the inversion. A workflow that is transparent and repeatable is important for addressing the complexities encountered in an inversion and make the iterative process more effective. In this presentation, we will show examples of a workflow for an inversion of MT data from the Krysuvik geothermal area in Iceland. The workflow is implemented as a Python script, where the SimPEG framework (http://simpeg.xyz) is used to implement an open-source code base to work with MT data . The SimPEG framework uses finite-volume discretizations for structured and semi-structured meshes and has built-in routines for optimization and regularization. By structuring and writing the simpegMT code in compliance with the modular, object-oriented framework of SimPEG, full compatibility is achieved and its functions can be leveraged directly. Further, the structure facilitates interoperability with other problem-types, streamlines testing and fosters cooperations between developers and users of the code. We will step through the workflow, where varying complexity leads to a recovery of 3D models of different detail. We show how we guide steps in the workflow based on the preceding results, iteratively
NASA Astrophysics Data System (ADS)
Bischoff, S. H.; Flesch, L. M.
2015-12-01
Piecing together the uplift and growth of the Tibetan Plateau requires a robust understanding of the present-day dynamics of the India-Eurasia collision zone. To aid in the understanding of mountain building and plateau growth, we developed a 3D finite element model of the Tibetan Plateau following Flesch and Bendick (2012). Our model is based on the vast collection of published geophysical data and employs COMSOL Multiphysics (www.comsol.com). We assume model material properties from the wide variety of published seismic and MT studies, incorporated with an updated, vertically averaged, effective viscosity distribution from Flesch et al. (2001). We test potential relationships between conductance/seismic velocity and strength (viscosity) by modeling strength difference contacts at imaged interfaces. We quantify fitness of candidate 3D viscosity functions by comparing solved model surface velocities to observed surface velocities inferred from GPS and Quaternary fault slip rates. The model geometry incorporates Earth curvature and extends eastward from 65° to 110°E, northward from 15° to 45°N, and vertically down to 100 km below sea level. The physics of deformation is governed by the Stokes equations describing incompressible Newtonian fluid flow. Boundary conditions consist of free slip across the bottom surface (representing the lithosphere-asthenosphere boundary) and moving edge walls constrained by a GPS-derived, continuous velocity field. Model results indicate a tradeoff between crust and mantle dominant strength. Best-fit models are achieved by a combination of strong crust/upper mantle with additional strain accommodation in localized weak zones.
High-resolution imaging and inversion of 3D GPR data for layered media
NASA Astrophysics Data System (ADS)
Slob, Evert
2013-04-01
Ground penetrating radar is increasingly being used to provide quantitative information of layered structures. For application in civil engineering these can be roads, highway pavements, airport runways, bridges, tunnels, or buildings. Monitoring is important for the management and safety of these structures. Standard imaging uses a modeled wavefield extrapolator to image the data and the quality of the image depends heavily on the quality of the modeled extrapolator. Usually, data inversion is implemented by minimizing a cost function involving the measured data and the modeled data. The model is modified such that data computed from the model fits to the measured data. The data itself is not used, except as a measure of the model data fit. A recently developed alternative method is to use results from inverse scattering theory to first construct an image while all multiple reflections are simultaneously eliminated from the data. This image can be constructed from surface reflection data if the data allows separating the subsurface reflection response from the down going emitted field. For 3D waves in a layered medium this requires knowledge of all horizontal electric and magnetic field components. If the data is properly sampled the solution is unique. In layered media the plane wave decomposition allows computing the image for each angle of incidence separately as a function of image time that is equal to the one-way intercept time. Once the image is constructed for all available angles of incidence a simple matrix inversion leads to the desired electric permittivity and magnetic permeability values in each layer. Finally these values provide interval velocities that can be used to convert image time to depth and the inverse problem is solved. The theory requires infinite bandwidth frequency domain data, which is equivalent to measuring the true impulse response. This is not possible in practice and numerical results show that data with finite bandwidths can be
Understanding how Fault-bounded Blocks Deform in 3D by Inverse Modelling
NASA Astrophysics Data System (ADS)
Jouen, G.; White, N.
2004-05-01
Normal faults play a crucial role in modifying basin stratigraphy. At the exploration scale, the internal deformation of tilted blocks is governed by the three-dimensional geometry of large-scale faults which bound these blocks. At the reservoir scale, the geometry and growth of normal faulting control the deformation of strata and the compartmentalisation of reservoir intervals. Despite their importance, large-scale normal faults are often difficult to image. The purpose of structural validation is two-fold: to determine the 3D shape of normal faults and to investigate the relationship between fault geometry and deformed stratigraphy including the intra-block faults. We have developed methods for tackling structural validation at a variety of scales in two and three dimensions. The cornerstone of our approach is the use of geophysical inverse theory to calculate optimal fault geometries from deformed strata. This approach allows us to focus on key questions: does a solution exist? Are there several possible solutions or just one unique one? In a complex normal fault system, which part of the fault controls the motion responsible for the deformation in the hanging-wall? Traditional forward modelling cannot answer these fundamental issues. We have applied the inversion on seismic data in particularly complex areas in the northern North Sea. The aims of this project are to determine the geometry of the basin-bounding fault, to assess the likelihood of out-of-plane motion as well as understanding the mode of deformation leading to the complexity of the present structure. Closely spaced inverse models show that the basin-bounding fault on the UK side is steeper and more planar than previously thought. This method also helped us to have a better view of what could have been the cause of the organisation and density of the intra-block faulting where it occurs. The North Cormorant study has shown how inverse modelling can yield important, quantitative, insights. Our
NASA Astrophysics Data System (ADS)
Meqbel, Naser M.; Egbert, Gary D.; Wannamaker, Philip E.; Kelbert, Anna; Schultz, Adam
2014-09-01
Long period (10-20,000 s) magnetotelluric (MT) data are being acquired across the continental USA on a quasi-regular grid of ˜70 km spacing as an electromagnetic component of the National Science Foundation EarthScope/USArray Program. These data are sensitive to fluids, melts, and other orogenic indicators, and thus provide a valuable complement to other components of EarthScope. We present and interpret results of 3-D MT data inversion from 325 sites acquired from 2006-2011 to provide a regional scale view of electrical resistivity from the middle crust to nearly the mantle transition zone, covering an area from NW Washington to NW Colorado. Beneath the active extensional subprovinces in the south-central region, on average we see a resistive upper crust, and then extensive areas of low resistivity in the lower crust and uppermost mantle. Further below, much of the upper half of the upper mantle appears moderately resistive, then subsequently the lower upper mantle becomes moderately conductive. This column suggests a dynamic process of moderately hydrated and fertile deeper upper mantle upwelling during extension, intersection of that material with the damp solidus causing dehydration and melting, and upward exodus of generated mafic melts to pond and exsolve saline fluids near Moho levels. Lithosphere here is very thin. To the east and northeast, thick sections of resistive lithosphere are imaged under the Wyoming and Medicine Hat Cratons. These are punctuated with numerous electrically conductive sutures presumably containing graphitic or sulfide-bearing meta-sediments deeply underthrust and emplaced during ancient collisions. Below Cascadia, the subducting Juan de Fuca and Gorda lithosphere appears highly resistive. Suspected oceanic lithosphere relicts in the central NW part of the model domain also are resistive, including the accreted “Siletzia” terrane beneath the Coast Ranges and Columbia Embayment, and the seismically fast “slab curtain” beneath
Xie, G.; Li, J.; Majer, E.; Zuo, D.
1998-07-01
This paper describes a new 3D parallel GILD electromagnetic (EM) modeling and nonlinear inversion algorithm. The algorithm consists of: (a) a new magnetic integral equation instead of the electric integral equation to solve the electromagnetic forward modeling and inverse problem; (b) a collocation finite element method for solving the magnetic integral and a Galerkin finite element method for the magnetic differential equations; (c) a nonlinear regularizing optimization method to make the inversion stable and of high resolution; and (d) a new parallel 3D modeling and inversion using a global integral and local differential domain decomposition technique (GILD). The new 3D nonlinear electromagnetic inversion has been tested with synthetic data and field data. The authors obtained very good imaging for the synthetic data and reasonable subsurface EM imaging for the field data. The parallel algorithm has high parallel efficiency over 90% and can be a parallel solver for elliptic, parabolic, and hyperbolic modeling and inversion. The parallel GILD algorithm can be extended to develop a high resolution and large scale seismic and hydrology modeling and inversion in the massively parallel computer.
Inks, T.L.; Agena, W.F.
2008-01-01
In February 2007, the Mt. Elbert Prospect stratigraphic test well, Milne Point, North Slope Alaska encountered thick methane gas hydrate intervals, as predicted by 3D seismic interpretation and modeling. Methane gas hydrate-saturated sediment was found in two intervals, totaling more than 100 ft., identified and mapped based on seismic character and wavelet modeling.
An MTBE plume in the Upper Glacial Aquifer of Long Island, NY was simulated by combining MODFLOW and MT3D with a semi-analytical model for a gasoline release. The first step was to develop and calibrate a 3-dimensional steady-state numerical ground water flow model of the aquife...
NASA Astrophysics Data System (ADS)
Preston, Leiph Alexander
We develop and apply a non-linear inversion of direct and wide-angle reflection travel times for 3-D P-wave velocity structure, earthquake hypocenters, and reflector geometry under NW Washington focusing on the structure of the subducting Juan de Fuca plate. The first-arrival travel times are derived from both active-source experiments and from local earthquakes. The reflection arrivals were picked from data collected during the 1998 Wet SHIPS active-source experiment, which consisted of air-gun sources within the inland water-ways of NW Washington and SW British Columbia to land-based stations. Our inversion procedure reduces the well-known trade-off between reflector position and the velocities above it by the combination of simultaneous inversion and adequate crossing paths. We interpret the wide-angle reflector as the Moho of the subducting Juan de Fuca slab. The relocated intraslab earthquakes separate into two groups: those located up-dip of the 45km reflector depth contour generally lie below the reflector in material whose velocity exceeds 7.7km/s, placing them within the subducting mantle, while those down-dip of this contour occur within material whose velocities are 6.8--7.5km/s, placing them within subducted oceanic crust. We interpret these groups of earthquakes as resulting from serpentine dehydration in the subducted mantle and the basalt to eclogite transformation in the subducted crust. We have performed velocity checkerboard, slab velocity resolution, and parameter sensitivity tests to estimate our ability to resolve the relationship among the reflector, intraslab hypocenters, and slab velocity structure. These tests indicate we have the necessary resolvability and can distinguish the relative locations among the velocities, reflector, and intraslab hypocenters within the subducting slab to +/-2km. The occurrence of events within the subducted mantle geometrically allows for larger magnitude earthquakes than could occur if they were confined to
Chern insulators without band inversion in Mo S2 monolayers with 3 d adatoms
NASA Astrophysics Data System (ADS)
Wei, Xinyuan; Zhao, Bao; Zhang, Jiayong; Xue, Yang; Li, Yun; Yang, Zhongqin
2017-02-01
Electronic and topological properties of Mo S2 monolayers endowed with 3 d transition metal (TM) adatoms (V-Fe) are explored by using ab initio methods and k .p models. Without the consideration of the Hubbard U interaction, the V, Cr, and Fe adatoms tend to locate on the top of the Mo atoms, while the most stable site for the Mn atom is at the hollow position of the Mo-S hexagon. After the Hubbard U is applied, the most stable sites of all the systems become the top of the Mo atoms. Chern insulators without band inversion are achieved in these systems. The V and Fe adsorption systems are the best candidates to produce the topological states. The k .p model calculations indicate that these topological states are determined by the TM magnetism, the C3 v crystal field from the Mo S2 substrate, and the TM atomic spin-orbit coupling (SOC). The special two-meron pseudospin texture is found to contribute to the topology. The apparent difference between the Berry curvatures for the V and Fe adsorption systems is also explored. Our results widen the understanding of the Chern insulators and are helpful for the applications of the Mo S2 monolayers in the future electronics and spintronics.
Statistical Inverse Ray Tracing for Image-Based 3D Modeling.
Liu, Shubao; Cooper, David B
2014-10-01
This paper proposes a new formulation and solution to image-based 3D modeling (aka "multi-view stereo") based on generative statistical modeling and inference. The proposed new approach, named statistical inverse ray tracing, models and estimates the occlusion relationship accurately through optimizing a physically sound image generation model based on volumetric ray tracing. Together with geometric priors, they are put together into a Bayesian formulation known as Markov random field (MRF) model. This MRF model is different from typical MRFs used in image analysis in the sense that the ray clique, which models the ray-tracing process, consists of thousands of random variables instead of two to dozens. To handle the computational challenges associated with large clique size, an algorithm with linear computational complexity is developed by exploiting, using dynamic programming, the recursive chain structure of the ray clique. We further demonstrate the benefit of exact modeling and accurate estimation of the occlusion relationship by evaluating the proposed algorithm on several challenging data sets.
Wang, G.L.; Chew, W.C.; Cui, T.J.; Aydiner, A.A.; Wright, D.L.; Smith, D.V.
2004-01-01
Three-dimensional (3D) subsurface imaging by using inversion of data obtained from the very early time electromagnetic system (VETEM) was discussed. The study was carried out by using the distorted Born iterative method to match the internal nonlinear property of the 3D inversion problem. The forward solver was based on the total-current formulation bi-conjugate gradient-fast Fourier transform (BCCG-FFT). It was found that the selection of regularization parameter follow a heuristic rule as used in the Levenberg-Marquardt algorithm so that the iteration is stable.
Bedekar, Vivek; Morway, Eric D.; Langevin, Christian D.; Tonkin, Matthew J.
2016-09-30
MT3D-USGS, a U.S. Geological Survey updated release of the groundwater solute transport code MT3DMS, includes new transport modeling capabilities to accommodate flow terms calculated by MODFLOW packages that were previously unsupported by MT3DMS and to provide greater flexibility in the simulation of solute transport and reactive solute transport. Unsaturated-zone transport and transport within streams and lakes, including solute exchange with connected groundwater, are among the new capabilities included in the MT3D-USGS code. MT3D-USGS also includes the capability to route a solute through dry cells that may occur in the Newton-Raphson formulation of MODFLOW (that is, MODFLOW-NWT). New chemical reaction Package options include the ability to simulate inter-species reactions and parent-daughter chain reactions. A new pump-and-treat recirculation package enables the simulation of dynamic recirculation with or without treatment for combinations of wells that are represented in the flow model, mimicking the above-ground treatment of extracted water. A reformulation of the treatment of transient mass storage improves conservation of mass and yields solutions for better agreement with analytical benchmarks. Several additional features of MT3D-USGS are (1) the separate specification of the partitioning coefficient (Kd) within mobile and immobile domains; (2) the capability to assign prescribed concentrations to the top-most active layer; (3) the change in mass storage owing to the change in water volume now appears as its own budget item in the global mass balance summary; (4) the ability to ignore cross-dispersion terms; (5) the definition of Hydrocarbon Spill-Source Package (HSS) mass loading zones using regular and irregular polygons, in addition to the currently supported circular zones; and (6) the ability to specify an absolute minimum thickness rather than the default percent minimum thickness in dry-cell circumstances.Benchmark problems that implement the new
NASA Astrophysics Data System (ADS)
Meqbel, N. M.; Egbert, G. D.; Kelbert, A.
2010-12-01
Long period (10-20,000 s) magnetotelluric (MT) data are being acquired in a series of temporary arrays deployed across the continental United States through EMScope, a component of EarthScope, a multidisciplinary decade-long project to study the structure and evolution of the North American Continent. MT deployments in 2006-2010 have so far acquired data at 237 sites on an approximately regular grid, with the same nominal spacing as the USArray broadband seismic transportable array (~70 km), covering the Northwestern US, from the Oregon-Washington coast across the Rocky Mountains, into Montana and Wyoming. Preliminary 3-D inversion results (Patro and Egbert; 2008), based on data from the 110 westernmost “Cascadia” sites collected in the first two years, revealed extensive areas of high conductivity in the lower crust beneath the Northwest Basin and Range (NBR), inferred to result from fluids (including possibly partial melt at depth) associated with magmatic underplating, and beneath the Cascade Mountains, probably due to fluids released by the subducting Juan de Fuca slab. Here we extend this study, refining and further testing the preliminary results from Cascadia, and extending the inversion domain to the East, to include all of the EarthScope data. Although site spacing is very broad, distinct regional structures are clearly evident even in simple maps of apparent resistivity, phase and induction vectors. For the 3-D inversion we are using the parallelized version of our recently developed Modular Code (ModEM), which supports Non-Linear Conjugate Gradient and several Gauss-Newton type schemes. Our initial 3-D inversion results using 212 MT sites, fitting impedances and vertical field transfer functions (together and separately) suggest several conductive and resistive structures which appear to be stable and required by the measured data. These include: - A conductive structure elongated in the N-S direction underneath the volcanic arc of the Cascadia
Burgess, Ian B; Aizenberg, Joanna; Lončar, Marko
2013-12-01
Structural hierarchy and complex 3D architecture are characteristics of biological photonic designs that are challenging to reproduce in synthetic materials. Top-down lithography allows for designer patterning of arbitrary shapes, but is largely restricted to planar 2D structures. Self-assembly techniques facilitate easy fabrication of 3D photonic crystals, but controllable defect-integration is difficult. In this paper we combine the advantages of top-down and bottom-up fabrication, developing two techniques to deposit 2D-lithographically-patterned planar layers on top of or in between inverse-opal 3D photonic crystals and creating hierarchical structures that resemble the architecture of the bright green wing scales of the butterfly, Parides sesostris. These fabrication procedures, combining advantages of both top-down and bottom-up fabrication, may prove useful in the development of omnidirectional coloration elements and 3D-2D photonic crystal devices.
A new model of the Arctic crustal thickness from 3D gravity inversion
NASA Astrophysics Data System (ADS)
Lebedeva-Ivanova, N. N.; Gaina, C.; Minakov, A.; Kashubin, S.
2015-12-01
The remarkable increase of new data collections and compilations for the Arctic region during the last decade motivate for a re-evaluation of our knowledge about the crustal structure and the tectonic evolution of the Arctic basins. 3D forward and inverse gravity modelling methods in the spectral domain (Minakov et al. 2012); lithosphere thermal gravity anomaly correction (Alvey et al., 2008); a vertical density variation for the sedimentary layer and lateral crustal variation density are integrated in the algorithm for derive the crustal thickness of the High Arctic region. Recently updated grids of bathymetry (Jakobsson et al., 2012), gravity anomaly (Gaina et al, 2011) and dynamic topography (Spasojevic & Gurnis, 2012) were used as input data for the algorithm. TeMAr sedimentary thickness grid (Petrov et al., 2015) was modified according to the most recent published seismic data, and was re-gridded and utilized as input data. Other input parameters for the algorithm were calibrated using seismic crustal scale profiles. Derived crustal thickness and Moho depth grids cover the area northward from 66° N and fit within a few kilometres with seismic crustal models for the most parts of the High Arctic region. Greater misfit in Moho depth between our results and seismic study (Chain & Lebedeva-Ivanova, 2015) under the northern Canada Basin suggest exceptional property of crust or/and mantel in this part of the Basin. Assumed mantle density of 3.25 kg/cm3provide the best fit for the region; it may indicate pervasive subcontinental lithospheric mantle (Goldstein et al., 2008) under the whole Arctic region. New results show a possible crustal connection between the Alpha and the Lomonosov ridges near the Canadian margin. The deepest Moho depth of c.34 km for Alpha-Mendeleev Ridge System is observed under the southern Mendeleev Ridge. The derived crustal thickness and Moho depth show a substantial improvement from the publicly available grids (CRUST1 (Laske et al., 2013
3D inversion of full gravity gradient tensor data using SL0 sparse recovery
NASA Astrophysics Data System (ADS)
Meng, Zhaohai
2016-04-01
We present a new method dedicated to the interpretation of full gravity gradient tensor data, based on SL0 sparse recovery inversion. The SL0 sparse recovery method aims to find out the minimum value of the objective function to fit the data function and to solve the non-zero solution to the objective function. Based on continuous iteration, we can easily obtain the final global minimum (namely the property and space attribute of the inversion target). We consider which type of tensor data combination produces the best inversion results based on the inversion results of different full gravity gradient tensor data combinations (separate tensor data and combined tensor data). We compare the recovered models obtained by inverting the different combinations of different gravity gradient tensor components to understand how different component combinations contribute to the resolution of the recovered model. Based on the comparison between the SL0 sparse recovery inversion results and the smoothest and focusing inversion results of the full gravity gradient tensor data, we show that SL0 sparse recovery inversion can obtain more stable and efficient inversion results with relatively sharp edge information, and that this method can also produce a stable solution of the inverse problem for complex geological structures. This new method to resolve very large full gravity gradient tensor datasets has the considerable advantage of being highly efficient; the full gravity gradient tensor inversion requires very little time. This new method is very effective in explaining the full gravity tensor which is very sensitive to small changes in local anomaly. The numerical simulation and inversion results of the compositional model indicates that including multiple components for inversion increases the resolution of the recovered density model and improves the structure delineation. We apply our inversion method to invert the gravity gradient tensor survey data from the Vinton salt
NASA Astrophysics Data System (ADS)
Meqbel, N. M.; Egbert, G. D.; Wannamaker, P. E.; Kelbert, A.; Schultz, A.
2013-12-01
Long period (10-20,000 s) magnetotelluric (MT) data are being acquired across the continental USA on a quasi-regular grid of ~70 km spacing as an electromagnetic component of the National Science Foundation EarthScope/USArray Program. These data are sensitive to fluids, melts, and other orogenic indicators, and thus provide a valuable complement to other components of EarthScope. We present and interpret results of 3-D MT data inversion from 325 sites acquired from 2006-2011 to provide a regional scale view of electrical resistivity from the middle crust to nearly the mantle transition zone, covering an area from NW Washington to NW Colorado. Extensive areas of low resistivity are imaged in the lower crust and uppermost mantle beneath the extensional provinces, most plausibly explained by underplated, hybridized magmas and associated exsolved highly saline fluids. These pervasive low resistivities show aligned or 'streaky' textures roughly parallel to seismic fast-axes, possibly reflecting widespread flow induced alignment of melt in this area. Thick sections of resistive lithosphere imaged in the eastern and northeastern part of the domain coincide spatially with the Wyoming and Medicine Hat Cratons. Sutures bounding these cratonic blocks are electrically conductive most likely due to meta-sediments emplaced during ancient collisions. Below the Cascadia forearc, the subducting Juan de Fuca and Gorda lithosphere appears highly resistive. Other resistive zones in the NW part of the domain may denote relict oceanic lithosphere: the accreted 'Siletzia' terrane beneath the Coast Ranges and Columbia Embayment, and the seismically fast 'slab curtain' beneath eastern Idaho interpreted by others as stranded Farallon lithosphere. Quasi-horizontal patches of low resistivity in the deep crust beneath the Cascade volcanic arc and fore-arc likely represent fluids evolved from breakdown of hydrous minerals in the down-going slab. In the backarc, low resistivities concentrate in
NASA Astrophysics Data System (ADS)
Torres-Verdin, C.
2007-05-01
This paper describes the successful implementation of a new 3D AVA stochastic inversion algorithm to quantitatively integrate pre-stack seismic amplitude data and well logs. The stochastic inversion algorithm is used to characterize flow units of a deepwater reservoir located in the central Gulf of Mexico. Conventional fluid/lithology sensitivity analysis indicates that the shale/sand interface represented by the top of the hydrocarbon-bearing turbidite deposits generates typical Class III AVA responses. On the other hand, layer- dependent Biot-Gassmann analysis shows significant sensitivity of the P-wave velocity and density to fluid substitution. Accordingly, AVA stochastic inversion, which combines the advantages of AVA analysis with those of geostatistical inversion, provided quantitative information about the lateral continuity of the turbidite reservoirs based on the interpretation of inverted acoustic properties (P-velocity, S-velocity, density), and lithotype (sand- shale) distributions. The quantitative use of rock/fluid information through AVA seismic amplitude data, coupled with the implementation of co-simulation via lithotype-dependent multidimensional joint probability distributions of acoustic/petrophysical properties, yields accurate 3D models of petrophysical properties such as porosity and permeability. Finally, by fully integrating pre-stack seismic amplitude data and well logs, the vertical resolution of inverted products is higher than that of deterministic inversions methods.
NASA Astrophysics Data System (ADS)
Zhu, Lupei; Zhou, Xiaofeng
2016-10-01
Source inversion of small-magnitude events such as aftershocks or mine collapses requires use of relatively high frequency seismic waveforms which are strongly affected by small-scale heterogeneities in the crust. In this study, we developed a new inversion method called gCAP3D for determining general moment tensor of a seismic source using Green's functions of 3D models. It inherits the advantageous features of the "Cut-and-Paste" (CAP) method to break a full seismogram into the Pnl and surface-wave segments and to allow time shift between observed and predicted waveforms. It uses grid search for 5 source parameters (relative strengths of the isotropic and compensated-linear-vector-dipole components and the strike, dip, and rake of the double-couple component) that minimize the waveform misfit. The scalar moment is estimated using the ratio of L2 norms of the data and synthetics. Focal depth can also be determined by repeating the inversion at different depths. We applied gCAP3D to the 2013 Ms 7.0 Lushan earthquake and its aftershocks using a 3D crustal-upper mantle velocity model derived from ambient noise tomography in the region. We first relocated the events using the double-difference method. We then used the finite-differences method and reciprocity principle to calculate Green's functions of the 3D model for 20 permanent broadband seismic stations within 200 km from the source region. We obtained moment tensors of the mainshock and 74 aftershocks ranging from Mw 5.2 to 3.4. The results show that the Lushan earthquake is a reverse faulting at a depth of 13-15 km on a plane dipping 40-47° to N46° W. Most of the aftershocks occurred off the main rupture plane and have similar focal mechanisms to the mainshock's, except in the proximity of the mainshock where the aftershocks' focal mechanisms display some variations.
Realistic 3D coherent transfer function inverse filtering of complex fields
Cotte, Yann; Toy, Fatih M.; Arfire, Cristian; Kou, Shan Shan; Boss, Daniel; Bergoënd, Isabelle; Depeursinge, Christian
2011-01-01
We present a novel technique for three-dimensional (3D) image processing of complex fields. It consists in inverting the coherent image formation by filtering the complex spectrum with a realistic 3D coherent transfer function (CTF) of a high-NA digital holographic microscope. By combining scattering theory and signal processing, the method is demonstrated to yield the reconstruction of a scattering object field. Experimental reconstructions in phase and amplitude are presented under non-design imaging conditions. The suggested technique is best suited for an implementation in high-resolution diffraction tomography based on sample or illumination rotation. PMID:21833359
Realistic 3D coherent transfer function inverse filtering of complex fields.
Cotte, Yann; Toy, Fatih M; Arfire, Cristian; Kou, Shan Shan; Boss, Daniel; Bergoënd, Isabelle; Depeursinge, Christian
2011-08-01
We present a novel technique for three-dimensional (3D) image processing of complex fields. It consists in inverting the coherent image formation by filtering the complex spectrum with a realistic 3D coherent transfer function (CTF) of a high-NA digital holographic microscope. By combining scattering theory and signal processing, the method is demonstrated to yield the reconstruction of a scattering object field. Experimental reconstructions in phase and amplitude are presented under non-design imaging conditions. The suggested technique is best suited for an implementation in high-resolution diffraction tomography based on sample or illumination rotation.
3D Gravity Inversion of Northern Sinai Peninsula: A Case Study
NASA Astrophysics Data System (ADS)
Khalil, Mohamed A.; Santos, Fernando M.
2014-07-01
The Sinai Peninsula has attracted the attention of many geological and geophysical studies as it is influenced and bounded by major tectonic events. Those are (1) the Mesozoic to Early Cenozoic tectonically active opening of the Tethys, (2) the Late Cretaceous to Early Tertiary (Laramide) Syrian arc system, due to closing of the Tethys (3) the Oligo-Miocene Gulf of Suez rifted basin, and (4) the Late Miocene to Recent transform Dead Sea-Gulf of Aqaba rift. Additionally, the shear zones inside Sinai such as the Ragabet El-Naam and Minsherah-Abu Kandu Shear Zones. Each of these major tectonic events has affected dramatically the structure evolution of the northern Sinai area. The present paper estimates the 3D density contrast model using the gravity data of northern Sinai. The estimated 3D density contrast model elucidated the peculiarities of the main structural elements in the region. The estimated 3D density contrast model showed the high and low gravity anomalies that form the main mountains and main valleys in northern Sinai. The estimated low density zones are in agreement with the inferred faults resulting from the first horizontal derivative. Comparing the 3D model with the tectonic history of the region and the results of the first horizontal derivative and least square separation increased the reliability of the model.
3-D Sound Propagation and Acoustic Inversions in Shallow Water Oceans
2012-12-19
fixed arc-length grid.] 10 Modeling comparisons Propagate over seamount , off center Source at 250 m, 100Hz 4 cases - (1) Nx2D, (2) Cartesian, (3...cylindrical PE. Figure 2. PE model comparisons I Sound propagation over a seamount are computed by different 3-D PE models, including (1) Nx2- D (2
NASA Astrophysics Data System (ADS)
Sladen, A.; Monteiller, V.
2014-12-01
Most large earthquakes are generated in subduction zones. To study the complexity of these events, teleseismic body waves offer many advantages over other types of data: they allow to study both the temporal and spatial evolution of slip during the rupture, they don't depend on the presence of nearby land and they allow to study earthquakes regardless of their location. Since the development of teleseismic finite-fault inversion in the 1980th, teleseismic body waves have been simulated using 1D velocity models to take into account propagation effects at the source. Yet, subduction zones are known to be highly heterogeneous: they are characterized by curved and dipping structures, strong seismic velocity contrasts, strong variations of topography and height of the water column. The main reason for relying on a 1D approximation is the computational cost of 3D simulations. And while forward simulations of teleseismic waves in a 3D Earth are only starting to be tractable on modern computers at the frequency range of interest (0.1Hz or shorter), finite-fault source studies require a large number of these simulations. In this work, we present a new and efficient approach to compute 3D teleseismic body waves, in which the full 3D propagation is only computed in a regional domain using discontinuous Galerkin finite-element method, while the rest of the seismic wave field is propagated in a background axisymmetric Earth. The regional and global wave fields are matched using the so-called Total-Field/Scattered-Field technique. This new simulation approach allows us to study the waveform complexities resulting from 3D propagation and investigate how they could improve the resolution and reduce the non-uniqueness of finite-fault inversions.
Electromagnetic Response Inversion for a 3D Distribution of Conductivity/Dielect
Newman, Gregory
2001-10-24
NLCGCS inverts electromagnetic responses for a 3D distribution of electrical conductivity and dielectric permittivity within the earth for geophysical applications using single processor computers. The software comes bundled with a graphical user interface to aid in model construction and analysis and viewing of earth images. The solution employs both dipole and finite size source configurations for harmonic oscillatory sources. A new nonlinear preconditioner is included in the solution to speed up solution convergence.
NASA Astrophysics Data System (ADS)
Santhanam, Anand P.; Min, Yugang; Mudur, Sudhir P.; Rastogi, Abhinav; Ruddy, Bari H.; Shah, Amish; Divo, Eduardo; Kassab, Alain; Rolland, Jannick P.; Kupelian, Patrick
2010-07-01
A method to estimate the deformation operator for the 3D volumetric lung dynamics of human subjects is described in this paper. For known values of air flow and volumetric displacement, the deformation operator and subsequently the elastic properties of the lung are estimated in terms of a Green's function. A Hyper-Spherical Harmonic (HSH) transformation is employed to compute the deformation operator. The hyper-spherical coordinate transformation method discussed in this paper facilitates accounting for the heterogeneity of the deformation operator using a finite number of frequency coefficients. Spirometry measurements are used to provide values for the airflow inside the lung. Using a 3D optical flow-based method, the 3D volumetric displacement of the left and right lungs, which represents the local anatomy and deformation of a human subject, was estimated from 4D-CT dataset. Results from an implementation of the method show the estimation of the deformation operator for the left and right lungs of a human subject with non-small cell lung cancer. Validation of the proposed method shows that we can estimate the Young's modulus of each voxel within a 2% error level.
Structural results for La Palma island using 3-D gravity inversion
NASA Astrophysics Data System (ADS)
Camacho, A. G.; FernáNdez, J.; GonzáLez, P. J.; Rundle, J. B.; Prieto, J. F.; Arjona, A.
2009-05-01
A recent gravity survey composed of 317 bench marks all over the island of La Palma (Canary Islands) is used, in combination with satellite data for regional aspects, to obtain results about structural properties of the island connected with the tectonic environment and local volcanism. To that end, a nonlinear three-dimensional gravity inversion approach is considered. The inversion scheme provides, in a nonsubjective form, the geometry of the anomalous bodies constructed in a random growth process. Results from the inversion can be interpreted in the framework of the geologic evolution of this ocean island volcano as a complex composite volcano with a large central body with high-density corresponding to the older intrusive part of the basalt complex. New unexpected features are enlightened, such as large thermal anomalies in the upper mantle southward of La Palma, as well as fracture en echelon zones associable to a slow active process of dislocation related to the recent volcanism in the southern half of the island. The results obtained for La Palma as a test site testify to the usefulness of the developed gravity inversion methodology for structural studies on islands in general.
NASA Astrophysics Data System (ADS)
Timur, Emre
2016-04-01
There are numerous geophysical methods used to investigate geothermal areas. The major purpose of this magnetic survey is to locate the boudaries of active hydrothermal system in the South of Gediz Graben in Salihli (Manisa/Turkey). The presence of the hydrothermal system had already been inferred from surface evidence of hydrothermal activity and drillings. Firstly, 3-D prismatic models were theoretically investigated and edge detection methods were utilized with an iterative inversion method to define the boundaries and the parameters of the structure. In the first step of the application, it was necessary to convert the total field anomaly into a pseudo-gravity anomaly map. Then the geometric boudaries of the structures were determined by applying a MATLAB based software with 3 different edge detection algorithms. The exact location of the structures were obtained by using these boundary coordinates as initial geometric parameters in the inversion process. In addition to these methods, reduction to pole and horizontal gradient methods were applied to the data to achieve more information about the location and shape of the possible reservoir. As a result, the edge detection methods were found to be successful, both in the field and as theoretical data sets for delineating the boundaries of the possible geothermal reservoir structure. The depth of the geothermal reservoir was determined as 2,4 km from 3-D inversion and 2,1 km from power spectrum methods.
NASA Astrophysics Data System (ADS)
Hara, Tatsuhiko
2004-08-01
We implement the Direct Solution Method (DSM) on a vector-parallel supercomputer and show that it is possible to significantly improve its computational efficiency through parallel computing. We apply the parallel DSM calculation to waveform inversion of long period (250-500 s) surface wave data for three-dimensional (3-D) S-wave velocity structure in the upper and uppermost lower mantle. We use a spherical harmonic expansion to represent lateral variation with the maximum angular degree 16. We find significant low velocities under south Pacific hot spots in the transition zone. This is consistent with other seismological studies conducted in the Superplume project, which suggests deep roots of these hot spots. We also perform simultaneous waveform inversion for 3-D S-wave velocity and Q structure. Since resolution for Q is not good, we develop a new technique in which power spectra are used as data for inversion. We find good correlation between long wavelength patterns of Vs and Q in the transition zone such as high Vs and high Q under the western Pacific.
Full-wave Moment Tensor and Tomographic Inversions Based on 3D Strain Green Tensor
2010-01-31
G. Jahnke, Wave propagation in 3D spherical sections: effects of subduction zones , Phys. Earth Planet. Inter., 132, 219-234, 2002. Komastitsch, D...is at scales smaller than the Fresnel zone . For example, a 1-Hz P/Pn wave recorded by a receiver ~1000 km from the source has a Fresnel zone width...approach, Eos Trans. AGU, 89(53), Fall Meet. Suppl., abstract T11E-06 Invited, 2008b. Sigloch, K., N. McQuarrie, G. Nolet, Two-stage subduction
Development of direct-inverse 3-D methods for applied transonic aerodynamic wing design and analysis
NASA Technical Reports Server (NTRS)
Carlson, Leland A.
1989-01-01
An inverse wing design method was developed around an existing transonic wing analysis code. The original analysis code, TAWFIVE, has as its core the numerical potential flow solver, FLO30, developed by Jameson and Caughey. Features of the analysis code include a finite-volume formulation; wing and fuselage fitted, curvilinear grid mesh; and a viscous boundary layer correction that also accounts for viscous wake thickness and curvature. The development of the inverse methods as an extension of previous methods existing for design in Cartesian coordinates is presented. Results are shown for inviscid wing design cases in super-critical flow regimes. The test cases selected also demonstrate the versatility of the design method in designing an entire wing or discontinuous sections of a wing.
Development and Tuning of a 3-D Stochastic Inversion Methodology for the European Arctic
2008-09-01
Norway and is subdivided into four tectonic nappes. Obduction started in the Vendian to Middle Cambrian and lasted until the Silurian . 2008...10 and 150 s period were combined with existing data provided by the University of Colorado at Boulder. This new data set was inverted for maps...showing the 2D group-velocity distribution of Love and Rayleigh waves for specific periods . Using a Monte Carlo inversion technique (Shapiro and
Optimization of computations for adjoint field and Jacobian needed in 3D CSEM inversion
NASA Astrophysics Data System (ADS)
Dehiya, Rahul; Singh, Arun; Gupta, Pravin K.; Israil, M.
2017-01-01
We present the features and results of a newly developed code, based on Gauss-Newton optimization technique, for solving three-dimensional Controlled-Source Electromagnetic inverse problem. In this code a special emphasis has been put on representing the operations by block matrices for conjugate gradient iteration. We show how in the computation of Jacobian, the matrix formed by differentiation of system matrix can be made independent of frequency to optimize the operations at conjugate gradient step. The coarse level parallel computing, using OpenMP framework, is used primarily due to its simplicity in implementation and accessibility of shared memory multi-core computing machine to almost anyone. We demonstrate how the coarseness of modeling grid in comparison to source (comp`utational receivers) spacing can be exploited for efficient computing, without compromising the quality of the inverted model, by reducing the number of adjoint calls. It is also demonstrated that the adjoint field can even be computed on a grid coarser than the modeling grid without affecting the inversion outcome. These observations were reconfirmed using an experiment design where the deviation of source from straight tow line is considered. Finally, a real field data inversion experiment is presented to demonstrate robustness of the code.
NASA Astrophysics Data System (ADS)
Ekinci, Yunus Levent; Ertekin, Can
2015-04-01
Concern about sedimentary basins is generally related to their genetic and economic significance. Analysis of sedimentary basins requires the acquisition of data through outcrop studies and subsurface investigations that encompass drilling and geophysics. These data are commonly analysed by computer-assisted techniques. One of these methods is based on analysing gravity anomalies to compute the depth of sedimentary basin-basement rock interface. Sedimentary basins produce negative gravity anomalies, because they have mostly lower densities than that of the surrounding basement rocks. Density variations in a sedimentary fill increase rapidly at shallower depths then gradually reach the density of surrounding basement rocks due to the geostatic pressure i.e. compaction. The decrease of the density contrast can be easily estimated by a quadratic function. Hence, if the densities are chosen properly and the regional background is removed correctly, the topographical relief of the sedimentary basin-basement rock interface might be estimated by the inversion of the gravity data using an exponential density-depth relation. Three dimensional forward modelling procedure can be carried out by introducing a Cartesian coordinate system, and placing vertical prisms just below observation points on the grid plane. Depth to the basement, namely depths to the bottom of the vertical prisms are adjusted in an iterative manner by minimizing the differences between measured and calculated residual gravity anomalies. In this study, we present a MATLAB-based inversion code for the interpretation of sedimentary basins by approximating the topographical relief of sedimentary basin-basement rock interfaces. For a given gridded residual gravity anomaly map, the procedure estimates the bottom depths of vertical prisms by considering some published formulas and assumptions. The utility of the developed inversion code was successfully tested on theoretically produced gridded gravity data set
Aydmer, A.A.; Chew, W.C.; Cui, T.J.; Wright, D.L.; Smith, D.V.; Abraham, J.D.
2001-01-01
A simple and efficient method for large scale three-dimensional (3-D) subsurface imaging of inhomogeneous background is presented. One-dimensional (1-D) multifrequency distorted Born iterative method (DBIM) is employed in the inversion. Simulation results utilizing synthetic scattering data are given. Calibration of the very early time electromagnetic (VETEM) experimental waveforms is detailed along with major problems encountered in practice and their solutions. This discussion is followed by the results of a large scale application of the method to the experimental data provided by the VETEM system of the U.S. Geological Survey. The method is shown to have a computational complexity that is promising for on-site inversion.
Bhattacharya, Jishnu; Wolverton, C
2013-05-07
Spinel oxides represent an important class of cathode materials for Li-ion batteries. Two major variants of the spinel crystal structure are normal and inverse. The relative stability of normal and inverse ordering at different stages of lithiation has important consequences in lithium diffusivity, voltage, capacity retention and battery life. In this paper, we investigate the relative structural stability of normal and inverse structures of the 3d transition metal oxide spinels with first-principles DFT calculations. We have considered ternary spinel oxides LixM2O4 with M = Ti, V, Cr, Mn, Fe, Co and Ni in both lithiated (x = 1) and delithiated (x = 0) conditions. We find that for all lithiated spinels, the normal structure is preferred regardless of the metal. We observe that the normal structure for all these oxides has a lower size mismatch between octahedral cations compared to the inverse structure. With delithiation, many of the oxides undergo a change in stability with vanadium in particular, showing a tendency to occupy tetrahedral sites. We find that in the delithiated oxide, only vanadium ions can access a +5 oxidation state which prefers tetrahedral coordination. We have also calculated the average voltage of lithiation for these spinels. The calculated voltages agree well with the previously measured and calculated values, wherever available. For the yet to be characterized spinels, our calculation provides voltage values which can motivate further experimental attention. Lastly, we observe that all the normal spinel oxides of the 3d transition metal series have a driving force for a transformation to the non-spinel structure upon delithiation.
Dirac Circles and Quantum Hall Effect in 3D Inversion-Symmetric Crystals
NASA Astrophysics Data System (ADS)
Wieder, Benjamin J.; Kim, Youngkuk; Kane, C. L.
2015-03-01
In the presence of inversion and time-reversal symmetries, materials with weak spin-orbit coupling may host topologically protected Dirac line nodes. A band inversion transition in these systems can produce a line node which closes on itself and forms a protected Dirac circle. The surfaces parallel to this circle host zero-energy puddles in momentum space which are flat if the inverting bands have the same effective mass. In cases with differing effective masses, the surface modes disperse, but the bulk Dirac circle remains gapless. Adding an external magnetic field perpendicular to this circle creates surface Landau levels, whose number can be controlled by tuning the field strength. When a new level is created or destroyed, the bulk becomes gapless and the zero-temperature bulk conductivity displays a sharp peak. The sequence of conductivity peaks describes an unusual manifestation of the integer quantum hall effect. We characterize surface and bulk transport as a function of magnetic field strength and in the presence of disorder.
A numerical method for the inverse problem of cell traction in 3D
NASA Astrophysics Data System (ADS)
Vitale, G.; Preziosi, L.; Ambrosi, D.
2012-09-01
Force traction microscopy is an inversion method that allows us to obtain the stress field applied by a living cell on the environment on the basis of a pointwise knowledge of the displacement produced by the cell itself. This classical biophysical problem, usually addressed in terms of Green’s functions, can be alternatively tackled in a variational framework. In such a case, a variation of the error functional under suitable regularization is operated in view of its minimization. This setting naturally suggests the introduction of a new equation, based on the adjoint operator of the elasticity problem. In this paper, we illustrate a numerical strategy of the inversion method that discretizes the partial differential equations associated with the optimal control problem by finite elements. A detailed discussion of the numerical approximation of a test problem (with known solution) that contains most of the mathematical difficulties of the real one allows a precise evaluation of the degree of confidence that one can achieve in the numerical results.
MT2DInvMatlab—A program in MATLAB and FORTRAN for two-dimensional magnetotelluric inversion
NASA Astrophysics Data System (ADS)
Lee, Seong Kon; Kim, Hee Joon; Song, Yoonho; Lee, Choon-Ki
2009-08-01
MT2DInvMatlab is an open-source MATLAB® software package for two-dimensional (2D) inversion of magnetotelluric (MT) data; it is written in mixed languages of MATLAB and FORTRAN. MT2DInvMatlab uses the finite element method (FEM) to compute 2D MT model responses, and smoothness-constrained least-squares inversion with a spatially variable regularization parameter algorithm to stabilize the inversion process and provide a high-resolution optimal earth model. It is also able to include terrain effects in inversion by incorporating topography into a forward model. This program runs under the MATLAB environment so that users can utilize the existing general interface of MATLAB, while some specific functions are written in FORTRAN 90 to speed up computation and reuse pre-existing FORTRAN code in the MATLAB environment with minimal modification. This program has been tested using synthetic models, including one with variable topography, and on field data. The results were assessed by comparing inverse models obtained with MT2DInvMatlab and with a non-linear conjugate gradient (NLCG) algorithm. In both tests the new inversion software reconstructs the subsurface resistivity structure very closely and provides an improvement in both resolution and stability.
NASA Astrophysics Data System (ADS)
Götze, Hans-Jürgen; Choi, Sungchan
2015-04-01
We combined the global gravity dataset EGM2008 and a local terrestrial gravity data survey to conduct constrained 3-D crustal density modeling of a strato-volcanic complex and the surrounding area located close to the border of North Korea and China. The independent geophysical (seismic, seismology, geochemistry) and petrological constraints will be presented together with the preprocessing of data base by curvature analysis and Euler deconvolution. The multiple data base is used to assist a general interpretation of the investigated area, and the 3D density model (modelled by the in-house IGMAS+ software). Mt. Paekdu is characterized by a low of Bouguer anomaly of some -110 × 10-5 m/s2, which is caused by the combined gravity effects of (1) Moho depth of about 40 km, (2) a zone with both lower P-wave velocity and density than the surrounding, (3) low density volcanic rocks at the surface, and (4) the presence of a magma chamber that has not previously been identified. The terrestrial gravity field measured along the seismic profile shows a remarkable anomaly descending from the southern- to the northern flank of the Mt. Paekdu volcano, which should be a typical anomaly pattern generally observed over the active volcanic area in the world (e.g. the Yellow Stone volcano). The trend is interpreted to be caused by a prominent density difference between a serious of high density mid crustal sill beneath the southern flank and a predicted partial melted zone locating in the northern flank. With the help of several geoscientific observations (seismic, electromagnetic, gravity and geochemistry) and the 3D density model we conclude that a high density sill was formed in Pliocene and early Pleistocene after pre-shield plateau-forming eruption. Since the Pliocene, volcanic activity in the Mt. Paekdu region might be migrated from the southeastern of North Korea to the northwest, following the path of NW-SE-trending faults. Recently observed seismic tremors can be explained
NASA Astrophysics Data System (ADS)
Goetze, H. J.; Choi, S.
2014-12-01
In the presentation we get use of the global gravity dataset EGM2008 and a local terrestrial gravity data survey for a constrained 3-D crustal density modeling of a stratovolcano and its surrounding area located close to the border of North Korea and China. The independent geophysical (seismic, seismology, geochemistry) and petrological constraints will be presented together with the preprocessing of data base by curvature analysis and Euler deconvolution. The multiple data base is used to assist a general interpretation of the investigated area in time, and the 3D density model (modelled by the inhouse IGMAS+ software). Mt. Paekdu is characterized by a low of Bouguer anomaly of some -110 ´ 10-5 m/s2, which is caused by the combined gravity effects of (1) Moho depth of about 40 km, (2) a zone with both lower P-wave velocity and density than the surrounding, (3) low density volcanic rocks at the surface, and (4) the presence of a magma chamber that has not previously been identified. The terrestrial gravity field measured along the seismic profile shows a remarkable anomaly descending from the southern- to the northern flank of the Mt. Paekdu volcano, which should be a typical anomaly pattern generally obsered over the active volcanic area in the world (e.g. the Yellow Stone volcano). The trend is interpreted to be caused by a prominent density difference between a serious of high density mid crustal sill beneath the southern flank and a predicted partial melted zone locating in the northern flank. With the help of several geoscientific observations (seismic, electromagnetic, gravity and geochemistry) and the 3D density model we conclude that a high density sill was formed in Pliocene and early Pleistocene after pre-shield plateau-forming eruption. Since the Pliocene, volcanic activity in the Mt. Paekdu region might be migrated from the southeastern of North Korea to the northwest, following the path of NW-SE-trending faults. Recently observed seismic tremors can
NASA Astrophysics Data System (ADS)
Bertrand, E. A.; Caldwell, G.; Bannister, S. C.; Hill, G.; Bennie, S.
2013-12-01
The Taupo Volcanic Zone (TVZ), located in the central North Island of New Zealand, is a rifted arc that contains more than 20 liquid-dominated high-temperature geothermal systems, which together discharge ~4.2 GW of heat at the surface. The shallow (upper ~500 m) extent of these geothermal systems is marked by low-resistivity, mapped by tens-of-thousands of DC resistivity measurements collected throughout the 1970's and 80's. Conceptual models of heat transport through the brittle crust of the TVZ link these low-resistivity anomalies to the tops of vertically ascending plumes of convecting hydrothermal fluid. Recently, data from a 40-site array of broadband seismometers with ~4 km station spacing, and an array of 270 broadband magnetotelluric (MT) measurements with ~2 km station spacing, have been collected in the south-eastern part of the TVZ in an experiment to image the deep structure (or roots) of the geothermal systems in this region. Unlike DC resistivity, these MT measurements are capable of resolving the resistivity structure of the Earth to depths of 10 km or more. 2-D and 3-D models of subsets of these MT data have been used to provide the first-ever images of quasi-vertical low-resistivity zones (at depths of 3-7 km) that connect with the near-surface geothermal fields. These low-resistivity zones are interpreted to represent convection plumes of high-temperature fluids ascending within fractures, which supply heat to the overlying geothermal fields. At the Rotokawa, Ngatamariki and Ohaaki geothermal fields, these plumes extend to a broad layer of low-resistivity, inferred to represent a magmatic, basal heat source located below the seismogenic zone (at ~7-8 km depth) that drives convection in the brittle crust above. Little is known about the mechanisms that transfer heat into the hydrothermal regime. However, at Rotokawa, new 3-D resistivity models image a vertical low-resistivity zone that lies directly beneath the geothermal field. The top of this
Barlebo, H.C.; Rosbjerg, D.; Hill, M.C.
1996-01-01
An extensive amount of data including hydraulic heads, hydraulic conductivities and concentrations of several solutes from controlled injections have been collected during the MADE 1 and MADE 2 experiments at a heterogeneous site near Columbus, Mississippi. In this paper the use of three-dimensional inverse groundwater models including simultaneous estimation of flow and transport parameters is proposed to help identify the dominant characteristics at the site. Simulations show that using a hydraulic conductivity distribution obtained from 2187 borehole flowmeter tests directly in the model produces poor matches to the measured hydraulic heads and tritium concentrations. Alternatively, time averaged hydraulic head maps are used to define zones of constant hydraulic conductivity to be estimated. Preliminary simulations suggest that in the case of conservative transport many, but not all, of the major plume characteristics can be explained by large-scale heterogeneity in recharge and hydraulic conductivity.
Development of direct-inverse 3-D methods for applied transonic aerodynamic wing design and analysis
NASA Technical Reports Server (NTRS)
Carlson, Leland A.
1989-01-01
Progress in the direct-inverse wing design method in curvilinear coordinates has been made. This includes the remedying of a spanwise oscillation problem and the assessment of grid skewness, viscous interaction, and the initial airfoil section on the final design. It was found that, in response to the spanwise oscillation problem that designing at every other spanwise station produced the best results for the cases presented, a smoothly varying grid is especially needed for the accurate design at the wing tip, the boundary layer displacement thicknesses must be included in a successful wing design, the design of high and medium aspect ratio wings is possible with this code, and the final airfoil section designed is fairly independent of the initial section.
Hall-Effect Sign Inversion in a Realizable 3D Metamaterial
NASA Astrophysics Data System (ADS)
Kadic, Muamer; Schittny, Robert; Bückmann, Tiemo; Kern, Christian; Wegener, Martin
2015-04-01
In 2009, Briane and Milton proved mathematically the existence of three-dimensional isotropic metamaterials with a classical Hall coefficient that is negative with respect to that of all of the metamaterial constituents. Here, we significantly simplify their blueprint towards an architecture composed of only a single-constituent material in vacuum or air, which can be seen as a special type of porosity. We show numerically that the sign of the Hall voltage is determined by a separation parameter between adjacent tori. This qualitative behavior is robust even for only a small number of metamaterial unit cells. The combination of simplification and robustness brings experimental verification of this striking sign inversion into reach. Furthermore, we provide a simple intuitive explanation of the underlying physical mechanism.
Brosten, T.R.; Day-Lewis, F. D.; Schultz, G.M.; Curtis, G.P.; Lane, J.W.
2011-01-01
Electromagnetic induction (EMI) instruments provide rapid, noninvasive, and spatially dense data for characterization of soil and groundwater properties. Data from multi-frequency EMI tools can be inverted to provide quantitative electrical conductivity estimates as a function of depth. In this study, multi-frequency EMI data collected across an abandoned uranium mill site near Naturita, Colorado, USA, are inverted to produce vertical distribution of electrical conductivity (EC) across the site. The relation between measured apparent electrical conductivity (ECa) and hydraulic conductivity (K) is weak (correlation coefficient of 0.20), whereas the correlation between the depth dependent EC obtained from the inversions, and K is sufficiently strong to be used for hydrologic estimation (correlation coefficient of -0.62). Depth-specific EC values were correlated with co-located K measurements to develop a site-specific ln(EC)-ln(K) relation. This petrophysical relation was applied to produce a spatially detailed map of K across the study area. A synthetic example based on ECa values at the site was used to assess model resolution and correlation loss given variations in depth and/or measurement error. Results from synthetic modeling indicate that optimum correlation with K occurs at ~0.5m followed by a gradual correlation loss of 90% at 2.3m. These results are consistent with an analysis of depth of investigation (DOI) given the range of frequencies, transmitter-receiver separation, and measurement errors for the field data. DOIs were estimated at 2.0??0.5m depending on the soil conductivities. A 4-layer model, with varying thicknesses, was used to invert the ECa to maximize available information within the aquifer region for improved correlations with K. Results show improved correlation between K and the corresponding inverted EC at similar depths, underscoring the importance of inversion in using multi-frequency EMI data for hydrologic estimation. ?? 2011.
NASA Astrophysics Data System (ADS)
Balkaya, Çağlayan; Ekinci, Yunus Levent; Göktürkler, Gökhan; Turan, Seçil
2017-01-01
3D non-linear inversion of total field magnetic anomalies caused by vertical-sided prismatic bodies has been achieved by differential evolution (DE), which is one of the population-based evolutionary algorithms. We have demonstrated the efficiency of the algorithm on both synthetic and field magnetic anomalies by estimating horizontal distances from the origin in both north and east directions, depths to the top and bottom of the bodies, inclination and declination angles of the magnetization, and intensity of magnetization of the causative bodies. In the synthetic anomaly case, we have considered both noise-free and noisy data sets due to two vertical-sided prismatic bodies in a non-magnetic medium. For the field case, airborne magnetic anomalies originated from intrusive granitoids at the eastern part of the Biga Peninsula (NW Turkey) which is composed of various kinds of sedimentary, metamorphic and igneous rocks, have been inverted and interpreted. Since the granitoids are the outcropped rocks in the field, the estimations for the top depths of two prisms representing the magnetic bodies were excluded during inversion studies. Estimated bottom depths are in good agreement with the ones obtained by a different approach based on 3D modelling of pseudogravity anomalies. Accuracy of the estimated parameters from both cases has been also investigated via probability density functions. Based on the tests in the present study, it can be concluded that DE is a useful tool for the parameter estimation of source bodies using magnetic anomalies.
Enhanced imaging of CO2 at the Ketzin storage site: Inversion of 3D time-lapse seismic data
NASA Astrophysics Data System (ADS)
Gil, M.; Götz, J.; Ivanova, A.; Juhlin, C.; Krawczyk, C. M.; Lüth, S.; Yang, C.
2012-04-01
The Ketzin test site, located near Berlin, is Europe's longest-operating on-shore CO2 storage site. As of December 2011, more than 56,000 tons of food grade CO2 has been injected since June 2008 in an anticlinal structure of the Northeast German Basin. The target reservoir consists of porous, brine bearing sandstone units of the Upper Triassic Stuttgart Formation at approximately 630 to 650 m depth. In order to enhance the understanding of the structural geometry of the site and to investigate the extension of the CO2-plume, several geophysical monitoring methods are being applied at Ketzin, among these are active seismic measurements, geoelectrics and borehole measurements. Among the various seismic techniques (e.g. 2D reflection surveys, crosshole tomography, Vertical Seismic Profiling, 2D- and 3D-Moving Source Profiling) employed at this pilot site, 3D time-lapse reflection surveys are an important component. The baseline 3D survey was acquired in 2005 and the first repeat measurements were performed in 2009 after injection of about 22,000 tons of CO2. The second repeat survey is planned to be carried out in fall 2012. These measurements allow the time-lapse signature of the injected CO2 to be imaged. The time-lapse amplitude variation attributed to the injected CO2 in the reservoir matches, considering detection limits of seismic surface measurements, the expected distribution of the CO2 plume derived from reservoir simulations. Previous attempts towards a quantitative interpretation were based on integrative considerations of different types of geophysical measurements using strict assumptions and characterized by large error bars. In order to increase the resolution and reliability of the data and to improve estimation of rock properties and especially to enhance the imaging resolution of the CO2-plume, the time-lapse 3D seismic data have now been inverted for seismic impedances with different methods, which is the focus of this presentation. One difficulty
NASA Astrophysics Data System (ADS)
Mair, H. D.; Ciorau, P.; Owen, D.; Hazelton, T.; Dunning, G.
2000-05-01
Two ultrasonic simulation packages: Imagine 3D and SIMSCAN have specifically been developed to solve the inverse problem for blade root and rotor steeple of low-pressure turbine. The software was integrated with the 3D drawing of the inspected parts, and with the dimensions of linear phased-array probes. SIMSCAN simulates the inspection scenario in both optional conditions: defect location and probe movement/refracted angle range. The results are displayed into Imagine 3-D, with a variety of options: rendering, display 1:1, grid, generated UT beam. The results are very useful for procedure developer, training and to optimize the phased-array probe inspection sequence. A spreadsheet is generated to correlate the defect coordinates with UT data (probe position, skew and refracted angle, UT path, and probe movement). The simulation models were validated during experimental work with phased-array systems. The accuracy in probe position is ±1 mm, and the refracted/skew angle is within ±0.5°. Representative examples of phased array focal laws/probe movement for a specific defect location, are also included.
NASA Astrophysics Data System (ADS)
Petrov, P.; Newman, G. A.
2011-12-01
Recent developments in high resolution imaging technology of subsurface objects involves a combination of different geophysical measurements (gravity, EM and seismic). A joint image of the subsurface geophysical attributes (velocity, electrical conductivity and density) requires the consistent treatment of the different geophysical data due to their differing physical nature. For example, in conducting media, which is typical of the Earth's interior, EM energy propagation is defined by a diffusive mechanism and may be characterized by two specific length scales: wavelength and skin depth. However, the propagation of seismic signals is a multiwave process and is characterized by a set of wavelengths. Thus, to consistently treat seismic and electromagnetic data an additional length scale is needed for seismic data that does not directly depend on a wavelength and describes a diffusive process, similar to EM wave propagation in the subsurface. Works by Brown et al.(2005), Shin and Cha(2008), and Shin and Ha(2008) suggest that an artificial damping of seismic wave fields via Laplace-Fourier transformation can be an effective approach to obtain a seismic data that have similar spatial resolution to EM data. The key benefit of such transformation is that diffusive wave-field inversion works well for both data sets: seismic (Brown et al.,2005; Shin and Cha,2008) and electromagnetic (Commer and Newman,2008; Newman et al.,2010). With the recent interest in the Laplace-Fourier domain full waveform inversion, 3D fourth and second-order finite-difference schemes for modeling of seismic wave propagation have been developed (Petrov and Newman, 2010). Incorporation of attenuation and anisotropy into a velocity model is a necessary step for a more realistic description of subsurface media. Here we consider the extension of our method which includes attenuation and VTI anisotropy. Our approach is based on the integro-interpolation technique for velocity-stress formulation. Seven
NASA Astrophysics Data System (ADS)
Jesús Moral García, Francisco; Rebollo Castillo, Francisco Javier; Monteiro Santos, Fernando
2016-04-01
Maps of apparent electrical conductivity of the soil are commonly used in precision agriculture to indirectly characterize some important properties like salinity, water, and clay content. Traditionally, these studies are made through an empirical relationship between apparent electrical conductivity and properties measured in soil samples collected at a few locations in the experimental area and at a few selected depths. Recently, some authors have used not the apparent conductivity values but the soil bulk conductivity (in 2D or 3D) calculated from measured apparent electrical conductivity through the application of an inversion method. All the published works used data collected with electromagnetic (EM) instruments. We present a new software to invert the apparent electrical conductivity data collected with VERIS 3100 and 3150 (or the more recent version with three pairs of electrodes) using the 1D spatially constrained inversion method (1D SCI). The software allows the calculation of the distribution of the bulk electrical conductivity in the survey area till a depth of 1 m. The algorithm is applied to experimental data and correlations with clay and water content have been established using soil samples collected at some boreholes. Keywords: Digital soil mapping; inversion modelling; VERIS; soil apparent electrical conductivity.
NASA Astrophysics Data System (ADS)
Oh, Ju-Won; Alkhalifah, Tariq
2016-09-01
Multiparameter full waveform inversion (FWI) applied to an elastic orthorhombic model description of the subsurface requires in theory a nine-parameter representation of each pixel of the model. Even with optimal acquisition on the Earth surface that includes large offsets, full azimuth, and multicomponent sensors, the potential for trade-off between the elastic orthorhombic parameters are large. The first step to understanding such trade-off is analysing the scattering potential of each parameter, and specifically, its scattering radiation patterns. We investigate such radiation patterns for diffraction and for scattering from a horizontal reflector considering a background isotropic model. The radiation patterns show considerable potential for trade-off between the parameters and the potentially limited resolution in their recovery. The radiation patterns of C11, C22, and C33 are well separated so that we expect to recover these parameters with limited trade-offs. However, the resolution of their recovery represented by recovered range of model wavenumbers varies between these parameters. We can only invert for the short wavelength components (reflection) of C33 while we can mainly invert for the long wavelength components (transmission) of the elastic coefficients C11 and C22 if we have large enough offsets. The elastic coefficients C13, C23, and C12 suffer from strong trade-offs with C55, C44, and C66, respectively. The trade-offs between C13 and C55, as well as C23 and C44, can be partially mitigated if we acquire P-SV and SV-SV waves. However, to reduce the trade-offs between C12 and C66, we require credible SH-SH waves. The analytical radiation patterns of the elastic constants are supported by numerical gradients of these parameters.
Star-Lack, J; Nelson, S J; Kurhanewicz, J; Huang, L R; Vigneron, D B
1997-08-01
A T1 insensitive solvent suppression technique-band selective inversion with gradient dephasing (BASING)-was developed to suppress water and lipids for 1H magnetic resonance spectroscopy (MRS). BASING, which consists of a frequency selective RF inversion pulse surrounded by spoiler gradient pulses of opposite signs, was used to dephase stopband resonances and minimally impact passband metabolites. Passband phase linearity was achieved with a dual BASING scheme. Using the Shinnar-Le Roux algorithm, a highpass filter was designed to suppress water and rephase the lactate methyl doublet independently of TE, and water/lipid bandstop filters were designed for the brain and prostate. Phantom and in vivo experimental 3D PRESS CSI data were acquired at 1.5 T to compare BASING with CHESS and STIR suppression. With BASING, the measured suppression factor was over 100 times higher than with CHESS or STIR causing baseline distortions to be removed. It was shown that BASING can be incorporated into a variety of sequences to offer improved suppression in the presence of B1 and T1 inhomogeneites.
NASA Astrophysics Data System (ADS)
Windhari, Ayuty; Handayani, Gunawan
2015-04-01
The 3D inversion gravity anomaly to estimate topographical density using a matlab source code from gridded data provided by Parker Oldenburg algorithm based on fast Fourier transform was computed. We extend and improved the source code of 3DINVERT.M invented by Gomez Ortiz and Agarwal (2005) using the relationship between Fourier transform of the gravity anomaly and the sum of the Fourier transform from the topography density. We gave density contrast between the two media to apply the inversion. FFT routine was implemented to construct amplitude spectrum to the given mean depth. The results were presented as new graphics of inverted topography density, the gravity anomaly due to the inverted topography and the difference between the input gravity data and the computed ones. It terminates when the RMS error is lower than pre-assigned value used as convergence criterion or until maximum of iterations is reached. As an example, we used the matlab program on gravity data of Banten region, Indonesia.
Łeski, Szymon; Wójcik, Daniel K; Tereszczuk, Joanna; Swiejkowski, Daniel A; Kublik, Ewa; Wróbel, Andrzej
2007-01-01
Estimation of the continuous current-source density in bulk tissue from a finite set of electrode measurements is a daunting task. Here we present a methodology which allows such a reconstruction by generalizing the one-dimensional inverse CSD method. The idea is to assume a particular plausible form of CSD within a class described by a number of parameters which can be estimated from available data, for example a set of cubic splines in 3D spanned on a fixed grid of the same size as the set of measurements. To avoid specificity of particular choice of reconstruction grid we add random jitter to the points positions and show that it leads to a correct reconstruction. We propose different ways of improving the quality of reconstruction which take into account the sources located outside the recording region through appropriate boundary treatment. The efficiency of the traditional CSD and variants of inverse CSD methods is compared using several fidelity measures on different test data to investigate when one of the methods is superior to the others. The methods are illustrated with reconstructions of CSD from potentials evoked by stimulation of a bunch of whiskers recorded in a slab of the rat forebrain on a grid of 4x5x7 positions.
NASA Astrophysics Data System (ADS)
Afonso, J. C.; Fullea, J.; Yang, Y.; Griffin, W. L.; Jones, A. G.; Connolly, J.; Lebedev, S.; O'Reilly, S. Y.
2011-12-01
High-resolution imaging and characterization of the thermal and compositional structure of the lithospheric and sublithospheric upper mantle are the basis for understanding the formation and evolution of the lithosphere and the interaction between the crust-mantle and lithosphere-asthenosphere systems. Unfortunately, such imaging and characterization using available geophysical-geochemical methods still present unsolved and technically challenging problems. In this contribution we present a new full-3D multi-observable inversion method particularly designed for high-resolution (regional) thermal and compositional mapping of the lithosphere and sublithospheric upper mantle. Ambient noise tomography, multiple plane wave earthquake tomography, magnetotelluric, thermal, thermodynamic, and potential field modelling are all combined in a single thermodynamic-geophysical framework and appraised within a general probabilistic (Bayesian) formulation. This circumvents the problems of strong non-linearity involved in traditional inversions, provides highly refined seismic information, minimizes the problem of trade-off between temperature and composition in wave speeds, offers critical insights into incompatibilities between traditional stand-alone methods, and takes advantage of a priori local geochemical information. Both synthetic models and preliminary results in real-case examples will be used to discuss the benefits, robustness, and limitations of this method.
NASA Astrophysics Data System (ADS)
Choi, S.; Kim, C.; Kim, H. R.; Park, C.; Park, H. Y.
2015-12-01
We performed the marine magnetic and the bathymetry survey in the Lau basin for finding the submarine hydrothermal deposits in October 2009. We acquired magnetic and bathymetry datasets by using Overhouser Proton Magnetometer SeaSPY(Marine Magnetics Co.) and Multi-Beam Echo Sounder EM120(Kongsberg Co.). We conducted the data processing to obtain detailed seabed topography, magnetic anomaly and reduction to the pole(RTP). The Lau basin is one of the youngest back-arc basins in the Southwest Pacific. This region was a lot of hydrothermal activities and hydrothermal deposits. In particular, Tofua Arc(TA) in the Lau basin consists of various and complex stratovolcanos(from Massoth et al., 2007).), We calculated the magnetic susceptibility distribution of the TA19-1 seamount(longitude:176°23.5'W, latitude: 22°42.5'W)area using the RTP data by 3-D magnetic inversion from Jung's previous study(2013). Based on 2D 'compact gravity inversion' by Last & Kubik(1983), we expend it to the 3D algorithm using iterative reweighted least squares method with some weight matrices. The used weight matrices are two types: 1) the minimum gradient support(MGS) that controls the spatial distribution of the solution from Porniaguine and Zhdanov(1999); 2) the depth weight that are used according to the shape of subsurface structures. From the modeling, we derived the appropriate scale factor for the use of depth weight and setting magnetic susceptibility. Furthermore, we have to enter a very small error value to control the computation of the singular point of the inversion model that was able to be easily calculated for modeling. In addition, we applied separately weighted value for the correct shape and depth of the magnetic source. We selected the best results model by change to converge of RMS. Compared between the final modeled result and RTP values in this study, they are generally similar to the each other. But the input values and the modeled values have slightly little difference
NASA Astrophysics Data System (ADS)
Bignardi, S.; Mantovani, A.; Abu Zeid, N.
2016-08-01
OpenHVSR is a computer program developed in the Matlab environment, designed for the simultaneous modeling and inversion of large Horizontal-to-Vertical Spectral Ratio (HVSR or H/V) datasets in order to construct 2D/3D subsurface models (topography included). The program is designed to provide a high level of interactive experience to the user and still to be of intuitive use. It implements several effective and established tools already present in the code ModelHVSR by Herak (2008), and many novel features such as: -confidence evaluation on lateral heterogeneity -evaluation of frequency dependent single parameter impact on the misfit function -relaxation of Vp/Vs bounds to allow for water table inclusion -a new cost function formulation which include a slope dependent term for fast matching of peaks, which greatly enhances convergence in case of low quality HVSR curves inversion -capability for the user of editing the subsurface model at any time during the inversion and capability to test the changes before acceptance. In what follows, we shall present many features of the program and we shall show its capabilities on both simulated and real data. We aim to supply a powerful tool to the scientific and professional community capable of handling large sets of HSVR curves, to retrieve the most from their microtremor data within a reduced amount of time and allowing the experienced scientist the necessary flexibility to integrate into the model their own geological knowledge of the sites under investigation. This is especially desirable now that microtremor testing has become routinely used. After testing the code over different datasets, both simulated and real, we finally decided to make it available in an open source format. The program is available by contacting the authors.
Pavei, Gaspare; Seminati, Elena; Cazzola, Dario; Minetti, Alberto E.
2017-01-01
The dynamics of body center of mass (BCoM) 3D trajectory during locomotion is crucial to the mechanical understanding of the different gaits. Forward Dynamics (FD) obtains BCoM motion from ground reaction forces while Inverse Dynamics (ID) estimates BCoM position and speed from motion capture of body segments. These two techniques are widely used by the literature on the estimation of BCoM. Despite the specific pros and cons of both methods, FD is less biased and considered as the golden standard, while ID estimates strongly depend on the segmental model adopted to schematically represent the moving body. In these experiments a single subject walked, ran, (uni- and bi-laterally) skipped, and race-walked at a wide range of speeds on a treadmill with force sensors underneath. In all conditions a simultaneous motion capture (8 cameras, 36 markers) took place. 3D BCoM trajectories computed according to five marker set models of ID have been compared to the one obtained by FD on the same (about 2,700) strides. Such a comparison aims to check the validity of the investigated models to capture the “true” dynamics of gaits in terms of distance between paths, mechanical external work and energy recovery. Results allow to conclude that: (1) among gaits, race walking is the most critical in being described by ID, (2) among the investigated segmental models, those capturing the motion of four limbs and trunk more closely reproduce the subtle temporal and spatial changes of BCoM trajectory within the strides of most gaits, (3) FD-ID discrepancy in external work is speed dependent within a gait in the most unsuccessful models, and (4) the internal work is not affected by the difference in BCoM estimates. PMID:28337148
Modeling Coastal Salinity in Quasi 2D and 3D Using a DUALEM-421 and Inversion Software.
Davies, Gareth; Huang, Jingyi; Monteiro Santos, Fernando Acacio; Triantafilis, John
2015-01-01
Rising sea levels, owing to climate change, are a threat to fresh water coastal aquifers. This is because saline intrusions are caused by increases and intensification of medium-large scale influences including sea level rise, wave climate, tidal cycles, and shifts in beach morphology. Methods are therefore required to understand the dynamics of these interactions. While traditional borehole and galvanic contact resistivity (GCR) techniques have been successful they are time-consuming. Alternatively, frequency-domain electromagnetic (FEM) induction is potentially useful as physical contact with the ground is not required. A DUALEM-421 and EM4Soil inversion software package are used to develop a quasi two- (2D) and quasi three-dimensional (3D) electromagnetic conductivity images (EMCI) across Long Reef Beach located north of Sydney Harbour, New South Wales, Australia. The quasi 2D models discern: the dry sand (<10 mS/m) associated with the incipient dune; sand with fresh water (10 to 20 mS/m); mixing of fresh and saline water (20 to 500 mS/m), and; saline sand of varying moisture (more than 500 mS/m). The quasi 3D EMCIs generated for low and high tides suggest that daily tidal cycles do not have a significant effect on local groundwater salinity. Instead, the saline intrusion is most likely influenced by medium-large scale drivers including local wave climate and morphology along this wave-dominated beach. Further research is required to elucidate the influence of spring-neap tidal cycles, contrasting beach morphological states and sea level rise.
Pavei, Gaspare; Seminati, Elena; Cazzola, Dario; Minetti, Alberto E
2017-01-01
The dynamics of body center of mass (BCoM) 3D trajectory during locomotion is crucial to the mechanical understanding of the different gaits. Forward Dynamics (FD) obtains BCoM motion from ground reaction forces while Inverse Dynamics (ID) estimates BCoM position and speed from motion capture of body segments. These two techniques are widely used by the literature on the estimation of BCoM. Despite the specific pros and cons of both methods, FD is less biased and considered as the golden standard, while ID estimates strongly depend on the segmental model adopted to schematically represent the moving body. In these experiments a single subject walked, ran, (uni- and bi-laterally) skipped, and race-walked at a wide range of speeds on a treadmill with force sensors underneath. In all conditions a simultaneous motion capture (8 cameras, 36 markers) took place. 3D BCoM trajectories computed according to five marker set models of ID have been compared to the one obtained by FD on the same (about 2,700) strides. Such a comparison aims to check the validity of the investigated models to capture the "true" dynamics of gaits in terms of distance between paths, mechanical external work and energy recovery. Results allow to conclude that: (1) among gaits, race walking is the most critical in being described by ID, (2) among the investigated segmental models, those capturing the motion of four limbs and trunk more closely reproduce the subtle temporal and spatial changes of BCoM trajectory within the strides of most gaits, (3) FD-ID discrepancy in external work is speed dependent within a gait in the most unsuccessful models, and (4) the internal work is not affected by the difference in BCoM estimates.
NASA Astrophysics Data System (ADS)
Trovato, C.; Lokmer, I.; De Martin, F.; Aochi, H.
2016-11-01
Since a few decades volcanic long period (LP) events have been recorded on many active volcanoes and their study has been recognized as an important tool to characterize volcanic activity. LP event analyses through moment tensor (MT) inversions have led to kinematic descriptions of various source mechanisms. The main challenge in these inversions is to `strip out' the propagation effect in order to isolate the source; hence the velocity model used controls the accuracy of the retrieved source mechanism. We first carry out several synthetic tests of inversions on Mt Etna volcano (Italy). Four geological models with topography are considered with increasing complexity: the most complex model is used to generate synthetic data, while the other three models are used to calculate the Green's functions for inversions. The retrieved solutions from the three velocity models are similar. The MT solutions for a deeper source are well retrieved, while a shallower source test suffers from high uncertainties and strong misinterpretation of the source orientation. The homogeneous model gives the lowest misfit value, but source location and mechanism decomposition are inaccurate. When a complex model different from the true one is used, a high misfit value and a wrong solution is obtained. We then incorporate our findings into the MT inversion of an LP event recorded on Mt Etna in 2008. We obtain very different solutions among the three models in terms of source location and mechanism decomposition. The overall shape of the retrieved source time functions are similar, but some amplitude differences arise, especially for the homogeneous model. Our work highlights the importance of including the unconsolidated surface materials in the computation of Green's functions especially when dealing with shallow sources.
NASA Astrophysics Data System (ADS)
Camacho, Antonio G.; Carmona, Enrique; García-Jerez, Antonio; Sánchez-Martos, Francisco; Prieto, Juan F.; Fernández, José; Luzón, Francisco
2015-11-01
This paper presents a gravimetric study (based on 382 gravimetric stations in an area about 32 km2) of a nearly flat basin: the Low Andarax valley. This alluvial basin, close to its river mouth, is located in the extreme south of the province of Almería and coincides with one of the existing depressions in the Betic Cordillera. The paper presents new methodological work to adapt a published inversion approach (GROWTH method) to the case of an alluvial valley (sedimentary stratification, with density increase downward). The adjusted 3D density model reveals several features in the topography of the discontinuity layers between the calcareous basement (2,700 kg/m3) and two sedimentary layers (2,400 and 2,250 kg/m3). We interpret several low density alignments as corresponding to SE faults striking about N140-145°E. Some detected basement elevations (such as the one, previously known by boreholes, in Viator village) are apparently connected with the fault pattern. The outcomes of this work are: (1) new gravimetric data, (2) new methodological options, and (3) the resulting structural conclusions.
NASA Astrophysics Data System (ADS)
Kaban, Mikhail K.; Stolk, Ward; Tesauro, Magdala; El Khrepy, Sami; Al-Arifi, Nassir; Beekman, Fred; Cloetingh, Sierd A. P. L.
2016-11-01
We construct a new-generation 3D density model of the upper mantle of Asia and its surrounding areas based on a joint interpretation of several data sets. A recent model of the crust combining nearly all available seismic data is employed to calculate the impact of the crust on the gravity anomalies and observed topography and to estimate the residual mantle anomalies and residual topography. These fields are jointly inverted to calculate the density variations in the lithosphere and upper mantle down to 325 km. As an initial approximation, we estimate density variations using a seismic tomography model. Seismic velocity variations are converted into temperatures and then to density variations based on mineral physics constraints. In the Occam-type inversion, we fit both the residual mantle gravity anomalies and residual topography by finding deviations to the initial model. The obtained corrections improve the resolution of the initial model and reflect important features of the mantle structure that are not well resolved by the seismic tomography. The most significant negative corrections of the upper mantle density, found in the Siberian and East European cratons, can be associated with depleted mantle material. The most pronounced positive density anomalies are found beneath the Tarim and South Caspian basins, Barents Sea, and Bay of Bengal. We attribute these anomalies to eclogites in the uppermost mantle, which have substantially affected the evolution of the basins. Furthermore, the obtained results provide evidence for the presence of eclogites in the oceanic subducting mantle lithosphere.
3D structural cartography based on magnetic and gravity data inversion - Case of South-West Algeria
NASA Astrophysics Data System (ADS)
Hichem, Boubekri; Mohamed, Hamoudi; Abderrahmane, Bendaoud; Ivan, Priezzhev; Karim, Allek
2015-12-01
This article presents the results of 3D aeromagnetic and gravity data inversion across the West African Craton (WAC) in South West Algeria. Although the used data have different origins and resolutions, the performed manual and automatic interpretation for each dataset shows a good correlation with some earlier geological studies of the region, major structural aspects of the locality, as well as other new structural features. Many curved faults parallel to the suture zone indicate the presence of terranes or the metacratonization of the WAC and a related fault network of great importance with NE-SW and NW-SE directions. The mega shear zones from north to south, which are visible at the surface in the Hoggar, are also observed along the Saharan Platform. The fact that these faults are observed since the Cambro-Ordovician in all crust (including the Saharan Basins) indicates that this area, which is situated on the border of the WAC, remained active during the entire period of time.
NASA Astrophysics Data System (ADS)
Petrov, P.; Newman, G. A.
2010-12-01
-Fourier domain we had developed 3D code for full-wave field simulation in the elastic media which take into account nonlinearity introduced by free-surface effects. Our approach is based on the velocity-stress formulation. In the contrast to conventional formulation we defined the material properties such as density and Lame constants not at nodal points but within cells. This second order finite differences method formulated in the cell-based grid, generate numerical solutions compatible with analytical ones within the range errors determinate by dispersion analysis. Our simulator will be embedded in an inversion scheme for joint seismic- electromagnetic imaging. It also offers possibilities for preconditioning the seismic wave propagation problems in the frequency domain. References. Shin, C. & Cha, Y. (2009), Waveform inversion in the Laplace-Fourier domain, Geophys. J. Int. 177(3), 1067- 1079. Shin, C. & Cha, Y. H. (2008), Waveform inversion in the Laplace domain, Geophys. J. Int. 173(3), 922-931. Commer, M. & Newman, G. (2008), New advances in three-dimensional controlled-source electromagnetic inversion, Geophys. J. Int. 172(2), 513-535. Newman, G. A., Commer, M. & Carazzone, J. J. (2010), Imaging CSEM data in the presence of electrical anisotropy, Geophysics, in press.
NASA Astrophysics Data System (ADS)
Kuvshinov, Alexey; Semenov, Alexey
2012-06-01
We present a novel frequency-domain inverse solution to recover the 3-D electrical conductivity distribution in the mantle. The solution is based on analysis of local C-responses. It exploits an iterative gradient-type method - limited-memory quasi-Newton method - for minimizing the penalty function consisting of data misfit and regularization terms. The integral equation code is used as a forward engine to calculate responses and data misfit gradients during inversion. An adjoint approach is implemented to compute misfit gradients efficiently. Further improvements in computational load come from parallelizing the scheme with respect to frequencies, and from setting the most time-consuming part of the forward calculations - calculation of Green's tensors - apart from the inversion loop. Convergence, performance, and accuracy of our 3-D inverse solution are demonstrated with a synthetic numerical example. A companion paper applies the strategy set forth here to real data.
Weng Cho Chew
2004-10-27
The project aim was the improvement, evaluation, and application of one dimensional (1D) inversion and development and application of three dimensional (3D) inversion to processing of data collected at waste pits at the Idaho National Engineering and Environmental Laboratory. The inversion methods were intended mainly for the Very Early Time Electromagnetic (VETEM) system which was designed to improve the state-of-the-art of electromagnetic imaging of the shallow (0 to about 5m) subsurface through electrically conductive soils.
Commer, Michael; Helwig, Stefan, L.; Hordt, Andreas; Scholl,Carsten; Tezkan, Bulent
2006-06-14
Three long-offset transient electromagnetic (LOTEM) surveyswerecarried out at the active volcano Merapi in Central Java (Indonesia)during the years 1998, 2000, and 2001. The measurements focused on thegeneral resistivity structure of the volcanic edifice at depths of 0.5-2km and the further investigation of a southside anomaly. The measurementswere insufficient for a full 3D inversion scheme, which could enable theimaging of finely discretized resistivity distributions. Therefore, astable, damped least-squares joint-inversion approach is used to optimize3D models with a limited number of parameters. The mode ls feature therealistic simulation of topography, a layered background structure, andadditional coarse 3D blocks representing conductivity anomalies.Twenty-eight LOTEM transients, comprising both horizontal and verticalcomponents of the magnetic induction time derivative, were analyzed. Inview of the few unknowns, we were able to achieve reasonable data fits.The inversion results indicate an upwelling conductor below the summit,suggesting hydrothermal activity in the central volcanic complex. Ashallow conductor due to a magma-filled chamber, at depths down to 1 kmbelow the summit, suggested by earlier seismic studies, is not indicatedby the inversion results. In conjunction with an anomalous-density model,derived from arecent gravity study, our inversion results provideinformation about the southern geological structure resulting from amajor sector collapse during the Middle Merapi period. The density modelallows to assess a porosity range andthus an estimated vertical salinityprofile to explain the high conductivities on a larger scale, extendingbeyond the foothills of Merapi.
NASA Astrophysics Data System (ADS)
Lebedeva-Ivanova, Nina; Gaina, Carmen; Minakov, Alexander; Kashubin, Sergey
2016-04-01
We derived Moho depth and crustal thickness for the High Arctic region by 3D forward and inverse gravity modelling method in the spectral domain (Minakov et al. 2012) using lithosphere thermal gravity anomaly correction (Alvey et al., 2008); a vertical density variation for the sedimentary layer and lateral crustal variation density. Recently updated grids of bathymetry (Jakobsson et al., 2012), gravity anomaly (Gaina et al, 2011) and dynamic topography (Spasojevic & Gurnis, 2012) were used as input data for the algorithm. TeMAr sedimentary thickness grid (Petrov et al., 2013) was modified according to the most recently published seismic data, and was re-gridded and utilized as input data. Other input parameters for the algorithm were calibrated using seismic crustal scale profiles. The results are numerically compared with publically available grids of the Moho depth and crustal thickness for the High Arctic region (CRUST 1 and GEMMA global grids; the deep Arctic Ocean grids by Glebovsky et al., 2013) and seismic crustal scale profiles. The global grids provide coarser resolution of 0.5-1.0 geographic degrees and not focused on the High Arctic region. Our grids better capture all main features of the region and show smaller error in relation to the seismic crustal profiles compare to CRUST 1 and GEMMA grids. Results of 3D gravity modelling by Glebovsky et al. (2013) with separated geostructures approach show also good fit with seismic profiles; however these grids cover the deep part of the Arctic Ocean only. Alvey A, Gaina C, Kusznir NJ, Torsvik TH (2008). Integrated crustal thickness mapping and plate recon-structions for the high Arctic. Earth Planet Sci Lett 274:310-321. Gaina C, Werner SC, Saltus R, Maus S (2011). Circum-Arctic mapping project: new magnetic and gravity anomaly maps of the Arctic. Geol Soc Lond Mem 35, 39-48. Glebovsky V.Yu., Astafurova E.G., Chernykh A.A., Korneva M.A., Kaminsky V.D., Poselov V.A. (2013). Thickness of the Earth's crust in the
NASA Astrophysics Data System (ADS)
Miller, Craig A.; Williams-Jones, Glyn; Fournier, Dominique; Witter, Jeff
2017-02-01
Active, large volume, silicic magma systems are potentially the most hazardous form of volcanism on Earth. Knowledge of the location, size, and physical properties of silicic magma reservoirs, is therefore important for providing context in which to accurately interpret monitoring data and make informed hazard assessments. Accordingly, we present the first geophysical image of the Laguna del Maule volcanic field magmatic system, using a novel 3D inversion of gravity data constrained by thermodynamic modelling. The joint analysis of gravity and thermodynamic data allows for a rich interpretation of the magma system, and highlights the importance of considering the full thermodynamic effects on melt density, when interpreting gravity models of active magmatic systems. We image a 30 km3, low density, volatile rich magma reservoir, at around 2 km depth, containing at least 85% melt, hosted within a broader 115 km3 body interpreted as wholly or partially crystallised (>70% crystal) cumulate mush. Our model suggests a magmatic system with shallow, crystal poor magma, overlying deeper, crystal rich magma. Even though a large density contrast (-600 kg/m3) with the surrounding crust exists, the lithostatic load is 50% greater than the magma buoyancy force, suggesting buoyancy alone is insufficient to trigger an eruption. The reservoir is adjacent to the inferred extension of the Troncoso fault and overlies the location of an intruding sill, driving present day deformation. The reservoir is in close proximity to the 2.0 km3 Nieblas (rln) eruption at 2-3 ka, which we calculate tapped approximately 7% of the magma reservoir. However, we suggest that the present day magma system is not large enough to have fed all post-glacial eruptions, and that the location, or size of the system may have migrated or varied over time, with each eruption tapping only a small aliquot of the available magma. The presence of a shallow reservoir of volatile rich, near liquidus magma, in close
Theoretical assessment of 3-D magnetotelluric method for oil and gas exploration: Synthetic examples
NASA Astrophysics Data System (ADS)
Zhang, Kun; Wei, Wenbo; Lu, Qingtian; Dong, Hao; Li, Yanqing
2014-07-01
In petroleum explorations, seismic reflection technique has been almost always the preferred method for its high exploration depth and resolution. However, with the development of three dimensional (3D) inversion and interpretation schemes, much potential has been shown in MT method dealing with complex geological structures as in oil and gas exploration. In this study, synthetic geophysical models of petroleum reservoir structures are modeled and utilized to demonstrate that feasibility of 3-D MT technique for hydrocarbon exploration. A series of typical reservoir structure models are constructed and used to generate synthetic MT and seismic data to test the capabilities of 2-D/3-D MT and 2-D seismic inversion techniques. According to the inversion comparison, in addition to correctly retrieve the original forward model, the 3-D MT method also has some advantages over the reflective seismology method, which suffered from the lack of reflection wave and multiple wave problems. With the presented 3-D high resolution MT inversion method, MT techniques should be employed as one of the first choices for petroleum explorations.
NASA Astrophysics Data System (ADS)
Morgan, Joanna; Warner, Michael; Arnoux, Gillean; Hooft, Emilie; Toomey, Douglas; VanderBeek, Brandon; Wilcock, William
2016-02-01
3-D full-waveform inversion (FWI) is an advanced seismic imaging technique that has been widely adopted by the oil and gas industry to obtain high-fidelity models of P-wave velocity that lead to improvements in migrated images of the reservoir. Most industrial applications of 3-D FWI model the acoustic wavefield, often account for the kinematic effect of anisotropy, and focus on matching the low-frequency component of the early arriving refractions that are most sensitive to P-wave velocity structure. Here, we have adopted the same approach in an application of 3-D acoustic, anisotropic FWI to an ocean-bottom-seismometer (OBS) field data set acquired across the Endeavour oceanic spreading centre in the northeastern Pacific. Starting models for P-wave velocity and anisotropy were obtained from traveltime tomography; during FWI, velocity is updated whereas anisotropy is kept fixed. We demonstrate that, for the Endeavour field data set, 3-D FWI is able to recover fine-scale velocity structure with a resolution that is 2-4 times better than conventional traveltime tomography. Quality assurance procedures have been employed to monitor each step of the workflow; these are time consuming but critical to the development of a successful inversion strategy. Finally, a suite of checkerboard tests has been performed which shows that the full potential resolution of FWI can be obtained if we acquire a 3-D survey with a slightly denser shot and receiver spacing than is usual for an academic experiment. We anticipate that this exciting development will encourage future seismic investigations of earth science targets that would benefit from the superior resolution offered by 3-D FWI.
NASA Astrophysics Data System (ADS)
Gao, J.; Zhang, H.
2015-12-01
Near surface geophysical exploration for the purpose of engineering design or construction For this reason, geophysical imaging demands a higher resolution and a better quantitative interpretation. Seismic travel time tomography and direct current resistivity tomography are two main methods for the near surface survey. Because of the limited coverage of observation system and the complex physical relationship between physical parameters and observations, individual geophysical method suffers issues of non-uniqueness and resolution limitation to some degree. We have developed a joint inversion method to combine seismic travel time tomography and full channel resistivity tomography. For the full channel resistivity survey, it uses two electrodes for power supply and all the other electrodes for recording. Compared with the traditional resistivity method, it collects more data and has a better model converge. Our joint inversion strategy relies on the structure constraint enforced through minimizing cross gradients between seismic velocity and resistivity models (Gallardo, 2003). For resistivity tomography, sensitivity kernels are obtained through the adjoint method by solving the electrostatic field equation with the finite-difference method. For seismic travel time tomography, ray paths and travel times are calculated using the fast marching method. We have tested our joint inversion method for a 2D cross-hole problem where two small zones with high and low velocity/resistivity anomalies. Seismic/electrical sources/receivers are installed in two boreholes. For separate seismic inversion, the smearing effect is evident and two anomaly zones are distorted and misplaced. For separate electric resistivity inversion, although two anomaly zones are positioned correctly their values are not accurate. By joint inversion, two velocity anomaly zones are clearly imaged and the smearing effect is greatly reduced. In comparison, for the resistivity model, the two anomaly zones
Design optimization of axial flow hydraulic turbine runner: Part I - an improved Q3D inverse method
NASA Astrophysics Data System (ADS)
Peng, Guoyi; Cao, Shuliang; Ishizuka, Masaru; Hayama, Shinji
2002-06-01
With the aim of constructing a comprehensive design optimization procedure of axial flow hydraulic turbine, an improved quasi-three-dimensional inverse method has been proposed from the viewpoint of system and a set of rotational flow governing equations as well as a blade geometry design equation has been derived. The computation domain is firstly taken from the inlet of guide vane to the far outlet of runner blade in the inverse method and flows in different regions are solved simultaneously. So the influence of wicket gate parameters on the runner blade design can be considered and the difficulty to define the flow condition at the runner blade inlet is surmounted. As a pre-computation of initial blade design on S2m surface is newly adopted, the iteration of S1 and S2m surfaces has been reduced greatly and the convergence of inverse computation has been improved. The present model has been applied to the inverse computation of a Kaplan turbine runner. Experimental results and the direct flow analysis have proved the validation of inverse computation. Numerical investigations show that a proper enlargement of guide vane distribution diameter is advantageous to improve the performance of axial hydraulic turbine runner. Copyright
NASA Astrophysics Data System (ADS)
McNulty, David; Geaney, Hugh; Carroll, Elaine; Garvey, Shane; Lonergan, Alex; O’Dwyer, Colm
2017-02-01
Engineering Co3O4 nanoparticles into highly ordered, 3D inverse opal (IO) structures is shown to significantly improve their performance as more efficient conversion mode Li-ion anode materials. By comparison with Co3O4 microparticles, the advantages of the porous anode architecture are clearly shown. The inverse opal material markedly enhances specific capacity and capacity retention. The impact of various C rates on the rate of the initial charge demonstrates that higher rate charging (10 C) was much less destructive to the inverse opal structure than charging at a slow rate (0.1 C). Slower C rates that affect the IO structure resulted in higher specific capacities (more Li2O) as well as improved capacity retention. The IO structures cycle as CoO, which improves Coulombic efficiency and limits volumetric changes, allowing rate changes more efficiently. This work demonstrates how 3D IOs improve conversion mode anode material performance in the absence of additive or binders, thus enhancing mass transport of Li2O charge–discharge product through the open structure. This effect mitigates clogging by structural changes at slow rates (high capacity) and is beneficial to the overall electrochemical performance.
Devavry, Séverine; Legros, Céline; Brasseur, Chantal; Delagrange, Philippe; Spadoni, Gilberto; Cohen, William; Malpaux, Benoît; Boutin, Jean A; Nosjean, Olivier
2012-08-01
Melatonin receptors have been described to activate different G protein-dependent signaling pathways, both in laboratory, heterologous, cellular models and in physiological conditions. Furthermore, the constitutive activity of G protein-coupled receptors has been shown to be key in physiological and pathological conditions. In the case of melatonin receptors, information is rather scare and concerns only MT1 receptors. In the present report, we show that the G protein-coupled melatonin receptors do have a constitutive, nonmelatonin-induced signaling activity using two cellular models of different origins, the Chinese hamster ovary cell line and Neuro2A, a neuroblastoma cell line. Furthermore, we show that this constitutive activity involves mainly Gi proteins, which is consistent with the common knowledge on the melatonin receptors. Importantly, we also describe, for the first time, inverse agonist properties for melatonin ligands. Although it is clear than more in-depth, biochemistry-based studies will be required to better understand by which pathway(s) the constitutively active melatonin receptors transfer melatonin information into intracellular biochemical events; our data open interesting perspectives for understanding the importance of the constitutive activity of melatonin receptors in physiological conditions.
Chajon, Enrique; Dumas, Isabelle; Touleimat, Mahmoud B.Sc.; Magne, Nicolas; Coulot, Jeremy; Verstraet, Rodolfe; Lefkopoulos, Dimitri; Haie-Meder, Christine
2007-11-01
Purpose: The purpose of this study was to evaluate the inverse planning simulated annealing (IPSA) software for the optimization of dose distribution in patients with cervix carcinoma treated with MRI-based pulsed-dose rate intracavitary brachytherapy. Methods and Materials: Thirty patients treated with a technique using a customized vaginal mold were selected. Dose-volume parameters obtained using the IPSA method were compared with the classic manual optimization method (MOM). Target volumes and organs at risk were delineated according to the Gynecological Brachytherapy Group/European Society for Therapeutic Radiology and Oncology recommendations. Because the pulsed dose rate program was based on clinical experience with low dose rate, dwell time values were required to be as homogeneous as possible. To achieve this goal, different modifications of the IPSA program were applied. Results: The first dose distribution calculated by the IPSA algorithm proposed a heterogeneous distribution of dwell time positions. The mean D90, D100, and V100 calculated with both methods did not differ significantly when the constraints were applied. For the bladder, doses calculated at the ICRU reference point derived from the MOM differed significantly from the doses calculated by the IPSA method (mean, 58.4 vs. 55 Gy respectively; p = 0.0001). For the rectum, the doses calculated at the ICRU reference point were also significantly lower with the IPSA method. Conclusions: The inverse planning method provided fast and automatic solutions for the optimization of dose distribution. However, the straightforward use of IPSA generated significant heterogeneity in dwell time values. Caution is therefore recommended in the use of inverse optimization tools with clinical relevance study of new dosimetric rules.
Solid-liquid phase transitions in 3D systems with the inverse-power and Yukawa potentials
NASA Astrophysics Data System (ADS)
Vaulina, O. S.; Koss, X. G.
2016-03-01
The melting of face-centered cubic (fcc) and body-centered cubic (bcc) crystal lattices was studied analytically and numerically for the systems of particles interacting via the inverse-power-law and Yukawa potentials. New approach is proposed for determination of the solid-liquid phase transitions in these systems. The suggested approach takes into account a nonlinearity (anharmonicity) of pair interaction forces and allows to correctly predict the conditions of melting of the systems with various isotropic pair interaction potentials. The obtained results are compared with the existing theoretical and numerical data.
NASA Astrophysics Data System (ADS)
Martin, Roland; Monteiller, Vadim; Chevrot, Sébastien; Wang, Yi; Komatitsch, Dimitri; Dufréchou, Grégory
2015-04-01
We describe here a method of inversion applied to seismic data sets constrained by gravity data at the regional scale. This will allow us to obtain robust models of P and S wave velocities but also of density, providing key constraints on the composition and thermal state of the lithosphere. Our approach relies on teleseimic waves, which illuminate the medium from below. We have developped a hybrid method in which a wave propagation method at the global scale (DSM/Direct solution method) is coupled with a spectral element method at the regional scale (Monteiller et al. 2013). With the spectral element method, we are able to model the 3D wave propagation effects in a computational domain of 400km long x 400km wide and 200 km deep, for an incident teleseismic wavefront introduced at the boundaries of this domain with periods as short as 2 s. The DSM global method allows to compute this incident field for a spherical Earth model. We use a multi-scale joint inversion of both gravity and seismic waveform data, accounting for the long wavelengths of the gravity field taken from a global model. In terms of inversion technique, we have validated an adjoint method for the inversion of seismic waveforms. An optimized BFGS inversion technique is used to minimize the difference between observed and computed full waveforms. The gradient of the misfit function gives the direction over which the model must be perturbed to minimize this difference. At each step of the inversion procedure we choose an optimal step length that accelerates the minimization. This is the crucial ingredient that allows us to build an efficient iterative full waveform inversion. We have extended this method by incorporating gravity data provided by the BGI/Bureau Gravimétrique International into the inversion. If the waveforms allow us to constrain the seismic velocities, they are less sensitive to the structure in density, which gives independent and crucial information to constrain the nature of rocks
Speidel, M; Hatt, C; Tomkowiak, M; Raval, A; Funk, T
2014-06-15
Purpose: To develop a method for the fusion of 3D echocardiography and Scanning-Beam Digital X-ray (SBDX) fluoroscopy to assist with catheter device and soft tissue visualization during interventional procedures. Methods: SBDX is a technology for low-dose inverse geometry x-ray fluoroscopy that performs digital tomosynthesis at multiple planes in real time. In this study, transesophageal echocardiography (TEE) images were fused with SBDX images by estimating the 3D position and orientation (the “pose”) of the TEE probe within the x-ray coordinate system and then spatially transforming the TEE image data to match this pose. An initial pose estimate was obtained through tomosynthesis-based 3D localization of points along the probe perimeter. Position and angle estimates were then iteratively refined by comparing simulated projections of a 3D probe model against SBDX x-ray images. Algorithm performance was quantified by imaging a TEE probe in different known orientations and locations within the x-ray field (0-30 degree tilt angle, up to 50 mm translation). Fused 3D TEE/SBDX imaging was demonstrated by imaging a tissue-mimicking polyvinyl alcohol cylindrical cavity as a catheter was navigated along the cavity axis. Results: Detected changes in probe tilt angle agreed with the known changes to within 1.2 degrees. For a 50 mm translation along the source-detector axis, the detected translation was 50.3 mm. Errors for in-plane translations ranged from 0.1 to 0.9 mm. In a fused 3D TEE/SBDX display, the catheter device was well visualized and coincident with the device shadow in the TEE images. The TEE images portrayed phantom boundaries that were not evident under x-ray. Conclusion: Registration of soft tissue anatomy derived from TEE imaging and device imaging from SBDX x-ray fluoroscopy is feasible. The simultaneous 3D visualization of these two modalities may be useful in interventional procedures involving the navigation of devices to soft tissue anatomy.
NASA Astrophysics Data System (ADS)
An, Zhiguo; Di, Qingyun
2016-12-01
The Alxa area in Inner Mongolia has been selected as a possible site for geological disposal of high-level radioactive waste (HLRW). Based on results of a previous study on crustal stability, the Tamusu rock mass has been chosen as the target. To determine the geological structure of this rock mass, aeromagnetic and gravity data are collected and inverted. Three-dimensional (3D) inversion horizontal slices show that the internal density of the rock mass and the distribution of magnetic properties are not uniform, with fractures and fragmentation being present. To confirm this result, the controlled source audio-frequency magnetotelluric method (CSAMT) was applied to explore the geological structures, the typical CSAMT sounding curve was analyzed, and the response characteristics of the geological structure and surrounding rock are distinguished. The original data were processed and interpreted in combination with data from surface geology and drilling and logging data. It is found that the CSAMT results were consistent with those from 3D inversion of the gravity and magnetic data, confirming the existence of fractures and fragmentation in the exploration area.
NASA Astrophysics Data System (ADS)
Uzkeda, Hodei; Bulnes, Mayte; Poblet, Josep; García-Ramos, José Carlos; Piñuela, Laura
2016-09-01
We constructed a geological map, a 3D model and cross-sections, carried out a structural analysis, determined the stress fields and tectonic transport vectors, restored a cross section and performed a subsidence analysis to unravel the kinematic evolution of the NE emerged portion of the Asturian Basin (NW Iberian Peninsula), where Jurassic rocks crop out. The major folds run NW-SE, normal faults exhibit three dominant orientations: NW-SE, NE-SW and E-W, and thrusts display E-W strikes. After Upper Triassic-Lower Jurassic thermal subsidence, Middle Jurassic doming occurred, accompanied by normal faulting, high heat flow and basin uplift, followed by Upper Jurassic high-rate basin subsidence. Another extensional event, possibly during Late Jurassic-Early Cretaceous, caused an increment in the normal faults displacement. A contractional event, probably of Cenozoic age, led to selective and irregularly distributed buttressing and fault reactivation as reverse or strike-slip faults, and folding and/or offset of some previous faults by new generation folds and thrusts. The Middle Jurassic event could be a precursor of the Bay of Biscay and North Atlantic opening that occurred from Late Jurassic to Early Cretaceous, whereas the Cenozoic event would be responsible for the Pyrenean and Cantabrian ranges and the partial closure of the Bay of Biscay.
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)
Ross, Z.; Ben-Zion, Y.; Zhu, L.; Graves, R. W.
2015-12-01
We perform a full source tensor inversion of several M > 4 earthquakes that occurred in the San Jacinto fault zone in southern California, with an emphasis on resolving signatures of volumetric source changes. A previous study on these events with Green's functions based on a 1D velocity model identified statistically significant explosive isotropic components (Ross et al. 2015). Here we use the SCEC 3D Community Velocity Model to derive Green's functions with source-receiver reciprocity and finite-difference calculations based on the code of Graves (1996). About 50 stations are used at epicentral distances of up to 55 km. The inversions are performed using the 'generalized Cut and Paste' method, which includes CLVD and isotropic components (Zhu and Ben-Zion 2013). The derived source tensors are compared to the results of the previous study based on the simplified 1D velocity model. The results are analyzed with bootstrap analysis to estimate uncertainties involved. Additional tests are performed using synthetic waveforms to study the effects of neglecting various features on the source inversions.
Wang, Aijun; Liu, Wenfang; Tang, Junjie; Chen, Sheng-Li; Dong, Peng
2014-04-15
A photonic bandgap (PBG) extension of surface-disordered 3D photonic crystals (PCs) based on the TiO2 inverse opal (TiO2-IO) architecture has been demonstrated. By using a liquid phase deposition (LPD) process based on the controlled hydrolysis of ammonium hexafluorotitanate and boric acid, an extra layer of TiO2 nanoparticles were deposited onto the internal surface of the air voids in the TiO2-IOs to increase their surface roughness, thereby introducing surface disorder in the 3D order structures. The PBG relative width of surface-disordered TiO2-IOs has been broadened significantly, and, compared to the original TiO2-IO, its largest rate of increase (27%) has been obtained. It was found that the PBG relative width increased rapidly at first and then to a much slower rate of change with increase of the duration of the LPD time. A possible cause for this finding is discussed in this Letter.
Oldenburg, Amy L.; Yu, Xiao; Gilliss, Thomas; Alabi, Oluwafemi; Taylor, Russell M.; Troester, Melissa A.
2015-01-01
The progression of breast cancer is known to be affected by stromal cells within the local microenvironment. Here we study the effect of stromal fibroblasts on the in-place motions (motility) of mammary epithelial cells within organoids in 3D co-culture, inferred from the speckle fluctuation spectrum using optical coherence tomography (OCT). In contrast to Brownian motion, mammary cell motions exhibit an inverse power-law fluctuation spectrum. We introduce two complementary metrics for quantifying fluctuation spectra: the power-law exponent and a novel definition of the motility amplitude, both of which are signal- and position-independent. We find that the power-law exponent and motility amplitude are positively (p<0.001) and negatively (p<0.01) correlated with the density of stromal cells in 3D co-culture, respectively. We also show how the hyperspectral data can be visualized using these metrics to observe heterogeneity within organoids. This constitutes a simple and powerful tool for detecting and imaging cellular functional changes with OCT. PMID:26973862
NASA Astrophysics Data System (ADS)
Mashayekhi, Parisa; Ghorbani-Dashtaki, Shoja; Mosaddeghi, Mohammad Reza; Shirani, Hossein; Nodoushan, Ali Reza Mohammadi
2016-04-01
In this study, HYDRUS-2D/3D was used to simulate ponded infiltration through double-ring infiltrometers into a hypothetical loamy soil profile. Twelve scenarios of inverse modelling (divided into three groups) were considered for estimation of Mualem-van Genuchten hydraulic parameters. In the first group, simulation was carried out solely using cumulative infiltration data. In the second group, cumulative infiltration data plus water content at h = -330 cm (field capacity) were used as inputs. In the third group, cumulative infiltration data plus water contents at h = -330 cm (field capacity) and h = -15 000 cm (permanent wilting point) were used simultaneously as predictors. The results showed that numerical inverse modelling of the double-ring infiltrometer data provided a reliable alternative method for determining soil hydraulic parameters. The results also indicated that by reducing the number of hydraulic parameters involved in the optimization process, the simulation error is reduced. The best one in infiltration simulation which parameters α, n, and Ks were optimized using the infiltration data and field capacity as inputs. Including field capacity as additional data was important for better optimization/definition of soil hydraulic functions, but using field capacity and permanent wilting point simultaneously as additional data increased the simulation error.
NASA Astrophysics Data System (ADS)
Martin, Roland; Chevrot, Sébastien; Komatitsch, Dimitri; Seoane, Lucia; Spangenberg, Hannah; Wang, Yi; Dufréchou, Grégory; Bonvalot, Sylvain; Bruinsma, Sean
2017-01-01
We image the internal density structure of the Pyrenees by inverting gravity data using an a priori density model derived by scaling a Vp model obtained by full waveform inversion of teleseismic P-waves. Gravity anomalies are computed via a 3D high-order finite-element integration in the same high-order spectral-element grid as the one used to solve the wave equation and thus to obtain the velocity model. The curvature of the Earth and surface topography are taken into account in order to obtain a density model as accurate as possible. The method is validated through comparisons with exact semi-analytical solutions. We show that the spectral element method drastically accelerates the computations when compared to other more classical methods. Different scaling relations between compressional velocity and density are tested, and the Nafe-Drake relation is the one that leads to the best agreement between computed and observed gravity anomalies. Gravity data inversion is then performed and the results allow us to put more constraints on the density structure of the shallow crust and on the deep architecture of the mountain range.
West Flank Coso FORGE Magnetotelluric 3D Data
Doug Blankenship
2016-01-01
This is the 3D version of the MT data for the West Flank FORGE area.The Coso geothermal field has had three Magnetotelluric (MT) datasets collected including surveys in 2003, 2006, and 2011. The final collection, in 2011, expanded the survey to the west and covers the West Flank of FORGE area.This most recent data set was collected by Schlumberger/WesternGeco and inverted by the WesternGeco GeoSolutions Integrated EM Center of Excellence in Milan, Italy; the 2003 and 2006 data were integrated for these inversions in the present study.
NASA Technical Reports Server (NTRS)
Pazmany, Andrew L.; Reehorst, Andrew (Technical Monitor)
2001-01-01
The Mount Washington Icing Sensors Project (MWISP) was a multi-investigator experiment with participants from Quadrant Engineering, NOAA Environmental Technology Laboratory (NOAA/ETL), the Microwave Remote Sensing Laboratory (MIRSL) of the University of Massachusetts (UMass), and others. Radar systems from UMass and NOAA/ETL were used to measure X-, Ka-, and W-band backscatter data from the base of Mt. Washington, while simultaneous in-situ particle measurements were made from aircraft and from the observatory at the summit. This report presents range and time profiles of liquid water content and particle size parameters derived from range profiles of radar reflectivity as measured at X-, Ka-, and W-band (9.3, 33.1, and 94.9 GHz) using an artificial neural network inversion algorithm. In this report, we provide a brief description of the experiment configuration, radar systems, and a review of the artificial neural network used to extract cloud parameters from the radar data. Time histories of liquid water content (LWC), mean volume diameter (MVD) and mean Z diameter (MZD) are plotted at 300 m range intervals for slant ranges between 1.1 and 4 km. Appendix A provides details on the extraction of radar reflectivity from measured radar power, and Appendix B provides summary logs of the weather conditions for each day in which we processed data.
NASA Astrophysics Data System (ADS)
Kardell, Dominik A.
The two end-member concept of mantle plume-driven versus far field stress-driven continental rifting anticipates high volumes of magma emplaced close to the rift-initiating plume, whereas relatively low magmatic volumes are predicted at large distances from the plume where the rifting is thought to be driven by far field stresses. We test this concept at the Guinea Plateau, which represents the last area of separation between Africa and South America, by investigating for rift-related volumes of magmatism using borehole, 3D seismic, and gravity data to run structural 3D inversions in two different data areas. Despite our interpretation of igneous rocks spanning large areas of continental shelf covered by the available seismic surveys, the calculated volumes in the Guinea Plateau barely match the magmatic volumes of other magma-poor margins and thus endorse the aforementioned concept. While the volcanic units on the shelf seem to be characterized more dominantly by horizontally deposited extrusive volcanic flows distributed over larger areas, numerous paleo-seamounts pierce complexly deformed pre and syn-rift sedimentary units on the slope. As non-uniqueness is an omnipresent issue when using potential field data to model geologic features, our method faced some challenges in the areas exhibiting complicated geology. In this situation less rigid constraints were applied in the modeling process. The misfit issues were successfully addressed by filtering the frequency content of the gravity data according to the depth of the investigated geology. In this work, we classify and compare our volume estimates for rift-related magmatism between the Guinea Fracture Zone (FZ) and the Saint Paul's FZ while presenting the refinements applied to our modeling technique.
NASA Astrophysics Data System (ADS)
Afonso, Juan Carlos; Rawlinson, Nicholas; Yang, Yingjie; Schutt, Derek L.; Jones, Alan G.; Fullea, Javier; Griffin, William L.
2016-10-01
We apply a novel 3-D multiobservable probabilistic tomography method that we have recently developed and benchmarked, to directly image the thermochemical structure of the Colorado Plateau and surrounding areas by jointly inverting P wave and S wave teleseismic arrival times, Rayleigh wave dispersion data, Bouguer anomalies, satellite-derived gravity gradients, geoid height, absolute (local and dynamic) elevation, and surface heat flow data. The temperature and compositional structures recovered by our inversion reveal a high level of correlation between recent basaltic magmatism and zones of high temperature and low Mg# (i.e., refertilized mantle) in the lithosphere, consistent with independent geochemical data. However, the lithospheric mantle is overall characterized by a highly heterogeneous thermochemical structure, with only some features correlating well with either Proterozoic and/or Cenozoic crustal structures. This suggests that most of the present-day deep lithospheric architecture reflects the superposition of numerous geodynamic events of different scale and nature to those that created major crustal structures. This is consistent with the complex lithosphere-asthenosphere system that we image, which exhibits a variety of multiscale feedback mechanisms (e.g., small-scale convection, magmatic intrusion, delamination, etc.) driving surface processes. Our results also suggest that most of the present-day elevation in the Colorado Plateau and surrounding regions is the result of thermochemical buoyancy sources within the lithosphere, with dynamic effects (from sublithospheric mantle flow) contributing only locally up to ˜15-35%.
NASA Astrophysics Data System (ADS)
Suzuki, W.; Aoi, S.; Maeda, T.; Sekiguchi, H.; Kunugi, T.
2013-12-01
Source inversion analysis using near-source strong-motion records with an assumption of 1-D underground structure models has revealed the overall characteristics of the rupture process of the 2011 Tohoku-Oki mega-thrust earthquake. This assumption for the structure model is acceptable because the seismic waves radiated during the Tohoku-Oki event were rich in the very-low-frequency contents lower than 0.05 Hz, which are less affected by the small-scale heterogeneous structure. The analysis using more reliable Green's functions even in the higher-frequency range considering complex structure of the subduction zone will illuminate more detailed rupture process in space and time and the transition of the frequency dependence of the wave radiation for the Tohoku-Oki earthquake. In this study, we calculate the near-source Green's functions using a 3-D underground structure model and perform the source inversion analysis using them. The 3-D underground structure model used in this study is the Japan Integrated Velocity Structure Model (Headquarters for Earthquake Research Promotion, 2012). A curved fault model on the Pacific plate interface is discretized into 287 subfaults at ~20 km interval. The Green's functions are calculated using GMS (Aoi et al., 2004), which is a simulation program package for the seismic wave field by the finite difference method using discontinuous grids (Aoi and Fujiwara, 1999). Computational region is 136-146.2E in longitude, 34-41.6N in latitude, and 0-100 km in depth. The horizontal and vertical grid intervals are 200 m and 100 m, respectively, for the shallower region and those for the deeper region are tripled. The number of the total grids is 2.1 billion. We derive 300-s records by calculating 36,000 steps with a time interval of 0.0083 second (120 Hz sampling). It takes nearly one hour to compute one case using 48 Graphics Processing Units (GPU) on TSUBAME2.0 supercomputer owned by Tokyo Institute of Technology. In total, 574 cases are
NASA Astrophysics Data System (ADS)
Macquet, Marie; Paul, Anne; Pedersen, Helle A.; Villaseñor, Antonio; Chevrot, Sébastien; Sylvander, Matthieu; Wolyniec, David; Pyrope Working Group
2014-10-01
The lithospheric architecture of the Pyrenees is still uncertain and highly debated. Here, we provide new constraints from a high-resolution 3-D S-wave velocity model of the Pyrenees and the adjacent foreland basins. This model is obtained from ambient noise tomography on records of temporary and permanent seismic arrays installed in southwestern France and northern Spain. We first computed group velocity maps for Rayleigh waves in the 5 to 55 s period range using noise correlation stacks at 1500-8500 station pairs. As the crust is very heterogeneous, poor results were obtained using a single starting model in a linearized inversion of group velocity dispersion curves for the shear wave structure. We therefore built a starting model for each grid node by full exploration of the model space. The resulting 3-D shear wave velocity model is compared to data from previous geophysical studies as a validation test. Despite the poor sensitivity of surface waves to seismic discontinuities, the geometry of the top of the basement and the Moho depth are retrieved well, except along the Cantabrian coast. Major reflectors of the ECORS deep seismic sounding profiles in the central and western Pyrenees coincide with sharp velocity gradients in our velocity model. We retrieve the difference between the thicker Iberian crust and the thinner European crust, the presence of low-velocity material of the Iberian crust underthrust beneath the European crust in the central Pyrenees, and the structural dissymmetry between the South Pyrenean Zone and the North Pyrenean Zone at the shallow crustal level. In the Labourd-Mauléon-Arzacq region (western Pyrenees), there is a high S-wave velocity anomaly at 20-30 km in depth, which might explain the positive Bouguer anomaly of the Labourd Massif. This high-velocity lower crust, which is also detected beneath the Parentis area, might be an imprint of the Albian-Aptian rifting phase. The southeastern part of the Massif Central has an unusual
NASA Astrophysics Data System (ADS)
Suzuki, Yuki; Kawai, Kenji; Geller, Robert J.; Borgeaud, Anselme F. E.; Konishi, Kensuke
2016-12-01
We conduct waveform inversion to infer the three-dimensional (3-D) S-velocity structure in the lowermost 400 km of the mantle (the D'' region) beneath the Northern Pacific region. Our dataset consists of about 20,000 transverse component broadband body-wave seismograms observed at North American stations for 131 intermediate and deep earthquakes which occurred beneath the western Pacific subduction region. We use S, ScS, and other phases that arrive between them. Resolution tests indicate that our methods and dataset can resolve the velocity structure in the target region with a horizontal scale of about 150 km and a vertical scale of about 50 km. The 3-D S-velocity model obtained in this study shows three prominent features: (1) prominent sheet-like lateral high-velocity anomalies up to ˜3% faster than the Preliminary Reference Earth Model (PREM) with a thickness of ˜200 km, whose lower boundary is ˜150 km above the core-mantle boundary (CMB). (2) A prominent low-velocity anomaly block located to the west of the Kamchatka peninsula, which is ˜2.5% slower than PREM, immediately above the CMB beneath the high-velocity anomalies. (3) A relatively thin (˜300 km) low-velocity structure continuous from the low-velocity anomaly "(2)" to at least 400 km above the CMB. We also detect a continuous low-velocity anomaly from the east of the Kamchatka peninsula at an altitude of 50 km above the CMB to the far east of the Kuril islands at an altitude of 400 km above the CMB. We interpret these features respectively as: (1) remnants of slab material where the bridgmanite to Mg-post-perovskite phase transition may have occurred within the slab, (2, 3) large amounts of hot and less dense materials beneath the cold Kula or Pacific slab remnants just above the CMB which ascend and form a passive plume upwelling at the edge of the slab remnants.[Figure not available: see fulltext.
NASA Astrophysics Data System (ADS)
Widiwijayanti, C.; Tiberi, C.; Diament, M.; Deplus, C.; Mikhailov, V.; Louat, R.; Tikhotsky, S.; Gvishiani, A.
2003-04-01
The Molucca Sea extending from northeastern Indonesia to southern Philippines islands, is a zone of oceanic basin closure between two opposite-facing subduction zones. This convergence results in the collision of two subduction zones, which style evolves from the southern to the northern parts of the Molucca Sea. In order to provide new insights into the present-day lithospheric structures in the Molucca Sea area, we inverted satellite and sea-surface gravity data into an iterative scheme including a priori seismological data. The seismological data were collected from two networks of Ocean Bottom Seismometer (OBS). These data allowed us to locate local seismic events and to build 3D tomographic images. We relate these results to the different stages of collision. The gravity data consists of combined sea-surface and satellite derived gravity. We used Kolmogorov-Wiener optimal (mean-square) filter to extract the gravity signal associated with lithospheric structures, then analyzed it to determine main regional features of lithospheric structure. For this purpose we employed a selection of Euler solutions based on a new clustering technique. To identify the geometry and nature of lithospheric structures, we also performed a 3-D gravity inversion for the northern Molucca Sea data, introducing our tomographic model as an independent constraint. The combination of both methods permits us to obtain a coherent image of the lithospheric structure. The results of this study illustrate the heterogeneity of lithospheric units in the northern Molucca Sea, which results from the collision between the Sangihe margin and lithospheric fragments from the Phillipine plate such as the Snellius plateau or the Halmahera volcanic arc. Three phenomena can explain the observed lithospheric structure: 1) the rupture of the Molucca Sea plate, accompanied by the appearance at the surface of slices of oceanic crust, favoring the development of suture zones as the collision evolves, 2) the
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 Astrophysics Data System (ADS)
Balvert, Marleen; van Hoof, Stefan J.; Granton, Patrick V.; Trani, Daniela; den Hertog, Dick; Hoffmann, Aswin L.; Verhaegen, Frank
2015-07-01
Advances in precision small animal radiotherapy hardware enable the delivery of increasingly complicated dose distributions on the millimeter scale. Manual creation and evaluation of treatment plans becomes difficult or even infeasible with an increasing number of degrees of freedom for dose delivery and available image data. The goal of this work is to develop an optimisation model that determines beam-on times for a given beam configuration, and to assess the feasibility and benefits of an automated treatment planning system for small animal radiotherapy. The developed model determines a Pareto optimal solution using operator-defined weights for a multiple-objective treatment planning problem. An interactive approach allows the planner to navigate towards, and to select the Pareto optimal treatment plan that yields the most preferred trade-off of the conflicting objectives. This model was evaluated using four small animal cases based on cone-beam computed tomography images. Resulting treatment plan quality was compared to the quality of manually optimised treatment plans using dose-volume histograms and metrics. Results show that the developed framework is well capable of optimising beam-on times for 3D dose distributions and offers several advantages over manual treatment plan optimisation. For all cases but the simple flank tumour case, a similar amount of time was needed for manual and automated beam-on time optimisation. In this time frame, manual optimisation generates a single treatment plan, while the inverse planning system yields a set of Pareto optimal solutions which provides quantitative insight on the sensitivity of conflicting objectives. Treatment planning automation decreases the dependence on operator experience and allows for the use of class solutions for similar treatment scenarios. This can shorten the time required for treatment planning and therefore increase animal throughput. In addition, this can improve treatment standardisation and
NASA Astrophysics Data System (ADS)
Tietze, Kristina; Ritter, Oliver
2013-10-01
3-D inversion techniques have become a widely used tool in magnetotelluric (MT) data interpretation. However, with real data sets, many of the controlling factors for the outcome of 3-D inversion are little explored, such as alignment of the coordinate system, handling and influence of data errors and model regularization. Here we present 3-D inversion results of 169 MT sites from the central San Andreas Fault in California. Previous extensive 2-D inversion and 3-D forward modelling of the data set revealed significant along-strike variation of the electrical conductivity structure. 3-D inversion can recover these features but only if the inversion parameters are tuned in accordance with the particularities of the data set. Based on synthetic 3-D data we explore the model space and test the impacts of a wide range of inversion settings. The tests showed that the recovery of a pronounced regional 2-D structure in inversion of the complete impedance tensor depends on the coordinate system. As interdependencies between data components are not considered in standard 3-D MT inversion codes, 2-D subsurface structures can vanish if data are not aligned with the regional strike direction. A priori models and data weighting, that is, how strongly individual components of the impedance tensor and/or vertical magnetic field transfer functions dominate the solution, are crucial controls for the outcome of 3-D inversion. If deviations from a prior model are heavily penalized, regularization is prone to result in erroneous and misleading 3-D inversion models, particularly in the presence of strong conductivity contrasts. A `good' overall rms misfit is often meaningless or misleading as a huge range of 3-D inversion results exist, all with similarly `acceptable' misfits but producing significantly differing images of the conductivity structures. Reliable and meaningful 3-D inversion models can only be recovered if data misfit is assessed systematically in the frequency
3-D Magnetotelluric studies of Pre-Cambrian basement beneath southern Alberta
NASA Astrophysics Data System (ADS)
Nieuwenhuis, G.; Unsworth, M.; Pana, D.; Craven, J.
2012-12-01
The Pre-Cambrian basement rocks beneath Alberta record the tectonic events that led to the assembly of Laurentia in the Proterozoic. Since these rocks are covered with younger sedimentary rocks, they must be investigated with geophysical methods. In the 1990s, these basement rocks were studied with a number of long-period magnetotelluric (MT) profiles collected by the Lithoprobe project. Dimensionality analysis of these data show that they appear to be two dimensional (2-D) in the period band 1-1000 s. However 2-D inversion models were unable to reproduce these MT data with a realistic resistivity model. The inversion models were very rough and characterized by many closely spaced conductors. Since the Lithoprobe data gave indications of 3-D resistivity structure, especially in the Archean Loverna block, additional MT data were collected by the University of Alberta from 2006-2010 using NIMS instruments. The goal was to develop an array that would constrain a fully 3-D model of crustal and upper mantle resistivity. The data at periods 1-10,000 s were inverted using a 3-D inversion algorithm. Comparisons between 2-D and 3-D inversions show that both models fit the measured MT data equally well. The 3-D model shows that the structure is dominated by an upper mantle conductor beneath the Loverna Block (the Loverna conductor). This conductor was previously imaged by the 2-D inversion of the Lithoprobe data. Our 3-D model shows that the Loverna conductor extends throughout the Archean Loverna block (part of the Hearne Domain) and is bounded to the south by a potential field anomaly known as the Vulcan Structure. Initial interpretations of the Vulcan Structure explained it as an intracontinental rift zone, while more recent studies show that it is more likely a north dipping subduction zone between two Archean blocks. This interpretation is supported by our 3-D resistivity model, which shows a good correlation between north dipping reflectors and the top of conductivity
NASA Astrophysics Data System (ADS)
Chen, Y. L.; Wu, J. C.; Guo, L. Y.; Wang, X. Y.; Tan, H. B.; Shen, C. Y.
2015-08-01
Conventional D-InSAR (Differential SAR Interferometry) can only monitor 1-D surface deformation along LOS (line of sight) direction. In order to overcome this limitation and extract 3-D coseismic displacement, we combine the LOS displacement derived from D-InSAR technology, the OKADA elastic half space dislocation model theory, jointly the surface rupture distribution by field investigations and the fault model inverted by GPS, level data and gravity survey to retrieve the directions of surface co-seismic displacement, and then have got Wenchuan Ms8.0 Earthquake 3D displacement. Firstly, thirty six L-band PALSAR images of six adjacent ascending tracks are processed with D-InSAR technology to obtain the coseismic displacements along LOS direction. According to the OKADA model and the thrust fault movement model of the Long-Men-Shan Fault , we specify the three directions of surface coseismic displacements. And thus the 3D coseismic displacement field is then recovered by using LOS displacement and relevant geometric projection formulas, obviously including horizontal displacements field and vertical deformation contour maps. By comparing with the 3D displacement estimated from OKADA dislocation model and fault model, the displacement retrieved in this study can give more detail, and reflect seismic characteristics more truly.
NASA Astrophysics Data System (ADS)
Kreeprasertkul, Kritti
Joint PP and SS seismic inversion is a geophysical analysis method that uses the variation in reflectivity of PP and SS seismic data to derive elastic properties of the subsurface and specifically reservoir properties. The major objective of this research is to characterize the elastic properties of the Bakken petroleum system at Banner Field, Mountrail County, North Dakota using joint PP and SS seismic inversion. The combination of PP and SS data into the inversion process helped improve the accuracy of inversion results by increasing the correlation coefficient by 10 percent on both P- and S- impedance. Results from this analysis show a variation in elastic properties from the western part to the eastern portion of the study area. The western part is interpreted to be higher in shale volume, porosity and maturity. Additionally, the natural fracture network is an integral part of a tight oil reservoir. The existence and delineation of natural fractures coincides with faults that have been delineated within the reservoir interval. The elastic parameters and structural detailing from this study has been used to delineate a potential exploration opportunity in the Three Forks Formation on the west side of the survey.
ERIC Educational Resources Information Center
Brown, Malcolm
2009-01-01
Inversions are fascinating phenomena. They are reversals of the normal or expected order. They occur across a wide variety of contexts. What do inversions have to do with learning spaces? The author suggests that they are a useful metaphor for the process that is unfolding in higher education with respect to education. On the basis of…
Rutz, Tobias; Piccini, Davide; Coppo, Simone; Chaptinel, Jerome; Ginami, Giulia; Vincenti, Gabriella; Stuber, Matthias; Schwitter, Juerg
2016-12-01
The border zone of post-infarction myocardial scar as identified by late gadolinium enhancement (LGE) has been identified as a substrate for arrhythmias and consequently, high-resolution 3D scar information is potentially useful for planning of electrophysiological interventions. This study evaluates the performance of a novel high-resolution 3D self-navigated free-breathing inversion recovery magnetic resonance pulse sequence (3D-SN-LGE) vs. conventional 2D breath-hold LGE (2D-LGE) with regard to sharpness of borders (SBorder) of post-infarction scar. Patients with post-infarction scar underwent two magnetic resonance examinations for conventional 2D-LGE and high-resolution 3D-SN-LGE acquisitions (both 15 min after 0.2 mmol/kg Gadobutrol IV) at 1.5T. In the prototype 3D-SN-LGE sequence, each ECG-triggered radial steady-state-free-precession read-out segment is preceded by a non-slice-selective inversion pulse. Scar volume and SBorder were assessed on 2D-LGE and matching reconstructed high-resolution 3D-SN-LGE short-axis slices. In 16 patients (four females, 58 ± 10y) all scars visualized by 2D-LGE could be identified on 3D-SN-LGE (time between 2D-LGE and 3D-SN-LGE 48 ± 53 days). A good agreement of scar volume by 3D-SN-LGE vs. 2D-LGE was found (Bland-Altman: -3.7 ± 3.4 ml, correlation: r = 0.987, p < 0.001) with a small difference in scar volume (20.5 (15.8, 35.2) ml vs. 24.5 (20.0, 41.9)) ml, respectively, p = 0.002] and a good intra- and interobserver variability (1.1 ± 4.1 and -1.1 ± 11.9 ml, respectively). SBorder of border "scar to non-infarcted myocardium" was superior on 3D-SN-LGE vs. 2D-LGE: 0.180 ± 0.044 vs. 0.083 ± 0.038, p < 0.001. Detection and quantification of myocardial scar by 3D-SN-LGE is feasible and accurate in comparison to 2D-LGE. The high spatial resolution of the 3D sequence improves delineation of scar borders.
NASA Astrophysics Data System (ADS)
Marinin, I. V.; Kabanikhin, S. I.; Krivorotko, O. I.; Karas, A.; Khidasheli, D. G.
2012-04-01
We consider new techniques and methods for earthquake and tsunami related problems, particularly - inverse problems for the determination of tsunami source parameters, numerical simulation of long wave propagation in soil and water and tsunami risk estimations. In addition, we will touch upon the issue of database management and destruction scenario visualization. New approaches and strategies, as well as mathematical tools and software are to be shown. The long joint investigations by researchers of the Institute of Mathematical Geophysics and Computational Mathematics SB RAS and specialists from WAPMERR and Informap have produced special theoretical approaches, numerical methods, and software tsunami and earthquake modeling (modeling of propagation and run-up of tsunami waves on coastal areas), visualization, risk estimation of tsunami, and earthquakes. Algorithms are developed for the operational definition of the origin and forms of the tsunami source. The system TSS numerically simulates the source of tsunami and/or earthquakes and includes the possibility to solve the direct and the inverse problem. It becomes possible to involve advanced mathematical results to improve models and to increase the resolution of inverse problems. Via TSS one can construct maps of risks, the online scenario of disasters, estimation of potential damage to buildings and roads. One of the main tools for the numerical modeling is the finite volume method (FVM), which allows us to achieve stability with respect to possible input errors, as well as to achieve optimum computing speed. Our approach to the inverse problem of tsunami and earthquake determination is based on recent theoretical results concerning the Dirichlet problem for the wave equation. This problem is intrinsically ill-posed. We use the optimization approach to solve this problem and SVD-analysis to estimate the degree of ill-posedness and to find the quasi-solution. The software system we developed is intended to
NASA Astrophysics Data System (ADS)
Zhou, Li; Zhang, Wei; Shen, Yang; Chen, Xiaofei; Zhang, Jie
2016-06-01
With dense seismic arrays and advanced imaging methods, regional three-dimensional (3D) Earth models have become more accurate. It is now increasingly feasible and advantageous to use a 3D Earth model to better locate earthquakes and invert their source mechanisms by fitting synthetics to observed waveforms. In this study, we develop an approach to determine both the earthquake location and source mechanism from waveform information. The observed waveforms are filtered in different frequency bands and separated into windows for the individual phases. Instead of picking the arrival times, the traveltime differences are measured by cross-correlation between synthetic waveforms based on the 3D Earth model and observed waveforms. The earthquake location is determined by minimizing the cross-correlation traveltime differences. We then fix the horizontal location of the earthquake and perform a grid search in depth to determine the source mechanism at each point by fitting the synthetic and observed waveforms. This new method is verified by a synthetic test with noise added to the synthetic waveforms and a realistic station distribution. We apply this method to a series of M W3.4-5.6 earthquakes in the Longmenshan fault (LMSF) zone, a region with rugged topography between the eastern margin of the Tibetan plateau and the western part of the Sichuan basin. The results show that our solutions result in improved waveform fits compared to the source parameters from the catalogs we used and the location can be better constrained than the amplitude-only approach. Furthermore, the source solutions with realistic topography provide a better fit to the observed waveforms than those without the topography, indicating the need to take the topography into account in regions with rugged topography.
2008-09-01
nonunique properties of inversion methods, we may often find a solution for one data type, but we must acknowledge that, although it can predict behavior...of density prisms and a suite of 1D fundamental mode group velocities. (a) A single cell with its input geographic coordinate system . (b) For a...H. K., H. Kanamori, P. C. Jennings, and C. Kissling (Eds.) (2002). International Handbook of Earthquake and Engineering Seismology (CD-ROM
NASA Astrophysics Data System (ADS)
Pilz, Marco; Parolai, Stefano; Woith, Heiko
2017-01-01