Search for: All records

Abstract Drought intensity is commonly characterized using meteorologically‐based metrics that do not provide insight into water deficits within deeper hydrologic systems. In contrast, global positioning system (GPS) displacements are sensitive to both local and regional hydrologic‐storage fluctuations. While a few studies have leveraged this sensitivity to produce geodetic drought indices, hydrologic drought characterization using GPS is not commonly accounted for in drought assessment and management. To motivate this application, we produce a new geodetic drought index (GDI) and quantify its ability to characterize hydrologic drought conditions in key surface and sub‐surface hydrologic reservoirs/pools across California. In northern California, the GDI exhibits a strong regional association with surface‐reservoir storage at the 1‐month time scale (correlation coefficient: 0.83) and groundwater levels at the 3‐month time scale (correlation coefficient: 0.87), along with moderate associations with stream discharge at the daily (instantaneous) time scale (correlation coefficient: 0.50). Groundwater in southern California is best characterized with a 12‐month GDI (correlation coefficient: 0.77), and surface‐reservoir storage is optimized with the 3‐month GDI (correlation coefficient: 0.72). Two sigma uncertainties are ±0.03. Differences between northern and southern California reveal that the GDI is sensitive to unique aquifer and drainage basin characteristics. In addition to capturing long‐term hydrologic trends, rapid changes in the GDI initiate during clusters of large atmospheric river events that closely mirror fluctuations in traditional hydrologic and meteorological observations. We show that GPS‐based hydrologic drought indices provide a significant opportunity to improve drought assessment, in California and beyond, by improving our understanding of the hydrologic cycle. 

The data provided here accompany the publication "Drought Characterization with GPS: Insights into Groundwater and Reservoir Storage in California" [Young et al., (2024)] which is currently under review with Water Resources Research. (as of 28 May 2024)Please refer to the manuscript and its supplemental materials for full details. (A link will be appended following publication)File formatting information is listed below, followed by a sub-section of the text describing the Geodetic Drought Index Calculation. The longitude, latitude, and label for grid points are provided in the file "loading_grid_lon_lat".Time series for each Geodetic Drought Index (GDI) time scale are provided within "GDI_time_series.zip".The included time scales are for 00- (daily), 1-, 3-, 6-, 12- 18- 24-, 36-, and 48-month GDI solutions.Files are formatted following...Title: "grid point label L****"_"time scale"_monthFile Format: ["decimal date" "GDI value"]Gridded, epoch-by-epoch, solutions for each time scale are provided within "GDI_grids.zip".Files are formatted following...Title: GDI_"decimal date"_"time scale"_monthFile Format: ["longitude" "latitude" "GDI value" "grid point label L****"]2.2 GEODETIC DROUGHT INDEX CALCULATION We develop the GDI following Vicente-Serrano et al. (2010) and Tang et al. (2023), such that the GDI mimics the derivation of the SPEI, and utilize the log-logistic distribution (further details below). While we apply hydrologic load estimates derived from GPS displacements as the input for this GDI (Figure 1a-d), we note that alternate geodetic drought indices could be derived using other types of geodetic observations, such as InSAR, gravity, strain, or a combination thereof. Therefore, the GDI is a generalizable drought index framework. A key benefit of the SPEI is that it is a multi-scale index, allowing the identification of droughts which occur across different time scales. For example, flash droughts (Otkin et al., 2018), which may develop over the period of a few weeks, and persistent droughts (>18 months), may not be observed or fully quantified in a uni-scale drought index framework. However, by adopting a multi-scale approach these signals can be better identified (Vicente-Serrano et al., 2010). Similarly, in the case of this GPS-based GDI, hydrologic drought signals are expected to develop at time scales that are both characteristic to the drought, as well as the source of the load variation (i.e., groundwater versus surface water and their respective drainage basin/aquifer characteristics). Thus, to test a range of time scales, the TWS time series are summarized with a retrospective rolling average window of D (daily with no averaging), 1, 3, 6, 12, 18, 24, 36, and 48-months width (where one month equals 30.44 days). From these time-scale averaged time series, representative compilation window load distributions are identified for each epoch. The compilation window distributions include all dates that range ±15 days from the epoch in question per year. This allows a characterization of the estimated loads for each day relative to all past/future loads near that day, in order to bolster the sample size and provide more robust parametric estimates [similar to Ford et al., (2016)]; this is a key difference between our GDI derivation and that presented by Tang et al. (2023). Figure 1d illustrates the representative distribution for 01 December of each year at the grid cell co-located with GPS station P349 for the daily TWS solution. Here all epochs between between 16 November and 16 December of each year (red dots), are compiled to form the distribution presented in Figure 1e. This approach allows inter-annual variability in the phase and amplitude of the signal to be retained (which is largely driven by variation in the hydrologic cycle), while removing the primary annual and semi-annual signals. Solutions converge for compilation windows >±5 days, and show a minor increase in scatter of the GDI time series for windows of ±3-4 days (below which instability becomes more prevalent). To ensure robust characterization of drought characteristics, we opt for an extended ±15-day compilation window. While Tang et al. (2023) found the log-logistic distribution to be unstable and opted for a normal distribution, we find that, by using the extended compiled distribution, the solutions are stable with negligible differences compared to the use of a normal distribution. Thus, to remain aligned with the SPEI solution, we retain the three-parameter log-logistic distribution to characterize the anomalies. Probability weighted moments for the log-logistic distribution are calculated following Singh et al., (1993) and Vicente-Serrano et al., (2010). The individual moments are calculated following Equation 3. These are then used to calculate the L-moments for shape (), scale (), and location () of the three-parameter log-logistic distribution (Equations 4 – 6). The probability density function (PDF) and the cumulative distribution function (CDF) are then calculated following Equations 7 and 8, respectively. The inverse Gaussian function is used to transform the CDF from estimates of the parametric sample quantiles to standard normal index values that represent the magnitude of the standardized anomaly. Here, positive/negative values represent greater/lower than normal hydrologic storage. Thus, an index value of -1 indicates that the estimated load is approximately one standard deviation dryer than the expected average load on that epoch. *Equations can be found in the main text. 

Abstract Atmospheric rivers (ARs) deliver significant and essential precipitation to the western United States (US) with consequential interannual variability. The intensity and frequency of ARs strongly influence reservoir levels, mountain snowpack, and groundwater recharge, which are key drivers of water‐resource availability and natural hazards. Between October 2022 and April 2023, western states experienced exceptionally heavy precipitation from several families of powerful ARs. Using observations of surface‐loading deformation from Global Navigation Satellite Systems, we find that terrestrial water‐storage gains exceeded 100% of normal within vital California watersheds. Independent water‐storage solutions derived from different data‐analysis and inversion methods provide an important measure of precision. The sustained storage increases, which we show are closely associated with ARs at daily‐to‐weekly timescales, alleviated both meteorological and hydrological drought conditions in the region, with a lag in hydrological‐drought improvements. Quantifying water‐storage recovery associated with extreme precipitation after drought advances understanding of an increasingly variable hydrologic cycle. 

Abstract Groundwater flowing within mountainous regions is a critical water source for mountain and adjacent low‐elevation communities and ecosystems. Despite their importance, quantifying storage changes within these groundwater systems remains a challenge. To evaluate the fate of mountain groundwater and the role of extremely wet winters in replenishing groundwater storage, we use high‐precision Global Navigation Satellite Systems surface displacements to constrain changes in groundwater storage within the Sierra Nevada and Cascade mountain ranges, two important mountainous regions spanning much of the westernmost US. We find that groundwater storage has significantly declined over the past two decades, comprising over 90% of the total water storage lost in these regions, associated with recurrent, severe, multi‐year droughts. However, extremely wet winters can replenish groundwater by more than twice the average annual flux over short periods, driving these systems from historical lows to above normal conditions. Furthermore, discharge from these groundwater systems is 2–3 times less variable than recharge, indicating they release a relatively constant amount of water year over year, allowing for surplus storage following extremely wet winter periods. These findings show that extreme precipitation (snow + rain) rapidly recharges mountain groundwater and can offset multi‐year losses. Moreover, mountain groundwater systems can retain storage gains for at least 1 year post‐event, providing a durable source of streamflow, agricultural, and municipal water supply. As extremely wet winters are predicted to become more intense and frequent, we hypothesize they may help offset groundwater loss due to aridification, sustaining ecosystems and buffering against a new, drier climate normal. 

Abstract We installed a purpose‐built network of co‐located Global Navigation Satellite System (GNSS) stations and meteorological instrumentation to investigate water storage in a high‐mountain watershed along the Idaho‐Montana border. Twelve GNSS stations are distributed across the Selway‐Lochsa watersheds at approximately 30–40 km spacing, filling a critical observational gap between localized point measurements and regional geodetic and satellite data sets. The unique coupling of geodetic and hydrologic observations in this network enables direct comparison between co‐located GNSS measurements of the elastic response of the solid Earth and local changes in measured water storage. This network is specifically designed to address questions of hydrologic storage and movement at the mountain watershed scale. Here, we describe technical details of the network and its deployment; introduce new hydrologic, meteorologic, and geodetic data sets recorded by the network; process and analyze the source data (e.g., time series of daily three‐dimensional GNSS site positions, removal of non‐hydrologic signals); and characterize basic empirical relationships between water storage, water movement, and GNSS‐inferred surface displacement. The network shows preliminary evidence for spatial differences in displacement resulting from a range of snow loads across elevations, but longer and more complete data records are needed to support these initial findings. We also provide examples of additional scientific applications of this network, including estimations of snow depth and snow water equivalent from GNSS multipath reflectometry. Finally, we consider the challenges, limitations, and opportunities of deploying GNSS and weather stations at high elevations with heavy snowpack and offer ideas for technical improvements. 

We investigate deformation mechanics of fracture networks in unsaturated fractured rocks from subsurface conventional detonation using dynamic noble gas measurements and changes in air permeability. We dynamically measured the noble gas isotopic composition and helium exhalation of downhole gas before and after a large subsurface conventional detonation. These noble gas measurements were combined with measurements of the subsurface permeability field from 64 discrete sampling intervals before and after the detonation and subsurface mapping of fractures in borehole walls before well completion. We saw no observable increase in radiogenic noble gas release from either an isotopic composition or a helium exhalation point of view. Large increases in permeability were observed in 13 of 64 discrete sampling intervals. Of the sampling intervals which saw large increases in flow, only two locations did not have preexisting fractures mapped at the site. Given the lack of noble gas release and a clear increase in permeability, we infer that most of the strain accommodation of the fractured media occurred along previously existing fractures, rather than the creation of new fractures, even for a high strain rate event. These results have significant implications for how we conceptualize the deformation of rocks with fracture networks above the percolation threshold, with application to a variety of geologic and geological engineering problems.