Printer Friendly

An alternative method to improve gravity field models by incorporating GOCE gradient data/Metodo alternativo para mejorar los modelos de campo gravitacional al incorporar informacion del satelite Explorador de la Circulacion Oceanica y de Gravedad.


The Gravity field and steady-state Ocean Circulation Explorer (GOCE) satellite (European Space Agency, 1999) was launched on March 17th, 2009 and was terminated on November 11th, 2013. By combining the gradient and orbit data from GOCE, three groups of gravity field models with increasingly high accuracies have been derived by the ESA (European Space Agency) GOCE level 2 processing facility. Indeed, in addition to GOCE data, a number of other kinds of observations can also be used to derive gravity field models, including observations from altimeter satellites (Marchenko, 2003a; Flechtner et al., 2006), the CHAllenging Minisatellite Payload (CHAMP) (Reigber et al., 1996) and the Gravity Recovery and Climate Experiment (GRACE) (Tapley et al., 2004), as well as some ground observations (Marchenko, 2003a; Pavlis et al., 2012). The use of such data has enabled a range of institutes to resolve more than 100 gravity field models, released by the International Centre for Global Earth Models operated by the German Research Centre for Geosciences. In comparison to other kinds of gravity information, GOCE data are of higher accuracy in some special frequency bands, and the spatial consistency of the accuracy is superior to that of ground observations (Marchenko et al., 2016). Nevertheless, deriving an accurate gravity field model using GOCE gravity gradient observations in combination with other data remains an important issue (Gilardoni et al., 2013).

One possible approach to this question is to initially construct a series of normal equations based on all the observations (Pail et al., 2010) and then resolve them using least square estimation method (LSM). However, as the large quantity of data makes the computational cost of this process prohibitive, an alternative method is to utilize GOCE gravity gradient data to improve existing gravity field models that have been derived from other kinds of observations. If this approach is adopted, there is no need to process the other (non-GOCE) observation data over again.

In this study, we adopt the latter method in which the existing gravity field models are selected as reference models to be improved. To resolve the gravity field model, two approaches can be used, i.e., LSM and spherical harmonic analysis. Currently, the former is more commonly applied for processing GOCE gradient data for gravity field recovery, for example, the resolving of time-wise solution (Pail et al., 2011) and direct solution of the GOCE gravity field models (Bruinsma et al, 2010) by ESA DGM-1S (Hashemi et al, 2013), etc. One of the main advantages of the LSM approach is the ease with which a gravity field model can be resolved using gradient and orbit data in combination. Considering that the existing gravity field models are improved mainly by gradient data in this study, we adopt spherical harmonic analysis to recover the gravity field model.

Resolving of a gravity field model via the use of spherical harmonics can be treated as a Boundary Value Problem (BVP) (Rummel, et al, 1993), the first step of which is to construct a so-called boundary value condition based on the observation equation. A great deal of research has been conducted in this area. Rummel and van Gelderen (1992) analyzed the spectral characteristics of different gravity gradient components and generated a series of two-dimensional Fourier expressions, van Gelderen and Rummel (2001) summarized the solution of a general geodetic BVP under gravity gradient boundary conditions. Martinec (2003) discussed Green function to gradiometric BVP. In addition, invariants of gradient tensors (Rummel, 1986, Sacerdote & Sanso, 1989, Vermeer, 1990, Marchenko, 2003b, Baur et al, 2008) were discussed to recover gravity field models by LSM. Holota (1989) and Yu and Zhao (2010) constructed the boundary conditions using invariants of a gravity gradient tensor. Those conditions contain few attitude errors but create linear errors. Marchenko et al.(2016) investigated a special version of the space-wise method based on the second degree radial derivatives of GOCE missions. Their results show that the second order derivatives of EGG_TRF_2 data can improve the middle and long waves of the geo-potential (Marchenko et al, 2016).

The radial gravity gradients are utilized in this study to construct the boundary condition. It is noteworthy that Fuchs et al. (2013) proposed a superior and highly accurate method for vertical gravity gradient computation, which they successfully applied to a discussion of the Tokyo earthquake. By this method, Bouman et al.(2016) calculated satellite gravity gradient grids for geophysics.

Meng et al.(2016) proposed an alternative combination of [V.sub.xx], [V.sub.w] and [V.sub.zz] to compute a new [V.sub.zz] which then is used for computing the radial gravity gradient. Since accurate radial gravity gradients can be obtained, we believe that they can be used as a boundary condition to improve some existing gravity field model and to obtain new ones. This is the starting point of this paper.


Radial gravity gradient computation

GOCE measurements include the [V.sub.xx], [V.sub.yy], and [V.sub.zz] components of the gravity gradient tensors (Stummer et al, 2012). Thus, if no noise is observed, the Laplace equation can be applied, as follows:

[V'.sub.zz] = -([V.sub.zz] + [V.sub.yy]) (1)

This leads us to obtain two vertical gravity gradient values (i.e., [V.sub.zz] and [V'.sub.zz]) that are, however, unequal to one another because of noise. Fuchs et al. (2013) addressed this by using mean values of [V.sub.zz] and [V'.sub.zz] to improve the accuracy of vertical gravity gradients. Meng et al. (2016) proposed the following combination to compute new vertical gravity gradients,

[mathematical expression not reproducible] (2)

Compared to [V.sub.zz], noises of [V.sub.zz,c] are reduced about 42.3% (Meng et al, 2016). However, [V.sub.zz,c] are the vertical components of gravity gradient tensors in Gradiometer Reference Frame (GRF) and not the radial gravity gradients, hence the coordinate transformation is needed. Although Fuchs et al. (2013) and Meng et al. (2016) proposed new methods for computing new vertical gravity gradients, horizontal gravity gradiens (ie, [V.sub.xx] and [V.sub.yy]) are also needed for the data rotation. However, no literature discussed how to compute high accurate [V.sub.xx] and [V.sub.yy]. Indeed, according to the Laplace equation, we can construct the following equations,

[mathematical expression not reproducible] (3)

This leads us to obtain two x-components (i.e., [V.sub.xx] and [V.sub.xx]) and two y-components ( i.e., [V.sub.yy] and [V'.sub.yy]) of a gravity gradient tensor. We provide combinations for these values by solving the following equations:

[mathematical expression not reproducible] (4)

[mathematical expression not reproducible] (5)

In these expressions, a and b are the combination factors used for computing the new x-component of the gravity gradient tensor, c and d are the combination factors used for computing the new y-component of the gravity gradient tensor, while [[delta].sup.2.sub.ii] (i = x, y, z) are the variance of [V.sub.ii] noises, and [[delta].sup.2.sub.ii,jj](i = x, y, z) refer to the co-variances between different [V.sub.ii] noises. The solutions of these expressions are therefore as follows:

[mathematical expression not reproducible] (6)

[mathematical expression not reproducible] (7)

Thus, as long as [[delta].sup.2.sub.ii,jj] (i = x, y, z) are known, it is possible to obtain optimal factors from equation (6) and (7). Indeed, if the components of the gradient tensor are independent, equation (6) and (7) can be simplified as follows:

[mathematical expression not reproducible] (8)

[mathematical expression not reproducible] (9)

The error distribution of the GOCE gradient tensor applied by Fuchs et al. (2013) is used here, as follows:

[[[delta].sub.zz] = [[delta].sub.xz] = 2[[delta].sub.xx] = 2[[delta].sub.yy]=2[[sigma].sub.0]] (10)

Solutions for a and b, c and d can therefore be derived from equation (8) and equation (9), as follows:

[mathematical expression not reproducible] (11)

[mathematical expression not reproducible] (12)

[mathematical expression not reproducible] (13)

Hence, [mathematical expression not reproducible] and [mathematical expression not reproducible]. The noises of horizontal gravity gradients are reduced by about 8.9%. In this paper, besides [V.sub.xx,c] and [V.sub.yy,c] expressed by equation (13), [V.sub.zz,c] expressed by equation (2) are used to compute the radial gravity gradients by the following coordinate transformation,

[mathematical expression not reproducible] (14)

where R denotes the transformation matrix from the GRF coordinate system to the Local North Oriented Frame (LNOF) (Gruber et al., 2010) coordinate system. In LNOF the Z direction is same as the radial direction of the Earth. It should be noted that, owing to the poor accuracy of [V.sub.xy] and [V.sub.yz], we replace them with [U.sub.xy] and [U.sub.yz] respectively, that were calculated using the reference model (i.e., the model to be modified). Since accuracies [V.sub.xx,c], [V.sub.yy,c], and [V.sub.zz,c] are higher than those of [V.sub.xx], [V.sub.yy], and [V.sub.zz], this means that the accuracy of [V.sub.rr] can be improved via the use of [V.sub.xx,c], [V.sub.yy,c], and [V.sub.zz,c].

Residual gradient data were used in this study for data rotation as follows:

[mathematical expression not reproducible] (15)

where f and [T.sub.xz] denote the respective differences between [V.sub.xx],[V.sub.yy], [V.sub.zz], [V.sub.xz] and [U.sub.xx], [U.sub.yy], [U.sub.zz], [U.sub.xz] which are obtained from the reference model. [T.sub.xy] and [T.sub.xz] are both equal to zero because [V.sub.xy] and [V.sub.yz] are both replaced by [U.sub.xy] and [U.sub.yz], respectively. The use of this rotation approach enables us to obtain [T.sub.rr].

The variable [T.sub.rr], obtained from equation (15), was then used to construct the following boundary value problem (Wan & Yu, 2013), as follows:

[mathematical expression not reproducible] (16)

where/denotes [T.sub.rr], and S refers to the mean surface of GOCE orbits. Thus, the solution of the boundary value problem (16) provides the spherical harmonic coefficients of T at satellite altitude. By solving this BVP, we can obtain solutions of ([DELTA][C.sub.nm], [DELTA][S.sub.nm]), as follows:

[mathematical expression not reproducible] (17)

where GM denotes the gravitational constant of the Earth, a is the semi-major axis of the Earth, (r, [theta], [lambda]) denote the spherical coordinates, and [P.sub.nm] (cos [theta]) is the fully normalized Associated Legendre Functions with degree n and order m. The final solutions of the gravity field model are as follows:

[mathematical expression not reproducible] (18)

where ([C.sub.nm,0], [S.sub.nm,0]) denote the coefficients of the initial model, i.e., EIGEN5C (Forste et al., 2008) applied in this study. ([C.sub.nm], [S.sub.nm]) are the coefficients of the final model.

Noise processing of GOCE gradient data

Before using GOCE gradient observations to recover a gravity field model, the observed noise should be processed first (Schuh, 2003; Reguzzoni & Tselfes, 2009; Bruinsma et al, 2010; Pail et al, 2011; Gatti et al, 2014). Three types of noise should be considered. The first one is low frequency noise. As implied by the name, this type of noise is located in the low frequency band (i.e., lower than 0.005 Hz in this study) and constitutes the main type of noise in GOCE gradient data. This low frequency noise is large in magnitude, and within a known frequency band. The second type of noise is random noise. This noise is smaller in magnitude and is distributed throughout the full frequency band. The combinations presented in equations (2) and (13) are used to reduce this random noise. Hence, the techniques of white noise processing are not discussed here. The third one is gross error. Their magnitude is very large and its distribution is scattered. Gross error is found in all forms of observational data used in the gravity field recovery, including gradient, orbit, and attitude data. Gross error can also be created by a kind of computing, e.g., the warm effect of filtering.

a. Time domain processing

Forward and backward Finite Impulse Response (FIR) filtering was initially applied to the residual gravity gradient data, which are the difference values of the actual GOCE gradient data and the simulated gradient data from the reference gravity field model (i.e., the reference model to be modified). As noted, EIGEN5C model was selected here as the reference model, which was derived by GFZ and CNES (National Centre for Space Studies) using five years of GRACE data as well as some ground observations. We only selected 0-300 degrees6 / orders of EIGEN5C as the reference model, considering the limits of measurement band wise (MBW) and the accuracy of the gradient observations.

To show the filtering effect, we processed the gradient data provided by EGG_NOM2 from ESA over the short time interval between 0:0:1.6 on October 17th, 2010 and 23:59:59.6 on October 26th, 2010, using a band-pass FIR filter. The pass band frequency is 0.005~0.1 Hz. Fig. 1 shows the Power Spectral Density (PSD) of the data after filtering. It is clearly that signals in the pass band remained unchanged by this filtering processing, while those beyond the pass band were largely removed and only occupy a very small proportion in the final data obtained from the filtering. This filtering step therefore enabled us to remove low frequency noise while saving the signals during the pass band. It is also noteworthy that no phase drift occurred (e.g., Yu & Wan, 2011); thus, the data after filtering can still be considered as gravity gradient signals.

In order to verify the effects of FIR processing, we perform a statistical analysis to the absolute values of Laplace Equation. Table 1 gives the results. Smaller values of this statistic indicate a higher accuracy; thus, according to Table 1, the accuracy was improved by a factor of 159172 by FIR filtering in terms of the mean value.

b. Space domain processing

Although the dominant noises have been removed in the time domain via FIR filtering, we still need to remove the residual gross errors. These errors were processed in the space domain. When the radial gradients of the disturbing gravity were obtained using equation (15), we were able to perform grid processing on the satellite orbit surface. We first divided the surface into girds, such as 20' x 20' and then conducted statistical processing. Mean values, variances, and standard errors of the data in each grid were then obtained, and we removed the data that deviated from the mean value by more than twice the standard deviation. This processing led to the removal of approximately 3.84% of the data.

Processing was carried out in the space domain because gross errors are sometimes stable over short time periods, perhaps because of an abnormal satellite situation, making these errors difficult to remove in the time domain. However, even though the data in the same grid are from different points in time, their values should be the same or close to each other. If some abnormal values exist, there must be gross errors that can be isolated.

Results and Discussion

All the gradient data obtained by GOCE throughout the lifetime of the project are processed in this paper. As described above, we first processed the noises in both time and space domains. Radial disturbing gravity gradients were then calculated using equation (15). In this section, the EIGEN5C is taken as an example to be improved with the proposed method. The final solution is derived by integration based on the orthogonality of spherical harmonic functions given by equation (16) and equation (17). The new model is named as the Earth Gravitational Model of Gravity field and Ocean Circulation Explorer (EGMGOCE).

Comparisons in spherical domain

The degree variances of EGMGOCE and EGM08 (Pavlis et al, 2012), are compared in Fig. 2 and the geoid heights are compared in Fig. 3.

In Fig. 2, Variance_EIGEN5C denotes the degree variance differences between EGM08 and EIGEN5C, and Variance_EGMGOCE denotes the differences between EGM08 and EGMGOCE. In Fig. 3, EIGEN5C-EGM08 and EGMGOCE-EGM08 indicate the cumulative geoid differences between EIGEN5C, EGMGOCE and EGM08 respectively. These two figures illustrate that the initial model was mainly improved in the degrees from 100 to 200. If EGM08 is considered as the true model, the geoid accuracy of the new model is higher than that of the initial model by 4 cm up to 200 degree.

For a comparison with GOCE models, EGMD5 from ESA is selected, which is the newest GOCE gravity field model derived by the direct method. The degree variance differences and cumulative geoid differences between EGMGOCE, EIGEN5C and EGMD5 are shown in Fig. 4 and Fig. 5, respectively. The new model is clearly much closer to EGMD5 than EIGEN5C, particularly at degrees from 80 to 230. If EGMD5 is considered as the true model, the geoid accuracy of the new model is higher than that of the initial model by 7 cm up to 200 degree. These results indicated that the extraction of the GOCE gravity gradient signals and modification of EIGEN5C were successful.

Comparisons in space domain

The differences between EGMGOCE and EIGEN5C can be seen as the contribution from GOCE gradient observation. In order to present GOCE contribution in space domain, we selected 0-200 degrees / orders of EIGEN5C, EGMGOCE, EGM08 and EGMD5 to compute geoid values at mean surface of the Earth. The data were calculated at grids of 1[degrees] x 1[degrees]. Statistics of their differences are given in Table 2 and Table 3. According to Table 2, if EGM08 is the true model, the accuracy of EIGEN5C is improved by 4.6 cm because of the use of GOCE gradient data. According to Table 3, if EGMD5 is the true model, the accuracy of EIGEN5C is improved by 7 cm. In order to show the differences of the geoids between EIGEN5C and EGMGOCE, Fig. 6 and Fig. 7 show the geoids calculated by EIGEN5C and EGM08 respectively, and Fig. 8 shows their differences. Table 4 gives the statistics of absolute values of the differences. According to this table, mean value of the difference is about 5 cm and the standard deviation is about 9 cm. According to Fig. 8, large differences mainly occur in Africa, South America, Antarctica and Himalaya (Yi, 2001; Rummel &Gruber, 2012). This indicates that the data used to derive EIGEN5C has a little poorer accuracy in these areas, since accuracy of GOCE gradient data has a high space consistence, so that accuracy of EIGEN5C can be improved by GOCE gradient observations.


This paper discusses the method for modifying existing gravity field models using GOCE gradient observations. We used combinations of [V.sub.xx]. [V.sub.yy] and [V.sub.zz] to compute new horizontal and vertical gravity gradients, which can reduce white noises of them by 8.9%, 8.9% and 42.3% respectively. Compared with the conventional method, accuracies of the radial gravity gradients can be improved correspondingly.

Strategies for processing noise within GOCE gradient observation data in both time and space domains are also discussed in this study. Band pass FIR filtering is adopted in the time domain and grid processing based on the statistical properties of each grid is implemented in the space domain. These strategies enable the removal of low frequency noises and gross errors. The combinations given by equations (2) and (13) can reduce white noises.

The EIGEN5C model was improved as an example in this study using all available GOCE gradient data, and a new gravity field model named as EGMGOCE was derived. Compared with EGM08, modifications using GOCE gradient data are mainly reflected in the degrees from 100 to 200. Up to the 200th degree, the cumulative geoid difference is reduced by 4cm. Compared with EGMD5, modifications from GOCE data are mainly reflected in degrees from 80 to 230, and up to the 200th degree, the cumulative geoid difference is reduced by 7 cm. These results showed that the GOCE gradient observations can contribute meaningfully to the development of EGMGOCE model and demonstrate the effectiveness of the methods adopted in this paper. Other gravity field models can also be processed similarly to derive new models.

Since the radial gravity gradients as well as the horizontal components of gravity gradients can be derived with high accuracy, they could not only be utilized to recover gravity field models, but also be used to extract and interpret geophysical signals directly, such as Moho geometry (Ye et al., 2016), earthquake (Fuchs, et al., 2013), etc. Hence, the methods proposed in this paper can also be powerful tools for researches in geophysics.


We thank the ESA for providing GOCE data. We also thank the reviewers for their helpful suggestions. This paper is funded by the National Natural Science Foundation of China (Nos. 41674026,41404019) and the open fund of Key Laboratory of Space Utilization, Chinese Academy of Sciences (CSU-WX-A-KJ-2016-044).


Baur, O., Sneeuw, N., & Grafarend, E.W. (2008). Methodology and use of tensor invariants for satellite gravity gradiometry. Journal of Geodesy, 82(4-5), 279-293

Bouman, J., Ebbing, J., Fuchs, M., Sebera, J., Lieb, V., Szwillus, W.. Haagmans, R., & Novak, P. (2016). Satellite gravity gradient grids for geophysics. Scientific Reports, 6, 21050. DOI:10.1038/srep21050

Bruinsma, S. L., Marty, J. C., Balmino, G., Biancale, R., Foerste, C., Abrikosov, O., & Neumayer, H. (2010). GOCE Gravity Field Recovery by Means of the Direct Numerical Method. ESA Living Planet Symposium, 27th June - 2nd July 2010, Bergen, Norway.

European Space Agency. (1999). Gravity Field and Steady-State Ocean Circulation Mission--The Four Candidate Earth Explorer Core Missions. Technical Report, ESA SP-1233(1), European Space Agency Publications Division, Noordwijk, The Netherlands.

Forste, C., Flechtner, F., & Schmidt, R. (2008). EIGEN-GL05C--a new global combined high-resolution GRACE-based gravity field model of the GFZ-GRGS cooperation. In: General Assembly European Geosciences Union, Vienna, Austria 2008 (Geophysical Research Abstracts, vol. 10, Abstract No. EGU2008-A-06944).

Flechtner, F., Schmidt, R., Meyer, U., Schoene, T., Esselborn, S., Foerste, C., & Rothacher, M. (2006). The benefit of EIGEN gravity field models for altimetry and vice versa. Proceedings of the Symposium on 15 Years of Progress in Radar Altimetry, by Danesy, D. ISBN:92-9092-925-l. Noordwijk, Netherlands: European Space Agency.

Fuchs, M., Bouman, J., Broerse, T., Visser, P., & Vermmersen, B. (2013). Observing coseismic gravity change from the Japan Tohoku-Oki 2011 earthquake with GOCE gravity gradiometry. Journal of Geophysical Research, 118, 5712-5721. DOI:10.1002/jgrb.50381

Gatti, A., Reguzzoni, M., Migliaccio, F., & Sanso, F. (2014). Space-wise grids of gravity gradients from GOCE data at nominal satellite altitude. Presented at the 5th International GOCE User Workshop, 25-28 November 2014, UNESCO, Paris, France

Gilardoni, A., Reguzzoni, M., Sampietro, D., & Sanso, F. (2013). Combing EGM2008 with GOCE gravity field models. Bollettino di Geofisica Teorica ed Applicata, 54(4), 285-302.

Gruber, T., Rummel, R., & Abrikosov, O. (2010). GOCE Level 2 Product Data Handbook, ESA.

Hashemi, F., Ditmar, P., & Klees, R. (2013). The static gravity field model DGM-1S from GRACE and GOCE data: computation, validation and an analysis of GOCE mission's added value. Journal of Geodesy, 87(9), 843-867.

Holota, P. (1989). Boundary Value Problems and Invariants of the Gravitational Tensor in Satellite Gradiometry. Lecture Notes in Earth Sciences, 25, 447-457.

Marchenko, A. N. (2003a). Regional gravity field model from satellite altimetry, gravimetry, and GPS-leveling data in the Black Sea area. Bolletino di Geodesia et Scienze Affini, LXII(4), 246-259.

Marchenko, A. N. (2003b). A note on the eigenvalue--eigenvector problem. In: N. Kiihtreiber (Eds,) Festschrift dedicated to Helmut Moritz on the occasion of his 70th birthday (143-154). Institute for Geodesy, Graz University of Technology.

Marchenko, A. N., Marchenko, D. A., & Lopushansky, A. N. (2016). Gravity field models derived from the second degree radial derivatives of the GOCE mission: a case study. Annals of Geophysics, 59(6), 1-11.

Martinec, Z. (2003) Green's function solution to spherical gradiometric boundary-value problems. Journal of Geodesy, 77(1-2), 41-49, DOI: 10.1007/s00190-002-0288-z.

Meng, X. C., Wan, X. Y., Yu, J. H., Zhu, Y C., & Feng, W. (2016). Optimization method for computing radial gravity gradient using gravity gradient observations. Acta Geodaetica et Cartographica Sinica,45(7), 775-781.

Pail, R., Bruinsma, S., Migliaccio, F., Forste, C., Goiginger, H., Schuh, W. D., Hock, E., Reguzzoni, M., Brockmann, J. M., Abrikosov, O., Veicherts, M., Fecher, T., Mayrhofer, R., Krasbutter, I., Sanso, F., & Tscherning, C. C. (2011). First GOCE gravity field models derived by three different approaches. Journal of Geodesy, 85, 819-843. DOI: 10.1007/s00190-011-0467-x.

Pail, R., Goiginger, H., Schuh, W. D., Hock, E., Brockmann, J. M., Fecher, T., Gruber, T., Mayer-Gurr, T., Kusche, J., Jaggi, A., & Rieser, D. (2010). Combined satellite gravity field model GOCO01S derived from GOCE and GRACE. Geophysical Research Letters, 37, L20314. DOI: 10.1029/2010GL044906.

Pavlis, N. K., Holmes, S. A., Kenyon, S. C., & Factor, J. K. (2012). The development and evaluation of the earth gravitational model 2008 (EGM2008). Journal of Geophysical Research-Solid Earth, 111, B04406. DOI: 10.1029/2011JB008916.

Reguzzoni, M., & Tselfes, N. (2009). Optimal multi-step collocation: application to the space-wise approach for GOCE data analysis. Journal of Geodesy, 83(1), 13-29. DOI: 10.1007/s00190-008-0225-x.

Reigber, C., Bock, R., Forste, C., Grunwaldt, L., Jakowski, N., Liihr, H., Schwintzer, P., & Tilgner, C. (1996). CHAMP Phase B: Executive Summary; Scientific Technical Report STR96/13; GeoForschungsZentrum Potsdam.

Rummel, R. (1986). Satellite gradiometry. Lecture Notes in Earth Sciences, 7, 317-363.

Rummel, R., & van Gelderen, M. (1992). Spectral-analysis of the full gravity tensor. Geophysical Journal International, 111(1), 159-169.

Rummel, R., van Gelderen, M., Koop, R., Schrama, E., Sanso., F., Brovelli, M., Miggliaccio, F., & Sacerdote, F. (1993). Spherical harmonic analysis of satellite gradiometry. Report no. 39, Publications on Geodesy, New Series, Netherlands Geodetic Commission, Delft, The Netherlands.

Rummel, R., & Gruber, T. (2012). Product Acceptance Review (PAR)--Level 2 Product ReportRep., ESA.

Sacerdote, F., & Sanso, F. (1989). Some problems related to satellite gradiometry. Bulletin of Geodesy, 63, 405-415.

Schuh, W. D. (2003). The processing of band-limited measurements: filtering techniques in the least squares context and in the presence of data gaps. Space Science Reviews, 108(1-2), 67-78.

Stummer, C., Siemes, C., Pail, R., Frommknecht, B., & Floberghagen, R. (2012). Upgrade of the GOCE level lb gradiometer processor. Advances in Space Research, 49(4), 739-752. DOI: 10.1016/j.asr.2011.11.027.

Tapley, B. D., & Bettadpur, S. (2004). The gravity recovery and climate experiment: Mission overview and early results. Geophysical Research Letters, 31(9), 1-4.

Vermeer, M. (1990). Observable quantities in satellite gradiometry. Journal of Geodesy, 64(4), 347-361.

Van Gelderen, M., & Rummel, R. (2001). The solution of the general geodetic boundary value problem by least squares. Journal of Geodesy, 75(1), 1-11.

Wan, X. Y., & Yu, J. H. (2013). Derivation of the radial gradient of the gravity based on non-full tensor satellite gravity gradients. Journal of Geodynamics, 66, 59-64.

Ye, Z., Tenzer, R., Sneeuw, N., Liu, L., & Wild-Pfeiffer, F. (2016). Generalized model for a Moho inversion from gravity and vertical gravity-gradient data. Geophysical Journal International, 207(1), 111-128. DOI:

Yi, W. Y (2012). An alternative computation of a gravity field model from GOCE. Advances in Space Research, 50(3), 371-384.

Yu, J. H., & Wan, X. Y (2011). The frequency analysis of gravity gradients and the methods of filtering processing. In: Proceedings of '4th International GOCE User Workshop', Munich, Germany (ESA SP-696, July 2011).

Yu, J. H., & Zhao, D.M. (2010). The gravitational gradient tensor's invariants and the related boundary conditions. Science in China Series D-Earth Sciences 53(5), 781-790.

Xiaoyun Wan (1) & Jiangjun Ran (2*)

(1) School of Land Science And Technology, China University of Geosciences (Beijing)

(2*) Delft University of Technology, Stevinweg 1, 2628 CN, Delft, The Netherlands

(*) Corresponding author:


Manuscript received: 07/05/2017

Accepted for publication: 25/05/2018

How to cite item

Wan, X., & Ran, J. (2018). An alternative method to improve gravity field models by incorporating GOCE gradient data, Colombia. Earth Sciences Research Journal, 22(3), 187-193.

Table 1. Laplace Equation statistics

              Raw data      FIR

Number        863999        858500
Min(mE)       789524.4           2.036E-05
Max(mE)       800881.4         957.125
Mean (mE)     795860.75          5.617
Std (mE)        2128.44          5.334

Table 2. Statistical geoid differences between EIGEN5C, EGMGOCE and

               EIGEN5C-08        EGMGOCE-08

Num            64800              64800
Sum(m)            72.4087           122.1064
Minimum(m)        -3.9080            -2.3716
Maximum(m)         3.2736             3.0994
Mean(m)            0.0011             0.0019
Std(m)             0.2000             0.1538

Table 3. Statistical geoid differences between EIGEN5C, EGMGOCE and

             EIGEN5C-D5   EGMGOCE-D5

Num          64800        64800
Sum(m)        -906.3277    -856.6300
Minimum(m)      -3.6455      -1.7668
Maximum(m)       3.8129       2.0641
Mean(m)         -0.0140      -0.0132
Std(m)           0.2272       0.1570

Table 4. Statistical absolute values of geoid differences between

Num     Min (m)    Max (m)   Mean (m)   Std (m)

64800   4.40E-07   2.2011    0.0523     0.0911
COPYRIGHT 2018 Universidad Nacional de Colombia, Departamento de Geociencias
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2018 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:GEODESY; Gravity field and steady-state Ocean Circulation Explorer
Author:Wan, Xiaoyun; Ran, Jiangjun
Publication:Earth Sciences Research Journal
Article Type:Report
Date:Sep 1, 2018
Previous Article:Geospatial assessment of eco-environmental changes in desertification area of the Brazilian semi-arid region/Evaluacion geoespacial de los cambios...
Next Article:Empirical Correlation between Geotechnical and Geophysical Parameters in a Landslide Zone (Case Study: Nargeschal Landslide)/Correlacion empirica...

Terms of use | Privacy policy | Copyright © 2021 Farlex, Inc. | Feedback | For webmasters |