Observed Changes in Climate and Streamflow in the Upper Rio Grande Basin

Observed streamflow and climate data are used to test the hypothesis that climate change is already affecting Rio Grande streamflow volume derived from snowmelt runoff in ways consistent with model‐based projections of 21st‐Century streamflow. Annual and monthly changes in streamflow volume and surface climate variables on the Upper Rio Grande, near its headwaters in southern Colorado, are assessed for water years 1958–2015. Results indicate winter and spring season temperatures in the basin have increased significantly, April 1 snow water equivalent (SWE) has decreased by approximately 25%, and streamflow has declined slightly in the April–July snowmelt runoff season. Small increases in precipitation have reduced the impact of declining snowpack on trends in streamflow. Changes in the snowpack–runoff relationship are noticeable in hydrographs of mean monthly streamflow, but are most apparent in the changing ratios of precipitation (rain + snow, and SWE) to streamflow and in the declining fraction of runoff attributable to snowpack or winter precipitation. The observed changes provide observational confirmation for model projections of decreasing runoff attributable to snowpack, and demonstrate the decreasing utility of snowpack for predicting subsequent streamflow on a seasonal basis in the Upper Rio Grande Basin.     

FOREST GROWTH IMPACT ON PEAK SNOW PACK MEASUREMENTS AT LONG TERM SNOW COURSES1

ABSTRACT: Snow course surveys in late winter provide stream‐flow forecasters with their best information for making water supply and flood forecasts for the subsequent spring and summer runoff period in mountainous regions of western North America. Snow survey data analyses are generally based on a 30‐year “normal” period. It is well documented that forest cover changes over time will affect snow accumulation on the ground within forests. This paper seeks to determine if forest cover changes over decades at long term snow courses decrease measured peak snow water equivalent (SWE) enough to affect runoff prediction. Annual peak SWE records were analyzed at four snow courses in two different forest types having at least 25 years of snowpack data to detect any decreases in SWE due to forest growth. No statistically significant decreases in annual peak SWE over time were found at any of these four snow courses. The wide range of annual winter precipitation and correspondingly highly variable peak snowpack accumulation, as well as many other weather and site variables, masked any minor trends in the data.     

EFFECTS OF CLIMATE CHANGE ON HYDROLOGY AND WATER RESOURCES IN THE COLUMBIA RIVER BASIN1

ABSTRACT: As part of the National Assessment of Climate Change, the implications of future climate predictions derived from four global climate models (GCMs) were used to evaluate possible future changes to Pacific Northwest climate, the surface water response of the Columbia River basin, and the ability of the Columbia River reservoir system to meet regional water resources objectives. Two representative GCM simulations from the Hadley Centre (HC) and Max Planck Institute (MPI) were selected from a group of GCM simulations made available via the National Assessment for climate change. From these simulations, quasi-stationary, decadal mean temperature and precipitation changes were used to perturb historical records of precipitation and temperature data to create inferred conditions for 2025, 2045, and 2095. These perturbed records, which represent future climate in the experiments, were used to drive a macro-scale hydrology model of the Columbia River at 1/8 degree resolution. The altered streamflows simulated for each scenario were, in turn, used to drive a reservoir model, from which the ability of the system to meet water resources objectives was determined relative to a simulated hydrologic base case (current climate). Although the two GCM simulations showed somewhat different seasonal patterns for temperature change, in general the simulations show reasonably consistent basin average increases in temperature of about 1.8–2.1°C for 2025, and about 2.3–2.9°C for 2045. The HC simulations predict an annual average temperature increase of about 4.5°C for 2095. Changes in basin averaged winter precipitation range from -1 percent to + 20 percent for the HC and MPI scenarios, and summer precipitation is also variously affected. These changes in climate result in significant increases in winter runoff volumes due to increased winter precipitation and warmer winter temperatures, with resulting reductions in snowpack. Average March 1 basin average snow water equivalents are 75 to 85 percent of the base case for 2025, and 55 to 65 percent of the base case by 2045. By 2045 the reduced snowpack and earlier snow melt, coupled with higher evapotranspiration in early summer, would lead to earlier spring peak flows and reduced runoff volumes from April-September ranging from about 75 percent to 90 percent of the base case. Annual runoff volumes range from 85 percent to 110 percent of the base case in the simulations for 2045. These changes in streamflow create increased competition for water during the spring, summer, and early fall between non-firm energy production, irrigation, instream flow, and recreation. Flood control effectiveness is moderately reduced for most of the scenarios examined, and desirable navigation conditions on the Snake are generally enhanced or unchanged. Current levels of winter-dominated firm energy production are only significantly impacted for the MPI 2045 simulations.     

Use of observed hydroclimatic trends to constrain projections of snowmelt season runoff in the Rio Grande headwaters

A method is developed for choosing 21st Century streamflow projections among widely varying results from a large ensemble of climate model-driven simulations. We quantify observed trends in climate–streamflow relationships in the Rio Grande headwaters, which has experienced warming temperature and declining snowpack since the mid-20th Century. Prominent trends in the snowmelt runoff season are used to assess corresponding statistics in downscaled global climate model projections. We define “Observationally Consistent (OC)” simulations as those that reproduce historical changes to linear statistics of diminished snowpack–streamflow coupling in the headwaters and an associated increase in the contribution of spring season (post-peak snowpack) precipitation to streamflow. Only a modest fraction of the ensemble of simulations meets these consistency metrics. The subset of OC simulations projects significant decreases in headwaters flow, whereas the simulations that poorly replicate historical trends exhibit a much wider range of projected changes. These results bolster confidence in model-based projections of declining runoff in the Rio Grande headwaters in the snowmelt runoff season and offer an example of a methodology for evaluating model-based projections in basins with similar hydroclimates that have experienced pronounced climate changes in the recent historical record.     

Defining the Diurnal Pattern of Snowmelt Using a Beta Distribution Function

Snow is an important component of the hydrologic cycle for many regions worldwide. In addition to vital water resources, snowmelt can be important for forest ecosystem dynamics and flood risk. However, standard design events in the United States lack a design snowmelt event, including only precipitation events, though snowmelt has been shown to be larger than rainfall. In this article, we present a method using hourly snow water equivalent data to develop and test a function for representing the diurnal pattern of snowmelt. A two‐parameter beta distribution function is modified for the purposes of this study and found to fit the pattern of snowmelt well with a root mean squared error of 0.008. Soil moisture sensors were additionally utilized to assess the timing of the snowmelt water outflow from the base of the snowpack that supports the shape of the function, but suggests that the timing of losses recorded on snow pillows lag as much as 3 h. Further testing of the function showed the shape of the function to be accurate. The methods developed and tested in this paper can be applied for design purposes comparing snowmelt and rainfall events or to improve hydrological models investigating processes such as streamflow or groundwater recharge.     

Hydrologic Calibration and Validation of SWAT in a Snow‐Dominated Rocky Mountain Watershed,Montana, U.S.A.1

Abstract: The Soil and Water Assessment Tool (SWAT) has been applied successfully in temperate environments but little is known about its performance in the snow‐dominated, forested, mountainous watersheds that provide much of the water supply in western North America. To address this knowledge gap, we configured SWAT to simulate the streamflow of Tenderfoot Creek (TCSWAT). Located in central Montana, TCSWAT represents a high‐elevation watershed with ～85% coniferous forest cover where more than 70% of the annual precipitation falls as snow, and runoff comes primarily from spring snowmelt. Model calibration using four years of measured daily streamflow, temperature, and precipitation data resulted in a relative error (RE) of 2% for annual water yield estimates, and mean paired deviations (Dv) of 36 and 31% and Nash‐Sutcliffe (NS) efficiencies of 0.90 and 0.86 for monthly and daily streamflow, respectively. Model validation was conducted using an additional four years of data and the performance was similar to the calibration period, with RE of 4% for annual water yields, Dv of 43% and 32%, and NS efficiencies of 0.90 and 0.76 for monthly and daily streamflow, respectively. An objective, regression‐based model invalidation procedure also indicated that the model was validated for the overall simulation period. Seasonally, SWAT performed well during the spring and early summer snowmelt runoff period, but was a poor predictor of late summer and winter base flow. The calibrated model was most sensitive to snowmelt parameters, followed in decreasing order of influence by the surface runoff lag, ground water, soil, and SCS Curve Number parameter sets. Model sensitivity to the surface runoff lag parameter reflected the influence of frozen soils on runoff processes. Results indicated that SWAT can provide reasonable predictions of annual, monthly, and daily streamflow from forested montane watersheds, but further model refinements could improve representation of snowmelt runoff processes and performance during the base flow period in this environment.     

Trends in Western U.S. Snowpack and Related Upper Colorado River Basin Streamflow1

Miller, W. Paul and Thomas C. Piechota, 2011. Trends in Western U.S. Snowpack and Related Upper Colorado River Basin Streamflow. Journal of the American Water Resources Association (JAWRA) 47(6):1197–1210. DOI: 10.1111/j.1752‐1688.2011.00565.x Abstract: Water resource managers in the Western United States (U.S.) are currently faced with the challenge of adapting to unprecedented drought and uncertain impacts of climate change. Recent research has indicated increasing regional temperature and changes to precipitation and streamflow characteristics throughout the Western U.S. As such, there is increased uncertainty in hydroclimatological forecasts, which impact reservoir operations and water availability throughout the Western U.S., particularly in the Colorado River Basin. Previous research by the authors hypothesized a change in the character of precipitation (i.e., the frequency and amount of rainfall and snowfall events) throughout the Colorado River Basin. In the current study, 398 snowpack telemetry stations were investigated for trends in cumulative precipitation, snow water equivalent, and precipitation events. Observations of snow water equivalent characteristics were compared to observations in streamflow characteristics. Results indicate that the timing of the last day of the snow season corresponds well to the volume of runoff observed over the traditional peak flow season (April through July); conversely, the timing of the first day of the snow season does not correspond well to the volume of runoff observed over the peak flow season. This is significant to water resource managers and river forecasters, as snowpack characteristics may be indicative of a productive or unproductive runoff season.     

Spatiotemporal Variability of Snow Depletion Curves Derived from SNODAS for the Conterminous United States, 2004‐2013

Assessment of water resources at a national scale is critical for understanding their vulnerability to future change in policy and climate. Representation of the spatiotemporal variability in snowmelt processes in continental‐scale hydrologic models is critical for assessment of water resource response to continued climate change. Continental‐extent hydrologic models such as the U.S. Geological Survey National Hydrologic Model (NHM) represent snowmelt processes through the application of snow depletion curves (SDCs). SDCs relate normalized snow water equivalent (SWE) to normalized snow covered area (SCA) over a snowmelt season for a given modeling unit. SDCs were derived using output from the operational Snow Data Assimilation System (SNODAS) snow model as daily 1‐km gridded SWE over the conterminous United States. Daily SNODAS output were aggregated to a predefined watershed‐scale geospatial fabric and used to also calculate SCA from October 1, 2004 to September 30, 2013. The spatiotemporal variability in SNODAS output at the watershed scale was evaluated through the spatial distribution of the median and standard deviation for the time period. Representative SDCs for each watershed‐scale modeling unit over the conterminous United States (n = 54,104) were selected using a consistent methodology and used to create categories of snowmelt based on SDC shape. The relation of SDC categories to the topographic and climatic variables allow for national‐scale categorization of snowmelt processes.     

Projections of 21st Century Sierra Nevada Local Hydrologic Flow Components Using an Ensemble of General Circulation Models1

Abstract: Sierra Nevada snowmelt and runoff is a key source of water for many of California’s 38 million residents and nearly the entire population of western Nevada. The purpose of this study was to assess the impacts of expected 21st Century climatic changes in the Sierra Nevada at the subwatershed scale, for all hydrologic flow components, and for a suite of 16 General Circulation Models (GCMs) with two emission scenarios. The Soil and Water Assessment Tool (SWAT) was calibrated and validated at 35 unimpaired streamflow sites. Results show that temperatures are projected to increase throughout the Sierra Nevada, whereas precipitation projections vary between GCMs. These climatic changes drive a decrease in average annual streamflow and an advance of snowmelt and runoff by several weeks. The largest streamflow reductions were found in the mid‐range elevations due to less snow accumulation, whereas the higher elevation watersheds were more resilient due to colder temperatures. Simulation results showed that decreases in snowmelt affects not only streamflow, but evapotranspiration, surface, and subsurface flows, such that less water is available in spring and summer, thus potentially affecting aquatic and terrestrial ecosystems. Declining spring and summer flows did not equally affect all subwatersheds in the region, and the subwatershed perspective allowed for identification for the most sensitive basins throughout the Sierra Nevada.     

GENERAL-CIRCULATION-MODEL SIMULATIONS OF FUTURE SNOWPACK IN THE WESTERN UNITED STATES1

ABSTRACT: April 1 snowpack accumulations measured at 311 snow courses in the western United States (U.S.) are grouped using a correlation-based cluster analysis. A conceptual snow accumulation and melt model and monthly temperature and precipitation for each cluster are used to estimate cluster-average April 1 snowpack. The conceptual snow model is subsequently used to estimate future snowpack by using changes in monthly temperature and precipitation simulated by the Canadian Centre for Climate Modeling and Analysis (CCC) and the Hadley Centre for Climate Prediction and Research (HADLEY) general circulation models (GCMs). Results for the CCC model indicate that although winter precipitation is estimated to increase in the future, increases in temperatures will result in large decreases in April 1 snowpack for the entire western U.S. Results for the HADLEY model also indicate large decreases in April 1 snowpack for most of the western US, but the decreases are not as severe as those estimated using the CCC simulations. Although snowpack conditions are estimated to decrease for most areas of the western US, both GCMs estimate a general increase in winter precipitation toward the latter half of the next century. Thus, water quantity may be increased in the western US; however, the timing of runoff will be altered because precipitation will more frequently occur as rain rather than as snow.     

REVISITING THE DEGREE-DAY METHOD FOR SNOWMELT COMPUTATIONS1

ABSTRACT: The simple, empirical degree-day approach for calculating snowmelt and runoff from mountain basins has been in use for more than 60 years. It is frequently suggested that the degree-day method be replaced by the more physically-based energy balance approach. The degree-day approach, however, maintains its popularity, applicability, and effectiveness. It is shown that the degree-day method is reliable for computing total snowmelt depths for periods of a week to the entire snowmelt season. It can also be used for daily snowmelt depths when utilized in connection with an adequate snowmelt runoff model for computing the basin runoff. The degree-day ratio is shown to vary seasonally as opposed to being constant as is often assumed. Additionally, in order to evaluate the degree-day ratio correctly, the changing snow cover extent in a basin during the snowmelt season must be taken into account. It is also possible to combine the degree-day approach with a radiation component so that short time interval (<24 hours) computations of snowmelt depth can be made. When snowmelt input is transformed to basin output (runoff) by a snowmelt runoff model, there is little difference between the degree-day approach and a radiation-based approach. This is fortuitous because the physically-based energy balance models will not soon displace the degree-day methods because of their excessive data requirements.     

Management Implications of Snowpack Sensitivity to Temperature and Atmospheric Moisture Changes in Yosemite National Park,CA

In order to investigate snowpack sensitivity to temperature increases and end‐member atmospheric moisture conditions, we applied a well‐constrained energy‐ and mass‐balance snow model across the full elevation range of seasonal snowpack using forcing data from recent wet and dry years. Humidity scenarios examined were constant relative humidity (high) and constant vapor pressure between storms (low). With minimum calibration, model results captured the observed magnitude and timing of snowmelt. April 1 snow water equivalent (SWE) losses of 38%, 73%, and 90% with temperature increases of 2, 4, and 6°C in a dry year centered on areas of greatest SWE accumulation. Each 2°C increment of warming also resulted in seasonal snowline moving upslope by 300 m. The zone of maximum melt was compressed upward 100–500 m with 6°C warming, with the range reflecting differences in basin hypsometry. Melt contribution by elevations below 2,000 m disappeared with 4°C warming. The constant‐relative‐humidity scenario resulted in 0–100 mm less snowpack in late spring vs. the constant‐vapor‐pressure scenario in a wet year, a difference driven by increased thermal radiation (+1.2 W/m²) and turbulent energy fluxes (+1.2 W/m²) to the snowpack for the constant‐relative‐humidity case. Loss of snowpack storage and potential increases in forest evapotranspiration due to warming will result in a substantial shift in forest water balance and present major challenges to land management in this mountainous region.     

Snowmelt Runoff and Water Yield Along Elevation and Temperature Gradients in California’s Southern Sierra Nevada1

Hunsaker, Carolyn T., Thomas W. Whitaker, and Roger C. Bales, 2012. Snowmelt Runoff and Water Yield Along Elevation and Temperature Gradients in California’s Southern Sierra Nevada. Journal of the American Water Resources Association (JAWRA) 48(4): 667‐678. DOI: 10.1111/j.1752‐1688.2012.00641.x Abstract: Differences in hydrologic response across the rain‐snow transition in the southern Sierra Nevada were studied in eight headwater catchments – the Kings River Experimental Watersheds – using continuous precipitation, snowpack, and streamflow measurements. The annual runoff ratio (discharge divided by precipitation) increased about 0.1 per 300 m of mean catchment elevation over the range 1,800‐2,400 m. Higher‐elevation catchments have lower vegetation density, shallow soils with rapid permeability, and a shorter growing season when compared with those at lower elevations. Average annual temperatures ranged from 6.8°C at 2,400 m to 8.6 at 1,950 m elevation, with annual precipitation being 75‐95% snow at the highest elevations vs. 20‐50% at the lowest. Peak discharge lagged peak snow accumulation on the order of 60 days at the higher elevations and 20 to 30 days at the lower elevations. Snowmelt dominated the daily streamflow cycle over a period of about 30 days in higher elevation catchments, followed by a 15‐day transition to evapotranspiration dominating the daily streamflow cycle. Discharge from lower elevation catchments was rainfall dominated in spring, with the transition to evapotranspiration dominance being less distinct. Climate warming that results in a longer growing season and a shift from snow to rain would result in earlier runoff and a lower runoff ratio.     

EFFECT OF SIMULATED CLIMATE CHANGE ON SNOWMELT RUNOFF MODELING IN SELECTED BASINS1

ABSTRACT: The projected increase in the concentration of CO₂ and other greenhouse gases in the atmosphere is likely to result in a global temperature increase. This paper reports on the probable effects of a temperature increase and changes in transpiration on basin discharge in two different mountain snowmelt regions of the western United States. The hydrological effects of the climate changes are modeled with a relatively simple conceptual, semi-distributed snowmelt runoff model. Based on the model results, it may be concluded that increased air temperatures will result in a shift of snowmelt runoff to earlier in the snowmelt season. Furthermore, it is shown that it is very important to include the expected change in climate-related basin conditions resulting from the modeled temperature increase in the runoff simulation. The effect of adapting the model parameters to reflect the changed basin conditions resulted in a further shift of streamflow to April and an even more significant decrease of snowmelt runoff in June and July. If the air temperatures increase by approximately 5°C and precipitation and accumulated snow amounts remain about the same, runoff in April and May, averaged for the two basins, is expected to increase by 185 percent and 26 percent, respectively. The runoff in June and July will decrease by about 60 percent each month. Overall, the total seasonal runoff decreases by about 6 percent. If increased CO₂ concentrations further change basin conditions by reducing transpiration by the maximum amounts reported in the literature, then, combined with the 5°C temperature increase, the April, May, June, and July changes would average +230 percent, +40 percent, ?55 percent, and ?45 percent, respectively. The total seasonal runoff change would be +11 percent.     

VERIFICATION OF THE RHEA-OROGRAPHIC-PRECIPITATION MODEL1

ABSTRACT: Observed April 1 snowpack accumulations within and near the Gunnison River basin in southwestern Colorado are compared with simulations from the Rhea-orographic-precipitation model to determine if the model simulates reliable magnitudes and temporal and spatial variability in winter precipitation for the basin. Twenty simulations of the Rhea model were performed using‘optimal’parameter sets determined for 10-kilometer (km) grids (10-km by 10-km grid cells) through stochastic calibration. Comparisons of Rhea-model simulations of winter precipitation with April 1 snowpack accumulations at 32 snowcourse stations were performed for the years 1972–1990. For most stations and most years the Rhea model reliably simulates the temporal and spatial variability in April 1 snowpack accumulations. However, in general, the Rhea-model underestimates April 1 snowpack accumulations in the Gunnison River basin area, and the underestimation is greatest for locations that receive the largest amount of snow. A significant portion of the error in Rhea-model simulations is due to the calibration of the Rhea model using gauge-catch precipitation measurements which can be as much as 50 percent below actual snowfall accumulations. Additional error in the Rhea-model simulations is a result of the comparison of gridded precipitation values to observed values measured at points.     

POTENTIAL IMPACTS OF CLIMATE CHANGE ON CALIFORNIA HYDROLOGY1

ABSTRACT: Previous reports based on climate change scenarios have suggested that California will be subjected to increased wintertime and decreased summertime streamflow. Due to the uncertainty of projections in future climate, a new range of potential climatological future temperature shifts and precipitation ratios is applied to the Sacramento Soil Moisture Accounting Model and Anderson Snow Model in order to determine hydrologic sensitivities. Two general circulation models (GCMs) were used in this analysis: one that is warm and wet (HadCM2 run 1) and one that is cool and dry (PCM run B06.06), relative to the GCM projections for California that were part of the Third Assessment Report of the Intergovernmental Panel on Climate Change. A set of specified incremental temperature shifts from 1.5°C to 5.0°C and precipitation ratios from 0.70 to 1.30 were also used as input to the snow and soil moisture accounting models, providing for additional scenarios (e.g., warm/dry, cool/wet). Hydrologic calculations were performed for a set of California river basins that extend from the coastal mountains and Sierra Nevada northern region to the southern Sierra Nevada region; these were applied to a water allocation analysis in a companion paper. Results indicate that for all snow‐producing cases, a larger proportion of the streamflow volume will occur earlier in the year. The amount and timing is dependent on the characteristics of each basin, particularly the elevation. Increased temperatures lead to a higher freezing line, therefore less snow accumulation and increased melting below the freezing height. The hydrologic response varies for each scenario, and the resulting solution set provides bounds to the range of possible change in streamflow, snowmelt, snow water equivalent, and the change in the magnitude of annual high flows. An important result that appears for all snowmelt driven runoff basins, is that late winter snow accumulation decreases by 50 percent toward the end of this century.     

Modeling Discharge Rates Using a Coupled Modeled Approach for the Merced River in Yosemite National Park

This study describes the application of the NASA version of the Carnegie‐Ames‐Stanford Approach (CASA) ecosystem model coupled with a surface hydrologic routing scheme previously called the Hydrological Routing Algorithm (HYDRA) to model monthly discharge rates from 2000 to 2007 on the Merced River drainage in Yosemite National Park, California. To assess CASA‐HYDRA's capability to estimate actual water flows in extreme precipitation years, the focus of this study is the 2007 water year, which was very dry, and the 2005 water year, which was a moderately wet year in the historical record. Prior to comparisons to gauge records, CASA‐HYDRA snowmelt algorithms were modified with equations from the U.S. Department of Agriculture Snowmelt‐Runoff Model (SRM), which has been designed to predict daily streamflow in mountain basins where snowmelt is a major runoff factor. Results show that model predictions closely matched monthly flow rates at the Pohono Bridge gauge station (USGS#11266500), with R² = 0.67 and Nash‐Sutcliffe (E) = 0.65. By subdividing the upper Merced River basin into subbasins with high spatial resolution in the gridded modeling approach, we were able to determine which biophysical characteristics in the Sierra differed to the largest degree in extreme low‐flow and high‐flow years. Average elevation and snowpack accumulation were found to be the most important explanatory variables to understand subbasin contributions to monthly discharge rates.     

Impacts of 21st‐Century Climate Change on Hydrologic Extremes in the Pacific Northwest Region of North America

Climate change projections for the Pacific Northwest (PNW) region of North America include warmer temperatures (T), reduced precipitation (P) in summer months, and increased P during all other seasons. Using a physically based hydrologic model and an ensemble of statistically downscaled global climate model scenarios produced by the Columbia Basin Climate Change Scenarios Project, we examine the nature of changing hydrologic extremes (floods and low flows) under natural conditions for about 300 river locations in the PNW. The combination of warming, and shifts in seasonal P regimes, results in increased flooding and more intense low flows for most of the basins in the PNW. Flood responses depend on average midwinter T and basin type. Mixed rain and snow basins, with average winter temperatures near freezing, typically show the largest increases in flood risk because of the combined effects of warming (increasing contributing basin area) and more winter P. Decreases in low flows are driven by loss of snowpack, drier summers, and increasing evapotranspiration in the simulations. Energy‐limited basins on the west side of the Cascades show the strongest declines in low flows, whereas more arid, water‐limited basins on the east side of the Cascades show smaller reductions in low flows. A fine‐scale analysis of hydrologic extremes over the Olympic Peninsula echoes the results for the larger rivers discussed above, but provides additional detail about topographic gradients.     

Winter Hydrology and NO3− Concentrations in a Forested Watershed: A Detailed Field Study in the Adirondack Mountains of New York1

Abstract: More than 85% of NO₃^? losses from watersheds in the northeastern United States are exported during winter months (October 1 to May 30). Interannual variability in NO₃^? loads to individual streams is closely related to interannual climatic variations, particularly during the winter. The objective of our study was to understand how climatic and hydrogeological factors influence NO₃^? dynamics in small watersheds during the winter. Physical parameters including snow depth, soil temperature, stream discharge, and water table elevation were monitored during the 2007‐2008 winter in two small catchments in the Adirondack Mountains, New York State. Snowpack persisted from mid‐December to mid‐April, insulating soils such that only two isolated instances of soil frost were observed during the study period. NO₃^? export during a mid‐winter rain‐on‐snowmelt event comprised between 8 and 16% of the total stream NO₃^? load for the four‐month winter study period. This can be compared with the NO₃^? exported during the final spring melt, which comprised between 38 and 45% of the total four‐month winter NO₃^? load. Our findings indicate that minor melt events were detectable with changes in soil temperature, streamflow, groundwater level, and snow depth. But, based on loading, these events were relatively minor contributors to winter NO₃^? loss. A warmer climate and fluctuating snowpack may result in more major mid‐winter melt events and greater NO₃^? export to surface waters.     

Influence of manure application on surface energy and snow cover: model development and sensitivities

Winter landspreading is an important part of manure management in the U.S. Upper Midwest. Although the practice is thought to lead to excessive P runoff losses, surprisingly little has been learned from field experiments or current water quality models. We captured knowledge gained from winter manure landspreading experiments by modifying a mechanistic snow ablation model to include manure. The physically based, modified model simulated the observed delay in snow cover disappearance and surface energy balance changes caused by application of the manure. Additional model simulations of surface energy balance estimates of radiation and turbulent fluxes showed that during intense melting events the manure on top of snow significantly reduced the energy available for melt of the snow underneath, slowing melt. The effect was most pronounced when snowmelt was driven by both relatively high solar radiation and turbulent heat fluxes. High absorbed shortwave radiation caused significant warming of the manure, which led to substantial losses in turbulent fluxes and longwave radiation. Simulations of snowmelt also showed that manure applications between 45 and 100 Mg ha(-1) significantly reduced peak snowmelt rates, in proportion to the manure applied. Lower snowmelt rates beneath manure may allow more infiltration of meltwater compared with bare snow. This infiltration and attenuated snowmelt runoff may partially mitigate the enhanced likelihood of P runoff from unincorporated winter-spread manure.