Spatial and Temporal Evaluation of Hydrological Response to Climate and Land Use Change in Three South Dakota Watersheds

This study analyzed changes in hydrology between two recent decades (1980s and 2010s) with the Soil and Water Assessment Tool (SWAT) in three representative watersheds in South Dakota: Bad River, Skunk Creek, and Upper Big Sioux River watersheds. Two SWAT models were created over two discrete time periods (1981‐1990 and 2005‐2014) for each watershed. National Land Cover Datasets 1992 and 2011 were, respectively, ingested into 1981‐1990 and 2005‐2014 models, along with corresponding weather data, to enable comparison of annual and seasonal runoff, soil water content, evapotranspiration (ET), water yield, and percolation between these two decades. Simulation results based on the calibrated models showed that surface runoff, soil water content, water yield, and percolation increased in all three watersheds. Elevated ET was also apparent, except in Skunk Creek watershed. Differences in annual water balance components appeared to follow changes in land use more closely than variation in precipitation amounts, although seasonal variation in precipitation was reflected in seasonal surface runoff. Subbasin‐scale spatial analyses revealed noticeable increases in water balance components mostly in downstream parts of Bad River and Skunk Creek watersheds, and the western part of Upper Big Sioux River watershed. Results presented in this study provide some insight into recent changes in hydrological processes in South Dakota watersheds. Editor's note: This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.     

Impact of Changing Climate and Land Cover on Flood Magnitudes in the Delaware River Basin,USA

Changing climate and land cover are expected to impact flood hydrology in the Delaware River Basin over the 21st Century. HEC‐HMS models (U.S. Army Corps of Engineers Hydrologic Engineering Center‐Hydrologic Modeling System) were developed for five case study watersheds selected to represent a range of scale, soil types, climate, and land cover. Model results indicate that climate change alone could affect peak flood discharges by ?6% to +58% a wide range that reflects regional variation in projected rainfall and snowmelt and local watershed conditions. Land cover changes could increase peak flood discharges up to 10% in four of the five watersheds. In those watersheds, the combination of climate and land cover change increase modeled peak flood discharges by up to 66% and runoff volumes by up to 44%. Precipitation projections are a key source of uncertainty, but there is a high likelihood of greater precipitation falling on a more urbanized landscape that produces larger floods. The influence of climate and land cover changes on flood hydrology for the modeled watersheds varies according to future time period, climate scenario, watershed land cover and soil conditions, and flood frequency. The impacts of climate change alone are typically greater than land cover change but there is substantial geographic variation, with urbanization the greater influence on some small, developing watersheds.     

Agricultural practices influence flow regimes of headwater streams in western Iowa

Agricultural tillage influences runoff and infiltration, but consequent effects on watershed hydrology are poorly documented. This study evaluated 25 yr (1971-1995) hydrologic records from four first-order watersheds in Iowa's loess hills. Two watersheds were under conventional tillage and two were under conservation (ridge) tillage, one of which was terraced. All four watersheds grew corn (Zea mays L.) every year. Flow-frequency statistics and autoregressive modeling were used to determine how conservation treatments influenced stream hydrology. The autoregressive modeling characterized variations in discharge, baseflow, and runoff at multi-year, annual, and shorter time scales. The ridge-tilled watershed (nonterraced) had 47% less runoff and 36% more baseflow than the conventional watershed of similar landform and slope. Recovery of baseflow after drought was quicker in the conservation watersheds, as evidenced by 365-d moving average plots, and 67% greater baseflow during the driest 2 yr. The two conventional watersheds were similar, except the steeper watershed discharged more runoff and baseflow during short (<30 d), wet periods. Significant multi-year and annual cycles occurred in all variables. Under ridge-till, seasonal (annual-cycle) variations in baseflow had greater amplitude, showing the seasonality of subsurface contaminant movement could increase under conservation practices. However, deviations from the modeled cycles of baseflow were also more persistent under conservation practices, indicating baseflow was more stable. Indeed, flow-frequency curves showed wet-weather discharge decreased and dry-weather discharge increased under conservation practices. Although mean discharge increased in the conservation watersheds, variance and skewness of daily values were smaller. Ridge tillage with or without terraces increased stream discharge but reduced its variability.     

Impacts of Climate Variability and Land Use Alterations on Frequency Distributions of Terrestrial Runoff Loading to Coastal Waters in Southern California1

Abstract: The transport of water, sediment, dissolved and particulate chemicals, and bacteria from coastal watersheds affects the nearshore marine and estuarine waters. In southern California, coastal watersheds deliver water and associated constituents to the nearshore system in discrete pulses. To better understand the pulsed nature of these watersheds, frequency distributions of simulated runoff events are presented for: (1) three land use conditions (1929, 1998, 2050); (2) three time periods (all water years 1989‐2002), only El Nino years (1992, 1993, 1995, 1998); and only non‐El Nino years; and (3) three regions (watershed, uplands, and lowlands). At the watershed scale, there was a significant increase (>200%) in mean event runoff from 1929 to 2050 (0.4‐1.3 cm) due to localized urbanization, which shifted the dominant sources of runoff from the mountains in 1929 (78% of watershed runoff) to the coastal plane for 2050 conditions (51% of watershed runoff). Inter‐annual climate variability was strong in the rainfall and runoff frequency distributions, with mean event rainfall and runoff 66 and 60% larger in El Nino relative to non‐El Nino years. Combining urbanization and climate variability, 2050 land conditions resulted in El Nino years being five times more likely to produce large (>3.0 cm) runoff events relative to non‐El Nino years. Combining frequency distributions of event runoff with regional nutrient export relationships, we show that in El Nino years, one in five events produced runoff ≥2.5 cm and temporary nearshore nitrate and phosphate concentrations of 12 and 1.4 μM, respectively, or approximately 5‐10 times above ambient conditions.     

Modeling Hydrology in a Small Rocky Mountain Watershed Serving Large Urban Populations1

Abstract: This article describes the development of a calibrated hydrologic model for the Blue River watershed (867 km²) in Summit County, Colorado. This watershed provides drinking water to over a third of Colorado’s population. However, more research on model calibration and development for small mountain watersheds is needed. This work required integration of subsurface and surface hydrology using GIS data, and included aspects unique to mountain watersheds such as snow hydrology, high ground‐water gradients, and large differences in climate between the headwaters and outlet. Given the importance of this particular watershed as a major urban drinking‐water source, the rapid development occurring in small mountain watersheds, and the importance of Rocky Mountain water in the arid and semiarid West, it is useful to describe calibrated watershed modeling efforts in this watershed. The model used was Soil and Water Assessment Tool (SWAT). An accurate model of the hydrologic cycle required incorporation of mountain hydrology‐specific processes. Snowmelt and snow formation parameters, as well as several ground‐water parameters, were the most important calibration factors. Comparison of simulated and observed streamflow hydrographs at two U.S. Geological Survey gaging stations resulted in good fits to average monthly values (0.71 Nash‐Sutcliffe coefficient). With this capability, future assessments of point‐source and nonpoint‐source pollutant transport are possible.     

Temporal Changes in Streamflow and Attribution of Changes to Climate and Landuse in Wisconsin Watersheds

Previous historic trends analyses on 21st Century hydrologic data in the United States generally focus on annual flow statistics and have continued to use USGS hydro‐climatic data network (HCDN) stations, although post‐1988 diversions and runoff regulations are not reflected in the HCDN. Using a more recent dataset, Geospatial Attributes of Gages for Evaluating Streamflow, version II (GAGES II), compiled by Falcone (2012), which includes more watersheds with reference conditions, a comprehensive analysis of changes in seasonal, and annual streamflow in Wisconsin watersheds is demonstrated. Given the pronounced influence of seasonal hydrology in Wisconsin watersheds, the objective of this study is to elucidate the nature of temporal (annual, seasonal, and monthly) changes in runoff. Considerable temporal and regional variability was found in annual and seasonal streamflow changes between the two historic periods 1951‐1980 and 1981‐2010 considered in the study. For example, the northern watersheds show relatively small changes in streamflow discharge ranging from ?6.0 to 4.2%, while the southern watersheds show relatively large increases in streamflow discharge ranging from 13.1 to 18.2%. To apportion streamflow changes to climate and nonclimatic factors, a method based on potential evapotranspiration changes is demonstrated. Results show that nonclimatic factors account for more than 60% of changes in annual runoff in Wisconsin watersheds considered in the study.     

A Hydrologic Data Screening Procedure for Exploring Monotonic Trends and Shifts in Rainfall and Runoff Patterns

A thorough understanding of past and present hydrologic responses to changes in precipitation patterns is crucial for predicting future conditions. The main objectives of this study were to determine temporal changes in rainfall‐runoff relationship and to identify significant trends and abrupt shifts in rainfall and runoff time series. Ninety‐year rainfall and runoff time series datasets from the Gasconade and Meramec watersheds in east‐central Missouri were used to develop data screening procedure to assess changes in the rainfall and runoff temporal patterns. A statistically significant change in mean and variance was detected in 1980 in the rainfall and runoff time series within both watersheds. In addition, both the rainfall and runoff time series indicated the presence of nonstationary attributes such as statistically significant monotonic trends and/or change in mean and variance, which should be taken into consideration when using the time series to predict future scenarios. The annual peak runoff and the annual low flow in the Meramec watershed showed significant temporal changes compared to that in the Gasconade watershed. Water loss in both watersheds was found to be significantly increasing which is potentially due to the increase in groundwater pumping for water supply purposes.     

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.     

Hydrology and Water Budget for a Forested Atlantic Coastal Plain Watershed,South Carolina1

Abstract: Increases in timber demand and urban development in the Atlantic Coastal Plain over the past decade have motivated studies on the hydrology, water quality, and sustainable management of coastal plain watersheds. However, studies on baseline water budgets are limited for the low‐lying, forested watersheds of the Atlantic Coastal Plain. The purpose of this study was to document the hydrology and a method to quantify the water budget of a first‐order forested watershed, WS80, located within the USDA Forest Service Santee Experimental Forest northeast of Charleston, South Carolina. Annual Rainfall for the 2003 and 2004 periods were 1,671 mm (300 mm above normal) and 962 mm (over 400 mm below normal), respectively. Runoff coefficients (outflow as a fraction of total rainfall) for the 2003 and 2004 periods were 0.47 and 0.08, respectively, indicating a wide variability of outflows as affected by antecedent conditions. A spreadsheet‐based Thornthwaite monthly water balance model was tested on WS80 using three different potential evapotranspiration estimators [Hamon, Thornthwaite, and Penman‐Monteith (P‐M)]. The Hamon and P‐M‐based methods performed reasonably well with average absolute monthly deviations of 12.6 and 13.9 mm, respectively, between the measured and predicted outflows. Estimated closure errors were all within 9% for the 2003, 2004, and seasonal water budgets. These results may have implications on forest management practices and provide necessary baseline or reference information for Atlantic Coastal Plain watersheds.     

DIFFERENCES BETWEEN MODELED SURPLUS AND USGS-MEASURED DISCHARGE IN LAKE PONTCHARTRAIN BASIN,LOUISIANA1

ABSTRACT: A one-layer decreasing-availability monthly water balance model is used to estimate monthly surplus that flows into the Lake Pontchartrain Basin from the Amite, Tickfaw, Natalbany, Tangipahoa, and Tchefuncte Rivers for water years 1949 through 1990. The modeled annual surplus for each drainage basin is compared to gauged annual discharge obtained from the United States Geological Survey. This provides an estimate of the differential success of the model over watersheds of various sizes, and also suggests appropriate adjustment factors to be used in future water balance analyses of similar basins in humid subtropical climate regions. Results show that annual surplus values agree well with the USGS values, after an annual adjustment of about 140 mm (11 to 28 percent of the basin surplus) is subtracted from the annual modeled totals to compensate for overestimation by the model. However, inter-annual variability is high in the annual cycles. Winter and spring discharges can also be modeled successfully.     

Evaluating agricultural best management practices in tile-drained subwatersheds of the Mackinaw River, Illinois

Best management practices (BMPs) are widely promoted in agricultural watersheds as a means of improving water quality and ameliorating altered hydrology. We used a paired watershed approach to evaluate whether focused outreach could increase BMP implementation rates and whether BMPs could induce watershed-scale (4000 ha) changes in nutrients, suspended sediment concentrations, or hydrology in an agricultural watershed in central Illinois. Land use was >90% row crop agriculture with extensive subsurface tile drainage. Outreach successfully increased BMP implementation rates for grassed waterways, stream buffers, and strip-tillage within the treatment watershed, which are designed to reduce surface runoff and soil erosion. No significant changes in nitrate-nitrogen (NO-N), total phosphorus (TP), dissolved reactive phosphorus, total suspended sediment (TSS), or hydrology were observed after implementation of these BMPs over 7 yr of monitoring. Annual NO-N export (39-299 Mg) in the two watersheds was equally exported during baseflow and stormflow. Mean annual TP export was similar between the watersheds (3.8 Mg) and was greater for TSS in the treatment (1626 ± 497 Mg) than in the reference (940 ± 327 Mg) watershed. Export of TP and TSS was primarily due to stormflow (>85%). Results suggest that the BMPs established during this study were not adequate to override nutrient export from subsurface drainage tiles. Conservation planning in tile-drained agricultural watersheds will require a combination of surface-water BMPs and conservation practices that intercept and retain subsurface agricultural runoff. Our study emphasizes the need to measure conservation outcomes and not just implementation rates of conservation practices.     

Hydrologic Landscape Classification to Estimate Bristol Bay,Alaska Watershed Hydrology

Hydrologic landscapes (HLs) have proven to be a useful tool for broad scale assessment and classification of landscapes across the United States as they help organize larger geographical areas into areas of similar hydrologic characteristics. We developed a HL classification for the Bristol Bay watershed of southwest Alaska that incorporates indices of annual climate and seasonality, terrain, geology, and the influences of large lakes and glaciers. A HL classification is particularly useful in this large watershed because of its hydrologic and landscape variability, important salmon fishery, variety of environmental and potential anthropogenic stressors, and lack of widespread hydrologic data. Following creation of Bristol Bay basin‐wide HL classes, we compared the HL distributions within watersheds grouped by two calculated runoff parameters derived from available long‐term streamflow records and found HL distributions within these groups provided predictive insight on hydrologic behavior. Using these developed runoff groups, we estimated expected hydrologic behavior in watersheds across the larger Bristol Bay watershed that lacked gauged streamflow records. The HL approach provides a scientific basis for estimating the first‐order hydrologic behavior of watersheds and landscapes that lack detailed hydrologic information.     

The Influence of Forest Regrowth on the Stream Discharge in the North Carolina Piedmont Watersheds

This study focuses on the relationships of watershed runoff with historical land use/land cover (LULC) and climate trends. Over the 20th Century, LULC in the Southeast United States, particularly the North Carolina Piedmont, has evolved from an agriculture dominated to an extensively forested landscape with more recent localized urbanization. The regrowth of forest has an important influence on the hydrology of the region as it enhances ecosystem interaction with recent climate change. During 1920‐2009, the amount of precipitation in some parts of the North Carolina Piedmont forest regrowth area showed increasing trends without corresponding increments in runoff. We employed the Soil and Water Assessment Tool (SWAT) to backcast long‐term hydrologic behavior of watersheds in North Carolina with different LULC conditions: (1) LULC conversion from agricultural to forested area and (2) long‐term stable forested area. Comparing U.S. Geological Survey‐measured stream discharge with SWAT‐simulated stream discharge under the assumption of constant 2006 LULC, we found significant stream discharge underprediction by SWAT in two LULC conversion watersheds during the early simulation period (1920s) with differences gradually decreasing by the mid‐1970s. This model bias suggests that forest regrowth on abandoned agricultural land was a key factor contributing to mitigate the impact of increased precipitation on runoff due to increasing water consumption driven by changes in vegetation.     

RUNOFF RESPONSE FROM TWO RECLAIMED WATERSHEDS1

ABSTRACT: Detailed studies of the surface hydrology of reclaimed surface-mined watersheds for both rainfall and snowmelt events are non-existent for central Alberta yet this information is crucial for design of runoff conveyance and storage structures. A study was initiated in 1992 with principal objectives of quantifying surface runoff for both summer rainfall and spring snowmelt events and identifying the dominant flow processes occurring in two reclaimed watersheds. Snowmelt accounted for 86 and 100% of annual watershed runoff in 1993 and 1994, respectively. The highest instantaneous peak flow was recorded during a summer rainfall event with a return period of greater than 50 years. Infiltration-excess overland flow was identified as the dominant flow process occurring within the Sandy Subsoil Watershed, whereas saturation overland flow was the principal runoff process occurring within the West Watershed.     

Estimating Forest Ecosystem Evapotranspiration at Multiple Temporal Scales With a Dimension Analysis Approach1

Abstract:  It is critical that evapotranspiration (ET) be quantified accurately so that scientists can evaluate the effects of land management and global change on water availability, streamflow, nutrient and sediment loading, and ecosystem productivity in watersheds. The objective of this study was to derive a new semi‐empirical ET modeled using a dimension analysis method that could be used to estimate forest ET effectively at multiple temporal scales. The model developed describes ET as a function of water availability for evaporation and transpiration, potential ET demand, air humidity, and land surface characteristics. The model was tested with long‐term hydrometeorological data from five research sites with distinct forest hydrology in the United States and China. Averaged simulation error for daily ET was within 0.5 mm/day. The annual ET at each of the five study sites were within 7% of measured values. Results suggest that the model can accurately capture the temporal dynamics of ET in forest ecosystems at daily, monthly, and annual scales. The model is climate‐driven and is sensitive to topography and vegetation characteristics and thus has potential to be used to examine the compounding hydrologic responses to land cover and climate changes at multiple temporal scales.     

Flood Risk Reduction from Agricultural Best Management Practices

Best management practices (BMPs) play an important role in improving impaired water quality from conventional row crop agriculture. In addition to reducing nutrient and sediment loads, BMPs such as fertilizer management, reduced tillage, and cover crops could alter the hydrology of agricultural systems and reduce surface water runoff. While attention is devoted to the water quality benefits of BMPs, the potential co‐benefits of flood loss reduction are often overlooked. This study quantifies the effects of selected commonly applied BMPs on expected flood loss to agricultural and urban areas in four Iowa watersheds. The analysis combines a watershed hydrologic model, hydraulic model outputs, and a loss estimation model to determine relationships between hydrologic changes from BMP implementations and annual economic flood loss. The results indicate a modest reduction in peak discharge and economic loss, although loss reduction is substantial when urban centers or other high‐value assets are located downstream in the watershed. Among the BMPs, wetlands, and cover crops reduce losses the most. The research demonstrates that watershed‐scale implementation of agricultural BMPs could provide benefits of flood loss reduction in addition to water quality improvements.     

MODEL PREDICTIONS OF WATERSHED HYDROLOGIC COMPONENTS: COMPARISON AND VERIFICATIION1

ABSTRACT: Model predictions of the relatively simple soil compartment model SESOIL are compared with those of the more data-intensive terrestrial ecosystem hydrology model AGTEHM. Comparisons were performed using data from a deciduous forest stand watershed, a grassland watershed, and two agricultural field plots. Good agreement was obtained between model predictions for annual values of infiltration, evapotranspiration, surface runoff, and groundwater runoff. SESOIL model predictions also compare well with empirical measurements at the forest stand and the grassland watersheds.     

Effects of Watershed Impervious Cover on Dissolved Silica Loading in Storm Flow1

Abstract: Dissolved silica (DSi) availability is a factor that affects the composition of algal populations in aquatic ecosystems. DSi cycling is tightly linked to the hydrological cycle, which is affected by human alterations of the landscape. Development activities that increase impervious cover change watershed hydrology and may increase the discharge of DSi‐poor rainwater and decrease the discharge of DSi‐rich ground water into aquatic ecosystems, possibly shifting algal community composition toward less desirable assemblages. In this study, DSi loadings from two adjacent coastal watersheds with different percent impervious cover were compared during four rain and five nonrain events. Loadings in the more impervious watershed contained a significantly larger proportion of surface runoff than base flow (ground‐water discharge) and had lower [DSi] water during rain events than the less impervious watershed. Application of the Soil Conservation Service Curve Number (CN) method showed that the minimum rainfall height necessary to yield runoff was significantly lower for the more impervious watershed, implying that runoff volumes increase with impervious cover as well as the frequency of runoff‐yielding events. Empirical data collected during this study and estimates derived from the CN method suggest that impervious cover may be responsible for both short‐term DSi limitation during rain events as well as long‐term reduction of DSi inputs into aquatic ecosystems.     

Assessing Critical Source Areas in Watersheds for Conservation Buffer Planning and Riparian Restoration

A science-based geographic information system (GIS) approach is presented to target critical source areas in watersheds for conservation buffer placement. Critical source areas are the intersection of hydrologically sensitive areas and pollutant source areas in watersheds. Hydrologically sensitive areas are areas that actively generate runoff in the watershed and are derived using a modified topographic index approach based on variable source area hydrology. Pollutant source areas are the areas in watersheds that are actively and intensively used for such activities as agricultural production. The method is applied to the Neshanic River watershed in Hunterdon County, New Jersey. The capacity of the topographic index in predicting the spatial pattern of runoff generation and the runoff contribution to stream flow in the watershed is evaluated. A simple cost-effectiveness assessment is conducted to compare the conservation buffer placement scenario based on this GIS method to conventional riparian buffer scenarios for placing conservation buffers in agricultural lands in the watershed. The results show that the topographic index reasonably predicts the runoff generation in the watershed. The GIS-based conservation buffer scenario appears to be more cost-effective than the conventional riparian buffer scenarios.     

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.