scholarly journals LOCAL EVALUATION OF EARTH GRAVITATIONAL MODELS, CASE STUDY: IRAN

2017 ◽  
Vol 43 (1) ◽  
pp. 1-13 ◽  
Author(s):  
Ismael FOROUGHI ◽  
Yosra AFRASTEH ◽  
Sabah RAMOUZ ◽  
Abdolreza SAFARI

Global gravity models are being developed according to new data sets available from satellite gravity missions and terrestrial/marine gravity data which are provided by different countries. Some countries do not provide all their available data and the global gravity models have many vague computational methods. Therefore, the models need to be evaluated locally before using. It is generally understood that the accuracy of global gravity models is enough for local (civil, mining, construction, etc.) projects, however, our results in Iran show that the differences between synthesized values and observation data reach up to ∼300 mGal for gravity anomalies and ∼2 m for geoid heights. Even by applying the residual topographical correction to synthetized gravity anomalies, the differences are still notable. The accuracy of global gravity models for predicting marine gravity anomalies is also investigated in Persian Gulf and the results show differences of ∼140 mGal in coastal areas. The results of evaluating selected global gravity models in Iran indicate that the EIGEN-6C4 achieves the lowest RMS for estimating the geoid heights. EGM08 predicts the closest results to terrestrial gravity anomalies. DIR-R5 GOCE satellite-only model estimates the low-frequency part of gravity field more accurately. The best prediction of marine gravity anomalies is also achieved by EGM08.

Geophysics ◽  
1999 ◽  
Vol 64 (6) ◽  
pp. 1748-1755 ◽  
Author(s):  
Mara M. Yale ◽  
D. T. Sandwell

Gravity field recovery from satellite altimetry provides global marine coverage but lacks the accuracy and resolution needed for many exploration geophysics studies. The repeating ground tracks of the ERS-1/2, Geosat, and Topex/Poseidon altimeters offer the possibility of improving the accuracy and resolution of gravity anomalies along widely spaced (∼40-km spacing) tracks. However, complete ocean coverage is usually needed to convert the sea‐surface height (or along‐track slope) measurements into gravity anomalies. Here we develop and test a method for constructing stacked gravity profiles by using a published global gravity grid (Sandwell and Smith, 1997), V7.2, as a reference model for the slope‐to‐gravity anomaly conversion. The method is applied to stacks (averages) of Geosat/ERM (up to 62 cycles), ERS-1/2 (up to 43 cycles), and Topex (up to 142 cycles) satellite altimeter profiles. We assess the accuracies of the ERS-1/2 profiles through a comparison with a gravity model of the northern Gulf of Mexico (profiles provided by EDCON Inc.). The 40 ERS profiles evaluated have a mean rms difference of 3.77 mGal and full wavelength resolution (0.5 coherence) of 24 km. Our processing retains wavelengths as short as 10 km so smaller, large‐amplitude features can be resolved, especially in shallow ocean areas (<1000 m deep). We provide an example of combining these higher resolution profiles with lower resolution gravity data in the Caspian Sea.


2021 ◽  
Vol 9 ◽  
Author(s):  
Xiaoyun Wan ◽  
Weipeng Han ◽  
Jiangjun Ran ◽  
Wenjie Ma ◽  
Richard Fiifi Annan ◽  
...  

Marine gravity data from altimetry satellites are often used to derive bathymetry; however, the seafloor density contrast must be known. Therefore, if the ocean water depths are known, the density contrast can be derived. This study experimented the total least squares algorithm to derive seafloor density contrast using satellite derived gravity and shipborne depth observations. Numerical tests are conducted in a local area of the Atlantic Ocean, i.e., 34°∼32°W, 3.5°∼4.5°N, and the derived results are compared with CRUST1.0 values. The results show that large differences exist if the gravity and shipborne depth data are used directly, with mean difference exceeding 0.4 g/cm3. However, with a band-pass filtering applied to the gravity and shipborne depths to ensure a high correlation between the two data sets, the differences between the derived results and those of CRUST1.0 are reduced largely and the mean difference is smaller than 0.12 g/cm3. Since the spatial resolution of CRUST1.0 is not high and in many ocean areas the shipborne depths and gravity anomalies are much denser, the method of this study can be an alternative method for providing seafloor density variation information.


2021 ◽  
Vol 56 (3) ◽  
pp. 78-100
Author(s):  
Eyasu Alemu

Abstract To estimate Moho depth, geoid, gravity anomaly, and other geopotential functionals, gravity data is needed. But, gravity survey was not collected in equal distribution in Ethiopia, as the data forming part of the survey were mainly collected on accessible roads. To determine accurate Moho depth using Global Gravity Models (GGMs) for the study area, evaluation of GGMs is needed based on the available terrestrial gravity data. Moho depth lies between 28 km and 32 km in Afar. Gravity disturbances (GDs) were calculated for the terrestrial gravity data and the recent GGMs for the study area. The model-based GDs were compared with the corresponding GD obtained from the terrestrial gravity data and their differences in terms of statistical comparison parameters for determining the best fit GGM at a local scale in Afar. The largest standard deviation (SD) (36.10 mGal) and root mean square error (RMSE) (39.00 mGal) for residual GD and the lowest correlation with the terrestrial gravity (0.61 mGal) were obtained by the satellite-only model (GO_CONS_GCF_2_DIR_R6). The next largest SD (21.27 mGal) and RMSE (25.65 mGal) for residual GD were obtained by the combined gravity model (XGM2019e_2159), which indicates that it is not the best fit model for the study area as compared with the other two GGMs. In general, the result showed that the combined models are more useful tools for modeling the gravity field in Afar than the satellite-only GGMs. But, the study clearly revealed that for the study area, the best model in comparison with the others is the EGM2008, while the second best model is the EIGEN6C4.


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.


1999 ◽  
Vol 36 (1) ◽  
pp. 75-89 ◽  
Author(s):  
Hamid Telmat ◽  
Jean-Claude Mareschal ◽  
Clément Gariépy

Gravity data were obtained along two transects on the southern coast of Ungava Bay, which provide continuous gravity coverage between Leaf Bay and George River. The transects and the derived gravity profiles extend from the Superior craton to the Rae Province across the New Quebec Orogen (NQO). Interpretation of the transect along the southwestern coast of Ungava Bay suggests crustal thickening beneath the NQO and crustal thinning beneath the Kuujjuaq Terrane, east of the NQO. Two alternative interpretations are proposed for the transect along the southeastern coast of the bay. The first model shows crustal thickening beneath the George River Shear Zone (GRSZ) and two shallow bodies correlated with the northern extensions of the GRSZ and the De Pas batholith. The second model shows constant crustal thickness and bodies more deeply rooted than in the first model. The gravity models are consistent with the easterly dipping reflections imaged along a Lithoprobe seismic line crossing Ungava Bay and suggest westward thrusting of the Rae Province over the NQO. Because no gravity data have been collected in Ungava Bay, satellite altimetry data have been used as a means to fill the gap in data collected at sea. The satellite-derived gravity data and standard Bouguer gravity data were combined in a composite map for the Ungava Bay region. The new land-based gravity measurements were used to verify and calibrate the satellite data and to ensure that offshore gravity anomalies merge with those determined by the land surveys in a reasonable fashion. Three parallel east-west gravity profiles were extracted: across Ungava Bay (59.9°N), on the southern shore of the bay (58.5°N), and onshore ~200 km south of Ungava Bay (57.1°N). The gravity signature of some major structures, such as the GRSZ, can be identified on each profile.


2013 ◽  
Vol 7 (4) ◽  
pp. 3497-3541 ◽  
Author(s):  
B. C. Gunter ◽  
O. Didova ◽  
R. E. M. Riva ◽  
S. R. M. Ligtenberg ◽  
J. T. M. Lenaerts ◽  
...  

Abstract. This study explores an approach that simultaneously estimates Antarctic mass balance and glacial isostatic adjustment (GIA) through the combination of satellite gravity and altimetry data sets. The results improve upon previous efforts by incorporating reprocessed data sets over a longer period of time, and now include a firn densification model to account for firn compaction and surface processes. A range of different GRACE gravity models were evaluated, as well as a new ICESat surface height trend map computed using an overlapping footprint approach. When the GIA models created from the combination approach were compared to in-situ GPS ground station displacements, the vertical rates estimated showed consistently better agreement than existing GIA models. In addition, the new empirically derived GIA rates suggest the presence of strong uplift in the Amundsen Sea and Philippi/Denman sectors, as well as subsidence in large parts of East Antarctica. The total GIA mass change estimates for the entire Antarctic ice sheet ranged from 53 to 100 Gt yr−1, depending on the GRACE solution used, and with an estimated uncertainty of ±40 Gt yr−1. Over the time frame February 2003–October 2009, the corresponding ice mass change showed an average value of −100 ± 44 Gt yr−1 (EA: 5 ± 38, WA: −105 ± 22), consistent with other recent estimates in the literature, with the mass loss mostly concentrated in West Antarctica. The refined approach presented in this study shows the contribution that such data combinations can make towards improving estimates of present day GIA and ice mass change, particularly with respect to determining more reliable uncertainties.


2016 ◽  
Author(s):  
Godfred Osukuku ◽  
Abiud Masinde ◽  
Bernard Adero ◽  
Edmond Wanjala ◽  
John Ego

Abstract This research work attempts to map out the stratigraphic sequence of the Kerio Valley Basin using magnetic, gravity and seismic data sets. Regional gravity data consisting of isotactic, free-air and Bouguer anomaly grids were obtained from the International Gravity Bureau (BGI). Magnetic data sets were sourced from the Earth Magnetic Anomaly grid (EMAG2). The seismic reflection data was acquired in 1989 using a vibrating source shot into inline geophones. Gravity Isostacy data shows low gravity anomalies that depict a deeper basement. Magnetic tilt and seismic profiles show sediment thickness of 2.5-3.5 Km above the basement. The Kerio Valley Basin towards the western side is underlain by a deeper basement which are overlain by succession of sandstones/shales and volcanoes. At the very top are the mid Miocene phonolites (Uasin Gishu) underlain by mid Miocene sandstones/shales (Tambach Formation). There are high gravity anomalies in the western and southern parts of the basin with the sedimentation being constrained by two normal faults. The Kerio Valley Basin is bounded to the west by the North-South easterly dipping fault system. Gravity data was significantly of help in delineating the basement, scanning the lithosphere and the upper mantle according to the relative densities. The basement rocks as well as the upper cover of volcanoes have distinctively higher densities than the infilled sedimentary sections within the basin. From the seismic profiles, the frequency of the shaley rocks and compact sandstones increases with depths. The western side of the basin is characterized by the absence of reflections and relatively higher frequency content. The termination of reflectors and the westward dip of reflectors represent a fault (Elgeyo fault). The reflectors dip towards the west, marking the basin as an asymmetrical syncline, indicating that the extension was towards the east. The basin floor is characterized by a nearly vertical fault which runs parallel to the Elgeyo fault. The seismic reflectors show marked discontinuities which may be due to lava flows. The deepest reflector shows deep sedimentation in the basin and is in reasonable agreement with basement depths delineated from potential methods (gravity and magnetic). Basement rocks are deeper at the top of the uplift footwall of the Elgeyo Escarpment. The sediments are likely of a thickness of about 800 M which is an interbed of sandstones and shales above the basement.


2021 ◽  
Author(s):  
Yan Ming Wang ◽  
Xiaopeng Li ◽  
Kevin Ahlgren ◽  
Jordan Krcmaric ◽  
Ryan Hardy ◽  
...  

&lt;p&gt;For the upcoming North American-Pacific Geopotential Datum of 2022, the National Geodetic Survey (NGS), the Canadian Geodetic Survey (CGS) and the&amp;#160;National Institute of Statistics and Geography of Mexico (INEGI) computed the first joint experimental gravimetric geoid model (xGEOID) on 1&amp;#8217;x1&amp;#8217; grids that covers a region bordered by latitude 0 to 85 degree, longitude 180 to 350 degree east.&amp;#160;xGEOID20 models are computed using terrestrial gravity data, the latest satellite gravity model GOCO06S, altimetric gravity data DTU15, and an additional nine airborne gravity blocks of the GRAV-D project, for a total of 63 blocks. In addition, a digital elevation model in a 3&amp;#8221; grid was produced by combining MERIT, TanDEM-X, and USGS-NED and used for the topographic/gravimetric reductions. The geoid models computed from the height anomalies (NGS) and from the Helmert-Stokes scheme (CGS) were combined using two different weighting schemes, then evaluated against the independent GPS/leveling data sets. The models perform in a very similar way, and the geoid comparisons with the most accurate Geoid Slope Validation Surveys (GSVS) from 2011, 2014 and 2017 indicate that the&amp;#160;relative geoid accuracy&amp;#160;could be around 1-2 cm baseline lengths up to 300 km for these GSVS lines in the United States. The xGEOID20 A/B models were selected from the combined models based on the validation results. The geoid accuracies were also estimated using the forward modeling.&lt;/p&gt;


2018 ◽  
Vol 8 (1) ◽  
pp. 145-153 ◽  
Author(s):  
O.I. Apeh ◽  
E.C. Moka ◽  
V.N. Uzodinma

Abstract Spherical harmonic expansion is a commonly applied mathematical representation of the earth’s gravity field. This representation is implied by the potential coeffcients determined by using elements/parameters of the field observed on the surface of the earth and/or in space outside the earth in the spherical harmonic expansion of the field. International Centre for Gravity Earth Models (ICGEM) publishes, from time to time, Global Gravity Field Models (GGMs) that have been developed. These GGMs need evaluation with terrestrial data of different locations to ascertain their accuracy for application in those locations. In this study, Bouguer gravity anomalies derived from a total of eleven (11) recent GGMs, using sixty sample points, were evaluated by means of Root-Mean-Square difference and correlation coeficient. The Root-Mean-Square differences of the computed Bouguer anomalies from ICGEMwebsite compared to their positionally corresponding terrestrial Bouguer anomalies range from 9.530mgal to 37.113mgal. Additionally, the correlation coe_cients of the structure of the signal of the terrestrial and GGM-derived Bouguer anomalies range from 0.480 to 0.879. It was observed that GECO derived Bouguer gravity anomalies have the best signal structure relationship with the terrestrial data than the other ten GGMs. We also discovered that EIGEN-6C4 and GECO derived Bouguer anomalies have enormous potential to be used as supplements to the terrestrial Bouguer anomalies for Enugu State, Nigeria.


Sign in / Sign up

Export Citation Format

Share Document