Measuring crustal strain and seismic moment accumulation, is crucial for understanding the growth and distribution of seismic hazards along major fault systems. Here, we develop a methodology to integrate 4.5 years (2015–2019.5) of Sentinel‐1 Interferometric Synthetic Aperture Radar (InSAR) and continuous Global Navigation Satellite System (GNSS) time series to achieve 6 to 12‐day sampling of surface displacements at ∼500 m spatial resolution over the entire San Andreas fault system. Numerous interesting deformation signals are identified with this product (video link:
Global Navigation Satellite System (GNSS) vertical displacements measuring the elastic response of Earth's crust to changes in hydrologic mass have been used to produce terrestrial water storage change (∆TWS) estimates for studying both annual ∆TWS as well as multi‐year trends. However, these estimates require a high observation station density and minimal contamination by nonhydrologic deformation sources. The Gravity Recovery and Climate Experiment (GRACE) is another satellite‐based measurement system that can be used to measure regional TWS fluctuations. The satellites provide highly accurate ∆TWS estimates with global coverage but have a low spatial resolution of ∼400 km. Here, we put forward the mathematical framework for a joint inversion of GNSS vertical displacement time series with GRACE ∆TWS to produce more accurate spatiotemporal maps of ∆TWS, accounting for the observation errors, data gaps, and nonhydrologic signals. We aim to utilize the regional sensitivity to ∆TWS provided by GRACE mascon solutions with higher spatial resolution provided by GNSS observations. Our approach utilizes a continuous wavelet transform to decompose signals into their building blocks and separately invert for long‐term and short‐term mass variations. This allows us to preserve trends, annual, interannual, and multi‐year changes in TWS that were previously challenging to capture by satellite‐based measurement systems or hydrological models, alone. We focus our study in California, USA, which has a dense GNSS network and where recurrent, intense droughts put pressure on freshwater supplies. We highlight the advantages of our joint inversion results for a tectonically active study region by comparing them against inversion results that use only GNSS vertical deformation as well as with maps of ∆TWS from hydrological models and other GRACE solutions. We find that our joint inversion framework results in a solution that is regionally consistent with the GRACE ∆TWS solutions at different temporal scales but has an increased spatial resolution that allows us to differentiate between regions of high and low mass change better than using GRACE alone.
more » « less- Award ID(s):
- 2054160
- NSF-PAR ID:
- 10374466
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 127
- Issue:
- 3
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract https://www.youtube.com/watch?v=SxNLQKmHWpY ). We decompose the line‐of‐sight InSAR displacements into three dimensions by combining the deformation azimuth from a GNSS‐derived interseismic fault model. We then construct strain rate maps using a smoothing interpolator with constraints from elasticity. The resulting deformation field reveals a wide array of crustal deformation processes including, on‐ and off‐fault secular and transient tectonic deformation, creep rates on all the major faults, and vertical signals associated with hydrological processes. The strain rate maps show significant off‐fault components that were not captured by GNSS‐only models. These results are important in assessing the seismic hazard in the region. -
more » « less
Abstract Accurate and timely inland waterbody extent and location data are foundational information to support a variety of hydrological applications and water resources management. Recently, the Cyclone Global Navigation Satellite System (CYGNSS) has emerged as a promising tool for delineating inland water due to distinct surface reflectivity characteristics over dry versus wet land which are observable by CYGNSS’s eight microsatellites with passive bistatic radars that acquire reflected L-band signals from the Global Positioning System (GPS) (i.e., signals of opportunity). This study conducts a baseline 1-km comparison of water masks for the contiguous United States between latitudes of 24°N-37°N for 2019 using three Earth observation systems: CYGNSS (i.e., our baseline water mask data), the Moderate Resolution Imaging Spectroradiometer (MODIS) (i.e., land water mask data), and the Landsat Global Surface Water product (i.e., Pekel data). Spatial performance of the 1-km comparison water mask was assessed using confusion matrix statistics and optical high-resolution commercial satellite imagery. When a mosaic of binary thresholds for 8 sub-basins for CYGNSS data were employed, confusion matrix statistics were improved such as up to a 34% increase in F1-score. Further, a performance metric of ratio of inland water to catchment area showed that inland water area estimates from CYGNSS, MODIS, and Landsat were within 2.3% of each other regardless of the sub-basin observed. Overall, this study provides valuable insight into the spatial similarities and discrepancies of inland water masks derived from optical (visible) versus radar (Global Navigation Satellite System Reflectometry, GNSS-R) based satellite Earth observations.
-
Abstract The West Antarctic Ice Sheet (WAIS) overlies a thin, variable-thickness lithosphere and a shallow upper-mantle region of laterally varying and, in some regions, very low (~1018 Pa s) viscosity. We explore the extent to which viscous effects may affect predictions of present-day geoid and crustal deformation rates resulting from Antarctic ice mass flux over the last quarter century and project these calculations into the next half century, using viscoelastic Earth models of varying complexity. Peak deformation rates at the end of a 25-yr simulation predicted with an elastic model underestimate analogous predictions that are based on a 3D viscoelastic Earth model (with minimum viscosity below West Antarctica of 1018 Pa s) by ~15 and ~3 mm yr−1 in the vertical and horizontal directions, respectively, at sites overlying low-viscosity mantle and close to high rates of ice mass flux. The discrepancy in uplift rate can be reduced by adopting 1D Earth models tuned to the regional average viscosity profile beneath West Antarctica. In the case of horizontal crustal rates, adopting 1D regional viscosity models is no more accurate in recovering predictions that are based on 3D viscosity models than calculations that assume a purely elastic Earth. The magnitude and relative contribution of viscous relaxation to crustal deformation rates will likely increase significantly in the next several decades, and the adoption of 3D viscoelastic Earth models in analyses of geodetic datasets [e.g., Global Navigation Satellite System (GNSS); Gravity Recovery and Climate Experiment (GRACE)] will be required to accurately estimate the magnitude of Antarctic modern ice mass flux in the progressively warming world.more » « less
-
Zhang, Kai (Ed.)Realistic representation of hydrological drought events is increasingly important in world facing decreased freshwater availability. Index-based drought monitoring systems are often adopted to represent the evolution and distribution of hydrological droughts, which mainly rely on hydrological model simulations to compute these indices. Recent studies, however, indicate that model derived water storage estimates might have difficulties in adequately representing reality. Here, a novel Markov Chain Monte Carlo - Data Assimilation (MCMC-DA) approach is implemented to merge global Terrestrial Water Storage (TWS) changes from the Gravity Recovery And Climate Experiment (GRACE) and its Follow On mission (GRACE-FO) with the water storage estimations derived from the W3RA water balance model. The modified MCMC-DA derived summation of deep rooted soil and groundwater storage estimates is then used to compute 0.5∘ standardized groundwater drought indices globally to show the impact of GRACE/GRACE-FO DA on a global index-based hydrological drought monitoring system. Our numerical assessment covers the period of 2003–2021, and shows that integrating GRACE/GRACE-FO data modifies the seasonality and inter-annual trends of water storage estimations. Considerable increases in the length and severity of extreme droughts are found in basins that exhibited multiyear water storage fluctuations and those affected by climate teleconnections.more » « less
-
more » « less
Abstract Geodetic methods can monitor changes in terrestrial water storage (TWS) across large regions in near real‐time. Here, we investigate the effect of assumed Earth structure on TWS estimates derived from Global Navigation Satellite System (GNSS) displacement time series. Through a series of synthetic tests, we systematically explore how the spatial wavelength of water load affects the error of TWS estimates. Large loads (e.g., >1,000 km) are well recovered regardless of the assumed Earth model. For small loads (e.g., <10 km), however, errors can exceed 75% when an incorrect model for the Earth is chosen. As a case study, we consider the sensitivity of seasonal TWS estimates within mountainous watersheds of the western U.S., finding estimates that differ by over 13% for a collection of common global and regional structural models. Errors in the recovered water load generally scale with the total weight of the load; thus, long‐term changes in storage can produce significant uplift (subsidence), enhancing errors. We demonstrate that regions experiencing systematic and large‐scale variations in water storage, such as the Greenland ice sheet, exhibit significant differences in predicted displacement (over 20 mm) depending on the choice of Earth model. Since the discrepancies exceed GNSS observational precision, an appropriate Earth model must be adopted when inverting GNSS observations for mass changes in these regions. Furthermore, regions with large‐scale mass changes that can be quantified using independent data (e.g., altimetry, gravity) present opportunities to use geodetic observations to refine structural properties of seismologically derived models for the Earth's interior structure.
