skip to main content


Title: Subsurface Water Dominates Sierra Nevada Seasonal Hydrologic Storage
Abstract

Vertical displacements (dz) in permanent Global Positioning System (GPS) station positions enable estimation of water storage changes (ΔS), which historically have been impossible to measure directly. We use dz from 924 GPS stations in the western United States to estimate daily ΔS in California's Sierra Nevada and compare it to seasonal snow accumulation and melt over water years 2008–2017. Seasonal variations in GPS‐based ΔS are ~1,000 mm. Typically, only ~30% of ΔS is attributable to snow water equivalent (SWE). ΔS lags the snow cycle, peaking after maximum SWE and remaining positive when all snow has melted (SWE = 0). We conclude that seasonal ΔS fluctuations are not primarily driven by SWE but by rainfall and snowmelt stored in the shallow subsurface (as soil moisture and/or groundwater) and released predominantly through evapotranspiration. Seasonal peak GPS ΔS is larger than accumulated precipitation from the Parameter‐elevation Relationships on Independent Slopes Model and North American Land Data Assimilation System, indicating that these standard inputs underestimate mountain precipitation.

 
more » « less
NSF-PAR ID:
10448932
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
46
Issue:
21
ISSN:
0094-8276
Page Range / eLocation ID:
p. 11993-12001
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Accurate precipitation estimates are critical to simulating seasonal snowpack evolution. We conduct and evaluate high‐resolution (4‐km) snowpack simulations over the western United States (WUS) mountains in Water Year 2013 using the Noah with multi‐parameterization (Noah‐MP) land surface model driven by precipitation forcing from convection‐permitting (4‐km) Weather Research and Forecasting (WRF) modeling and four widely used high‐resolution datasets that are derived from statistical interpolation based on in situ measurements. Substantial differences in the precipitation amount among these five datasets, particularly over the western and northern portions of WUS mountains, significantly affect simulated snow water equivalent (SWE) and snow depth (SD) but have relatively limited effects on snow cover fraction (SCF) and surface albedo. WRF generally captures observed precipitation patterns and results in an overall best‐performed SWE and SD in the western and northern portions of WUS mountains, where the statistically interpolated datasets lead to underpredicted precipitation, SWE, and SD. Over the interior WUS mountains, all the datasets consistently underestimate precipitation, causing significant negative biases in SWE and SD, among which the results driven by the WRF precipitation show an average performance. Further analysis reveals systematic positive biases in SCF and surface albedo across the WUS mountains, with similar bias patterns and magnitudes for simulations driven by different precipitation datasets, suggesting an urgent need to improve the Noah‐MP snowpack physics. This study highlights that convection‐permitting modeling with proper configurations can have added values in providing decent precipitation for high‐resolution snowpack simulations over the WUS mountains in a typical ENSO‐neutral year, particularly over observation‐scarce regions.

     
    more » « less
  2. 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
  3. Abstract

    Climate change is altering biogeochemical, metabolic, and ecological functions in lakes across the globe. Historically, mountain lakes in temperate regions have been unproductive because of brief ice‐free seasons, a snowmelt‐driven hydrograph, cold temperatures, and steep topography with low vegetation and soil cover. We tested the relative importance of winter and summer weather, watershed characteristics, and water chemistry as drivers of phytoplankton dynamics. Using boosted regression tree models for 28 mountain lakes in Colorado, we examined regional, intraseasonal, and interannual drivers of variability in chlorophyllaas a proxy for lake phytoplankton. Phytoplankton biomass was inversely related to the maximum snow water equivalent (SWE) of the previous winter, as others have found. However, even in years with average SWE, summer precipitation extremes and warming enhanced phytoplankton biomass. Peak seasonal phytoplankton biomass coincided with the warmest water temperatures and lowest nitrogen‐to‐phosphorus ratios. Although links between snowpack, lake temperature, nutrients, and organic‐matter dynamics are increasingly recognized as critical drivers of change in high‐elevation lakes, our results highlight the additional influence of summer conditions on lake productivity in response to ongoing changes in climate. Continued changes in the timing, type, and magnitude of precipitation in combination with other global‐change drivers (e.g., nutrient deposition) will affect production in mountain lakes, potentially shifting these historically oligotrophic lakes toward new ecosystem states. Ultimately, a deeper understanding of these drivers and pattern at multiple scales will allow us to anticipate ecological consequences of global change better.

     
    more » « less
  4. Abstract

    The impacts of warming temperatures and declining snowpack on seasonal water yields in the Missouri River Headwaters are not well understood, revealing a gap in our understanding of regional hydroclimate and drivers of streamflow within the Upper Missouri River basin. This study presents the first annually resolved tree‐ring reconstruction of spring precipitation for the Missouri River Headwaters. This reconstruction along with existing tree‐ring reconstructions of 1 April snow water equivalence (SWE) and water year streamflow are used to detect variable influences of winter and spring precipitation on streamflow over past centuries, and relative to the modern period. The results suggest that spring precipitation has been a more consistent influence on water year streamflow in the Missouri River Headwaters over past centuries than winter snowpack. The strong relationship between 1 April SWE values and water year streamflow in the Missouri River Headwaters observed over much of the 20th century is not found to be a consistent feature of these multicentury paleorecords. Instead, the climatic influences on streamflow within the Missouri River Headwaters are likely more variable than 20th‐century instrumental records indicate.

     
    more » « less
  5. Abstract. A critical component of hydrologic modeling in cold andtemperate regions is partitioning precipitation into snow and rain, yetlittle is known about how uncertainty in precipitation phase propagates intovariability in simulated snow accumulation and melt. Given the wide varietyof methods for distinguishing between snow and rain, it is imperative toevaluate the sensitivity of snowpack model output to precipitation phasedetermination methods, especially considering the potential of snow-to-rainshifts associated with climate warming to fundamentally change the hydrologyof snow-dominated areas. To address these needs we quantified thesensitivity of simulated snow accumulation and melt to rain–snowpartitioning methods at sites in the western United States using theSNOWPACK model without the canopy module activated. The methods in thisstudy included different permutations of air, wet bulb and dew pointtemperature thresholds, air temperature ranges, and binary logisticregression models. Compared to observations of snow depth and snow water equivalent (SWE), thebinary logistic regression models produced the lowest mean biases, whilehigh and low air temperature thresholds tended to overpredict andunderpredict snow accumulation, respectively. Relative differences betweenthe minimum and maximum annual snowfall fractions predicted by the differentmethods sometimes exceeded 100 % at elevations less than 2000 m in theOregon Cascades and California's Sierra Nevada. This led to rangesin annual peak SWE typically greater than 200 mm,exceeding 400 mm in certain years. At the warmer sites, ranges in snowmelttiming predicted by the different methods were generally larger than 2 weeks, while ranges in snow cover duration approached 1 month and greater.Conversely, the three coldest sites in this work were relatively insensitiveto the choice of a precipitation phase method, with average ranges in annualsnowfall fraction, peak SWE, snowmelt timing, and snow cover duration of lessthan 18 %, 62 mm, 10 d, and 15 d, respectively. Average ranges in snowmeltrate were typically less than 4 mm d−1 and exhibited a smallrelationship to seasonal climate. Overall, sites with a greater proportionof precipitation falling at air temperatures between 0 and4 ∘C exhibited the greatest sensitivity to method selection,suggesting that the identification and use of an optimal precipitation phasemethod is most important at the warmer fringes of the seasonal snow zone. 
    more » « less