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.     

Weighting Nitrogen and Phosphorus Pixel Pollutant Loads to Represent Runoff and Buffering Likelihoods

Watershed models often estimate annual nitrogen (N) or phosphorus (P) pollutant loads in rural areas with export coefficient (EC) (kg/ha/yr) values based on land cover, and in urban areas as the product of spatially uniform event mean concentration (EMC) (mg/L) values and runoff volume. Actual N and P nonpoint source (NPS) pollutant loading has more spatial complexity due to watershed variation in runoff likelihood and buffering likelihood along surface and subsurface pathways, which can be represented in a contributing area dispersal area (CADA) NPS model. This research develops a CADA NPS model to simulate how watershed properties of elevation, land cover, and soils upslope and downslope of each watershed pixel influence nutrient loading. The model uses both surface and subsurface runoff indices (RI), and surface and subsurface buffer indices (BI), to quantify the runoff and buffering likelihood for each watershed pixel, and generate maps of weighted EC and EMC values that identify NPS pollutant loading hotspots. The research illustrates how CADA NPS model maps and pixel loading values are sensitive to the spatial resolution and accuracy of elevation and land cover data, and model predictions can represent the lower and upper bounds of NPS loading. The model provides managers with a tool to rapidly visualize, rank, and investigate likely areas of high nutrient export.     

From agricultural intensification to conservation: sediment transport in the Raccoon River, Iowa, 1916-2009

Fluvial sediment is a ubiquitous pollutant that negatively affects surface water quality and municipal water supply treatment. As part of its routine water supply monitoring, the Des Moines Water Works (DMWW) has been measuring turbidity daily in the Raccoon River since 1916. For this study, we calibrated daily turbidity readings to modern total suspended solid (TSS) concentrations to develop an estimation of daily sediment concentrations in the river from 1916 to 2009. Our objectives were to evaluate long-term TSS patterns and trends, and relate these to changes in climate, land use, and agricultural practices that occurred during the 93-yr monitoring period. Results showed that while TSS concentrations and estimated sediment loads varied greatly from year to year, TSS concentrations were much greater in the early 20th century despite drier conditions and less discharge, and declined throughout the century. Against a backdrop of increasing discharge in the Raccoon River and widespread agricultural adaptations by farmers, sediment loads increased and peaked in the early 1970s, and then have slowly declined or remained steady throughout the 1980s to present. With annual sediment load concentrated during extreme events in the spring and early summer, continued sediment reductions in the Raccoon River watershed should be focused on conservation practices to reduce rainfall impacts and sediment mobilization. Overall, results from this study suggest that efforts to reduce sediment load from the watershed appear to be working.     

WATER QUALITY IN AGRICULTURAL,URBAN, AND MIXED LAND USE WATERSHEDS1

ABSTRACT: Water quality and nonpoint source (NPS) pollution are important issues in many areas of the world, including the Inner Bluegrass Region of Kentucky where urban development is changing formerly rural watersheds into urban and mixed use watersheds. In watersheds where land use is mixed, the relative contributions of NPS pollution from rural and urban land uses can be difficult to separate. To better understand NPS pollution sources in mixed use watersheds, surface water samples were taken at three sites that varied in land use to examine the effect of land use on water quality. Within the group of three watersheds, one was predominately agriculture (Agricultural), one was predominately urban (Urban), and a third had relatively equal representation of both types of land uses (Mixed). Nitrogen (N), phosphorus (P), total suspended solids (TSS), turbidity, pH, temperature, and streamflow were measured for one year. Comparisons are made among watersheds for concentration and fluxes of water quality parameters. Nitrate and orthophosphate concentrations were found to be significantly higher in the Agricultural watershed. Total suspended solids, turbidity, temperature, and pH, were found to be generally higher in the Urban and Mixed watersheds. No differences were found for streamflow (per unit area), total phosphorus, and ammonium concentrations among watersheds. Fluxes of orthophosphate were greater in the Agricultural watershed that in the Urban watershed while fluxes of TSS were greater in the Mixed watershed when compared to the Agricultural watershed. Fluxes of nitrate, ammonium, and total phosphorus did not vary among watersheds. It is apparent from the data that Agricultural land uses are generally a greater source of nutrients than the Urban land uses while Urban land uses are generally a greater source of suspended sediment.     

On the Role of Spatial,Temporal, and Climatic Forces on Stream Sediment Loading from Rural and Urban Ecosystems

In this study, we characterize the greatest sediment loading events by their sediment delivery behavior; dominant climate, watershed, and antecedent conditions; and their seasonal distribution for rural and urban land uses. The study area is Paradise Creek Watershed, a mixed land use watershed in northern Idaho dominated by saturation excess processes in the upstream rural area and infiltration excess in the downstream urban area. We analyzed 12 years of continuous streamflow, precipitation, and watershed data at two monitoring stations. We identified 137 sediment loading events in the upstream rural section of the watershed and 191 events in the downstream urban section. During the majority of these events conditions were transport limited and the sediment flush occurred early in the event, generally in the first 20% of elapsed event time. Statistical analysis including two dozen explanatory variables showed peak discharge, event duration, and antecedent baseflow explained most of the variation in event sediment load at both stations and for the watershed as a whole (R² = 0.73‐0.78). In the rural area, saturated soils combined with spring snowmelt in March led to the greatest loading events. The urban area load contribution peaked in January, which could be a re‐suspension of streambed sediments from the previous water year. Throughout the study period, one event contributed, on average, 33% of the annual sediment load but only accounted for 2% of the time in a year.     

BUILDUP,WASHOFF, AND EVENT MEAN CONCENTRATIONS1

ABSTRACT: Most watershed water quality simulation models require the user to specify pollutant buildup and washoff rate parameters for pollutants, by land use. Buildup and washoff rates are difficult to measure directly, and only limited guidance and few observed data are available from the literature. Many studies, however, report storm event mean concentrations (EMCs). These EMCs must arise as a result of the buildup and washoff processes, but typically represent the net contribution from a variety of pervious and impervious surfaces. This paper explores the relationship between EMCs and buildup/washoff parameters. An assumption of the mathematical form of the buildup/washoff relationship gives an algebraic expression for the EMC consistent with model assumptions. This yields techniques to separate observed EMCs into contributions from different land uses and from pervious and impervious surfaces. Given this relationship, numerical optimization may be used to estimate site specific values of buildup and washoff parameters from observed storm EMCs for use in modeling. Use of this approach helps ensure that model parameters are consistent with observed data, providing a rational starting point for final model calibration. Several site examples demonstrate use of the method.     

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.     

LAND USE AND AQUATIC BIOINTEGRITY IN THE BLACKFOOT RWER WATERSHED,MONTANA1

ABSTRACT: Benthic macroinvertebrate samples representing 151 taxa were collected in August 1995 to examine the linkage between land use, water quality, and aquatic biointegrity in seven tributaries of the Blackfoot River watershed, Montana. The tributaries represent silvicultural (timber harvesting), agricultural (irrigated alfalfa and hay and livestock grazing), and wilderness land uses. A 2.4 km (1.5 mile) reach of a recently restored tributary also was sampled for comparison with the other six sites. A geographic information system (GIS) was used to characterize the seven subwatersheds and estimate soil erosion, using the Modified Universal Soil Loss Equation, and sediment delivery. The wilderness stream had the highest aquatic biointegrity. Two agricultural streams had the largest estimated soil erosion and sediment delivery rates, the greatest habitat impairment from nonpoint source pollution, and the most impoverished macroinvertebrate communities. The silvicultural subwatersheds had greater rates of estimated soil erosion and sediment delivery and lower aquatic biointegrity than the wilderness reference site but evinced better conditions than the agricultural sites. A multiple-use (forestry, grazing, and wildlife management) watershed and the restored site ranked between the silvicultural and agricultural sites. This spectrum of land use and aquatic biointegrity illustrates both the challenges and opportunities that define watershed management.     

Effects of calibration on L-THIA GIS runoff and pollutant estimation

Urbanization can result in alteration of a watershed's hydrologic response and water quality. To simulate hydrologic and water quality impacts of land use changes, the Long-Term Hydrologic Impact Assessment (L-THIA) system has been used. The L-THIA system estimates pollutant loading based on direct runoff quantity and land use based pollutant coefficients. The accurate estimation of direct runoff is important in assessing water quality impacts of land use changes. An automated program was developed to calibrate the L-THIA model using the millions of curve number (CN) combinations associated with land uses and hydrologic soil groups. L-THIA calibration for the Little Eagle Creek (LEC) watershed near Indianapolis, Indiana was performed using land use data for 1991 and daily rainfall data for six months of 1991 (January 1-June 30) to minimize errors associated with use of different temporal land use data and rainfall data. For the calibration period, the Nash-Sutcliffe coefficient was 0.60 for estimated and observed direct runoff. The calibrated CN values were used for validation of the model for the same year (July 1-December 31), and the Nash-Sutcliffe coefficient was 0.60 for estimated and observed direct runoff. The Nash-Sutcliffe coefficient was 0.52 for January 1, 1991 to December 31, 1991 using uncalibrated CN values. As shown in this study, the use of better input parameters for the L-THIA model can improve accuracy. The effects on direct runoff and pollutant estimation of the calibrated CN values in the L-THIA model were investigated for the LEC. Following calibration, the estimated average annual direct runoff for the LEC watershed increased by 34%, total nitrogen by 24%, total phosphorus by 22%, and total lead by 43%. This study demonstrates that the L-THIA model should be calibrated and validated prior to application in a particular watershed to more accurately assess the effects of land use changes on hydrology and water quality.     

NUTRIENT CONTRIBUTION OF NONPOINT SOURCE RUNOFF IN THE LAS VEGAS VALLEY1

ABSTRACT: A Geographic Information System (GIS) based non‐point source runoff model is developed for the Las Vegas Valley, Nevada, to estimate the nutrient loads during the years 2000 and 2001. The estimated nonpoint source loads are compared with current wastewater treatment facilities loads to determine the non‐point source contribution of total phosphorus (TP), total nitrogen (TN), and total suspended solids (TSS) on a monthly and annual time scale. An innovative calibration procedure is used to estimate the pollutant concentrations for different land uses based on available water quality data at the outlet. Results indicate that the pollutant concentrations are higher for the Las Vegas Valley than previous published values for semi‐arid and arid regions. The total TP and TN loads from nonpoint sources are approximately 15 percent and 4 percent, respectively, of the total load to the receiving water body, Lake Mead. The TP loads during wet periods approach the permitted loads from the wastewater treatment plants that discharge into Las Vegas Wash. In addition, the GIS model is used to track pollutant loads in the stream channels for one of the subwatersheds. This is useful for planning the location of Best Management Practices to control nonpoint pollutant loads.     

WATERSHED OPTIMIZATION OF AGRICULTURAL BEST MANAGEMENT PRACTICES: CONTINUOUS SIMULATION VERSUS DESIGN STORMS1

ABSTRACT: Nonpoint source (NPS) models and expert opinions are often used to prescribe best management practices (BMPs) for controlling NPS pollution. An optimization algorithm (e.g., a genetic algorithm, or GA) linked with a NPS model (e.g., Annualized AGricultural Nonpoint Source pollution model, or AnnAGNPS), can be used to more objectively prescribe BMPs and to optimize NPS pollution control measures by maximizing pollutant reduction and net monetary return from a watershed. Pollutant loads from design storms and annual loads from a continuous simulation can both be used for optimizing BMP schemes. However, which strategy results in a better solution (in terms of providing water quality protection) for a watershed is not clear. The specific objective of the study was to determine the differences in watershed pollutant loads, in an experimental watershed in Pennsylvania, resulting from optimization analyses performed using pollutant loads from a series of five 2‐yr 24‐hr storm events, a series of five 5‐yr 24‐hr storm events, and cumulative pollutant loads from a continuous simulation of five years of weather data. For each of these three different event alternatives, 100 near optimal solutions (BMP schemes) were generated. Sediment (Sed), sediment nitrogen (SedN), dissolved N (SolN), sediment organic carbon (SedOC), and sediment phosphorus (SedP) loads from a different five‐year period (an evaluation period) suggest that the optimal BMP schemes resulting from the use of annual cumulative pollutant loads from a continuous simulation of five years of weather data provide smaller cumulative NPS pollutant loads at the watershed outlet.     

Modeling Changes to Streamflow,Sediment, and Nutrient Loading from Land Use Changes Due to Potential Natural Gas Development

Natural gas development using hydraulic fracturing has many potential environmental impacts, but among the most certain is the land disturbance required to build the well pads and other infrastructure required to drill and extract the gas. We used the Soil and Water Assessment Tool (SWAT) model to investigate how natural gas development could impact streamflow and sediment, total nitrogen (TN), and total phosphorous (TP) loadings in the upper Delaware River Basin (DRB), a relatively undeveloped watershed of 7,950 km² that lies above the Marcellus Shale formation. If gas development was permitted, our projections show the DRB could experience development of over 600 well pads to extract natural gas at build out, which, with supporting infrastructure (roads, gathering pipelines), could convert over 5,000 ha from existing land uses in the study area. In subbasins with development activity we found sediment, TN, and TP yields could increase by an average of 15, 0.08, and 0.03 kg/ha/yr, respectively (an increase of 2, 3, and 15%, respectively) for each one percent of subbasin land area converted into natural gas infrastructure. At the study area outlet on the Delaware River at Port Jervis, New York, we found increases in the annual average streamflow and sediment, nitrogen, and phosphorus loads of up to 0.01, 0.2, 0.2, and 1%, respectively, for a rapid development year, and 0.08, 1.3, 2.0, and 11%, respectively, for the full development scenario. 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.     

DRY AND WET WEATHER FLOW NUTRIENT LOADS FROM A LOS ANGELES WATERSHED1

Effective watershed management requires an accurate assessment of the pollutant loads from the associated point and nonpoint sources. The importance of wet weather flow (WWF) pollutant loads is well known, but in semi‐arid regions where urbanization is significant the pollutant load in dry weather flow (DWF) may also be important. This research compares the relative contributions of potential contaminants discharged in DWF and WWF from the Ballona Creek Watershed in Los Angeles, California. Models to predict DWF and WWF loads of total suspended solids, biochemical oxygen demand, nitrate‐nitrogen, nitrite‐nitrogen, ammonia‐nitrogen, total Kjeldahl nitrogen, and total phosphorus from the Ballona Creek Watershed for six water years dating from 1991 to 1996 were developed. The contaminants studied were selected based on data availability and their potential importance in the degradation of Ballona Creek and Santa Monica Bay beneficial uses. Wet weather flow was found to contribute approximately 75 percent to 90 percent of the total annual flow volume discharged by the Ballona Creek Watershed. Pollutant loads are also predominantly due to WWF, but during the dry season, DWF is a more significant contributor. Wet weather flow accounts for 67 to 98 percent of the annual load of the constituents studied. During the dry season, however, the portion attributable to DWF increases to greater than 40 percent for all constituents except biochemical oxygen demand and total suspended solids. When individual catchments within the watershed are considered, the DWF pollutant load from the largest catchment is similar to the WWF pollutant load in two other major catchments. This research indicates WWF is the most significant source of nonpoint source pollution load on an annual basis, but management of the effects of the nonpoint source pollutant load should consider the seasonal importance of DWF.     

EFFECT OF WATERSHED SUBDIVISION ON SWAT FLOW,SEDIMENT, AND NUTRIENT PREDICTIONS1

ABSTRACT: The size, scale, and number of subwatersheds can affect a watershed modeling process and subsequent results. The objective of this study was to determine the appropriate level of subwatershed division for simulating flow, sediment, and nutrients over 30 years for four Iowa watersheds ranging in size from 2,000 to 18,000 km² with the Soil and Water Assessment Tool (SWAT) model. The results of the analysis indicated that variation in the total number of subwatersheds had very little effect on streamflow. However, the opposite result was found for sediment, nitrate, and inorganic P; the optimal threshold subwatershed sizes, relative to the total drainage area for each watershed, required to adequately predict these three indicators were found to be around 3, 2, and 5 percent, respectively. Decreasing the size of the subwatersheds below these threshold levels does not significantly affect the predicted levels of these environmental indicators. These threshold subwatershed sizes can be used to optimize input data preparation requirements for SWAT analyses of other watersheds, especially those within a similar size range. The fact that different thresholds emerged for the different indicators also indicates the need for SWAT users to assess which indicators should have the highest priority in their analyses.     

Influence of Sampling Frequency on Estimation of Annual Total Phosphorus and Total Suspended Solids Loads1

Abstract: The determination of sediment and nutrient loads is typically based on the collection and analysis of grab samples. The frequency and regularity of traditional sampling may not provide representation of constituent loading, particularly in systems with flashy hydrology. At two sites in the Little Bear River, Utah, continuous, high‐frequency turbidity was used with surrogate relationships to generate estimates of total phosphorus and total suspended solids concentrations, which were paired with discharge to estimate annual loads. The high frequency records were randomly subsampled to represent hourly, daily, weekly, and monthly sampling frequencies and to examine the effects of timing, and resulting annual load estimates were compared to the reference loads. Higher frequency sampling resulted in load estimates that better approximated the reference loads. The degree of bias was greater at the more hydrologically responsive site in the upper watershed, which required a higher sampling frequency than the lower watershed site to achieve the same level of accuracy in estimating the reference load. The hour of day and day of week of sampling impacted load estimation, depending on site and hydrologic conditions. The effects of sampling frequency on the determination of compliance with a water quality criterion were also examined. These techniques can be helpful in determining necessary sampling frequency to meet the objectives of a water quality monitoring program.     

Development of a SWAT Patch for Better Estimation of Sediment Yield in Steep Sloping Watersheds1

Abstract: The watershed scale Soil and Water Assessment Tool (SWAT) model divides watersheds into smaller subwatersheds for simulation of rainfall‐runoff and sediment loading at the field level and routing through stream networks. Typically, the SWAT model first needs to be calibrated and validated for accurate estimation through adjustment of sensitive input parameters (i.e., Curve Number values, USLE P, slope and slope‐length, and so on). However, in some instances, SWAT‐simulated results are greatly affected by the watershed delineation and Digital Elevation Models (DEM) cell size. In this study, the SWAT ArcView GIS Patch II was developed for steep sloping watersheds, and its performance was evaluated for various threshold values and DEM cell size scenarios when delineating subwatersheds using the SWAT model. The SWAT ArcView GIS Patch II was developed using the ArcView GIS Avenue program and Spatial Analyst libraries. The SWAT ArcView GIS Patch II improves upon the SWAT ArcView GIS Patch I because it reflects the topographic factor in calculating the field slope‐length of Hydrologic Response Units in the SWAT model. The simulated sediment value for 321 subwatersheds (watershed delineation threshold value of 25 ha) is greater than that for 43 subwatersheds (watershed delineation threshold value of 200 ha) by 201% without applying the SWAT ArcView GIS Patch II. However, when the SWAT ArcView GIS Patch II was applied, the difference in simulated sediment yield decreases for the same scenario (i.e., difference in simulated sediment with 321 subwatersheds and 43 subwatersheds) was 12%. The simulated sediment value for DEM cell size of 50 m is greater than that for DEM cell size of 10 m by 19.8% without the SWAT ArcView GIS Patch II. However, the difference becomes smaller (3.4% difference) between 50 and 10 m with the SWAT ArcView GIS Patch II for the DEM scenarios. As shown in this study, the SWAT ArcView GIS Patch II can reduce differences in simulated sediment values for various watershed delineation and DEM cell size scenarios. Without the SWAT ArcView GIS Patch II, variations in the SWAT‐simulated results using various watershed delineation and DEM cell size scenarios could be greater than those from input parameter calibration. Thus, the results obtained in this study show that the SWAT ArcView GIS Patch II should be used when simulating hydrology and sediment yield for steep sloping watersheds (especially if average slope of the subwatershed is >25%) for more accurate simulation of hydrology and sediment using the SWAT model. The SWAT ArcView GIS Patch II is available at http://www.EnvSys.co.kr/~swat for free download.     

SPATIAL AND TEMPORAL VARIABILITY IN SUSPENDED SEDIMENT YIELDS OF NORTH ISLAND BASINS,NEW ZEALAND1

ABSTRACT: Specific annual suspended sediment yields and their standard deviations are presented for 47 basins of North Island, New Zealand. Most of the variance in yields is explained by catchment mean rainfall. Rivers with similar flow range have similar suspended sediment concentration ratings, independent of differing watershed lithology and regolith, except for six basins having an abundance of soft fine sediments. Prediction equations for yield and its standard deviation are derived for four essentially arbitrary regions. AU feature rainfall as the independent variable. Differences between regions may owe to variations in intensity, frequency, and duration patterns of storms and, in one area, to bed material size as well. The temporal distribution of annual yields from a basin m be modeled by a two-parameter lognormal function: the prediction equations above may be used to evaluate this function at a site for which suspended sediment data are unavailable.     

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.