Upward continuation of Dome-C airborne gravity and comparison with GOCE gradients at orbit altitude in east Antarctica

An airborne gravity campaign was carried out at the Dome-C survey area in East Antarctica between the 17th and 22nd of January 2013, in order to provide data for an experiment to validate GOCE satellite gravity gradients. After typical filtering for airborne gravity data, the cross-over error statistics for the few crossing points are 11.3 mGal root mean square (rms) error, corresponding to an rms line error of 8.0 mGal. This number is relatively large due to the rough flight conditions, short lines and field handling procedures used. Comparison of the airborne gravity data with GOCE RL4 spherical harmonic models confirmed the quality of the airborne data and that they contain more high-frequency signal than the global models. First, the airborne gravity data were upward continued to GOCE altitude to predict gravity gradients in the local North-East-Up reference frame. In this step, the least squares collocation using the ITGGRACE2010S field to degree and order 90 as reference field, which is subtracted from both the airborne gravity and GOCE gravity gradients, was applied. Then, the predicted gradients were rotated to the gradiometer reference frame using level 1 attitude quaternion data. The validation with the airborne gravity data was limited to the accurate gradient anomalies (T_XX, T_YY, T_ZZ and T_XZ) where the long-wavelength information of the GOCE gradients has been replaced with GOCO03s signal to avoid contamination with GOCE gradient errors at these wavelengths. The comparison shows standard deviations between the predicted and GOCE gradient anomalies T_XX, T_YY, T_ZZ and T_XZ of 9.9, 11.5, 11.6 and 10.4 mE, respectively. A more precise airborne gravity survey of the southern polar gap which is not observed by GOCE would thus provide gradient predictions at a better accuracy, complementing the GOCE coverage in this region.     

扰动重力梯度张量单分量和组合分量最小二乘配置法模型的建立

建立了利用扰动重力梯度张量T_zz分量和T_xx+T_yy、T_zz-T_xx-T_yy组合分量确定地球重力场的调和分析法模型,进一步推导了扰动重力梯度张量对角线三分量的自协方差和互协方差函数的级数展开式,推导了单分量、组合分量与重力位系数之间协方差函数的实用计算公式,给出了利用单分量和组合分量解算地球重力场模型的最小二乘配置法基本原理公式.结果表明,最小二乘配置法具有一定的抗差能力,随着观测数据误差的不断增大,其恢复的重力场模型有效阶次不断降低,精度也不断下降;T_zz-T_xx-T_yy组合分量解算重力场模型的精度最高,其次为T_zz分量,T_xx+T_yy组合分量最差.     

Using ground gravity to improve ice mass change estimation from GOCE gravity gradients in mid-west Greenland

Vertical gravity gradient anomalies from the Gravity and steady-state Ocean Circulation Explorer (GOCE) DIR-3 model have been used to determine gravity anomalies in mid-west Greenland by using Least-Squares Collocation (LSC) and the Reduced Point Mass (RPM) method. The two methods give nearly identical results. However, compared to LSC, the RPM method needs less computational time as the number of equations to be solved in LSC equals the number of observations. The advantage of the LSC, however, is the acquired error estimates. The observation periods are winter 2009 and summer 2012. In order to enhance the accuracy of the calculated gravity anomalies, ground gravity data from West Greenland is used over locations where the gravity change resulting from ice mass changes is negligible, i.e. over solid rock. In the period considered, the gravity anomaly change due to changes in ice mass varies from ?5 mGal to 4 mGal. It is negative over the outlet glacier Jacobshavn Isbræ, where the mass loss corresponds to a gravity change of approximately ?4 mGal. When using only GOCE vertical gravity gradients, the error estimates range from 5 mGal at the coast to 17 mGal over the ice sheet. Introducing the ground gravity data from West Greenland in the prediction reduces the errors to range from 2 to 10 mGal.     

An estimation of the height system bias parameter <Emphasis Type="Italic">N</Emphasis><Subscript>0</Subscript> using least squares collocation from observed gravity and GPS-levelling data

This paper deals with the analysis of gravity anomaly and precise levelling in conjunction with GPS-Levelling data for the computation of a gravimetric geoid and an estimate of the height system bias parameter N_o for the vertical datum in Pakistan by means of least squares collocation technique. The long term objective is to obtain a regional geoid (or quasi-geoid) modeling using a combination of local data with a high degree and order Earth gravity model (EGM) and to determine a bias (if there is one) with respect to a global mean sea surface. An application of collocation with the optimal covariance parameters has facilitated to achieve gravimetric height anomalies in a global geocentric datum. Residual terrain modeling (RTM) technique has been used in combination with the EGM96 for the reduction and smoothing of the gravity data. A value for the bias parameter N_o has been estimated with reference to the local GPS-Levelling datum that appears to be 0.705 m with 0.07 m mean square error. The gravimetric height anomalies were compared with height anomalies obtained from GPS-Levelling stations using least square collocation with and without bias adjustment. The bias adjustment minimizes the difference between the gravimetric height anomalies with respect to residual GPS-Levelling data and the standard deviation of the differences drops from 35 cm to 2.6 cm. The results of this study suggest that N_o adjustment may be a good alternative for the fitting of the final gravimetric geoid as is generally done when using FFT methods.     

Regional gravimetric quasi-geoid model and transformation surface to national height system for Turkey (THG-09)

Turkish regional geoid models have been developed by employing a reference earth gravitational model, surface gravity observations and digital terrain models. The gravimetric geoid models provide a ready transformation from ellipsoidal heights to the orthometric heights through the use of GPS/leveling geoid heights determined through the national geodetic networks. The recent gravimetric models for Turkish territory were computed depending on OSU91 (TG-91) and EGM96 (TG-03) earth gravitational models. The release of the Earth Gravitational Model 2008 (EGM08), the collection of new surface gravity observations, the advanced satellite altimetry-derived gravity over the sea, and the availability of the high resolution digital terrain model have encouraged us to compute a new geoid model for Turkey. We used the Remove-Restore procedure based on EGM08 and applied Residual Terrain Model (RTM) reduction of the surface gravity data. Fast Fourier Transformation (FFT) was then used to obtain the residual quasigeoid from the reduced gravity. We restored the individual contributions of EGM08 and RTM to the whole quasi-geoid height (TQG-09). Since the Helmert orthometric height system is adopted in Turkey, the quasi-geoid model (TQG-09) was then converted to the geoid model (TG-09) by making use of Bouguer gravity anomalies and digital terrain model. After all we combined a gravimetric geoid model with GPS/leveling geoid heights in order to obtain a hybrid geoid model (THG-09) (or a transformation surface) to be used in GPS applications. The RMS of the post-fit residuals after the combination was found to be ± 0.95 cm, which represents the internal precision of the final combination. And finally, we tested the hybrid geoid model with GPS/leveling data, which were not used in the combination, to assess the external accuracy. Results show that the external accuracy of the THG-09 model is ± 8.4 cm, a precision previously not achieved in Turkey until this study.     

Error sources and data limitations for the prediction of surface gravity: a case study using benchmarks

Gravity-based heights require gravity values at levelled benchmarks (BMs), which sometimes have to be predicted from surrounding observations. We use the Earth Gravitational Model 2008 (EGM2008) and the Australian National Gravity Database (ANGD) as examples of model and terrestrial observed data respectively to predict gravity at Australian National Levelling Network (ANLN) BMs. The aim is to quantify errors that may propagate into the predicted BM gravity values and then into gravimetric height corrections (HCs). Our results indicate that an approximate ±1 arc-min horizontal position error of the BMs causes maximum errors in EGM2008 BM gravity of ~22 mGal (~55 mm in the HC at ~2200 m elevation) and ~18 mGal for ANGD BM gravity because the values are not computed at the true location of the BM. We use RTM (residual terrain modelling) techniques to show that ~50% of EGM2008 BM gravity error in a moderately mountainous region can be accounted for by signal omission. Non-representative sampling of ANGD gravity in this region may cause errors of up to 50 mGals (~120 mm for the Helmert orthometric correction at ~2200 m elevation). For modelled gravity at BMs to be viable, levelling networks need horizontal BM positions accurate to a few metres, while RTM techniques can be used to reduce signal omission error. Unrepresentative gravity sampling in mountains can be remedied by denser and more representative re-surveys, and/or gravity can be forward modelled into regions of sparser gravity.     

Investigation of topographic reductions and aliasing effects on gravity and the geoid over Greece based on various digital terrain models

The reduction of gravity-field related quantities (e.g., gravity anomalies, geoid heights) due to the topography plays a crucial role in both geodetic and geophysical applications, since in the former it is an intermediate step towards geoid prediction and in the latter it reveals lateral as well as radial density contrasts and infers the geology of the area under study. The computations are usually carried out by employing a DTM and/or a DBM, which describe the topography and bathymetry, respectively. Errors in these DTMs/DBMs will introduce errors in the computed topographic effects, while poor spatial resolution of the topography and bathymetry models will result in aliasing effects to both gravity anomalies and geoid heights, both influencing the accuracy of the estimated solutions. The scope of this work is twofold. First, a validation and accuracy assessment of the SRTM 3″ (90 m) DTM over Greece is performed through comparisons with existing global models as well as with the Greek 450 m national DTMs. Whenever a misrepresentation of the topography is identified in the SRTM data, it is “corrected” using the local 450 m DTM. This process resulted in an improved SRTM DTM called SRTMGr, which was then used to determine terrain effects to gravity field quantities. From the fine-resolution SRTMGr DTMs, coarser models of 15″, 30″, 1′, 2′ and 5′ have been determined in order to investigate aliasing effects on both gravity anomalies and geoid heights by computing terrain effects at variable spatial resolutions. From the results acquired in two test areas, it was concluded that SRTMGr provides similar results to the local DTM making the use of other older global DTMs obsolete. The study for terrain aliasing effects proved that when high-resolution and accuracy gravity and geoid models are needed, then the highest possible resolution DTM should be employed to compute the respective terrain effects. Based on the results acquired from two the test areas a corrected SRTMGr DTM has been compiled for the entire Greek territory towards the development of a new gravimetric geoid model. Results from that analysis are presented based on the well-known remove-compute-restore method, employing land and marine gravity data, EGM08 as a reference geopotential model and the SRTMGr DTM for the computation of the RTM effects.     

Determination of gravity anomaly at sea level from inversion of satellite gravity gradiometric data

Gravity gradients can be used to determine the local gravity field of the Earth. This paper investigates downward continuation of all elements of the disturbing gravitational tensor at satellite level using the second-order partial derivatives of the extended Stokes formula in the local-north oriented frame to determine the gravity anomaly at sea level. It considers the inversion of each gradient separately as well as their joint inversion. Numerical studies show that the gradients T_zz, T_xx, T_yy and T_xz have similar capability of being continued downward to sea level in the presence of white noise, while the gradient T_yz is considerably worse than the others. The bias-corrected joint inversion process shows the possibility of recovering the gravity anomaly with 1 mGal accuracy. Variance component estimation is also tested to update the observation weights in the joint inversion.     

Experiences with the use of mass-density maps in residual gravity forward modelling

In many modern local and regional gravity field modelling concepts, the short-wavelength gravitational signal modeled by the residual terrain modelling (RTM) technique is used to augment global geopotential models, or to smooth observed gravity prior to data gridding. In practice, the evaluation of RTM effects mostly relies on a constant density assumption, because of the difficulty and complexity of obtaining information on the actual distribution of density of topographic masses. Where the actual density of topographic masses deviates from the adopted value, errors are present in the RTM mass-model, and hence, in the forward-modelled residual gravity field. In this paper we attempt to overcome this problem by combining the RTM technique with a high-resolution mass-density model. We compute RTM gravity quantities over New Zealand, with different combinations of elevation models and mass-density assumptions using gravity and GPS/levelling measurements, precise terrain and bathymetry models, a high-resolution mass-density model and constant density assumptions as main input databases. Based on gravity observations and the RTM technique, optimum densities are detected for North Island of ~2500 kg m^?3, South Island of ~2600 kg m^?3, and the whole New Zealand of ~2590 kg m^?3. Comparison among the three sets of residual gravity disturbances computed from different mass-density assumptions show that, together with a global potential model, the high-resolution New Zealand density model explains ~89.5% of gravitational signals, a constant density assumption of 2670 kg m^?3 explains ~90.2%, while a regionally optimum mass-density explains ~90.3%. Detailed comparison shows that the New Zealand density model works best over areas with small residual heights. Over areas with larger residual heights, subsurface density variations appear to affect the residual gravity disturbance. This effect is found to reach about 30 mGal over Southern Alpine Fault. In order to improve the RTM modelling with mass-density maps, a higher-quality mass-density model that provides radially varying mass-density data would be desirable.     

RTM-based omission error corrections for global geopotential models: Case study in Central Europe

The aim of this paper is to evaluate the effects of residual terrain model (RTM) on potential and on gravity and to point out how significant can the omission error of global geopotential models (GGMs) be and how it can influence their testing. The RTM for Central Europe is computed in the spherical approximation. The topography is modelled by spherical tesseroids and the gravitational effect of the topography is obtained as a sum of their partial gravitational effects. A detailed picture of RTM in Slovakia is shown. The testing of GOCE (Gravity Field and Steady-State Ocean Circulation Explorer) global geopotential models in Central Europe published earlier is re-evaluated with the more rigorous omission error estimation. Experimental results show significantly better agreement between the gravity anomalies computed from global geopotential models with the omission-error estimation and gravity anomalies obtained from the direct measurements. On the other hand, for height anomalies such an improvement is not observed. The results are discussed in context of the other previously published studies.     

On High-Frequency Topography-Implied Gravity Signals for a Height System Unification Using GOCE-Based Global Geopotential Models

National height reference systems have conventionally been linked to the local mean sea level, observed at individual tide gauges. Due to variations in the sea surface topography, the reference levels of these systems are inconsistent, causing height datum offsets of up to ±1–2 m. For the unification of height systems, a satellite-based method is presented that utilizes global geopotential models (GGMs) derived from ESA’s satellite mission Gravity field and steady-state Ocean Circulation Explorer (GOCE). In this context, height datum offsets are estimated within a least squares adjustment by comparing the GGM information with measured GNSS/leveling data. While the GNSS/leveling data comprises the full spectral information, GOCE GGMs are restricted to long wavelengths according to the maximum degree of their spherical harmonic representation. To provide accurate height datum offsets, it is indispensable to account for the remaining signal above this maximum degree, known as the omission error of the GGM. Therefore, a combination of the GOCE information with the high-resolution Earth Gravitational Model 2008 (EGM2008) is performed. The main contribution of this paper is to analyze the benefit, when high-frequency topography-implied gravity signals are additionally used to reduce the remaining omission error of EGM2008. In terms of a spectral extension, a new method is proposed that does not rely on an assumed spectral consistency of topographic heights and implied gravity as is the case for the residual terrain modeling (RTM) technique. In the first step of this new approach, gravity forward modeling based on tesseroid mass bodies is performed according to the Rock–Water–Ice (RWI) approach. In a second step, the resulting full spectral RWI-based topographic potential values are reduced by the effect of the topographic gravity field model RWI_TOPO_2015, thus, removing the long to medium wavelengths. By using the latest GOCE GGMs, the impact of topography-implied gravity signals on the estimation of height datum offsets is analyzed in detail for representative GNSS/leveling data sets in Germany, Austria, and Brazil. Besides considerable changes in the estimated offset of up to 3 cm, the conducted analyses show that significant improvements of 30–40% can be achieved in terms of a reduced standard deviation and range of the least squares adjusted residuals.     

Transformation between gravimetric and GPS/levelling-derived geoids using additional gravity information

The transformation from the gravimetric to the GPS/levelling-derived geoid using additional gravity information for the covariance function of geoid height differences has been investigated in a test area in south-western Canada. A “corrector surface” model, which accounts for datum inconsistencies, long-wavelength geoid errors, vertical network distortions and GPS errors, has been constructed using least-squares collocation. The local covariance function of geoid height differences is usually obtained from residual values between the GPS/levelling and gravimetric geoid heights after the elimination of all known systematic distortions. If additional gravity data (in the form of gravity anomalies) are available, the covariance function of geoid height differences can be determined by the following steps: (1) transforming the GPS/levelling-derived geoid heights into gravity anomalies; (2) forming differences between the computed in step 1 and given gravity anomalies; (3) determining the parameters of the local covariance function of the gravity anomaly differences; (4) constructing an analytical covariance model for the geoid height differences from the covariance function of the gravity anomaly differences using the parameters derived in step 3. The advantage of the proposed method stems from the great number of gravity data used to derive the empirical covariance function. A comparison with the least-squares adjustment shows that the standard deviation of the residuals of the predicted geoid height differences with respect to the control point values decreases by 2.4 cm.     

我国重力场新的统计特征参数的计算分析

基于地形信息,构造了完全布格异常和完全空间异常的协方差和代表误差模型参数等新统计量,利用某试验区密集的重力和地形数据进行了统计分析,并将计算方案推广到了全国大范围、多区域内的统计计算,给出了6种不同地形类别区域（平原、丘陵、小山区、中山区、大山区、特大山区）的完全空间异常和完全布格异常的方差、协方差以及代表误差模型参数.试验结果表明：依据本文的统计模型、算法与思路,在实际测量数据的支撑下,可以给出全国范围的统计参数的网格值、等值线图等,提供后续重力场相关研究工作使用.     

A comparison of GOCE gravitational models with EGM2008

Four new gravity field models from GOCE, two of them combined with GRACE, are compared here with EGM2008. The objectives are to look into the differences in consecutive ranges of the spherical harmonic expansion globally as well as in selected geographical regions and in the regions of the various data sources used for EGM2008. In general, GOCE is able to contribute to improved global gravity models in the spherical harmonic range between 120 and 200 (and above). The agreement between EGM2008 and the GOCE models is very good in well-surveyed regions such as North America, Europe and Australia, with geoid RMS-differences on the order of 4–6 cm. In other regions, where the surface gravity data available for the development of EGM2008 were poor, such as South America, Africa, South-East Asia or China the RMS-differences are on a level of 30 cm. Here GOCE leads to a significant improvement. These findings are confirmed by the analysis of the areas of the various EGM2008 data sources. In the regions of the so-called “fill-in” data of EGM2008 RMS-geoid height differences are high. In Antarctica GOCE also gives important improvements in terms of spatial resolution and accuracy. In general, the agreement between EGM2008 and the GOCE-models up to degree and order (d/o) 200 is good, with a global (excluding the polar gaps of GOCE orbits, throughout) geoid difference RMS of 11 cm, in the ocean areas 8 cm and 20 cm in the continental areas. GOCE models are better suited for ocean circulation studies because no prior ocean information enters into the data reduction process, as it is the case when deducing gravity anomalies from an altimetric mean sea surface. On the other hand, the good consistency between GOCE-models and EGM2008 in ocean areas very likely indicates that the influence of ocean circulation information on EGM2008 is rather small. The four tested GOCE models behave similarly except at the highest latitudes where GOCE lacks data due to its orbit inclination of 96.5° and some form of regularization which has to be applied.     

GOCE, Satellite Gravimetry and Antarctic Mass Transports

In 2009 the European Space Agency satellite mission GOCE (Gravity Field and Steady-State Ocean Circulation Explorer) was launched. Its objectives are the precise and detailed determination of the Earth’s gravity field and geoid. Its core instrument, a three axis gravitational gradiometer, measures the gravity gradient components V _xx , V _yy , V _zz and V _xz (second-order derivatives of the gravity potential V) with high precision and V _xy , V _yz with low precision, all in the instrument reference frame. The long wavelength gravity field is recovered from the orbit, measured by GPS (Global Positioning System). Characteristic elements of the mission are precise star tracking, a Sun-synchronous and very low (260 km) orbit, angular control by magnetic torquing and an extremely stiff and thermally stable instrument environment. GOCE is complementary to GRACE (Gravity Recovery and Climate Experiment), another satellite gravity mission, launched in 2002. While GRACE is designed to measure temporal gravity variations, albeit with limited spatial resolution, GOCE is aiming at maximum spatial resolution, at the expense of accuracy at large spatial scales. Thus, GOCE will not provide temporal variations but is tailored to the recovery of the fine scales of the stationary field. GRACE is very successful in delivering time series of large-scale mass changes of the Antarctic ice sheet, among other things. Currently, emphasis of respective GRACE analyses is on regional refinement and on changes of temporal trends. One of the challenges is the separation of ice mass changes from glacial isostatic adjustment. Already from a few months of GOCE data, detailed gravity gradients can be recovered. They are presented here for the area of Antarctica. As one application, GOCE gravity gradients are an important addition to the sparse gravity data of Antarctica. They will help studies of the crustal and lithospheric field. A second area of application is ocean circulation. The geoid surface from the gravity field model GOCO01S allows us now to generate rather detailed maps of the mean dynamic ocean topography and of geostrophic flow velocities in the region of the Antarctic Circumpolar Current.     

Gravity field contribution analysis of GOCE gravitational gradient components

A gravity field model is computed from the four accurate gravitational gradient components of GOCE (Gravity field and steady-state Ocean Circulation Explorer), combined with the analysis of the kinematic orbits, and some moderate constraint (or stabilization) in the polar areas where no observation from GOCE is available due to the orbit geometry. The normal matrix of each component is computed individually in order to study its contribution to the combined solution. The results show that the contribution of V_zz is the largest, with an average value of 32.74% of the total solution; the second and the third largest are V_zz and V_yy, with average values of 28.04% and 26.08%, respectively; the component V_xz contributes 11.81%. Validation with external data shows that each component has its characteristic value and that the information content of the component V_xz is not negligible and should be included for gravity field recovery. The orbit part as derived from high-low satellite-to-satellite tracking (SST-hl) to the GPS contributes mostly to the coefficients below degree and order (d/o) 20, and to non-zonal coefficients from d/o 20 to 80. The mean value of the contribution of the polar stabilization is the smallest with a value of 0.22%, nevertheless it is important. In addition to the contribution analysis in terms of the normal matrices, each individual component of the gradiometer has been combined with SST and polar stabilization, to give a set of single component gravity field models. These partially combined solutions are compared to the fully combined solution in terms of geoid differences. They show that the partially combined solution with V_zz is closest to the complete solution. Even closer is a combination with V_xx and V_yy. In addition to the GOCE-only solution, a GOCE-GRACE (Gravity Recovery And Climate Experiment) combined gravity field model is derived and the information content of GOCE and an available set of normal equations of GRACE are investigated. Results show that, as expected, GRACE dominates the solution below degree 90 and GOCE above degree 140.     

单加速度计模式下的GOCE卫星重力场建模方法研究

GOCE卫星由于加速度计的特殊安装方式,其非保守力主要由普通模式的组合加速度提供,使得单个加速度计的特征更难提取.本文首次采用实测数据,研究了单加速度计模式下的高低跟踪数据处理.利用GOCE任务2009年(2009-11—2009-12)的实测数据,分别以GOCE卫星梯度仪坐标系三个坐标轴正向的加速度计为研究对象,利用1s间隔的高采样轨道数据,采用动力法同时进行卫星重力场建模和加速度计的精密校准.为了克服两极地区的数据缺失对重力场模型低次系数的影响,即所谓的极空白问题,引入同期GRACE卫星的观测数据,采用方差分量估计方法,建立了GRACE/GOCE卫星跟踪卫星重力场模型WHU-GRGO-SST.该模型完全到100阶次,经6169个美国GPS水准点数据检验,在同阶次上与EGM2008和GGM05S的精度水平相同.分析发现,GOCE卫星的加速度计偏差参数存在显著的漂移,也显示了单加速度计模式处理GOCE高低跟踪数据的优势.本文的研究成果为建立静态高分辨率、高精度的GRACE/GOCE重力场模型提供了更严密的模型与技术方案,同时也为GOCE卫星梯度仪校准,以及梯度数据的深入分析提供了重要的参考信息.     

Interpretation of tensor gravity data using an adaptive tilt angle method

Full Tensor Gravity Gradiometry (FTG) data are routinely used in exploration programmes to evaluate and explore geological complexities hosting hydrocarbon and mineral resources. FTG data are typically used to map a host structure and locate target responses of interest using a myriad of imaging techniques. Identified anomalies of interest are then examined using 2D and 3D forward and inverse modelling methods for depth estimation. However, such methods tend to be time consuming and reliant on an independent constraint for clarification. This paper presents a semi‐automatic method to interpret FTG data using an adaptive tilt angle approach. The present method uses only the three vertical tensor components of the FTG data (T_zx, T_zy and T_zz) with a scale value that is related to the nature of the source (point anomaly or linear anomaly). With this adaptation, it is possible to estimate the location and depth of simple buried gravity sources such as point masses, line masses and vertical and horizontal thin sheets, provided that these sources exist in isolation and that the FTG data have been sufficiently filtered to minimize the influence of noise. Computation times are fast, producing plausible results of single solution depth estimates t hat relate directly to anomalies. For thick sheets, the method can resolve the thickness of these layers assuming the depth to the top is known from drilling or other independent geophysical data. We demonstrate the practical utility of the method using examples of FTG data acquired over the Vinton Salt Dome, Louisiana, USA and basalt flows in the Faeroe‐Shetland Basin, UK. A major benefit of the method is the ability to quickly construct depth maps. Such results are used to produce best estimate initial depth to source maps that can act as initial models for any detailed quantitative modelling exercises using 2D/3D forward/inverse modelling techniques.