Hydrological Angular Momentum Research Articles

Analysis of combined series of hydrological angular momentum developed based on climate models

Mass changes in the hydrosphere represent an important contributor to polar motion (PM) variations, especially at seasonal time scales (i.e., annual and semiannual). Although well studied, hydrological angular momentum (HAM) remains a major source of uncertainty in estimating PM excitation. In this work, we use a large number of climate models from the sixth phase of the Coupled Model Intercomparison Project (CMIP6) to determine HAM series both from individual models and their combination, formed with a multi-model mean, a weighted mean, and a three-cornered hat (TCH) method. The CMIP6-based HAM series are analysed in several spectral bands and evaluated with a reference hydrological signal in geodetically observed PM excitation (GAO). HAM determined from CMIP6 were also compared to HAM calculated from Gravity Recovery and Climate Experiment (GRACE) measurements. We find that while climate models do not allow for reliable estimation of non-seasonal changes in HAM, they can help interpret seasonal variability. For annual prograde and semiannual retrograde oscillations, several combined CMIP6-based series exhibit higher amplitude and phase consistency with GAO than the corresponding series computed from GRACE data. Whether one uses a simple average of the models, a weighted average, or a combination of models from the TCH method has little impact on the resulting HAM series and their level of agreement with GAO. Our study advances the understanding of hydrological signal in Earth’s rotation at seasonal time scales.

Analysis of The Itsg-Grace Daily Models in The Determination of Polar Motion Excitation Function

ABSTRACT The main aim of this study is to evaluate the usefulness of Institute of Geodesy at Graz University of Technology (ITSG) daily gravity field models in the determination of hydrological angular momentum (HAM) at nonseasonal time scales. We compared the equatorial components (χ1 and χ2) of HAM calculated with the ITSG daily gravity field models (ITSG-Gravity Recovery and Climate Experiment [ITSG-GRACE] 2016 and ITSG-GRACE 2018) with HAM and sea-level angular momentum (SLAM) from hydrological land surface discharge model (LSDM) and the hydrological signal in the polar motion excitation (known as geodetic residuals [GAO]). Data from ITSG have a daily temporal resolution and allow us to determine oscillations with higher frequencies than the more commonly used monthly data. We limited our study to the period between 2004 and 2011 because of the gaps in GRACE observations before and after this period. We evaluated HAM obtained from ITSG GRACE models in spectral and time domains and determined the amplitude spectra of the analyzed series in the spectral range from 2 to 120 days. Our analyses confirm the existence of a sub-monthly signal in the HAM series determined from ITSG daily data. We observed a similar signal in LSDM-based HAM, but with notably weaker amplitudes. We also observed common peaks around 14 days in the amplitude spectra for the GAO- and ITSG-based series, which may be related to the Earth’s tides. ITSG daily gravity field models can be useful to determine the equatorial components of HAM at nonseasonal time scales.

Geodetic residual time series: A combined series by minimization of their internal noise level

This study aims to assess the hydrological effects of polar motion calculated from different combinations of geophysical excitations at decadal, seasonal, and non-seasonal periods. The geodetic residuals GAO, being a difference between observed geodetic excitation function of polar motion Geodetic Angular Momentum (GAM) and atmospheric (Atmospheric Angular Momentum—AAM) plus oceanic excitation functions (Oceanic Angular Momentum—OAM), are compared. Estimating hydrological effects on Earth’s rotation differs significantly, especially when using various oceanic models. Up to now, studies of geophysical excitations of polar motion containing AAM, OAM, and hydrological angular momentum (HAM) have not achieved entire agreement between geophysical (sum of AAM, OAM, and HAM obtained from the models) and geodetic excitation. Many geophysical models of the atmosphere, oceans, and land hydrology can be used to compute polar motion excitation. However, these models are very complex and still have uncertainties in the process descriptions, parametrization, and forcing. This work aims to show differences between various GAO solutions calculated from different mass and motion terms of various AAM and OAM models. Justifying to use one combination of GAO to estimate geodetic residuals is comparing those time series to hydrological excitations computed from Gravity Recovery and Climate Experiment (GRACE) data and the Land Surface Discharge Model (LSDM) model. Especially the quality of each geodetic residual time series is determined by estimating their noise level using a generalized formulation of the “three-cornered hat method” (3CH). This study presents a combined series of geodetic residuals GAO in polar motion (PM), wherein the internal noise level is shortened to a minimum by using the 3CH method. The combined GAO time series are compared with results obtained from the GRACE/GRACE Follow-On (GRACE-FO) solution provided by International Combination Service for Time variable Gravity Fields (COST-G) and the single solution elaborated by the Center for Space Research (CSR) and from the HAM LSDM hydrological model. The results show that higher consistency between GAO and HAM excitations can be achieved by minimizing the internal noise level in the GAO combined excitation series using the 3CH method, especially for the overall broadband and seasonal oscillations. For seasonal spectral bands, an agreement between combined GAO and the best-correlated series of GRACE CSR achieve correlations as high as 0.97 and 0.83 for the χ1 and χ2 equatorial components of PM excitation, respectively. This study’s combined geodetic residual time series slightly improved consistency between observed geodetic polar motion excitations and geophysical ones.

Preliminary study on hydrological angular momentum determined from CMIP6 historical simulations

Polar motion (PM) is an essential parameter needed to transform coordinates between celestial and terrestrial reference frames, thus playing a crucial role in precise positioning and navigation. The role of hydrological signals in PM excitation is not yet fully understood, which is largely because of the lack of agreement between estimates of hydrological angular momentum (HAM) computed from different data sources. In this study, we used data obtained from the latest, sixth phase of the Coupled Model Intercomparison Project (CMIP6) to assess the impact of the continental hydrosphere on PM excitation. To do so, we exploited soil moisture and snow water variables obtained from historical simulations of CMIP6 to estimate climate-based HAM. The HAM series were computed, then we analysed their variability in terms of trends, seasonal and non-seasonal oscillations. An important part of this study is the validation of HAM estimates based on comparison with the hydrological signal in geodetically observed PM excitation (geodetic residuals, GAO). In addition, HAM series based on climate models were compared with those determined from global gravimetric data provided by the Gravity Recovery and Climate Experiment (GRACE) mission, and from the Land Surface Discharge Model (LSDM). This study also aimed to identify the most appropriate CMIP6 models for interpretation of PM variations. Overall, the correspondence between GAO and HAM received from CMIP6 was lower than the previously obtained consistency with GRACE results, and the level of agreement was dependent on the oscillation considered and the model used. However, it may be possible to identify several CMIP6 models from among the almost 100 available that provides a HAM series more compatible with GAO than HAM from GRACE or LSDM, especially in annual oscillations. The GISS-E2-1-G_historical_r10i1p1f1 model was found to provide the highest consistency with GAO for annual prograde amplitudes, GFDL-CM4_historical_r1i1p1f1 for annual retrograde amplitudes, BCC-ESM1_historical_r3i1p1f1 for the annual prograde phase, and MIROC-ES2L_historical_r2i1p1f2 for the annual retrograde phase. Because of their length, the CMIP6 data allow for analysis of the past and future changes in HAM from 1850 to 2100, which is of particular importance in the exploration of the impact of climate change on PM excitation.Graphical

Prograde and Retrograde Terms of Gravimetric Polar Motion Excitation Estimates from the GRACE Monthly Gravity Field Models

From 2002 to 2017, the Gravity Recovery and Climate Experiment (GRACE) mission’s twin satellites measured variations in the mass redistribution of Earth’s superficial fluids, which disturb polar motion (PM). In this study, the PM excitation estimates were computed from two recent releases of GRACE monthly gravity field models, RL05 and RL06, and converted into prograde and retrograde circular terms by applying the complex Fourier transform. This is the first such analysis of circular parts in GRACE-based excitations. The obtained series were validated by comparison with the residuals of observed polar motion excitation (geodetic angular momentum (GAM)–atmospheric angular momentum (AAM)–oceanic angular momentum (OAM) (GAO)) determined from precise geodetic measurements of the pole coordinates. We examined temporal variations of hydrological excitation function series (or hydrological angular momentum, HAM) in four spectral bands: seasonal, non-seasonal, non-seasonal short-term, and non-seasonal long-term. The general conclusions arising from the conducted analyses of prograde and retrograde terms were consistent with the findings from the equatorial components of PM excitation studies drawn in previous research. In particular, we showed that the new GRACE RL06 data increased the consistency between different solutions and improved the agreement between GRACE-based excitation series and reference data. The level of agreement between HAM and GAO was dependent on the oscillation considered and was higher for long-term than short-term variations. For most of the oscillations considered, the highest agreement with GAO was obtained for CSR RL06 and ITSG-Grace2018 solutions. This study revealed that both prograde and retrograde circular terms of PM excitation can be determined by GRACE with similar levels of accuracy. The findings from this study may help in choosing the most appropriate GRACE solution for PM investigations and can be useful in future improvements to GRACE data processing.

Determining and Evaluating the Hydrological Signal in Polar Motion Excitation from Gravity Field Models Obtained from Kinematic Orbits of LEO Satellites

This study evaluates the gravity field solutions based on high-low satellite-to-satellite tracking (hl-SST) of low-Earth-orbit (LEO) satellites: GRACE, Swarm, TerraSAR-X, TanDEM-X, MetOp-A, MetOp-B, and Jason 2, by converting them into hydrological polar motion excitation functions (or hydrological angular momentum (HAM)). The resulting HAM series are compared with the residuals of observed polar motion excitation (geodetic residuals, GAO) derived from precise geodetic measurements, and the HAM obtained from the GRACE ITSG 2018 solution. The findings indicate a large impact of orbital altitude and inclination on the accuracy of derived HAM. The HAM series obtained from Swarm data are found to be the most consistent with GAO. Visible differences are found in HAM obtained from GRACE and Swarm orbits and provided by different processing centres. The main reasons for such differences are likely to be different processing approaches and background models. The findings of this study provide important information on alternative data sets that may be used to provide continuous polar motion excitation observations, of which the Swarm solution provided by the Astronomical Institute, Czech Academy of Sciences, is the most accurate. However, further analysis is needed to determine which processing algorithms are most appropriate to obtain the best correspondence with GAO.

Terrestrial water storage variations and their effect on polar motion

The role of continental water in polar motion excitation can be illustrated by determining Hydrological Angular Momentum calculated from terrestrial water storage (TWS). In this paper we compare global and regional changes in TWS computed using Coupled Model Intercomparison Project Phase 5 climate models, Global Land Data Assimilation System (GLDAS) land hydrology models and observations from the Gravity Recovery and Climate Experiment (GRACE) satellite mission. We also compare hydrological excitation functions derived from models with those obtained from the GRACE mission and the hydrological signal in observed polar motion excitation (the so-called geodetic residuals). The results confirm that GLDAS models of seasonal and non-seasonal TWS change are more consistent with GRACE data than climate models; on the other hand, none of the considered models are fully consistent with GRACE data or geodetic residuals. In turn, GRACE observations are most consistent with the non-seasonal hydrological signal in observed excitation. A detailed study of the contribution of different TWS components to the hydrological excitation function shows that soil moisture dominates.

Hydrological excitation of polar motion by different variables from the GLDAS models

Continental hydrological loading by land water, snow and ice is a process that is important for the full understanding of the excitation of polar motion. In this study, we compute different estimations of hydrological excitation functions of polar motion (as hydrological angular momentum, HAM) using various variables from the Global Land Data Assimilation System (GLDAS) models of the land-based hydrosphere. The main aim of this study is to show the influence of variables from different hydrological processes including evapotranspiration, runoff, snowmelt and soil moisture, on polar motion excitations at annual and short-term timescales. Hydrological excitation functions of polar motion are determined using selected variables of these GLDAS realizations. Furthermore, we use time-variable gravity field solutions from the Gravity Recovery and Climate Experiment (GRACE) to determine the hydrological mass effects on polar motion excitation. We first conduct an intercomparison of the maps of variations of regional hydrological excitation functions, timing and phase diagrams of different regional and global HAMs. Next, we estimate the hydrological signal in geodetically observed polar motion excitation as a residual by subtracting the contributions of atmospheric angular momentum and oceanic angular momentum. Finally, the hydrological excitations are compared with those hydrological signals determined from residuals of the observed polar motion excitation series. The results will help us understand the relative importance of polar motion excitation within the individual hydrological processes, based on hydrological modeling. This method will allow us to estimate how well the polar motion excitation budget in the seasonal and inter-annual spectral ranges can be closed.

Hydrological Excitations of Polar Motion Derived from Different Variables of Fgoals − g2 Climate Model

Abstract The hydrological contribution to decadal, inter-annual and multi-annual suppress polar motion derived from climate model as well as from GRACE (Gravity Recovery and Climate Experiment) data is discussed here for the period 2002.3-2016.0. The data set used here are Earth Orientation Parameters Combined 04 (EOP C04), Flexible Global Ocean-Atmosphere-Land System Model: Grid-point Version 2 (FGOAL-g2) and Global Land Data Assimilation System (GLDAS) climate models and GRACE CSR RL05 data for polar motion, hydrological and gravimetric excitation, respectively. Several Hydrological Angular Momentum (HAM) functions are calculated here from the selected variables: precipitation, evaporation, runoff, soil moisture, accumulated snow of the FGOALS and GLDAS climate models as well as from the global mass change fields from GRACE data provided by the International Earth Rotation and Reference System Service (IERS) Global Geophysical Fluids Center (GGFC). The contribution of different HAM excitation functions to achieve the full agreement between geodetic observations and geophysical excitation functions of polar motion is studied here.

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.

High frequency variations of Earth Rotation Parameters from GPS and GLONASS observations.

The Earth's rotation undergoes changes with the influence of geophysical factors, such as Earth's surface fluid mass redistribution of the atmosphere, ocean and hydrology. However, variations of Earth Rotation Parameters (ERP) are still not well understood, particularly the short-period variations (e.g., diurnal and semi-diurnal variations) and their causes. In this paper, the hourly time series of Earth Rotation Parameters are estimated using Global Positioning System (GPS), Global Navigation Satellite System (GLONASS), and combining GPS and GLONASS data collected from nearly 80 sites from 1 November 2012 to 10 April 2014. These new observations with combining different satellite systems can help to decorrelate orbit biases and ERP, which improve estimation of ERP. The high frequency variations of ERP are analyzed using a de-trending method. The maximum of total diurnal and semidiurnal variations are within one milli-arcseconds (mas) in Polar Motion (PM) and 0.5 milli-seconds (ms) in UT1-UTC. The semidiurnal and diurnal variations are mainly related to the ocean tides. Furthermore, the impacts of satellite orbit and time interval used to determinate ERP on the amplitudes of tidal terms are analyzed. We obtain some small terms that are not described in the ocean tide model of the IERS Conventions 2010, which may be caused by the strategies and models we used or the signal noises as well as artifacts. In addition, there are also small differences on the amplitudes between our results and IERS convention. This might be a result of other geophysical excitations, such as the high-frequency variations in atmospheric angular momentum (AAM) and hydrological angular momentum (HAM), which needs more detailed analysis with more geophysical data in the future.

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.

Hydrological and oceanic effects on polar motion from GRACE and models

Terrestrial water storage (TWS) and ocean bottom pressure (OBP) are major contributors to the observed polar motion excitations, second only to atmospheric mass movement. However, quantitative assessment of the hydrological and oceanic effects on polar motion remains unclear because of the lack of global observations. In this paper, hydrological and oceanic mass excitations to polar motion are investigated using monthly TWS and OBP derived from the Gravity Recovery and Climate Experiment (GRACE) for January 2003 until December 2008. The results from this analysis are compared with hydrological model excitations from the European Center for Medium‐Range Weather Forecasts (ECMWF) and oceanic model excitations obtained from the Jet Propulsion Laboratory (JPL) using Estimating the Circulation and Climate of the Ocean (ECCO). Results show that the GRACE‐derived OBP and TWS better explain the geodetic residual polar motion excitations for the Px component at the annual period, while the GRACE OBP and ECMWF hydrological angular momentum agree better with the geodetic residuals for the annual Py excitation. GRACE ocean and hydrology excitations better explain the geodetic residuals for the semiannual Py excitation. However, the JPL ECCO and ECMWF models better explain the intraseasonal geodetic residual of polar motion excitation in the Px and Py components. The GRACE data demonstrate much higher intraseasonal variability than either the models or the geodetic observations.

Atmospheric, oceanic and hydrological contributions to seasonal variations in length of day

The annual and semiannual residuals derived in the axial angular momentum budget of the solid Earth–atmosphere system reflect significant signals. They must be caused by further excitation sources. Since, in particular, the contribution for the wind term from the atmospheric layer between the 10 and 0.3 hPa levels to the seasonal variations in length of day (LOD) is still missing, it is necessary to extend the top level into the upper stratosphere up to 0.3 hPa. Under the conservation of the total angular momentum of the entire Earth, variations in the oceanic angular momentum (OAM) and the hydrological angular momentum (HAM) are further significant excitation sources at seasonal time scales. Focusing on other contributions to the Earth's axial angular momentum budget, the following data are used in this study: axial atmospheric angular momentum (AAM) data derived for the 10–0.3 hPa layer from 1991 to 1997 for computing the missing wind effects; axial OAM functions as generated by oceanic general circulation models (GCMs), namely for the ECHAM3 and the MICOM models, available from 1975 to 1994 and from 1992 to 1994, respectively, for computing the oceanic contributions to LOD changes, and, concerning the HAM variations, the seasonal estimates of the hydrological contribution as derived by Chao and O'Connor [(1988) Geophys J 94: 263–270]. Using vector representation, it is shown that the vectors achieve a close balance in the global axial angular momentum budget within the estimated uncertainties of the momentum quantities on seasonal time scales.

Hydrological and oceanic excitations to polar motion andlength-of-day variation

SUMMARY Water mass redistributions in the global hydrosphere, including continental water storage change and non-steric sea level change, introduce variations in the hydrological angular momentum (HAM) and the oceanic angular momentum (OAM). Under the conservation of angular momentum, HAM and OAM variations are signi¢cant excitation sources of the Earth rotational variations at a wide range of timescales. In this paper, we estimate HAM and OAM variations and their excitations to polar motion and length-of-day variation using soil moisture and snow estimates and non-steric sea level change determined by TOPEX/Poseidon satellite radar altimeter observations and a simpli¢ed steric sea level change model. The results are compared with the variations of polar motion and LOD that are not accounted for by the atmosphere. This study indicates that seasonal continental water storage change provides signi¢cant contributions to both polar motion and LOD variation, especially to polar motion X, and the non-steric sea level change is responsible for a major part of the remaining excitations at both seasonal scale and high frequencies, particularly in polar motion Y and LOD. The good correlation between OAM contributions and the remaining excitations shows that large-scale non-tidal mass variation exists in the oceans and can be detected by TOPEX/Poseidon altimeter observations.