冰川均衡调整对南极冰质量平衡监测的影响及其不确定性

本文研究了新的全球冰川均衡调整(GIA)模型对南极冰盖质量平衡监测的影响,考虑现有冰川负荷模型和地幔黏滞度模型的差异,完整评估了结果的不确定性,最后结合GRACE和卫星测高的结果进行了对比分析.结果表明,GIA对GRACE监测的等效水柱变化有重大影响,较大的GIA影响出现在西南极,沿罗斯冰架-卡姆布冰流-罗尼冰架-南极...     

Assessment of terrestrial water contributions to polar motion from GRACE and hydrological models

The hydrological contribution to polar motion is a major challenge in explaining the observed geodetic residual of non-atmospheric and non-oceanic excitations since hydrological models have limited input of comprehensive global direct observations. Although global terrestrial water storage (TWS) estimated from the Gravity Recovery and Climate Experiment (GRACE) provides a new opportunity to study the hydrological excitation of polar motion, the GRACE gridded data are subject to the post-processing de-striping algorithm, spatial gridded mapping and filter smoothing effects as well as aliasing errors. In this paper, the hydrological contributions to polar motion are investigated and evaluated at seasonal and intra-seasonal time scales using the recovered degree-2 harmonic coefficients from all GRACE spherical harmonic coefficients and hydrological models data with the same filter smoothing and recovering methods, including the Global Land Data Assimilation Systems (GLDAS) model, Climate Prediction Center (CPC) model, the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis products and European Center for Medium-Range Weather Forecasts (ECMWF) operational model (opECMWF). It is shown that GRACE is better in explaining the geodetic residual of non-atmospheric and non-oceanic polar motion excitations at the annual period, while the models give worse estimates with a larger phase shift or amplitude bias. At the semi-annual period, the GRACE estimates are also generally closer to the geodetic residual, but with some biases in phase or amplitude due mainly to some aliasing errors at near semi-annual period from geophysical models. For periods less than 1-year, the hydrological models and GRACE are generally worse in explaining the intraseasonal polar motion excitations.     

Assessment of the Global and Regional Land Hydrosphere and Its Impact on the Balance of the Geophysical Excitation Function of Polar Motion

The impact of continental hydrological loading from land water, snow and ice on polar motion excitation, calculated as hydrological angular momentum (HAM), is difficult to estimate, and not as much is known about it as about atmospheric angular momentum (AAM) and oceanic angular momentum (OAM). In this paper, regional hydrological excitations to polar motion are investigated using monthly terrestrial water storage data derived from the Gravity Recovery and Climate Experiment (GRACE) mission and from the five models of land hydrology. The results show that the areas where the variance shows large variability are similar for the different models of land hydrology and for the GRACE data. Areas which have a small amplitude on the maps make an important contribution to the global hydrological excitation function of polar motion. The comparison of geodetic residuals and global hydrological excitation functions of polar motion shows that none of the hydrological excitation has enough energy to significantly improve the agreement between the observed geodetic excitation and geophysical ones.     

Improved Constraints on Models of Glacial Isostatic Adjustment: A Review of the Contribution of Ground-Based Geodetic Observations

The provision of accurate models of Glacial Isostatic Adjustment (GIA) is presently a priority need in climate studies, largely due to the potential of the Gravity Recovery and Climate Experiment (GRACE) data to be used to determine accurate and continent-wide assessments of ice mass change and hydrology. However, modelled GIA is uncertain due to insufficient constraints on our knowledge of past glacial changes and to large simplifications in the underlying Earth models. Consequently, we show differences between models that exceed several mm/year in terms of surface displacement for the two major ice sheets: Greenland and Antarctica. Geodetic measurements of surface displacement offer the potential for new constraints to be made on GIA models, especially when they are used to improve structural features of the Earth’s interior as to allow for a more realistic reconstruction of the glaciation history. We present the distribution of presently available campaign and continuous geodetic measurements in Greenland and Antarctica and summarise surface velocities published to date, showing substantial disagreement between techniques and GIA models alike. We review the current state-of-the-art in ground-based geodesy (GPS, VLBI, DORIS, SLR) in determining accurate and precise surface velocities. In particular, we focus on known areas of need in GPS observation level models and the terrestrial reference frame in order to advance geodetic observation precision/accuracy toward 0.1 mm/year and therefore further constrain models of GIA and subsequent present-day ice mass change estimates.     

GRACE terrestrial water storage data assimilation based on the ensemble four-dimensional variational method PODEn4DVar:Method and validation

Seasonal and interannual changes in the Earth's gravity field are mainly due to mass exchange among the atmosphere,ocean,and continental water sources.The terrestrial water storage changes,detected as gravity changes by the Gravity Recovery and Climate Experiment(GRACE) satellites,are mainly caused by precipitation,evapotranspiration,river transportation and downward infiltration processes.In this study,a land data assimilation system LDAS-G was developed to assimilate the GRACE terrestrial water storage(TWS) data into the Community Land Model(CLM3.5) using the POD-based ensemble four-dimensional variational assimilation method PODEn4 DVar,disaggregating the GRACE large-scale terrestrial water storage changes vertically and in time,and placing constraints on the simulation of vertical hydrological variables to improve land surface hydrological simulations.The ideal experiments conducted at a single point and assimilation experiments carried out over China by the LDAS-G data assimilation system showed that the system developed in this study improved the simulation of land surface hydrological variables,indicating the potential of GRACE data assimilation in large-scale land surface hydrological research and applications.     

Significance of secular trends of mass variations determined from GRACE solutions

Since 2002 the Earth’s gravity field is globally observed by the Gravity Recovery and Climate Experiment (GRACE) satellite mission. The GRACE monthly gravity field solutions, available from several analysis centres, reflect mass variations in the atmosphere, hydrosphere and geosphere. Due to correlated noise contained in these solutions, it is, however, first necessary to apply an appropriate filtering technique. The resulting, smoothed time series are applied not only to determine variations with different periodic signatures (e.g., seasonal, short and medium-term), but to derive long-periodic mass variations and secular trends as well. As the GRACE monthly solutions always show the integral effect of all mass variations, for separation of single processes, like the GIA (Glacial isostatic adjustment)-related mass increase in Fennoscandia, appropriate reduction models (e.g. from hydrology) are necessary.In this study we show for the example of the Fennoscandian uplift area that GRACE solutions from different analysis centres yield considerably different secular trends. Furthermore, it turns out that the inevitable filtering of the monthly gravity field models affects not only the amplitudes of the signals, but also their spatial resolution and distribution such as the spatial form of the detected signals. It also becomes evident that the determination of trends has to be performed together with the determination of periodic components. All periodic terms which are really contained in the data, and only such, have to be included. The restricted time span of the available GRACE measurements, however, limits the separation of long-periodic and secular signals. It is shown that varying the analysis time span affects the results considerably. Finally, a reduction of hydrological signals from the detected integral secular trends using global hydrological models (WGHM, LaDWorld, GLDAS) is attempted. The differences among the trends resulting from different models illustrate that the state-of-the-art hydrology models are not suitable for this purpose as yet. Consequently, taking the GRACE monthly gravity field solutions from one centre, choosing a single filter and applying an insufficiently reliable reduction model leads sometimes to a misinterpretation of considered geophysical processes. Therefore, one has to be cautious with the final interpretation of the results.     

青藏高原GRACE卫星重力长期变化

本文采用最新的GRACE(Gravity Recovery and Climate Experiment)(RL05)数据,通过水文模型(Global Land Data Assimilation System,GLDAS与Climate Prediction Center,CPC)扣除土壤水及雪水的影响,利用Paulson提供的冰川模型结果扣除GIA(Glacial Isostatic Adjustment)的影响,采用尺度因子的方法减少数据处理过程中误差的影响,最终基于最小二乘计算方法得到2003—2013中国及周边地区长期性重力异常变化情况.结果发现青藏高原有较为明显的重力上升信号,我们认为该信号可能由印度板块俯冲欧亚板块导致青藏高原地壳增厚所引起.接着依据GPS观测结果和艾黎均衡假说构建了地壳形变模型并通过直立长方体模型予以正演模拟分析.以班公湖—怒江断裂带为界将青藏高原划分为南北两大区块,结果显示青藏高原重力异常大致以0.2μGal·a-1的速率在递增,小于GRACE得到的0.3±0.08μGal·a-1的增长速率(对应于地壳增厚速率约3mm·a-1),剩余未解释部分可能与湖水、冰川因素、冻土因素等有关.该结果对于认识青藏高原隆升动力学有一定参考意义.     

Detection of Continental Hydrology and Glaciology Signals from GRACE: A Review

Since its launch in March 2002, the Gravity Recovery and Climate Experiment (GRACE) has provided a global mapping of the time-variations of the Earth’s gravity field. Tiny variations of gravity from monthly to decadal time scales are mainly due to redistributions of water mass inside the surface fluid envelops of our planet (i.e., atmosphere, ocean and water storage on continents). In this article, we present a review of the major contributions of GRACE satellite gravimetry in global and regional hydrology. To date, many studies have focused on the ability of GRACE to detect, for the very first time, the time-variations of continental water storage (including surface waters, soil moisture, groundwater, as well as snow pack at high latitudes) at the unprecedented resolution of ~400–500 km. As no global complete network of surface hydrological observations exists, the advances of satellite gravimetry to monitor terrestrial water storage are significant and unique for determining changes in total water storage and water balance closure at regional and continental scales.     

Comparison of the geophysical excitations of polar motion from the period: 1980.0–2009.0

In this study we compared contributions to polar motion excitation determined separately from each of three kinds of geophysical data: atmospheric pressure, equivalent water height estimated from hydrological models, and harmonic coefficients of the Earth gravity field obtained from Gravity Recovery and Climate Experiment (GRACE). Hydrological excitation function (Hydrological Angular Momentum — HAM) has been estimated from models of global hydrology, based on the observed distribution of surfacewater, snow, ice, and soil moisture. In our considerationwe used several global models of land hydrosphere and models ofAtmospheric Angular Momentum (AAM) and Oceanic Angular Momentum (OAM). All of themwere compared with observed Geodetic Angular Momentum (GAM). The spectra of the following excitation functions of polar motion: GAM, AAM+OAM, AAM+OAM+HAM, GAM-AAM-OAM residual geodetic excitation function, and HAM were computed too. The time variable spectra of geodetic, gravimetric, and the sum of atmospheric, oceanic, and hydrological excitation functions are also presented. Phasor diagrams of the seasonal components of polar motion excitation functions of all HAM excitation functions as well as of two GRACE solutions: Center for Space Research (CSR), Centre National d’Etudes Spatiales/Groupe de Recherche en Geodesie Spatiale (CNES/GRGS) were determined and discussed.     

Glacial isostatic adjustment in Fennoscandia from GRACE data and comparison with geodynamical models

The Earth’s gravity field observed by the Gravity Recovery and Climate Experiment (GRACE) satellite mission shows variations due to the integral effect of mass variations in the atmosphere, hydrosphere and geosphere. Several institutions, such as the GeoForschungsZentrum (GFZ) Potsdam, the University of Texas at Austin, Center for Space Research (CSR) and the Jet Propulsion Laboratory (JPL), Pasadena, provide GRACE monthly solutions, which differ slightly due to the application of different reduction models and centre-specific processing schemes. The GRACE data are used to investigate the mass variations in Fennoscandia, an area which is strongly influenced by glacial isostatic adjustment (GIA). Hence the focus is set on the computation of secular trends. Different filters (e.g. isotropic and non-isotropic filters) are discussed for the removal of high frequency noise to permit the extraction of the GIA signal. The resulting GRACE based mass variations are compared to global hydrology models (WGHM, LaDWorld) in order to (a) separate possible hydrological signals and (b) validate the hydrology models with regard to long period and secular components. In addition, a pattern matching algorithm is applied to localise the uplift centre, and finally the GRACE signal is compared with the results from a geodynamical modelling. The GRACE data clearly show temporal gravity variations in Fennoscandia. The secular variations are in good agreement with former studies and other independent data. The uplift centre is located over the Bothnian Bay, and the whole uplift area comprises the Scandinavian Peninsula and Finland. The secular variations derived from the GFZ, CSR and JPL monthly solutions differ up to 20%, which is not statistically significant, and the largest signal of about 1.2 Gal/year is obtained from the GFZ solution. Besides the GIA signal, two peaks with positive trend values of about 0.8 Gal/year exist in central eastern Europe, which are not GIA-induced, and also not explainable by the hydrology models. This may indicate that the recent global hydrology models have to be revised with respect to long period and secular components. Finally, the GRACE uplift signal is also in quite good agreement with the results from a simple geodynamical modelling.     

Seismologic applications of GRACE time-variable gravity measurements

The Gravity Recovery and Climate Experiment (GRACE) has been measuring temporal and spatial variations of mass redistribution within the Earth system since 2002. As large earthquakes cause significant mass changes on and under the Earth’s surface, GRACE provides a new means from space to observe mass redistribution due to earthquake deformations. GRACE serves as a good complement to other earthquake measurements because of its extensive spatial coverage and being free from terrestrial restriction. During its over 10 years mission, GRACE has successfully detected seismic gravitational changes of several giant earthquakes, which include the 2004 Sumatra–Andaman earthquake, 2010 Maule (Chile) earthquake, and 2011 Tohoku-Oki (Japan) earthquake. In this review, we describe by examples how to process GRACE time-variable gravity data to retrieve seismic signals, and summarize the results of recent studies that apply GRACE observations to detect co- and post-seismic signals and constrain fault slip models and viscous lithospheric structures. We also discuss major problems and give an outlook in this field of GRACE application.     

Use of GRACE determined secular gravity rates for glacial isostatic adjustment studies in North-America

Monthly geopotential spherical harmonic coefficients from the GRACE satellite mission are used to determine their usefulness and limitations for studying glacial isostatic adjustment (GIA) in North-America. Secular gravity rates are estimated by unweighted least-squares estimation using release 4 coefficients from August 2002 to August 2007 provided by the Center for Space Research (CSR), University of Texas. Smoothing is required to suppress short wavelength noise, in addition to filtering to diminish geographically correlated errors, as shown in previous studies. Optimal cut-off degrees and orders are determined for the destriping filter to maximize the signal to noise ratio. The halfwidth of the Gaussian filter is shown to significantly affect the sensitivity of the GRACE data (with respect to upper mantle viscosity and ice loading history). Therefore, the halfwidth should be selected based on the desired sensitivity.It is shown that increase in water storage in an area south west of Hudson Bay, from the summer of 2003 to the summer of 2006, contributes up to half of the maximum estimated gravity rate. Hydrology models differ in the predictions of the secular change in water storage, therefore even 4-year trend estimates are influenced by the uncertainty in water storage changes. Land ice melting in Greenland and Alaska has a non-negligible contribution, up to one-fourth of the maximum gravity rate.The estimated secular gravity rate shows two distinct peaks that can possibly be due to two domes in the former Pleistocene ice cover: west and south east of Hudson Bay. With a limited number of models, a better fit is obtained with models that use the ICE-3G model compared to the ICE-5G model. However, the uncertainty in interannual variations in hydrology models is too large to constrain the ice loading history with the current data span. For future work in which GRACE will be used to constrain ice loading history and the Earth's radial viscosity profile, it is important to include realistic uncertainty estimates for hydrology models and land ice melting in addition to the effects of lateral heterogeneity.     

Improvement of Global Hydrological Models Using GRACE Data

After about 6 years of GRACE (Gravity Recovery and Climate Experiment) satellite mission operation, an unprecedented global data set on the spatio-temporal variations of the Earth’s water storage is available. The data allow for a better understanding of the water cycle at the global scale and for large river basins. This review summarizes the experiences that have been made when comparing GRACE data with simulation results of global hydrological models and it points out the prerequisites and perspectives for model improvements by combination with GRACE data. When evaluated qualitatively at the global scale, water storage variations on the continents from GRACE agreed reasonably well with model predictions in terms of their general seasonal dynamics and continental-scale spatial patterns. Differences in amplitudes and phases of water storage dynamics revealed in more detailed analyses were mainly attributed to deficiencies in the meteorological model forcing data, to missing water storage compartments in the model, but also to limitations and errors of the GRACE data. Studies that transformed previously identified model deficiencies into adequate modifications of the model structure or parameters are still rare. Prerequisites for a comprehensive improvement of large-scale hydrological models are in particular the consistency of GRACE observation and model variables in terms of filtering, reliable error estimates, and a full assessment of the water balance. Using improvements in GRACE processing techniques, complementary observation data, multi-model evaluations and advanced methods of multi-objective calibration and data assimilation, considerable progress in large-scale hydrological modelling by integration of GRACE data can be expected.     

联合GRACE和ICESat数据分离南极冰川均衡调整(GIA)信号

2002年发射的GRACE重力卫星为南极冰盖质量平衡提供了一种新的测量方式,但由于南极GIA模型的不确定较大,进而影响GRACE结果的可靠性.本文联合2003—2009年的GRACE和ICESat等数据实现了南极GIA信号的分离,联合方法所分离的GIA不依赖于不确定性很大的冰负荷等假设模型,而是直接基于卫星观测数据估算而来的,具有更大的可靠性.在分离过程中,本文提出了冰流速度加权改正法和GPS球谐拟合改正法对GIA结果进行精化,同时引入了南极GPS观测站的位移数据对分离的GIA进行详细的评估和验证,GPS验证表明经过冰流速度加权和GPS球谐拟合双改正后的GIA结果精度明显得到提高.最后本文利用所分离的GIA对GRACE和ICESat结果进行了改正,得到2003—2009年南极冰盖质量变化的趋势为-66.7±54.5 Gt/a(GRACE)和-77.2±21.5Gt/a(ICESat),相比采用其他的GIA模型,本文的GIA结果使GRACE和ICESat这两种不同观测技术得到的南极冰盖质量变化结果更加趋于一致.     

Using postglacial sea level,crustal velocities and gravity-rate-of-change to constrain the influence of thermal effects on mantle lateral heterogeneities

Lateral heterogeneities in the mantle can be caused by thermal, chemical and non-isotropic pre-stress effects. Here, we investigate the possibility of using observations of the glacial isostatic adjustment (GIA) process to constrain the thermal contribution to lateral variations in mantle viscosity. In particular, global historic relative sea level, GPS in Laurentide and Fennoscandia, altimetry together with tide-gauge data in the Great Lakes area, and GRACE data in Laurentide are used. The lateral viscosity perturbations are inferred from the seismic tomography model S20A by inserting the scaling factor β to determine the contribution of thermal effects versus compositional heterogeneity and non-isotropic pre-stress effects on lateral heterogeneity in mantle viscosity. When β = 1, lateral velocity variations are caused by thermal effects alone. With β < 1, the contribution of thermal effect decreases, so that for β = 0, there is no lateral viscosity variation and the Earth is laterally homogeneous. These lateral viscosity variations are superposed on four different reference models which differ significantly in the lower mantle viscosity. The Coupled Laplace Finite Element method is used to predict the GIA response on a spherical, self-gravitating, compressible, viscoelastic Earth with self-gravitating oceans, induced by the ICE-4G deglaciation model.Results show that the effect of β on uplift rates and gravity rate-of-change is not simple and involves the trade-off between the contribution of lateral viscosity variations in the transition zone and in the lower mantle. Models with small viscosity contrast in the lower mantle cannot explain the observed uplift rates in Laurentide and Fennoscandia. However, the RF3S20 model with a reference viscosity profile simplified from Peltier's VM2 with the value of β around 0.2–0.4 is found to explain most of the global RSL data, the uplift rates in Laurentide and Fennoscandia and the BIFROST horizontal velocity data. In addition, the changes in GIA signals caused by changes in the value of β are large enough to be detected by the data, although uncertainty in other parameters in the GIA models still exists. This may encourage us to further utilize GIA observations to constrain the thermal effect on mantle lateral heterogeneity as geodetic and satellite gravity measurements are improved.     

Using postglacial sea level, crustal velocities and gravity-rate-of-change to constrain the influence of thermal effects on mantle lateral heterogeneities

Lateral heterogeneities in the mantle can be caused by thermal, chemical and non-isotropic pre-stress effects. Here, we investigate the possibility of using observations of the glacial isostatic adjustment (GIA) process to constrain the thermal contribution to lateral variations in mantle viscosity. In particular, global historic relative sea level, GPS in Laurentide and Fennoscandia, altimetry together with tide-gauge data in the Great Lakes area, and GRACE data in Laurentide are used. The lateral viscosity perturbations are inferred from the seismic tomography model S20A by inserting the scaling factor β to determine the contribution of thermal effects versus compositional heterogeneity and non-isotropic pre-stress effects on lateral heterogeneity in mantle viscosity. When β = 1, lateral velocity variations are caused by thermal effects alone. With β < 1, the contribution of thermal effect decreases, so that for β = 0, there is no lateral viscosity variation and the Earth is laterally homogeneous. These lateral viscosity variations are superposed on four different reference models which differ significantly in the lower mantle viscosity. The Coupled Laplace Finite Element method is used to predict the GIA response on a spherical, self-gravitating, compressible, viscoelastic Earth with self-gravitating oceans, induced by the ICE-4G deglaciation model.Results show that the effect of β on uplift rates and gravity rate-of-change is not simple and involves the trade-off between the contribution of lateral viscosity variations in the transition zone and in the lower mantle. Models with small viscosity contrast in the lower mantle cannot explain the observed uplift rates in Laurentide and Fennoscandia. However, the RF3S20 model with a reference viscosity profile simplified from Peltier's VM2 with the value of β around 0.2–0.4 is found to explain most of the global RSL data, the uplift rates in Laurentide and Fennoscandia and the BIFROST horizontal velocity data. In addition, the changes in GIA signals caused by changes in the value of β are large enough to be detected by the data, although uncertainty in other parameters in the GIA models still exists. This may encourage us to further utilize GIA observations to constrain the thermal effect on mantle lateral heterogeneity as geodetic and satellite gravity measurements are improved.     

冰川均衡调整对东亚重力和海平面变化的影响

新的全球冰川均衡调整(GIA)模型RF3L20(β=0.4)+ICE-4G考虑了地幔黏滞度沿横向的变化,其黏滞度参数得到大地测量、历史相对海平面变化观测和地震剪切波层析模型的较好约束.本文利用该模型预测了东亚现今重力变化和海平面变化,根据当前末次冰川时空变化和黏滞度参考模型中下地幔下部黏滞度认识的差异,评估了预测的不确定性.结果表明,GIA对东亚地区重力场和海平面长期变化有显著的影响:例如,在哈尔滨、长春、泰安、蓟县、郑州、武汉等测站,GIA重力影响达几十纳伽,可用超导重力仪和未来原子重力仪观测出来;在东亚大陆GIA对GRACE监测的等效水柱长期变化的影响为3%~10%,其中青藏高原西部、华北和三峡地区的影响较大.在东海-太平洋区,GIA的相对影响高达20%~40%;GIA使东亚海域绝对海平面以0.27~0.37 mm/a的速率在长期下降,在黄海、东海卫星测高监测的绝对海平面长期变化中,GIA的相对影响分别达6.9%和7.5%;在58个验潮站,平均相对海平面长期上升速率为2.22 mm/a,GIA影响为-0.17 mm/a,其中14个测站GIA的影响达-0.3~-0.4 mm/a.本文GIA预测的结果,对在东亚地区发现弱的地球动力学过程信号、监测水质量长期变化、监测海平面长期变化和分析其机制,提供精密的改正模型.     

Integration of GRACE mass variations into a global hydrological model

Time-variable gravity data of the GRACE (Gravity Recovery And Climate Experiment) satellite mission provide global information on temporal variations of continental water storage. In this study, we incorporate GRACE data for the first time directly into the tuning process of a global hydrological model to improve simulations of the continental water cycle. For the WaterGAP Global Hydrology Model (WGHM), we adopt a multi-objective calibration framework to constrain model predictions by both measured river discharge and water storage variations from GRACE and illustrate it on the example of three large river basins: Amazon, Mississippi and Congo. The approach leads to improved simulation results with regard to both objectives. In case of monthly total water storage variations we obtained a RMSE reduction of about 25 mm for the Amazon, 6 mm for the Mississippi and 1 mm for the Congo river basin. The results highlight the valuable nature of GRACE data when merged into large-scale hydrological modeling. Furthermore, they reveal the utility of the multi-objective calibration framework for the integration of remote sensing data into hydrological models.     

Correction of Hydrological and Oceanic Effects from GRACE Data by Combination of the Steric Sea Level,Altimetry Data and GLDAS Model

In this study, a scheme to estimate oceanic and hydrological effects in the GRACE (Gravity Recovery and Climate Experiment) data is presented. The aim is to reveal tectonic signals for the case of the Sumatra earthquake on 26 December 2004. The variations of hydrological and oceanic effects are estimated with the aid of data set of GRACE, altimetry, World Ocean Atlas, and the GLDAS model for a period of January 2003 to December 2006. The time series of computed gravity changes over Sumatra region show some correlations to the deformation resulting from the earthquake occurred in December 2004. The maximum and minimum impacts of hydrological and oceanic effects on gravity changes are about 3 μGal in radial direction and–5 μGal in northward direction. The maximum and minimum amounts of gravitational gradient changes after the correction are 0.2 and–0.25 mE, which indicates the significant influences of hydrological and oceanic sources on the desired signal.