Assessing Natural and Anthropogenic Variability in Wetland Structure for Two Hydrogeomorphic Riverine Wetland Subclasses

The hydrogeomorphic approach (HGM) to wetland classification and functional assessment has been applied regionally throughout the United States, but the ability of HGM functional assessment models to reflect wetland condition has limited verification. Our objective was to determine how variability derived from anthropogenic effects and natural variability impacted site assessment variables within regional wetland subclasses in central Oklahoma. We collected data for nine potential assessment variables including vegetation physiognomy (e.g., tree basal area, herbaceous cover, canopy cover, etc.) and soil organic matter at wetlands of two HGM riverine subclasses (oxbow and riparian) in May and June, 2010. Using Akaike Information Criteria, we identified limited relationships between landscape disturbance metrics and assessment variables within subclasses. The high degree of natural variability from climatic and hydrologic factors within both subclasses may be masking the impact of landscape disturbance on the other measured assessment variables. Precipitation had significant effects on assessment variables within each of the subclasses. To reduce natural climatic variability, the reference domain may need to be further subdivided. The approach used in this study provides fairly rapid and quantitative methods for evaluating the effectiveness of using HGM assessment variables in assessing wetland condition regionally.     

Losing Function Through Wetland Mitigation in Central Pennsylvania,USA

In the United States, the Clean Water Act requires mitigation for wetlands that are negatively impacted by dredging and filling activities. During the mitigation process, there generally is little effort to assess function for mitigation sites and function is usually inferred based on vegetative cover and acreage. In our study, hydrogeomorphic (HGM) functional assessment models were used to compare predicted and potential levels of functional capacity in created and natural reference wetlands. HGM models assess potential function by measurement of a suite of structural variables and these modeled functions can then be compared to those in natural, reference wetlands. The created wetlands were built in a floodplain setting of a valley in central Pennsylvania to replace natural ridge-side slope wetlands. Functional assessment models indicated that the created sites differed significantly from natural wetlands that represented the impacted sites for seven of the ten functions assessed. This was expected because the created wetlands were located in a different geomorphic setting than the impacted sites, which would affect the type and degree of functions that occur. However, functional differences were still observed when the created sites were compared with a second set of reference wetlands that were located in a similar geomorphic setting (floodplain). Most of the differences observed in both comparisons were related to unnatural hydrologic regimes and to the characteristics of the surrounding landscape. As a result, the created wetlands are not fulfilling the criteria for successful wetland mitigation.     

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

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

Rapid Assessment of Urban Wetlands: Do Hydrogeomorphic Classification and Reference Criteria Work?

The Hydrogeomorphic (HGM) functional assessment method is predicated on the ability of hydrogeomorphic wetland classification and visual assessment of alteration to provide reference standards against which functions in individual wetlands can be evaluated. The effectiveness of this approach was tested by measuring nitrogen cycling functions in forested wetlands in an urbanized region in New Jersey, USA. Fourteen sites represented three HGM classes and were characterized as “least disturbed reference” or “non-reference” based on initial visual assessment. Water table levels and in situ rates of net nitrogen mineralization, net nitrification, and denitrification were measured over one year in each site. Hydrological alterations, resulting in consistently low or flashy water table levels, were not correlated with a priori designations as reference and non-reference. Although the flat-riverine wetland class had lower net nitrification and higher denitrification rates than riverine or mineral flat wetland classes, this difference was attributable to the lack of hydrologically-altered wetlands in the flat-riverine class, and thus more consistently wet conditions. Within all HGM classes, a classification based on the long-term hydrological record that separated sites with “normal,” saturated hydrology from those with “altered,” drier hydrology, clearly distinguished sites with different nitrogen cycling function. Based on these findings, current practices for designating reference standard sites to judge wetland functions, at least in urbanized regions, are ineffective and potentially misleading. At least one year of hydrological monitoring data is suggested to classify wetlands into groups that have different nutrient cycling functions, particularly in urban landscapes.     

Groundwaters at Risk: Wetland Loss Changes Sources,Lengthens Pathways,and Decelerates Rejuvenation of Groundwater Resources

Wetland loss alters the hydrology of wetlandscapes in poorly understood ways. To quantify the effects of wetland loss on subsurface hydrology, a physically based hydrologic model that simulates the timing and pathways of subsurface hydrologic connections was coupled with wetland inventories over a 50‐year period during which substantial wetland loss occurred. The model revealed, based on vertical variations in saturated hydraulic conductivities, wetland loss of different degrees led to a contraction of catchment contributing areas to local surface waters but an expansion of contributing areas to the regional surface water body. This shift in groundwater contributing areas reflected (1) a decrease in baseflow contribution to the local surface water bodies, and (2) an increase in the transit time and length of subsurface hydrologic connections with an associated increase in the magnitude and age of baseflow discharging to the regional surface water body. The model also showed regions with thick permeable aquifers were particularly sensitive to the loss of wetlands. Our ability to predict these changes in hydrology of the watershed provides important support for designing science‐based policies to promote sustainable water resource management.     

Testing the Basic Assumption of the Hydrogeomorphic Approach to AssessingWetland Functions

The hydrogeomorphic (HGM) approach for developing "rapid" wetland function assessment methods stipulates that the variables used are to be scaled based on data collected at sites judged to be the best at performing the wetland functions (reference standard sites). A critical step in the process is to choose the least altered wetlands in a hydrogeomorphic subclass to use as a reference standard against which other wetlands are compared. The basic assumption made in this approach is that wetlands judged to have had the least human impact have the highest level of sustainable performance for all functions. The levels at which functions are performed in these least altered wetlands are assumed to be "characteristic" for the subclass and "sustainable." Results from data collected in wetlands in the lowlands of western Washington suggest that the assumption may not be appropriate for this region. Teams developing methods for assessing wetland functions did not find that the least altered wetlands in a subclass had a range of performance levels that could be identified as "characteristic" or "sustainable." Forty-four wetlands in four hydrogeomorphic subclasses (two depressional subclasses and two riverine subclasses) were rated by teams of experts on the severity of their human alterations and on the level of performance of 15 wetland functions. An ordinal scale of 1-5 was used to quantify alterations in water regime, soils, vegetation, buffers, and contributing basin. Performance of functions was judged on an ordinal scale of 1-7. Relatively unaltered wetlands were judged to perform individual functions at levels that spanned all of the seven possible ratings in all four subclasses. The basic assumption of the HGM approach, that the least altered wetlands represent "characteristic" and "sustainable" levels of functioning that are different from those found in altered wetlands, was not confirmed. Although the intent of the HGM approach is to use level of functioning as a metric to assess the ecological integrity or "health" of the wetland ecosystem, the metric does not seem to work in western Washington for that purpose.     

Hydrogeomorphic Reference Condition and Its Relationship with Macroinvertebrate Assemblages in Southeastern U.S. Sand Hills Streams

Defining stream reference conditions is integral to providing benchmarks to ecological perturbation. We quantified channel geometry, hydrologic and environmental variables, and macroinvertebrates in 62 low‐gradient, SE United States (U.S.) Sand Hills (Level IV ecoregion) sand‐bed streams. To identify hydrogeomorphic reference condition (HGM), we clustered channel geometry deviation from expectations given watershed area (Aws), resulting in two HGM groups discriminated by area at the top of bank (Atob) residuals <0.6 m² and >0.6 m² predicted to be HGM reference/nonreference streams, respectively. Two independent partial least squares discriminate analyses used (1) hydrologic/environmental variables and (2) macroinvertebrate mean trait values (mT) on 10 reference/nonreference stream pairs of similar Aws for classification validation. Nonreference streams had flashier hydrographs and altered flow magnitudes, lower organic matter, coarser substrate, higher pH/specific conductivity compared with reference streams. Macroinvertebrate assemblages corresponded to HGM groupings, with mT indicative of multivoltinism, collector‐gatherer functional feeding groups, fast current‐preference taxa, and lower Ephemeroptera, Plecoptera, and Trichoptera richness and biotic integrity in nonreference streams. HGM classifications in Sand Hills, sand‐bed streams were determined from channel geometry. This easily implemented classification is indicative of contemporary hydrologic disturbance resulting in contrasting macroinvertebrate assemblages.     

Comparison of Function of Created Wetlands of Two Age Classes in Central Pennsylvania

Hydrogeomorphic (HGM) functional assessment models were used to assess whether function in created wetlands of two ages (1 year old and >12 years old) was equivalent to that of natural (reference) mainstem floodplain wetlands. Reference wetlands scored higher than both created age classes for providing energy dissipation and short-term surface water storage. Reference wetlands scored higher in maintaining native plant community and structure than 1-year-old sites, and 12-year-old wetlands scored higher than reference sites for providing vertebrate habitat structure. Analysis of individual model variables showed that reference wetlands had greater vegetative biomass and higher soil organic matter content than both created wetland age classes. Created wetlands were farther from natural wetlands and had smaller mean forest patch sizes within a 1-km-radius circle around the site than did the reference sites, indicating less hydrologic connectivity. Created wetlands also had less microtopographic variation than reference wetlands. The 1-year-old created sites were placed in landscape settings with greater land use diversity and road density than reference sites. The 12-year-old sites had a higher gradient and a higher percentage of their surrounding area in urban land use. These results show that the created wetlands were significantly structurally different (if not functionally so) from reference wetlands even after 12 years. The most profound differences were in hydrology and the characteristics of the surrounding landscape. More attention needs to be focused on placing created wetlands in appropriate settings to encourage proper hydrodynamics, eliminate habitat fragmentation, and minimize the effects of stressors to the site.     

Evaluating cumulative effects of disturbance on the hydrologic function of bogs,fens, and mires

Few quantitative studies have been done on the hydrology of fens, bogs, and mires. Consequently predicting the cumulative impacts of disturbances on their hydrologic functions is extremely difficult. For example, few data are available on the role of bogs and fens with respect to flood desynchronization and shoreline anchoring. However, recent studies suggest that very small amounts of groundwater discharge are sufficient to radically modify mire surface-water chemistry, and consequently, vegetation communities and their associated surface-water hydrology. Bogs and fens are, in a sense, hydrobiologic systems, and any evaluation of cumulative impacts will have to (1) consider the complicated and little understood interactions among wetland hydrology, water chemistry, and biota, and (2) place the effect of individual wetland impacts within the context of the cumulative impacts contributed to the watershed from other geomorphic areas and land uses.It is difficult to evaluate the potential cumulative impacts on wetland hydrology because geologic settings of wetlands are often complex and the methods used to measure wetland streamflow, groundwater flow, and evapotranspiration are inexact (Winter 1988). This is especially so for bogs, fens, and mires underlain by thick organic soils. These wetlands, found in the circumboreal areas of North America, Europe, and Asia, are major physiographic features in eastern North America, northern Europe, and Siberia (Kivenen and Pakarinen 1981, Gore 1983, Glaser and Janssens 1986). Their very scale makes it difficult to quantify the hydrologic function accurately. The hydrology of small bogs and fens found elsewhere is just as poorly understood because of conflicting conceptual models of pertinent hydrologic processes.This article (1) reviews our current understanding of the hydrologic function of bogs, fens, and mires at different scales and in different physiographic settings and (2) presents hypotheses on potential cumulative impacts on the hydrologic function that might occur with multiple disturbances.     

Evaluation of US EPA Environmental Monitoring and Assessment Program's (EMAP)-Wetlands sampling design and classification

The United States Environmental Protection Agency's Environmental Monitoring and Assessment Program (EMAP) will monitor the nation's resources by evaluating the status and trends of selected indicators of condition using a probability-based sampling design. The EMAP-Wetlands program will monitor the condition of the nation's wetlands. The EMAP classification system is an aggregation of the many subclasses of the US Fish and Wildlife Service's National Wetlands Inventory (NWI) classification system. This aggregation results in fewer wetland classes with more wetlands per class than the NWI system. Aggregation of the NWI classification was based primarily on dominant vegetation cover, flooding regimes, dominant water source, and adjacency to rivers and lakes. We evaluated the EMAP classification system and sampling design using NWI digital wetlands data for portions of Illinois, Washington, North Dakota, and South Dakata. Relative numbers of wetlands, total areas, average areas, and common versus rare classes were compared between the EMAP and NWI classification systems. As expected, the EMAP classification provided fewer wetland polygons, each with larger areas, without altering total wetland area. Summary statistics comparing sample estimates to true population parameters (represented by the NWI data) demonstrated the effectiveness of the EMAP sampling design with the exception of rare EMAP classes in the selected regions. Although simple random sampling is inadequate for both large and small wetlands, the EMAP sampling design is readily adapted to provide better estimates for these categories. Aggregating the NWI classification to the EMAP classification provides fewer wetland classes, with more wetlands per class, for EMAP's annual reports and statistical summaries. The research in this report has been funded by the United States Environmental Protection Agency (EPA) under contracts 68-C8-0006 to ManTech Environmental Technology, Inc. and 68-03-3532 to The Bionetics Corporation. Mention of trade names does not constitute endorsement or recommendation for use.     

DESIGNATION OF WETLANDS BY WEIGHTED AVERAGES OF VEGETATION DATA: A PRELIMINARY EVALUATION1

ABSTRACT: Weighted averages (WA) was investigated as a vegetation-based method for wetland designation, to be used in conjunction with the wetland indicator status of plants from Wetland Plants of the United States of America 1986 (Reed, 1986). Ecological indices were assigned to indicator groups and were used to compute weighted averages for quantitative data obtained from four studies of wetland vegetation conducted in various regions of the United States. Weighted averages of vegetation data proved to be a useful tool for assessing wetland status of the vegetation types included in our study: (1) rankings of vegetation stands or types by WA correlated well with their positions on environmental moisture gradients; and (2) the results of WA could be used, together with a wetland/upland break-point, to designate vegetation types as wetland or upland in a way that agreed well, in three of the four studies, with an alternative classification of wetland habitats. The variation of weighted averages among the sampling units representing a vegetation type was generally small relative to the range of ecological indices assigned. However, designations based on weighted average scores close to the break-point should be considered provisional and must be verified with supplementary data on soils and hydrology.     

GIS‐BASED HYIROLOGIC MODELING OF RIPARIAN AREAS: IMPLICATIONS FOR STREAM WATER QUALITY1

ABSTRACT: Riparian buffers have potential for reducing excess nutrient levels in surface water. Spatial variation in riparian buffer effectiveness is well recognized, yet researchers and managers still lack effective general tools for understanding the relevance of different hydrologic settings. We present several terrain‐based GIS models to predict spatial patterns of shallow, subsurface hydrologic flux and riparian hydrology. We then link predictions of riparian hydrology to patterns of nutrient export in order to demonstrate potential for augmenting the predictive power of land use/land cover (LU/LC) maps. Using predicted hydrology in addition to LUILC, we observed increases in the explained variation of nutrient exports from 290 sites across Lower Michigan. The results suggest that our hydrologic predictions relate more strongly to patterns of nutrient export than the presence or absence of wetland vegetation, and that in fact the influence of vegetative structure largely depends on its hydrologic context. Such GIS models are useful and complimentary tools for exploring the role of hydrologic routing in riparian ecosystem function and stream water quality. Modeling efforts that take a similar GIS approach to material transport might be used to further explore the causal implications of riparian buffers in heterogeneous watersheds.     

WETLAND ECOSYSTEM STUDIES FROM A HYDROLOGIC PERSPECTIVE1

ABSTRACT: Selected studies from the literature were reviewed to determine the extent of knowledge about the relationship between hydrology and wetland ecosystem studies. Wetland studies of chemical input-output relationships have been the most dependent on hydrologic data of all wetland investigations; yet, very few of these studies have attempted to measure all components of a wetland's water balance. Usually, unmeasured components were calculated as the difference between measured inputs and outputs. Ground water frequently was overlooked. Chemical input-output investigations primarily were concerned with determining the amount of input retained in the wetlands. Few studies also included direct measurement of biogeochemical processes within wetlands of elements that were part of simultaneous input-output investigations. The importance of uncertainties in chemical budgets that are due to uncertainties in hydrologic budgets has been addressed in very few wetland investigations. Although many studies have emphasized the importance of hydrology to wetland ecosystem research, few studies have documented this, so that hydrology remains one of the least understood components of wetland ecosystems.     

LOS ALAMOS FOREST FIRE IMPACT MODELING1

ABSTRACT: A controlled burn at Bandelier National Monument got out of control and burned about 43,000 acres (17,400 hectares) near Los Alamos, New Mexico, in May 2000. The wildfire caused dramatic changes in infiltration capacity and wettability of soils in many of the watersheds above Los Alamos National Laboratory (LANL) and destroyed the duff layer, dramatically reducing the interception and infiltration capacity of the formerly forested watersheds. These sudden changes in basin hydrology necessitated a rapid assessment of hydrology and hydraulics for the canyons running through LANL property to evaluate flood risk, plan emergency flood protection measures, and assess potential sediment and actinide transport. This paper presents the results of hydrologic and hydraulic modeling of Los Alamos Canyon following the wildfire. The large scale modeling effort, with over 13,000 cross sections for the hydraulic model (5,000 for Los Alamos Canyon, 8,000 for Guaje Canyon), relied heavily on a geographic information system (GIS) for model input and floodplain delineation. The HEC‐geoRAS model provided good integration between the hydraulic model (HEC‐RAS, Version 3.0.1) and the GIS (ArcView, v. 3.3). These modeling results are being used in drainage master planning efforts at LANL and in the development of sediment transport models using HEC‐6T. Sediment transport modeling results will be used to develop actinide transport models for the canyons at LANL.     

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.     

Bioaccumulation of selenium from natural geologic sources in western states and its potential consequences

Ecological impacts of water-quality problems have developed in the western United States resulting from the disposal of seleniferous agricultural wastewater in wetland areas. Overt effects of selenium toxicosis occurred at five areas where deformities of wild aquatic birds were similar to those first observed at Kesterson National Wildlife Refuge in the west-central San Joaquin Valley of California. These areas are: Tulare Lake Bed Area, California, Middle Green River Basin, Utah, Kendrick Reclamation Project Area, Wyoming, Sun River Basin, Montana, and Stillwater Wildlife Management Area, Nevada. Potential for ecological damage is indicated at six more sites in Oregon, Colorado, the Colorado/Kansas border, and South Dakota out of 16 areas in 11 states where biological tissue data were collected. This conclusion is based on the fact that selenium bioaccumulated in bird livers to median levels that had exceeded or were in the range associated with adverse reproductive effects. Selenium concentrations in samples of fish and bird eggs support these conclusions at a majority of these areas. Reason for concern is also given for the lower Colorado River Valley, although this is not exclusively a conclusion from these reconnaissance data. Biogeochemical conditions and the extent of selenium contamination of water, bottom sediment, and biota from which this assessment was made are given here. In a companion paper, the biogeochemical pathway postulated for selenium contamination to take place from natural geologic sources to aquatic wildlife is defined.     

GEOGRAPHIC VARIATION OF THE STATISTICAL PARAMETERS OF HYDROMETEOROLOGICAL DATA OF THE UNITED STATES1

ABSTRACT: Having determined various statistical parameters for five mean monthly hydrometeorological time series of the United States, directional variation of the autocorrelation numbers, their spatial distribution over the United States and statistical significances are presented in this paper. Two conceptually different approaches (one using directional strips of large combined watersheds and the other with analysis of variances) are employed to explore the geographic variation of the statistical parameters (autocorrelations and explained variances) in question. Results adequately indicate the correspondence between these variations and the observed topographic, climatologic and hydrologic characteristics over the United States.     

Assessment of cropland inundation due to the operation of the Reelfoot Lake spillway in West Tennessee

Running Reelfoot Bayou (RRB) is the outlet canal of Reelfoot Lake, the largest natural lake in Tennessee. RRB is not able to contain discharge from Reelfoot Lake greater than the bankfull discharge of 28 m³/s (1000 ft³/s), which typically occurs at the beginning of the growing season (April–June). Historically, the planting of crops has been delayed until flooding subsides and cropland has drained. The objective of this study is a preliminary quantification of cropland inundation to determine its spatial distribution in the RRB floodplain. Inundated croplands in the RRB floodplain were delineated over a range of spillway discharges from 2 to 57 m³/s (70–2000 ft³/s), using one-dimensional–two-dimensional hydrodynamic modeling and multispectral satellite images (Landsat 8 and Sentinel-2). The composite maps made by combining the simulated and image-derived flood maps were overlaid on the United States Department of Agriculture CropScape layer to determine the inundation of individual summer crops during the growing season. About 25% of the inundated croplands are flooded at discharges of RRB less than 28 m³/s, implying wetland hydrology. The results of this analysis can be used to inform operational management of the Reelfoot Lake spillway.     

HYDROGEOLOGIC EVALUATION OF WETLAND BASINS FOR LAND USE PLANNING1

ABSTRACT: Wetlands occur in a variety of geologic, hydrologic, and topographic settings and exhibit diverse hydrogeologic characteristics. A wetland is more than an organic mat - it is an element in a larger hydrogeologic system. Three distinct but related effects of wetlands are: modifying the character of runoff from a basin, influencing the discharge/recharge relationship with the underlying aquifer, and effecting the potential for ground water development in a wetland dominated basin. An important goal of wetland research is to define the diverse roles that wetlands play in the regional hydrology and to define the geologic, hydrologic, and topographic factors that will allow meaningful distinctions among wetlands. Geologic and hydrologic factors include character and thickness of surficial materials; bedrock type; hydrologic position; permeability of organic layer; transmissivity, quality, and hydrologic connection of wetland related aquifers; ground water outflow; and depth of standing water. Topographic factors are position in the drainage basin, relative size, and absolute size of wetlands. A wetland classification to aid hydrologists and land use planners is proposed using selected factors involving hydrologic position, topographic position, and geologic type.     

“The Kesterson effect”

Hypothesized to be derived from Cretaceous marine sedimentary rocks, selenium contamination of the Kesterson National Wildlife Refuge is traced through irrigation drainage to the source bedrock of the California Coast Ranges. This biogeochemical pathway of selenium is defined here as the “Kesterson effect.” At the refuge ponds, this effect culminated in 1983 in a 64% rate of deformity and death of embryos and hatchlings of wild aquatic birds. From the previous companion paper on irrigation drainage, the Kesterson effect has been implicated in nine of 11 reconnaissance areas studied in the western United States. Deformities have resulted in at least five of these sites. Climatic, geologic, hydrologic, and soil conditions in these reconnaissance areas are similar to those in the area surrounding Kesterson National Wildlife Refuge in the west-central San Joaquin Valley of California. In California, selenium, as selenate, was ultimately found weathered with sulfur from marine sources in soluble sodium and magnesium sulfate salts, which are concentrated by evaporation on farmland soils. The Se, mobilized by irrigation drainage, is bioaccumulated to toxic levels in refuge wetland ponds that are located mainly in hydrologically closed basins and thus act as concentrating disposal points. The depositional environment of the ponds may be similar to that of the nutrient-rich continental shelf edge and slope in which Cretaceous, Eocene, and Miocene sediments found to be seleniferous in the California Coast Ranges were deposited. Bioaccumulation may be therefore a primary mechanism of selenium enrichment in ancient sediments in addition to that of the formerly suggested Cretaceous volcanic pathway.