scholarly journals Inner zone terrain correction calculation using interpolated heights

2015 ◽  
Vol 45 (3) ◽  
pp. 219-235 ◽  
Author(s):  
Pavol Zahorec

Abstract The discrepancy between real heights of gravity points and the elevation model has a significant impact on the terrain corrections calculation especially within the inner zone. The concept of interpolated heights of calculation points used instead of measured ones within the specified inner zone can considerably decrease the resulting errors. The choice of appropriate radius of the inner zone for use of interpolated heights is analysed on synthetic topography model as well as real data. The tests with synthetic models showed the appropriate radius of this zone is proportional to the deformation wavelength of the model. Simple statistical analysis of a particular elevation model can give an estimate of the appropriate radius for the calculation using interpolated heights. A concept with interpolated heights in the zone 0–250 m is used in actual practice in Slovakia. The analysis of regional gravity data from the Tatry Mountains test area indicates the searched radius should be about 100 m. Detailed gravity measurements from different areas showed the searched radius does not play so important role but the use of interpolated heights instead of measured ones is still relevant. The more reasonable method instead of using interpolated heights is also presented when calculating the topographic effect.

1965 ◽  
Vol 5 (41) ◽  
pp. 637-650 ◽  
Author(s):  
Charles E. Corbató

AbstractGravity measurements at 146 stations on lower Blue Glacier were used to determine the subglacial bedrock configuration. The gravity values, station elevations and density contrast were carefully measured, and terrain corrections thoroughly evaluated to insure accuracy of the Bottguer anomalies. A series of successive approximations results in evaluation of the regional gravity field and a three-dimensional model of the glacier whose gravimetric effects fit the range of the observational and computational errors. Comparison with bore holes and seismic reflections indicates no significant errors in the model and accuracies of 5–10 per cent in the calculated thicknesses of the glacier.


2019 ◽  
Vol 49 (2) ◽  
pp. 207-227 ◽  
Author(s):  
Pavol Zahorec ◽  
Juraj Papčo ◽  
Peter Vajda ◽  
Stanislav Szabó

Abstract Results from a detailed gravity survey realized along the planned highway tunnel in the karstic area of Slovak Karst in the eastern Slovakia are presented. Detailed gravity profiles crossed an area of rugged topography, therefore the terrain corrections played a crucial role in the gravity data processing. The airborne laser scanning technique (LiDAR) was used in order to compile a high-resolution digital terrain model (DTM) of the surrounding area and to calculate terrain corrections properly. The difference between the Bouguer anomalies calculated with an available nationwide DTM and those with new LiDAR-based model can be significant in some places as it is presented in the paper. A new method for Bouguer correction density analysis based on surface data is presented. Special underground gravity measurements in the existing nearby railway tunnel were also conducted in order to determine the mean density of the topographic rocks. The Bouguer anomalies were used to interpret lithological contacts and tectonic/karstic discontinuities.


1965 ◽  
Vol 5 (41) ◽  
pp. 637-650 ◽  
Author(s):  
Charles E. Corbató

AbstractGravity measurements at 146 stations on lower Blue Glacier were used to determine the subglacial bedrock configuration. The gravity values, station elevations and density contrast were carefully measured, and terrain corrections thoroughly evaluated to insure accuracy of the Bottguer anomalies. A series of successive approximations results in evaluation of the regional gravity field and a three-dimensional model of the glacier whose gravimetric effects fit the range of the observational and computational errors. Comparison with bore holes and seismic reflections indicates no significant errors in the model and accuracies of 5–10 per cent in the calculated thicknesses of the glacier.


2021 ◽  
Vol 95 (2) ◽  
Author(s):  
Mirjam Bilker-Koivula ◽  
Jaakko Mäkinen ◽  
Hannu Ruotsalainen ◽  
Jyri Näränen ◽  
Timo Saari

AbstractPostglacial rebound in Fennoscandia causes striking trends in gravity measurements of the area. We present time series of absolute gravity data collected between 1976 and 2019 on 12 stations in Finland with different types of instruments. First, we determine the trends at each station and analyse the effect of the instrument types. We estimate, for example, an offset of 6.8 μgal for the JILAg-5 instrument with respect to the FG5-type instruments. Applying the offsets in the trend analysis strengthens the trends being in good agreement with the NKG2016LU_gdot model of gravity change. Trends of seven stations were found robust and were used to analyse the stabilization of the trends in time and to determine the relationship between gravity change rates and land uplift rates as measured with global navigation satellite systems (GNSS) as well as from the NKG2016LU_abs land uplift model. Trends calculated from combined and offset-corrected measurements of JILAg-5- and FG5-type instruments stabilized in 15 to 20 years and at some stations even faster. The trends of FG5-type instrument data alone stabilized generally within 10 years. The ratio between gravity change rates and vertical rates from different data sets yields values between − 0.206 ± 0.017 and − 0.227 ± 0.024 µGal/mm and axis intercept values between 0.248 ± 0.089 and 0.335 ± 0.136 µGal/yr. These values are larger than previous estimates for Fennoscandia.


2012 ◽  
Vol 2 (1) ◽  
pp. 53-64 ◽  
Author(s):  
H. Yildiz ◽  
R. Forsberg ◽  
J. Ågren ◽  
C. Tscherning ◽  
L. Sjöberg

Comparison of remove-compute-restore and least squares modification of Stokes' formula techniques to quasi-geoid determination over the Auvergne test areaThe remove-compute-restore (RCR) technique for regional geoid determination implies that both topography and low-degree global geopotential model signals are removed before computation and restored after Stokes' integration or Least Squares Collocation (LSC) solution. The Least Squares Modification of Stokes' Formula (LSMS) technique not requiring gravity reductions is implemented here with a Residual Terrain Modelling based interpolation of gravity data. The 2-D Spherical Fast Fourier Transform (FFT) and the LSC methods applying the RCR technique and the LSMS method are tested over the Auvergne test area. All methods showed a reasonable agreement with GPS-levelling data, in the order of a 3-3.5 cm in the central region having relatively smooth topography, which is consistent with the accuracies of GPS and levelling. When a 1-parameter fit is used, the FFT method using kernel modification performs best with 3.0 cm r.m.s difference with GPS-levelling while the LSMS method gives the best agreement with GPS-levelling with 2.4 cm r.m.s after a 4-parameter fit is used. However, the quasi-geoid models derived using two techniques differed from each other up to 33 cm in the high mountains near the Alps. Comparison of quasi-geoid models with EGM2008 showed that the LSMS method agreed best in term of r.m.s.


2017 ◽  
Vol 5 (2) ◽  
pp. T243-T257 ◽  
Author(s):  
Martin Landrø ◽  
Mark Zumberge

We have developed a calibrated, simple time-lapse seismic method for estimating saturation changes from the [Formula: see text]-storage project at Sleipner offshore Norway. This seismic method works well to map changes when [Formula: see text] is migrating laterally away from the injection point. However, it is challenging to detect changes occurring below [Formula: see text] layers that have already been charged by some [Formula: see text]. Not only is this partly caused by the seismic shadow effects, but also by the fact that the velocity sensitivity for [Formula: see text] change in saturation from 0.3 to 1.0 is significantly less than saturation changes from zero to 0.3. To circumvent the seismic shadow zone problem, we combine the time-lapse seismic method with time-lapse gravity measurements. This is done by a simple forward modeling of gravity changes based on the seismically derived saturation changes, letting these saturation changes be scaled by an arbitrary constant and then by minimizing the least-squares error to obtain the best fit between the scaled saturation changes and the measured time-lapse gravity data. In this way, we are able to exploit the complementary properties of time-lapse seismic and gravity data.


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.


Geophysics ◽  
2013 ◽  
Vol 78 (2) ◽  
pp. G15-G24 ◽  
Author(s):  
Pejman Shamsipour ◽  
Denis Marcotte ◽  
Michel Chouteau ◽  
Martine Rivest ◽  
Abderrezak Bouchedda

The flexibility of geostatistical inversions in geophysics is limited by the use of stationary covariances, which, implicitly and mostly for mathematical convenience, assumes statistical homogeneity of the studied field. For fields showing sharp contrasts due, for example, to faults or folds, an approach based on the use of nonstationary covariances for cokriging inversion was developed. The approach was tested on two synthetic cases and one real data set. Inversion results based on the nonstationary covariance were compared to the results from the stationary covariance for two synthetic models. The nonstationary covariance better recovered the known synthetic models. With the real data set, the nonstationary assumption resulted in a better match with the known surface geology.


2021 ◽  
Vol 6 (24) ◽  
pp. 213-225
Author(s):  
Shazad Jamal Jalal ◽  
Tajul Ariffin Musa ◽  
Ami Hassan Md Din ◽  
Wan Anom Wan Aris

Gravity data and computing gravity anomalies are regarded as vital for both geophysics and physical geodesy fields. The mountainous areas of Iraq are characterized by the lack of regional gravity data because gravity surveys are rarely performed in the past four decades due to the Iraq-Iran war and the internal unstable political situation of this particular region. In addition, the formal map of the available terrestrial gravity which was published by the French Database of Bureau Gravimetrique International (International Gravimetric Bureau-in English) (BGI), introduces Iraq and the study area as a remote area and in white color because of the unavailability of gravity data. However, a dense and local (not regional) gravity data is available which was conducted by geophysics researchers 13 years ago. Therefore, the regional gravity survey of 160 gravity points was performed by the authors at an average 11 km apart, which was covers the whole area of Sulaymaniyah Governorate (part of the mountainous areas of Iraq). In spite of Although the risk of mine fields within the study area, suitable safe routes as well as a helicopter was used for the gravity survey of several points on the top of mountains. The survey was conducted via Lacoste and Romberg geodetic gravimeter and GPS handheld. The objective of the study is to determine and map the gravity anomalies for the entire study area, the data of which would assist different geosciences applications.


Geosciences ◽  
2021 ◽  
Vol 11 (11) ◽  
pp. 482
Author(s):  
Dharmendra Kumar ◽  
Arun Singh ◽  
Mohammad Israil

The magnetotelluric (MT) method is one of the useful geophysical techniques to investigate deep crustal structures. However, in hilly terrains, e.g., the Garhwal Himalayan region, due to the highly undulating topography, MT responses are distorted. Such responses, if not corrected, may lead to the incorrect interpretation of geoelectric structures. In the present paper, we implemented terrain corrections in MT data recorded from the Garhwal Himalayan Corridor (GHC). We used AP3DMT, a 3D MT data modeling and inversion code written in the MATLAB environment. Terrain corrections in the MT impedance responses for 39 sites along the Roorkee–Gangotri profile in the period range of 0.01 s to 1000 s were first estimated using a synthetic model by recording the topography and locations of MT sites. Based on this study, we established the general character of the terrain and established where terrain corrections were necessary. The distortion introduced by topography was computed for each site using homogenous and heterogeneous models with actual topographic variations. Period-dependent, galvanic and inductive distortions were observed at different sites. We further applied terrain corrections to the real data recorded from the GHC. The corrected data were inverted, and the inverted model was compared with the corresponding inverted model obtained with uncorrected data. The modification in electrical resistivity features in the model obtained from the terrain-corrected response suggests the necessity of terrain correction in MT data recorded from the Himalayan region.


Sign in / Sign up

Export Citation Format

Share Document