Spatial distribution of DDT in sediments from estuarine rivers of central Florida

Sediments may act as both a carrier for and a potential source of contaminants such as toxic organics in aquatic environments. This study investigated the spatial distribution of the pesticide DDT [1,1, 1-trichloro-2,2-bis(p-chlorophenyl)ethane] in sediments from the Cedar and Ortega Rivers located in the lower St. Johns River basin, Florida, USA, using field measurements and three-dimensional kriging analysis. High DDT concentrations were found near the junction of the Cedar and Ortega Rivers and at the north end of the Ortega River in the upper 0.5 m of the sediments, indicating that the sediment was enriched with DDT in the top layer although use of this chlorinated compound was banned in 1972. Further study revealed that the influence of sediment grain size or texture on DDT contamination was negligible in this river system and no linear correlations existed among DDT and its metabolites such as DDD [1,1-dichloro-2,2-bis(p-chlorophenyl)ethane] and DDE [1,1-dichloro-2,2-bis(p-chlorophenyl)ethylene]. Comparison of three-dimensional distribution of DDT content to the Florida sediment quality assessment guideline or probable effect level (PEL) showed several "hot spots" in the Ortega River sediments, where DDT contents exceeded the PEL value of 4.78 microg kg(-1). Such contamination may pose a significant hazard to aquatic life.     

Temporal and spatial distributions of sediment total organic carbon in an estuary river

Understanding temporal and spatial distributions of naturally occurring total organic carbon (TOC) in sediments is critical because TOC is an important feature of surface water quality. This study investigated temporal and spatial distributions of sediment TOC and its relationships to sediment contaminants in the Cedar and Ortega Rivers, Florida, USA, using three-dimensional kriging analysis and field measurement. Analysis of field data showed that large temporal changes in sediment TOC concentrations occurred in the rivers, which reflected changes in the characteristics and magnitude of inputs into the rivers during approximately the last 100 yr. The average concentration of TOC in sediments from the Cedar and Ortega Rivers was 12.7% with a maximum of 22.6% and a minimum of 2.3%. In general, more TOC accumulated at the upper 1.0 m of the sediment in the southern part of the Ortega River although the TOC sedimentation varied with locations and depths. In contrast, high concentrations of sediment contaminants, that is, total polycyclic aromatic hydrocarbons (PAHs) and polychlorinated biphenyls (PCBs), were found in sediments from the Cedar River. There was no correlation between TOC and PAHs or PCBs in these river sediments. This finding is in contradiction to some other studies which reported that the sorption of hydrocarbons is highly related to the organic matter content of sediments. This discrepancy occurred because of the differences in TOC and hydrocarbon source input locations. It was found that more TOC loaded into the southern part of the Ortega River, while almost all of the hydrocarbons entered into the Cedar River. This study suggested that the locations of their input sources as well as the land use patterns should also be considered when relating hydrocarbons to sediment TOC.     

Characterization of the pesticide chlordane in estuarine river sediments

Sediments are increasingly recognized as both carrier and potential source of contaminants in aquatic environments. This study investigated the characteristics and spatial distribution of total chlordane and its three most abundant compounds, including alpha-chlordane, gamma-chlordane, and trans-nonachlor, in sediments from the Cedar and Ortega rivers, Florida, USA, using geographic information system (GIS)-based kriging analyses and field measurements. Kriging analysis showed that two areas, one from the Cedar River area and the other from the northern end of the Ortega River area, were contaminated. The maximum concentrations of total chlordane, gamma-chlordane, alpha-chlordane, and trans-nonachlor in the sediments were, respectively, 101.8, 20.1, 26.3, and 19.2 microg/kg. A plot of total organic carbon (TOC)-normalized chlordane concentrations showed that effects of grain size on sediment chlordane contamination were negligible. A principal axis analysis further revealed that a linear correlation existed between alpha-chlordane and total chlordane as well as between gamma-chlordane and total chlordane, whereas no correlation existed between trans-nonachlor and total chlordane. Comparison of total chlordane concentration with Florida Sediment Assessment Guidelines showed that the Cedar River and the northern end of the Ortega River had total chlordane concentrations above the probable effect level (4.79 microg/kg), which could pose a potential risk to aquatic life.     

GIS-based spatial regression and prediction of water quality in river networks: A case study in Iowa

Nonpoint source pollution is the leading cause of the U.S.’s water quality problems. One important component of nonpoint source pollution control is an understanding of what and how watershed-scale conditions influence ambient water quality. This paper investigated the use of spatial regression to evaluate the impacts of watershed characteristics on stream NO₃NO₂-N concentration in the Cedar River Watershed, Iowa. An Arc Hydro geodatabase was constructed to organize various datasets on the watershed. Spatial regression models were developed to evaluate the impacts of watershed characteristics on stream NO₃NO₂-N concentration and predict NO₃NO₂-N concentration at unmonitored locations. Unlike the traditional ordinary least square (OLS) method, the spatial regression method incorporates the potential spatial correlation among the observations in its coefficient estimation. Study results show that NO₃NO₂-N observations in the Cedar River Watershed are spatially correlated, and by ignoring the spatial correlation, the OLS method tends to over-estimate the impacts of watershed characteristics on stream NO₃NO₂-N concentration. In conjunction with kriging, the spatial regression method not only makes better stream NO₃NO₂-N concentration predictions than the OLS method, but also gives estimates of the uncertainty of the predictions, which provides useful information for optimizing the design of stream monitoring network. It is a promising tool for better managing and controlling nonpoint source pollution.     

SURFACE WATER RESOURCES OF ST. JOHNS RIVER,FLORIDA1

ABSTRACT: The St. Johns River basin is the largest watershed entirely within the State of Florida, and is one of the few northward flowing rivers in the United States. The river basin contains 11,431 square miles, of which 9,430 square miles are drained by the river and its tributaries. The remainder drains into the Atlantic Ocean or the Intracoastal Waterway. Its largest sub-basin is the Oklawaha River basin, which has a drainage area of 2,870 square miles. Ground elevations range from sea level to 200 feet above mean sea level in the main river basin and as high as 300 feet above mean sea level in the Oklawaha River basin. This study was designed to investigate the surface water resources of the St. Johns River and the existing consumptive uses. The analysis revealed that the river is an extremely large and valuable resource which has been under-utilized and could play a much stronger role in serving the needs of the people in the basin.     

HYSTAR Sediment Model: Distributed Two‐Dimensional Simulation of Watershed Erosion and Sediment Transport Using Time‐Area Routing

An erosion and sediment transport component incorporated in the HYdrology Simulation using Time‐ARea method (HYSTAR) upland watershed model provides grid‐based prediction of erosion, transport and deposition of sediment in a dynamic, continuous, and fully distributed framework. The model represents the spatiotemporally varied flow in sediment transport simulation by coupling the time‐area routing method and sediment transport capacity approach within a grid‐based spatial data model. This avoids the common, and simplistic, approach of using the Universal Soil Loss Equation (USLE) to estimate erosion rates with a delivery ratio to relate gross soil erosion to sediment yield of a watershed, while enabling us to simulate two‐dimensional sediment transport processes without the complexity of numerical solution of the partial differential governing equations. In using the time‐area method for routing sediment, the model offers a novel alternative to watershed‐scale sediment transport simulation that provides detailed spatial representation. In predicting four‐year sediment hydrographs of a watershed in Virginia, the model provided good performance with R² of 0.82 and 0.78 and relative error of ?35% and 11% using the Yalin and Yang's sediment transport capacity equations, respectively. Prediction of spatiotemporal variation in sediment transport processes was evaluated using maps of sediment transport rates, concentrations, and erosion and deposition mass, which compare well with expected behavior of flow hydraulics and sediment transport processes.     

Investigating the Fate and Transport of Escherichia coli in the Charles River,Boston, Using High‐Resolution Observation and Modeling1

Abstract: The processes affecting the fate and transport of Escherichia coli in surface waters were investigated using high‐resolution observation and modeling. The concentration patterns in Boston’s Charles River were observed during four sampling events with a total of 757 samples, including two spatial surveys with two along‐river (1,500 m length) and three across‐river (600 m length) transects at approximately 25‐m intervals, and two temporal surveys at a fixed location (Community Boating) over seven days at hourly intervals. The data reveal significant spatial and temporal structure at scales not resolved by typical monitoring programs. A mechanistic, time‐variable, three‐dimensional coupled hydrodynamic and water quality model was developed using the ECOMSED and RCA modeling frameworks. The computational grid consists of 3,066 grid cells with average length dimension of 25 m. Forcing functions include upstream and downstream boundary conditions, Stony Brook, and Muddy River (major tributaries) combined sewer overflow (CSO) and non‐CSO discharge and wind. The model generally reproduces the observed spatial and temporal patterns. This includes the presence and absence of a plume in the study area under similar loading, but different hydrodynamic conditions caused by operation of the New Charles River Dam (downstream) and wind. The model also correctly predicts an episode of high concentrations at the time‐series station following seven days of no rainfall. The model has an overall root mean square error (RMSE) of 250 CFU/100 ml and an error rate (above or below the USEPA‐recommended single sample criteria value of 235 CFU/100 ml) of 9.4%. At the time series station, the model has an RMSE of 370 CFU/100 ml and an error rate of 15%.     

ASSESSING CHANGES IN WATERSHED FLOW REGIMES WITH SPATIALLY EXPLICIT HYDRAULIC MODELS1

ABSTRACT: Urbanization, farming, and other watershed activities can significantly alter storm hydrographs and sediment erosion rates within a watershed. These changes routinely cause severe economic and ecological problems manifested in the form of increased flooding and significant changes in channel morphology. As the activities within a watershed influence the hydrologic, hydraulic, and ecological conditions within a river, interdisciplinary approaches to predict and assess the impacts that different land uses have on streams need to be developed. An important component of this process is ascertaining how hydrologic changes induced by specific watershed activities will affect hydraulic conditions and the accompanying flood levels, sediment transport rates, and habitat conditions within a stream. A conceptual model for using spatially explicit (two‐dimensional) hydraulic models to help evaluate the impacts that changes in flow regime might have on a river is presented. This framework proposes that reproducing and quantifying flow complexity allows one to compare the hydraulic conditions within urban, urbanizing, and non‐urban streams in a more biologically and economically meaningful way. The justification, advantage, and need for such a method is argued through the results of one‐ and two‐dimensional hydraulic model studies. The implementation of this methodology in watershed urbanization studies is described.     

SIMULATING SEDIMENT TRANSPORT IN THE CARSON RIVER AND LAHONTAN RESERVOIR,NEVADA, USA1

ABSTRACT: The discovery of the Comstock Lode in western Nevada in 1859 initiated the use of liquid mercury (Hg) or “quicksilver” to remove gold (Au) and silver (Ag) from crushed ore. Today, Hg is present in historic mill tailings piles, in alluvial deposits adjacent to the Carson River, and in Lahontan Reservoir. Mercury concentrations in Carson River water have been reported as high as 61 μg/L by the U.S. Geological Survey. Fish from Lahontan Reservoir have methylmorcury (MeHg) concentrations as much as four times the 1.0 μg/g limit for human consumption. Since more than 95 percent of total Hg in water can be associated with particulates, the transport of sediment must be quantified to understand the fate of Hg in the system. By linking the RIVMOD hydrodynamic model with the WASP5 water quality model, and using suspended sediment rating curves along with bedload transport equations, reliable predictions of sediment transport can be made. Measured suspended sediment data from the Carson River, and an estimate of annual sediment loading to Lahontan Reservoir, were used to create a calibrated sediment transport model. Model simulations predicted the long term transport of sediment into Lahontan Reservoir, the transport of sediment into Lahontan Reservoir during a flood year (1986 water year), and concentrations of total Hg in the Carson River using an estimate of sediment Hg concentrations. This research will eventually be used with an Hg model to predict the fate of Hg in the river and reservoir system.     

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.     

Trace Metal Concentrations in an Intensive Agricultural Watershed in British Columbia,Canada1

Abstract: Effects of agricultural intensification and a naturally occurring landslide of asbestos material upon water and sediment quality in a transboundary watershed were investigated. The water and sediments of the Sumas River watershed were analyzed for copper (Cu), zinc (Zn), chromium (Cr), and nickel (Ni) concentrations in 1993/1994 and 2003/2004 and differences within sites over time were examined. Based upon a review of the literature, Cu and Zn were used as indicators of agricultural impacts while Cr and Ni were used as indicators of impacts from an asbestos landslide. Animal unit equivalents (AUEs) were calculated on a per area basis as an indicator of livestock density using detailed statistical census data. Whatman #42 filtered metals, bioavailable metals, and sediment‐bound metals (in the <63 μm fraction) were determined at 22 sites along the mainstem and tributaries, including two reference sites. Temperature, pH, and dissolved organic carbon (DOC) were also measured. The bioavailable metal fraction was determined using the diffusive gradient thin film technique (DGT). Sediment‐bound results were compared with British Columbia’s Interim Sediment Quality Guidelines (ISQGs) and Severe Effects Levels (SELs). A Wilcoxon signed rank test was used to determine if the concentrations of metals changed significantly within sites between 1993/1994 and 2003/2004. Spearman rank correlation analysis was used to determine relationships between trace metals, water quality parameters, and AUEs/hectare. The results indicate that Cu and Zn levels in sediments increased significantly to concentrations above the ISQGs of 35.7 mg/kg and 123 mg/kg, respectively from 1993/1994 to 2003/2004 in streams, where associated land use was dominated by intensive agriculture. Higher AUEs/hectare were significantly correlated with greater bioavailable levels of Zn as well as higher sediment‐bound Zn concentrations. Neither Cu nor Cr were detected by the DGTs on any of the sampling occasions. The Cr and Ni sediment concentrations were highest in Swift Creek, the headwater tributary affected by the natural landslide of asbestos material, and decreased in the Sumas River downstream from the point of input. Cr and Ni concentrations have increased in the mid‐region of the Sumas River since 1993/1994, suggesting downstream movement of the asbestos material over time. DGT results indicated that bioavailable Zn is significantly positively correlated to sediment‐bound Zn and livestock density, and bioavailable Ni is significantly correlated to sediment‐bound Ni.     

流域污染物通量测算方法研究

流域水系内污染物通量不仅能够用于评价各类污染源的水污染物入河负荷,也是对流域污染特征,水污染物在河流水体中复杂迁移、转化过程的最直观反应。准确测算流域水系内污染物跨界通量及其时空分布是进行流域水环境风险预警和风险管理的重要前提之一。针对目前多种污染物通量测算方法在进行污染物年通量估算时,结果不确定性大这一突出问题,以流域水质监测站年内逐日流量、悬浮颗粒物监测数据作为悬浮颗粒物年通量参考值,基于以月、半月、周为周期的监测策略,将逐日流量、悬浮颗粒物监测数据重新筛选抽样构造,由此,系统分析了不同流域集水面积、污染通量监测频次和目前常用通量估算方法对污染物年通量估算不确定性的影响。所得方法和结论可为进一步制定流域污染物通量的测算规范提供方法指引和技术支持。     

River Channel Geometry and Rating Curve Estimation Using Height above the Nearest Drainage

River channel geometry is an important input to hydraulic and hydrologic models. Traditional approaches to quantify river geometry have involved surveyed river cross sections, which cannot be extended to ungaged basins. In this paper, we describe a method for developing a synthetic rating curve to relate flow to water level in a stream reach based on reach‐averaged channel geometry properties developed using the Height above Nearest Drainage (HAND) method. HAND uses a digital elevation model (DEM) of the terrain and computes the elevation difference between each land surface cell and the stream bed cell to which it drains. Taking increments in water level in the stream, HAND defines the inundation zone and a water depth grid within this zone, and the channel characteristics are defined from this water depth grid. We apply our method to the Blanco River (Texas) and the Tar River (North Carolina) using 10‐m terrain data from the United States Geological Survey (USGS) 3D Elevation Program (3DEP) dataset. We evaluate the method's performance by comparing the reach‐average stage‐river geometry relationships and rating curves to those from calibrated Hydrologic Engineering Center's River Analysis System (HEC‐RAS) models and USGS gage observations. The results demonstrate that after some adjustment, the river geometry information and rating curves derived from HAND using national‐coverage datasets are comparable to those obtained from hydraulic models or gage measurements. We evaluate the inundation extent and show our approach is able to capture the majority of the Federal Emergency Management Agency (FEMA) 100‐year floodplain.     

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

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

Modeling the Highly Dynamic Loading of Mercury Species in the Carson River and Lahontan Reservoir System,Nevada

The semiarid Carson River — Lahontan Reservoir system in Nevada, United States is highly contaminated with mercury (Hg) from historic mining with contamination dispersed throughout channel and floodplain deposits. Work builds on previous research using a fully dynamic numerical model to outline a complete conceptualization of the system that includes transport and fate of both sorbed and dissolved constituents. Flow regimes are defined to capture significant mechanisms of Hg loading that include diffusion, channel pore water advective flux, bank erosion, and overbank deposition. Advective flux of pore water is required to reduce dilution and likely represents colloidal‐mediated transport. Fluvial concentrations span several orders of magnitude with spatial and temporal trends simulated within 10‐24% error for all modeled species. Over the simulation period, 1991‐2008, simulated loads are 582 kg/yr (THg²⁺), 4.72 kg/yr (DHg²⁺), 0.54 kg/yr (TMeHg), and 0.07 kg/yr (DMeHg) with bank erosion processes the principal mechanism of loading for both total and dissolved species. Prediction error in the reservoir is within one‐order of magnitude and considered qualitative; however, simulated results indicate internal cycling within the receiving reservoir accounts for only 1% of the reservoir's water column contamination, with river channel sediment sources more influential in the upper reservoir and bank erosion processes having greater influence in the lower reservoir.     

Freshwater Mussel Population Status and Habitat Quality in the Clinch River,Virginia and Tennessee,USA: A Featured Collection

The Clinch River of southwestern Virginia and northeastern Tennessee is arguably the most important river for freshwater mussel conservation in the United States. This featured collection presents investigations of mussel population status and habitat quality in the Clinch River. Analyses of historic water‐ and sediment‐quality data suggest that water column ammonia and water column and sediment metals, including Cu and Zn, may have contributed historically to declining densities and extirpations of mussels in the river's Virginia sections. These studies also reveal increasing temporal trends for dissolved solids concentrations throughout much of the river's extent. Current mussel abundance patterns do not correspond spatially with physical habitat quality, but they do correspond with specific conductance, dissolved major ions, and water column metals, suggesting these and/or associated constituents as factors contributing to mussel declines. Mussels are sensitive to metals. Native mussels and hatchery‐raised mussels held in cages in situ accumulated metals in their body tissues in river sections where mussels are declining. Organic compound and bed‐sediment contaminant analyses did not reveal spatial correspondences with mussel status metrics, although potentially toxic levels were found. Collectively, these studies identify major ions and metals as water‐ and sediment‐quality concerns for mussel conservation in the Clinch River.     

WATER QUALITY MODEL OUTPUT UNCERTAINTY AS AFFECTED BY SPATIAL RESOLUTION OF INPUT DATA1

ABSTRACT: Resolution of the input GIS data used to parameterize distributed‐parameter hydrologic/water quality models may affect uncertainty in model outputs and impact the subsequent application of model results in watershed management. In this study we evaluated the impact of varying spatial resolutions of DEM, land use, and soil data (30 × 30 m, 100 × 100 m, 150 × 150 m, 200 × 200 m, 300 × 300 m, 500 × 500 m, and 1,000 × 1,000 m) on the uncertainty of SWAT predicted flow, sediment, NO₃‐N, and TP transport. Inputs included measured hydrologic, meteorological, and watershed characteristics as well as water quality data from the Moores Creek watershed in Washington County, Arkansas. The SWAT model output was most affected by input DEM data resolution. A coarser DEM data resolution resulted in decreased representation of watershed area and slope and increased slope length. Distribution of pasture, forest, and urban areas within the watershed was significantly affected at coarser resolution of land use and resulted in significant uncertainty in predicted sediment, NO₃‐N, and TP output. Soils data resolution had no significant effect on flow and NO₃‐N predictions; however, sediment was overpredicted by 26 percent, and TP was underpredicted by 26 percent at 1,000 m resolution. This may be due to change in relative distribution of various hydrologic soils groups (HSGs) in the watershed. Minimum resolution for input GIS data to achieve less than 10 percent model output error depended upon the output variable of interest. For flow, sediment, NO₃‐N, and TP predictions, minimum DEM data resolution should range from 30 to 300 m, whereas minimum land use and soils data resolution should range from 300 to 500 m.     

INVESTIGATION OF TURBIDITY MAXIMUM IN A MESOTIDAL ESTUARY,TAIWAN1

ABSTRACT: To comprehend the distributions of salinity, temperature, and suspended sediment in the Danshuei River estuary in Taiwan, monthly field surveys were conducted in 2003. These included several high and low slackwater surveys and intensive surveys. The results show that the Danshuei River estuary is predominately a partially mixed estuary. The highest concentration of suspended sediment is typically observed at the Chung‐Hsin Bridge, the most upstream sampling station. The suspended sediment concentration exhibits a general decreasing trend in the downstream direction. It may be concluded that the sediments mostly come from the upstream reach. A locally high concentration of suspended sediment is found at the Kuan‐Du station because of the local deep channel bathymetry and two‐layered estuarine circulation. A vertical two‐dimensional hydrodynamic and sediment transport model is applied to investigate the tidally averaged salinity distribution, residual circulation, and suspended sediment concentration. The modeling results reveal that, under the Q₇₅ flow condition (i.e., low flow), a turbidity maximum occurs at the Kuan‐Du station due to the strong estuarine circulation. The model simulation with a much higher river flow condition results in a weaker residual circulation and weaker turbidity maximum.     

Effects of Impervious Cover at Multiple Spatial Scales on Coastal Watershed Streams1

Abstract: The spatial scale and location of land whose development has the strongest influence on aquatic ecosystems must be known to support land use decisions that protect water resources in urbanizing watersheds. We explored impacts of urbanization on streams in the West River watershed, New Haven, Connecticut, to identify the spatial scale of watershed imperviousness that was most strongly related to water chemistry, macroinvertebrates, and physical habitat. A multiparameter water quality index was used to characterize regional urban nonpoint source pollution levels. We identified a critical level of 5% impervious cover, above which stream health declined. Conditions declined with increasing imperviousness and leveled off in a constant state of impairment at 10%. Instream variables were most correlated (0.77 ≤ |r| ≤ 0.92, p < 0.0125) to total impervious area (TIA) in the 100‐m buffer of local contributing areas (～5‐km² drainage area immediately upstream of each study site). Water and habitat quality had a relatively consistent strong relationship with TIA across each of the spatial scales of investigation, whereas macroinvertebrate metrics produced noticeably weaker relationships at the larger scales. Our findings illustrate the need for multiscale watershed management of aquatic ecosystems in small streams flowing through the spatial hierarchies that comprise watersheds with forest‐urban land use gradients.     

CO-KRIGING TO INCORPORATE SCREENING DATA: HUDSON RIVER SEDIMENT PCBs1

ABSTRACT: Kriging methods of geostatistical analysis provide valuable techniques for analysis of sediment contamination problems, including interpolation of concentration maps from point data and estimation of global mean concentrations. Sample collection efforts frequently include preliminary screening data of considerably more extensive coverage than the laboratory analyses on which estimation is usually based. How should these be incorporated in kriging? Screening and laboratory analysis constitute two separate estimates of the same spatial field but of very different characteristics. A modified version of co-kriging is developed to include the imprecise screening information in the analysis of contaminant distribution. Use of the method is demonstrated on a data set of sediment PCB samples from the Upper Hudson River, for which preliminary categorical mass spectrometry screening was used to select a smaller set of samples for gas chromatograph analysis. The method is widely applicable to many situations of contaminant and natural resource estimation.