Joint inversion of ground gravity data and satellite gravity gradients between Nepal and Bhutan: New insights on structural and seismic segmentation of the Himalayan arc

Author(s):  
Rodolphe Cattin ◽  
Théo Berthet ◽  
György Hetényi ◽  
Anita Saraswati ◽  
Isabelle Panet ◽  
...  
2020 ◽  
Author(s):  
Yu Tian ◽  
Yong Wang

<p>The North China Craton (NCC) is one of the oldest cratons in the world. Currently, the destruction mechanism and geodynamics of the NCC still remain controversial. All of the proposed views regarding the issues involve studying the internal density structure of the NCC lithosphere. Gravity field data are one of the most important data in regard to investigating the lithospheric density structure, the gravity gradient data and the gravity data possess their own advantages. Given the inconsistency of the on orbit GOCE satellite gravity gradient and surface gravity observation plane height, also effects of the initial density model upon of the inversion results, the joint inversion of gravity gradient and gravity are divided into two integrated processes. By using the preconditioned conjugate gradient (PCG) inversion algorithm, the density data are calculated using the preprocessed remaining gravity anomaly data. The newly obtained high resolution density data are then used as the initial density model, which can be served as the constraints for the subsequent gravity gradient inversion. Downward continuation, terrain correction, interface undulation correction and long wavelength correction are performed for the four gravity gradient tensor data(<strong>T</strong><sub>xx</sub>,<strong>T</strong><sub>xz</sub>,<strong>T</strong><sub>yy</sub>,<strong>T</strong><sub>zz</sub>)of the Gravity field and steady-state Ocean Circulation Explorer (GOCE) satellite,  after which the remaining gravity gradient anomaly data(<strong>T</strong>'<sub>xx</sub>,<strong>T</strong>'<sub>xz</sub>,<strong>T</strong>'<sub>yy</sub>,<strong>T</strong>'<sub>zz</sub>) are used as the new observation quantity. Finally, the ultimate lithospheric density distribution within the depth range of 0–180 km in the NCC is obtained using the same PCG algorithm.</p>


2021 ◽  
Author(s):  
Wolfgang Szwillus

<p>Commonly, the physical properties of the Earth (e.g., velocity, density) are parameterized as continuous fields. The most popular representation are grids and basis functions like spherical harmonics or splines. In an inversion context it is quite common that not all the parameters are fully constrained by the available inputdata. This relates to the common issues of insufficient resolution, incomplete coverage, and trade-offs due tonon-uniqueness. By applying some form of regularization to the inverse problem, a well-behaved and unique solution can be obtained, but this solution depends on the details of the chosen regularization.</p><p>Transdimensional approaches address the regularization problem by using a model representation with a variable number of parameters. The number of parameters is adjusted according to the requirements of the input data using the reversible jump Monte Carlo Markov Chain (rj-MCMC) algorithm. The output is an ensemble of variable resolution models that provides insight into the required model complexity and trade-offbetween parameters.</p><p>Here, I present synthetic tests from a joint inversion of satellite gravity gradients and normal modes for the Earth's velocity and density structure. The mantle's seismic velocity and density inside a 2-D spherical annulus are described by a variable number of discrete anomalous volumes, each with a variable size, shape, location and strength of velocity and density anomaly. The discrete anomalies are adjusted using the transdimensional approach in order to fit the gravity and normal mode data. This synthetic example shows promising results, because the synthetic model can recovered reasonably well.</p>


2020 ◽  
Author(s):  
Bernhard Weise ◽  
Max Moorkamp ◽  
Stewart Fishwick

<p>The EarthScope USArray project provides high quality magnetotelluric and seismic observations, which have been used to identify tectonic boundaries of the USA. Combining these data sets together with satellite gravity observations, we investigate how the different data sets can complement each other in order to find a consistent model of the subsurface. Using a cross-gradient constraint, we first invert the magnetotelluric and gravity data sets in order to demonstrate the feasibility of our approach and to identify any difficulties. Once a joint conductivity and density model is found, we perform a full joint inversion of all three data sets. By comparison with models derived from separate inversions of the individual observables we can show how the different data sets interact. Examining the magnitude of the cross-gradient lets us distinguish parts of the model where a good agreement of the recovered structures has been achieved from those where differing patterns are necessary in order to achieve an acceptable data fit. In this presentation we will give an overview of our approach, highlight our strategy and show results from individual and joint inversions.</p>


2020 ◽  
Author(s):  
Yu Tian ◽  
Yong Wang

Abstract. The North China Craton (NCC) is one of the oldest cratons in the world. Currently, the destruction mechanism and geodynamics of the NCC still remain controversial. All of the proposed views regarding the issues involve studying the internal density structure of the NCC lithosphere. Gravity field data are one of the most important data in regard to investigating the lithospheric density structure, the gravity gradient data and the gravity data possess their own advantages. Given the inconsistency of the on orbit GOCE satellite gravity gradient and surface gravity observation plane height, also effects of the initial density model upon of the inversion results, the joint inversion of gravity gradient and gravity are divided into two integrated processes. By using the preconditioned conjugate gradient (PCG) inversion algorithm, the density data are calculated using the preprocessed remaining gravity anomaly data. The newly obtained high resolution density data are then used as the initial density model, which can be served as the constraints for the subsequent gravity gradient inversion. Several essential corrections are performed for the four gravity gradient tensor (Txx,Txz,Tyy,Tzz) of the GOCE satellite, after which the remaining gravity gradient anomaly (T'xx,


2021 ◽  
Vol 3 (5) ◽  
Author(s):  
Dostdar Hussain ◽  
Aftab Ahmed Khan ◽  
Syed Najam Ul Hassan ◽  
Syed Ali Asad Naqvi ◽  
Akhtar Jamil

AbstractMountains regions like Gilgit-Baltistan (GB) province of Pakistan are solely dependent on seasonal snow and glacier melt. In Indus basin which forms in GB, there is a need to manage water in a sustainable way for the livelihood and economic activities of the downstream population. It is important to monitor water resources that include glaciers, snow-covered area, lakes, etc., besides traditional hydrological (point-based measurements by using the gauging station) and remote sensing-based studies (traditional satellite-based observations provide terrestrial water storage (TWS) change within few centimeters from the earth’s surface); the TWS anomalies (TWSA) for the GB region are not investigated. In this study, the TWSA in GB region is considered for the period of 13 years (from January 2003 to December 2016). Gravity Recovery and Climate Experiment (GRACE) level 2 monthly data from three processing centers, namely Centre for Space Research (CSR), German Research Center for Geosciences (GFZ), and Jet Propulsion Laboratory (JPL), System Global Land Data Assimilation System (GLDAS)-driven Noah model, and in situ precipitation data from weather stations, were used for the study investigation. GRACE can help to forecast the possible trends of increasing or decreasing TWS with high accuracy as compared to the past studies, which do not use satellite gravity data. Our results indicate that TWS shows a decreasing trend estimated by GRACE (CSR, GFZ, and JPL) and GLDAS-Noah model, but the trend is not significant statistically. The annual amplitude of GLDAS-Noah is greater than GRACE signal. Mean monthly analysis of TWSA indicates that TWS reaches its maximum in April, while it reaches its minimum in October. Furthermore, Spearman’s rank correlation is determined between GRACE estimated TWS with precipitation, soil moisture (SM) and snow water equivalent (SWE). We also assess the factors, SM and SWE which are the most efficient parameters producing GRACE TWS signal in the study area. In future, our results with the support of more in situ data can be helpful for conservation of natural resources and to manage flood hazards, droughts, and water distribution for the mountain regions.


2020 ◽  
Vol 221 (3) ◽  
pp. 1542-1554 ◽  
Author(s):  
B C Root

SUMMARY Current seismic tomography models show a complex environment underneath the crust, corroborated by high-precision satellite gravity observations. Both data sets are used to independently explore the density structure of the upper mantle. However, combining these two data sets proves to be challenging. The gravity-data has an inherent insensitivity in the radial direction and seismic tomography has a heterogeneous data acquisition, resulting in smoothed tomography models with de-correlation between different models for the mid-to-small wavelength features. Therefore, this study aims to assess and quantify the effect of regularization on a seismic tomography model by exploiting the high lateral sensitivity of gravity data. Seismic tomography models, SL2013sv, SAVANI, SMEAN2 and S40RTS are compared to a gravity-based density model of the upper mantle. In order to obtain similar density solutions compared to the seismic-derived models, the gravity-based model needs to be smoothed with a Gaussian filter. Different smoothening characteristics are observed for the variety of seismic tomography models, relating to the regularization approach in the inversions. Various S40RTS models with similar seismic data but different regularization settings show that the smoothening effect is stronger with increasing regularization. The type of regularization has a dominant effect on the final tomography solution. To reduce the effect of regularization on the tomography models, an enhancement procedure is proposed. This enhancement should be performed within the spectral domain of the actual resolution of the seismic tomography model. The enhanced seismic tomography models show improved spatial correlation with each other and with the gravity-based model. The variation of the density anomalies have similar peak-to-peak magnitudes and clear correlation to geological structures. The resolvement of the spectral misalignment between tomographic models and gravity-based solutions is the first step in the improvement of multidata inversion studies of the upper mantle and benefit from the advantages in both data sets.


2017 ◽  
Vol 7 (1) ◽  
Author(s):  
K. Breili ◽  
R. Hougen ◽  
D. I. Lysaker ◽  
O. C. D. Omang ◽  
B. Tangen

AbstractThe Norwegian Mapping Authority (NMA) has recently established a new gravity laboratory in Ny-Ålesund at Svalbard, Norway. The laboratory consists of three independent pillars and is part of the geodetic core station that is presently under construction at Brandal, approximately 1.5 km north of NMA’s old station. In anticipation of future use of the new gravity laboratory, we present benchmark gravity values, gravity gradients, and final coordinates of all new pillars. Test measurements indicate a higher noise level at Brandal compared to the old station. The increased noise level is attributed to higher sensitivity to wind.We have also investigated possible consequences of moving to Brandal when it comes to the gravitational signal of present-day ice mass changes and ocean tide loading. Plausible models representing ice mass changes at the Svalbard archipelago indicate that the gravitational signal at Brandal may differ from that at the old site with a size detectable with modern gravimeters. Users of gravity data from Ny-Ålesund should, therefore, be cautious if future observations from the new observatory are used to extend the existing gravity record. Due to its lower elevation, Brandal is significantly less sensitive to gravitational ocean tide loading. In the future, Brandal will be the prime site for gravimetry in Ny-Ålesund. This ensures gravity measurements collocated with space geodetic techniques like VLBI, SLR, and GNSS.


2012 ◽  
Vol 18 (4) ◽  
pp. 549-563 ◽  
Author(s):  
Vagner G. Ferreira ◽  
Zheng Gong ◽  
Samuel A. Andam-Akorful

GRACE satellite gravity data was used to estimate mass changes within the Volta River basin in West African for the period of January, 2005 to December, 2010. We also used the precipitation data from the Tropical Rainfall Measurement Mission (TRMM) to determine relative contributions source to the seasonal hydrological balance within the Volta River basin. We found out that the seasonal mass change tends to be detected by GRACE for periods from 1 month in the south to 4 months in the north of the basin after the rainfall events. The results suggested a significant gain in water storage in the basin at reference epoch 2007.5 and a dominant annual cycle for the period under consideration for both in the mass changes and rainfall time series. However, there was a low correlation between mass changes and rainfall implying that there must be other processes which cause mass changes without rainfall in the upstream of the Volta River basin.


Sign in / Sign up

Export Citation Format

Share Document