An official website of the United States government
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.
more »
« less
Identification and Extraction of Seasonal Geodetic Signals Due to Surface Load Variations
Abstract Deformation of the Earth's surface associated with redistributions of continental water mass explains, to first order, the seasonal signals observed in geodetic position time series. Discriminating these seasonal signals from other sources of deformation in geodetic measurements is essential to isolate tectonic signals and to monitor spatio‐temporal variations in continental water storage. We propose a new methodology to identify and extract these seasonal signals. The approach uses a variational Bayesian Independent Component Analysis (vbICA) to extract the seasonal signals and a gravity‐based deformation model to identify which of these signals are caused by surface loading. We test the procedure on two study areas, the Arabian Peninsula and the Nepal Himalaya, and find that the technique successfully extracts the seasonal signals with one or two independent components, depending on whether the load is stationary or moving. The approach is robust to spatial heterogeneities inherent to geodetic measurements and can help extract systematic errors in geodetic products (e.g., draconitic errors). We also discuss how to handle the degree‐1 deformation field present in the geodetic data set but not captured by the gravity‐based model.
more »
« less
- Award ID(s):
- 1821853
- PAR ID:
- 10370790
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 123
- Issue:
- 12
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Fractures within ice shelves are zones of weakness, which can deform on timescales from seconds to decades. Icequakes produced during the fracturing process show a higherb‐valuein the Gutenberg‐Richter scaling relationship than continental earthquakes. We investigate icequakes on the east side of rift WR4 in the Ross Ice Shelf, Antarctica. Our model suggests a maximum icequake slip depth that is ∼7.8 m below the rift mélange, where the slip area can only grow laterally along the fracture planes. We propose ductile deformation below this depth, potentially due to the saturation of unfrozen water. We use remote sensing and geodetic tools to quantify surface movement on different timescales and find that the majority of icequakes occur during falling tides. The total seismic moment is <1% of the estimated geodetic moment during a tidal cycle. This study demonstrates the feasibility of using seismology and geodesy to investigate ice rift zone rheology.more » « less
-
Abstract Africa's continental crust hosts a variety of geologic terrains and is crucial for understanding the evolution of its longest‐lived cratons. However, few of its seismological models are yet to incorporate the largest continent‐wide noise dispersion data sets. Here, we report on new insights into Africa's crustal architecture obtained using a new data set and model assessment product, ADAMA, which comprises a large ensemble of short‐period surface wave dispersion measurements: 5–40 s. We construct a continent‐wide model ofAfrica'sCrustEvaluated with ADAMA'sRayleighPhase maps (ACE‐ADAMA‐RP). Dispersion maps, and uncertainties, are obtained with a probabilistic approach. This model update, and a crustal taxonomy derived from unsupervised machine learning, reveals that the architecture of Africa's crust can be classified into two main types:primitive(C1: faster velocities with little gradients) andmodified(C2–C4: slower velocities in the shallow crust with more pronounced gradients). The Archean shields are “primitive,” showing little variation or secular evolution. The basins, orogens, and continental margins are “modified” and retain imprints of surface deformation. The crustal taxonomy is obtained without a‐priori geological information and differs from previous classification schemes. While most of our reported features are robust, probabilistic modeling suggests caution in the quantitative interpretations where illumination is compromised by low‐quality measurements, sparse coverage or both. Future extension of our approach to other complementary seismological and geophysical data sets—for example, multimode earthquake dispersion, receiver functions, gravity, and mineral physics, will enable continent‐wide lithospheric modeling that extends resolution to the upper mantle.more » « less
-
Abstract Heat storage change (HSC) is a crucial component of lake's thermal energy budget. Conventional temperature profile based models of HSC require location specific parameters such as lakebed topography. Based on the half‐order time‐derivative formula of heat fluxes, an analytical model was formulated for estimating HSC from water surface temperature and solar radiation without using geography dependent parameters. The proposed model was tested against field measurements at Poyang Lake, a shallow inland lake, which has pronounced seasonal variations in water level and lake area. Our analysis indicates that the model accurately simulates diurnal HSC with a coefficient of determination of 0.94 and a root mean squared error (RMSE) of 77.5 ± 21.6 Wm−2for the study period. Larger nighttime RMSE (75.0 ± 26.8 Wm−2) than the daytime value (55.1 ± 19.7 W m−2) is attributable to larger measurement errors of nighttime turbulent fluxes. The estimation of HSC independent of temperature profile and lake‐specific parameters by the proposed model facilitates remote sensing monitoring the HSC of global water bodies.more » « less
-
ABSTRACT As part of the 2022 revision of the Aotearoa New Zealand National Seismic Hazard Model (NZ NSHM 2022), deformation models were constructed for the upper plate faults and subduction interfaces that impact ground-shaking hazard in New Zealand. These models provide the locations, geometries, and slip rates of the earthquake-producing faults in the NZ NSHM 2022. For upper plate faults, two deformation models were developed: a geologic model derived directly from the fault geometries and geologic slip rates in the NZ Community Fault Model version 1.0 (NZ CFM v.1.0); and a geodetic model that uses the same faults and fault geometries and derives fault slip-deficit rates by inverting geodetic strain rates for back slip on those specified faults. The two upper plate deformation models have similar total moment rates, but the geodetic model has higher slip rates on low-slip-rate faults, and the geologic model has higher slip rates on higher-slip-rate faults. Two deformation models are developed for the Hikurangi–Kermadec subduction interface. The Hikurangi–Kermadec geometry is a linear blend of the previously published interface models. Slip-deficit rates on the Hikurangi portion of the deformation model are updated from the previously published block models, and two end member models are developed to represent the alternate hypotheses that the interface is either frictionally locked or creeping at the trench. The locking state in the Kermadec portion is less well constrained, and a single slip-deficit rate model is developed based on plate convergence rate and coupling considerations. This single Kermadec realization is blended with each of the two Hikurangi slip-deficit rate models to yield two overall Hikurangi–Kermadec deformation models. The Puysegur subduction interface deformation model is based on geometry taken directly from the NZ CFM v.1.0, and a slip-deficit rate derived from published geodetic plate convergence rate and interface coupling estimates.more » « less
