Search USGSSearch

SEARCH · Search USGS

Results for “Environmental and Sustainability Indicators”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

209 records · Page 12Linked to original sources

Ground-water recharge in the arid and semiarid southwestern United States

Ground-water recharge in the arid and semiarid southwestern United States results from the complex interplay of climate, geology, and vegetation across widely ranging spatial and temporal scales. Present-day recharge tends to be narrowly focused in time and space. Widespread water-table declines accompanied agricultural development during the twentieth century, demonstrating that sustainable ground-water supplies are not guaranteed when part of the extracted resource represents paleorecharge. Climatic controls on ground-water recharge range from seasonal cycles of summer monsoonal and winter frontal storms to multimillennial cycles of glacial and interglacial periods. Precipitation patterns reflect global-scale interactions among the oceans, atmosphere, and continents. Large-scale climatic influences associated with El Niño and Pacific Decadal Oscillations strongly, but irregularly, control weather in the study area, so that year-to-year variations in precipitation and ground-water recharge are large and difficult to predict. Proxy data indicate geologically recent periods of naturally occurring multidecadal droughts unlike any in the modern instrumental record. Any anthropogenically induced climate change will likely reduce ground-water recharge through diminished snowpack at higher elevations. Future changes in El Niño and monsoonal patterns, both crucial to precipitation in the study area, are highly uncertain in current models. Current land-use modifications influence ground-water recharge through vegetation, irrigation, and impermeable area. High mountain ranges bounding the study area—the San Bernadino Mountains and Sierra Nevada to the west, and the Wasatch and southern Colorado Rocky Mountains to the east—provide external geologic controls on ground-water recharge. Internal geologic controls stem from tectonic processes that led to numerous, variably connected alluvial-filled basins, exposure of extensive Paleozoic aquifers in mountainous recharge areas, and distinct modes of recharge in the Colorado Plateau and Basin and Range subregions. The chapters in this professional paper present (first) an overview of climatic and hydrogeologic framework (chapter A), followed by a regional analysis of ground-water recharge across the entire study area (chapter B). These are followed by an overview of site-specific case studies representing different subareas of the geographically diverse arid and semiarid southwestern United States (chapter C); the case studies themselves follow in chapters D–K. The regional analysis includes detailed hydrologic modeling within the framework of a high-resolution geographic-information system (GIS). Results from the regional analysis are used to explore both the distribution of ground-water recharge for mean climatic conditions as well as the influence of two climatic patterns—the El Niño-Southern Oscillation and Pacific Decadal Oscillation—that impart a high degree of variability to the hydrologic cycle. Individual case studies employ a variety of geophysical and geochemical techniques to investigate recharge processes and relate the processes to local geologic and climatic conditions. All of the case studies made use of naturally occurring tracers to quantify recharge. Thermal and geophysical techniques that were developed in the course of the studies are presented in appendices. The quantification of ground-water recharge in arid settings is inherently difficult due to the generally low amount of recharge, its spatially and temporally spotty nature, and the absence of techniques for directly measuring fluxes entering the saturated zone from the unsaturated zone. Deep water tables in arid alluvial basins correspond to thick unsaturated zones that produce up to millennial time lags between changes in hydrologic conditions at the land surface and subsequent changes in recharge to underlying ground water. Recent advances in physical, chemical, isotopic, and modeling techniques have fostered new types of recharge assessments. Chemical and isotopic techniques include an increasing variety of environmental tracers that are useful and robust. Physically based techniques include the use of heat as a tracer and computationally intensive geophysical imaging tools for characterizing hydrologic conditions in the unsaturated zone. Modeling-based techniques include spatially distributed water-budget computations using high-resolution remotely sensed and ground-based geographic data. Application of these techniques to arid and semiarid settings in the southwestern United States reveals distinct patterns of recharge corresponding to geologic setting, climatic and vegetative history, and land use. Analysis of recharge patterns shows that large expanses of alluvial basin floors are drying out under current climatic conditions, with little to no recharge to underlying ground water. Ground-water recharge occurs mainly beneath upland catchments in which thin soils overlie permeable bedrock, ephemeral channels in which flow may average only several hours per year, and active agricultural areas. The chapters in this professional paper represent a coordinated attempt to develop a better understanding of one of the Nation's most critical yet difficult-to-quantify renewable resources.

Professional Paper

Monitoring to assess progress toward meeting the Assabet River, Massachusetts, phosphorus total maximum daily load - Aquatic macrophyte biomass and sediment-phosphorus flux

In 2004, the Total Maximum Daily Load (TMDL) for Total Phosphorus in the Assabet River, Massachusetts, was approved by the U.S. Environmental Protection Agency. The goal of the TMDL was to decrease the concentrations of the nutrient phosphorus to mitigate some of the instream ecological effects of eutrophication on the river; these effects were, for the most part, direct consequences of the excessive growth of aquatic macrophytes. The primary instrument effecting lower concentrations of phosphorus was to be strict control of phosphorus releases from four major wastewatertreatment plants in Westborough, Marlborough, Hudson, and Maynard, Massachusetts. The improvements to be achieved from implementing this control were lower concentrations of total and dissolved phosphorus in the river, a 50-percent reduction in aquatic-plant biomass, a 30-percent reduction in episodes of dissolved oxygen supersaturation, no low-flow dissolved oxygen concentrations less than 5.0 milligrams per liter, and a 90-percent reduction in sediment releases of phosphorus to the overlying water. In 2007, the U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, initiated studies to evaluate conditions in the Assabet River prior to the upgrading of wastewater-treatment plants to remove more phosphorus from their effluents. The studies, completed in 2008, implemented a visual monitoring plan to evaluate the extent and biomass of the floating macrophyte Lemna minor (commonly known as lesser duckweed) in five impoundments and evaluated the potential for phosphorus flux from sediments in impounded and free-flowing reaches of the river. Hydrologically, the two study years 2007 and 2008 were quite different. In 2007, summer streamflows, although low, were higher than average, and in 2008, the flows were generally higher than in 2007. Visually, the effects of these streamflow differences on the distribution of Lemna were obvious. In 2007, large amounts of floating macrophytes accumulated behind bridge constrictions and dams; in 2008, high flows during the early part of the growing season carried floating macrophytes past bridges and over dams, minimizing accumulations. Samples of Lemna were collected and weighed to provide an estimate of Lemna biomass based on areal coverage during the summer growing seasons at eight sites in the five impoundments. Average estimated biomass during 2007 was approximately twice the 2008 biomass in each of the areas monitored. In 2007, in situ hyperspectral and high-resolution, multispectral data from the IKONOS satellite were obtained to evaluate the feasibility of using remote sensing to monitor the extent of aquatic plant growth in Assabet River impoundments. Three vegetation indices based on light reflectance were used to develop metrics with which the hyperspectral and satellite data were compared. The results of the comparisons confirmed that the high-resolution satellite imagery could differentiate among the common aquatic-plant associations found in the impoundments. The use of satellite imagery could counterbalance emphasis on the subjective judgment of a human observer, and airborne hyperspectral data can provide higher resolution imagery than multispectral satellite data. In 2007 and 2008, the potential for sediment flux of phosphorus was examined in free-flowing reaches of the river and in the two largest impoundments-Hudson and Ben Smith. These studies were undertaken to determine in situ flux rates prior to the implementation of the Assabet River Total Maximum Daily Load (TMDL) for phosphorus and to compare these rates with those used in the development and evaluation of the TMDL. Water samples collected from a chamber placed on the river bottom were analyzed for total phosphorus and orthophosphorus. Ambient dissolved oxygen concentrations and seasonal temperature differences appeared to affect the rates of sequestration and sediment release of phosphorus. When dissolved oxygen concentrations remained relatively high in the chambers and when the temperature was relatively low, the tendency was for phosphorus concentrations to decrease in the chambers, indicating sediment sequestration of phosphorus; when dissolved oxygen concentrations dropped to near zero and temperatures were warmest, phosphorus concentrations increased in the chambers, indicating phosphorus flux from the sediment. The rates of release and sequestration in the in situ studies were generally comparable with the rates determined in laboratory studies of Assabet River sediment cores for State and Federal agencies. Sediment-core and chamber studies produced substantial sediment fluxes to the water column only under extremely low-DO or anaerobic conditions rarely found in the Assabet River impoundments; thus, sediment is not likely to be a major phosphorus source, especially when compared to the wastewater effluent, which sustains higher ambient concentrations. The regulatory agencies now (2011) have substantial laboratory and field data with which to determine the required 90-percent reduction in phosphorus flux after the completion of upgrades to the wastewater-treatment plants that discharge to the Assabet River.

Massachusetts

Revised groundwater-flow model of the glacial aquifer system north of Aberdeen, South Dakota, through water year 2015

The city of Aberdeen, in northeastern South Dakota, requires an expanded and sustainable supply of water to meet current and future demands. Conceptual and numerical models of the glacial aquifer system in the area north of Aberdeen were developed by the U.S. Geological Survey in cooperation with the City of Aberdeen in 2012. The U.S. Geological Survey, in cooperation with the City of Aberdeen, completed a study to revise the original numerical groundwater-flow model using data through water year (WY) 2015 to aid the City of Aberdeen in their development of plans and strategies for a sustainable water supply and to increase understanding of the glacial aquifer system and groundwater-flow system near Aberdeen. The original model was revised to improve the fit between model-simulated values and observed (measured or estimated) data, provide greater insight into surface-water interactions, and improve the usefulness of the model for water-supply planning. The revised groundwater-flow model (hereafter referred to as the “revised model”) presented in this report supersedes the original model. The purpose of this report is to describe a revised groundwater-flow model including data collection, model calibration, and model results for the glacial aquifer system including the Elm, Middle James, and Deep James aquifers north of Aberdeen, South Dakota, using updated hydrologic data through WY 2015. The original numerical model was revised in several ways. The model was modified by adding four new layers, which included a surficial layer, two intervening confining layers, and a shale bedrock layer. The revised model provides an improved understanding of the groundwater-flow system in comparison to the original model. The principal aquifers of the model area include portions of the Elm, Middle James, and Deep James aquifers. The lithologic information used to define and describe the aquifers in the model area was unaltered; however, aquifer properties and boundary conditions were reviewed and updated using geological information reported by the South Dakota Department of Environmental and Natural Resources and information obtained from geophysical investigations for this study. The horizontal extent of the Elm, Middle James, and Deep James aquifers was unaltered from the original model. The thickness of the Deep James aquifer was modified based on interpretations from the geophysical investigations. In general, groundwater in the Elm aquifer flowed from northwest to southeast and locally towards rivers and streams. Similarly, in the Middle James and Deep James aquifers, groundwater also typically flowed southeast. The revisions made to the original model include use of the following MODFLOW stress packages: Recharge, Evapotranspiration, Time-Variant Specified Head, Wells, Drains, and Stream Flow Routing, all of which were updated from the original model except for the Stream Flow Routing Package, which replaced the River Package used in the original model. Model calibration is the process of estimating model parameters to minimize the differences, or residuals, between observed data and simulated values; therefore, Parameter ESTimation (PEST) software was used to optimize model input parameters by matching model-simulated values to observed data. Calibration parameters included horizontal hydraulic conductivity, vertical hydraulic conductivity, specific yield, specific storage, and vertical streambed conductance for stream and drain cells. Multipliers were used to calibrate the recharge and evapotranspiration stresses. Evapotranspiration extinction depth also was adjusted during model calibration. Comparisons to the original model are described to highlight the changes made in the revised model. In general, the revised model adequately simulates the natural system and compares favorably with observed hydrologic data. Simulated water levels were evaluated by comparing them to single water-level observations at selected well locations. The selected wells were the same wells used in the original model. The coefficient of determination value between simulated and observed water levels for the revised model was 0.89 and included simulated and observed values from October 1, 1974 (WY 1975), through September 30, 2015 (WY 2015). The coefficient of determination value for the original model was 0.94 and included simulated and observed values from October 1, 1974, through September 30, 2009. The difference may indicate that the original model could have been overfit to hydraulic head observations because base flow was not simulated. The additional data used in the revised model included some climatically wetter, more extreme periods, such as 2011, in which annual precipitation was 30.9 inches. Average annual precipitation for the original model timeframe, which included data from WYs 1975–2009, was 20.26 inches. Additional precipitation data for WYs 2010–15, included in the revised model timeframe, resulted in an average annual precipitation for WYs 1975–2015 in the model area of 20.6 inches. The larger variability in climate data coupled with the additional water-level data could explain the lower coefficient of determination for water levels in the revised model. The revised model was used to calculate various groundwater-budget components for steady-state and transient conditions for WYs 1975–2015. The time-variant specified-head cells in the revised model had the largest change when compared to the original steady-state model for inflows and outflows. Comparing the transient budget components between the original and the revised models indicated that inflow from recharge and time-variant specified-head cells had the greatest effect on groundwater inflows, and outflow from storage had the greatest effect on groundwater outflows. The simulated potentiometric contours from the revised model were compared with (1) the observed (interpreted) potentiometric surface (layer 2) and the hydraulic head values (layers 4 and 6) and (2) the simulated contours from the original model. The simulated hydraulic gradients and general direction of groundwater flow in the Elm aquifer in the revised model generally matched the observed potentiometric contours, the simulated potentiometric contours from the original model, and general flow directions interpreted to be perpendicular to the contours. Minor discrepancies between simulated potentiometric contours from the revised model and the observed potentiometric contours may be due to the lack of observed data in the model area. The revised model was designed to reduce the limitations of the original model. The revisions were validated by comparing the results of the original model with the revised model. A primary benefit of the revised model is the inclusion of the surficial deposits and the confining units as explicit layers in the model. The addition of the surficial layer was beneficial for three primary reasons: (1) more accurate representation of recharge from precipitation, (2) more accurate representation of groundwater evapotranspiration, and (3) more accurate representation of groundwater and surface-water interactions. The groundwater model is a numeric approximation of a complex physical hydrologic system, and the revised model data were interpolated in regions with sparse data. Additionally, model discretization included averaged and interpolated values for water use, withdrawal rates, and hydraulic conductivity. The revised model provides a useful estimate for hydraulic gradients, groundwater-flow directions, and aquifer response to groundwater withdrawals.

South Dakota

Groundwater discharge characteristics for selected streams within the Loup River Basin, Nebraska, 2014–16

Streams in the Loup River Basin are sensitive to groundwater withdrawals because of the close hydrologic connection between groundwater and surface water. Groundwater discharge is the primary component of streamflow in the Loup River Basin and constitutes more than 90 percent of streamflow in the central part of the Sand Hills. To improve the understanding of geologic controls and various climatic and land-use changes on groundwater discharge, the U.S. Geological Survey (USGS), in cooperation with the Upper Loup Natural Resources District (NRD), the Lower Loup NRD, and the Nebraska Environmental Trust, studied the spatial and temporal characteristics of groundwater discharge within the Loup River Basin. This report documents the methods of data collection and analysis, which include the collection of approximately 350 river miles of aerial thermal infrared imagery and continuous groundwater-level and temperature data from six streamflow-gaging stations within the Loup River Basin. The results from the stream reconnaissance and examination of aerial thermal infrared imagery demonstrated the influence of the surficial and subsurface geology on the spatial characteristics of groundwater discharge to streams in the Loup River Basin. At the headwaters of the South Loup River, streamflow is sustained and increased from focused groundwater discharge emanating from Quaternary deposits at many small (less than 0.1 cubic foot per second) focused points. The volume of water produced from this dense network of focused groundwater discharge points along the North Fork South Loup River is sufficient to provide approximately 40 percent of the flow measured at the South Loup River at Arnold, Nebraska streamflow-gaging station (USGS station 06781600) during the irrigation season. Approximately 5 miles downstream from the South Loup River at Arnold, Nebr., streamflow-gaging station, the river incises into Pliocene-age sand and gravel deposits, which provide additional groundwater discharge to the stream. The streamflow of the South Loup River increases by a factor of 5 across a 62-mile reach of the middle South Loup River. Increases in streamflow along the upper Dismal River result from a dense network of focused groundwater discharge points within semiconsolidated Pliocene-age deposits. Below the Dismal River near Thedford, Nebr., streamflow-gaging station (USGS station 06775900), the Dismal River incises into the Ogallala Formation over a short reach before flowing over coarser, more permeable Quaternary-age alluvial deposits. Diffuse groundwater discharge sustains and increases the streamflow of the lower Dismal River in this reach. Groundwater sapping was evident on some stream reaches and has increased the size and flow of focused groundwater discharge points. Previous researchers have documented streambed incision and groundwater sapping on the upper Dismal River that have created and enlarged focused groundwater discharge points capturing additional groundwater. Similar processes appear to have played a role in the formation of larger focused groundwater discharge points, which sustain the flow of the middle South Loup River. The constant flow of groundwater into the South Loup River has removed finer-grained Quaternary sediments and further exposed Pliocene-age gravel deposits. Headward erosion is evident where some of the large focused groundwater discharge points have incised their own draws and terminate in bowl-like depressions away from the stream. Within the Loup River NRDs, the percentage of groundwater-irrigated land in a stream basin is one factor that affects groundwater discharge to streams. A striking example was at the South Loup River at Saint Michael, Nebr., groundwater and streamflow-gaging station (USGS station 06784000) where the shallow groundwater levels declined below the level of the stream during the middle to late part of the growing season (July to September) when consumptive groundwater use was at its peak. The South Loup River Basin above the South Loup River at Saint Michael, Nebr., streamflow-gaging station has the highest percentage of groundwater-irrigated row crops of all the basins examined in this study. Continuous groundwater and surface-water levels measured at the North Loup River at the Taylor, Nebr., streamflow-gaging station (USGS station 06786000) indicate that the stream is receiving groundwater throughout the year; however, when consumptive groundwater use peaks during the middle to late part of the growing season (July to September), the difference in elevation between the groundwater level and the stream elevation decreases, which indicates a reduction in the amount of groundwater discharge received.

Nebraska

In-situ arsenic removal during groundwater recharge through unsaturated alluvium

OBJECTIVES The purpose of this study was to determine the feasibility and sustainability of in-situ removal of arsenic from water infiltrated through unsaturated alluvium. BACKGROUND Arsenic is naturally present in aquifers throughout the southwestern United States and elsewhere. In January 2006, the U.S. Environmental Protection Agency (EPA) lowered the Maximum Contaminant Level (MCL) for arsenic from 50 to 10 micrograms per liter (g/L). This raised concerns about naturally-occurring arsenic in groundwater. Although commercially available systems using sorbent iron or aluminum oxide resins are available to treat high-arsenic water, these systems are expensive to build and operate, and may generate hazardous waste. Iron and aluminum oxides occur naturally on the surfaces of mineral grains that compose alluvial aquifers. In areas where alluvial deposits are unsaturated, these oxides may sorb arsenic in the same manner as commercial resins, potentially providing an effective low-cost alternative to commercially engineered treatment systems. APPROACH The Antelope Valley within the Mojave Desert of southern California contains a shallow water-table aquifer with arsenic concentrations of 5 g/L, and a deeper aquifer with arsenic concentrations of 30 g/L. Water was pumped from the deep aquifer into a pond and infiltrated through an 80 m-thick unsaturated zone as part of field-scale and laboratory experiments to treat high-arsenic groundwater and recharge the shallow water table aquifer at the site. The field-scale recharge experiment included the following steps: 1) construction of a recharge pond 2) test drilling for sample collection and instrument installation adjacent to the pond 3) monitoring downward migration of water infiltrated from the pond 4) monitoring changes in selected trace-element concentrations as water infiltrated through the unsaturated zone Data from instruments within the borehole adjacent to the pond were supplemented with borehole and surface geophysical data to evaluate the lateral spreading of water as it moved downward through the unsaturated zone. Three laboratory studies were undertaken. Sequential extraction was used to evaluate the abundance of iron, aluminum, and manganese oxides and selected trace elements on operationally defined sites on the surfaces of mineral grains collected before and after infiltration from the pond. Secondly, radio-labeled arsenic-73 microcosm experiments evaluated the potential for incorporation of arsenic sorbed to exchange sites on mineral grains into less reactive crystalline mineral structures with time. Finally, column studies evaluated arsenic sorption and the pH dependence of sorption for selected unsaturated zone materials. RESULTS/CONCLUSIONS Between December 2010 and July 2012, more than 120,000 cubic meters (m3 ) (about 97 acre-feet) of high-arsenic groundwater was pumped from the deep aquifer into a 0.11 hectare (about 0.27 acres) pond and infiltrated though an 80-meter (about 260 feet) thick unsaturated zone to recharge a water-table aquifer. Arsenic concentrations were lowered from 30 to 2 g/L as water infiltrated though the unsaturated zone at the site. Some uranium, possibly associated with past agricultural land use at the site, was mobilized to concentrations as high as 66 g/L within the unsaturated zone during the experiment. Uranium was resorbed and the high uranium concentrations did not reach the water table at the site. Concentrations of other trace elements, including antimony, chromium, vanadium, and selenium were low throughout the study. Infiltration rates from the pond were as high as 0.4 meters per day (1.1 feet per day, ft/d), and the wetting front moved downward about 25 centimeters per day (cm/d) (0.8 ft/d) to a depth of about 50 m (about 165 feet). Clay layers at that depth slowed the downward movement of the wetting front to about 5 cm/d (0.16 ft/d). Lateral movement of the wetting front was monitored using sequential direct-current (DC) surface and sequential electromagnetic (EM) and DC borehole resistivity. Most lateral movement occurred on a clay layer about 50 m (about 165 feet) below land surface. Infiltrated water reached the water table in January 2013. At the water table, the “wetted footprint” of water infiltrated from the pond, indicated by surface resistivity data, was about 13 hectares (about 32 acres). On the basis of data collected at the site, there is enough sorbent material to operate this pond and treat groundwater having an arsenic concentration of 30 g/L to 2 g/L for about 500 years. Toxicity Characteristic Leaching Procedure (TCLP) data showed arsenic concentrations to be below hazardous levels beneath the pond after the experiment. Pond maintenance may be required to keep infiltration rates high, and prevent accumulation of organic material on the pond bottom, although organic material on the pond bottom may increase removal of other trace elements in infiltrated water including chromium, selenium, and vanadium. Laboratory results are consistent with the field data and show sorption of arsenic in 10 cm (0.3 feet) columns to about 2 g/L over a pH range of 6 to 8, and at influent arsenic concentrations as high as 300 g/L, without breakthrough in 50 pore volumes. Column results suggest that the insitu treatment may remove arsenic in a range of hydrogeologic settings, and would not necessarily be restricted to alkaline alluvial aquifers common throughout the southwestern United States. Radiolabeled arsenic-73 experiments show that although arsenic is initially weakly sorbed (and potentially mobile), with time arsenic is incorporated into amorphous materials. One year after sorption onto surface exchange sites, most sorbed arsenic is incorporated into crystalline oxide minerals on the surfaces of primary mineral grains and is less mobile. Results of the study suggest that long-term land use restrictions on sites used for in-situ treatment of arsenic may not be needed to control water applied to surface materials. This minimizes some regulatory concerns about future land use at sites used for in-situ arsenic treatment. However, future land uses that may alter reduction-oxidation conditions in the subsurface should be avoided, such as infiltration of stormwater recharge or recharge with other water having high organic carbon concentrations (including unsewered residential land use, dairy or other confined animal operations).

Report

Selenium hazards in the Salton Sea environment—Summary of current knowledge to inform future wetland management

Quaternary marine and continental shales in the western United States are sources of selenium that can be loaded into the aquatic environment through mining, agricultural, and energy production processes. The mobilization of selenium from shales through agricultural irrigation has been recognized since the 1930s; however, discovery of deformities in birds and other wildlife using agricultural habitats during the 1980s spurred studies to determine the extent and effects of the contamination. Through these early studies, researchers determined that biota in the Salton Sea drainage basin was at risk from legacy selenium contamination in the Colorado River watershed. The Salton Sea and its surrounding managed and unmanaged wetlands provide vital inland habitat and trophic support for diverse assemblages of resident and migratory wildlife, and understanding regional selenium hazards for these trust species is a priority for many Federal and State agencies. The modern Salton Sea is a shallow, landlocked saline lake in Riverside and Imperial Counties (not shown) of California that is sustained by irrigation return and perennial river inflow. Changes in water transfer agreements under the 2003 Quantification Settlement Agreement (QSA) have resulted in reduced irrigation flow, declining lake levels, and the evolution of unmanaged wetlands in areas where drains and rivers no longer reach the Salton Sea. These wetlands provide additional habitat for some species of concern, but their potential to increase selenium hazards for trust species is largely unknown. From the 1980s to 2020, efforts to document selenium contamination and effects throughout the region have resulted in a considerable amount of selenium data from the Salton Sea and its surrounding drainage basin; however, no long-term (greater than 20 years), consistent sampling program has been established, and all data have been collected by different entities using a variety of protocols and analytical techniques. This lack of coordination has been previously documented in regional management plans and has led to difficulty in reliably assessing selenium hazards in the Salton Sea environment. This report provides a summary of the available disparate selenium information collected from water, sediment, and biota in the Salton Sea region since the 1980s and to identify data gaps that need to be filled to understand the potential effects of selenium on species of concern, including federally endangered desert pupfish ( Cyprinodon macularius ) and Yuma Ridgway’s Rail ( Rallus obsoletus yumanensis; formerly Yuma Clapper Rail, Rallus longirostris yumanensis ). Available data from the Salton Sea drainage basin show that water from the Colorado River has the lowest selenium concentration of all surface water sources. All other surface water flowing into the Salton Sea has elevated selenium concentrations due to evaporation and evapotranspiration that occurs in agricultural fields and associated water delivery infrastructure or leaching of selenium from irrigated farmland soils. The Salton Sea has lower selenium concentrations because of various biogeochemical processes that recycle selenium into the sediment or volatilize it to the atmosphere; however, these mechanisms are not well defined, and it is not clear if selenium cycling will change in response to possible changes in the oxidation state of the Salton Sea bottom waters as water levels decline. Agricultural drains have the highest average selenium concentrations, but few drains have been sampled since changes in irrigation practices have occurred (due to the 2003 QSA). Groundwater selenium concentrations are variable; some wells south of the Salton Sea have selenium concentrations as high as 300 micrograms per liter (µg/L), whereas selenium concentrations are below detection in other wells. Groundwater and surface-water geothermal discharge zones around the margins of the Salton Sea and in unmanaged wetlands have not been studied in detail, and published selenium measurements are not available for these surface features. Selenium concentrations in the sediment of the Salton Sea drainage basin are highest in wetland particulate organic matter and the Salton Sea lakebed, indicating that removal of selenium from the water to the sediment has been a primary mechanism for keeping selenium concentrations low in the water column. Sediment selenium concentrations in wetlands are lower than in the Salton Sea but higher than inflowing drains and rivers, indicating the lentic wetland sites also may be important sinks for selenium because of biogeochemical processes. Sediment selenium data have not been collected in agricultural drains since changes in irrigation practices occurred (due to the 2003 QSA), and it is unknown if selenium sequestration from the water column has changed in these systems. We divided biological data into broad taxonomic categories, including primary producers, invertebrates, herpetofauna, mammals, fishes, and birds to facilitate evaluation of selenium concentrations and spatiotemporal trends observed in the Salton Sea. Overall, selenium concentrations were substantially greater in algae samples compared to all vascular plant samples combined. Median selenium concentrations in several invertebrate taxa (Chironomidae, Formicidae, Corixidae, Corbiculidae and Nereididae, and Decapoda) exceeded the maximum suggested dietary threshold of 3.0–4.0 micrograms per gram (µg/g) dry weight (dw) for predators consuming invertebrates in aquatic food webs. The greatest number of samples were collected from fish, and selenium distributions among species and locations showed that the range for most samples was lower than the U.S. Environmental Protection Agency selenium criterion for aquatic life (8.5 µg/g dw whole body, 11.3 µg/g dw fillets). The median selenium concentrations for whole body fish were below the selenium criterion in most locations, except for bairdiella ( Bairdiella icistia ) from the Salton Sea and irrigation drains, a few individual tilapia spp. (family Cichlidae, including genera Tilapia, Oreochromis , and their hybrids) from the river and river outlets, and several western mosquitofish ( Gambusia affinis ) and sailfin molly ( Poecilia latipinna ) from irrigation drain outlets. For avian samples combined among years and locations, median selenium concentrations in livers from all families except waders and Ibis (family Threskiornithidae) were higher than levels expected to cause selenium toxicosis (10–20 µg/g dw), and all median egg concentrations were above or near 6.0 μg/g dw, which is a conservative threshold value for reproductive impairment. Most knowledge gaps we identified for water, sediment, and biota were interrelated, and the use of integrated approaches to address knowledge gaps can provide greater insight into the drivers behind selenium hazards. Integrated water, sediment, and biota studies could help identify cost-effective management solutions that serve multiple purposes. A comprehensive analysis of the hydrology, biogeochemistry, and food-web processes in wetlands and other habitats can inform predictive models to identify drivers of selenium bioavailability, uptake from the environment and subsequent trophic transfer, ultimately forming the basis for experimental habitat management manipulations to minimize selenium hazards to wildlife. Furthermore, a comprehensive, long-term sampling and analytical laboratory plan would enable comparison of data among different entities that are sampling at the Salton Sea. Such efforts are well suited to help fill knowledge gaps that preclude understanding of selenium hazards and future management options for biota using Salton Sea habitats, including newly formed wetlands throughout the region. All data compiled for this report are available in two U.S. Geological Survey data releases: Groover and others (2022) for water and sediment samples and De La Cruz and others (2022) for biological samples. The data releases include all publicly available data for selenium concentrations in water, sediment, and biological samples collected in and around the Salton Sea, including the Coachella and Imperial Valleys. The data releases also include previously unpublished data.

California

Grassland birds wintering at U.S. Navy facilities in southern Texas

Grassland birds have undergone widespread decline throughout North America during the past several decades. Causes of this decline include habitat loss and fragmentation because of conversion of grasslands to cropland, afforestation in the East, brush and shrub invasion in the Southwest and western United States, and planting of exotic grass species to enhance forage production. A large number of exotic plant species, including grasses, have been introduced in North America, but most research on the effects of these invasions on birds has been limited to breeding birds, primarily those in northern latitudes. Research on the effects of exotic grasses on birds in winter has been extremely limited. This is the first study in southern Texas to examine and compare winter bird responses to native and exotic grasslands. This study was conducted during a period of six years (2003–2009) on United States Navy facilities in southern Texas including Naval Air Station–Corpus Christi, Naval Air Station–Kingsville, Naval Auxiliary Landing Field Waldron, Naval Auxiliary Landing Field Orange Grove, and Escondido Ranch, all of which contained examples of native grasslands, exotic grasslands, or both. Data from native and exotic grasslands were collected and compared for bird abundance and diversity; ground cover, vegetation density, and floristic diversity; bird and vegetation relationships; diversity of insects and arachnids; and seed abundance and diversity. Effects of management treatments in exotic grasslands were evaluated by comparing numbers and diversity of birds and small mammals in mowed, burned, and control areas. To determine bird abundance and bird species richness, birds were surveyed monthly (December–February) during the winters of 2003–2008 in transects (100 meter × 20 meter) located in native and exotic grasslands distributed at all five U.S. Navy facilities. To compare vegetation in native and exotic grasslands, vegetation characteristics were measured during 2003–2008 in the same transects used for bird surveys and included five measures of ground cover, plus estimates of plant species richness, vegetation density (visual obstruction) at two different heights, and shrub numbers. These data, plus seasonal rainfall, were then used to evaluate components of variation in native and exotic grasslands. Relations between total bird numbers and bird species richness with environmental variation in native and exotic grasslands were compared. To compare diversity of arthropods in native and exotic grasslands, insects and arachnids were collected using three different methodologies (standardized sweep-net, random sweep-net, and pitfall traps) during four seasons, (2005–2006), at Naval Air Station–Corpus Christi, Naval Auxiliary Landing Field Waldron, and Naval Air Station–Kingsville. To compare seed abundance and diversity between native and exotic grasslands, seeds were collected for two winters (2004–2006) at Naval Air Station–Corpus Christi and Naval Air Station–Kingsville. To evaluate effects of management on grassland vertebrates, abundance and diversity of birds and small mammals were estimated and compared in exotic grasses subjected to mowing, burning, or no active management (control) for one full year (2008–2009). Observations were made of 1,044 birds of 30 species in grassland transects during five winters. The Savannah Sparrow ( Passerculus sandwichensis ) was the most common bird, which, with 644 detections, accounted for 63 percent of all individuals identified to species. Meadowlarks ( Sturnella spp. ) and Le Conte’s Sparrows ( Ammodramus leconteii ) were the second (10 percent) and third (7 percent) most abundant bird species, respectively. Six of the seven most abundant species detected in grasslands were grassland species, and their numbers accounted for 87 percent of all birds, but 20 of the 30 species (67 percent) that used grasslands were not grassland species. Seven species observed in grassland transects during the study were Species of Conservation Concern: Le Conte’s Sparrow, Sedge Wren ( Cistothorus platensis ), Grasshopper Sparrow ( Ammodramus savannarum ), Long-billed Curlew ( Numenius americanus ), Sprague’s Pipit ( Anthus spragueii ), Cassin’s Sparrow ( Aimophila cassinii ), and Loggerhead Shrike ( Lanius ludovicianus ). Native grasslands consistently supported greater bird species richness than exotic grasslands. In one winter, exotic grasslands supported more birds than native grasslands. Native grasslands were determined to have more forb cover, more bare ground, and greater plant species richness than exotic grasslands, whereas exotic grasslands were characterized by more grass cover and relatively greater vegetation density during dry years. Not only did these individual measures differ between native and exotic grasslands, but components of variation also differed. In native grasslands, grass density and cover contributed more to variation, whereas in exotic grasslands, non-grass vegetation was a greater component of variation. Total bird numbers and bird species richness in native grasslands were related to the principal component that contained a measure of litter cover. Total bird numbers and bird species richness in exotic grasslands indicated no significant relationships with any of the principal components of variation. The two most common insect orders in native grasslands were Hymenoptera and Coleoptera, which accounted for 42 percent of all insects. The two most common insect orders in exotic grasslands were Hemiptera and Homoptera, which accounted for about 80 percent of all insects. Insect family richness was greater in exotic grasslands than in native grasslands in two of four seasons. Proportions of arachnid families were similar in native and exotic grasslands, but arachnid family richness was greater in exotic grasslands than in native grasslands. Abundance of seeds was greater in exotic than in native grasslands. However, seed diversity was greater in native grasslands than in exotic grasslands. Among the three types of management (mowed, burned, and control) applied to exotic grasses, birds were most abundant in the mowed area. Sedge Wrens, however, were never encountered in mowed sites. Meadowlarks were similarly abundant in all treatments, but Le Conte’s Sparrows were detected only in the control (unmanaged) area. Hispid cotton rats ( Sigmodon hispidus ) accounted for 93 percent of all rodent captures, with the number of captures peaking December through February. Hispid cotton rat numbers and total rodent numbers were greatest in control and pre-burn areas, and lowest in the mowed area. Mammal diversity, however, was greatest in the mowed habitat. Native and exotic grasslands differed essentially in all categories (bird numbers and diversity, vegetation characteristics, components of variation, diversity of insects and arachnids, and seed abundance and diversity) used to measure and compare them. This indicates that fundamental ecosystem processes have been altered after native grasslands have undergone invasion and ultimate domination by exotic grass species. Future research in Texas grassland ecosystems is essential because: 1) Texas sustains more area in grasslands than any other state or province in the Central Flyway; 2) Texas serves as the winter destination or migration pathway for hundreds of species of birds, including winter residents and Neotropical migrants; 3) ecology, distribution, and numbers of grassland birds wintering in southern latitudes of the United States remains poorly understood; and 4) climate change threatens to further accelerate advances of invading grass species.

Open-File Report

Hydrogeology and hydrology of the Punta Cabullones wetland area, Ponce, southern Puerto Rico, 2007-08

The U.S. Geological Survey, in cooperation with the Municipio Autónomo de Ponce and the Puerto Rico Department of Natural and Environmental Resources, conducted a study of the hydrogeology and hydrology of the Punta Cabullones area in Ponce, southern Puerto Rico. (Punta Cabullones is also referred to as Punta Cabullón.) The Punta Cabullones area is about 9 square miles and is an ecological system made up of a wetland, tidal flats, saltflats, mangrove forests, and a small fringing reef located a short distance offshore. The swales or depressions between successive beach ridges became development avenues for saline to hypersaline wetlands. The Punta Cabullones area was designated by the U.S. Fish and Wildlife Service as a coastal barrier in the 1980s because of its capacity to act as a buffer zone to ameliorate the impacts of natural phenomenon such as storm surges. Since 2003, Punta Cabullones has been set aside for preservation as part of the mitigation effort mandated by Federal and State laws to compensate for the potential environmental effects that might be caused by the construction of the Las Américas Transshipment Port. Total rainfall measured during 2008 within the Punta Cabullones area was 36 inches, which is slightly greater than the long-term annual average of 32 inches for the coastal plain near Ponce. Two evapotranspiration estimates, 29 and 37 inches, were obtained for the subarea of the Punta Cabullones area that is underlain by fan-delta and alluvial deposits by using two variants of the Penman semi-empirical equation. The long-term water stage and chemical character of the wetland in Punta Cabullones are highly dependent on the seasonal and annual variations of both rainfall and sea-wave activity. Also, unseasonal short-term above-normal rainfall and sea-wave events resulting from passing storms may induce substantial changes in the water stage and the chemical character of the wetland. In general, tidal fluctuations exert a minor role in modifying the water quality and stage of the wetland in Punta Cabullones. The role of the tidal fluctuations becomes important during those times when the outlets/inlets to the sea are not blocked by a sand bar and is allowed to freely flow into the wetland interior. The salinity of the wetland varies from brackish to hypersaline. The hypersaline conditions, including the occurrence of saltflats, within the Punta Cabullones wetland area result from a high evapotranspiration rate. The hypersaline conditions are further enhanced by a sand bar that blocks the inlet/outlet of the wetland’s easternmost channel, particularly during the dry season. Groundwater in Punta Cabullones mostly is present within beds of silisiclastic sand and gravel. During the study period, the depth to groundwater did not exceed 4 feet below land surface. The movement and direction of the groundwater flow in Punta Cabullones are driven by density variations that in turn result from the wide range of salinities in the groundwater. The salinity of the groundwater decreases within the first 60 to 100 feet of depth and decreases outward from a mound of hypersaline groundwater centered on piezometer nest PN2. The main groundwater types within the Punta Cabullones area vary from calcium-bicarbonate type in the northernmost part of the study area to a predominantly sodium-potassium-chloride groundwater type southward. According to stable-isotope data, groundwater within the study area is both modern meteoric water and seawater highly affected by evaporation. The chemical and stable-isotopic character of local groundwater is highly influenced by evapotranspiration because of its shallow depth. Equivalent freshwater heads indicate groundwater moves away from a mound centered on piezometer nest PN2, in a pattern similar to the spatial distribution of groundwater salinity. Vertical groundwater flow occurs in Punta Cabullones due to local differences in density. In the wetland subarea of Punta Cabullones, groundwater and surface water are hydraulically coupled. Locally, surface-hypersaline water sinks into the aquifer, providing recharge and serving as a mechanism to redistribute salt throughout the study area. The evapotranspiration in the wetland subarea is estimated at about 11 million gallons per day (Mgal/d) that is equivalent to about 12,586 acre-feet per year. The balance of evapotranspiration, in excess of the about 0.5 Mgal/d of groundwater flow within the wetland, is supplied by saline to hypersaline surface water that may include seawater and meteoric water highly affected by evaporation with dissolved salts. In one of the extreme scenarios in which no groundwater is intercepted by pumpage at the Restaurada well field, the amount of saline to hypersaline water in the wetland consumed by evapotranspiration is about 10.5 Mgal/d. In the opposite extreme in which the entire regional groundwater flow is intercepted by pumpage in the Restaurada well field, the entire evapotranpiration requirement is met by saline to hypersaline water. Hydrologic, isotopic, and chemical data indicate that all of, or a large portion of, the historical groundwater flow to Punta Cabullones is being captured by the Puerto Rico Aqueducts and Sewer Authority pumpage at the Restaurada well field at a rate of about 2 Mgal/d. As a consequence, seawater intrusion into the aquifer at the Punta Cabullones area seems to be occurring, while the current pumpage at the Restaurada well field is sustained by storage depletion of the aquifer.

Puerto Rico

Sequence stratigraphy, seismic stratigraphy, and seismic structures of the lower intermediate confining unit and most of the Floridan aquifer system, Broward County, Florida

Deep well injection and disposal of treated wastewater into the highly transmissive saline Boulder Zone in the lower part of the Floridan aquifer system began in 1971. The zone of injection is a highly transmissive hydrogeologic unit, the Boulder Zone, in the lower part of the Floridan aquifer system. Since the 1990s, however, treated wastewater injection into the Boulder Zone in southeastern Florida has been detected at three treated wastewater injection utilities in the brackish upper part of the Floridan aquifer system designated for potential use as drinking water. At a time when usage of the Boulder Zone for treated wastewater disposal is increasing and the utilization of the upper part of the Floridan aquifer system for drinking water is intensifying, there is an urgency to understand the nature of cross-formational fluid flow and identify possible fluid pathways from the lower to upper zones of the Floridan aquifer system. To better understand the hydrogeologic controls on groundwater movement through the Floridan aquifer system in southeastern Florida, the U.S. Geological Survey and the Broward County Environmental Planning and Community Resilience Division conducted a 3.5-year cooperative study from July 2012 to December 2015. The study characterizes the sequence stratigraphy, seismic stratigraphy, and seismic structures of the lower part of the intermediate confining unit aquifer and most of the Floridan aquifer system. Data obtained to meet the study objective include 80 miles of high-resolution, two-dimensional (2D), seismic-reflection profiles acquired from canals in eastern Broward County. These profiles have been used to characterize the sequence stratigraphy, seismic stratigraphy, and seismic structures in a 425-square-mile study area. Horizon mapping of the seismic-reflection profiles and additional data collection from well logs and cores or cuttings from 44 wells were focused on construction of three-dimensional (3D) visualizations of eight sequence stratigraphic cycles that compose the Eocene to Miocene Oldsmar, Avon Park, and Arcadia Formations. The mapping of these seismic-reflection and well data has produced a refined Cenozoic sequence stratigraphic, seismic stratigraphic, and hydrogeologic framework of southeastern Florida. The upward transition from the Oldsmar Formation to the Avon Park Formation and the Arcadia Formation embodies the evolution from (1) a tropical to subtropical, shallow-marine, carbonate platform, represented by the Oldsmar and Avon Park Formations, to (2) a broad, temperate, mixed carbonate-siliciclastic shallow marine shelf, represented by the lower part of the Arcadia Formation, and to (3) a temperate, distally steepened carbonate ramp represented by the upper part of the Arcadia Formation. In the study area, the depositional sequences and seismic sequences have a direct correlation with hydrogeologic units. The approximate upper boundary of four principal permeable units of the Floridan aquifer system (Upper Floridan aquifer, Avon Park permeable zone, uppermost major permeable zone of the Lower Floridan aquifer, and Boulder Zone) have sequence stratigraphic and seismic-reflection signatures that were identified on cross sections, mapped, or both, and therefore the sequence stratigraphy and seismic stratigraphy were used to guide the development of a refined spatial representation of these hydrogeologic units. In all cases, the permeability of the four permeable units is related to stratiform megaporosity generated by ancient dissolution of carbonate rock associated with subaerial exposure and unconformities at the upper surfaces of carbonate depositional cycles of several hierarchical scales ranging from high-frequency cycles to depositional sequences. Additionally, interparticle porosity also contributes substantially to the stratiform permeability in much of the Upper Floridan aquifer. Information from seismic stratigraphy allowed 3D geomodeling of hydrogeologic units—an approach never before applied to this area. Notably, the 3D geomodeling provided 3D visualizations and geocellular models of the depositional sequences, hydrostratigraphy, and structural features. The geocellular data could be used to update the hydrogeologic structure inherent to groundwater flow simulations that are designed to address the sustainability of the water resources of the Floridan aquifer system. Two kinds of pathways that could enable upward cross-formational flow of injected treated wastewater from the Boulder Zone have been identified in the 80 miles of high-resolution seismic data collected for this study: a near-vertical reverse fault and karst collapse structures. The single reverse fault, inferred to be of tectonic origin, is in extreme northeastern Broward County and has an offset of about 19 feet at the level of the Arcadia Formation. Most of the 17 karst collapse structures identified manifest as columniform, vertically stacked sagging seismic reflections that span early Eocene to Miocene age rocks equivalent to much of the Floridan aquifer system and the lower part of the overlying intermediate confining unit. In some cases, the seismic-sag structures extend upward into strata of Pliocene age. The seismic-sag structures are interpreted to have a semicircular shape in plan view on the basis of comparison to (1) other seismic-sag structures in southeastern Florida mapped with two 2D seismic cross lines or 3D data, (2) comparison to these structures located in other carbonate provinces, and (3) plausible extensional ring faults detected with multi-attribute analysis. The seismic-sag structures in the study area have heights as great as 2,500 vertical feet, though importantly, one spans about 7,800 feet. Both multi-attribute analysis and visual detection of offset of seismic reflections within the seismic-sag structures indicate faults and fractures are associated with many of the structures. Multi-attribute analysis highlighting chimney fluid pathways also indicates that the seismic-sag structures have a high probability for potential vertical cross-formational fluid flow along the faulted and fractured structures. A collapse of the seismic-sag structures within a deep burial setting evokes an origin related to hypogenic karst processes by ascending flow of subsurface fluids. In addition, paleo-epigenic karst related to major regional subaerial unconformities within the Florida Platform generated collapse structures (paleo-sinkholes) that are much smaller in scale than the cross-formational seismic-sag structures.

Florida

Elk monitoring in Lewis and Clark National Historical Park: 2008-2012 synthesis report

Maintaining elk (Cervus elaphus roosevelti) herds that frequent Lewis and Clark National Historical Park (NHP) is central to the park’s purpose of preserving the historic, cultural, scenic, and natural resources associated with the winter encampment of the Lewis and Clark expedition. Elk were critically important to the Lewis and Clark expedition in providing food and hides that sustained the expedition during the winter of 1805-06 and supplied them for their return east during 1806. Today, elk remain a key component of interpreting the Lewis and Clark story to over 200,000 park visitors each year at the Fort Clatsop visitor center. In 2008, the US Geological Survey (USGS) began collaborating with Lewis and Clark NHP and the NPS North Coast and Cascades Network to develop a protocol for monitoring long-term changes in the magnitude and spatial patterns of elk use within and adjacent to Lewis and Clark NHP (Griffin et al. 2011). Specific objectives of the monitoring program were to measure trends in (1) relative use of the Fort Clatsop unit by elk during winter; (2) the proportion of areas where elk sign is present in the Fort Clatsop unit in winter; and (3) the frequency of elk sightings from roads in and around the Fort Clatsop unit. This report synthesizes the results of the first four years of monitoring elk distribution and use in Lewis and Clark NHP from 2008-2012. We also present data from FY2012 (Appendix 1), in lieu of an annual report for that year. We used fecal pellet group surveys as the cornerstone for monitoring trends in both relative use of the Fort Clatsop Unit by elk and the proportion of areas where elk sign was present at the end of winter. We estimated pellet group density based on data collected from a network of fecal pellet plots distributed systematically throughout the unit. We developed a double observer sampling scheme that enabled us to estimate detection biases and improve the accuracy of pellet group density estimates. We computed the estimated detection probability for any pellet group observed; this probability was a function of the pellet group size and stage of decay, as well as lighting and vegetation conditions, and the number of observers (one or two) searching for pellets in that subplot. We then used these estimated detection probabilities to adjust the raw counts of the detected pellet groups to account for groups that likely went undetected under similar pellet and environmental conditions (each observed pellet group was weighted by the inverse of its estimated detection probability). We also used results from the late winter fecal pellet surveys to quantify the proportion of areas where elk pellets occurred (PAO), which was based on the presence of fecal pellet groups and estimation of detection biases (i.e., accounting for pellet groups that likely went undetected by both observers). In this synthesis, we report temporal trends in both pellet group density and PAO from 2008-2012, based on weighted linear regression analyses as well as spatial variation of pellet group densities over time. We completed late winter fecal pellet surveys at 61-66 plots annually, depending on yearly variation in access. We cleared fecal pellets at survey points in late October / early November each year and returned in late February / early March to count pellet groups left by elk over the winter. The estimated probability that a pellet group was detected by any one observer during late winter was affected most by the pellet group size and was less affected by decay class and lighting conditions. Per-observer detection probabilities ranged from as low as ~10-15% for single pellets to ~85-90% for pellet groups with 50 pellets. Average pellet group density in the Fort Clatsop unit ranged annually from 0.58 (+/- 1.43 standard error [SE]) to 0.93 (+/- 2.25 SE) pellet groups per 3-m radius subplot. Pellet group density declined over time, at approximately 8.8% per year (+/- 2.5% SE), but that slope was not statistically distinguishable from zero (2- tailed P=0.16). Following correction for detection biases, the proportion of surveyed points used by elk (i.e., PAO) ranged from 0.44 (+/- 0.07 SE) to 0.53 (+/- 0.07 SE) during the 4 winters. The estimated proportion of areas where elk pellets occurred (PAO) declined at a rate of 2.6% per year (+/- 1.2% per year SE), but that trend also was not statistically distinguishable from zero (2- tailed P=0.17). Statistical significance of a measure’s trend depends on both the magnitude (i.e., slope) of the observed trend and the number of years the trend continues in the same increasing or decreasing direction. Through simulation modeling we determined how many additional years of surveys would be required to reveal a statistically significant trend, based on the same trends in pellet group density and PAO, and associated variation, observed from 2009-2012. Assuming the same trends persist in the future, simulations indicated that there is a 70% probability that a statistically significant trend would be detected after two more years of conducting pellet group surveys. Relative use by elk during winter, as indexed by elk pellet group density, was generally greatest in the southeast region of the Fort Clatsop unit in or near the large freshwater marsh at the mouth of Colewort Creek and adjacent upland areas. Pellet group density was also higher than average in the north-central forested area, not far from a privately-owned pasture north of the park boundary. This spatial pattern in pellet group densities across the Fort Clatsop unit was consistent across all four years, although specific pellet group densities varied from year to year. Pellet group density declined significantly over time at two points in the southeast of the Fort Clatsop unit, even though pellet group density at those points remained higher than the unit average. Pellet group density increased significantly over time at one point in the north-central region, and at one point in the south-central region of the unit, indicating a slight shift in the distribution of elk use within the Fort Clatsop Unit over the four years. As an index of visitors’ opportunities to see elk in and around the Fort Clatsop Unit, we conducted replicated roadside elk surveys 3-5 times monthly during February, April, June, August, October and December 2008-2012. During each morning of survey, we searched for elk along four routes that totaled 32 km. We examined bimonthly trends in the numbers of elk groups seen, the total number of elk seen, and the observed composition ratios for those six months of the year. The average number of elk groups seen per survey ranged from 0.75 (+/- 0.32 SE) during February to a peak of 1.95 (+/- 0.36 SE) during June. Despite this seasonal variation in numbers of elk groups seen, the average total number of elk seen per morning was less variable. The average ratios of antlered elk to antlerless adult elk (i.e., bulls:cows) and calves to antlerless adult elk (i.e. calves:cows) varied seasonally, with the highest of both average ratios observed in August. We detected no significant trends in the average number of elk groups and total numbers of elk seen per survey from 2008-2012. Similarly, ratios of calves and antlered elk per antlerless elk did not differ over time. Elk groups were frequently seen from January to August in the southeast region of the Fort Clatsop unit, in the vicinity of Colewort Creek. Outside of NPS lands, we observed elk most frequently in open areas near the Astoria regional airport, in the pastures and forests immediately north of the Fort Clatsop unit and, prior to the construction of a residential development, in a pasture northwest of the Fort Clatsop unit. Elk monitoring at Lewis and Clark NHP is still in its initial years and additional monitoring will be required to verify trends that appear to be emerging. For example, the initial monitoring suggested incipient declining trends in both pellet group density and proportion of plots with pellets present, as well as, potentially, a small shift in elk distribution away from a new trail that was recently constructed in the southeast portion of the Fort Clatsop unit. Continued monitoring will aid in determining whether this local change in distribution persists (or, alternatively, resulted from short-term random variation), and whether there will be any positive or negative effect in the northern portion of the unit where a new trail has been constructed. High variability in road counts prevented our ability to find any clear trend in numbers or composition of elk observed in and near Fort Clatsop, but changes in the patterns of observations of elk from roadways suggest that residential development outside the park has reduced the available habitat for elk in some of the areas surveyed, and may have affected spatial use patterns of elk adjacent to some areas of the park. In addition to monitoring future effects of land use changes outside the park, continued monitoring may also prove useful for assessing elk responses to natural succession in forests disturbed by windthrow in December 2007 and to NPS vegetation management activities such as variable density thinning in the forest, trail development, and restoration at Otter Point tidal area and Colewort Creek Slough.

Oregon;Washington

Coupled arsenotrophy in a hot spring photosynthetic biofilm at Mono Lake, California

Red-pigmented biofilms grow on rock and cobble surfaces present in anoxic hot springs located on Paoha Island in Mono Lake. The bacterial community was dominated (∼ 85% of 16S rRNA gene clones) by sequences from the photosynthetic Ectothiorhodospira genus. Scraped biofilm materials incubated under anoxic conditions rapidly oxidized As(III) to As(V) in the light via anoxygenic photosynthesis but could also readily reduce As(V) to As(III) in the dark at comparable rates. Back-labeling experiments with 73 As(V) demonstrated that reduction to 73 As(III) also occurred in the light, thereby illustrating the cooccurrence of these two anaerobic processes as an example of closely coupled arsenotrophy. Oxic biofilms also oxidized As(III) to As(V). Biofilms incubated with [ 14 C]acetate oxidized the radiolabel to 14 CO 2 in the light but not the dark, indicating a capacity for photoheterotrophy but not chemoheterotrophy. Anoxic, dark-incubated samples demonstrated As(V) reduction linked to additions of hydrogen or sulfide but not acetate. Chemoautotrophy linked to As(V) as measured by dark fixation of [ 14 C]bicarbonate into cell material was stimulated by either H 2 or HS − . Functional genes for the arsenate respiratory reductase ( arrA ) and arsenic resistance ( arsB ) were detected in sequenced amplicons of extracted DNA, with about half of the arrA sequences closely related (∼98% translated amino acid identity) to those from the family Ectothiorhodospiraceae . Surprisingly, no authentic PCR products for arsenite oxidase ( aoxB ) were obtained, despite observing aerobic arsenite oxidation activity. Collectively, these results demonstrate close linkages of these arsenic redox processes occurring within these biofilms. Oxyanions of the group 15 element arsenic, arsenate [As(V)] and arsenite [As(III)], have been known for millennia to be potent poisons. Despite its well-established toxicity to life, the phenomenon of arsenic resistance was discovered whereby some microorganisms maintain an otherwise “normal” existence in the presence of high concentrations of As(V) or As(III) (17, 29, 31). More recently it has become recognized that certain representatives from the bacterial and archaeal domains can actually exploit the electrochemical potential of the As(V)/As(III) redox couple (+130 mV) to gain energy for growth. This can be achieved either by employing As(III) as an autotrophic electron donor or by using As(V) as a respiratory electron acceptor (18, 21, 34). The latter phenomenon, although most commonly associated with chemoheterotrophy, can also employ inorganic substances like sulfide or H 2 . Indeed, As(V)-respiring anaerobes displaying a capacity for chemoautotrophy with these electron donors have been isolated and described (5, 7, 16). We recently reported that photoautotrophy is supported by As(III) in anoxic biofilms located in hot springs on Paoha Island in Mono Lake, CA (15). This process represented a novel means of As(III) oxidation achieved via anoxygenic photosynthesis occurring in certain photosynthetic bacteria (i.e., Ectothiorhodospira ) and possibly within some cyanobacteria as well (e.g., “ Oscillatoria ”). Whether or not a microbial habitat is overtly oxic or anoxic, or temporally shifts between these two states over a diel cycle, critical energy linkages between aerobes and anaerobes have long been known for the biogeochemical cycles of key elements, such as sulfur, iron, and nitrogen. Most prominently studied is the case of nitrogen, whereby an ecological coupling exists between the processes of nitrification and denitrification (9, 10, 28). The former process provides energy to aerobic nitrifiers, while the latter process consumes the nitrate produced by this reaction, thereby meeting the energy needs of the denitrifiers. For arsenic, the detection of both As(III) oxidation and As(V) reduction in oxic and anoxic incubations of freshly collected periphyton suggested that an analogous coupled process may also occur for this element (12). Similarly, several uncontaminated soils in Japan displayed a capacity for either As(V) reduction or As(III) oxidation upon arsenic oxyanion amendment and whether they were incubated under oxic or anoxic conditions (39). A defined coculture consisting of an aerobic As(III) oxidizer (strain OL1) and an anaerobic As(V) respirer (strain Y5) was shown to function in this fashion under manipulated laboratory conditions of oxygen tension (26). We pursued the phenomenon of coupled arsenic metabolism further by using materials collected from the hot spring biofilms in Mono Lake, but we focused on examination of the cycling of arsenic under anoxic conditions. In this paper we report results obtained by manipulated incubations of red-pigmented biofilms found in the hot springs of Paoha Island. Preliminary community characterizations of these biofilms show that they are dominated by B acteria from the genus Ectothiorhodospira but also harbor an assemblage of Archaea related to the Halobacteriacaea . Incubation results have demonstrated the presence of the following arsenic metabolic activities: respiratory As(V) reduction, photosynthetic anaerobic As(III) oxidation, and aerobic As(III) oxidation, along with the ecophysiological conditions under which they occur. Surprisingly, we were unable to obtain authentic PCR products for arsenite oxidase genes ( aoxB ), despite observing aerobic As(III) oxidation activity. These biofilms serve as a model system for how anaerobic cycling of arsenic can be sustained with oxidation of As(III) by anoxygenic photosynthesis coupled to regeneration of this electron donor via dissimilatory As(V) reduction. The significance that such a light-driven anaerobic ecosystem may have played in the Archean Earth is discussed.

California