Abstract Extreme floods and landslides in high‐latitude watersheds have been associated with rain‐on‐snow (ROS) events. Yet, the risks of changing precipitation phases on a declining snowpack under a warming climate remain unclear. Normalizing the total annual duration of ROS with that of the seasonal snowpack, the ERA5 data (1941–2023) show that the frequency of high‐runoff ROS events is a characteristic feature of high‐latitude coastal zones, particularly over the coasts of south‐central Alaska and southern Newfoundland. Total rainfall accumulation per seasonal snowpack duration has increased across western mountain ranges, with the Olympic Mountains experiencing more than 40 mm of additional rainfall over the snowpack in the past eight decades, followed by the Sierra Nevada. These trends could drive an 8% increase in rainfall extremes (e.g., more than 10 mm for 6 hr storm with a 15‐year return period), highlighting the need for resilient flood control systems in high‐latitude coastal cities. 
                        more » 
                        « less   
                    
                            
                            The sensitivity of modeled snow accumulation and melt to precipitation phase methods across a climatic gradient
                        
                    
    
            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   
        
    
                            - Award ID(s):
- 1637686
- PAR ID:
- 10126821
- Date Published:
- Journal Name:
- Hydrology and Earth System Sciences
- Volume:
- 23
- Issue:
- 9
- ISSN:
- 1607-7938
- Page Range / eLocation ID:
- 3765 to 3786
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
- 
            
- 
            Intermittent snow depth observations can be leveraged with data assimilation (DA) to improve model estimates of snow water equivalent (SWE) at the point scale. A key consideration for scaling a DA system to the basin scale is its performance at locations with forest cover – where canopy-snow interactions affect snow accumulation and melt, yet are difficult to model and parameterize. We implement a particle filter (PF) assimilation technique to assimilate intermittent depth observations into the Flexible Snow Model (FSM2), and validate the output against snow density and SWE measurements across paired forest and open sites, at two locations with different climates and forest structures. Assimilation reduces depth error by 70-90%, density error by 5-30%, and SWE error by 50-70% at forest locations (relative to control model runs) and brings errors in-line with adjacent open sites. The PF correctly simulates the seasonal evolution of the snowpack under forest canopy, including cases where interception lowers SWE in the forest during accumulation, and shading reduces melt during the ablation season (relative to open sites). The snow model outputs are sensitive to canopy-related parameters, but DA reduces the range in depth and SWE estimates resulting from spatial variations or uncertainties in these parameters by more than 50%. The results demonstrate that the challenge of accurately measuring, estimating, or calibrating canopy-related parameters is reduced when snow depth observations are assimilated to improve SWE estimates.</p>more » « less
- 
            Abstract. Data from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition allowed us to investigate the temporal dynamics of snowfall, snow accumulation and erosion in great detail for almost the whole accumulation season (November 2019 to May 2020). We computed cumulative snow water equivalent (SWE) over the sea ice based on snow depth and density retrievals from a SnowMicroPen and approximately weekly measured snow depths along fixed transect paths. We used the derived SWE from the snow cover to compare with precipitation sensors installed during MOSAiC. The data were also compared with ERA5 reanalysis snowfall rates for the drift track. We found an accumulated snow mass of 38 mm SWE between the end of October 2019 and end of April 2020. The initial SWE over first-year ice relative to second-year ice increased from 50 % to 90 % by end of the investigation period. Further, we found that the Vaisala Present Weather Detector 22, an optical precipitation sensor, and installed on a railing on the top deck of research vessel Polarstern, was least affected by blowing snow and showed good agreements with SWE retrievals along the transect. On the contrary, the OTT Pluvio2 pluviometer and the OTT Parsivel2 laser disdrometer were largely affected by wind and blowing snow, leading to too high measured precipitation rates. These are largely reduced when eliminating drifting snow periods in the comparison. ERA5 reveals good timing of the snowfall events and good agreement with ground measurements with an overestimation tendency. Retrieved snowfall from the ship-based Ka-band ARM zenith radar shows good agreements with SWE of the snow cover and differences comparable to those of ERA5. Based on the results, we suggest the Ka-band radar-derived snowfall as an upper limit and the present weather detector on RV Polarstern as a lower limit of a cumulative snowfall range. Based on these findings, we suggest a cumulative snowfall of 72 to 107 mm and a precipitation mass loss of the snow cover due to erosion and sublimation as between 47 % and 68 %, for the time period between 31 October 2019 and 26 April 2020. Extending this period beyond available snow cover measurements, we suggest a cumulative snowfall of 98–114 mm.more » « less
- 
            null (Ed.)Abstract When compared with differences in snow accumulation predicted by widely used hydrological models, there is a much greater divergence among otherwise “good” models in their simulation of the snow ablation process. Here, we explore differences in the performance of the Variable Infiltration Capacity model (VIC), Noah land surface model with multiparameterization options (Noah-MP), the Catchment model, and the third-generation Simplified Simple Biosphere model (SiB3) in their ability to reproduce observed snow water equivalent (SWE) during the ablation season at 10 Snowpack Telemetry (SNOTEL) stations over 1992–2012. During the ablation period, net radiation generally has stronger correlations with observed melt rates than does air temperature. Average ablation rates tend to be higher (in both model predictions and observations) at stations with a large accumulation of SWE. The differences in the dates of last snow between models and observations range from several days to approximately a month (on average 5.1 days earlier than in observations). If the surface cover in the models is changed from observed vegetation to bare soil in all of the models, only the melt rate of the VIC model increases. The differences in responses of models to canopy removal are directly related to snowpack energy inputs, which are further affected by different algorithms for surface albedo and energy allocation across the models. We also find that the melt rates become higher in VIC and lower in Noah-MP if the shrub/grass present at the observation sites is switched to trees.more » « less
- 
            Abstract In the western United States, water supplies largely originate as snowmelt from forested land. Forests impact the water balance of these headwater streams, yet most predictive runoff models do not explicitly account for changing snow‐vegetation dynamics. Here, we present a case study showing how warmer temperatures and changing forests in the Henrys Fork of the Snake River, a seasonally snow‐covered headwater basin in the Greater Yellowstone Ecosystem, have altered the relationship between April 1st snow water equivalent (SWE) and summer streamflow. Since the onset and recovery of severe drought in the early 2000s, predictive models based on pre‐drought relationships over‐predict summer runoff in all three headwater tributaries of the Henrys Fork, despite minimal changes in precipitation or snow accumulation. Compared with the pre‐drought period, late springs and summers (May–September) are warmer and vegetation is greener with denser forests due to recovery from multiple historical disturbances. Shifts in the alignment of snowmelt and energy availability due to warmer temperatures may reduce runoff efficiency by changing the amount of precipitation that goes to evapotranspiration versus runoff and recharge. To quantify the alignment between snowmelt and energy on a timeframe needed for predictive models, we propose a new metric, the Vegetation‐Water Alignment Index (VWA), to characterize the synchrony of vegetation greenness and snowmelt and rain inputs. New predictive models show that in addition to April 1st SWE, the previous year's VWA and summer reference evapotranspiration are the most significant predictors of runoff in each watershed and provide more predictive power than traditionally used metrics. These results suggest that the timing of snowmelt relative to the start of the growing season affects not only annual partitioning of streamflow, but can also determine the groundwater storage state that dictates runoff efficiency the following spring.more » « less
 An official website of the United States government
An official website of the United States government 
				
			 
					 
					
 
                                    