Vertical and horizontal components of GNSS displacements in the Long Valley Caldera and adjacent Sierra Nevada range show a clear correlation with hydrological trends at both multiyear and seasonal time scales. We observe a clear vertical and horizontal seasonal deformation pattern primarily attributable to the solid earth response to hydrological surface loading at large-to-regional (Sierra Nevada range) scales. Several GNSS sites, located at the eastern edge of the Sierra Nevada along the southwestern rim of Long Valley Caldera, also show significant horizontal deformation that cannot be explained by elastic deformation from surface loading. Due to the location of these sites and the strong correlation between their horizontal displacements and spring discharge, we hypothesize that this deformation reflects poroelastic processes related to snowmelt runoff water infiltrating into the Sierra Nevada slopes around Long Valley Caldera. Interestingly, this is also an area where water infiltrates to feed the local hydrothermal system, and where snowmelt-induced earthquake swarms have been recently detected.

Increased climate variability is driving changes in water storage across the contiguous United States (CONUS). Observational estimates of these storage changes are important for validation of hydrological models and predicting future water availability. We estimate CONUS terrestrial water storage anomalies (TWSA) from 2007–2017 using Global Positioning System (GPS) displacements, constrained by lower-resolution TWSA observations from Gravity Recovery and Climate Experiment (GRACE) satellite gravity—a combination that provides higher spatiotemporal resolution than previous estimates. The relative contribution of seasonal, interannual, and subseasonal TWSA varies widely across CONUS watersheds, with implications for regional water security. Separately, we find positive correlation between TWSA and the El Niño/Southern Oscillation in the southeastern Texas-Gulf and South Atlantic-Gulf watersheds and an unexpected negative correlation in the southwest. In the western United States, atmospheric rivers (ARs) drive a large fraction of subseasonal TWSA, with the top 5% of ARs contributing 73% of total AR-related TWSA increases.

Increasingly accurate, and spatio-temporally dense, measurements of Earth surface movements enable us to identify multiple deformation patterns and highlight the need to properly characterize the related source processes. This is particularly important in tectonically active areas, where deformation measurement is crucial for monitoring ongoing processes and assessing future hazard. Long Valley Caldera, California, USA, is a volcanic area where frequent episodes of unrest involve inflation and increased seismicity. Ground- and satellite-based instruments show that volcanic inflation renewed in 2011, and is continuing as of early 2021. Additionally, Long Valley Caldera is affected by the large, but spatially and temporally variable, amounts of precipitation falling on the adjacent Sierra Nevada Range. The density and long duration of deformation measurements at Long Valley Caldera provide an excellent collection of data to decompose time-series and separate multiple superimposed deformation sources. We analyze Global Navigation Satellite System (GNSS) time-series and apply variational Bayesian Independent Component Analysis (vbICA) decomposition method to isolate inflation-related signals from other processes. We show that hydrological forcing causes significant horizontal and vertical deformation at different temporal (seasonal and multiyear) and spatial (few to hundreds of km) scales that cannot be ignored while analyzing and modeling the tectonic signal. Focusing on the last inflation episode, we then improve on prior simplistic models of the inflation reservoir by including heterogeneous subsurface material properties and topography. Our results suggest the persistence and stability of the reservoir (prolate ellipsoid at about 8 km beneath the resurgent dome) and indicate a 40-50% reduction of the inflation rate after about 3 years from the inflation onset. The onset of the reduced inflation rate corresponded in time with the occurrence of a strong seismic swarm in the Caldera, but also to the temporal variation of climatic conditions in the area.