A Comparison of Statistical Approaches for Predicting Stream Temperatures Across Heterogeneous Landscapes1

Abstract: Estimating stream temperatures across broad spatial extents is important for regional conservation of running waters. Although statistical models can be useful in this endeavor, little information exists to aid in the selection of a particular statistical approach. Our objective was to compare the accuracy of ordinary least‐squares multiple linear regression, generalized additive modeling, ordinary kriging, and linear mixed modeling (LMM) using July mean stream temperatures in Michigan and Wisconsin. Although LMM using low‐rank thin‐plate smoothing splines to measure the spatial autocorrelation in stream temperatures was the most accurate modeling approach; overall, there were only slight differences in prediction accuracy among the evaluated approaches. This suggests that managers and researchers can select a stream temperature modeling approach that meets their level of expertise without sacrificing substantial amounts of prediction accuracy. The most accurate models for Michigan and Wisconsin had root mean square errors of 2.0‐2.3°C, suggesting that only relatively coarse predictions can be produced from landscape‐based statistical models at regional scales. Explaining substantially more variability in stream temperatures likely will require the collection of finer‐scale hydrologic and physiographic data, which may be cost prohibitive for monitoring and assessing stream temperatures at regional scales.     

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.     

INFLUENCES OF SOIL DATASET RESOLUTION ON HYDROLOGIC MODELING1

ABSTRACT: Soils represent a fundamental abiotic parameter in defining the characteristics of an ecosystem. The U.S. Department of Agriculture (USDA) Natural Resources Conservation Service (NRCS) produces the most detailed digital spatial soil datasets that are publicly available. The Soil Survey Geographic (SSURGO) database contains basic attributes for the continuous coverage of soils across the United States. In its standard format, the SSURGO database is incompatible for use within the ArcView Soil and Water Assessment Tool (SWAT). A modified version of the State Soil and Geographic (STATSGO) database is the template soils dataset used by ArcView SWAT. This paper presents the methodology and development of a SSURGO database preprocessor extension for the ArcView SWAT model. A case study for the Upper Sabinal River Watershed near Uvalde, Texas, is given. Results indicate that hydro‐logic output parameter differences occur when comparing the STATSGO and SSURGO database information in the ArcView SWAT model under identical modeling conditions. Specifically, the SSURGO model produced a greater daily mean water yield with evapotranspiration and surface runoff being found consistently lower across the watershed. The most likely causes assigned to this phenomenon were higher percolation and resulting ground water return flow values due to significantly larger saturated hydraulic conductivity values associated with the SSURGO 2.x database.     

EVALUATION OF THE APPLICABILITY OF THE SWAT MODEL FOR COASTAL WATERSHEDS IN SOUTHEASTERN LOUISIANA1

ABSTRACT: The Soil and Water Assessment Tool (SWAT) has been used for hydrologic analyses at various watershed scales. However, little is known about the model's performance in coastal watersheds. In this study SWAT was evaluated for its applicability in three Louisiana coastal watersheds: the Amite, Tickfaw, and Tangipahoa River watersheds. The model was calibrated with daily discharge from 1976 to 1977 and validated from 1979 to 1999 for the Amite and Tangipahoa and with daily discharge from 1979 to 1989 for the Tickfaw. Deviation of mean discharge and the Nash‐Sutcliffe model efficiency were used to evaluate model behavior. The study found that Manning's roughness coefficient for the main channel, SCS curve number, and soil evaporation compensation factor were the most sensitive parameters for these coastal watersheds. The Manning's roughness coefficient showed the greatest effect on the response time of surface runoff, suggesting the critical role of channel routing in hydrologic modeling for lowland watersheds. The SWAT model demonstrated an excellent performance, with Nash‐Sutcliffe efficiencies of 0.935, 0.940, and 0.960 for calibrations of the Amite, Tickfaw, and Tangipahoa watersheds, respectively, and of 0.851, 0.811, and 0.867 for validations. The modeling results demonstrate that SWAT is capable of simulating hydrologic processes for medium scale to large scale coastal lowland watersheds in Louisiana.     

PRISM vs. CFSR Precipitation Data Effects on Calibration and Validation of SWAT Models

Precipitation is one of the most important drivers in watershed models. Our objective was to compare two sources of interpolated precipitation data in terms of their effect on calibration and validation of two Soil and Water Assessment Tool (SWAT) models. One model was a suburban watershed in metropolitan Atlanta, Georgia. The precipitation sources were Parameter‐elevation Relationships on Independent Slopes Model (PRISM) data on a 4‐km grid and climate forecast system reanalysis (CFSR) data on a 38‐km grid. The PRISM data resulted in a better fit to the calibration data (Nash Sutcliffe efficiency [NSE] = 0.64, Kling‐Gupta efficiency [KGE] = 0.74, p‐factor = 0.84, and r‐factor = 0.43) than the CFSR data (NSE = 0.47, KGE = 0.53, p‐factor = 0.67, and r‐factor = 0.39). Validation results were similar. Sensitive parameters were similar in both the PRISM and CFSR models, but fitted values indicated more rapid groundwater flow to the streams with the PRISM data. The same comparison was made in the Big Creek watershed located approximately 1,000 km away, in central Louisiana. Results were similar with a more responsive groundwater system indicating PRISM data may produce better predictions of streamflow because of a more accurate estimate of rainfall within a watershed or because of a denser grid. Our study implies PRISM is providing a better estimate than CFSR of precipitation within a watershed when rain gauge data are not available, resulting in more accurate simulations of streamflows at the watershed outlet. 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.     

Using GIS to Delineate Headwater Stream Origins in the Appalachian Coalfields of Kentucky

Headwater streams have a significant nexus or physical, chemical, and/or biological connection to downstream reaches. Generally, defined as 1st‐3rd order with ephemeral, intermittent, or perennial flow regimes, these streams account for a substantial portion of the total stream network particularly in mountainous terrain. Due to their often remote locations, small size, and large numbers, conducting field inventories of headwater streams is challenging. A means of estimating headwater stream location and extent according to flow regime type using publicly available spatial data is needed to simplify this complex process. Using field‐collected headwater point of origin data from three control watersheds, streams were characterized according to a set of spatial parameters related to topography, geology, and soils. These parameters were (1) compared to field‐collected point of origin data listed in three nearby Jurisdictional Determinations, (2) used to develop a geographic information system (GIS)‐based stream network for identifying ephemeral, intermittent, and perennial streams, and (3) applied to a larger watershed and compared to values obtained using the high‐resolution National Hydrography Dataset (NHD). The parameters drainage area and local valley slope were the most reliable predictors of flow regime type. Results showed the high‐resolution NHD identified no ephemeral streams and 9 and 65% fewer intermittent and perennial streams, respectively, than the GIS model.     

WATER RESOURCE MANAGEMENT AND DROUGHT FORECASTS IN SOUTH CAROLINA1

ABSTRACT: Despite potential benefits for resource planning, community water systems managers have not used seasonal climate forecasts extensively. Obstacles to forecast use include a lack of awareness of their existence, distrust of their accuracy, perceived irrelevance to management decisions, and competition from other technological innovations. In this paper, ways in which seasonal forecasts might be extended to address more directly some concerns of South Carolina community water systems managers are explored. From May 1998 through August 2002, this group experienced drought conditions that threatened water quality and supply and required restrictions on water consumption. Methods for incorporating long lead forecasts with joint probabilities of monthly temperature and precipitation to produce drought forecasts are demonstrated. When tailored to specific places, such forecasts show the likelihood of exceeding drought thresholds that would trigger water use restrictions. The methods illustrate how long lead forecasts can be extended and customized into secondary products that address issues of greater relevance to water resource managers.     

Reducing Fertilizer‐Nitrogen Losses from Rowcrop Landscapes: Insights and Implications from a Spatially Explicit Watershed Model

We present conceptual and quantitative models that predict changes in fertilizer‐derived nitrogen delivery from rowcrop landscapes caused by agricultural conservation efforts implemented to reduce nutrient inputs and transport and increase nutrient retention in the landscape. To evaluate the relative importance of changes in the sources, transport, and sinks of fertilizer‐derived nitrogen across a region, we use the spatially explicit SPAtially Referenced Regression On Watershed attributes watershed model to map the distribution, at the small watershed scale within the Upper Mississippi‐Ohio River Basin (UMORB), of: (1) fertilizer inputs; (2) nutrient attenuation during delivery of those inputs to the UMORB outlet; and (3) nitrogen export from the UMORB outlet. Comparing these spatial distributions suggests that the amount of fertilizer input and degree of nutrient attenuation are both important in determining the extent of nitrogen export. From a management perspective, this means that agricultural conservation efforts to reduce nitrogen export would benefit by: (1) expanding their focus to include activities that restore and enhance nutrient processing in these highly altered landscapes; and (2) targeting specific types of best management practices to watersheds where they will be most valuable. Doing so successfully may result in a shift in current approaches to conservation planning, outreach, and funding.     

VIRGINIA USA WATER QUALITY, 1978 TO 1995: REGIONAL INTERPRETATION1

ABSTRACT: Nine surface water‐quality variables were analyzed for trend at 180 Virginia locations over the 1978 to 1995 period. Median values and seasonal Kendall's tau, a trend indicator statistic, were generated for dissolved oxygen saturation (DO), biochemical oxygen demand (BOD), pH (PH), total residue (TR), nonfilterable residue (NFR), nitrate‐nitrite nitrogen (NN), total Kjeldahl nitrogen (TKN), total phosphorus (TP), and fecal coliform (FC) at each location. Each location was assigned to one of four physiographic regions, and mean state and regional medians and taus were calculated. Widespread BOD and NFR improvements were detected and FC improvements occurred in the state's western regions. TR and TKN exhibited predominantly increasing trends at locations throughout the state. BOD, TKN, NFR, and TR medians were higher at coastal locations than in other regions. NN, TKN, and TR exhibited predominantly increasing trends in regions with high median concentrations, while declining trends predominated in regions with relatively high BOD, FC, and NFR medians. Appalachian locations exhibited the greatest regional water‐quality improvements for BOD, FC, NFR, and TKN. Factors responsible for regional differences appear to include geology, land use, and landscape features; these factors vary regionally.     

Targeted Application of Seasonal Load Duration Curves Using Multivariate Analysis in Two Watersheds Flowing Into Lake Houston1

Teague, Aarin, Philip B. Bedient, and Birnur Guven, 2011. Targeted Application of Seasonal Load Duration Curves Using Multivariate Analysis in Two Watersheds Flowing Into Lake Houston. Journal of the American Water Resources Association (JAWRA) 47(3):620‐634. DOI: 10.1111/j.1752‐1688.2011.00529.x Abstract: Water quality is a problem in Lake Houston, the primary source of drinking water for the City of Houston, Texas, due to pollutant loads coming from the influent watersheds, including Spring Creek and Cypress Creek. Statistical analysis of the historic water quality data was developed to understand the source characterization and seasonality of the watershed. Multivariate analysis including principal component, cluster, and discriminant analysis provided a custom seasonal assessment of the watersheds so that loading curves may be targeted for season specific pollutant source characterization. The load duration curves have been analyzed using data collected by the U.S. Geologic Survey with corresponding City of Houston water quality data at the sites to characterize the behavior of the pollutant sources and watersheds. Custom seasons were determined for Spring Creek and Cypress Creek watersheds and pollutant source characterization compared between the seasons and watersheds.     

Locating Existing Best Management Practices Within a Watershed: The Value of Multiple Methods

There is an increasing need to document the impacts of conservation‐related best management practices (BMPs) on water quality within a watershed. However, this impact analysis depends upon accurate geospatial locations of existing practices, which are difficult to obtain. This study demonstrates and evaluates three different methods for obtaining geospatial information for BMPs. This study was focused on the Eagle Creek Watershed, a mixed use watershed in central Indiana. We obtained geospatial information for BMPs through government records, producer interviews, and remote‐sensing aerial photo interpretation. Aerial photos were also used to validate the government records and producer interviews. This study shows the variation in results obtained from the three sources of information as well as the benefits and drawbacks of each method. Using only one method for obtaining BMP information can be incomplete, and this study demonstrates how multiple methods can be used for the most accurate picture.     

CLIMATE CHANGE AND WATERSHED PLANNING IN WASHINGTON STATE1

ABSTRACT: This paper draws on interviews with Washington State Watershed Planning Leads (Planning Leads) and interactions with local watershed planning units to identify factors that may influence the inclusion of climate change in watershed planning efforts in Washington State. These factors include the interest of individual planning unit members in climate change; Planning Lead familiarity with climate impacts; the influence of trust, leadership, and “genetic knowledge” on planning units; and perceptions of strategic gain. The research also identifies aspects of the planning process that may create opportunities for addressing climate impacts in future planning. These aspects include continuation of watershed planning units after plans are developed; commitment to updating watershed plans; recognition of climate impacts in planning documentation; dedicated incentive funding; and the availability of hydrologic modeling tools for assessing hydrologic impacts. Additional types of technical assistance that could support integration of climate impacts are also identified. It is hoped that the insight provided by this analysis will help individuals involved in stakeholder‐based watershed planning recognize the various dynamics potentially affecting the inclusion of climate change in watershed planning and in doing so, contribute to the development of planning approaches and tools that will support local efforts to adapt to climate impacts.     

WATERSHED SCALE INVENTORY OF EXISTING RIPARIAN BUFFERS IN NORTHEAST MISSOURI USING GIS1

An observational study was conducted at the watershed scale using land cover (vegetation) data to assess the absence or presence of riparian buffers in three northeastern Missouri watersheds. Forests and grasslands lying within a 61 m (200 ft) parallel band directly adjacent to streams were considered “buffers” for improving or protecting water quality and were characterized according to their length, width, and vegetation type. Results indicated that riparian buffers were abundant throughout the watersheds but were typically narrow along first‐order and second‐order streams; in many cases they may not have been wide enough to provide adequate stream protection. At least 90 percent of all streams had buffer vegetation immediately adjacent to the streambanks, but as few as 31 percent of first‐order streams had buffers extending to 61 m from the stream on at least one side. On‐site evaluations are needed to determine the condition of these forests and grasslands and their ability to process nonpoint source pollutants. The results will be useful for providing natural resource managers with knowledge of current watershed conditions as well as in identifying specific locations for future conservation efforts within each watershed.     

A Comparison of DEM‐Based Indexes for Targeting the Placement of Vegetative Buffers in Agricultural Watersheds

Targeted placement of vegetative buffers may increase their effectiveness for improving water quality in agricultural watersheds. The use of digital elevation models (DEMs) enables precise mapping of runoff pathways for identifying where greater runoff loads can be intercepted and treated with buffers. Five different DEM‐based targeting indexes were compared and contrasted for the degree to which they identify similar locations in watersheds: Flow Accumulation [S.K. Jenson and J.O. Domingue (1988). Photogrammetric Engineering and Remote Sensing 54:1593], Wetness Index [I.D. Moore, R.B. Grayson, and A.R. Ladson (1991). Hydrological Processes 5:3], Topographic Index [M.T. Walter, T.S. Steenhuis, V.K. Mehta, D. Thongs, M. Zion, and E. Schneiderman (2002). Hydrological Processes 16:2041], and the Water Inflow and Sediment Retention Indexes [M.G. Dosskey, Z. Qiu, M.J. Helmers, and D.E. Eisenhauer (2011b). Journal of Soil and Water Conservation 66:362]. The indexes were applied in two different watersheds, one in New Jersey and one in Missouri. Results showed that they all tend to target similar locations in both watersheds which traces to the importance of larger contributing area to the rankings by each index. Disagreement among indexes traces to other variables which enable more accurate targeting under particular hydrologic circumstances. Effective use of these indexes poses special challenges, including selecting an index that better describes the hydrologic circumstances in a watershed and is simple enough to use, ensuring the accuracy of the DEM, and determining a maximum index value for the appropriateness of vegetative buffers. When properly applied, each index can provide a standardized basis and effective spatial resolution for targeting buffer placement in watersheds.     

The Regionalization of National‐Scale SPARROW Models for Stream Nutrients1

Schwarz, Gregory E., Richard B. Alexander, Richard A. Smith, and Stephen D. Preston, 2011. The Regionalization of National‐Scale SPARROW Models for Stream Nutrients. Journal of the American Water Resources Association (JAWRA) 47(5):1151‐1172. DOI: 10.1111/j.1752‐1688.2011.00581.x Abstract: This analysis modifies the parsimonious specification of recently published total nitrogen (TN) and total phosphorus (TP) national‐scale SPAtially Referenced Regressions On Watershed attributes models to allow each model coefficient to vary geographically among three major river basins of the conterminous United States. Regionalization of the national models reduces the standard errors in the prediction of TN and TP loads, expressed as a percentage of the predicted load, by about 6 and 7%. We develop and apply a method for combining national‐scale and regional‐scale information to estimate a hybrid model that imposes cross‐region constraints that limit regional variation in model coefficients, effectively reducing the number of free model parameters as compared to a collection of independent regional models. The hybrid TN and TP regional models have improved model fit relative to the respective national models, reducing the standard error in the prediction of loads, expressed as a percentage of load, by about 5 and 4%. Only 19% of the TN hybrid model coefficients and just 2% of the TP hybrid model coefficients show evidence of substantial regional specificity (more than ±100% deviation from the national model estimate). The hybrid models have much greater precision in the estimated coefficients than do the unconstrained regional models, demonstrating the efficacy of pooling information across regions to improve regional models.     

Effect of the Spatial Variability of Land Use,Soil Type,and Precipitation on Streamflows in Small Watersheds1

Abstract: The spatial variability of the data used in models includes the spatial discretization of the system into subsystems, the data resolution, and the spatial distribution of hydrologic features and parameters. In this study, we investigate the effect of the spatial distribution of land use, soil type, and precipitation on the simulated flows at the outlet of “small watersheds” (i.e., watersheds with times of concentration shorter than the model computational time step). The Soil and Water Assessment Tool model was used to estimate runoff and hydrographs. Different representations of the spatial data resulted in comparable model performances and even the use of uniform land use and soil type maps, instead of spatially distributed, was not noticeable. It was found that, although spatially distributed data help understand the characteristics of the watershed and provide valuable information to distributed hydrologic models, when the watershed is small, realistic representations of the spatial data do not necessarily improve the model performance. The results obtained from this study provide insights on the relevance of taking into account the spatial distribution of land use, soil type, and precipitation when modeling small watersheds.