- Split View
-
Views
-
Cite
Cite
Zdeněk Martinec, Mass-density Green's functions for the gravitational gradient tensor at different heights, Geophysical Journal International, Volume 196, Issue 3, March 2014, Pages 1455–1465, https://doi.org/10.1093/gji/ggt495
- Share Icon Share
Abstract
Four different forms of the tensor Green's function for the gravitational gradient tensor, derived in this article, give a theoretical basis for geophysical interpretations of the GOCE-based gravitational gradients in terms of the Earth's mass-density structure. The first form is an invariant expression of the tensor Green's function that can be used to evaluate numerically the gravitational gradients in different coordinate systems (e.g. Cartesian). The second form expresses the gravitational gradients in spherical coordinates (ϑ, φ) with the origin at the north pole as a series of tensor spherical harmonics. This form is convenient to apply when the GOCE data are represented in terms of the gravitational potential as a scalar spherical harmonic series, such as the GOCO03S satellite gravity model. The third form expresses gravitational gradients in spherical coordinates (ψ, α) with the pole at the computation point. The fourth form then expresses the corresponding isotropic kernels in a closed form. The last two forms are used to analyse the sensitivity of the gravitational gradients with respect to lateral distribution of the Earth's mass-density anomalies. They additionally provide a tool for evaluating the omission error of geophysically modelled gravitational gradients and its amplification when the bandwidth-limited GOCE-based gravitational gradients are interpreted at different heights above the Earth's surface. We show that the omission error of the bandwidth-limited mass-density Green's functions for gradiometric data at the GOCE satellite's altitude does not exceed 1 per cent in amplitude when compared to the full-spectrum Green's functions. However, when evaluating the bandwidth-limited Green's functions at lower altitudes, their omission errors are significantly amplified. In this case, we show that the short-wavelength content of the forward-modelled gravitational gradients generated by an a priori density structure of the Earth must be filtered out such that the omission error of the GOCE-based gravitational gradients (i.e. the signal that has not been modelled from the GOCE data) is equal to the omission error of the forward-modelled gravitational gradients. Only after performing such filtering can the GOCE-based gravitational gradients at low altitudes be interpreted in terms of the Earth's density structure. Both the closed and spherical-harmonic forms of the Green's functions allow a direct interpretation in terms of the minimum lateral extent of the Earth that needs to be considered in a regional model constrained by gravitational gradients if the full information provided by the Green's functions is to be retained. We show that this extent (i) is smallest for the vertical–vertical Green's function and (ii) linearly increases with increasing computation-point height. The spectral forms of the gravitational gradients are further used to calculate the sensitivity of the GOCE-based gravitational gradients to the depth of density anomalies, expressed in terms of the harmonic degree of the internal mass-density anomaly. We show that the largest gravitational gradient response is obtained for shallow mass anomalies, and is further amplified as the harmonic degree increases.
1 INTRODUCTION
The gravity field and steady-state ocean circulation explorer (GOCE), launched in 2009, provides gravitational gradients at a perigee height of 255 km with a spatial resolution of 90 km and a precision of 1 mGal (e.g. Rummel et al.2011). The GOCE gravitational gradient data are delivered by its on-board gradiometer. The gradients themselves are not directly measured, but derived from the differences between the accelerations of three pairs of accelerometers (Frommknecht et al.2011). The gravitational gradients are thus determined in the so-called gradiometer reference frame, which co-rotates with the satellite. The gradiometer-reference-frame configuration is such that four gravitational gradients are determined with a high level of accuracy, whereas the other two are less accurate. The GOCE gravitational gradients in the gradiometer reference frame are combined with complementary data from the gravity field missions CHAMP and GRACE (and also with SLR, terrestrial gravity data and satellite altimetry) and the combined global gravity field (GOCO) models are calculated (e.g. Pail et al.2013). In this text, we refer to the GOCE-based gravitational gradients as the gradients calculated from the spherical harmonic expansion of a GOCO model, for example GOCO03S (Mayer-Gürr et al.2012), with full spectral content up to degree and order 220.
The latest GOCO models have improved upon the available information about the terrestrial gravity field in comparison to the EGM2008 model (Pavlis et al.2012) in regions where the quantity and quality of terrestrial gravity data included in the EGM2008 model is poor (Bouman et al.2011; Bouman and Fuchs 2012; Hirt et al.2012). This is particularly evident in less well surveyed parts of the world, for example central Africa, where the noise in the terrestrial gravity data incorporated into the EGM2008 model has been substantially reduced. A number of studies have exploited the capabilities of the new GOCE data in constraining the structure of the Earth's crust and oceans (e.g. Bingham et al.2011; Álvarez et al.2012; Hirt et al.2012; Köther et al.2012; Mariani et al.2013), the lithosphere (Bouman et al.2013) and the upper mantle structure (Fullea et al.2013).
The representation of the gravitational potential of a GOCO model in terms of solid spherical harmonics allows the evaluation of the GOCE-based gravitational gradients at different heights above the Earth's surface, as well as the possibility of interpreting the GOCE-based gravitational gradients in terms of the solid Earth's structure, again at different altitudes. Two different approaches based on the height where the GOCE-based gravitational gradients are interpreted by modelled quantities have been considered. First, for the identification of geological units in unexplored parts of the world, such as the central part of Africa, it is useful to evaluate the GOCE-based gravitational gradients at altitudes close to the Earth's surface, since the gradiometric signal is amplified and better reflects the near-surface geological structure. However, the signal induced by the topographic masses is enhanced at low altitudes and the sensitivity to the uncertainties in topographic-mass density is increased. The standard value of 2670 kg m−3 for topographic-mass density reductions may introduce an undesired error to geophysically modelled gravitational gradients, for instance, in areas of elevated, low-density sediments (e.g. the Congo basin). In addition, the evaluation of the GOCE-based gravitational gradients at low altitudes amplifies not only the signal, but also their omission errors (in general, the omission error is the signal that has not been modelled, Losch et al.2002).
The second approach is based on the fact that a short-wavelength gravitational signal due to near-surface geology is dampened at satellite altitudes more efficiently than a long-wavelength gravitational signal with lithospheric and sublithospheric origins. Hence, studies on lithospheric or upper-mantle structure interpret the GOCE-based gravitational gradients at the satellite's altitudes (Bouman et al.2013; Fullea et al.2013). In the same way, the omission error of forward-modelled gravitational gradients is dampened efficiently at GOCE satellite's altitudes and these gravitational gradients are then easier to compare with the bandwidth-limited GOCE-based gravitational gradients.
The omission error of the GOCE-based gravitational gradients can be mathematically described and its amplification at the computation point located near the Earth's surface numerically estimated. This motivates the work of this paper. We present a detailed and systematic derivation of the mass-density Green's functions for satellite gravitational gradients in both the spherical-harmonic and closed forms. We demonstrate that these two alternative forms provide a mathematical tool for calculating the omission error of the bandwidth-limited gradiometric data at satellite altitude and its amplification when the GOCE-based satellite gravitational gradients are considered at altitudes near the Earth's surface. The spectral forms of the gravitational gradients are further used to calculate the sensitivity of the satellite gravitational gradients to the depth of density anomalies.
There are, in particular, three main questions addressed in this paper. (i) How important is the omission error of the bandwidth-limited Green's functions when the GOCE-based gravitational gradients are interpreted at different heights above the Earth's surface? (ii) How do the gravitational gradients depend on the spatial distribution of the Earth's density structure? (iii) What is the minimum lateral extent of the portion of the Earth that needs to be considered for regional geophysical models constrained by the GOCE-based satellite gravitational gradients if the full information provided by the Green's functions is to be retained?
It is a known fact that the inverse gravimetric problem, that is the determination of the Earth's density structure from knowledge about the external gravitational potential or vertical gravity, is an ill-posed problem since it violates three Hadamard's criteria on the well-posedness of a mathematical model (Hadamard 1923). There is an extensive literature on dealing with ill-posed inverse problems (e.g. Menke 1989; Parker 1994; Xu 1998; Tarantola 2005). In particular, the inverse gravimetric problem has been studied by for example Pěč & Martinec (1984), Sansò et al. (1986), Matyska (1987), Ballani et al. (1993a,b). The ill-posedness of the inverse gradiometric problem, that is the case where the gravitational gradients are interpreted in terms of the Earth's density structure, is of a similar character as that of the inverse gravimetric problem. However, no theoretical study has been published to investigate the differences in these two inverse problems.
It is not the aim of this paper to provide a comprehensive theoretical study of the inverse gradiometric problem, nor to solve it while determining the Earth's density structure over regional scales. For the latter case, we refer to a number of recent efforts to interpret the GOCE-based gravitational gradients in terms of a regional density model (Braitenberg et al.2011; Álvarez et al.2012; Bouman et al.2013; Fullea et al.2013; Martinec & Fullea 2013). In fact, the difficulties when carrying out these interpretations have motivated the theoretical work presented hereafter.
2 GREEN'S FUNCTIONS FOR THE GRAVITATIONAL POTENTIAL, VECTOR AND GRADIENT TENSOR
3 GREEN'S FUNCTIONS IN SPHERICAL-HARMONIC FORM
3.1 Gravitational potential
3.2 Gravitational vector
3.3 Gravitational gradient tensor
4 CLOSED FORM OF THE ISOTROPIC KERNELS OF THE GRAVITATIONAL VECTOR AND GRADIENT TENSOR
4.1 Gravitational vector
4.2 Gravitational gradient tensor
We continue deriving the closed form of the mass-density Green's function for the gravitational gradient tensor. Eq. (28) implies that to express the kernels Krr(t, x), KrΩ(t, x) and KΩΩ in a closed form, it is sufficient to find the closed form of the rr, rϑ and ϑφ components of the |${\rm grad\, \displaystyle }{\rm grad\, \displaystyle }(1/L)$| tensor.
5 DECOMPOSITION OF THE GRAVITATIONAL GRADIENT TENSOR
6 APPLICATIONS
6.1 Omission error and spatial dependency of the gravitational gradient tensor
Let us first investigate the sensitivity of the gravitational gradients to the height above the Earth's surface where they are computed. The left-hand panel of Fig. 1 shows the isotropic kernels Krr(t, cos ψ) (red solid line), KrΩ(t, cos ψ) (green solid line) and KΩΩ(t, cos ψ) (blue solid line) evaluated by the closed forms (47) for angular distances 0° ≤ ψ ≤ 6° and t = 0.96224, which corresponds to the GOCE perigee height of 255 km. The graphs of the kernels are bell-shaped (terminology used by Eshagh 2011) but with different positions of their maximum values. Function Krr has its maximum at point ψ = 0° and decreases monotonically to zero with increasing angular distance ψ, whereas functions KrΩ and KΩΩ vanish at ψ = 0°, increase their amplitudes for increasing distance ψ, before reaching their maximum amplitudes at ψmax ≈ 1.1° and ψmax ≈ 1.8°, respectively, and then monotonically decreasing with distance ψ. Quantitatively, the amplitude of function Krr decreases to <1 per cent of its maximum value for distances ψ > 3°, while functions KrΩ and KΩΩ still have significant non-zero amplitudes from that distance, while their amplitudes are less than 10−2 × KrΩ(t, ψmax) and 10−2 × KΩΩ(t, ψmax) at distances ψ > 10° and ψ > 18°, respectively.
Projecting these facts into the integrals (50) tells us that the mass density below the computation point, that is the density distribution around the point ψ = 0°, has the largest effect on the |${\bf {\Gamma }}_{rr}$| gravitational gradient , although a negligible effect on the |${\bf {\Gamma }}_{r\Omega }$| and |${\bf {\Gamma }}_{\Omega \Omega }$| gravitational gradients. The density contribution to |${\bf {\Gamma }}_{rr}$| gradually decreases with increasing distance ψ, while the density contribution to |${\bf {\Gamma }}_{r\Omega }$| and |${\bf {\Gamma }}_{\Omega \Omega }$| first increases with increasing distance ψ, reaches its largest contribution at ψmax and then gradually decreases with ψ. The density structure at distances greater than 3° from the computation point may still contribute to |${\bf {\Gamma }}_{rr}$| gravitational gradient by 1 per cent of the contribution from the density structure around the point ψ = 0°. On the contrary, the density structure at distances ψ > 3° up to ψ = 10° and ψ = 18° may contribute to the |${\bf {\Gamma }}_{r\Omega }$| and |${\bf {\Gamma }}_{\Omega \Omega }$| gravitational gradients by amounts that are greater than 1 per cent of the contribution from the density structure around the point ψmax, respectively,
In terms of solving the inverse gradiometric problem for the mass density distribution, the conclusions derived above imply that the |${\bf {\Gamma }}_{rr}$| gravitational gradient provides more localized information on the density distribution than the |${\bf {\Gamma }}_{r\Omega }$| and |${\bf {\Gamma }}_{\Omega \Omega }$| gravitational gradients. Hence, the |${\bf {\Gamma }}_{rr}$| gravitational gradient may be more suitable for solving the inverse problem for density than the other gravitational gradients. In addition, for a chosen computational point, the |${\bf {\Gamma }}_{rr}$| gravitational gradient component contains information about the density structure from rather different regions than the |${\bf {\Gamma }}_{r\Omega }$| and |${\bf {\Gamma }}_{\Omega \Omega }$| components.
The left-hand panel of Fig. 1 also shows the isotropic kernels Krr(t, cos ψ), KrΩ(t, cos ψ) and KΩΩ(t, cos ψ) (dashed lines), computed by summing the series of Legendre polynomials (30) up to cut-off degree jmax = 220, which is equal to the cut-off degree of the GOCO03S satellite gravity model (Mayer-Gürr et al.2012). We can see that the graphs of the truncated isotropic kernels differ very slightly from the graphs of the full-spectrum (that is, the closed form) isotropic kernels. In fact, they are almost indistinguishable within the thickness of the lines. Quantitatively, the largest relative error of the truncated Krr with respect to the closed-form Krr is 0.90 per cent at ψ = 0°. Likewise, the largest relative errors of the truncated KrΩ and KΩΩ are 0.53 per cent and 0.44 per cent at ψmax = 1.1° and 1.8°, respectively.
In geodesy, the term omission error refers to the unresolved or unmodelled part of the Earth's gravity field due to the truncation of the spherical harmonic series used to represent the gravity field. For the following discussion, we make a more generalized use of this term by saying that, the omission error is the signal that is not modelled. Within this terminology, the omission errors of the truncated isotropic kernels Krr, KrΩ and KΩΩ at an altitude of 255 km are , 0.53 and 0.44 per cent, respectively.
When using the GOCE-based gravitational gradients in lithospheric modelling, the isotropic kernels are often not represented in spectral forms (30), but rather the closed form (8) of the mass-density Greens’ function is applied and the volume integral (30) for gravitational gradient tensor |${\bf {\Gamma }}$| is discretized by various geometrical objects (Fullea et al.2008; Braitenberg et al.2011; Uieda et al.2011; Álvarez et al.2012; Hirt et al.2012; Tsoulis 2012; Grombein et al.2013). If the discretization of the integral is performed in a mesh with a fine discretization step Δ such that the Nyquist frequency jN = π/Δ is much higher than the cut-off degree jmax of the GOCE-based model, that is when jN ≫ jmax, then the forward-modelled gravitational gradients will contain a high-frequency (short-wavelength) part that is missing in the GOCE-based gravitational gradients, namely the omission error evaluated in the previous paragraph. Since this error is relatively small at the GOCE satellite's altitude, it can be tolerated for certain types of geophysical applications (Fullea et al.2013). Alternatively, this high-frequency part can be excluded from the modelling of gravitational gradients by replacing the closed-form Green's functions (47) by their bandwidth-limited forms (30).
Complementary to the left-hand panel of Fig. 1, the right-hand panel shows the isotropic kernels for t = 0.99221, which corresponds to a computation-point height of 50 km. The graphs of the full-spectrum isotropic kernels (solid lines) can be described in a similar way as those in Fig. 1 for a computational point at 255-km height. However, the functions decrease faster with increasing angular distance ψ than those in the left-hand panel, which is a well-known fact from potential field theory. Quantitatively, Krr, KrΩ and KΩΩ have amplitudes less than 1 per cent of their maximum amplitudes at distances ψ > 0.6°, ψ > 1.5° and ψ > 2.5°, respectively. Moreover, the amplitudes of the isotropic kernels for a computational point at 50-km height are significantly larger than those for a computational point at 255-km height, meaning that the gravitational gradient signals are amplified when the computation point moves from the GOCE satellite's altitude towards the Earth's surface. This fact can also be deduced from the Meissl scheme that relates surface spherical harmonic coefficients of the Earth's external gravitational potential with coefficients of its first- and second-order derivatives on spheres at different altitudes (Meissl 1971; Rummel & van Gelderen 1995).
What differs substantially are the graphs of the truncated isotropic kernels. Whereas at the computational height of 255 km the full-spectrum and truncated isotropic kernels are coincident within an omission error not larger than 0.9 per cent of the maximum amplitudes, at the computation height of 50 km, the full-spectrum and truncated kernels differ significantly. An advantage of the vertical-vertical Krr kernel with respect to the KrΩ and KΩΩ kernels is that both the full-spectrum and truncated Krr have their maximum at ψ = 0°, but the maximum value is reduced by four times when the kernel Krr is truncated at degree jmax = 220. The difference between the full-spectrum and truncated Krr kernels therefore gives the omission error of the truncated Krr, with the omission error having an amplitude comparable to that of the full-spectrum Krr.
Similarly to Krr, the kernels KrΩ and KΩΩ are reduced in amplitude when omitting their short-wavelength part, but in contrast to Krr, the positions of the maximum amplitudes (denoted by inverted triangles) are shifted towards greater distances ψ from the observer. Shifting the maximum values of the truncated KrΩ and KΩΩ kernels means that the sensitivity of the |${\bf {\Gamma }}_{r\Omega }$| and |${\bf {\Gamma }}_{\Omega \Omega }$| gravitational gradients is transferred to the density structure at locations different from those for the original, full-spectrum isotropic kernels. The difference between the full-spectrum and truncated KrΩ and KΩΩ gives the omission error of the truncated kernels. We can therefore see that the omission error of the bandwidth-limited gravitational gradients increases when the computation point moves from GOCE satellite's altitude towards the Earth's surface.
Both the closed and spherical-harmonic forms of the Green's functions allow a direct interpretation in terms of the minimum lateral extent of the Earth that needs to be considered when a regional model is constrained by gravitational gradients, if the full information provided by the Green's functions is to be retained. For instance, the left-hand panel of Fig. 1 shows that, at a height of 255 km above the Earth's surface, Krr could be safely modelled in a region with a lateral extent of 3°, that is Krr has been reduced to almost zero by this distance, but clearly that would be insufficient for the vertical-horizontal, and horizontal-horizontal Green's functions, where a substantial proportion of both KrΩ and KΩΩ are still present by this distance.
To develop and quantify this idea, we define an angular distance ψ0 at which the ratio between the amplitude of the Green's functions to the corresponding maximum amplitude is equal to δ. Fig. 2 shows the angular distance ψ0 as a function of the height at which the gradients are being considered for the closed-form isotropic kernels expressed by eq. (47) (solid lines), and the bandwidth-limited isotropic kernels computed by summing the series of Legendre polynomials (30) truncated at degree jmax = 220 (dashed lines). Two examples of ratios between the amplitude of an isotropic kernel at angular distance ψ0 to its maximum amplitude are given: δ = 1 per cent (left-hand panel) and δ = 10 per cent (right-hand panel). We can make three observations. First, the lateral size ψ0 for the closed-form kernels (nearly exactly) linearly increases with increasing computation-point height. This hidden property of closed forms (47) may be used as a rule of thumb for approximately choosing the minimum lateral size of a regional model when it is constrained by the gravitational gradients. Due to the truncation, the lateral size ψ0 for the truncated kernels slightly oscillates around, or is nearby, the closed-form kernels, meaning that there is no significant difference in choosing ψ0 for the truncated kernels. However, the positions of the maximum values of the truncated kernels differ significantly from those of the closed-form kernels for computation-point heights near the Earth's surface, as demonstrated for the height of 50 km in the right-hand panel of Fig. 1. The second observation confirms the results shown in Fig. 1. For a given computation-point height, the Krr kernel requires the smallest lateral extent. This extent is greater for kernel KrΩ and is the greatest for kernel KΩΩ. For instance, at a height of 100 km and for δ = 1 per cent, the lateral size ψ0 for the closed-form Krr, KrΩ and KΩΩ are equal to 1.1°, 4.0° and 7.2°, respectively. Finally, comparing the two panels of Fig. 2, we can see that the minimum lateral size ψ0 is smaller if the full information provided by the kernels is allowed to be reduced by 10 per cent rather than by 1 per cent. For instance, at a height of 100 km and for δ = 10 per cent, the lateral extent ψ0 for the closed-form Krr, KrΩ and KΩΩ are equal to 0.8°, 2.0° and 3.1°, respectively.
For the identification of geological units in unexplored parts of the world, it is advantageous to interpret the GOCE-based gravitational gradients at lower altitudes since the gravitational gradient signals are amplified and better reflect near-surface geological structures. Since the omission error for a near-surface computation point is also significantly amplified, the forward geophysical modelling must ensure that the omission error of the forward-modelled gravitational gradients is the same as the omission error of the GOCE-based gravitational gradients. Only after this are the GOCE-based gravitational gradients able to be interpreted in terms of geological structure. One way of performing this step is to pass the forward-modelled gravitational gradients through the bandpass filter with the same bandwidth as the GOCE-based gravitational gradients (e.g. Martinec 1991). This approach has been applied by Martinec & Fullea (2013) when interpreting the GOCO03S gravity model to refine a model of sedimentary rock cover in the southeastern part of the Congo basin.
6.2 Sensitivity of the gravitational gradients to the depth of the mass anomalies
7 SUMMARY AND CONCLUSIONS
This paper has been motivated by an effort to create an adequate mathematical tool for analysing the sensitivity of satellite gradiometric data, as provided by the GOCE mission, to the Earth's internal density structure and to the height where the GOCE-based gravitational gradients are considered. It has been shown that an approach based on Green's functions is highly convenient for carrying out such an analysis. We express the Green's functions for gravitational gradients in the spectral form as a series of tensor spherical harmonics. This form of the Green's functions can be used for representing the forward-modelled gravitational gradients when the GOCE-based gravitational gradients are interpreted at different heights above the Earth's surface. Alternatively, by means of the addition theorems for tensor spherical harmonics, the spectral forms are converted to closed spatial forms. These two alternative forms provide a powerful tool for analysing the omission error of the bandwidth-limited GOCE-based gravitational gradients at the satellite's altitude and its amplification at lower altitudes. We show that the omission error of the bandwidth-limited Green's functions for gradiometric data at the GOCE satellite altitude does not exceed 1 per cent in amplitude when compared to the full-spectrum Green's functions. Such an error may be tolerated for certain types of geophysical applications (Bouman et al.2013; Fullea et al.2013). However, when computing the bandwidth-limited gravitational gradients close to the Earth's surface, the omission error is significantly enhanced. We conclude that the GOCE-based gravitational gradients at altitudes close to the Earth's surface could be interpreted in terms of the Earth's density structure only after filtering the short-wavelength content of the forward-modelled gravitatinal gradients, such that the omission error of the GOCE-based gravitational gradients is equal to the omission error of the forward-modelled gravitational gradients.
The spectral forms of the gravitational gradients additionally allow us to calculate the sensitivity of the satellite gravitational gradients to the depth of a density anomaly, expressed in terms of the harmonic degree of the internal density anomaly. We show that the largest gravitational gradient response is obtained for shallow mass anomalies, and is further amplified as the harmonic degree increases. That is why the GOCE gravitational gradients are able to improve the modelling of the Earth's crustal and lithospheric structure. In addition, if they are evaluated near the Earth's surface, the gravitational gradient signal is amplified and better reflects the near-surface geological structure, which may be an advantage for geological mapping in less well-surveyed regions.
The author would like to thank Roger Haagmans for directing attention to the problem studied in this paper, Kevin Fleming and two anonymous reviewers for their comments on the manuscript that have helped to improve the manuscript. This work has been carried out within the framework of the ESA sponsored GOCE-GDC study and the research program 11/RFP.1/GEO/3309 financed by the Science Foundation of Ireland. The author acknowledges this support.