The Great Lakes Hydrography Dataset: Consistent,Binational Watersheds for the Laurentian Great Lakes Basin

Ecosystem‐based management of the Laurentian Great Lakes, which spans both the United States and Canada, is hampered by the lack of consistent binational watersheds for the entire Basin. Using comparable data sources and consistent methods, we developed spatially equivalent watershed boundaries for the binational extent of the Basin to create the Great Lakes Hydrography Dataset (GLHD). The GLHD consists of 5,589 watersheds for the entire Basin, covering a total area of approximately 547,967 km², or about twice the 247,003 km² surface water area of the Great Lakes. The GLHD improves upon existing watershed efforts by delineating watersheds for the entire Basin using consistent methods; enhancing the precision of watershed delineation using recently developed flow direction grids that have been hydrologically enforced and vetted by provincial and federal water resource agencies; and increasing the accuracy of watershed boundaries by enforcing embayments, delineating watersheds on islands, and delineating watersheds for all tributaries draining to connecting channels. In addition, the GLHD is packaged in a publically available geodatabase that includes synthetic stream networks, reach catchments, watershed boundaries, a broad set of attribute data for each tributary, and metadata documenting methodology. The GLHD provides a common set of watersheds and associated hydrography data for the Basin that will enhance binational efforts to protect and restore the Great Lakes.     

Applications of Numerical Sediment Quality Targets for Assessing Sediment Quality Conditions in a US Great Lakes Area of Concern

Contaminated sediments are receiving increasing recognition around the world, leading to the development of various sediment quality indicators for assessment, management, remediation, and restoration efforts. Sediment chemistry represents an important indicator of ecosystem health, with the concentrations of contaminants of potential concern (COPCs) providing measurable characteristics for this indicator. The St. Louis River Area of Concern (AOC), located in the western arm of Lake Superior, provides a case study for how numerical sediment quality targets (SQTs) for the protection of sediment-dwelling organisms can be used to support the interpretation of sediment chemistry data. Two types of SQTs have been established for 33 COPCs in the St. Louis River AOC. The Level I SQTs define the concentrations of contaminants below which sediment toxicity is unlikely to occur, whereas the Level II SQTs represent the concentrations that, if exceeded, are likely to be associated with sediment toxicity. The numerical SQTs provide useful tools for making sediment management decisions, especially when considered as part of a weight-of-evidence approach that includes other sediment quality indicators, such as sediment contaminant chemistry and geochemical characteristics, sediment toxicity, and benthic macroinvertebrate community structure. The recommended applications of using the numerical SQTs in the St. Louis River AOC include: designing monitoring programs, interpreting sediment chemistry data, conducting ecological risk assessments, and developing site-specific sediment quality remediation targets for small, simple sites where adverse biological effects are likely. Other jurisdictions may benefit from using these recommended applications in their own sediment quality programs.     

Nutrient Inputs to the Laurentian Great Lakes by Source and Watershed Estimated Using SPARROW Watershed Models1

Robertson, Dale M. and David A. Saad, 2011. Nutrient Inputs to the Laurentian Great Lakes by Source and Watershed Estimated Using SPARROW Watershed Models. Journal of the American Water Resources Association (JAWRA) 47(5):1011‐1033. DOI: 10.1111/j.1752‐1688.2011.00574.x Abstract: Nutrient input to the Laurentian Great Lakes continues to cause problems with eutrophication. To reduce the extent and severity of these problems, target nutrient loads were established and Total Maximum Daily Loads are being developed for many tributaries. Without detailed loading information it is difficult to determine if the targets are being met and how to prioritize rehabilitation efforts. To help address these issues, SPAtially Referenced Regressions On Watershed attributes (SPARROW) models were developed for estimating loads and sources of phosphorus (P) and nitrogen (N) from the United States (U.S.) portion of the Great Lakes, Upper Mississippi, Ohio, and Red River Basins. Results indicated that recent U.S. loadings to Lakes Michigan and Ontario are similar to those in the 1980s, whereas loadings to Lakes Superior, Huron, and Erie decreased. Highest loads were from tributaries with the largest watersheds, whereas highest yields were from areas with intense agriculture and large point sources of nutrients. Tributaries were ranked based on their relative loads and yields to each lake. Input from agricultural areas was a significant source of nutrients, contributing ～33‐44% of the P and ～33‐58% of the N, except for areas around Superior with little agriculture. Point sources were also significant, contributing ～14‐44% of the P and 13‐34% of the N. Watersheds around Lake Erie contributed nutrients at the highest rate (similar to intensively farmed areas in the Midwest) because they have the largest nutrient inputs and highest delivery ratio.     

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.     

Assessing Risk in Operational Decisions Using Great Lakes Probabilistic Water Level Forecasts

/ A method adapted from the National Weather Service's Extended Streamflow Prediction technique is applied retrospectively to three Great Lakes case studies to show how risk assessment using probabilistic monthly water level forecasts could have contributed to the decision-mak-ing process. The first case study examines the 1985 International Joint Commission (IJC) decision to store water in Lake Superior to reduce high levels on the downstream lakes. Probabilistic forecasts are generated for Lake Superior and Lakes Michigan-Huron and used with riparian inundation value functions to assess the relative impacts of the IJC's decision on riparian interests for both lakes. The second case study evaluates the risk of flooding at Milwaukee, Wisconsin, and the need to implement flood-control projects if Lake Michigan levels were to continue to rise above the October 1986 record. The third case study quantifies the risks of impaired municipal water works operation during the 1964-1965 period of extreme low water levels on Lakes Huron, St. Clair, Erie, and Ontario. Further refinements and other potential applications of the probabilistic forecast technique are discussed.KEY WORDS: Great Lakes; Water levels; Forecasting; Risk; Decision making     

Direct and Indirect Effects of Forest Harvesting on Sediment Yield in Forested Watersheds of the United States

Managed forests generally produce high water quality, but degradation is possible via sedimentation if proper management is not implemented during forest harvesting. To mitigate harvesting effects on total watershed sediment yield, it is necessary to understand all processes that contribute to these effects. Forest harvesting best management practices (BMPs) focus almost exclusively on overland sediment sources, whereas in‐and‐near stream sources go unaddressed although they can contribute substantially to sediment yield. Thus, we propose a new framework to classify forest harvesting effects on stream sediment yield according to their direct and indirect processes. Direct effects are those caused by erosion and sediment delivery to surface water from overland sources (e.g., forest roads). Indirect effects are those caused by a shift in hydrologic processes due to tree removal that accounts for increases in subsurface and surface flows to the stream such that alterations in water quality are not predicated upon overland sediment delivery to the stream, but rather in‐stream processes. Although the direct/indirect distinction is often implicit in forest hydrology studies, we have formalized it as a conceptual model to help identify primary drivers of sediment yield after forest harvesting in different landscapes. Based on a literature review, we identify drivers of these effects in five regions of the United States, discuss current forest management BMPs, and identify research needs.     

Integration of SWAT and HSPF for Simulation of Sediment Sources in Legacy Sediment‐Impacted Agricultural Watersheds

A total maximum daily load for the Chesapeake Bay requires reduction in pollutant load from sources within the Bay watersheds. The Conestoga River watershed has been identified as a major source of sediment load to the Bay. Upland loads of sediment from agriculture are a concern; however, a large proportion of the sediment load in the Conestoga River has been linked to scour of legacy sediment associated with historic millpond sites. Clarifying this distinction and identifying specific segments associated with upland vs. channel sources has important implications for future management. In order to address this important question, we combined the strengths of two widely accepted watershed management models — Soil and Water Assessment Tool (SWAT) for upland agricultural processes, and Hydrologic Simulation Program FORTRAN (HSPF) for instream fate and transport — to create a novel linked modeling system to predict sediment loading from critical sources in the watershed including upland and channel sources, and to aid in targeted implementation of management practices. The model indicates approximately 66% of the total sediment load is derived from instream sources, in agreement with other studies in the region and can be used to support identification of these channel source segments vs. upland source segments, further improving targeted management. The innovated linked SWAT‐HSPF model implemented in this study is useful for other watersheds where both upland agriculture and instream processes are important sources of sediment load.     

Integrated Measures of Anthropogenic Stress in the U.S. Great Lakes Basin

Integrated, quantitative expressions of anthropogenic stress over large geographic regions can be valuable tools in environmental research and management. Despite the fundamental appeal of a regional approach, development of regional stress measures remains one of the most important current challenges in environmental science. Using publicly available, pre-existing spatial datasets, we developed a geographic information system database of 86 variables related to five classes of anthropogenic stress in the U.S. Great Lakes basin: agriculture, atmospheric deposition, human population, land cover, and point source pollution. The original variables were quantified by a variety of data types over a broad range of spatial and classification resolutions. We summarized the original data for 762 watershed-based units that comprise the U.S. portion of the basin and then used principal components analysis to develop overall stress measures within each stress category. We developed a cumulative stress index by combining the first principal component from each of the five stress categories. Maps of the stress measures illustrate strong spatial patterns across the basin, with the greatest amount of stress occurring on the western shore of Lake Michigan, southwest Lake Erie, and southeastern Lake Ontario. We found strong relationships between the stress measures and characteristics of bird communities, fish communities, and water chemistry measurements from the coastal region. The stress measures are taken to represent the major threats to coastal ecosystems in the U.S. Great Lakes. Such regional-scale efforts are critical for understanding relationships between human disturbance and ecosystem response, and can be used to guide environmental decision-making at both regional and local scales.     

Impacts of Land‐Cover Change on Suspended Sediment Transport in Two Agricultural Watersheds1

Schilling, Keith E., Thomas M. Isenhart, Jason A. Palmer, Calvin F. Wolter, and Jean Spooner, 2011. Impacts of Land‐Cover Change on Suspended Sediment Transport in Two Agricultural Watersheds. Journal of the American Water Resources Association (JAWRA) 47(4):672‐686. DOI: 10.1111/j.1752‐1688.2011.00533.x Abstract: Suspended sediment is a major water quality problem, yet few monitoring studies have been of sufficient scale and duration to assess the effectiveness of land‐use change or conservation practice implementation at a watershed scale. Daily discharge and suspended sediment export from two 5,000‐ha watersheds in central Iowa were monitored over a 10‐year period (water years 1996‐2005). In Walnut Creek watershed, a large portion of land was converted from row crop to native prairie, whereas in Squaw Creek land use remained predominantly row crop agriculture. Suspended sediment loads were similar in both watersheds, exhibiting flashy behavior typical of incised channels. Modeling suggested that expected total soil erosion in Walnut Creek should have been reduced 46% relative to Squaw Creek due to changes in land use, yet measured suspended sediment loads showed no significant differences. Stream mapping indicated that Walnut Creek had three times more eroding streambank lengths than did Squaw Creek suggesting that streambank erosion dominated sediment sources in Walnut Creek and sheet and rill sources dominated sediment sources in Squaw Creek. Our results demonstrate that an accounting of all sources of sediment erosion and delivery is needed to characterize sediment reductions in watershed projects combined with long‐term, intensive monitoring and modeling to account for possible lag times in the manifestation of the benefits of conservation practices on water quality.     

Assessing Naturalness in Northern Great Lakes Forests Based on Historical Land-Cover and Vegetation Changes

The concept of naturalness was developed to assess to what degree landscapes represent a natural state. Protected areas are often regarded as the remnants of untouched landscapes although many landscapes commonly perceived as pristine have a long history of human impact. Here, we introduced a historical perspective into the concept of naturalness and the analysis of the effectiveness of protected areas by analyzing historical trajectories in land-cover and forest communities for the Pictured Rocks National Lakeshore on Michigan’s Upper Peninsula (USA). Distribution of land-cover and forest community types was reconstructed for pre-settlement time (around 1850), the height of agricultural expansion (1928), and modern conditions (2000). Naturalness of the landscape was assessed by analyzing similarity between pre-settlement and current conditions and by assessing landscape continuity (1850–1928–2000). We compared changes in the strictly protected park core zone with those in the inland buffer zone with ongoing sustainable logging, and a not protected area adjacent to the park. Forest was the dominant land-cover type over the entire study period. We detected a gradient in land-cover continuity from the core zone (81 % continuity) to the inland buffer zone (74 %) and the area outside the park (66 %). Northern hardwood was the dominating forest type in all time points with high continuity (76 %). In contrast, pine forests show a more dynamic pattern with more than 50 % of the initial forests switching to non-forest or early succession forest types by 1928. More than half of the study area was considered as “natural virgin” (no changes in land-cover and forest community type) with a higher portion within the park than in the adjacent area. In contrast, areas with low naturalness are more abundant outside the park. Our study demonstrates the value of integrating historical information into naturalness assessments and the results provide useful information for future park management. More broadly speaking, our study advances research on the effectiveness of protected areas, by going beyond simple measures of averted deforestation, and introducing approaches to directly measure naturalness.     

Agricultural Land Use Effects on Sediment Loading and Fish Assemblages in Two Minnesota (USA) Watersheds

We examined the relationship between water quality and fish communities within two agricultural areas using a computer simulation model. Our analyses focused on a coolwater stream, Wells Creek in southeastern Minnesota, and a warmwater stream, the Chippewa River in western Minnesota. We used the Agricultural Drainage and Pesticide Transport (ADAPT) model in relation to land use to calculate instream suspended sediment concentrations using estimates of sediment delivery, runoff, baseflow and streambank erosion, and quantified the effects of suspended sediment exposure on fish communities. We predicted the effects of agricultural practices on stream fish communities under several possible land use scenarios, with reference to current conditions. Land use changes led to reductions in sediment loading of up to 84% in Wells Creek and 49% in the Chippewa River. The reduction in sediment loading across scenarios may be directly related to a reduction in runoff by about 35% in both study areas. We found a 98% decrease in lethal concentrations of suspended sediment on fish in Wells Creek with an increase in conservation tillage, riparian buffers, and permanent vegetative cover. However, the effects of suspended sediment did not significantly decrease in the Chippewa River. This difference between study areas was likely due to differences in tolerance to suspended sediment between coolwater and warmwater fish communities and differences in topography, runoff and bank erosion between the two streams. The Minnesota Cooperative Fish and Wildlife Research Unit is jointly sponsored by the US Geological Survey, the University of Minnesota, the Minnesota Department of Natural Resources, and the Wildlife Management Institute.     

A VSA-Based Strategy for Placing Conservation Buffers in Agricultural Watersheds

Conservation buffers have the potential to reduce agricultural nonpoint source pollution and improve terrestrial wildlife habitat, landscape biodiversity, flood control, recreation, and aesthetics. Conservation buffers, streamside areas and riparian wetlands are being used or have been proposed to control agricultural nonpoint source pollution. This paper proposes an innovative strategy for placing conservation buffers based on the variable source area (VSA) hydrology. VSAs are small, variable but predictable portion of a watershed that regularly contributes to runoff generation. The VSA-based strategy involves the following three steps: first, identifying VSAs in landscapes based on natural characteristics such as hydrology, land use/cover, topography and soils; second, targeting areas within VSAs for conservation buffers; third, refining the size and location of conservation buffers based on other factors such as weather, environmental objectives, available funding and other best management practices. Building conservation buffers in VSAs allows agricultural runoff to more uniformly enter buffers and stay there longer, which increases the buffers capacity to remove sediments and nutrients. A field-scale example is presented to demonstrate the effectiveness and cost-effectiveness of the within-VSA conservation buffer scenario relative to a typical edge-of-field buffer scenario. The results enhance the understanding of hydrological processes and interactions between agricultural lands and conservation buffers in agricultural landscapes, and provide practical guidance for land resource managers and conservationists who use conservation buffers to improve water quality and amenity values of agricultural landscape.     

Quantifying Targets for Rehabilitating Degraded Areas of the Great Lakes

/ One attempt to quantify targets for rehabilitating degraded aquatic ecosystems has been through a United States-Canada program to develop and implement comprehensive remedial action plans (RAPs) to restore beneficial uses in 42 Great Lakes Areas of Concern. The International Joint Commission has facilitated agreement on listing/delisting guidelines for determining when use impairments exist in areas of concern and when uses have been restored, while federal/state/provincial governments and local stakeholders have provided leadership in establishing quantitative targets for restoring uses and in determining how to achieve them. The listing/delisting guidelines have been instrumental in helping reach agreement on problem definition (lack of agreement on problem definition has historically been used as a reason to delay action) and reaching agreement on quantitative targets for restoring uses. Quantitative, ecosystem-based targets are being used to drive the RAP process, help organizations pursue a common mission of restoring uses, and help achieve greater accountability. As a priority, the target-setting process must also recognize the importance of establishing both short- and long-term milestones in order to measure and celebrate incremental progress in restoring uses.KEY WORDS: Use impairments; Restoring uses; Quantitative targets     

Landscape Planning for Agricultural Nonpoint Source Pollution Reduction III: Assessing Phosphorus and Sediment Reduction Potential

Riparian buffers have the potential to improve stream water quality in agricultural landscapes. This potential may vary in response to landscape characteristics such as soils, topography, land use, and human activities, including legacies of historical land management. We built a predictive model to estimate the sediment and phosphorus load reduction that should be achievable following the implementation of riparian buffers; then we estimated load reduction potential for a set of 1598 watersheds (average 54 km²) in Wisconsin. Our results indicate that land cover is generally the most important driver of constituent loads in Wisconsin streams, but its influence varies among pollutants and according to the scale at which it is measured. Physiographic (drainage density) variation also influenced sediment and phosphorus loads. The effect of historical land use on present-day channel erosion and variation in soil texture are the most important sources of phosphorus and sediment that riparian buffers cannot attenuate. However, in most watersheds, a large proportion (approximately 70%) of these pollutants can be eliminated from streams with buffers. Cumulative frequency distributions of load reduction potential indicate that targeting pollution reduction in the highest 10% of Wisconsin watersheds would reduce total phosphorus and sediment loads in the entire state by approximately 20%. These results support our approach of geographically targeting nonpoint source pollution reduction at multiple scales, including the watershed scale.     

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.     

Development of plans to restore degraded areas in the Great Lakes

The International Joint Commission's Water Quality Board has identified 42 Areas of Concern in the Great Lakes ecosystem where Great Lakes Water Quality Agreement objectives or jurisdictional standards, criteria or guidelines, established to protect uses, have been exceeded and remedial actions are necessary to restore beneficial uses. As a result of the 1985 report of the Water Quality Board, the eight Great Lakes states and the Province of Ontario committed themselves to developing a remedial action plan (RAP) to restore all uses in each Area of Concern within their political boundaries. Each RAP must identify the specific measures necessary to control existing sources of pollution, abate existing contamination (e.g., contaminated sediments), and restore all beneficial uses. Points which must be explicitly addressed in each RAP include: geographic extent of problem, beneficial uses impaired, causes of problems, remedial measures and a schedule for implementation, responsible agencies, and surveillance and monitoring activities that will be used to track effectiveness of remedial actions. The jurisdictions are responsible for developing RAPs, and the International Joint Commission is responsible for evaluating the adequacy of each RAP and tracking progress in restoring beneficial uses.     

Comparison of Indexes for Prioritizing Placement of Water Quality Buffers in Agricultural Watersheds1

Dosskey, Michael G. and Zeyuan Qiu, 2011. Comparison of Indexes for Prioritizing Placement of Water Quality Buffers in Agricultural Watersheds. Journal of the American Water Resources Association (JAWRA) 47(4):662‐671. DOI: 10.1111/j.1752‐1688.2011.00532.x Abstract: Five physically based, spatially distributed, empirical indexes were compared for the degree to which they identified the same or different locations in watersheds where vegetative buffers would function better for reducing agricultural nonpoint source pollution. All five indexes were calculated on a 10 m × 10 m digital elevation grid on agricultural land in the 144‐km² Neshanic River watershed in New Jersey. The indexes included the topography‐based Wetness Index (WI) and Topographic Index (TI) and three soil survey‐based indexes (sediment trapping efficiency [STE], water trapping efficiency [WTE], and groundwater interaction [GI]). Results showed that each index associated higher pollution risk and mitigation potential to a different part of the landscape. The WI and TI identified swales and riparian areas where runoff converges, whereas STE and WTE identified upland sites. The STE and WTE lack the fine scale of slope resolution and the accounting for convergent runoff patterns that can be important for properly locating buffers in some watersheds. The GI index indicates the existence of a shallow water table but the correspondence with WI‐ and TI‐identified sites was only modest. For watersheds where pollutant loading is generated by both saturation‐excess (emphasized by TI and WI) and infiltration‐excess processes (emphasized by STE and WTE), the indexes could be complementary. However, techniques would be needed for properly apportioning priority among sites identified by each index.