USGS Science⌕ Search

SEARCH · USGS Science

Results for “Environmental Research Letters”

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.

504 records · Page 28Linked to original sources
USGS2025

Scientific opportunities in the National Landscape Conservation System

The National Landscape Conservation System consists of unique and beautiful places across America’s landscapes where identified resources and values are protected and science is highlighted. The mission of the National Landscape Conservation System (NLCS), which is managed by the Bureau of Land Management and is often referred to as the agency’s National Conservation Lands, is to conserve, protect, and restore nationally significant landscapes for their cultural, ecological, and scientific values. This clear inclusion of science in the NLCS mission sets the stage for individual units to serve as places of learning, teaching, discovery, and innovation. Science is an integral part of managing the National Conservation Lands, and science conducted within and across the more than 900 units that make up the NLCS can inform and influence conservation and public land management well beyond its boundaries. Here, we highlight seven core aspects of National Conservation Lands that present valuable science opportunities: (1) the scientific values for which individual units are designated; (2) the many other resources, objects, and values within units; (3) the value of units as “control” sites for understanding the effects of activities such as mineral extraction that commonly occur elsewhere on multiple-use public lands but are often prohibited within National Conservation Lands; (4) the value of units for studying the effects of activities such as recreation that regularly occur and may be intensified on National Conservation Lands; (5) the high visibility of units, which draws strong interest and engagement from scientists, partners, and the public; (6) the functioning of the units as a network managed for a common purpose, which provides an opportunity to explore cross-cutting science questions across widely varying contexts and geographies; and (7) the opportunities units provide to promote and apply Indigenous Knowledge to scientific research to manage natural and cultural resources. Because of all of these characteristics, National Conservation Lands can serve as hubs for basic and applied science that can inform management of all public lands and resources into the future. We highlight these science opportunities through examples from existing units and suggest two actions that could help further science activities and impact on National Conservation Lands.

Parks Stewardship Forum↗
USGS2024

Influence of inherited structure on flexural extension in foreland basin systems: Evidence from the northern Arkoma basin and southern Ozark dome, USA

Extensional faults are key components of foreland basin systems. They form within the upper crust in response to flexure of the lithosphere and accommodate subsidence within the foredeep and forebulge depozones. Such faults are excellent proxies for orogenic system evolution and control the distribution of natural resources and hazards. However, the spatiotemporal evolution of flexural extension has not been documented previously at a regional scale, thereby limiting our understanding of underlying geodynamic controls. Here, we resolve late Paleozoic flexural extension in the northern Arkoma basin and southern Ozark dome, USA. We synthesize a large database of previous mapping, existing research, subsurface data, and geophysical data into 3D geologic and 2D kinematic models. Mesh surfaces representing several key horizons from the Carboniferous Period (ca. 335-306 Ma) were constructed. These surfaces were built from oil and gas well tops (n = ∼10,000) and surface geologic map contacts using an advanced kriging method. The mesh surfaces are offset by a complex 3D fault network, allowing detailed analysis of along-strike and down-dip variations in fault displacement. Analysis of the 3D model reveals a regular and repeated fault segmentation pattern wherein E -W striking, left- and foreland-stepping en échelon normal faults are segmented by inherited NE striking basement faults. Maximum vertical separation along the E -W normal faults is generally focused between the inherited NE-trending faults. This suggests that the inherited basement faults delocalized extensional strain during late Paleozoic normal faulting. Maximum vertical separation and fault localization may correlate to areas with high-amplitude positive magnetic anomalies interpreted as Mesoproterozoic granitic rocks. Speculative covariance of magnetic anomalies and fault displacements implies that the relatively strong basement granite concentrated stress, leading to localized faulting within the relatively thin sedimentary cover. Lastly, we show that flexural extension migrated southeast to northwest from the Chesterian-Morrowan (ca. 335-319 Ma) to the Desmoinesian (ca. 306 Ma). The migratory flexural extension may be explained by diachronous loading during Pangean assembly, or by synchronous loading but variable load compensation due to inherent factors.

Arkansas, Oklahoma↗
USGS2025

Metal fingerprints of Eocene rhyolite magmas coincident with Carlin-type gold deposition in Nevada USA

Eocene magmatic systems contemporaneous with world-class Carlin-type Au deposits in Nevada (USA) have been proposed by some researchers as a key ingredient for Au mineralization, though evidence conclusively demonstrating their genetic relationship remains tenuous. This study provides the first direct evidence of the pre-eruptive metal budget of volatile- and metal-charged silicic magmas coincident in time (~41 to 34 Ma) and space (within 5 km) with Carlin-type Au deposits. We characterize the pre-eruptive metal fingerprints of these diverse magmatic systems to assess their potential as sources of metals for Carlin-type Au mineralization. Metal abundances from quartz-hosted melt inclusions (Au, Te, Ag, Sb, Tl, Mo, W, Sn, As, Pb, Co, Cu, Ni, and Zn) characterized in situ by SHRIMP-RG and LA-ICP-MS represent our best (and only) estimates for the pre-eruptive metal budget in these systems. Median metal concentrations are generally within one order of magnitude of average upper crust and average continental rhyolite values. But there are two notable exceptions, with median Au contents extending >1 order of magnitude higher than average upper crust and median Cu contents ranging >1 order of magnitude lower than upper crust. Despite this, melts contain lower Au/Cu (<0.1), Au/Ag (<5), and Au/Tl (<0.3) than most ore-grade Carlin-type rock samples and quartz-hosted fluid inclusions, regardless of their age and timing relative to nearby Carlin-type Au mineralization. The metal fingerprints of these magmatic systems, de-fined both by traditional and multivariate compositional data analysis techniques, are distinct from one another. Yet none are particularly specialized, e.g., high Au/Cu, in terms of being ideal ingredients as postulated by magmatic models for Carlin-type Au mineralization. Magmatic Au contents do not appear to be correlated with rhyolite “flavors” in the way that Cu, Sn, and Nb contents are. Fluid/melt partitioning modeling and magma volume estimates support the idea that a diverse array of non-specialized silicic magmas could feasibly contribute some or potentially all of the Au, Ag, and Cu in Carlin-type systems. The compositional diversity among contemporaneous magmatic systems could possibly contribute to some of the diversity observed across Carlin-type Au districts in Nevada.

Nevada↗
USGS2025

Water supply in the conterminous United States, Alaska, Hawaii, and Puerto Rico, water years 2010–20

We present an assessment of water supply across the conterminous United States (CONUS), Alaska, Hawaii, and Puerto Rico covering water years 2010–20. Our analysis drew on two national hydrologic models, the National Hydrologic Model Precipitation-Runoff Modeling System and the Weather Research and Forecasting model hydrologic modeling system. Both models produced estimates of streamflow, evapotranspiration, soil moisture, snow water equivalent, and other hydrologic states and fluxes. The models were driven by the bias-adjusted 4-kilometer-resolution, long-term regional hydroclimate simulation over the conterminous United States dataset (CONUS404). We assessed spatial and temporal error distributions by comparing monthly simulations at the 12-digit hydrologic unit code and regional scale from both models against external benchmarking datasets. Results showed that average annual rainfall across the CONUS was 857 millimeters per year for the period of analysis, with water year 2012 the driest year (729 millimeters) and water year 2019 the wettest year (995 millimeters). Key interannual variability results included the following: (1) the California–Nevada hydrologic region had the highest variability in precipitation and snow accumulation, and (2) the Texas hydrologic region was among hydrologic regions with the highest variability in precipitation. We related interannual variability in precipitation to storage volumes in soil moisture, snow water equivalent, and lakes and reservoirs to highlight areas with little storage and large year-to-year variability in precipitation. These areas included the Southern High Plains, Central High Plains, Texas, Souris–Red–Rainy, Mississippi Embayment, and Midwest regions. Our analysis of groundwater-level data showed that several of these areas overlap aquifers where groundwater levels were considerably lower than historical averages, including the Colorado Plateaus aquifers, the Rio Grande aquifer system, and the Central and Southern regions of the High Plains aquifer. Many of these lowered groundwater levels are continuations of decades-long declines from overpumping that started well before the assessment period. The resulting water budgets and their analyses provide a high-resolution foundational assessment of the mean state and variability of the terrestrial hydrologic cycle across the CONUS and Alaska, Hawaii, and Puerto Rico to support a wide range of water resource management applications.

Professional Paper↗
USGS2026

Comparative assessment of STIC sensors, streamflow and rain gauges for quantifying river connectivity in intermittent systems

In intermittent stream systems, including those occurring in Texas, USA, the severity of low-flow conditions, duration of seasonal disconnection, and frequency of no-flow events have been amplified by drought. Documentation of these no-flow events is necessary to evaluate ecosystem health. However, many intermittent reaches remain un-gauged given that perennial river sec-tions are often prioritized for gauge placement. Our objectives were to 1) document stream flow using Stream Temperature, Intermittency, and Conductivity (STIC) loggers to determine the frequency and duration of no-flow events in intermittent tributaries of the Colorado River, Texas and 2) compare logger data to publicly available data from streamflow discharge and precipitation gauge networks to understand differences among these data types for drying event characterization. We use these comparisons to summarize benefits and limitations of the application of in-stream data loggers. STIC loggers were deployed at 19 sites, one in each pool and riffle habitat of a stream reach. STIC loggers recorded a measurement of relative conductance every six hours from June 2022 to March 2024, which was used to determine the presence or absence of flow connectivity in a reach. No-flow duration among intermittent reaches varied between 37 and 270 days across tributaries during an ongoing drought in the study area. Overall, logger data was more precise than discharge data for characterizing no-flow events or precipitation data when documenting presence of water in the stream channel due to runoff. Lack of discharge gauges in intermittent tributaries left large sections of stream reaches undocumented and resulted in mischaracterization of flow patterns. Drought severity across the tributaries did not follow longitudinal patterns that would be expected by the climatic precipitation gradient of the study area. More research is needed to determine if factors such as population size affect severity. Likewise, precipitation data did not correlate well with logger water presence data, lacking consideration for groundwater recharge, soil hydrophobicity, and surface compaction. This study shows that to monitor no-flow events, detailed spatial datasets are necessary and that STIC loggers are useful tools that provide data to fill spatial information gaps and facilitate more accurate flow characterization and water presence data in intermittent systems.

Texas↗
USGS2025

Update of the 2008 provisional Enhanced Geothermal Systems (EGS) assessment for the Great Basin, USA

In response to the Energy Act of 2020, the U.S. Geological Survey (USGS) is updating the Enhanced Geothermal Systems (EGS) resource assessment for the Great Basin, USA. The previous 2008 provisional assessment estimated how much electricity could be generated from EGS resources of the western United States using models of electric-grade heat, models of heat extraction over time, and estimates of how much rock might be stimulated to produce viable amounts of heat. Herein, a similar conceptual strategy is applied, using updated models of heat extraction as a function of fracture spacing and well distance. Previously used reservoir heat delivery models are updated to have a dependence on fracture and well spacing, potentially improving future estimates of EGS resources as ongoing research provides a better understanding about the success of reservoir stimulation as a function of geology and location. For a range of well distances (250-1000 m) and fracture spacings (1-50 m), heat extraction efficiency ranges from 25-62%, demonstrating the importance of accounting for the most likely results of proven viable fracturing technologies. Although fracturing is important, the biggest uncertainty by far in estimating the EGS resource for the Great Basin is estimating which geologic units at what depths can be stimulated sufficiently to produce geothermal energy economically and efficiently. Uncertainties in these factors yield estimates that range over two orders of magnitude with an upper limit of ~174 terawatts-thermal (TWth) produced for 30 years from the upper 7 km of the crust. This upper limit would require significant technological advances to access most of the electric-grade resource across the Great Basin. Assuming that 1% of this estimate will be accessible in the next few decades gives a resource estimate similar to that made in the 2008 provisional assessment. These estimated EGS heat extraction rates far exceed ( greater than 100x) the natural geothermal heat production rate, thus geothermal electricity production at these rates might not be sustainable unless heat is also recharged from other sources (e.g., excess solar energy when supply exceeds demand). In addition to assessment maps and cumulative estimates, the new models of fractured reservoirs developed herein can be used to estimate steady power production given a set of fractures and well spacing, and estimates can be made for setback distances to ensure no thermal interference with nearby powerplants.

Arizona, California, Idaho, Nevada, Oregon, Utah↗
USGS2025

Wangyanite, PdNi8S8, a new Pd end-member mineral of the pentlandite group from the J-M reef, Stillwater Complex, Montana, USA

Wangyanite (IMA2024-008a), ideally PdNi 8 S 8 , is a Pd end-member mineral of the pentlandite group that was discovered in the J-M reef of the Stillwater Complex, Montana, USA. Wangyanite occurs as anhedral-subhedral granular crystals 200–400 µm in size, associated with isoferroplatinum, braggite, pentlandite, and chalcopyrite interstitial to plagioclase grains within anorthosite. Wangyanite exhibits a yellowish brown color with a black streak and a metallic luster. It is brittle with uneven fractures, and has a calculated density of 5.14 g/cm 3 . The mineral does not show discernible pleochroism, bireflectance, or anisotropy. It has an average composition of 9.95 wt.% Pd, 31.95 wt.% Ni, 25.02 wt.% Fe, 0.57 wt.% Co, 31.74 wt.% S, totaling 99.23 wt.%. The empirical formula, based on eight sulfur atoms per formula unit, is (Pd 0.76 Co 0.08 ) Σ0.84 (Ni 4.39 Fe 3.60 ) Σ7.99 S 8 . Wangyanite has a cubic cell with a space group of Fm -3 m (#225), having lattice parameters of a = 10.1167(12) Å, V = 1035.4(4) Å 3 , and Z = 4. Its crystal structure has been solved by single-crystal three-dimensional electron diffraction study. The strongest X-ray diffraction lines of wangyanite are claculated at [ d in Å ( I %)( hkl )]: 5.841(14.03)(111), 3.050(100)(311), 1.947(29.16)(115,333), 1.264(11.66)(800), 3.577(8.79)(220), 2.920(20.82)(222), and 2.321(9.34)(331). Wangyanite shares the same crystal structure as pentlandite, but the octahedrally coordinated site is mainly occupied by Pd in wangyanite. Based on the textural features and previous experimental Pd-Fe-Ni-S phase system, wangyanite could form by peritectic reaction between braggite, pentlandite and sulfide liquid. These mineral associations are stable in a Ni-Pd-rich sulfide melt system at about 550 °C or even lower temperature. Therefore, wangyanite can potentially serve as an indicator of the presence of Pd-rich residual melts. The mineral is named in honor of Prof. Christina Yan Wang, a well-known researcher on platinum-group element (PGE) occurrences and enrichment mechanisms in mafic-ultramafic intrusions, notably those deposits related to the Emeishan large igneous province in China.

Montana↗
USGS2026

Ungulate migrations of the Western United States, volume 6

This report, volume 6 in the “Ungulate Migrations of the Western United States” report series, showcases the migrations of 23 ungulate herds in the Western United States. The report series is produced by the Corridor Mapping Team (CMT). Led by the U.S. Geological Survey, the CMT is a collaboration among 11 State agencies, as well as regional and Federal partners, and an expanding number of Tribal wildlife agencies. The CMT was initiated in response to the U.S. Department of the Interior Secretarial Order 3362, which was signed in 2018 and provided Federal support to expand existing research efforts to study ungulate populations and conserve their migrations throughout the Western United States. Including this volume, the report series has detailed the migrations of 237 unique ungulate herds throughout the Western United States and continues to serve as a valuable resource to guide local and regional management, policy, and on-the-ground work necessary to maintain intact and functional ungulate migrations. This report highlights several guiding principles of the CMT that facilitate collaboration among the diverse set of partners and contribute to the program’s continued successes. Notably, raw global positioning system data are not shared among participating agencies and the U.S. Geological Survey, delineating migration corridors and seasonal ranges relies on empirical data, the CMT provides flexible approaches to participating State and Tribal partners, and regular CMT meetings create a framework for open communication among agency partners that supports transboundary mapping of migrations. The 237 ungulate migrations that have been included in the report series are an expanding inventory, which can help maintain ungulate migrations in perpetuity.

Arizona, California, Colorado, Idaho, Montana, Nev↗
USGS2025

Factors affecting the distribution of water-bearing fractures in the bedrock aquifers of West Virginia

Bedrock aquifers cover 23,601 square miles within the State of West Virginia and comprise 97.4 percent of the surficial area within the State; the remaining 2.6 percent (621 square miles) consists of alluvial sand-and-gravel and glacial outwash aquifers bordering the State’s major rivers. While West Virginia’s alluvial aquifers have been studied extensively, bedrock aquifers have only been characterized for studies completed in a few areas in Jefferson, McDowell, and Monroe Counties. Bedrock aquifers are water supplies for public supply, agriculture, industry, and residential homeowner use. In this study, the U.S. Geological Survey, in cooperation with the West Virginia Department of Environmental Protection Division of Water and Waste Management, provides a statewide assessment of the occurrence and distribution of fractures within bedrock aquifers of the State and the various topographic, physiographic, and lithologic influences controlling the occurrence and distribution of bedrock fractures. The results of this study provide an increased understanding of the distribution of fractures in bedrock aquifers in West Virginia and help to verify trends that have been suspected for many years but were never well documented or verified by data. The results confirmed that the density of fractures and those that were determined to be water bearing decrease significantly with depth. A statistically significant difference in the density of fractures was observed at a depth of 215 feet for wells in the Appalachian Plateaus Physiographic Province’s and in the Valley and Ridge Physiographic Province’s aquifers; a higher density of fractures and water-bearing fractures were above a depth of 215 feet than below that depth. This is an important consideration when drilling wells for residential, commercial, industrial, or agricultural water supply. Abandoned underground coal mines are commonly believed to form large pools of water in the interconnected mine entries in abandoned room and pillar coal mines. Such pools of water can and do exist in abandoned underground coal mines, but many mines lack open entries and are held up by overburden strata and pillars that can collapse and form aquifers comprised of vast interconnected rubble zones (gob), especially in older mines. Data assessed for this study showed that shale-corrected values of effective porosity for limestone aquifers in West Virginia had a median value of 2 percent and an average value of 4 percent and generally are mineralized with low effective porosity. Argillaceous or sandy limestone has a median shale-corrected porosity of 4 percent and an average shale-corrected porosity of 5 percent. The median and average shale-corrected porosity of sandstone aquifers was estimated to be 14 percent, but the median shale-corrected porosity for argillaceous or calcareous sandstone was 5 percent and the average shale-corrected porosity for argillaceous or calcareous sandstone was 6 percent. Even though shale has a relatively high total sonic porosity compared to other lithologies, shale and siltstone had relatively low shale-corrected porosity, ranging from 0 to 2 percent. Well yields were previously documented to be highest in valley settings, lowest on hilltops, and intermediate on hillsides. Transmissivity data provided by this study confirm this general pattern within the Appalachian Plateaus Province; however, the Valley and Ridge Province does not follow this pattern. While still lowest on hilltop settings, the highest well yields were in hillside settings. The trend for the Valley and Ridge Province was likely skewed because of 9 high-yield wells specifically targeting deeper thin limestone units, such as the Tonoloway and Helderberg Limestones, at depths with transmissivity in excess of 2,000 feet squared per day in Mineral County, West Virginia, or targeting karst aquifers in Berkeley, Jefferson, or Greenbrier Counties, West Virginia. Finally, water-bearing fractures have been hypothesized to comprise a small number of all fractures within a typical bedrock well in West Virginia. Data collected for this study support this theory. A total of 3,403 fractures were identified during this study; 3,151 (92.6 percent) of those fractures are low-transmissive, and only 252 (7.4 percent) fractures are water-bearing. Even though a well may contain many fractures, less than 8 percent are considered water-bearing fractures.

West Virginia↗
USGS2025

Benthic habitat map of Olowalu Reef, Maui, Hawaii—Geomorphological structure, biological cover, and geologic zonation determined with spectral, lidar, and acoustic data

The fringing coral reef off Olowalu, Maui, Hawaii, has been identified as a local conservation priority site. In 2007, the National Oceanic and Atmospheric Administration (NOAA) produced a benthic habitat map of the Hawaiian Islands that was used as a foundation for this study. To support place-based management of the reef in the future, the U.S. Geological Survey (USGS) mapped the geologic zone, major and dominant geomorphological structure, biological cover type, and percent of biological cover for 11 square kilometers (km 2 ) of Olowalu Reef at a minimum mapping unit (MMU) of 100 square meters (m 2 ) to create a benthic habitat map. Heads-up digitization was employed on 0.50-meter (m) natural color satellite orthoimagery with ancillary 1-m acoustic backscatter imagery from single-scan sonar (sound navigation and ranging). A 1-m, 4-m, and 8-m digital bathymetric model (DBM) was interpolated from bathymetric lidar (light detection and ranging), and various geomorphometric layers derived from the DBMs were used for habitat interpretation. Still-frame imagery of the seafloor extracted from vessel-towed underwater video transects on Olowalu Reef served as ground validation points ( n =870) during active mapping and accuracy assessment points ( n =216) for thematic accuracy assessment. Thematic accuracy was cross-validated by the Hawai‘i Department of Land and Natural Resources Division of Aquatic Resources. Final thematic accuracy was 88.8 percent for major structure, 85.6 percent for dominant structure, 86.0 percent for major biological cover, and 78.6 percent for type and percent of major biological cover. Reef and hardbottom constituted 52 percent of the total mapped habitat, comprising mostly aggregate reef (31 percent) and pavement (11 percent), with large swaths of spur-and-groove (9 percent). Of this hardbottom, 17 percent was covered with moderate (10 to <50 percent) coral and 27 percent with high coral cover (50 to <90 percent). High (50 to <90 percent) macroalgae cover dominated the continuous sand sheets in offshore bank/shelf zones. The map created in this study supplements the NOAA 2007 map and expands on the observations made by USGS sampling of the reef. The NOAA 2007 map and our map differed in total areal extent by a negligible 6 m 2 and were in general thematic agreement. Our map is intended to serve as a baseline for public access, general research, local-level management, and reef change for future studies.

Hawaii↗
USGS2026

Distributed volcanic fields—Patterns, products, processes, and plumbing

Distributed volcanic fields—also referred to as “monogenetic fields” or simply “volcanic fields”—are found on every continent, on the seafloors, and in every tectonic setting on Earth. Yet, our understanding of what constitutes a distributed volcanic field, and our knowledge of the common products, spatiotemporal patterns in activity, and plumbing systems of these fields, are commonly poor and (or) debated. Here we provide an overview of the current understanding of distributed volcanic fields (DVFs) through review and compilation of published data, with an emphasis on understanding the characteristics of fields that are essential to prepare for future hazardous volcanic activity. Distributed fields can erupt restricted to variable magma compositions. Activity in these fields produces nearly all types of volcanoes, including scoria and spatter cones, maars, low shields, lava fields, stratovolcanoes, and calderas, but monogenetic volcanoes are the defining characteristic of DVFs. The number of vents and lifespans of DVFs vary greatly; vent numbers range from 2 to greater than 1,000, and fields can remain active for 150,000 to 10 million years. Some characteristics of DVFs appear field dependent. Although the volumes of single monogenetic eruptions within DVFs globally are variable (from less than 0.1 to a few cubic kilometers), each field may have a characteristic range in eruption volumes. Eruptions within distributed volcanic fields can have a random distribution and timing, but spatial and (or) temporal vent clustering and alignment are common, and eruptions may be time or volume predictable. Average eruption recurrence intervals in DVFs span over an order of magnitude, from less than 3 to greater than 100 thousand years. Exposures of the shallow portions of older systems reveal that many magma batches end their ascent before reaching the surface, with some dikes stopping within only meters of the surface. Petrologic insights into the deeper plumbing systems of DVF volcanoes indicate variable depths of magma crystallization and (or) stalling, and that magma transport and storage systems may be unique to a particular DVF and its underlying crustal structure. Both seismic and petrologic data for mafic DVF eruptions corroborate intrusions of magma into storage systems in the months to years before eruption, as well as variable ascent rates through the crust (days to weeks). Few distributed fields are closely monitored, as they erupt infrequently, but these observations imply that the amount of time available to install monitoring instrumentation after the first felt seismicity may be short, as precursory activity is frequently on the order of days to months for modern DVF eruptions. Satellite-based geodesy may provide useful warning in such cases. This chapter concludes with outstanding questions and areas for future research that will aid in our understanding and preparedness for future DVF eruptions.

Professional Paper↗
USGS2025

Preserving and increasing water resources—Natural infrastructure in dryland streams in Baja California Sur, Mexico

The Los Planes watershed of Baja California Sur, Mexico, and its underlying aquifer are experiencing groundwater decline owing to low average annual rainfall (28.1 centimeters per year) and rising water demand from population growth and agricultural activities. This decline in water availability can lead to desertification—a process that changes arable land to desert by degrading soil and vegetation—and can pose serious challenges to livelihoods that depend on the land. To address these issues, a ranch in the Los Planes watershed has installed many natural infrastructures in dryland streams (NIDS) in channels for soil and water conservation. In 2022, the U.S. Geological Survey (USGS) began working with regional researchers and land managers to investigate the effects of NIDS on natural biological, geochemical, and physical processes and determine the efficacy of NIDS for water augmentation in the Los Planes watershed. The USGS also worked with local academic institutions and nonprofit organizations to create public educational opportunities focused on the area’s hydrogeology. These and other collaborative efforts with the U.S. Water Partnership and Innovaciones Alumbra aim at enhancing water resources in the Baja California Sur region and promoting water security and safeguarding community well-being. Resume La cuenca de Los Planes, ubicada en Baja California Sur, México, y su acuífero subyacente, están sufriendo una disminución de las aguas subterráneas debido a la baja precipitación media anual (28.1 centímetros por año) y la alta demanda de agua por parte de una población creciente y la actividad agrícola. Esta disminución de la disponibilidad de agua puede conducir a la desertificación—un proceso que por medio de la degradación del suelo y la vegetación convierte a la tierra cultivable en desierto—representando un serio desafío para los medios de vida de las personas. Para abordar estos problemas, un rancho en la cuenca de Los Planes ha instalado numerosas obras de Infraestructura Natural en Arroyos de Tierras Áridas (INATS) para conservación del suelo y del agua. En 2022, el Servicio Geológico de los Estados Unidos (USGS, por sus siglas en inglés) comenzó a trabajar con investigadores regionales y gestores de tierras para estudiar los efectos de INATS en los procesos biológicos, geoquímicos y físicos, y determinar su eficacia en el aumento de los recursos hídricos en la cuenca de Los Planes. El USGS se ha asociado con instituciones académicas y organizaciones locales sin fines de lucro para crear oportunidades educativas públicas centradas en la hidrogeología de la zona. Estos y otros esfuerzos colaborativos con la Asociación del Agua de Estados Unidos (U.S. Water Partnership) e Innovaciones Alumbra, tienen como objetivo mejorar el uso de los recursos hídricos en la región de Baja California Sur, promover la seguridad hídrica y proteger el bienestar de la comunidad.

Baja California Sur↗
USGS2024

Post-wildfire debris flows

Post-wildfire debris flows pose severe hazards to communities and infrastructure near and within recently burned mountainous terrain. Intense heat of wildfires changes the runoff characteristics of a watershed by combusting the vegetative canopy, litter, and duff, introducing ash into the soil and creating water repellant soils. Following wildfire, rainfall on bare ground is less able to infiltrate into the fire-altered soils and overland flow is less impeded by vegetation. Rainfall runoff in recently burned areas can erode hillslopes owing to the removal of soil binding organic matter near the soil surface by fire. In channels, loose, dry-ravel deposits composed of sand and gravel are readily entrained by concentrated runoff in channels. Entrainment of soil on hillslopes and in channels bulks up the sediment concentration of the rainfall runoff to generate debris flows capable of transporting boulders and large woody debris. Post-wildfire debris flows can be triggered by rainfall conditions that would typically produce little runoff during unburned conditions. The primary rainfall trigger for post-wildfire debris flows is high intensity rainfall during short duration convective rainstorms or periods of high rainfall intensity embedded within a long-duration frontal storm. Numerous observations of debris flows triggered by storms lasting less than an hour following periods of little to no rainfall indicate that antecedent rainfall is not a requirement for initiation of post-wildfire debris flows. Post-wildfire debris-flow hazard assessment entails estimating probability and magnitude of debris flows in the burned area, estimating debris-flow runout and intensity, and defining rainfall intensity-duration thresholds for debris-flow initiation. In the United States, probability and magnitude is estimated using empirically derived models largely based on data collected in southern California. The models provide maps to identify watersheds and drainage paths where post-wildfire hazards are most pronounced. Rainfall intensity-duration thresholds can be incorporated into flood hazard forecasting tools. Currently, work is underway to identify how to best implement debris-flow runout models in burned areas with efficiency and accuracy. Post-wildfire debris flows have been a long-recognized process in the Transverse Ranges of southern California; however, climate change is driving more frequent wildfires to burn more mountainous terrain throughout the western United States and worldwide. As a result, post-wildfire debris flows are becoming a more common threat in areas where they were once infrequent. As the threat of post-wildfire debris flow expands into new areas, evaluating the hazard becomes challenging because the degree to which wildfire increases debris-flow susceptibility varies from region to region. This chapter summarizes the knowledge to date for evaluating post-wildfire debris-flow susceptibility and hazard assessment. We summarize the characteristics of wildfire burn severity, topography, underlying soil and geology, and rainfall conditions that contribute to making a watershed most likely to produce post-wildfire debris flows. Methods for hazard assessment in the United States and other countries are summarized. We highlight knowledge gaps for how post-wildfire debris-flow susceptibility varies throughout the western United States and worldwide and identify research needs to improve hazard assessment methods in different geographies.

Book chapter↗
USGS2024

Diel temperature signals track seasonal shifts in localized groundwater contributions to headwater streamflow generation at network scale

Groundwater contributions to streamflow sustain aquatic ecosystem resilience; streams without significant groundwater inputs often have well-coupled air and water temperatures that degrade cold-water habitat during warm low flow periods. Widespread uncertainty in stream-groundwater connectivity across space and time has created disparate predictions of energy and nutrient fluxes across headwater networks, hindering predictions of cold-water habitat resilience under climate change scenarios. Recently, annual paired air and water temperature signals have been harnessed to indicate stream water thermal sensitivity and the dominance of deep versus shallow groundwater influence, although the utility of diel air–water temperature signal metrics for hydrologic inference has remained unexplored. Here we analyzed two consecutive years of locally paired, air–water temperature data from 47 headwater stream sites in the Catskill Mountains, New York, USA, and discovered characteristic seasonal patterns in diel temperature signal sinusoid metrics (amplitude ratio, phase lag, and mean ratio) driven by shifts in streamflow generation mechanisms and stream network position. Hydrologic interpretations of observed patterns were supported by stream heat budget model scenarios and additional analysis of paired air–water temperature data from two streams in Shenandoah National Park, Virginia, USA, with well characterized stream-groundwater connectivity. We found that within smaller tributaries, streamflow generation transitions from runoff to groundwater dominance were driven by hillslope drying during seasonal periods of lower precipitation. This was evidenced by significant correlations (p < 0.01) between daily water:air temperature signal amplitudes (non-linear decreases of ∼ 50 %) and derived base-flow index at 22 of the 28 sites, indicating enhanced local groundwater influence on streamflow promotes decoupling of diel air–water temperature signals. Additionally, ratios between daily water:air temperature signal means were lower in tributaries (∼0.68) when compared to main-stem (∼0.8) sites, increasing linearly throughout the observational period. In conceptual stream heat budget models, groundwater inflow had minimal effects on daily phase lags (∼0.2 hr), but increases in fractional groundwater discharge (0–50 %) depressed daily amplitude (∼20 % to 50 %) and mean ratios (∼15 %), supporting the sensitivity of daily metrics to interpreted changes in seasonal groundwater contributions to streamflow. During observational periods (i.e., April through October 2021 and 2022), significant differences (p < 0.01) between tributary and main-stem air–water metrics occurred when base-flow contributions were highest (∼0.93 vs. ∼ 0.68), as sites lower in the network had daily temperature metrics dominated by stream channel thermal inertia, rather than local groundwater connectivity, showing enhanced air–water diel signal coupling during warmer, drier periods. Divergent air temperature coupling across the network was interpreted as being driven by distance from local groundwater source zones, additional lateral groundwater inflows do not contribute a meaningful fraction to channel discharge lower in the network. Given the growing footprint of stream temperature observations, diel air–water temperature signals can provide distributed metrics sensitive to upstream groundwater discharge. Consequently, these metrics can support ongoing efforts by resource managers and researchers seeking to forecast the resilience of cold-water habitat to climate warming and changing precipitation regimes in mountain headwater streams.

New York, Virginia↗
USGS2024

Peak streamflow trends in South Dakota and their relation to changes in climate, water years 1921–2020

Peak-flow (flood) frequency analysis is essential to water-resources management applications, including the design of critical infrastructure such as bridges and culverts, and floodplain mapping. Federal guidelines for performing peak-flow flood frequency analyses are presented in a U.S. Geological Survey Techniques and Methods Report known as Bulletin 17C. A basic assumption within Bulletin 17C, which documents the guidelines for determining annual peak streamflow frequency, is that, for basins without major hydrologic alterations (for example, regulation, diversion, and urbanization), statistical properties of the distribution of annual peak streamflows are stationary; that is, the mean, variance, and skew are constant through time. Nonstationarity is a statistical property of a peak-flow series such that the long-term (on the order of decades) distributional properties change one or more times either gradually or abruptly through time. Individual nonstationarities may be attributed to one source such as flow regulation, land-use change, or climate but are often the result of a combination of sources, making detection and attribution of nonstationarities challenging. In response to a growing concern regarding nonstationarity in peak streamflows in the region, the U.S. Geological Survey, in cooperation with the Departments of Transportation of Illinois, Iowa, Michigan, Minnesota, Missouri, South Dakota, and Wisconsin; the Montana Department of Natural Resources and Conservation; and the North Dakota Department of Water Resources, assessed the potential nonstationarity in peak streamflows in the north-central United States. This chapter characterizes the effects of natural hydroclimatic shifts and potential climate change on annual peak streamflows in the State of South Dakota. Annual peak and daily streamflow as well as model-simulated gridded climatic data were examined for temporal monotonic trends, change points, and other statistical properties indicative of changing climatic and environmental conditions. Changes in annual peak and daily flows were evaluated among 13, 35, and 81 qualifying U.S. Geological Survey streamgages for the 75-, 50-, and 30-year trend periods through water year 2020 (the period from October 1, 2019, to September 30, 2020) in South Dakota, respectively. No qualifying streamgages were in the 100-year trend period in the State. Statistical tests for autocorrelation (independent and identically distributed assumption), monotonic trends, and change points in the median and scale are analyzed to evaluate potential stationarity violations (nonstationarity) for performing at-site peak-flow flood-frequency analysis. The trends are reported using a likelihood approach as an alternative to simply reporting significant trends with an arbitrary p -value cutoff point. A distinct east-west spatial pattern of likely upward and downward monotonic trends and change points, respectively, was detected in 75- and 50-year trend periods, but an inconsistent spatial pattern was detected in the 30-year trend period. Additionally, change points in the median annual peak streamflows were detected in the late 1970s and early 1980s in the western part of the State, but in the east, the change point was more commonly detected in 1992–93. A similar east-west spatial pattern of likely upward and downward trends was detected in the annual peak-flow timing, the day of the year of the annal peak streamflow. In the western part of the State, the annual peak streamflows are arriving earlier, but in the east, the annual peak streamflows are arriving later. A peaks-over-threshold (POT) analysis where, on average, there are two events per year (POT2) and four events per year (POT4) was also used to evaluate changes in the frequency (count) of daily streamflows exceeding the threshold. Similar to detected changes in the annual peak streamflow, an east-west likely upward or downward change corresponding to an increase or decrease, respectively, in the frequency of daily streamflow greater than a POT2 and POT4 threshold was detected. A monthly water-balance model was used to evaluate hydroclimatic variation in annual and seasonal precipitation, snowfall, potential evapotranspiration, and soil moisture storage for all qualifying streamgages in the 75-, 50-, and 30-year trend periods. Detected trends in the annual hydroclimatic metrics for the 75- and 50-year trend periods indicate a spatially consistent statewide increase in precipitation, decrease in snowfall, increase in potential evapotranspiration, and increase in soil moisture storage. Furthermore, detected trends in seasonal precipitation in the 75- and 50-year trend periods highlight a pronounced change in precipitation in winter and later into the summer season, especially in the 50-year trend period in the eastern part of the State. Statewide increases in seasonal soil moisture storage were also detected, highlighting year-round increasing flood magnitudes, particularly in the eastern part of the State. Based on the results of these stationarity tests for the qualifying streamgages in South Dakota among the 75-, 50-, and 30-year trend periods, consistent temporal and spatial patterns of nonstationarity were detected among the 75- and 50-year trend periods. Furthermore, when nonstationarity is detected in daily streamflow, increased streamflow and volume (increasing frequency in POT), as well as potentially bridge scour, may have implications on culvert and highway design in the eastern part of South Dakota. Thus, when performing at-site peak-flow flood-frequency analyses in South Dakota, potential nonstationarities and alternative approaches are important considerations.

South Dakota↗
USGS2024

Simulation of groundwater flow in the Long Island, New York regional aquifer system for pumping and recharge conditions from 1900 to 2019

The U.S. Geological Survey has developed a transient, groundwater-flow model that simulates hydrologic conditions in the Long Island aquifer system as part of an ongoing (since 2016) multiyear, cooperative investigation with the New York State Department of Environmental Conservation. The goals of this investigation are to assist stakeholders and resource managers to evaluate the response of the hydrologic system to changes in future hydraulic stresses. Responses in the hydrologic system include changes in water levels in the hydrogeologic units; discharge to streams, coastal waters, and subsurface infrastructure; and the extent of saline groundwater in the aquifers. Hydraulic stresses include future water-supply management and changes in land use and infrastructure. The numerical model synthesizes a diverse set of physiographic, geologic, climatic, land-use, and historical population, water use, and infrastructure data to physically represent the Long Island aquifer system from land surface to bedrock and to simulate annual hydrologic conditions between 1900 and 2019. A three-dimensional hydrogeologic framework was developed from existing and recently collected borehole geologic and geophysical data collected as part of a companion drilling program. Water-transmitting properties of the principal aquifer sediments were defined in three dimensions from new and existing lithologic logs. The distribution of recharge from precipitation was estimated from landscape characteristics and climate data. Anthropogenic recharge from wastewater, leaky infrastructure, and storm runoff were estimated from population, infrastructure, and pumping data. Water-use data, including well locations, depths, and pumping rates, were obtained from historical sources and records and used to estimate pumping stresses continuously in time and space, at an annual average time scale. The data were incorporated into a three-dimensional numerical model using the U.S. Geological Survey finite difference modeling code MODFLOW 6; the model encompassed all of Long Island and surrounding surface waters and simulated historical hydrologic conditions from 1900 to 2019. The calibration process involved trial and error adjustments using prior knowledge to improve general fit to observations followed by an inverse calibration to update and optimize input parameters, using an iterative ensemble smoother algorithm implemented in PEST++ version 5.0. This resulted in a model that generally was in good agreement with observed, dynamically varying hydrologic conditions from 1900 to 2019. The calibrated model was used to develop two base-case models for scenario testing of future, hypothetical conditions where one represented average-annual conditions, and one represented average-seasonal conditions from 2010 to 2019. The model representing average-annual conditions was modified further to represent an alternate sea-level position of 6 feet above the North American Vertical Datum of 1988, and the model representing average-seasonal conditions was modified to represent the average seasonal effects of a 5-year drought imposed upon current hydrologic conditions. Recharge is the sole source of water to the aquifer system; groundwater discharges to coastal water and streams and is withdrawn by pumped wells. Model-estimated annual recharge ranged from about 11 inches in 1965 to 41 inches in 1983. On average, from 2010 to 2019, about 23 percent of water was pumped from wells, and about 47 and 27 percent discharged to coastal waters and streams, respectively; the remaining 4 percent was water that moved into storage in the aquifer matrix. Water levels on Long Island vary naturally during time in response to changes in recharge; the amount of variation is largest in the interior of the island, in areas with highest water table altitudes near groundwater divides and lowest near streams and the coastal waters. The total range of water table altitudes on Long Island between 1900 and 2019 ranged from near 0 to more than 70 feet in western parts of Long Island. The largest range in altitudes is in New York City and is associated with areas of large historical withdrawals between the 1920s and the late 1980s. Water table altitudes generally varied by less than 10 feet in eastern Suffolk County, where the aquifer is under more natural conditions. Saltwater intrusion is of great concern on Long Island, particularly in western Long Island where both the unconfined and confined parts of the aquifer system have been intruded in response to large-scale groundwater withdrawals; however, the volume of freshwater in the islandwide aquifer system only has changed by about 5 percent between 1900 and 2019. The decadal change in the freshwater volume was largest during the early and mid-20th century, corresponding to the largest historical pumping, but that volume change did not exceed 1 percent. The negligible change in freshwater volume suggests that saltwater intrusion as of 2019 was limited at an islandwide scale but continues to occur in local areas of Queens and Nassau Counties, adversely affecting current water supplies and limiting future water supplies for affected communities. The regional groundwater model developed for this investigation is a tool that can be used to help determine the viability of current and future water supplies at a regional scale and can be used to support development of additional models at finer scale to support more focused assessments of groundwater sustainability.

New York↗
USGS2025

Reconnaissance of potential alternate water supply sources for the City of Gary, West Virginia

Seven potential sources of water, consisting of free-flowing discharge from abandoned coal mines at six locations and one abandoned flooded underground coal mine air shaft, were sampled for chemical analysis to assess the quality of the groundwater emanating from the seven mine sources. The six free-flowing mine discharge sources were also assessed for discharge by current-meter measurements on two separate occasions. The U.S. Geological Survey assessed these seven sources to provide information to the City of Gary, West Virginia (W. Va.), and the City of Gary’s consulting engineer with groundwater-quality and flow data to allow them to assess the seven sites as potential alternate sources of water for the City of Gary to augment its existing supply. For the six sites where discharge could be measured, discharge ranged from a minimum of 0.082 cubic feet per second (ft 3 /s) to a maximum of 3.685 ft 3 /s. Of the six sites measured, only two, Harmon Branch at Thorpe, W. Va. (USGS site 372201081303501) and the abandoned public-supply water wells near Havaco, W. Va. (USGS site 372358081344601), had discharge in excess of 1.00 ft 3 /s. Discharge from the abandoned public supply wells was 3.685 ft 3 /s on September 20, 2023, and 2.888 ft 3 /s on October 16, 2023, and discharge from Harmon Branch at Thorpe, W. Va., was 1.049 ft 3 /s on September 22, 2023, and 1.038 ft 3 /s on October 17, 2023. Discharge in the abandoned underground mine air shaft (USGS site 372224081340901) could not be assessed, but the air shaft drains an abandoned mine that likely contains water stored in approximately 1.7 square miles (mi 2 ) of abandoned underground coal mines in the Pocahontas No. 3 coal seam, and possibly an additional 0.9 mi 2 of leakage from the overlying Pocahontas No. 4 coal seam. Discharge for the six sites measured for the study was measured during a period between September 20 and October 18, 2023, and corresponded to the 12th to the 15th percentile of flow-duration statistics for the Tug Fork downstream of Elkhorn Creek at Welch, W. Va. streamgage (USGS site 03212750). Water-quality data for the seven sites sampled overall were acceptable with respect to drinking water standards. Of the 203 constituents analyzed, only a few failed to meet applicable U.S. Environmental Protection Agency (EPA) drinking water standards. Iron exceeded the 300 micrograms per liter (μg/L) secondary maximum contaminant level (SMCL) at only 1 of the 7 sites (14.3 percent) sampled. Iron concentrations ranged from a minimum of less than (<) 5.00 μg/L to a maximum of 724 μg/L with a median concentration of 7.62 μg/L. Manganese exceeded the 50.0 μg/L SMCL at 2 of the 7 sites (28.6 percent) sampled. Manganese concentrations ranged from a minimum of 1.93 μg/L to a maximum of 271 μg/L with a median concentration of 4.03 μg/L. No sites sampled exceeded the arsenic maximum contaminant level (MCL) of 10 μg/L. Arsenic concentrations ranged from a minimum of <0.100 μg/L to a maximum of 2.35 μg/L with a median arsenic concentration of 0.200 μg/L. None of the seven sites sampled for selenium for this study exceeded the EPA MCL of 50.0 μg/L. Selenium concentrations ranged from a minimum of <0.050 μg/L to a maximum of 5.26 μg/L with a median concentration of 3.21 μg/L. All seven sites were sampled for volatile organic compounds (VOCs), semivolatile organic compounds (SVOCs), and polychlorinated biphenyls (PCBs), but most had concentrations below the detection limit. Of the 10 PCB compounds analyzed for the seven sites sampled, none contained detectable concentrations of PCBs or Aroclor compounds. Of the 44 SVOCs analyzed at each of the seven sites sampled, only 1 SVOC, acenaphthene, was detected, at a concentration of 0.02 μg/L. Of the 96 VOCs analyzed, from each of the seven sites sampled, only two were found at detectable concentrations. Trichloromethane was detected only at 1 of the 7 (14.3 percent) sites sampled at a concentration of 0.027 μg/L, and benzene was detected at the same site and 3 additional sites (4 of the 7 sites or 57.1 percent of the sites sampled) at concentrations of 0.028, 0.029, 0.021, and 0.035 μg/L, but none exceeded the EPA MCL for benzene of 5.00 μg/L. Total coliform bacteria are ubiquitous in the environment, and their presence only suggests the potential for contamination by near-surface processes. Escherichia coli ( E. coli ) bacteria are derived from either human or animal fecal material and can be an indicator of potential contamination by pathogenic bacteria or viruses. Total coliform bacteria were detected at all 7 sites sampled at concentrations ranging from 17.5 to greater than (>) 2,420 most probable number per 100 mL (MPN/100 mL) of sample, with a median total coliform concentration of 1,553 MPN/100 mL. Escherichia coli bacteria were detected at 4 of the 7 sites sampled at concentrations ranging from <1 to 11.9 MPN/100 mL, with a median E. coli concentration of 5.1 MPN/100 mL.

West Virginia↗
USGS2000

Hydrology and water and sediment quality at James Campbell National Wildlife Refuge near Kahuku, Island of Oahu, Hawaii

The James Campbell National Wildlife Refuge occupies two lowland marsh and pond complexes on the northern coastal plain of Oahu: the mostly natural ponds and wetlands of the Punamano Unit and the constructed ponds of the Kii Unit. The U.S. Fish and Wildlife Service manages the Refuge primarily to protect and enhance habitat for four endangered species of Hawaiian waterbirds. Kii Unit is fed by artesian wells and rainfall, whereas Punamano Unit is fed naturally by rainfall, runoff, and ground-water seepage. Streams drain from the uplands into lowland ditches that pass through Kii Unit on their way to the ocean. A high-capacity pump transfers water from the inner ditch terminus at Kii to the ocean outlet channel. Stormwaters also exit the inner ditch system over flood-relief swales near the outlet pump and through a culvert with a one-way valve. A hydrologic investigation was done from November 1996 through February 1998 to identify and quantify principal inflows and outflows of water to and from the Refuge, identify hydraulic factors affecting flooding, document ground-water/surface-water interactions, determine the adequacy of the current freshwater supply, and determine water and sediment quality. These goals were accomplished by installing and operating a network of stream-gaging stations, meteorology stations, and shallow ground-water piezometers, by computing water budgets for the two Refuge units, and by sampling and analyzing water and pond-bottom sediments for major ions, trace metals, and organic compounds. Streamflow during the study was dominated by winter stormflows, followed by a gradual recession of flow into summer 1997, as water that had been stored in alluvial fans drained to lowland ditches. Outflow at the ditch terminus in 1997 was 125 million gallons greater than measured inflow to the coastal plain, mainly reflecting gains from ground water along the ditches between outlying gages and the ditch terminus. Of the measured 1997 outflow, 98 percent was through the Kii outlet pump, with the outlet culvert valve only opening for brief periods during storms. Large volumes of stormflow overflowed the flood-relief swales unmeasured. The largest storm of the study, in November 1996, was estimated to have a flood frequency of about 3 to 4 years. Streamflow exceeded culvert capacity and overtopped Kamehameha Highway at Kalaeokahipa Stream and Hospital ditch. Slight overbank flooding in Kii ditch resulted strictly from high discharge. Minor overbank flooding farther out on the coastal plain probably was caused mainly by the small hydraulic gradients available to convey stormflows along the lowland ditches. Stormwaters flooded Kii ponds and flowed back upstream along Punamano ditch into Punamano marsh, introducing suspended sediment and possibly other contaminants to the Refuge. Two smaller storms in January 1997 resulted in smaller flows and no overbank flooding. The Kii outlet pump ran continuously for 7 days during the November 1996 storm and for 1 to 2 days during the January 1997 storms. During all three storms, the outlet culvert valve opened and the inner ditches overtopped the flood-relief swales, allowing free outflow of water from the inner ditch. Backwater effects hindered drainage during the January 1997 storms at Hospital ditch at Kamehameha Highway, and at Punamano ditch at Nudist Camp Road (where the backflow into Punamano marsh in November 1996 constituted an extreme backwater effect). A probable marine backwater effect was imposed at the ocean outlet ditch during the November 1996 storm through a combination of high spring tides and wave setup from large surf. Whether this backwater effect propagated upstream in the ditches to affect inland sites could not be determined conclusively. A sand plug may have built up in the ocean outlet channel before the November 1996 storm, but if so, it probably washed out prior to, or early in the storm, and was not present at the time of peak stage at inland sites. A season-long buildup of the sand plug in late 1997 was inferred from rising water levels in the outlet ditch. Seawater flows up the outlet channel or over the sand berm and into the outer ditch system on most high tides, and particularly during spring high tides. Ponds and ditches of the Refuge and surrounding lowlands have mud- and clay-lined bottoms that form an effective confining unit and inhibit interaction with an underlying shallow limestone aquifer. At Kii Unit, pond levels are higher than adjacent ditch levels and underlying ground-water levels, establishing lateral and downward head gradients that could foster seepage losses from the ponds. Regional ground-water discharge from the Koolau aquifer to the coastal-plain sediments is mostly diffuse, but is concentrated where ridges of Koolau Basalt plunge beneath coastal-plain sediments near Punamano Unit and at the head of Hospital ditch. Kii ditch gains brackish ground water downstream of Kamehameha Highway. Wastewater disposal from the sewage treatment plant adjoining Kii Unit poses little or no threat to Refuge habitat. Disposal is at six injection wells located 0.45 mi away at Kahuku, and the wells inject into confined limestone aquifers that do not extend to Kii Unit. The natural freshwater supply to Punamano Unit is adequate for maintaining the wildlife habitat, judging from stable pond levels and low salinities there. A monthly water budget for Punamano showed an apparent annual deficit in measurable flows in 1997, requiring unmeasured ground-water gains equalling 51 inches of water. The freshwater supply to Kii Unit is inadequate according to Refuge managers, because there is not enough water to manipulate levels adequately in the ponds during most of the year, and particularly during the driest months. This is confirmed by monthly deficits in the water budget for the Kii ponds during summer months. However, the Kii budget showed an annual surplus in measurable flows for 1997 equalling 24 inches of water. Unmeasured losses are required to explain the apparent annual surplus, such as discharge to the ditches through pond water-control structures and downward and/or lateral ground-water seepage. The apparent surplus at Kii is strictly hydrologic and is not a surplus in a management sense; it cannot be stored or used to supply the Refuge, but instead reflects losses from the system that render this amount of water unavailable for use. The budget year, 1997, was drier than normal (24 percent below long-term mean rainfall) and so the measured potential evaporation for 1997 was probably higher than the long-term mean. Few metals or organic compounds of potential concern were detected in pond and ditch waters and in pond-bottom sediments. Detected pesticides were at trace levels or just above minimum reporting limits. Exceptions that exceeded quality guidelines for freshwater sediment were copper and zinc in sediment from Kii ponds C and D, and copper in sediment from Punamano north pond. Therefore, urban and agricultural runoff probably have contributed little in the way of harmful metals or organic compounds to the Refuge, although the potential for such contribution remains from periodic flooding of the ponds by ditch stormflows. Salinity was low throughout most Refuge waters, qualifying as fresh to slightly brackish and suitable for the environmental needs of Refuge fauna. Higher salinities have been observed in ditches during past periods of sugarcane cultivation and saltwater aquaculture, however. Resumption of saltwater aquaculture could raise ditch salinities if saltwater effluents are disposed directly into the ditches, as they were in the past.

Hawaii↗