Hydrogeodesy, a relatively new field within the earth sciences, is the analysis of the distribution and movement of terrestrial water at Earth's surface using measurements of Earth's shape, orientation, and gravitational field. In this paper, we review the current state of hydrogeodesy with a specific focus on Global Navigation Satellite System (GNSS)/Global Positioning System measurements of hydrologic loading. As water cycles through the hydrosphere, GNSS stations anchored to Earth's crust measure the associated movement of the land surface under the weight of changing hydrologic loads. Recent advances in GNSS‐based hydrogeodesy have led to exciting applications of hydrologic loading and subsequent terrestrial water storage (TWS) estimates. We describe how GNSS position time series respond to climatic drivers, can be used to estimate TWS across temporal scales, and can improve drought characterization. We aim to facilitate hydrologists' use of GNSS‐observed surface deformation as an emerging tool for investigating and quantifying water resources, propose methods to further strengthen collaborative research and exchange between geodesists and hydrologists, and offer ideas about pressing questions in hydrology that GNSS may help to answer.
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.
more » « less- Award ID(s):
- 2021637
- NSF-PAR ID:
- 10480752
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Earth and Space Science
- Volume:
- 10
- Issue:
- 10
- ISSN:
- 2333-5084
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract In snowmelt‐driven mountain watersheds, the hydrologic connectivity between meteoric waters and stream flow generation varies strongly with the season, reflecting variable connection to soil and groundwater storage within the watershed. This variable connectivity regulates how streamflow generation mechanisms transform the seasonal and elevational variation in oxygen and hydrogen isotopic composition (δ18O and δD) of meteoric precipitation. Thus, water isotopes in stream flow can signal immediate connectivity or more prolonged mixing, especially in high‐relief mountainous catchments. We characterized δ18O and δD values in stream water along an elevational gradient in a mountain headwater catchment in southwestern Montana. Stream water isotopic compositions related most strongly to elevation between February and March, exhibiting higher δ18O and δD values with decreasing elevation. These elevational isotopic lapse rates likely reflect increased connection between stream flow and proximal snow‐derived water sources heavily subject to elevational isotopic effects. These patterns disappeared during summer sampling, when consistently lower δ18O and δD values of stream water reflected contributions from snowmelt or colder rainfall, despite much higher δ18O and δD values expected in warmer seasonal rainfall. The consistently low isotopic values and absence of a trend with elevation during summer suggest lower connectivity between summer precipitation and stream flow generation as a consequence of drier soils and greater transpiration. As further evidence of intermittent seasonal connectivity between the stream and adjacent groundwaters, we observed a late‐winter flush of nitrate into the stream at higher elevations, consistent with increased connection to accumulating mineralized nitrogen in riparian wetlands. This pattern was distinct from mid‐summer patterns of nitrate loading at lower elevations that suggested heightened human recreational activity along the stream corridor. These observations provide insights linking stream flow generation and seasonal water storage in high elevation mountainous watersheds. Greater understanding of the connections between surface water, soil water and groundwater in these environments will help predict how the quality and quantity of mountain runoff will respond to changing climate and allow better informed water management decisions.
-
Abstract. Seasonal snowpack is an essential component in the hydrological cycle and plays a significant role in supplying water resources to downstream users. Yet the snow water equivalent (SWE) in seasonal snowpacks, and its space–time variation, remains highly uncertain, especially over mountainous areas with complex terrain and sparse observations, such as in High Mountain Asia (HMA). In this work, we assessed the spatiotemporal distribution of seasonal SWE, obtained from a new 18-year HMA Snow Reanalysis (HMASR) dataset, as part of the recent NASA High Mountain Asia Team (HiMAT) effort. A Bayesian snow reanalysis scheme previously developed to assimilate satellite-derived fractional snow-covered area (fSCA) products from Landsat and MODIS platforms has been applied to develop the HMASR dataset (at a spatial resolution of 16 arcsec (∼500 m) and daily temporal resolution) over the joint Landsat–MODIS period covering water years (WYs) 2000–2017. Based on the results, the HMA-wide total SWE volume is found to be around163 km3 on average and ranges from 114 km3 (WY2001) to 227 km3 (WY2005) when assessed over 18 WYs. The most abundant snowpacks are found in the northwestern basins (e.g., Indus, Syr Darya and Amu Darya) that are mainly affected by the westerlies, accounting for around 66 % of total seasonal SWE volume. Seasonal snowpack in HMA is depicted by snow accumulating through October to March and April, typically peaking around April and depleting in July–October, with variations across basins and WYs. When examining the elevational distribution over the HMA domain, seasonal SWE volume peaks at mid-elevations (around 3500 m), with over 50 % of the volume stored above 3500 m. Above-average amounts of precipitation causes significant overall increase in SWE volumes across all elevations, while an increase in air temperature (∼1.5 K) from cooler to normal conditions leads to an redistribution in snow storage from lower elevations to mid-elevations. This work brings new insight into understanding the climatology and variability of seasonal snowpack over HMA, with the regional snow reanalysis constrained by remote-sensing data, providing a new reference dataset for future studies of seasonal snow and how it contributes to the water cycle and climate over the HMA region.more » « less
-
more » « less
Abstract The seasonal snowmelt period is a critical component of the hydrologic cycle for many mountainous areas. Changes in the timing and rate of snowmelt as a result of physical hydrologic flow paths, such as longitudinal intra‐snowpack flow paths, can have strong implications on the partitioning of meltwater amongst streamflow, groundwater recharge, and soil moisture storage. However, intra‐snowpack flow paths are highly spatially and temporally variable and thus difficult to observe. This study utilizes new methods to non‐destructively observe spatio‐temporal changes in the liquid water content of snow in combination with plot experiments to address the research question: What is the scale of influence that intra‐snowpack flow paths have on the downslope movement of liquid water during snowmelt across an elevational gradient? This research took place in northern Colorado with study plots spanning from the rain‐snow transition zone up to the high alpine. Results indicate an increasing scale of influence from intra‐snowpack flow paths with elevation, showing higher hillslope connectivity producing larger intra‐snowpack contributing areas for meltwater accumulation, quantified as the upslope contributing area required to produce observed changes in liquid water content from melt rate estimates. The total effective intra‐snowpack contributing area of accumulating liquid water was found to be 17, 6, and 0 m2for the above tree line, near tree line, and below tree line plots, respectively. Dye tracer experiments show capillary and permeability barriers result in increased number and thickness of intra‐snowpack flow paths at higher elevations. We additionally utilized aerial photogrammetry in combination with ground penetrating radar surveys to investigate the role of this hydrologic process at the small watershed scale. Results here indicate that intra‐snowpack flow paths have influence beyond the plot scale, impacting the storage and transmission of liquid water within the snowpack at the small watershed scale.
-
more » « less
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.
