![]() Meyrath, Thierry ![]() ![]() in Geophysical Journal International (2017), 208(2), 1126-1138 We study fluctuations in the degree-2 zonal spherical harmonic coefficient of the Earth's gravity potential, $C_{20}$, over the period 2003-2015. This coefficient is related to the Earth's oblateness and ... [more ▼] We study fluctuations in the degree-2 zonal spherical harmonic coefficient of the Earth's gravity potential, $C_{20}$, over the period 2003-2015. This coefficient is related to the Earth's oblateness and studying its temporal variations, $\Delta C_{20}$, can be used to monitor large-scale mass movements between high and low latitude regions. We examine $\Delta C_{20}$ inferred from six different sources, including satellite laser ranging (SLR), GRACE and global geophysical fluids models. We further include estimates that we derive from measured variations in the length-of-day (LOD), from the inversion of global crustal displacements as measured by GPS, as well as from the combination of GRACE and the output of an ocean model as described by \cite{sunetal2016}. We apply a sequence of trend- and seasonal moving average filters to the different time series in order to decompose them into an interannual, a seasonal and an intraseasonal component. We then perform a comparison analysis for each component, and we further estimate the noise level contained in the different series using an extended version of the three-cornered-hat method. For the seasonal component, we generally obtain a very good agreement between the different sources, and except for the LOD-derived series, we find that over 90\% of the variance in the seasonal components can be explained by the sum of an annual and semiannual oscillation of constant amplitudes and phases, indicating that the seasonal pattern is stable over the considered time period. High consistency between the different estimates is also observed for the intraseasonal component, except for the solution from GRACE, which is known to be affected by a strong tide-like alias with a period of about 161 days. Estimated interannual components from the different sources are generally in agreement with each other, although estimates from GRACE and LOD present some discrepancies. Slight deviations are further observed for the estimate from the geophysical models, likely to be related to the omission of polar ice and groundwater changes in the model combination we use. On the other hand, these processes do not seem to play an important role at seasonal and shorter time scales, as the sum of modelled atmospheric, oceanic and hydrological effects effectively explains the observed $C_{20}$ variations at those scales. We generally obtain very good results for the solution from SLR, and we confirm that this well-established technique accurately tracks changes in $C_{20}$. Good agreement is further observed for the estimate from the GPS inversion, showing that this indirect method is successful in capturing fluctuations in $C_{20}$ on scales ranging from intra- to interannual. Obtaining accurate estimates from LOD, however, remains a challenging task and more reliable models of atmospheric wind fields are needed in order to obtain high-quality $\Delta C_{20}$, in particular at the seasonal scale. The combination of GRACE data and the output of an ocean model appears to be a promising approach, particularly since corresponding $\Delta C_{20}$ is not affected by tide-like aliases, and generally gives better results than the solution from GRACE, which still seems to be of rather poor quality. [less ▲] Detailed reference viewed: 324 (23 UL)![]() Meyrath, Thierry ![]() ![]() in Journal of Geodesy (2017), 91(11), 1329-1350 Large-scale mass redistribution in the terrestrial water storage (TWS) leads to changes in the low-degree spherical harmonic coefficients of the Earth's surface mass density field. Studying these low ... [more ▼] Large-scale mass redistribution in the terrestrial water storage (TWS) leads to changes in the low-degree spherical harmonic coefficients of the Earth's surface mass density field. Studying these low-degree fluctuations is an important task that contributes to our understanding of continental hydrology. In this study, we use global GNSS measurements of vertical and horizontal crustal displacements that we correct for atmospheric and oceanic effects, and use a set of modified basis functions similar to Clarke et al. (2007) to perform an inversion of the corrected measurements in order to recover changes in the coefficients of degree-0 (hydrological mass change), degree-1 (center of mass shift) and degree-2 (flattening of the Earth) caused by variations in the TWS over the period January 2003 - January 2015. We infer from the GNSS-derived degree-0 estimate an annual variation in total continental water mass with an amplitude of $(3.49 \pm 0.19) \times 10^{3}$ Gt and a phase of $70 \pm 3^{\circ}$ (implying a peak in early March), in excellent agreement with corresponding values derived from the Global Land Data Assimilation System (GLDAS) water storage model that amount to $(3.39 \pm 0.10) \times 10^{3}$ Gt and $71 \pm 2^{\circ}$, respectively. The degree-1 coefficients we recover from GNSS predict annual geocentre motion (i.e. the offset change between the center of common mass and the center of figure) caused by changes in TWS with amplitudes of $0.69 \pm 0.07$ mm for GX, $1.31 \pm 0.08$ mm for GY and $2.60 \pm 0.13$ mm for GZ. These values agree with GLDAS and estimates obtained from the combination of GRACE and the output of an ocean model using the approach of Swenson et al. (2008) at the level of about 0.5, 0.3 and 0.9 mm for GX, GY and GZ, respectively. Corresponding degree-1 coefficients from SLR, however, generally show higher variability and predict larger amplitudes for GX and GZ. The results we obtain for the degree-2 coefficients from GNSS are slightly mixed, and the level of agreement with the other sources heavily depends on the individual coefficient being investigated. The best agreement is observed for $T_{20}^C$ and $T_{22}^S$, which contain the most prominent annual signals among the degree-2 coefficients, with amplitudes amounting to $(5.47 \pm 0.44) \times 10^{-3}$ and $(4.52 \pm 0.31) \times 10^{-3}$ m of equivalent water height (EWH), respectively, as inferred from GNSS. Corresponding agreement with values from SLR and GRACE is at the level of or better than $0.4 \times 10^{-3}$ and $0.9 \times 10^{-3}$ m of EWH for $T_{20}^C$ and $T_{22}^S$, respectively, while for both coefficients, GLDAS predicts smaller amplitudes. Somewhat lower agreement is obtained for the order-1 coefficients, $T_{21}^C$ and $T_{21}^S$, while our GNSS inversion seems unable to reliably recover $T_{22}^C$. For all the coefficients we consider, the GNSS-derived estimates from the modified inversion approach are more consistent with the solutions from the other sources than corresponding estimates obtained from an unconstrained standard inversion. [less ▲] Detailed reference viewed: 247 (15 UL)![]() Li, Zhao ![]() ![]() in Willis, Pascal (Ed.) Proceedings of the 2013 IAG Scientific Assembly, Potsdam, Germany, 1-6 September, 2013 (2015) To remove continental water storage (CWS) signals from the GPS data, CWS mass models are needed to obtain predicted surface displacements. We compared weekly GPS height time series with five CWS models ... [more ▼] To remove continental water storage (CWS) signals from the GPS data, CWS mass models are needed to obtain predicted surface displacements. We compared weekly GPS height time series with five CWS models: (1) the monthly and (2) three-hourly Global Land Data Assimilation System (GLDAS); (3) the monthly and (4) one-hourly Modern- Era Retrospective Analysis for Research and Applications (MERRA); (5) the six-hourly National Centers for Environmental Prediction-Department of Energy (NCEP-DOE) global reanalysis products (NCEP-R-2). We find that of the 344 selected global IGS stations, more than 77% of stations have their weighted root mean square (WRMS) reduced in the weekly GPS height by using both the GLDAS and MERRA CWS products to model the surface displacement, and the best improvement concentrate mainly in North America and Eurasia.We find that the one-hourly MERRA-Land dataset is the most appropriate product for modeling weekly vertical surface displacement caused by CWS variations. The threehourly GLDAS data ranks the second, while the GLDAS and MERRA monthly products rank the third. The higher spatial resolution MERRA product improves the performance of the CWS model in reducing the scatter of the GPS height by about 2–6% compared with the GLDAS. Under the same spatial resolution, the higher temporal resolution could also improve the performance by almost the same magnitude. We also confirm that removing the ATML and NTOL effects from the weekly GPS height would remarkably improve the performance of CWS model in correcting the GPS height by at least 10%, especially for coastal and island stations. Since the GLDAS product has a much greater latency than the MERRA product, MERRA would be a better choice to model surface displacements from CWS. Finally, we find that the NCEP-R-2 data is not sufficiently precise to be used for this application. Further work is still required to determine the reason. [less ▲] Detailed reference viewed: 151 (11 UL)![]() ; van Dam, Tonie ![]() in Journal of Geodynamics (2013), 72 Seasonal signals in GPS time series are of great importance for understanding the evolution of regional mass fluctuations, i.e., ice, hydrology, and ocean mass. Conventionally these signals quasi-annual ... [more ▼] Seasonal signals in GPS time series are of great importance for understanding the evolution of regional mass fluctuations, i.e., ice, hydrology, and ocean mass. Conventionally these signals quasi-annual and semi-annual signals are modeled by least-squares fitting harmonic terms with a constant amplitude and phase. In reality, however, such seasonal signals are modulated, i.e., they will have a time-variable amplitude and phase. Recently, Davis et al.(2012) proposed a Kalman filter based approach to capture the stochastic seasonal behavior of geodetic time series. Singular Spectrum Analysis (SSA) is a non-parametric method, which uses time domain data to extract information from short and noisy time series without a priori knowledge of the dynamics affecting the time series. A prominent benefit is that trends obtained in this way are not necessarily linear. Further, true oscillations can be amplitude and phase modulated. In this work, we will assess the value of SSA for extracting time-variable seasonal signals from GPS time series. We compare our SSA-based results to those obtained using 1) least-squares analysis and 2) Kalman filtering. Our results demonstrate that SSA is a viable and complementary tool for extracting modulated oscillations from GPS time series. [less ▲] Detailed reference viewed: 548 (28 UL) |
||