Search

GRACE‐D accelerometer data show significant bias jumps since one month after the launch of the GRACE Follow‐On (GRACE‐FO) satellites in May 2018, making them inapplicable for correcting GRACE‐FO products for non‐gravitational accelerations. The GRACE‐FO Science Data System (SDS) compensated this issue by transplanting the GRACE‐C accelerometer data toward that of GRACE‐D. Recently, GRACE‐FO SDS implemented an updated transplant method, used in the latest release of GRACE‐FO data. Here, we examine the impact of updated accelerometer transplant data (ACH) on GRACE‐FO measurements at all levels: (a) Level‐1B inter‐satellite laser ranging residuals measured along satellite orbit, (b) Level‐2 time‐variable gravity solutions from all SDS centers (JPL, CSR, and GFZ), and (c) Level‐3 mascon solutions. We show that inter‐satellite laser ranging residuals are modified at low frequencies below 1 mHz, affecting the along‐orbit analysis of large‐scale time‐variable gravity signals. When mapped into monthly Level‐2 spherical harmonic coefficients of geopotential, the low‐frequency change in inter‐satellite ranging residuals leads to substantial improvement of coefficients associated with resonant orders (i.e., 15, 30, 45, etc.) and C30. We also present an improved SLR‐derived C30 which significantly improves the agreement with updated GRACE‐FO C30 at seasonal and interannual timescales. Moreover, we demonstrate the noise reduction in mass change estimates from new GRACE‐FO Level‐2 data over oceans, Greenland, and Antarctica for all SDS solutions. GRACE‐FO mascon solutions show a moderate change in the updated release. Our comprehensive analyses demonstrate high‐quality estimates of non‐gravitational accelerations by the updated transplant method, resulting in more accurate GRACE‐FO time‐variable gravity and mass change observations. Plain Language Summary: GRACE Follow‐On (GRACE‐FO) consists of two satellites that orbit around the Earth at ∼500 km altitude, following each other at a distance of ∼200 km. The inter‐satellite distance between the two satellites changes because of the gravitational as well as non‐gravitational forces (e.g., air drag), and it is measured by high‐precision sensors. To estimate time‐variable gravity caused by mass change in the Earth system from inter‐satellite ranging data, the effect of non‐gravitational forces is reduced using on‐board accelerometer measurements. The accelerometer on‐board the GRACE‐D satellite started to return erroneous measurements about 1 month after the launch. The GRACE‐FO Science Data System (SDS) compensated this issue by transplanting the GRACE‐C accelerometer data toward that of GRACE‐D. An updated transplant accelerometer data set for GRACE‐D was recently released by the GRACE‐FO SDS and, consequently, new time‐variable gravity and mass change solutions were released. In this paper, we demonstrate the improved accuracy of the new GRACE‐FO time‐variable gravity and mass change solutions caused by high‐quality estimates of non‐gravitational accelerations by the updated transplant method. Key Points: GRACE Follow‐On (GRACE‐FO) laser ranging residuals are modified at low frequencies by ACH, affecting the along‐orbit analysis of large‐scale mass changesGeopotential spherical harmonic coefficients associated with resonant orders (i.e., ∼15, 30, 45, etc.) show the largest improvement by ACHMonthly gravity fields show significantly lower noise for months when the angle between GRACE‐FO orbital plane and Sun‐Earth direction is ∼0 [ABSTRACT FROM AUTHOR]

Large‐scale subsidence due to aquifer‐overdraft is an ongoing hazard in the San Joaquin Valley. Subsidence continues to cause damage to infrastructure and increases the risk of extensional fissures.Here, we use InSAR‐derived vertical land motion (VLM) to model the volumetric strain rate due to groundwater storage change during the 2007–2010 drought in the San Joaquin Valley, Central Valley, California. We then use this volumetric strain rate model to calculate surface tensile stress in order to predict regions that are at the highest risk for hazardous tensile surface fissures. We find a maximum volumetric strain rate of −232 microstrain/yr at a depth of 0 to 200 m in Tulare and Kings County, California. The highest risk of tensile fissure development occurs at the periphery of the largest subsiding zones, particularly in Tulare County and Merced County. Finally, we assume that subsidence is likely due to aquifer pressure change, which is calculated using groundwater level changes observed at 300 wells during this drought. We combine pressure data from selected wells with our volumetric strain maps to estimate the quasi‐static bulk modulus, K, a poroelastic parameter applicable when pressure change within the aquifer is inducing volumetric strain. This parameter is reflective of a slow deformation process and is one to two orders of magnitude lower than typical values for the bulk modulus found using seismic velocity data. The results of this study highlight the importance of large‐scale, high‐resolution VLM measurements in evaluating aquifer system dynamics, hazards associated with overdraft, and in estimating important poroelastic parameters. Key Points: Models of volume strain rate are developed for the San Joaquin Valley during the 2007–2010 droughtHigh risk of fissure initiation occurs on the periphery of subsiding zones based on the ratio of tensile stress to soil tensile strengthThe quasi‐static bulk modulus is estimated to be 0.19 GPa [ABSTRACT FROM AUTHOR]

California's millennium drought of 2012–2015 severely impacted the Central Valley aquifer system and caused permanent loss of groundwater and aquifer storage capacity. To quantify these impacts within the southern San Joaquin Valley, we analyze various complementary measurements, including gravity changes from Gravity Recovery and Climate Experiment (GRACE) satellites; vertical land motion from Global Positioning System, interferometric synthetic aperture radar, and extensometer; and groundwater level records. The interferometric data set acquired by the Sentinel‐1 satellite only spans the period January 2015 and October 2017, while the other data sets span the entire drought period. Using GRACE observations, we find an average groundwater loss of 6.1 ± 2.3 km3/year as a lower bound estimate for the San Joaquin Valley, amounting to a total volume of 24.2 ± 9.3 km3 lost during the period October 2011 to September 2015. This is consistent with the total volume of 29.25 ± 8.7 km3, estimated using only Global Positioning System deformation data. Our results highlight the advantage of using vertical land motion data to evaluate groundwater loss and thus fill the gaps between GRACE and GRACE‐Follow‐On missions and complement their estimates. We further determine that 0.4–3.25% of the aquifer system storage capacity is permanently lost during this drought period. Comparing groundwater level and vertical land motion data following September 2015, we determine an equilibration time of 0.5–1.5 years for groundwater levels within aquitard and aquifer units, during which residual compaction of aquitard and land subsidence continues beyond the drought period. We suggest that such studies can advance the knowledge of evolving groundwater resources, enabling managers and decision makers to better assess water demand and supply during and in‐between drought periods. Key Points: GPS measurements of aquifer compaction in San Joaquin Valley are analyzed to quantify groundwater loss during 2012–2015The GPS‐based estimate of groundwater loss is consistent with that obtained from satellite gravimetry, given uncertaintiesThe San Joaquin Valley aquifer system permanently lost up to 3% of its storage capacity due to overdraft during drought [ABSTRACT FROM AUTHOR]