Search for: All records

Abstract 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. 

Abstract Imaging tectonic creep along active faults is critical for measuring strain accumulation and ultimately understanding the physical processes that guide creep and the potential for seismicity. We image tectonic deformation along the central creeping section of the San Andreas Fault at the Dry Lake Valley paleoseismic site (36.468°N, 121.055°W) using three data sets with varying spatial and temporal scales: (1) an Interferometric Synthetic Aperture Radar (InSAR) velocity field with an ~100‐km footprint produced from Sentinel‐1 satellite imagery, (2) light detection and ranging (lidar) and structure‐from‐motion 3‐D topographic differencing that resolves a decade of deformation over a 1‐km aperture, and (3) surface fractures that formed over the 3‐ to 4‐m wide fault zone during a drought from late 2012 to 2014. The InSAR velocity map shows that shallow deformation is localized to the San Andreas Fault. We demonstrate a novel approach for differencing airborne lidar and structure‐from‐motion topography that facilitates resolving deformation along and adjacent to the San Andreas Fault. The 40‐m resolution topographic differencing resolves a 2.5 ± 0.2 cm/yr slip rate localized to the fault. The opening‐mode fractures accommodate cm/yr of fault slip. A 90% ± 30% of the 1‐km aperture deformation is accommodated over the several meter‐wide surface trace of the San Andreas Fault. The extension direction inferred from the opening‐mode fractures and topographic differencing is 40°–48° from the local trend of the San Andreas Fault. The localization of deformation likely reflects the well‐oriented and mature fault. 

Abstract We develop finite element models of the coseismic displacement field accounting for the 3D elastic structures surrounding the epicentral area of the 2019 Ridgecrest earthquake sequence containing two major events of M_w7.1 and M_w6.4. The coseismic slip distribution is inferred from the surface displacement field recorded by interferometric synthetic aperture radar. The rupture dip geometry is further optimized using a novel nonlinear‐crossover‐linear inversion approach. It is found that accounting for elastic heterogeneity and fault along‐strike curvilinearity improves the fit to the observed displacement field and yields a more accurate estimate of geodetic moment and Coulomb stress changes. We observe spatial correlations among the locations of aftershocks and patches of high slip, and rock anomalous elastic properties, suggesting that the shallow crust's elastic structures possibly controlled the Ridgecrest earthquake sequence. Most of the coseismic slip with a peak slip of 7.4 m at 3.6 km depth occurred above a zone of reducedS‐wave velocity and significant post‐M_w7.1 afterslip. This implies that viscous materials or fluid presence might have contributed to the low rupture velocity of the mainshock. Moreover, the zone of high slip on the northwest‐trending fault segment is laterally bounded by two aftershock clusters, whose location is characterized by intermediate rock rigidity. Notably, some minor orthogonal faults consistently end above a subsurface rigid body. Overall, these observations of structural controls improve our understandings of the seismogenesis within incipient fault systems.