Modeling Streamflow and Water Quality Sensitivity to Climate Change and Urban Development in 20 U.S. Watersheds

Watershed modeling in 20 large, United States (U.S.) watersheds addresses gaps in our knowledge of streamflow, nutrient (nitrogen and phosphorus), and sediment loading sensitivity to mid‐21st Century climate change and urban/residential development scenarios. Use of a consistent methodology facilitates regional scale comparisons across the study watersheds. Simulations use the Soil and Water Assessment Tool. Climate change scenarios are from the North American Regional Climate Change Assessment Program dynamically downscaled climate model output. Urban and residential development scenarios are from U.S. Environmental Protection Agency's Integrated Climate and Land Use Scenarios project. Simulations provide a plausible set of streamflow and water quality responses to mid‐21st Century climate change across the U.S. Simulated changes show a general pattern of decreasing streamflow volume in the central Rockies and Southwest, and increases on the East Coast and Northern Plains. Changes in pollutant loads follow a similar pattern but with increased variability. Ensemble mean results suggest that by the mid‐21st Century, statistically significant changes in streamflow and total suspended solids loads (relative to baseline conditions) are possible in roughly 30‐40% of study watersheds. These proportions increase to around 60% for total phosphorus and total nitrogen loads. Projected urban/residential development, and watershed responses to development, are small at the large spatial scale of modeling in this study.     

ALAWAT: A spatially allocated watershed model for approximating stream,sediment, and pollutant flows in Hawaii,USA

The Ala Wai Canal Watershed Model (ALAWAT) is a planning-level watershed model for approximating direct runoff, streamflow, sediment loads, and loads for up to five pollutants. ALAWAT uses raster GIS data layers including land use, SCS soil hydrologic groups, annual rainfall, and subwatershed delineations as direct model parameter inputs and can use daily total rainfall from up to ten rain gauges and streamflow from up to ten stream gauges. ALAWAT uses a daily time step and can simulate flows for up to ten-year periods and for up to 50 subwatersheds. Pollutant loads are approximated using a user-defined combination of rating curve relationships, mean event concentrations, and loading/washoff parameters for specific subwatersheds, land uses, and times of year. Using ALAWAT, annual average streamflow and baseflow relationships and urban suspended sediment loads were approximated for the Ala Wai Canal watershed (about 10,400 acres) on the island of Oahu, Hawaii. Annual average urban suspended sediments were approximated using two methods: mean event concentrations and pollutant loading and washoff. Parameters for the pollutant loading and washoff method were then modified to simulate the effect of various street sweeping intervals on sediment loads.     

Climatic and agricultural contributions to changing loads in two watersheds in Ohio

Trends in climatic variables, streamflow, agricultural practices, and loads of nutrients and suspended solids were estimated for 1976-1995 in the Maumee and Sandusky watersheds, two large agricultural basins draining to Lake Erie. To understand the contributions that various factors may have made to the trends in loads, earlier results of models linking loads to explanatory variables were combined with estimated trends in those variables. The study period was characterized by increases in temperature, wintertime precipitation and streamflow, conservation farming, and loads of nitrate and total suspended solids; decreases in snowfall and snow cover, fertilizer, manure from livestock, and loads of soluble reactive phosphorus; and relatively steady exports of total phosphorus. After removing the effects of trends in streamflow, nitrate loads increased much less while total suspended solids and total phosphorus loads declined. The analysis suggests that the nitrate increases were due largely to climatic factors, particularly increases in winter streamflow, decreases in snowfall and snow cover, and declining annual precipitation. Decreases in soluble reactive phosphorus were associated with changes in agricultural practices, particularly declines in fertilizer deliveries and head of livestock.     

THE EFFECTS OF CLIMATE CHANGE ON STREAM FLOW AND NUTRIENT LOADING'

ABSTRACT: This study assesses the potential impact of climate change on stream flow and nutrient loading in six watersheds of the Susquehanna River Basin using the Generalized Watershed Loading Function (GWLF). The model was used to simulate changes in stream flow and nutrient loads under a transient climate change scenario for each watershed. Under an assumption of no change in land cover and land management, the model was used to predict monthly changes in stream flow and nutrient loads for future climate conditions. Mean annual stream flow and nutrient loads increased for most watersheds, but decreased in one watershed that was intensively cultivated. Nutrient loading slightly decreased in April and late summer for several watersheds as a result of early snowmelt and increasing evapotranspiration. Spatial and temporal variability of stream flow and nutrient loads under the transient climate scenario indicates that different approaches for future water resource management may be useful.     

MODELING THE HYDROCHEMISTRY OF THE CANNONSVILLE WATERSHED WITH GENERALIZED WATERSHED LOADING FUNCTIONS (GWLF)1

ABSTRACT: A previous modeling study used the Generalized Watershed Loading Functions (GWLF) model to simulate stream‐flow, and nutrient and sediment loads to Cannonsville Reservoir from the West Branch Delaware River (WBDR). We made several model revisions, calibrated key parameters, and tested the original GWLF model and a revised GWLF model using more recent data. Model revisions included: addition of unsaturated leakage between unsaturated and saturated subsurface reservoirs; revised timing of sediment export; inclusion of urban sediments and dissolved nutrients; tracking of particulate nutrients from point sources; and revised timing of septic system loads. The revision of sediment yield timing resulted in significant improvements in monthly sediment and particulate phosphorus predictions as compared to the original model. Addition of unsaturated leakage improved hydrologic predictions during low flow months. The other model changes improve realism without adding significant model complexity or data requirements. Goodness of fit of revised model predictions versus stream measurements, as measured by the Nash‐Sutcliff coefficient of model efficiency, exceeded 0.8 for streamflow‐0.7 for sediment yield and dissolved nitrogen (N) and 0.6 for particulate and dissolved phosphorus (P). The revised GWLF model, with limited calibration, provides reasonable estimates of monthly streamflow, and nutrient and sediment loads in the Cannonsville watershed.     

Factors Affecting Stream Nutrient Loads: A Synthesis of Regional SPARROW Model Results for the Continental United States1

Preston, Stephen D., Richard B. Alexander, Gregory E. Schwarz, and Charles G. Crawford, 2011. Factors Affecting Stream Nutrient Loads: A Synthesis of Regional SPARROW Model Results for the Continental United States. Journal of the American Water Resources Association (JAWRA) 47(5):891‐915. DOI: 10.1111/j.1752‐1688.2011.00577.x Abstract: We compared the results of 12 recently calibrated regional SPARROW (SPAtially Referenced Regressions On Watershed attributes) models covering most of the continental United States to evaluate the consistency and regional differences in factors affecting stream nutrient loads. The models – 6 for total nitrogen and 6 for total phosphorus – all provide similar levels of prediction accuracy, but those for major river basins in the eastern half of the country were somewhat more accurate. The models simulate long‐term mean annual stream nutrient loads as a function of a wide range of known sources and climatic (precipitation, temperature), landscape (e.g., soils, geology), and aquatic factors affecting nutrient fate and transport. The results confirm the dominant effects of urban and agricultural sources on stream nutrient loads nationally and regionally, but reveal considerable spatial variability in the specific types of sources that control water quality. These include regional differences in the relative importance of different types of urban (municipal and industrial point vs. diffuse urban runoff) and agriculture (crop cultivation vs. animal waste) sources, as well as the effects of atmospheric deposition, mining, and background (e.g., soil phosphorus) sources on stream nutrients. Overall, we found that the SPARROW model results provide a consistent set of information for identifying the major sources and environmental factors affecting nutrient fate and transport in United States watersheds at regional and subregional scales.     

GENERALIZED WATERSHED LOADING FUNCTIONS FOR STREAM FLOW NUTRIENTS1

ABSTRACT: Loading functions are proposed as a general model for estimating monthly nitrogen and phosphorus fluxes in stream flow. The functions have a simple mathematical structure, describe a wide range of rural and urban nonpoint sources, and couple surface runoff and ground water discharge. Rural runoff loads are computed from daily runoff and erosion and monthly sediment yield calculations. Urban runoff loads are based on daily nutrient accumulation rates and exponential wash off functions. Ground water discharge is determined by lumped parameter unsaturated and saturated zone soil moisture balances. Default values for model chemical parameters were estimated from literature values. Validation studies over a three-year period for an 850 km² watershed showed that the loading functions explained at least 90 percent of the observed monthly variation in dissolved and total nitrogen and phosphorus fluxes in stream flow. Errors in model predictions of mean monthly fluxes were: dissolved phosphorus - 4 percent; total phosphorus - 2 percent; dissolved nitrogen - 18 percent; and total nitrogen - 28 percent. These results were obtained without model calibration.     

Water Quality Signals from Rural Land Use and Exurbanization in a Mountain Landscape: What's Clear and What's Confounded?

In mountainous landscapes with high climatic and geomorphic variability, how do rural land uses and exurbanization alter hydrology and water quality? We evaluated effects of rural land use and exurbanization on streamflows, suspended sediment concentrations and loads, specific conductance, and summer water temperatures in 12 streams and rivers within the Upper Little Tennessee River basin in the southern Appalachian Mountains. Eleven streams featured low levels of development (>61% forest cover) but differed in land use patterning, basin size, annual precipitation, and watershed morphology. One urban stream, located within the largest town in the basin, provided the high development comparative endpoint. Even low levels of rural development and exurbanization were associated with substantial increases in suspended sediment concentrations, sediment loads, and summer stream temperature daily maxima and diurnal variation. Observed summer temperature increases were much larger than would be expected due to global climate change over the next century. Specific conductance was idiosyncratic among the smaller streams. These water quality changes were not accompanied by streamflow changes that were discernible amid the high natural variation in precipitation and geomorphology. The water quality findings suggest the need for applying the best management practices, including riparian buffers, to even low levels of rural development.     

THE SENSITWITY OF NORTHERN SIERRA NEVADA STREAMFLOW TO CLIMATE CHANGE1

ABSTRACT: The sensitivity of streamflow to climate change was investigated in the American, Carson, and Truckee River Basins, California and Nevada. Nine gaging stations were used to represent streamflow in the basins. Annual models were developed by regressing 1961–1991 streamflow data on temperature and precipitation. Climate-change scenarios were used as inputs to the models to determine streamflow sensitivities. Climate-change scenarios were generated from historical time series by modifying mean temperatures by a range of +4°C to—4°C and total precipitation by a range of +25 percent to -25 percent. Results show that streamflow on the warmer, lower west side of the Sierra Nevada generally is more sensitive to temperature and precipitation changes than is streamflow on the colder, higher east side. A 2°C rise in temperature and a 25-percent decrease in precipitation results in stream-flow decreases of 56 percent on the American River and 25 percent on the Carson River. A 2°C decline in temperature and a 25-percent increase in precipitation results in streamflow increases of 102 percent on the American River and 22 percent on the Carson River.     

Sources and Transport of Phosphorus to Rivers in California and Adjacent States,U.S., as Determined by SPARROW Modeling

The SPARROW (SPAtially Referenced Regression on Watershed attributes) model was used to simulate annual phosphorus loads and concentrations in unmonitored stream reaches in California, U.S., and portions of Nevada and Oregon. The model was calibrated using de‐trended streamflow and phosphorus concentration data at 80 locations. The model explained 91% of the variability in loads and 51% of the variability in yields for a base year of 2002. Point sources, geological background, and cultivated land were significant sources. Variables used to explain delivery of phosphorus from land to water were precipitation and soil clay content. Aquatic loss of phosphorus was significant in streams of all sizes, with the greatest decay predicted in small‐ and intermediate‐sized streams. Geological sources, including volcanic rocks and shales, were the principal control on concentrations and loads in many regions. Some localized formations such as the Monterey shale of southern California are important sources of phosphorus and may contribute to elevated stream concentrations. Many of the larger point source facilities were located in downstream areas, near the ocean, and do not affect inland streams except for a few locations. Large areas of cultivated land result in phosphorus load increases, but do not necessarily increase the loads above those of geological background in some cases because of local hydrology, which limits the potential of phosphorus transport from land to streams.     

Streamflow Responses to Climate Change: Analysis of Hydrologic Indicators in a New York City Water Supply Watershed

Recent works have indicated that climate change in the northeastern United States is already being observed in the form of shorter winters, higher annual average air temperature, and more frequent extreme heat and precipitation events. These changes could have profound effects on aquatic ecosystems, and the implications of such changes are less understood. The objective of this study was to examine how future changes in precipitation and temperature translate into changes in streamflow using a physically based semidistributed model, and subsequently how changes in streamflow could potentially impact stream ecology. Streamflow parameters were examined in a New York City water supply watershed for changes from model‐simulated baseline conditions to future climate scenarios (2081‐2100) for ecologically relevant factors of streamflow using the Indicators of Hydrologic Alterations tool. Results indicate that earlier snowmelt and reduced snowpack advance the timing and increase the magnitude of discharge in the winter and early spring (November‐March) and greatly decrease monthly streamflow later in the spring in April. Both the rise and fall rates of the hydrograph will increase resulting in increased flashiness and flow reversals primarily due to increased pulses during winter seasons. These shifts in timing of peak flows, changes in seasonal flow regimes, and changes in the magnitudes of low flow can all influence aquatic organisms and have the potential to impact stream ecology.     

Dry Weather Water Quality Loadings in Arid,Urban Watersheds of the Los Angeles Basin,California, USA1

Abstract: Dry weather runoff in arid, urban watersheds may consist entirely of treated wastewater effluent and/or urban nonpoint source runoff, which can be a source of bacteria, nutrients, and metals to receiving waters. Most studies of urban runoff focus on stormwater, and few have evaluated the relative contribution and sources of dry weather pollutant loading for a range of constituents across multiple watersheds. This study assessed dry weather loading of nutrients, metals, and bacteria in six urban watersheds in the Los Angeles region of southern California to estimate relative sources of each constituent class and the proportion of total annual load that can be attributed to dry weather discharge. In each watershed, flow and water quality were sampled from storm drain and treated wastewater inputs, as well as from in‐stream locations during at least two time periods. Data were used to calculate mean concentrations and loads for various sources. Dry weather loads were compared with modeled wet weather loads under a range of annual rainfall volumes to estimate the relative contribution of dry weather load. Mean storm drain flows were comparable between all watersheds, and in all cases, approximately 20% of the flowing storm drains accounted for 80% of the daily volume. Wastewater reclamation plants (WRP) were the main source of nutrients, storm drains accounted for almost all the bacteria, and metals sources varied by constituent. In‐stream concentrations reflected major sources, for example nutrient concentrations were highest downstream of WRP discharges, while in‐stream metals concentrations were highest downstream of the storm drains with high metals loads. Comparison of wet vs. dry weather loading indicates that dry weather loading can be a significant source of metals, ranging from less than 20% during wet years to greater than 50% during dry years.     

EXPORT COEFFICIENT MODELING TO ASSESS PHOSPHORUS LOADING IN AN URBAN WATERSHED1

ABSTRACT: An export coefficient modeling approach was used to assess the influence of land use on phosphorus loading to a Southern Ontario stream. A model was constructed for the 1995–1996 water year and calibrated within ± 3 percent of the observed mean concentration of total phosphorus. It was found that runoff from urban areas contributed most to the loading of phosphorus to the stream. When the model was assessed by running it for the 1977–1978 water year, using water quality and land use data collected independently, agreement within ± 7 percent was obtained. The model was then used to forecast the impact of future urban development proposed for the watershed, in terms of phosphorus loading, and to evaluate the reduction in loading resulting from several urban best management practices (BMP). It was determined that phosphorus removal will have to be applied to all the urban runoff from the watershed to appreciably reduce stream phosphorus concentration. Of the BMP designs assessed, an infiltration pond system resulted in the greatest phosphorus load reduction, 50 percent from the 1995–1996 baseline.     

Integrated Modular Modeling of Water and Nutrients From Point and Nonpoint Sources in the Patuxent River Watershed1

Abstract: We present a simple modular landscape simulation model that is based on a watershed modeling framework in which different sets of processes occurring in a watershed can be simulated separately with different models. The model consists of three loosely coupled submodels: a rainfall‐runoff model (TOPMODEL) for runoff generation in a subwatershed, a nutrient model for estimation of nutrients from nonpoint sources in a subwatershed, and a stream network model for integration of point and nonpoint sources in the routing process. The model performance was evaluated using monitoring data in the watershed of the Patuxent River, a tributary to the Chesapeake Bay in Maryland, from July 1997 through August 1999. Despite its simplicity, the landscape model predictions of streamflow, and sediment and nutrient loads were as good as or better than those of the Hydrological Simulation Program‐Fortran model, one of the most widely used comprehensive watershed models. The landscape model was applied to predict discharges of water, sediment, silicate, organic carbon, nitrate, ammonium, organic nitrogen, total nitrogen, organic phosphorus, phosphate, and total phosphorus from the Patuxent watershed to its estuary. The predicted annual water discharge to the estuary was very close to the measured annual total in terms of percent errors for both years of the study period (≤2%). The model predictions for loads of nutrients were also good (20‐30%) or very good (<20%) with exceptions of sediment (40%), phosphate (36%), and organic carbon (53%) for Year 1.     

Effects of stream channelization on exports of nitrogen and phosphorus from north Carolina coastal plain watersheds

Nitrogen and phosphorus exports from channelizedstream watersheds were elevated over those from nearby natural swamp-stream watersheds. Nitrate exports were significantly greater from channelized-stream watersheds, and higher exports were attributed to faster groundwater drawdown, continual streamflow, and transformation of former floodplain to croplands following channelization. Exports of total organic nitrogen and total nitrogen were also significantly greater from channelized-stream watersheds. Differences in the exports of ammonium, filterable reactive phosphorus, and filterable unreactive phosphorus between the two watershed types were not detectable. Particulate phosphorus exports were significantly higher from channelized-stream watersheds, presumably because of greater erosion potential of nearby croplands and steep channel banks in the altered watersheds. The presence of nonpoint sources of pollution increased watershed exports of nutrients regardless of stream morphology. Examination of nutrient budgets for a portion of swamp floodplain at the base of one natural-stream watershed revealed that changes in local groundwater hydrology and stream morphology associated with channelization appeared to have greater effect on nutrient exports than simply the loss of bordering forested floodplain.     

IMPACTS OF CLIMATE CHANGE ON WATER YIELD IN THE UPPER WIND RIVER BASIN1

ABSTRACT: The potential impacts of climate change on water yield are examined in the Upper Wind River Basin. This is a high‐elevation, mountain basin with a snowfall/snowmelt dominated stream‐flow hydrograph. A variety of physiographic conditions are represented in the rangeland, coniferous forests, and high‐elevation alpine regions. The Soil Water Assessment Tool (SWAT) is used to model the baseline input time series data and climate change scenarios. Five hydroclimatic variables (temperature, precipitation, CO₂, radiation, and humidity) are examined using sensitivity tests of individual and coupled variables with a constant change and coupled variables with a monthly change. Results indicate that the most influential variable on annual water yield is precipitation; and, the most influential variable on the timing of streamflow is temperature. Carbon dioxide, radiation, and humidity each noticeably impact water yield, but less significantly. The coupled variable analyses represent a more realistic climate change regime and reflect the combined response of the basin to each variable; for example, increased temperature offsets the effects of increased precipitation and magnifies the effects of decreased precipitation. This paper shows that the hydrologic response to climate change depends largely on the hydroclimatic variables examined and that each variable has a unique effect (e.g., magnitude, timing) on water yield.     

LOADING FUNCTIONS FOR PREDICTING NUTRIENT LOSSES FROM COMPLEX WATERSHEDS1

ABSTRACT: A loading function methodology is presented for predicting runoff, sediment, and nutrient losses from complex watersheds. Separate models are defined for cropland, forest, urban and barnyard sources, and procedures for estimating baseflow nutrients are provided. The loading functions are designed for use as a preliminary screening tool to isolate the major contributors in a watershed. Input data sources are readily available and the functions do not require costly calibrations. Data requirements include watershed land use and soil information, daily precipitation and temperature records and rainfall erosivities. Comparison of predicted and measured water, sediment, and nutrient runoff fluxes for the West Branch Deleware River in New York, indicated that runoff was underpredicted by about 14 percent while dissolved nutrients were within 30 percent of observed values. Sediment and solid-phase nutrients were overpredicted by about 50 percent. An annual nutrient budget for the West Branch Delaware River showed that cornland was the major source of sediment, solid phase nutrients, and total phosphorus. Waste water treatment plants and ground water discharge contributed the most dissolved phosphorus and dissolved nitrogen, respectively.     

NUTRIENT MASS BALANCE FOR THE ALBEMARLE-PAMLICO DRAINAGE BASIN,NORTH CAROLINA AND VIRGINIA, 19901

ABSTRACT: A 1990 nitrogen and phosphorus mass balance calculated for eight National Stream Quality Accounting Network (NASQAN) basins in the Albemarle-Pamlico Drainage Basin indicated the importance of agricultural nonpoint sources of nitrogen and phosphorus and watershed nitrogen retention and processing capabilities. Basin total nitrogen and phosphorus input estimates were calculated for atmospheric deposition (which averaged 27 percent of total nitrogen inputs and 22 percent of total phosphorus inputs); crop fertilizer (27 and 25 percent); animal-waste (22 and 50 percent, respectively); point sources (3 percent each of total nitrogen and total phosphorus inputs); and biological nitrogen fixation (21 percent of total nitrogen inputs). Highest in-stream nitrogen and phosphorus loads were measured in predominantly agricultural drainage areas. Intermediate loads were observed in mixed agricultural/urban drainage areas; the lowest loads were measured in mixed agricultural/forested drainage areas. The difference between the sum of the nutrient input categories and the sum of the in-stream nutrient loads and crop-harvest nutrient removal was assigned to a residual category for the basin. The residual category averaged 51 percent of total nitrogen inputs and 54 percent of total phosphorus inputs.     

Streamflow Changes in the South Atlantic,United States During the Mid‐ and Late 20th Century1

Abstract: Repeated severe droughts over the last decade in the South Atlantic have raised concern that streamflow may be systematically decreasing, possibly due to climate variability. We examined the monthly and annual trends of streamflow, precipitation, and temperature in the South Atlantic for the time periods: 1934‐2005, 1934‐1969, and 1970‐2005. Streamflow and climate (temperature and precipitation) trends transitioned ca. 1970. From 1934 to 1969, streamflow and precipitation increased in southern regions and decreased in northern regions; temperature decreased throughout the South Atlantic. From 1970 to 2005, streamflow decreased, precipitation decreased, and temperature increased throughout the South Atlantic. It is unclear whether these will be continuing trends or simply part of a long‐term climatic oscillation. Whether these streamflow trends have been driven by climatic or anthropogenic changes, water resources management faces challenging prospects to adapt to decadal‐scale persistently wet and dry hydrologic conditions.     

Sensitivity of Stream flow and Water Table Depth to Potential Climatic Variability in a Coastal Forested Watershed1

Dai, Zhaohua, Carl C. Trettin, Changsheng Li, Devendra M. Amatya, Ge Sun, and Harbin Li, 2010. Sensitivity of Streamflow and Water Table Depth to Potential Climatic Variability in a Coastal Forested Watershed. Journal of the American Water Resources Association (JAWRA) 1–13. DOI: 10.1111/j.1752-1688.2010.00474.x Abstract: A physically based distributed hydrological model, MIKE SHE, was used to evaluate the effects of altered temperature and precipitation regimes on the streamflow and water table in a forested watershed on the southeastern Atlantic coastal plain. The model calibration and validation against both streamflow and water table depth showed that the MIKE SHE was applicable for predicting the streamflow and water table dynamics for this watershed with an acceptable model efficiency (E > 0.5 for daily streamflow and >0.75 for monthly streamflow). The simulation results from changing temperature and precipitation scenarios indicate that climate change influences both streamflow and water table in the forested watershed. Compared to current climate conditions, the annual average streamflow increased or decreased by 2.4% with one percentage increase or decrease in precipitation; a quadratic polynomial relationship between changes in water table depth (cm) and precipitation (%) was found. The annual average water table depth and annual average streamflow linearly decreased with an increase in temperature within the range of temperature change scenarios (0-6°C). The simulation results from the potential climate change scenarios indicate that future climate change will substantially impact the hydrological regime of upland and wetland forests on the coastal plain with corresponding implications to altered ecosystem functions that are dependent on water.