Evaluating stream water quality through land use analysis in two grassland catchments: impact of wetlands on stream nitrogen concentration

We evaluated the impacts of natural wetlands and various land uses on stream nitrogen concentration in two grassland-dominated catchments in eastern Hokkaido, Japan. Analyzing land use types in drainage basins, measuring denitrification potential of its soil, and water sampling in all seasons of 2003 were performed. Results showed a highly significant positive correlation between the concentration of stream NO3-N and the proportion of upland area in drainage basins in both catchments. The regression slope, which we assumed to reflect the impact on water quality, was 24% lower for the Akkeshi catchment (0.012 +/- 0.001) than for the Shibetsu catchment (0.016 +/- 0.001). In the Akkeshi catchment, there was a significant negative correlation between the proportion of wetlands in the drainage basins and stream NO3-N concentration. Stream dissolved organic nitrogen (DON) and carbon (DOC) concentrations were significantly higher in the Akkeshi catchment. Upland and urban land uses were strongly linked to increases in in-stream N concentrations in both catchments, whereas wetlands and forests tended to mitigate water quality degradation. The denitrification potential of the soils was highest in wetlands, medium in riparian forests, and lowest in grasslands; and was significant in wetlands and riparian forests in the Akkeshi catchment. The solubility of soil organic carbon (SOC) and soil moisture tended to determine the denitrification potential. These results indicate that the water environment within the catchments, which influences denitrification potential and soil organic matter content, could have caused the difference in stream water quality between the two catchments.     

Dissolved organic nitrogen regulation in freshwaters

Dissolved organic nitrogen (DON) has been hypothesized to play a major role in N cycling in a variety of ecosystems. Our aim was to assess the seasonal and concentration relationships between dissolved organic carbon (DOC), DON, and NO3- within 102 streams and 16 lakes within catchments of differing complexity situated in Wales. Further, we aimed to assess whether patterns of land use, soil type, and vegetation gave consistent trends in DON and dissolved inorganic nitrogen (DIN) relationships over a diverse range of catchments. Our results reinforce that DON constitutes a significant component of the total dissolved N pool typically representing 40 to 50% of the total N in streams and lakes but sometimes representing greater than 85% of the total dissolved N. Generally, the levels of DON were inversely correlated with the concentration of DIN. In contrast to DIN concentrations, which showed distinct seasonality, DON showed no consistent seasonal trend. We hypothesize that this reflects differences in the bioavailability of these two N types. The amount of DON, DOC, and DIN was significantly related to soil type with higher DON export from Histosol-dominated catchments in comparison with Spodosol-dominated watersheds. Vegetation cover also had a significant effect on DON concentrations independent of soil type with a nearly twofold decrease in DON export from forested catchments in comparison with nonforested watersheds. Due to the diversity in catchment DON behavior, we speculate that this will limit the adoption of DON as a broad-scale indicator of catchment condition for use in monitoring and assessment programs.     

Differences in Dissolved Organic Carbon and Nitrogen Responses to Storm‐Event and Ground‐Water Conditions in a Forested,Glaciated Watershed in Western New York1

Abstract: Differences in the storm‐event responses of dissolved organic carbon (DOC) and nitrogen (DON) in streamflow and ground water were evaluated for a glaciated forested watershed in western New York. Eight to ten storm events with varying rainfall amounts, intensities, and antecedent moisture conditions were studied for three catchments (1.6, 3.4, and 696 ha) over a three‐year period (2003‐2005). Concentrations of DOC in streamflow exiting the catchments were significantly higher for storm events following a dry period, whereas no similar response was observed for DON. Highest DON concentrations in streamflow were typically associated with storm events following wet antecedent moisture conditions. In addition to antecedent moisture conditions, DOC concentrations were also positively correlated with precipitation amounts, while DON did not reveal a consistent pattern. Streamwater and ground‐water concentrations of DOC during storm events were also strongly correlated with riparian ground‐water depths but a similar relationship was not observed for DON. Ground‐water DON concentrations were also more variable than DOC. We hypothesized that the differences in DOC and DON responses stemmed from the differences in catchment sources of these solutes. This study suggests that while DOC and DON are intrinsically linked as components of dissolved organic matter, their dynamics and exports from watersheds may be regulated by a different set of mechanisms and factors. Identifying these differences is critical for developing more reliable and robust models for transport of dissolved organic matter.     

Controls on catchment-scale patterns of phosphorus in soil, streambed sediment, and stream water

Many models of phosphorus (P) transfer at the catchment scale rely on input from generic databases including, amongst others, soil and land use maps. Spatially detailed geochemical data sets have the potential to improve the accuracy of the input parameters of catchment-scale nutrient transfer models. Furthermore, they enable the assessment of the utility of available, generic spatial data sets for the modeling and prediction of soil nutrient status and nutrient transfer at the catchment scale. This study aims to quantify the unique and joint contribution of soil and sediment properties, land cover, and point-source emissions to the spatial variation of P concentrations in soil, streambed sediments, and stream water at the scale of a medium-sized catchment. Soil parent material and soil chemical properties were identified as major factors controlling the catchment-scale spatial variation in soil total P and Olsen P concentrations. Soil type and land cover as derived from the generic spatial database explain 33.7% of the variation in soil total P concentrations and 17.4% of the variation in Olsen P concentrations. Streambed P concentrations are principally related to the major element concentrations in streambed sediment and P delivery from the hillslopes due to sediment erosion. During base flow conditions, the total phosphorus (<0.45 microm) concentrations in stream water are mainly controlled by the concentrations of P and the major elements in the streambed sediment.     

Factors controlling sediment and phosphorus export from two Belgian agricultural catchments

Sediment and total phosphorus (TP) export vary through space and time. This study was conducted to determine the factors controlling sediment and TP export in two agricultural catchments situated in the Belgian Loess Belt. At the outlet of these catchments runoff discharge was continuously measured and suspended sediment samples were taken during rainfall events. Within the catchments vegetation type and cover, soil surface parameters, erosion features, sediment pathways, and rainfall characteristics were monitored. Total P content and sediment characteristics such as clay, organic carbon, and suspended sediment concentration were correlated. Total sediment and TP export differ significantly between the monitored catchments. Much of the difference is due to the occurrence of an extreme event in one catchment and the morphology and spatial organization of land use in the catchments. In one catchment, the direct connection between erosive areas and the catchment outlet by means of a road system contributed to a high sediment delivery ratio (SDR) at the outlet. In the other catchment, the presence of a wide valley in the center of the catchment caused sediment deposition. Vegetation also had an effect on sediment production and deposition. Thus, many factors control sediment and TP export from small agricultural catchments; some of these factors are related to the physical catchment characteristics such as morphology and landscape structure and are (semi)permanent, while others, such as vegetation cover and land use, are time dependent.     

Spatial variability of soil phosphorus in relation to the topographic index and critical source areas: sampling for assessing risk to water quality

A measure of soil P status in agricultural soils is generally required for assisting with prediction of potential P loss from agricultural catchments and assessing risk for water quality. The objectives of this paper are twofold: (i) investigating the soil P status, distribution, and variability, both spatially and with soil depth, of two different first-order catchments; and (ii) determining variation in soil P concentration in relation to catchment topography (quantified as the "topographic index") and critical source areas (CSAs). The soil P measurements showed large spatial variability, not only between fields and land uses, but also within individual fields and in part was thought to be strongly influenced by areas where cattle tended to congregate and areas where manure was most commonly spread. Topographic index alone was not related to the distribution of soil P, and does not seem to provide an adequate indicator for CSAs in the study catchments. However, CSAs may be used in conjunction with soil P data for help in determining a more "effective" catchment soil P status. The difficulties in defining CSAs a priori, particularly for modeling and prediction purposes, however, suggest that other more "integrated" measures of catchment soil P status, such as baseflow P concentrations or streambed sediment P concentrations, might be more useful. Since observed soil P distribution is variable and is also difficult to relate to nationally available soil P data, any assessment of soil P status for determining risk of P loss is uncertain and problematic, given other catchment physicochemical characteristics and the sampling strategy employed.     

Flow,Organic, and Inorganic Sediment Yields from a Channelized Watershed in the South Carolina Lower Coastal Plain

Many small streams in coastal watersheds in the southeastern United States are modified for agricultural, residential, and commercial development. In the South Carolina Lower Coastal Plain, low‐relief topography and a shallow water table make stream channelization ubiquitous. To quantify the impacts of urbanization and stream channelization, we measured flow and sediment from an urbanizing watershed and a small forested watershed. Flow and sediment export rates were used to infer specific yields from forested and nonforested regions of the urbanizing watershed. Study objectives were to: (1) quantify the range of runoff‐to‐rainfall ratios; (2) quantify the range of specific sediment yields; (3) characterize the quantity and quality of particulate matter exported; and (4) estimate sediment yield attributable to agriculture, development, and channelization activities in the urbanizing watershed. Our results showed that the urban watershed exported over five times more sediment per unit area compared with the forested watershed. Sediment concentration was related to flow flashiness in the urban watershed and to flow magnitude in the forested watershed. Sediments from the forested watershed were dominated by organic matter, whereas mineral matter dominated sediment from the urban stream. Our results indicated that a significant shift in sediment quality and quantity are likely to occur as forested watersheds are transformed by urbanization in coastal South Carolina.     

Historical Land-Use Influences the Long-Term Stream Turbidity Response to a Wildfire

Wildfires commonly result in an increase in stream turbidity. However, the influence of pre-fire land-use practices on post-fire stream turbidity is not well understood. The Lower Cotter Catchment (LCC) in south-eastern Australia is part of the main water supply catchment for Canberra with land in the catchment historically managed for a mix of conservation (native eucalypt forest) and pine (Pinus radiata) plantation. In January 2003, wildfires burned almost all of the native and pine forests in the LCC. A study was established in 2005 to determine stream post-fire turbidity recovery within the native and pine forest areas of the catchment. Turbidity data loggers were deployed in two creeks within burned native forest and burned pine forest areas to determine turbidity response to fire in these areas. As a part of the study, we also determined changes in bare soil in the native and pine forest areas since the fire. The results suggest that the time, it takes turbidity levels to decrease following wildfire, is dependent upon the preceding land-use. In the LCC, turbidity levels decreased more rapidly in areas previously with native vegetation compared to areas which were previously used for pine forestry. This is likely because of a higher percentage of bare soil areas for a longer period of time in the ex-pine forest estate and instream stores of fine sediment from catchment erosion during post-fire storm events. The results of our study show that the previous land-use may exert considerable control over on-going turbidity levels following a wildfire.     

Transport and metabolic fate of sewage particles in a recipient stream

Although the implementation of wastewater treatment plants (WWTP) has dramatically increased the quality of surface waters in urbanized areas, WWTPs can still discharge noticeable amounts of solutes and particles to recipient streams. Although the fate of WWTP nutrients has received considerable attention, transport and in-stream transformation of sewage-derived particulate organic matter (SDPOM) have not. To investigate the transport and transformation of SDPOM in recipient streams, we experimentally injected fluorescently labeled SDPOM into a headwater stream and tracked its downstream fate at baseflow. Most SDPOM disappeared from the streamwater within a 160-m long reach with an average deposition velocity of 0.14 mm s(-1). We further coupled hydrometric measurements of specific water fluxes through the streambed interface with a mixing model to estimate streambed oxygen removal, and found significantly higher oxygen removal in the deposition (0.75 g O2 m(-2) d(-1)) than in the downstream post-deposition (0.36 g O2 m(-2) d(-1)) subreach. Contrary to our expectations, we did not detect any apparent effect of SDPOM deposition on streambed clogging. Our results show the capacity of a recipient stream to retain SDPOM and to reduce its downstream export, and thus contribute to a better understanding of ecosystem services of human-altered streams.     

Stream Communities Along a Catchment Land-Use Gradient: Subsidy-Stress Responses to Pastoral Development

When native grassland catchments are converted to pasture, the main effects on stream physicochemistry are usually related to increased nutrient concentrations and fine-sediment input. We predicted that increasing nutrient concentrations would produce a subsidy-stress response (where several ecological metrics first increase and then decrease at higher concentrations) and that increasing sediment cover of the streambed would produce a linear decline in stream health. We predicted that the net effect of agricultural development, estimated as percentage pastoral land cover, would have a nonlinear subsidy-stress or threshold pattern. In our suite of 21 New Zealand streams, epilithic algal biomass and invertebrate density and biomass were higher in catchments with a higher proportion of pastoral land cover, responding mainly to increased nutrient concentration. Invertebrate species richness had a linear, negative relationship with fine-sediment cover but was unrelated to nutrients or pastoral land cover. In accord with our predictions, several invertebrate stream health metrics (Ephemeroptera–Plecoptera–Trichoptera density and richness, New Zealand Macroinvertebrate Community Index, and percent abundance of noninsect taxa) had nonlinear relationships with pastoral land cover and nutrients. Most invertebrate health metrics usually had linear negative relationships with fine-sediment cover. In this region, stream health, as indicated by macroinvertebrates, primarily followed a subsidy-stress pattern with increasing pastoral development; management of these streams should focus on limiting development beyond the point where negative effects are seen.     

Recent Changes in Stream Flashiness and Flooding,and Effects of Flood Management in North Carolina and Virginia

The southeastern United States has undergone anthropogenic changes in landscape structure, with the potential to increase (e.g., urbanization) and decrease (e.g., reservoir construction) stream flashiness and flooding. Assessment of the outcome of such change can provide insight into the efficacy of current strategies and policies to manage water resources. We (1) examined trends in precipitation, floods, and stream flashiness and (2) assessed the relative influence of land cover and flow‐regulating features (e.g., best management practices and artificial water bodies) on stream flashiness from 1991 to 2013. We found mean annual precipitation decreased, which coincided with decreasing trends in floods. In contrast, stream flashiness, overall, showed an increasing trend during the period of study. However, upon closer examination, 20 watersheds showed stable stream flashiness, whereas 5 increased and 6 decreased in flashiness. Urban watersheds were among those that increased or decreased in flashiness. Watersheds that increased in stream flashiness gained more urban cover, lost more forested cover and had fewer best management practices installed than urban watersheds that decreased in stream flashiness. We found best management practices are more effective than artificial water bodies in regulating flashy floods. Flashiness index is a valuable and straightforward metric to characterize changes in streamflow and help to assess the efficacy of management interventions.     

Dissolved organic carbon and disinfection by-product precursor release from managed peat soils

A wetland restoration demonstration project examined the effects of a permanently flooded wetland on subsidence of peat soils. The project, started in 1997, was done on Twitchell Island, in the Sacramento-San Joaquin Delta of California. Conversion of agricultural land to a wetland has changed many of the biogeochemical processes controlling dissolved organic carbon (DOC) release from the peat soils, relative to the previous land use. Dissolved organic C in delta waters is a concern because it reacts with chlorine, added as a disinfectant in municipal drinking waters, to form carcinogenic disinfection byproducts (DBPs), including trihalomethanes (THMs) and haloacetic acids (HAAs). This study explores the effects of peat soil biogeochemistry on DOC and DBP release under agricultural and wetland management. Results indicate that organic matter source, extent of soil organic matter decomposition, and decomposition pathways all are factors in THM formation. The results show that historical management practices dominate the release of DOC and THM precursors. However, within-site differences indicate that recent management decisions can contribute to changes in DOC quality and THM precursor formation. Not all aromatic forms of carbon are highly reactive and certain environmental conditions produce the specific carbon structures that form THMs. Both HAA and THM precursors are elevated in the DOC released under wetland conditions. The findings of this study emphasize the need to further investigate the roles of organic matter sources, microbial decomposition pathways, and decomposition status of soil organic matter in the release of DOC and DBP precursors from delta soils under varying land-use practices.     

Spatial and temporal modeling of microbial contaminants on grazing farmlands

This paper introduces an integrated spatial and temporal modeling system developed mathematically for assessing microbial contaminants on animal-grazed farmlands. The model uses fecal coliform, specifically Escherichia coli, as an indicator of fecal contamination and describes the sources, sinks, transport processes, and fate of E. coli contaminants in catchments and associated streams. Spatial features include grazing location, land topography, distance to a nearby stream, and distance through the stream network to the outlet. Temporal features are population dynamics on the land surface, in flow, and on streambeds. The model applies the principles of conservation of mass balance on two different types of pools: grid cells on land surfaces and networked stream segments. The model aims to improve the prediction of the effects of different land management strategies on the fecal contamination of waterways. This is achieved by characterizing the movement of fecal contaminants from land to streams and in-stream mobilization. Processes of attenuation, diffusion, and transport govern the movement. Our study site is a hill land catchment with an area of 140 ha and is used exclusively for animal grazing. The model was calibrated with previous research results, and then tested using the data collected at the outlet of the catchment. The sensitivity of the model predictions was analyzed for different scenarios: effect of stock rate, attenuation rate, and flow volumes. The similar pattern between monitored and predicted E. coli concentration proved that the model captures the key features that control the population dynamics of fecal contaminants. Further experiments are required to expand the model's functionality for covering more mitigation options.     

A CASE STUDY OF SHALLOW FLOW PATHS IN A STEEP ZERO-ORDER BASIN1

ABSTRACT: Soil water potentials, slope throughflow, runoff chemistry, and isotopic composition were monitored in a 97 m² zero-order basin within the Maimai 8 watershed on the South Island of New Zealand, for a natural rain storm and two artificial water applications. Contrary to results previously reported for other portions of the Maimai catchment, much of the runoff occurred as a shallow subsurface organic layer flow. For the 47 mm natural rain event, pre-storm soil matric potential ranged from ?60 to ?150 cm H₂O. No saturation was produced within the profile, and the majority of storm runoff emanated from flow within the organic horizon perched on the mineral soil surface. Hillslope applications corroborated this interpretation by showing >90 percent new water flushing with negligible mineral soil moisture response. Although the mechanisms cited in the text are not representative of the entire catchment, the study demonstrates: (1) the value of a combined physical-chemical-isotopic approach in quantifying slope processes, and (2) the heterogeneous nature and diversity of slope runoff pathways in a relatively homogeneous catchment.     

Nitrogen and phosphorus attenuation within the stream network of a coastal, agricultural watershed

Streams alter the concentration of nutrients they transport and thereby influence nutrient loading to estuaries downstream; however, the relationship between in-stream uptake, discharge variability, and subsequent nutrient export is poorly understood. In this study, in-stream N and P uptake were examined in the stream network draining a row-crop agricultural operation in coastal North Carolina. The effect of in-stream nutrient uptake on estuarine loading was examined using continuous measurements of watershed nutrient export. From August to December 2003, 52 and 83% of the NH4+ and PO4(3-) loads were exported during storms while concurrent storm flow volume was 34% of the total. Whole-ecosystem mass transfer velocities (Vf) of NH4+ and PO4(3-), measured using short-term additions of inorganic nutrients, ranged from 0.1 to 25 mm min(-1). Using a mass balance approach, this in-stream uptake was found to attenuate 65 to 98% of the NH4+ flux and 78 to 98% of the PO4(3-) flux in small, first-order drainage ditches. For the larger channel downstream, an empirical model based on Vf and discharge was developed to estimate the percentage of the nutrient load retained in-stream. The model predicted that all of the upstream NH4+ and PO4(3-) load was retained during base flow, while 65 and 37% of the NH4+ and PO4(3-) load was retained during storms. Remineralization from the streambed (vs. terrestrial sources) was the apparent source of NH4+ and PO4(3-) to the estuary during base flow. In-stream uptake reduced the dissolved inorganic N to dissolved inorganic P ratio of water exported to the N-limited estuary, thus limiting the potential for estuarine phytoplankton growth.     

Controls on Nutrients Across a Prairie Stream Watershed: Land Use and Riparian Cover Effects

Nutrient inputs generally are increased by human-induced land use changes and can lead to eutrophication and impairment of surface waters. Understanding the scale at which land use influences nutrient loading is necessary for the development of management practices and policies that improve water quality. The authors assessed the relationships between land use and stream nutrients in a prairie watershed dominated by intermittent stream flow in the first-order higher elevation reaches. Total nitrogen, nitrate, and phosphorus concentrations were greater in tributaries occupying the lower portions of the watershed, closely mirroring the increased density of row crop agriculture from headwaters to lower-elevation alluvial areas. Land cover classified at three spatial scales in each sub-basin above sampling sites (riparian in the entire catchment, catchment land cover, and riparian across the 2 km upstream) was highly correlated with variation in both total nitrogen (r² = 53%, 52%, and 49%, respectively) and nitrate (r² = 69%, 65%, and 56%, respectively) concentrations among sites. However, phosphorus concentrations were not significantly associated with riparian or catchment land cover classes at any spatial scale. Separating land use from riparian cover in the entire watershed was difficult, but riparian cover was most closely correlated with in-stream nutrient concentrations. By controlling for land cover, a significant correlation of riparian cover for the 2 km above the sampling site with in-stream nutrient concentrations could be established. Surprisingly, land use in the entire watershed, including small intermittent streams, had a large influence on average downstream water quality although the headwater streams were not flowing for a substantial portion of the year. This suggests that nutrient criteria may not be met only by managing permanently flowing streams.     

USE OF STAGE DATA TO CHARACTERIZE HYDROLOGIC CONDITIONS IN AN URBANIZING ENVIRONMENT1

ABSTRACT: This paper presents the results of a study on the use of continuous stage data to describe the relation between urban development and three aspects of hydrologic condition that are thought to influence stream ecosystems—overall stage variability, stream flashiness, and the duration of extreme‐stage conditions. This relation is examined using data from more than 70 watersheds in three contrasting environmental settings—the humid Northeast (the metropolitan Boston, Massachusetts, area); the very humid Southeast (the metropolitan Birmingham, Alabama, area); and the semiarid West (the metropolitan Salt Lake City, Utah, area). Results from the Birmingham and Boston studies provide evidence linking increased urbanization with stream flashiness. Fragmentation of developed land cover patches appears to ameliorate the effects of urbanization on overall variability and flashiness. There was less success in relating urbanization and streamflow conditions in the Salt Lake City study. A related investigation of six North Carolina sites with long term discharge and stage data indicated that hydrologic condition metrics developed using continuous stage data are comparable to flow based metrics, particularly for stream flashiness measures.     

FOREST HYDROLOGY IN THE PACIFIC NORTHWEST: ADDITIONAL RESEARCH NEEDS1

ABSTRACT: While much is known about the hydrology of forested mountain catchments in the Pacific Northwest, important research questions remain. For example, the dynamics of storm precipitation amounts and the modeling of catchment outflows represent a continuing research need. Without an improved understanding of the spatial and temporal aspects of storm precipitation patterns, our ability to evaluate and improve physically-based hydrologic models is limited. From a practical perspective, tens of thousands of kilometers of access roads have been constructed across forested catchments of the Pacific Northwest. Yet, few forestry research programs focus on road drainage (e.g., ditches, culverts, fords). The few studies that address this issue indicate road drainage systems need to function effectively over a wide range of flow events and terrain conditions. In addition, historical forest practices associated with hillslopes and riparian systems have altered the character of many Pacific Northwest aquatic ecosystems. If restoration of these systems is to be effective, research efforts are needed to better understand the linkages between riparian forests, geomorphic processes, and hydrologic disturbance regimes.     

Impact of land drainage on peatland hydrology

There is a long history of drainage of blanket peat but few studies of the long-term hydrological impact of drainage. This paper aims to test differences in runoff production processes between intact and drained blanket peat catchments and determine whether there have been any long-term changes in stream flow since drainage occurred. Hillslope runoff processes and stream discharge were measured in four blanket peat catchments. Two catchments were drained with open-cut ditches in the 1950s. Ditching originally resulted in shorter lag times and flashier storm hydrographs but no change in the annual catchment runoff efficiency. In the period between 2002 and 2004, the hydrographs in the drained catchments, while still flashy, were less sensitive to rainfall than in the 1950s and the runoff efficiency had significantly increased. Drains resulted in a distinctive spatial pattern of runoff production across the slopes. Overland flow was significantly lower in the drained catchments where throughflow was more dominant. In the intact peatlands, matrix throughflow produced by peat layers below 10 cm was rare and produced <1% of the runoff. However, in drained peatlands, matrix throughflow in deeper peat layers was common and provided around 23% of the runoff from gauged plots. Macropore flow, the density of soil piping, and pipeflow were significantly greater in drained peatlands than in intact basins. Gradual changes to peat structure could explain the long-term changes in river flow, which are in addition to those occurring in the immediate aftermath of peatland drainage.     

Landscape metrics and estuarine sediment contamination in the mid-Atlantic and southern New England regions

In a previously published study, quantitative relationships were developed between landscape metrics and sediment contamination for 25 small estuarine systems within Chesapeake Bay. These analyses have been extended to include 75 small estuarine systems across the mid-Atlantic and southern New England regions of the USA. Because of the different characteristics and dynamics of the estuaries across these regions, adjustment for differing hydrology, sediment characteristics, and sediment origins were included in the analysis. Multiple linear regression with stepwise selection was used to develop statistical models for sediment metals, organics, and total polycyclic aromatic hydrocarbons (PAHs). The landscape metrics important for explaining the variation in sediment metals levels (R2 = 0.72) were the percent area of nonforested wetlands (negative contribution), percent area of urban land, and point source effluent volume and metals input (positive contributions). The metrics important for sediment organics levels (R2 = 0.5) and total PAHs (R2 = 0.46) were percent area of urban land (positive contribution) and percent area of nonforested wetlands (negative contribution). These models included silt-clay content (metals) or total organic C (organics, total PAHs) of sediments and grouping by estuarine hydrology, suggesting the importance of sediment characteristics and hydrology in mitigating the influence of the landscape metrics on sediment contamination levels. The overall results from this study are indicative of how statistical models can be developed relating landscape metrics to estuarine sediment contamination for distributions of land cover and point source discharges.