Fate and movement of atrazine, cyanazine, metolachlor and selected degradation products in water resources of the deep Loess Hills of Southwestern Iowa, USA

The environmental fate and movement of herbicides widely used for weed control in corn are assessed for a deep loess soil in southwestern Iowa. Beginning in the early 1980s, the herbicide-based weed control program emphasized the application of atrazine (ATR) or cyanazine (CYN) and metolachlor (MET) for both broadleaf and grass control. Between 1992 and 1995, concentrations of ATR, desethylatrazine (DEA), desisopropylatrazine (DIA), CYN and MET were measured in rainwater, both shallow and deep vadose zone water, and well water. Results show that the frequency of herbicide detections and the range and distribution of occurrences are dependent upon both landscape position and temporal inputs of recharge water from rainfall. Generally, DIA was observed more frequently and in higher mean concentration in well water than DEA, while DEA was observed more frequently than DIA in vadose zone groundwater. A chromatographic analogy is suggested to explain the occurrence patterns observed for both parent herbicide and degradation products within the unsaturated zone water. Analysis of rainwater samples collected during this time also revealed low concentrations of ATR, CYN and MET, with the timing of the detections indicative of non-local transport. Results show that the deep loess soil conducts both water and agricultural chemicals relatively rapidly and as such represents a production system which is vulnerable to contamination of shallow groundwater by herbicide-derived chemicals. Results also illustrate the importance of including major herbicide degradation products in water resource impact assessment studies.     

STUMOD—a Tool for Predicting Fate and Transport of Nitrogen in Soil Treatment Units

A typical onsite wastewater treatment system consists of a septic tank and a soil treatment unit to treat wastewater before it is discharged through the vadose zone to an aquifer. A tool was developed for the purpose of predicting the fate and transport of nitrogen in soil treatment units (STUMOD or Soil Treatment Unit Model). STUMOD calculates nitrogen species concentrations and the fraction of total nitrogen reaching the aquifer or a specified soil depth. Input data include parameters for hydraulics and nutrient transport and transformation. An analytical solution is used to calculate the profile of pressure based on Darcy’s equation and the relationships between suction head, unsaturated hydraulic conductivity, and soil moisture. Chemical transport is based on simplification of the advection–dispersion equation. STUMOD is relatively simple to use but accounts for important processes such as ammonium sorption, nitrification, and denitrification. STUMOD accounts for the effect of soil moisture content (a surrogate for redox conditions) on nitrification and denitrification reactions. The model has provisions to handle the influence of temperature and organic carbon content on nitrogen transformation. Model outputs, generated based on input parameters obtained from extensive literature review, were compared to a numerical model and data from laboratory tests and field sites. Both measured data and STUMOD outputs show a relatively higher removal in clayey soils compared to sandy soils. Consistent with literature data for most soils, STUMOD predicted ammonium conversion to nitrate within the first foot below the trench infiltrative surface.     

Hurricane Katrina-induced forest damage in relation to ecological factors at landscape scale

Forest stand stability to strong winds such as hurricanes has been found to be associated with a number of forest, soil and topography factors. In this study, through applying geographic information system (GIS) and logit regression, we assessed effects of forest characteristics and site conditions on pattern, severity and probability of Hurricane Katrina disturbance to forests in the Lower Pearl River Valley, USA. The factors included forest type, forest coverage, stand density, soil great group, elevation, slope, aspect, and stream buffer zone. Results showed that Hurricane Katrina damaged 60% of the total forested land in the region. The distribution and intensity of the hurricane disturbance varied across the landscape, with the bottomland hardwood forests on river floodplains most severely affected. All these factors had a variety of effects on vulnerability of the forests to the hurricane disturbance and thereby spatial patterns of the disturbance. Soil groups and stand factors including forest types, forest coverage and stand density contributed to 85% of accuracy in modeling the probability of the hurricane disturbance to forests in this region. Besides assessment of Katrina's damage, this study elucidates the great usefulness of remote sensing and GIS techniques combined with statistics modeling in assessment of large-scale risks of hurricane damage to coastal forests.     

Relationships between watershed emergy flow and coastal New England salt marsh structure, function, and condition

This study evaluated the link between watershed activities and salt marsh structure, function, and condition using spatial emergy flow density (areal empower density) in the watershed and field data from 10 tidal salt marshes in Narragansett Bay, RI, USA. The field-collected data were obtained during several years of vegetation, invertebrate, soil, and water quality sampling. The use of emergy as an accounting mechanism allowed disparate factors (e.g., the amount of building construction and the consumption of electricity) to be combined into a single landscape index while retaining a uniform quantitative definition of the intensity of landscape development. It expanded upon typical land use percentage studies by weighting each category for the intensity of development. At the RI salt marsh sites, an impact index (watershed emergy flow normalized for marsh area) showed significant correlations with mudflat infauna species richness, mussel density, plant species richness, the extent and density of dominant plant species, and denitrification potential within the high salt marsh. Over the 4-year period examined, a loading index (watershed emergy flow normalized for watershed area) showed significant correlations with nitrite and nitrate concentrations, as well as with the nitrogen to phosphorus ratios in stream discharge into the marshes. Both the emergy impact and loading indices were significantly correlated with a salt marsh condition index derived from intensive field-based assessments. Comparison of the emergy indices to calculated nitrogen loading estimates for each watershed also produced significant positive correlations. These results suggest that watershed emergy flow is a robust index of human disturbance and a potential tool for rapid assessment of coastal wetland condition.     

新疆典型城-郊-乡梯度带土壤微量元素空间分布与来源

从新疆某地典型城-郊-乡梯度带采集了77个表层(0~20 cm)土壤样品,基于GIS技术与多元统计分析方法,研究各梯度带表层土壤中Hg、As、V、Co、Ni、Cu、Zn、Cd、Pb和Sb等10种微量元素的空间分布特征与主要来源。结果表明:Hg元素在城区、郊区和乡村表层土壤中的平均含量均超出研究区土壤背景值,As元素在城区和郊区表层土壤中的平均含量超出背景值,Zn和Pb元素在城区表层土壤中的平均含量超出背景值,其余元素在3个梯度带表层土壤中的平均含量均低于相应的背景值。研究区表层土壤中,V、Co、Ni、Zn、Cd、Pb和Sb等7种元素的空间分布格局基本相似,均呈现沿城市化梯度带分布的地带性格局;As、Cu和Hg等3种元素的空间分布呈现岛状格局。来源分析结果表明,各梯度带表层土壤中的微量元素的来源各不相同。     

Effects of an Oil Spill on Soil Physico-Chemical Properties of a Spill Site in a Typic Udipsamment of the Niger Delta Basin of Nigeria

Physico-chemical analysis of soil samples at an oil spill site in a Typic Udipsamment of the Niger Delta Basin of Nigeria showed that the total hydrocarbon content of top soil layers ranged from 18.6 to 23.6 ppm in the heavy impact zone and the oil had penetrated to a depth op 8.4 m. The concentration of hydrocarbons in the medium impact zone ranged from 10.04 to 10.38 ppm while hydrocarbons were not detected in 85% of samples from the unimpacted reference zone. Heavy metal concentration measurements in the soil revealed a significant build-up (P < 0.05) of lead, copper and zinc in the heavy impact zone. Other quality parameters including electrical conductity, exchangeable cations, total nitrogen and available phosphorus in impacted soils were relatively low, while the total organic carbon was high compared with the reference site. Textural class of soil from the different depths showed a predominantly brown sand at the topsoil, loamy sand and grey coarse sand at medium depths, and grey coarse sand and greyish sandy clay at greater depths.     

Effects of an oil spill on soil physico-chemical properties of a spill site in the Niger Delta Area of Nigeria

Physico-chemical analysis of soil samples at an oil spill site in the Niger Delta Area of Nigeria showed that the total hydrocarbon content of top soil layers ranged from 0.8 to 12.4 ppm in the heavy impact zone and the oil had penetrated to a depth of 7.2 m. Hydrocarbon concentration in the medium impact zone ranged from 0.02 to 0.40 ppm while hydrocarbons were not detected in 75% of samples from the unimpacted reference zone. Measurement of heavy metal concentrations in the soils revealed a significant build-up (p<0.05) of lead, iron and zinc in the heavy impact zone. Other parameters including electrical conductivity, exchangeable cations, available phosphorus and total nitrogen in impacted soils were comparatively low while the total organic carbon was high, compared with the reference site. Textural class of soil from the different depths showed a predominantly brown sand at the top soil, loamy sand at medium depths and grey coarse sand at greater depths.     

Effects of geomorphology, habitat, and spatial location on fish assemblages in a watershed in Ohio, USA

In this paper, we evaluate relationships between in-stream habitat, water chemistry, spatial distribution within a predominantly agricultural Midwestern watershed and geomorphic features and fish assemblage attributes and abundances. Our specific objectives were to: (1) identify and quantify key environmental variables at reach and system wide (watershed) scales; and (2) evaluate the relative influence of those environmental factors in structuring and explaining fish assemblage attributes at reach scales to help prioritize stream monitoring efforts and better incorporate all factors that influence aquatic biology in watershed management programs. The original combined data set consisted of 31 variables measured at 32 sites, which was reduced to 9 variables through correlation and linear regression analysis: stream order, percent wooded riparian zone, drainage area, in-stream cover quality, substrate quality, gradient, cross-sectional area, width of the flood prone area, and average substrate size. Canonical correspondence analysis (CCA) and variance partitioning were used to relate environmental variables to fish species abundance and assemblage attributes. Fish assemblages and abundances were explained best by stream size, gradient, substrate size and quality, and percent wooded riparian zone. Further data are needed to investigate why water chemistry variables had insignificant relationships with IBI scores. Results suggest that more quantifiable variables and consideration of spatial location of a stream reach within a watershed system should be standard data incorporated into stream monitoring programs to identify impairments that, while biologically limiting, are not fully captured or elucidated using current bioassessment methods.     

Monitoring and identification of spatiotemporal landscape changes in multiple remote sensing images by using a stratified conditional Latin hypercube sampling approach and geostatistical simulation

This study develops a stratified conditional Latin hypercube sampling (scLHS) approach for multiple, remotely sensed, normalized difference vegetation index (NDVI) images. The objective is to sample, monitor, and delineate spatiotemporal landscape changes, including spatial heterogeneity and variability, in a given area. The scLHS approach, which is based on the variance quadtree technique (VQT) and the conditional Latin hypercube sampling (cLHS) method, selects samples in order to delineate landscape changes from multiple NDVI images. The images are then mapped for calibration and validation by using sequential Gaussian simulation (SGS) with the scLHS selected samples. Spatial statistical results indicate that in terms of their statistical distribution, spatial distribution, and spatial variation, the statistics and variograms of the scLHS samples resemble those of multiple NDVI images more closely than those of cLHS and VQT samples. Moreover, the accuracy of simulated NDVI images based on SGS with scLHS samples is significantly better than that of simulated NDVI images based on SGS with cLHS samples and VQT samples, respectively. However, the proposed approach efficiently monitors the spatial characteristics of landscape changes, including the statistics, spatial variability, and heterogeneity of NDVI images. In addition, SGS with the scLHS samples effectively reproduces spatial patterns and landscape changes in multiple NDVI images.     

Spatial distribution characteristics of soil organic carbon in subtropical forests of mountain Lushan,China

The study on the spatial distribution of forest soil organic carbon (SOC) is of great significance for accurate assessment of carbon storage in forest ecosystems. In the present study, by taking eight kinds of forest soils of Mountain Lushan in the subtropical area as the research object, we studied the spatial distribution characteristics of SOC in this mountainous area. The results showed that the SOC content and SOC density (SOCD) of main forest types in the Mountain Lushan were lower than the national and the world average. The soil layer of Lushan forest was thinner, and the SOC and active SOC (ASOC) contents of different forest types and SOCDs are the highest in the surface soil. SOCD of the topsoil accounts for 32.64–54.03% of the total SOCD in the whole soil profile. Surface litter is an important source of SOC, and the different vegetation types are the important reason for the different spatial distribution of SOC in this area. Soil SOC contents in the high-altitude forest (bamboo forest, deciduous broadleaf forest, Pinus taiwanensis forest, evergreen-deciduous forest, and coniferous-broadleaved mixed forest) were higher than those in the low-altitude forest (evergreen broadleaf forest, shrub, and Pinus massoniana forest). However, the difference in SOC content exhibited at the altitude gradient is significantly lower than that in SOC in the soil profile. This indicates that both soil depth and elevation are the important factors that affected SOC distribution. However, the influence of soil depth on spatial distribution of SOC may be more complex than that of altitude. Vegetation types and soil properties are the main reasons for the large differences of reduction rate in the contents of SOC and ASOC.     

基于污水厂监测数据的反硝化脱氮影响因素分析

以污水厂长序列运行数据为基础,反推脱氮实际影响因素,开展反硝化脱氮试验.结果表明,在2019年4月前后TN出水平均浓度降低58％;受季节和污水厂升级改造影响,温度由15℃升高至25℃,混合液悬浮固体质量浓度(MLSS)由4300 mg/L增至8000 mg/L,活性污泥成分增加,脱氮效果提升;BOD5/TN和pH值等变...     

The use of spatial empirical models to estimate soil erosion in arid ecosystems

The central objective of this project was to utilize geographical information systems and remote sensing to compare soil erosion models, including Modified Pacific South-west Inter Agency Committee (MPSIAC), Erosion Potential Method (EPM), and Revised Universal Soil Loss Equation (RUSLE), and to determine their applicability for arid regions such as Kuwait. The northern portion of Umm Nigga, containing both coastal and desert ecosystems, falls within the boundaries of the de-militarized zone (DMZ) adjacent to Iraq and has been fenced off to restrict public access since 1994. Results showed that the MPSIAC and EPM models were similar in spatial distribution of erosion, though the MPSIAC had a more realistic spatial distribution of erosion and presented finer level details. The RUSLE presented unrealistic results. We then predicted the amount of soil loss between coastal and desert areas and fenced and unfenced sites for each model. In the MPSIAC and EPM models, soil loss was different between fenced and unfenced sites at the desert areas, which was higher at the unfenced due to the low vegetation cover. The overall results implied that vegetation cover played an important role in reducing soil erosion and that fencing is much more important in the desert ecosystems to protect against human activities such as overgrazing. We conclude that the MPSIAC model is best for predicting soil erosion for arid regions such as Kuwait. We also recommend the integration of field-based experiments with lab-based spatial analysis and modeling in future research.     

Spatial and temporal variation of moisture content in the soil profiles of two different agricultural fields of semi-arid region

Modeling spatio-temporal variation of soil moisture with depth in the soil profile plays an important role for semi-arid crop production from an agro-hydrological perspective. This study was performed in Guvenc Catchment. Two soil series that were called Tabyabayir (TaS) and Kervanpinari (KeS) and classified as Leptosol and Vertisol Soil Groups were used in this research. The TeS has a much shallower (0–34 cm) than the KeS (0–134 cm). At every sampling time, a total of geo-referenced 100 soil moisture samples were taken based on horizon depths. The results indicated that soil moisture content changed spatially and temporally with soil texture and profile depth significantly. In addition, land use was to be important factor when soil was shallow. When the soil conditions were towards to dry, higher values for the coefficient of variation (CV) were observed for TaS (58 and 43 % for A and C horizons, respectively); however, the profile CV values were rather stable at the KeS. Spatial variability range of TaS was always higher at both dry and wet soil conditions when compared to that of KeS. Excessive drying of soil prevented to describe any spatial model for surface horizon, additionally resulting in a high nugget variance in the subsurface horizon for the TaS. On the contrary to TaS, distribution maps were formed all horizons for the KeS at any measurement times. These maps, depicting both dry and wet soil conditions through the profile depth, are highly expected to reduce the uncertainty associated with spatially and temporally determining the hydraulic responses of the catchment soils.     

Probability-based nitrate contamination map of groundwater in Kinmen

Groundwater supplies over 50 % of drinking water in Kinmen. Approximately 16.8 % of groundwater samples in Kinmen exceed the drinking water quality standard (DWQS) of NO₃ ^?-N (10 mg/L). The residents drinking high nitrate-polluted groundwater pose a potential risk to health. To formulate effective water quality management plan and assure a safe drinking water in Kinmen, the detailed spatial distribution of nitrate–N in groundwater is a prerequisite. The aim of this study is to develop an efficient scheme for evaluating spatial distribution of nitrate–N in residential well water using logistic regression (LR) model. A probability-based nitrate–N contamination map in Kinmen is constructed. The LR model predicted the binary occurrence probability of groundwater nitrate–N concentrations exceeding DWQS by simple measurement variables as independent variables, including sampling season, soil type, water table depth, pH, EC, DO, and Eh. The analyzed results reveal that three statistically significant explanatory variables, soil type, pH, and EC, are selected for the forward stepwise LR analysis. The total ratio of correct classification reaches 92.7 %. The highest probability of nitrate–N contamination map presents in the central zone, indicating that groundwater in the central zone should not be used for drinking purposes. Furthermore, a handy EC–pH-probability curve of nitrate–N exceeding the threshold of DWQS was developed. This curve can be used for preliminary screening of nitrate–N contamination in Kinmen groundwater. This study recommended that the local agency should implement the best management practice strategies to control nonpoint nitrogen sources and carry out a systematic monitoring of groundwater quality in residential wells of the high nitrate–N contamination zones.     

Implication of two in-stream processes in the fate of nutrients discharged by sewage system into a temporary river

The aim of this study was to better understand the fate of nutrients discharged by sewage treatment plants into an intermittent Mediterranean river, during a low-flow period. Many pollutants stored in the riverbed during the low-flow period can be transferred to the downstream environments during flood events. The study focused on two processes that affect the fate and the transport of nutrients, a physical process (retention in the riverbed sediments) and a biological process (denitrification). A spatial campaign was carried out during a low-flow period to characterize the nutrient contents of both water and sediments in the Vène River. The results showed high nutrient concentrations in the water column downstream of the treated wastewater disposal (up to 13,315 μg N/L for ammonium and 2,901 μg P/L for total phosphorus). Nutrient concentrations decreased rapidly downstream of the disposal whereas nutrient contents in the sediments increased (up to 1,898 and 784 μg/g for total phosphorus and Kjeldahl nitrogen, respectively). According to an in situ experiment using sediment boxes placed in the riverbed for 85 days, we estimated that the proportion of nutrients trapped in the sediments represents 25% (respectively 10%) of phosphorus (respectively nitrogen) loads lost from the water column. In parallel, laboratory tests indicated that denitrification occurred in the Vène River, and we estimated that denitrification likely coupled to nitrification processes during the 85 days of the experiment was significantly involved in the removal of nitrogen loads (up to 38%) from the water column and was greater than accumulation processes.     

An ecological function and services approach to total maximum daily load (TMDL) prioritization

Prioritizing total maximum daily load (TMDL) development starts by considering the scope and severity of water pollution and risks to public health and aquatic life. Methodology using quantitative assessments of in-stream water quality is appropriate and effective for point source (PS) dominated discharge, but less so in watersheds with mostly nonpoint source (NPS) related impairments. For NPSs, prioritization in TMDL development and implementation of associated best management practices should focus on restoration of ecosystem physical functions, including how restoration effectiveness depends on design, maintenance and placement within the watershed. To refine the approach to TMDL development, regulators and stakeholders must first ask if the watershed, or ecosystem, is at risk of losing riparian or other ecologically based physical attributes and processes. If so, the next step is an assessment of the spatial arrangement of functionality with a focus on the at-risk areas that could be lost, or could, with some help, regain functions. Evaluating stream and wetland riparian function has advantages over the traditional means of water quality and biological assessments for NPS TMDL development. Understanding how an ecosystem functions enables stakeholders and regulators to determine the severity of problem(s), identify source(s) of impairment, and predict and avoid a decline in water quality. The Upper Reese River, Nevada, provides an example of water quality impairment caused by NPS pollution. In this river basin, stream and wetland riparian proper functioning condition (PFC) protocol, water quality data, and remote sensing imagery were used to identify sediment sources, transport, distribution, and its impact on water quality and aquatic resources. This study found that assessments of ecological function could be used to generate leading (early) indicators of water quality degradation for targeting pollution control measures, while traditional in-stream water quality monitoring lagged in response to the deterioration in ecological functions.     

Numerical Modeling on the Effect of Dissolved Oxygen on Nitrogen Transformation and Transport in Unsaturated Porous System

Nitrogen pollution in groundwater resulting from wastewater application to land is a common problem, and it causes a major threat to groundwater-based drinking water supplies. In this study, a numerical model is developed to study the nitrogen species transport and transformation in unsaturated porous media. Further, a new mass transfer module for dissolved oxygen (DO) is incorporated in the one-dimensional numerical model for nitrogen species transport to describe the fate and transport of nitrogen species, dissolved oxygen, dissolved organic carbon (DOC), and biomass. The spatial and temporal variation of dissolved oxygen is incorporated in the model through the mass transfer from gaseous phase to water phase in an unsaturated porous system. The numerical results of the water flow model and single species and multispecies transport model in an unsaturated zone developed for this purpose have been validated with the available analytical/numerical solution. The developed model is applied in clay loam, silt, and sand soils to analyze the transport behavior of nitrogen species under unsaturated condition. The numerical results suggest that the high rate of oxygen mass transfer from the air phase to the water phase positively increases the dissolved oxygen in the applied wastewater and enhances the nitrification process. Because of this high oxygen mass transfer, the nitrate nitrogen concentration significantly increases in the unsaturated zone and the same is transported to a larger depth at higher simulation period. On the other hand, the low rate of oxygen mass transfer implicitly enhances the denitrification process and finally reduces the nitrate nitrogen concentration in the unsaturated zone. The numerical results also show that the nitrate nitrogen transport is rapid in sandy soil when compared with clay loam and silty soils under high oxygen mass transfer rate. In essence, the high oxygen mass transfer rate significantly increases the nitrate nitrogen in the unsaturated zone, especially at a greater depth at larger time levels and eventually affects the groundwater quality.     

Testing differences in methods of preparing moss samples

The main objective of this research was to estimate the total mass of nitrogen discharged from various sources in Korea using the mass balance approach. Three different nitrogen mass balances were presented: (1) agricultural activities including raising crops and animal husbandry; (2) domestic activities, and (3) activities in forest and urban areas. These nitrogen balances were combined to estimate riverine discharge of nitrogen to the ocean in national scale. Nitrogen inputs include atmospheric deposition, biological nitrogen fixation, application of inorganic fertilizers/manures, animal feed/imported foodstuffs, and meat/fish. Nitrogen outputs include ammonia volatilization, denitrification, human/animal waste generation, crop/meat production, and riverine discharge to the ocean. The estimated total nitrogen input in Korea was 1,194.5 × 10³ tons N/year. Nitrogen discharged into rivers was estimated as 408–422 × 10³ tons N/year, of which 66–71% was diffuse in origin. The estimated diffuse discharges for land uses were estimated as 82 × 10³ tons N/year from agricultural areas, 7 × 10³ tons N/year from forestry and 75 × 10³ tons N/year from urban and industrial areas.     

Riparian vegetation assessment of Cauvery River Basin of South India

The Cauvery river basin of South India has a large phyto-floristic wealth, rightfully enough to constitute a separate phyto-geographic unit. Increase in the anthropogenic pressures within the river basin and surrounding landscapes have persistently stressed the riparian ecosystem structure adversely, besides altering its composition. The objective of this study was to examine the present status of riparian vegetation along the Cauvery river basin, in response to anthropogenic pressures. For vegetation analysis, the riparian forest coming in the middle stretch of Cauvery river was categorized into two zones, viz., forest zone covering ~54 km stretch and agroecosystem zone covering ~80 km stretch. In each zone, tree species were quantified using transects at 8-km interval. Overall tree species accounting for both forest and agroecosystem were recorded and compared. The results indicate that the mean density and basal area of trees per plot were higher in the forest zone than agroecosystem zone. The Shannon–Weiner diversity of forest zone is 4.6, which is higher than agroecosystem. In addition, species composition indicated a relatively low or poor similarity between the two zones. The vegetation density and site disturbance scores for all the study sites reveals that sand mining and grazing areas have exerted negative impact on riparian forest. The results of the present study clearly brought out the need for preparing and implementing site-specific conservation plans for riparian ecosystem.     

Drainage Effects on Stream Nitrate-N and Hydrology in South-Central Minnesota (USA)

Excessive nitrate-N in south-central Minnesota ditches and streams is related to land-use change, and may be contributing to the development of the zone of hypoxia in the Gulf of Mexico. Intensive land-use (agricultural management) has progressively increased as subsurface drainage has improved crop productivity over the past 25 years. We have examined water at varying scales for delta18O and, nitrate-N concentrations. Additionally, analysis of annual peak flows, and channel geomorphic features provided a measure of hydrologic change. Laboratory and field results indicate that agricultural drainage has influenced riverine source waters, concentrations of nitrate-N, channel dimensions and hydrology in the Blue Earth River (BER) Basin. At the mouth of the BER shallow ground water comprises the largest source water component. The highest nitrate-N concentrations in the BER and tributaries typically occurred in May and June and ranged from 7-34 mg L(-1). Peak flows for the 1.01-2-yr recurrence intervals increased by 20-to-206% over the past 25 years. Geomorphic data suggest that small channels (ditches) were entrenched by design, whereas, natural that are disconnected from an accessible riparian corridor. Frequent access to a functioning riparian zone is important for denitrification.