Simulation of the time-variable gravity field by means of coupled geophysical models

Abstract. Time variable gravity fields, reflecting variations of mass distribution in the system Earth is one of the key parameters to understand the changing Earth. Mass variations are caused either by redistribution of mass in, on or above the Earth's surface or by geophysical processes in the Earth's interior. The first set of observations of monthly variations of the Earth gravity field was provided by the US/German GRACE satellite mission beginning in 2002. This mission is still providing valuable information to the science community. However, as GRACE has outlived its expected lifetime, the geoscience community is currently seeking successor missions in order to maintain the long time series of climate change that was begun by GRACE. Several studies on science requirements and technical feasibility have been conducted in the recent years. These studies required a realistic model of the time variable gravity field in order to perform simulation studies on sensitivity of satellites and their instrumentation. This was the primary reason for the European Space Agency (ESA) to initiate a study on ''Monitoring and Modelling individual Sources of Mass Distribution and Transport in the Earth System by Means of Satellites''. The goal of this interdisciplinary study was to create as realistic as possible simulated time variable gravity fields based on coupled geophysical models, which could be used in the simulation processes in a controlled environment. For this purpose global atmosphere, ocean, continental hydrology and ice models were used. The coupling was performed by using consistent forcing throughout the models and by including water flow between the different domains of the Earth system. In addition gravity field changes due to solid Earth processes like continuous glacial isostatic adjustment (GIA) and a sudden earthquake with co-seismic and post-seismic signals were modelled. All individual model results were combined and converted to gravity field spherical harmonic series, which is the quantity commonly used to describe the Earth's global gravity field. The result of this study is a twelve-year time-series of 6-hourly time variable gravity field spherical harmonics up to degree and order 180 corresponding to a global spatial resolution of 1 degree in latitude and longitude. In this paper, we outline the input data sets and the process of combining these data sets into a coherent model of temporal gravity field changes. The resulting time series was used in some follow-on studies and is available to anybody interested.

[1]  W. Farrell.  Deformation of the Earth by surface loads , 1972 .

[2]  Jakob Flury,et al.  Future Satellite Gravimetry and Earth Dynamics , 2006 .

[3]  H. Dobslaw,et al.  Atmospheric induced oceanic tides from ECMWF forecasts , 2005 .

[4]  M. Tamisiea,et al.  GRACE Gravity Data Constrain Ancient Ice Geometries and Continental Dynamics over Laurentia , 2007, Science.

[5]  R. Greatbatch.  A note on the representation of steric sea level in models that conserve volume rather than mass , 1994 .

[6]  E. Maier‐Reimer,et al.  The Hamburg Ocean primitive equation model - HOPE , 1996 .

[7]  S. Levitus.  Climatological Atlas of the World Ocean , 1982 .

[8]  Mohamed Ali Mekouar,et al.  15. United Nations Food and Agriculture Organization (FAO) , 2011 .

[9]  W. Hibler.  A Dynamic Thermodynamic Sea Ice Model , 1979 .

[10]  E. Rignot,et al.  Changes in the Velocity Structure of the Greenland Ice Sheet , 2006, Science.

[11]  E. Maier‐Reimer,et al.  Mean Circulation and Internal Variability in an Ocean Primitive Equation Model , 1996 .

[12]  Henryk Dobslaw.  Modellierung der allgemeinen ozeanischen Dynamik zur Korrektur und Interpretation von Satellitendaten. , 2008 .

[13]  Chen Ji,et al.  Crustal Dilatation Observed by GRACE After the 2004 Sumatra-Andaman Earthquake , 2006, Science.

[14]  James L. Davis,et al.  Geodetic Constraints on Glacial Isostatic Adjustment , 2013 .

[15]  J. Kusche.  Approximate decorrelation and non-isotropic smoothing of time-variable GRACE-type gravity field models , 2007 .

[16]  J. Thepaut,et al.  Assimilation and Modeling of the Atmospheric Hydrological Cycle in the ECMWF Forecasting System , 2005 .

[17]  Sol Hellerman,et al.  Normal Monthly Wind Stress Over the World Ocean with Error Estimates , 1983 .

[18]  Bradford H. Hager,et al.  Localization of the gravity field and the signature of glacial rebound , 1997, Nature.

[19]  D. L. Anderson,et al.  Preliminary reference earth model , 1981 .

[20]  E. van Meijgaard,et al.  Reassessment of the Antarctic surface mass balance using calibrated output of a regional atmospheric climate model , 2006 .

[21]  Maik Thomas,et al.  Simulation and observation of global ocean mass anomalies , 2007 .

[22]  Jonathan L. Bamber,et al.  A surface mass balance model for the Greenland Ice Sheet , 2005 .

[23]  T. Gruber,et al.  Simulation of oceanic bottom pressure for gravity space missions , 2001 .

[24]  R. Sabadini.  Applications of Normal Mode Relaxation Theory to Solid Earth Geophysics , 2004 .

[25]  L. S. Pereira,et al.  Revised FAO Procedures for Calculating Evapotranspiration: Irrigation and Drainage Paper No. 56 with Testing in Idaho , 2001 .

[26]  W. Peltier.  GLOBAL GLACIAL ISOSTASY AND THE SURFACE OF THE ICE-AGE EARTH: The ICE-5G (VM2) Model and GRACE , 2004 .

[27]  M. Watkins,et al.  GRACE Measurements of Mass Variability in the Earth System , 2004, Science.

[28]  Rongjiang Wang,et al.  Gravity changes due to the Sumatra-Andaman and Nias earthquakes as detected by the GRACE satellites: a reexamination , 2010 .

[29]  Jürgen Sündermann,et al.  Consideration of ocean tides in an OGCM and impacts on subseasonal to decadal polar motion excitation , 2001 .

[30]  L. Vermeersen,et al.  High-harmonic geoid signatures related to glacial isostatic adjustment and their detectability by GOCE , 2006 .

[31]  Marc F. P. Bierkens,et al.  Seasonal Predictability of European Discharge: NAO and Hydrological Response Time , 2009 .

[32]  Nico Sneeuw,et al.  Space-borne gravimetric satellite constellations and ocean tides: aliasing effects , 2010 .

[33]  L. Vermeersen,et al.  Sensitivity of glacial isostatic adjustment models with shallow low-viscosity earth layers to the ice-load history in relation to the performance of GOCE and GRACE , 2005 .

[34]  Eric Rignot,et al.  Recent Antarctic ice mass loss from radar interferometry and regional climate modelling , 2008 .

[35]  I. Velicogna.  Increasing rates of ice mass loss from the Greenland and Antarctic ice sheets revealed by GRACE , 2009 .

[36]  R. Sabadini,et al.  Influence of lithospheric and mantle stratification on global post‐seismic deformation , 1997 .

[37]  Bob E. Schutz,et al.  Glacial Isostatic Adjustment over Antarctica from combined ICESat and GRACE satellite data , 2009 .

[38]  M. R. van den Broeke,et al.  Partitioning Recent Greenland Mass Loss , 2009, Science.

[39]  B. Vermeersen,et al.  Challenges From Solid Earth Dynamics for Satellite Gravity Field Missions in the Post-Goce Era , 2004 .

[40]  H. Dobslaw,et al.  Impact of river run-off on global ocean mass redistribution , 2007 .