THE CONCEPT OF HYDROLOGIC LANDSCAPES1

ABSTRACT: Hydrologic landscapes are multiples or variations of fundamental hydrologic landscape units. A fundamental hydrologic landscape unit is defined on the basis of land‐surface form, geology, and climate. The basic land‐surface form of a fundamental hydrologic landscape unit is an upland separated from a lowland by an intervening steeper slope. Fundamental hydrologic landscape units have a complete hydrologic system consisting of surface runoff, ground‐water flow, and interaction with atmospheric water. By describing actual landscapes in terms of land‐surface slope, hydraulic properties of soils and geologic framework, and the difference between precipitation and evapotranspiration, the hydrologic system of actual landscapes can be conceptualized in a uniform way. This conceptual framework can then be the foundation for design of studies and data networks, syntheses of information on local to national scales, and comparison of process research across small study units in a variety of settings. The Crow Wing River watershed in central Minnesota is used as an example of evaluating stream discharge in the context of hydrologic landscapes. Lake‐research watersheds in Wisconsin, Minnesota, North Dakota, and Nebraska are used as an example of using the hydrologic‐land‐scapes concept to evaluate the effect of ground water on the degree of mineralization and major‐ion chemistry of lakes that lie within ground‐water flow systems.     

Delineation and evaluation of hydrologic-landscape regions in the United States using geographic information system tools and multivariate statistical analyses

Hydrologic-landscape regions in the United States were delineated by using geographic information system (GIS) tools combined with principal components and cluster analyses. The GIS and statistical analyses were applied to land-surface form, geologic texture (permeability of the soil and bedrock), and climate variables that describe the physical and climatic setting of 43,931 small (approximately 200 km2) watersheds in the United States. (The term "watersheds" is defined in this paper as the drainage areas of tributary streams, headwater streams, and stream segments lying between two confluences.) The analyses grouped the watersheds into 20 noncontiguous regions based on similarities in land-surface form, geologic texture, and climate characteristics. The percentage of explained variance (R-squared value) in an analysis of variance was used to compare the hydrologic-landscape regions to 19 square geometric regions and the 21 U.S. Environmental Protection Agency level-II ecoregions. Hydrologic-landscape regions generally were better than ecoregions at delineating regions of distinct land-surface form and geologic texture. Hydrologic-landscape regions and ecoregions were equally effective at defining regions in terms of climate, land cover, and water-quality characteristics. For about half of the landscape, climate, and water-quality characteristics, the R-squared values of square geometric regions were as high as hydrologic-landscape regions or ecoregions.     

Multivariate Condition Assessment of Watersheds with Linked Micromaps

A challenge for statewide stream monitoring is visualizing the spatial and statistical characteristics of such data to compare the biotic condition of watersheds and relate that condition to watershed‐level estimates of instream variables. We used linked micromaps on stream survey data of 25 subbasins (766‐5,982 km²) for biotic condition, nine water quality, and two habitat variables. Subbasin biotic condition was negatively correlated with conductivity, magnesium and sulfate concentrations, and weakly positively correlated with habitat scores of sedimentation and embeddedness, with higher scores indicating better habitat. Positive spatial autocorrelation was detected among the subbasins in both habitat variables, iron concentration, pH, and exceedances of fecal coliform criteria as shown in linked micromaps. A spatial principal components analysis reduced the 11 environmental variables to two principal axes. The first axis synthesized a gradient of water quality and habitat scores among the subbasins. Subbasin biotic condition regressed on first axis subbasin scores had a significant, negative slope and accounted for 55% of the variation. Subbasins in degraded biotic condition had elevated conductivities and ion concentrations in northern and southern subbasins, and low habitat scores in western subbasins. Through linked micromaps, we compared the biotic condition among subbasins and identified stressors prevalent among subbasins that affected biotic condition.     

MODELING FRAMEWORK FOR MANAGING COPPER RUNOFF IN URBAN WATERSHEDS1

ABSTRACT: A modeling framework was developed for managing copper runoff in urban watersheds that incorporates water quality characterization, watershed land use areas, hydrologic data, a statistical simulator, a biotic ligand binding model to characterize acute toxicity, and a statistical method for setting a watershed specific copper loading. The modeling framework is driven by export coefficients derived from water quality parameters and hydrologic inputs measured in an urban watershed's storm water system. This framework was applied to a watershed containing a copper roof built in 1992. A series of simulations was run to predict the change in receiving stream water chemistry caused by roof aging and to determine the maximum copper loading (at the 99 percent confidence level) a watershed could accept without causing acute toxicity in the receiving stream. Forecasting the amount of copper flux responsible for exceeding the assimilation capacity of a watershed can be directly related to maximum copper loadings responsible for causing toxicity in the receiving streams. The framework developed in this study can be used to evaluate copper utilization in urban watersheds.     

RELATION BETWEEN STREAM-WATER QUALITY AND GEOHYDROLOGY DURING BASE-FLOW CONDITIONS,ROBERTS CREEK WATERSHED,CLAYTON COUNTY,IOWA1

ABSTRACT: An investigation to determine the relation between stream water quality and geohydrology in the Roberts Creek watershed, Clayton County, Iowa, was conducted during selected base-flow periods in 1988–90. Discharge measurements were made and water samples collected for analyses of nutrients and selected herbicides in 19 subbasins along the main stem and tributaries of Roberts Creek. The areal extent of unconsolidated and bedrock units subcropping in each subbasin was quantified. The hydrologic data were correlated statistically with the geologic data to determine relations. Roberts Creek generally gained water and had larger nitrogen concentrations in subbasins in which bess and alluvial material were underlain primarily by low-permeability till and shale units. Roberts Creek generally lost water and had lower nitrate concentrations in subbasins with subcroppmg karstic units. Nitrogen concentrations decreased in streams underlain by the karstic units because the nitrogen removed by biological processes was not replaced by ground-water inflow. Seepage from Roberts Creek to ground water in areas of subcropping karstic carbonate rocks reduced the flow, which reduced the velocity, causing increased residence time of water in the stream. The additional residence time may allow additional time for biological processes to remove nitrogen from solution. There was no significant relation between dissolved orthophosphate or atrazine and the underlying geology.     

EFFECTS OF MULTI‐SCALE ENVIRONMENTAL CHARACTERISTICS ON AGRICULTURAL STREAM BIOTA IN EASTERN WISCONSIN1

ABSTRACT: The U.S. Geological Survey examined 25 agricultural streams in eastern Wisconsin the determine relations between fish, invertebrate, and algal metrics and multiple spatial scales of land cover, geologic setting, hydrologic, aquatic habitat, and water chemistry data. Spearman correlation and redundancy analyses were used to examine relations among biotic metrics and environmental characteristics. Riparian vegetation, geologic, and hydrologic conditions affected the response of biotic metrics to watershed agricultural land cover but the relations were aquatic assemblage dependent. It was difficult to separate the interrelated effects of geologic setting, watershed and buffer land cover, and base flow. Watershed and buffer land cover, geologic setting, reach riparian vegetation width, and stream size affected the fish IBI, invertebrate diversity, diatom IBI, and number of algal taxa; however, the invertebrate FBI, percentage of EPT, and the diatom pollution index were more influenced by nutrient concentrations and flow variability. Fish IBI scores seemed most sensitive to land cover in the entire stream network buffer, more so than watershed‐scale land cover and segment or reach riparian vegetation width. All but one stream with more than approximately 10 percent buffer agriculture had fish IBI scores of fair or poor. In general, the invertebrate and algal metrics used in this study were not as sensitive to land cover effects as fish metrics. Some of the reach‐scale characteristics, such as width/depth ratios, velocity, and bank stability, could be related to watershed influences of both land cover and geologic setting. The Wisconsin habitat index was related to watershed geologic setting, watershed and buffer land cover, riparian vegetation width, and base flow, and appeared to be a good indicator of stream quality Results from this study emphasize the value of using more than one or two biotic metrics to assess water quality and the importance of environmental characteristics at multiple scales.     

Linking Nitrogen Export to Landscape Heterogeneity: The Role of Infrastructure and Storm Flows in a Mediterranean Urban System

Urban ecosystems are often sources of nonpoint source (NPS) nitrogen (N) pollution to aquatic ecosystems. However, N export from urban watersheds is highly variable. Examples of densely urbanized watersheds are not well studied, and these may have comparatively low export rates. Commonly used metrics of landscape heterogeneity may obscure our ability to discern relationships among landscape characteristics that can explain these lower export rates. We expected that differences not often captured by these metrics in the relative cover of vegetation, structures, and impervious surfaces would better explain observed variation in N export. We examined these relationships during storms in residential watersheds. Contrary to expectations, land cover did not directly predict variation in N or water export. Instead, N export was strongly linked to drainage infrastructure density. Our research highlights the role of fine‐scaled landscape attributes, mainly infrastructure, in explaining patterns of N export from densely urbanized watersheds. Changes to hydrologic flow paths by infrastructure explained more variation in N export than land cover. Our findings support further development of landscape ecological models of urban N export that focus on hydrologic modification by infrastructure rather than traditional landscape measures such as land use, as indicators for evaluating patterns of NPS nitrogen pollution in densely urbanized watersheds.     

Oregon Hydrologic Landscapes: A Classification Framework1

Wigington, Parker J., Jr., Scott G. Leibowitz, Randy L. Comeleo, and Joseph L. Ebersole, 2012. Oregon Hydrologic Landscapes: A Classification Framework. Journal of the American Water Resources Association (JAWRA) 1‐20. DOI: 10.1111/jawr.12009 Abstract: There is a growing need for hydrologic classification systems that can provide a basis for broad‐scale assessments of the hydrologic functions of landscapes and watersheds and their responses to stressors such as climate change. We developed a hydrologic landscape (HL) classification approach that describes factors of climate‐watershed systems that control the hydrologic characteristics of watersheds. Our assessment units are incremental watersheds (i.e., headwater watersheds or areas draining directly into stream reaches). Major components of the classification include indices of annual climate, climate seasonality, aquifer permeability, terrain, and soil permeability. To evaluate the usefulness of our approach, we identified 30 rivers with long‐term streamflow‐gauging records and without major diversions and impoundments. We used statistical clustering to group the streams based on the shapes of their annual hydrographs. Comparison of the streamflow clusters and HL distributions within river basin clusters shows that the Oregon HL approach has the ability to provide insights about the expected hydrologic behavior of HLs and larger river basins. The Oregon HL approach has potential to be a useful framework for comparing hydrologic attributes of streams and rivers in the Pacific Northwest.     

GIS-ASSISTED REGRESSION ANALYSIS TO IDENTIFY SOURCES OF SELENIUM IN STREAMS1

ABSTRACT: Using a geographic information system, a regression model has been developed to identify and to assess potential sources of selenium in the Kendrick Reclamation Project Area, Wyoming. A variety of spatially distributed factors was examined to determine which factors are most likely to affect selenium discharge in tributaries to the North Platte River. Areas of Upper Cretaceous Cody Shale and Quaternary alluvial deposits and irrigated land, length of irrigation canals, and boundaries of hydrologic subbasins of the major tributaries to the North Platte River were digitized and stored in a geographic information system. Selenium concentrations in samples of soil, plant material, ground water, and surface water were determined and evaluated. The location of all sampling sites was digitized and stored in the geographic information system, together with the selenium concentrations in samples. A regression model was developed using stepwise multiple regression of median selenium discharges on the physical and chemical characteristics of hydrologic subbasins. Results indicate that the intensity of irrigation in a hydrologic subbasin, as determined by area of irrigated land and length of irrigation delivery canals, accounts for the largest variation in median selenium discharges among subbasins. Tributaries draining hydrologic subbasins with greater intensity of irrigation result in greater selenium discharges to the North Platte River than do tributaries draining subbasins with lesser intensity of irrigation.     

HYDROLOGIC LANDSCAPES ON THE DELMARVA PENINSULA PART 2: ESTIMATES OF BASE-FLOW NITROGEN LOAD TO CHESAPEAKE BAY1

ABSTRACT: Base-flow samples were collected from 47 sampling sites for four seasons from 1990–91 on the Delmarva Peninsula in Delaware and Maryland to relate stream chemistry to a “hydrologic landscape” and season. Two hydrologic landscapes were determined: (1) a well-drained landscape, characterized by a combination of a low percentage of forest cover, a low percentage of poorly drained soil, and elevated channel slope; and (2) poorly drained landscape, characterized by a combination of an elevated percentage of forest cover, an elevated percentage of poorly drained soil, and low channel slope. Concentrations of nitrogen were significantly related to the hydrologic landscape. Nitrogen concentrations tended to be higher in well-drained landscapes than in poorly drained ones. The highest instantaneous nitrogen yields occurred in well-drained landscapes during the winter. These yields were extrapolated over the part of the study area draining to Chesapeake Bay in order to provide a rough estimate of nitrogen load from base flow to the Bay and its estuarine tributaries. This estimate was compared to an estimate made by extrapolating from an existing long-term monitoring station. The load estimate from the stream survey data was 5 ± 10⁶ kg of N per year, which was about four times the estimate, made from the existing long-term monitoring station. The stream-survey estimate of base flow represents about 40 percent of the total nitrogen load that enters the Bay and estuarine tributaries from all sources in the study area.     

In search of effective scales for stream management: does agroecoregion,watershed, or their intersection best explain the variance in stream macroinvertebrate communities?

Our lack of understanding of relationships between stream biotic communities and surrounding landscape conditions makes it difficult to determine the spatial scale at which management practices are best assessed. We investigated these relationships in the Minnesota River Basin, which is divided into major watersheds and agroecoregions which are based on soil type, geologic parent material, landscape slope steepness, and climatic factors affecting crop productivity. We collected macroinvertebrate and stream habitat data from 68 tributaries among three major watersheds and two agroecoregions. We tested the effectiveness of the two landscape classification systems (i.e., watershed, agroecoregion) in explaining variance in habitat and macroinvertebrate metrics, and analyzed the relative influence on macroinvertebrates of local habitat versus regional characteristics. Macroinvertebrate community composition was most strongly influenced by local habitat; the variance in habitat conditions was best explained at the scale of intersection of major watershed and agroecoregion (i.e., stream habitat conditions were most homogeneous within the physical regions of intersection of these two landscape classification systems). Our results are consistent with findings of other authors that most variation in macroinvertebrate community data from large agricultural catchments is attributable to local physical conditions. Our results are the first to test the hypothesis and demonstrate that the scale of intersection best explains these variances. The results suggest that management practices adjusted for both watershed and ecoregion characteristics, with the goal of improving physical habitat characteristics of local streams, may lead to better basin-wide water quality conditions and stream biological integrity.     

Spatial Relationships of Levees and Wetland Systems within Floodplains of the Wabash Basin,USA

Given the unique biogeochemical, physical, and hydrologic services provided by floodplain wetlands, proper management of river systems should include an understanding of how floodplain modifications influence wetland ecosystems. The construction of levees can reduce river–floodplain connectivity, yet it is unclear how levees affect wetlands within floodplains, let alone the cumulative impacts within an entire watershed. This paper explores spatial relationships between levee and floodplain wetland systems in the Wabash Basin, United States. We used a hydrogeomorphic floodplain delineation technique to map floodplain extents and identify wetlands that may be hydrologically connected to river networks. We then spatially examined the relationship between levee presence, wetland area, and other river network attributes within discrete subbasins. Our results show that cumulative wetland area is relatively constant in subbasins that contain levees, regardless of maximum stream order within the subbasin. In subbasins that do not contain levees, cumulative wetland area increases with maximum stream order. However, we found that wetland distributions around levees can be complex, and further studies on the influence of levees on wetland habitat may need to consider finer resolution spatial scales.     

Integrating Landscape Ecology and Geoinformatics to Decipher Landscape Dynamics for Regional Planning

We used remote sensing and GIS in conjunction with multivariate statistical methods to: (i) quantify landscape composition (land cover types) and configuration (patch density, diversity, fractal dimension, contagion) for five coastal watersheds of Kalloni gulf, Lesvos Island, Greece, in 1945, 1960, 1971, 1990 and 2002/2003, (ii) evaluate the relative importance of physical (slope, geologic substrate, stream order) and human (road network, population density) variables on landscape composition and configuration, and (iii) characterize processes that led to land cover changes through land cover transitions between these five successive periods in time. Distributions of land cover types did not differ among the five time periods at the five watersheds studied because the largest cumulative changes between 1945 and 2002/2003 did not take place at dominant land cover types. Landscape composition related primarily to the physical attributes of the landscape. Nevertheless, increase in population density and the road network were found to increase heterogeneity of the landscape mosaic (patchiness), complexity of patch shape (fractal dimension), and patch disaggregation (contagion). Increase in road network was also found to increase landscape diversity due to the creation of new patches. The main processes involved in land cover changes were plough-land abandonment and ecological succession. Landscape dynamics during the last 50 years corroborate the ecotouristic-agrotouristic model for regional development to reverse trends in agricultural land abandonment and human population decline and when combined with hypothetical regulatory approaches could predict how this landscape could develop in the future, thus, providing a valuable tool to regional planning.     

NEAR‐STREAM LANDUSE EFFECTS ON STREAMWATER NUTRIENT DISTRIBUTION IN AN URBANIZING WATERSHED1

ABSTRACT: We investigated spatial and temporal relationships among surface and subsurface watershed attributes and stream nutrient concentrations in urbanizing Johnson Creek watershed in northern Oregon. We sampled stream water at eight urban and five nonurban locations from March 1998 through December 1999. We sampled eight wells distributed over the two primary aquifers in the watershed. Using a Geographic Information System (GIS), percentages of landuse attributes within a radius of 30, 91, and 152 m from each sample site were quantified. We analyzed relationships between (1) nutrient concentrations and percentage cover of different landuse attributes, and (2) nutrient concentrations and underlying hydrologic units. We did not find a significant relationship between ground water chemistry and stream water chemistry. We found elevated levels of phosphorus (P) concentrations correlated with urban landuse, while higher nitrogen (N) concentrations were correlated with nonurban (primarily agricultural) landuse. We concluded that elevated levels of N in nonurban areas of Johnson Creek watershed were associated with agricultural practices. We further concluded that urban development factors such as increases in storm drains, dry wells, and impermeable surfaces may be responsible for higher input of P to the stream in urbanizing areas of the Johnson Creek watershed.     

Surface‐Water Nutrient Conditions and Sources in the United States Pacific Northwest1

Wise, Daniel R. and Henry M. Johnson, 2011. Surface‐Water Nutrient Conditions and Sources in the United States Pacific Northwest. Journal of the American Water Resources Association (JAWRA) 47(5):1110‐1135. DOI: 10.1111/j.1752‐1688.2011.00580.x Abstract: The SPAtially Referenced Regressions On Watershed attributes (SPARROW) model was used to perform an assessment of surface‐water nutrient conditions and to identify important nutrient sources in watersheds of the Pacific Northwest region of the United States (U.S.) for the year 2002. Our models included variables representing nutrient sources as well as landscape characteristics that affect nutrient delivery to streams. Annual nutrient yields were higher in watersheds on the wetter, west side of the Cascade Range compared to watersheds on the drier, east side. High nutrient enrichment (relative to the U.S. Environmental Protection Agency’s recommended nutrient criteria) was estimated in watersheds throughout the region. Forest land was generally the largest source of total nitrogen stream load and geologic material was generally the largest source of total phosphorus stream load generated within the 12,039 modeled watersheds. These results reflected the prevalence of these two natural sources and the low input from other nutrient sources across the region. However, the combined input from agriculture, point sources, and developed land, rather than natural nutrient sources, was responsible for most of the nutrient load discharged from many of the largest watersheds. Our results provided an understanding of the regional patterns in surface‐water nutrient conditions and should be useful to environmental managers in future water‐quality planning efforts.     

REGIONAL CHARACTERISTICS OF NUTRIENT CONCENTRATIONS IN STREAMS AND THEIR APPLICATION TO NUTRIENT CRITERIA DEVELOPMENT1

ABSTRACT: In order to establish meaningful nutrient criteria, consideration must be given to the spatial variations in geographic phenomena that cause or reflect differences in nutrient concentrations in streams. Regional differences in stream nutrient concentrations were illustrated using stream data collected from 928 nonpoint‐source watersheds distributed throughout the country and sampled as part of the U.S. EPA National Eutrophication Survey (NES). Spatial patterns in the differences were compared and found to correspond with an a priori regional classification system based on regional patterns in landscape attributes associated with variation in nutrient concentrations. The classification consists of 14 regions composed of aggregations of the 84 U.S. EPA Level III Ecoregions. The primary distinguishing characteristics of each region and the factors associated with variability in water quality characteristics are presented. The use of the NES and many other extant monitoring data sets to develop regional reference conditions for nutrient concentrations in streams is discouraged on the basis of sample representation. The necessity that all sites used in such an effort be regionally representative and consistently screened for least possible impact is emphasized. These sampling issues are rigorously addressed by the U.S. EPA Environmental Monitoring and Assessment Program (EMAP). A case‐study, using EMAP data collected from the Central and Eastern Forested Uplands, demonstrates how regional reference conditions and draft nutrient criteria could be developed.     

A graphical screening method for assessing stream water quality using specific conductivity and alkalinity data

In areas of varying geology, it is difficult to infer water quality from specific conductance or electrical conductivity (EC) data without an understanding of the expected range of EC values based on local bedrock composition. This paper describes a user-friendly graphical screening method that addresses this issue by plotting the EC against concurrent alkalinity data, which correlates well with the presence of carbonate bedrock under natural conditions, and thus serves as an index of bedrock type. The upper limit of EC vs. alkalinity expected in a stream is determined using regional groundwater quality data, based on the assumption that stream chemistry reflects groundwater under baseflow conditions. Stream samples with EC/alkalinity values that consistently plot above this limit are considered impacted by anthropogenic sources. The effect of dilution and runoff on the EC vs. alkalinity plot of stream samples is considered using a simple baseflow/storm runoff-mixing model. The graphical method's utility as a screening tool is demonstrated by application to stream chemistry data from watersheds of southeastern Pennsylvania and northwestern New Jersey in several distinct geologic settings; however the method is general and widely applicable to watersheds in humid temperate regions. Its use is intended for watershed stewards of both professional and nonprofessional qualification.     

Headwater Influences on Downstream Water Quality

We investigated the influence of riparian and whole watershed land use as a function of stream size on surface water chemistry and assessed regional variation in these relationships. Sixty-eight watersheds in four level III U.S. EPA ecoregions in eastern Kansas were selected as study sites. Riparian land cover and watershed land use were quantified for the entire watershed, and by Strahler order. Multiple regression analyses using riparian land cover classifications as independent variables explained among-site variation in water chemistry parameters, particularly total nitrogen (41%), nitrate (61%), and total phosphorus (63%) concentrations. Whole watershed land use explained slightly less variance, but riparian and whole watershed land use were so tightly correlated that it was difficult to separate their effects. Water chemistry parameters sampled in downstream reaches were most closely correlated with riparian land cover adjacent to the smallest (first-order) streams of watersheds or land use in the entire watershed, with riparian zones immediately upstream of sampling sites offering less explanatory power as stream size increased. Interestingly, headwater effects were evident even at times when these small streams were unlikely to be flowing. Relationships were similar among ecoregions, indicating that land use characteristics were most responsible for water quality variation among watersheds. These findings suggest that nonpoint pollution control strategies should consider the influence of small upland streams and protection of downstream riparian zones alone is not sufficient to protect water quality.     

CHANGING NEAR‐STREAM LAND USE AND RWER CHANNEL MORPHOLOGY IN THE VENEZUELAN ANDES1

ABSTRACT: The shape of a river channel is linked to surrounding land use through interacting hydrologic and geologic processes. This study analyzes the relationship between the change in near‐stream land use and the shape of the adjacent river channel over time. Three watersheds in the foothills of the Venezuelan Andes that have experienced differing degrees of development were studied to determine river channel width, sinuosity, and position relative to surrounding land use. Change in land use over time was obtained from multiple‐date aerial photographs (1946 and 1980) referenced to 1996 Landsat Thematic Mapper (TM) satellite imagery, and verified by field inspection. Measurements of land‐use type and amount and river channel morphology from the two dates were made using geographic information system (GIS) methods. The three watersheds differed in the extent of deforestation, the location of remaining forested land, and how much land‐use change had already occurred by 1946. Change in river channel morphology was greatest at the most deforested sites. Valley shape and channel constraint also had a discernible effect on change in channel morphology. This study introduces a method for analyzing change in coupled terrestrial‐aquatic systems based on multiple‐date, remotely sensed data and GIS analysis of spatial properties. The results document human impacts on river channels through a comparison of multiple watersheds over a 35‐year time interval.