USGS ScienceSearch

SEARCH · USGS Science

Results for “International Journal of Coal Geology”

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.

1,660 records · Page 57Linked to original sources

Sources of water and salts for the Zuni Salt Lake in west-central New Mexico

The Zuni Salt Lake is located in a maar in west-central New Mexico and contains hypersaline water that has long been used by Native Americans for religious purposes and the collection of salt. There have been several investigations suggesting different sources for the water and salt to the lake. Springs, seeps, and ephemeral streamflow have all been observed to contribute freshwater to the lake, and brackish to hypersaline seeps have been documented along the banks of the lake. This report summarizes the findings of a study that characterizes the lake’s hydrology, its water and salinity sources, and the hydrogeologic conceptual model. Regional groundwater levels indicate that each of the aquifers in the area have the potential to discharge groundwater to the lake. There is also evidence of vertical groundwater flow pathways at the maar that were likely created by the igneous intrusion that fractured the intersecting aquifers. A detailed water budget was constructed from continuous lake stage, precipitation, and evaporation data to estimate the groundwater inflow to the Zuni Salt Lake. It was determined that groundwater inflow to the lake is 441 ±94 acre-feet per year, which composes as much as 77 percent of the total inflows. The high sodium and chloride concentrations measured in two hypersaline samples collected near the lake indicate that the majority of the dissolved solids entering the lake are from a hypersaline groundwater source. The geochemical and isotopic compositions measured in the lake and surrounding features support the interpretation that hypersaline groundwater is the primary source of salts to the lake, which is likely sourced from the older (and deeper) Permian units. The hypersaline groundwater samples collected during this investigation have a unique aqueous chemistry relative to each of the mapped aquifers, and variability in groundwater compositions is interpreted to result from differences in minerology and residence time.

New Mexico

Spatial variation in effectiveness and costs of sagebrush restoration treatments across the western United States

Natural resource managers strive to improve restoration efficiency across the vast network of public lands of the United States, seeking the best return on investment through innovation, geographic prioritization, and adaptive management. This challenge is exemplified by ecosystems dominated by sagebrush ( Artemisia spp.), which occur across the western United States and have been degraded by increasing wildfire frequency driven by invasive plant spread, among other factors. However, both costs and effectiveness of sagebrush restoration can vary spatially due to management practices in addition to biotic and abiotic factors, and characterizing this variation could inform broad-scale planning. We leveraged published models of sagebrush recovery and treatment costs to predict and evaluate cost-effectiveness of aerial and drill seeding Artemisia spp. across 429 718 km 2 that overlapped 12 US states. Compared to natural recovery, effectiveness for both methods was generally low, yet effectiveness was often greater, and costs were always greater, for drill than aerial seeding, resulting in slightly higher cost-effectiveness on average for drill seeding. Cost-effectiveness for both treatments increased near major roads and in areas with more repeated burns. Cost-effectiveness also increased with predisturbance cover of sagebrush but decreased with herbaceous cover. However, we also identified areas where aerial seeding was more cost-effective than drill seeding, particularly with greater slope. Our results provide spatially explicit estimates and potential mechanisms of cost-effectiveness for two common seeding methods for sagebrush, which can help prioritize limited resources, guide land use, and improve restoration efficiency and effectiveness across public rangelands of the western United States.

western United States

A temporal look at the influence of topographic amplification on earthquake-triggered landslides in 3D seismic simulations

Earthquakes are a primary trigger for landslides, often leading to catastrophic consequences. While numerous studies have explored the spatial distribution of earthquake-triggered landslides, understanding the interaction between seismic waves and topography remains a critical challenge. Topographic irregularities can cause seismic wave amplification, altering ground shaking, and can trigger landslides that are challenging for predictive models to anticipate. This study investigates the spatial and temporal evolution of topographically amplified landslides, focusing on coseismic landslides triggered by the Mw 7.5 mainshock of the 2018 Papua New Guinea earthquake and post-seismic landslides associated with its four aftershocks, each exceeding Mw 6.0. We employ low-frequency, three-dimensional numerical ground shaking simulations and data-driven multivariate analyses to examine how landslides evolved from the coseismic to post-seismic periods. Our findings reveal a spatial shift in landslide distribution, in which the mainshock triggered slope failures predominantly on steep hillslopes, and the aftershocks triggered landslides on gentler slopes, often near geologic boundaries. We attribute this transition partly to the earthquake legacy effect of the mainshock, where the mainshock caused weakening of these hillslopes, making them more prone to failure when aftershocks occur. Additionally, the concentration of failures along geologic contacts in the post-seismic phase suggests that site amplification, stemming from contrasts in subsurface materials, exerts a key influence on landslide occurrence. Although not explicitly captured in our current numerical simulations, this mechanism warrants further investigation for more accurate hazard modeling.

JGR Solid Earth

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

Potential corrosivity of untreated groundwater in Louisiana

Corrosive groundwater can cause lead, copper, and other metals to leach from pipes and plumbing fixtures in water distribution systems. Metals, if ingested, could lead to serious health implications to the nearly 2.9 million people in Louisiana who obtain their drinking water from groundwater sources. Four indices—the Langelier Saturation Index (LSI), Ryznar Stability Index (RSI), Puckorius Scaling Index (PSI), and the Potential to Promote Galvanic Corrosion (PPGC)—in addition to an analysis which normalized the results from the existing indices, the Combined Index (CI), were used to assess the corrosivity of groundwater in Louisiana and identify areas within eight major aquifers and aquifer systems with moderate to high corrosivity potential. The purpose of this study is to provide State and local governments, public water system managers, and the nearly 500,000 private well owners in Louisiana with information needed to manage drinking-water supplies and mitigate potential health risks related to leaching of metals from water pipes and fixtures. The average scores of untreated groundwater samples from approximately 375 wells by index are as follows: LSI, −1.28; RSI, 9.78; PSI, 9.34; and CI, 4.14. The PPGC does not produce a numerical score, but the total percentage of class counts can be used to assign a classification; overall, samples in Louisiana were classified as significant concern. The percentages of groundwater samples from wells classified as potentially corrosive, by index, are as follows: LSI, 53 percent; RSI, 94 percent; PSI, 81 percent; PPGC, 98 percent; and CI, 81 percent. The percentages of samples classified as indeterminate, by index, are as follows: LSI, 46 percent; RSI, 5 percent; PSI, 12 percent; PPGC, 0 percent; and CI, 18 percent.

Louisiana

Transmissivity and water quality of water-producing zones in the intermediate aquifer system, Sarasota County, Florida

The intermediate aquifer system is an important water source in Sarasota County, Florida, because the quality of water in it is usually better than that in the underlying Upper Floridan aquifer. The intermediate aquifer system consists of a group of up to three water-producing zones separated by less-permeable units that restrict the vertical movement of ground water between zones. The diverse lithology, that makes up the intermediate aquifer system, reflects the variety of depositional environments that occurred during the late Oligocene and Miocene epochs. Slight changes in the depositional environment resulted in aquifer heterogeneity, creating both localized connection between water-producing zones and abrupt culmination of water-producing zones that are not well documented. Aquifer heterogeneity results in vertical and areal variability in hydraulic and water-quality properties. The uppermost water-producing zone is designated producing zone 1 but is not extensively used because of its limited production capability and limited areal extent. The second water-producing zone is designated producing zone 2, and most of the domestic- and irrigation-supply wells in the area are open to this zone. Additionally, producing zone 2 is utilized for public supply in southern coastal areas of Sarasota County. Producing zone 3 is the lowermost and most productive water-producing zone in the intermediate aquifer system. Public-supply well fields serving the cities of Sarasota and Venice, as well as the Plantation and Mabry Carlton Reserve well fields, utilize producing zone 3. Heads within the intermediate aquifer system generally increase with aquifer depth. However, localized head-gradient reversals occur in the study area, coinciding with sites of intense ground-water withdrawals. Heads in producing zones 1, 2, and 3 range from 1 to 23, 0.2 to 34, and 7 to 42 feet above sea level, respectively. Generally, an upward head gradient exists between producing zones 3 and 2. However, an upward head gradient between producing zones 2 and 1 does not consistently occur throughout Sarasota County, probably the result of greater ground-water withdrawals from producing zone 2 than from producing zone 1. The transmissivity of the intermediate aquifer system is spatially variable. Specific-capacity data from selected wells penetrating producing zones 2 and 3, were used to estimate transmissivity. Estimated transmissivity values for producing zones 2 and 3 range from about 100 to 26,000 feet squared per day and from about 1,300 to 6,200 feet squared per day, respectively. Because the capacity of specific water-producing zones is highly variable from site to site, estimating the performance of a specific water-producing zone as a water resource is difficult. Water samples collected during the study were analyzed for major-ion concentrations. Generally, bicarbonate type water from rock interaction occurs in northern Sarasota County; enriched calcium-magnesium-sulfate type water from deeper aquifers occurs in central Sarasota County; and sodium-chloride type water from saltwater mixing occurs in southern Sarasota County. In some areas of northern Sarasota County, the major-ion concentrations in water are lower in producing zone 2 than in producing zone 1. Major-ion concentrations in water are higher in producing zone 3 throughout the study area. A major objective of the study was to evaluate hydraulic and water-quality data to determine distinctions that could be used to characterize a particular producing zone. However, data indicate that both hydraulic and water-quality properties are highly variable within and between zones, and are more related to the degree of connection between and areal extent of water-producing zones than to aquifer depth and distance from the coast.

Florida

Pediment formation and subsequent erosion in Gale crater: Clues to the climate history of Mars

Evidence of paleo-rivers, fans, deltas, lakes, and channel networks across Mars has prompted much debate about what climate conditions would permit the formation of these surface water derived features. Pediments, gently sloping erosional surfaces of low relief developed in bedrock, have also been identified on Mars. On Earth, these erosional landforms, often thought to be created by overland flow and shallow channelized runoff, are typically capped by fluvial sediments, and thus in exceptionally arid regions, pediments are interpreted to record past wet periods. Here we document the Greenheugh pediment in Gale crater, exploiting the observational capability of the Curiosity rover. The absence of a fluvial cap suggests that the pediment was likely cut by wind erosion, not fluvial processes. The pediment was then buried by an aeolian deposit (Stimson sandstone) that mantled the lower footslopes of Aeolis Mons (informally known as Mt. Sharp). This burial terminated active wind erosion, preserving the pediment surface (as an angular unconformity). Groundwater was present prior-to, during, and shortly after Stimson deposition, perhaps contributing to lithification and certainly to early diagenesis. Post lithification, wind erosion first cut canyons in the northern most footslopes (north of Vera Rubin ridge). Unlithified gravels were deposited in these canyons, likely due to runoff from Mt. Sharp. Boulder-rich fluvial and debris flow deposits built a > 70 m thick sequence (Gediz Vallis ridge) on the southern Greenheugh pediment. Continued wind erosion left elevated patches of gravel on the northern footslopes, and exposure age dating shows that erosion essentially ceased before 1 Ga (but possibly much earlier). Erosion to the south led to emergence of Vera Rubin ridge, retreat of the Greenheugh pediment, and the formation of Glen Torridon valley. Hence, this footslope environment of Mt. Sharp records climate-driven periods of wind erosion, aeolian deposition (and groundwater activity), surface runoff and sediment deposition, followed by further significant wind erosion that declined to present very slow rates. This likely occurred during the late Hesperian and possibly into the Amazonian.

Icarus

Seismic moment and local magnitude scales in Ridgecrest, CA from the SCEC/USGS Community Stress Drop Validation Study

We illustrate the systematic difference between moment magnitude and local magnitude caused by underlying earthquake source physics, using seismic moments submitted to the Statewide California Earthquake Center/United States Geological Survey Community Stress Drop Validation Study 2019 Ridgecrest data set. While the relationship between seismic moment and moment magnitude ( M or M w ) of log 10 ( M 0 ) ~ 1.5* M is uniformly valid for all earthquake sizes by definition (Hanks and Kanamori, 1979), the relationship between local magnitude M L and moment is itself magnitude dependent. For moderate events, ~3< M < ~6, M and M L are coincident; for earthquakes smaller than ~3, M L ~ 1.0 log 10 M 0 (Hanks and Boore, 1984). This is a physical consequence of the corner frequency fc becoming larger than the upper frequency of observation and implies that M L and M differ systematically by a factor of 1.5 for these small events. While this idea is not new, we propose a new, continuous relationship between local magnitude and moment, for magnitudes 2 to 6 which extrapolates to smaller and larger magnitudes, applicable to southern California specific to the Ridgecrest region. We make use of the plethora of seismic moments as submitted by many participants of the Community Stress Drop study, compared to the Southern California Seismic Network (SCSN) catalog magnitudes. Overall, the seismic moments in the Community Study recover moment magnitude well, so we use our new M L - M 0 to convert M L to M , refining the SCSN operational M Lr scale. This systematic difference of 50% in slope between local and moment magnitude at small magnitudes has implications for spectral stress drop estimates, earthquake ground motion modeling, as well as other magnitude scales and earthquake occurrence statistics.

California

Helium detection as a guide for uranium exploration

Helium, a byproduct of radioactive decay, may prove to be a valuable indicator of the presence and distribution of uranium deposits. Recent technological advances permit the development of instrumentation not previously adapted for this purpose. Commercially available equipment can provide high sensitivity at low cost and allow reasonable mobility for field use. A truck-mounted mass-spectrometer, tuned for helium-4, permits immediate adjustment or modification of sampling patterns in response to accumulating data. The inlet system of the spectrometer has been designed to allow flexibility in gas analyses from various sample types--soil gas, atmosphere, or gases in water. Sensitivity of the instrument is better than 50 parts of helium per 10 9 parts of gas. Replicate samples and standards can be analyzed in only 3 minutes. Helium in soil gas is being used for the initial evaluation of the technique. A hollow probe, as much as 2 m long, is driven into the ground; and a 10-cc syringe is used to purge the probe and extract the gas. This sample is then injected into the inlet system and introduced into the spectrometer at constant pressure. The output signal is displayed on a chart recorder. Investigations can be performed on a qualitative, relative basis, or on a quantitative basis by comparison to calibrated helium standards. Preliminary field testing includes studies of the responses to variations in wind speed, temperature, barometric pressure, moisture, and sampling depth over extended time periods, as well as studies of geologic controls on the helium content in soil gas. Surveys over known uranium occurrences reveal some anomalous helium distributions. This report describes an updated technological approach to an old idea: that of using helium, a byproduct of uranium radioactive decay, as an exploration tool for uranium. Helium is the sixth most abundant gas found in the earth's atmosphere (Table 1). It is an inert gas and very mobile, mixing rapidly. When the alpha particles from the radioactive decay of uranium and thorium pick up two electrons, atoms of the isotope helium-4 are formed. Several helium-4 atoms are produced from each decay series of U-238, U-235, and Th-232 (fig. 1). The half-lives of the parent isotopes are5also shown in figure 1. One gram of uranium will produce ~10 5 atoms of helium-4 per second, and one gram of thorium produces ~2.5 x 10 4 atoms of helium-4 per second (fig. 2). Calculations of the helium produced from the crust and mantle reveal that 1,125 x 10 30 atoms per year are produced, but only 7 x 10 30 atoms per year escape from the earth (Damon and Kulp, 1958). More is being produced than is being lost; in fact, the total atmospheric content of helium could be produced in only 2 million years (MacDonald, 1963)! However, all the crustal and mantle helium-4 does not degas into the atmosphere as it is formed. It is trapped in crystal lattices and in pore spaces within the earth. What this means is that there is excess helium-4 in the earth--excess, that is, compared to the atmospheric concentration in equilibrium with the helium-4 escape rate into space and the helium-4 flux from the crust and mantle (Nicolet, 1957). Near-surface pockets of high helium-4 concentrations are known, and some natural gas fields are so enriched that they are the source for commercial production of helium. Work by Roberts and others (1975) has shown high helium concentration in the soil gas associated with geothermal areas; and work by Clark and Kugler (1973), Dyck (1975), and Goldak (1974), for example, has noted high helium concentrations in soils and waters in the vicinity of uranium deposits. These latter studies certainly indicate the potential of helium detection for use as an exploration tool for uranium.

Open-File Report

Towards a refined definition of “tidal influence” for the coastal wetland sciences

The astronomical tide relates to lunar (and to a lesser extent solar) influences on rising and falling water bodies, and is a key control of many coastal processes. However, tides manifest along coastal margins in different ways to muddle application beyond Sir Isaac Newton’s first descriptions. Here, we argue that tide, or specifically “tidal,” has not been adequately defined for the coastal wetland sciences to facilitate contemporary interdisciplinary distinction for rapidly changing coastlines and future climate adaptation planning. We provide a brief history of the tidal concept, and we describe tides in a diversity of uncommonly recognized marginal coastal environments. Common astronomical tide manifestations are described through hydrographic examples, then expanded to include those tidal environments exposed additionally to fluvial, meteorological (wind), or groundwater influences. Our expanded definition of “tidal” is intended to be a forward-looking construct as environments will inevitably change with sea-level rise, water extraction, and anthropogenic flow alteration. Our definition of wetland tidal influence, which we broaden, includes those “wetlands affected by astronomic water level fluctuations within coastally restricted environments causing frequent or infrequent surface inundation or groundwater tidal variation that affects the biogeochemistry of the soil (e.g., ion exchange capacity, oxygen state, mineralization processes)”. Astronomic water level fluctuations are often amplified or dampened by non-astronomical events, such as wind. Expanding the definition of what is considered tidal by the scientific community is critical for improved understanding of coastal wetland function and the suite of ecosystem services they provide people and society, including climate mitigation opportunity, management, and adaptation.

Wetlands

Computing flow-field distortion coefficients from well-construction and formation properties

Direct measurements of groundwater velocity made with borehole flowmeters in screened wells must be compensated for the effects of flow-field distortion (also known as borehole acceleration). A theoretical equation developed by Drost et al. (1968) and simple inputs describing hydraulic properties of well construction and geologic formation were programmed into an Excel workbook to facilitate computation by groundwater-flowmeter users. Tables describing the physical and hydraulic properties for well constructions and gravel pack media are provided with an example to facilitate use of the workbook. Groundwater flowlines converge or diverge as they pass from a geologic formation, through a gravel pack and well screen. The extent of flowline convergence or divergence and the value of the flow-field distortion coefficient is related to the relative changes in hydraulic conductivity of the well screen, gravel pack, and geologic formation. Convergence or divergence is accompanied by acceleration or deceleration of groundwater. Direct measurements of groundwater velocity at the center of the monitoring well can be adjusted to provide a more accurate estimate of velocity in the formation by applying a correction for flow-field distortion. Variables required to compute the flow-field distortion coefficient include the hydraulic conductivity of the gravel pack, well screen, and the geologic formation surrounding the well screen; the borehole radius, and the inside radius and outside radius of the well screen.

Groundwater

Base-flow sampling to enhance understanding of the groundwater flow component of nitrogen loading in small watersheds draining into Long Island Sound

Excessive nitrogen discharge is a major concern for the Long Island Sound. Programs have been implemented to reduce point sources of nitrogen to the sound, but little is known about the nonpoint sources. This study aims to better understand the current groundwater contributions of nitrogen from nonpoint sources in the Long Island Sound watershed. During the spring and summer of 2022, the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, collected water-quality samples to analyze nutrients (nitrogen and phosphorus), chloride, and bromide at 45 stations in the Long Island Sound watershed in Connecticut, New York, and Rhode Island. The stations were in small drainage watersheds (5 to 30 square kilometers) in the southern part of the Long Island Sound watershed. During two separate synoptic sampling events, water-quality samples and instantaneous streamflow measurements were collected under base-flow conditions (where the streamflow is dominated by groundwater inputs rather than overland flow or runoff flow). One sampling event was in the nongrowing season (April 24–25, 2022), and the other was in the growing season (June 30–July 1, 2022). To calculate instantaneous nitrogen loads and yields, streamflow was measured at the time of sample collection. Nitrogen concentrations, loads, and yields varied among sampling stations and by season. Total filtered nitrogen concentrations were generally lower in the nongrowing season (from less than 0.14 to 1.9 milligrams per liter) than in the growing season (from less than 0.23 to 3.0 milligrams per liter). Nitrate plus nitrite concentrations showed little variation between the nongrowing and growing seasons. Unfiltered ammonia plus organic nitrogen concentrations were generally lower in the nongrowing season (from less than 0.07 to 0.83 milligram per liter) than in the growing season (from 0.11 to 0.98 milligram per liter). In contrast, total filtered and unfiltered nitrogen loads and yields were higher in the nongrowing season than during the growing season, likely because streamflows were higher during the nongrowing season. Total unfiltered nitrogen yields during the nongrowing season ranged from less than 0.15 to 5.0 kilograms per square kilometer per day. Total unfiltered nitrogen yields during the growing season ranged from less than 0.12 to 2.5 kilograms per square kilometer per day. Total filtered nitrogen yields during the nongrowing season ranged from less than 0.13 to 5.2 kilograms per square kilometer per day. Total filtered nitrogen yields during the growing season ranged from less than 0.06 to 2.5 kilograms per square kilometer per day.

Connecticut, New York, Rhode Island

Migration of seismicity from the mantle to the upper crust beneath Harrat Lunayyir volcanic field, Saudi Arabia

Harrat Lunayyir is a volcanic field in Saudi Arabia that experienced a Mw~5.4 earthquake driven by an upper-crustal dike intrusion in May 2009. This volcanic field has exhibited numerous forms of volcanic seismicity both prior to and since the 2009 dike intrusion. Significantly, earthquakes within the lithospheric mantle and, rarely, the lower crust are present in the two-decade long seismicity catalog of Harrat Lunayyir. Here we analyze 24 years of volcanic seismicity at Harrat Lunayyir from 1998 to 2022. We find that: 1) precursory seismicity began at least eight years prior to the 2009 event, with a particularly notable seismic episode one year prior; 2) lithospheric mantle seismicity is highly localized in space and in time, largely occurring in discrete sequences lasting on the order of a few hours to a few days; 3) one seismic sequence clearly migrates upward from the lithospheric mantle to the upper crust, including seismicity within the nominally ductile lower crust; 4) crustal seismicity has been slowly declining over time; and 5) lithospheric-mantle seismicity does not show any apparent decline with time. From these observations we infer that the seismicity is driven by magmatic fluids or volatiles, and seismic monitoring of this volcanic field should continue into the future.

Harrat Lunayyir volcanic field

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

Shallow storage of the explosive Earthquake Flat Pyroclastics magma body, Okataina Volcanic Center, Taupo Volcanic Zone, New Zealand: Evidence from phase-equilibria experiments

Rhyolitic tuffs range widely in their crystal contents from nearly aphyric to crystal-rich, and their crystal cargoes inform concepts of upper crustal magma reservoirs. The Earthquake Flat pyroclastics (Okataina Volcanic Center, Taupo Volcanic Zone, New Zealand) are 10 km 3 of rhyolitic tuffs with abundant (~ 40 vol.%) plagioclase and quartz, minor biotite, hornblende, and orthopyroxene, and accessory Fe-Ti oxides, apatite, and zircon, set in high-silica rhyolitic glass. Major minerals form large, euhedral phenocrysts and abundant glomerocrysts with few disequilibrium textures excepting some faintly resorbed quartz. Plagioclase phenocrysts have thick rims of nearly constant composition near An 30 , and hornblende is weakly zoned or unzoned. The abundant and texturally complex mineral assemblage contrasts with the nearby (~ 25 km), nearly synchronous, but more voluminous and crystal-moderate rhyolite tuffs from Rotoiti caldera. New H 2 O-saturated phase-equilibria results on the erupted Earthquake Flat melt (glass) determine its co-saturation with the partial phenocryst assemblage of plagioclase, quartz, biotite, and Fe-Ti oxides at: 140 MPa, 755 ºC. These closely approximate the conditions of the pre-eruptive magma body assuming it was saturated with nearly pure H 2 O and at an f O 2 of ~ Ni–NiO. Absence of hornblende and orthopyroxene from the synthesized assemblages may result from those minerals being in a peritectic reaction relation with melt to produce biotite, so they would not grow from the liquid used as starting material. Experimental results on Rotoiti rhyolite (Nicholls et al. 1992) show that the two bodies resided at similar pressures, temperatures, and f O 2 s. Lower crystal abundance of the Rotoiti tuffs may result from slight compositional differences. We interpret that the Earthquake Flat pyroclastics were sourced from the crystal-rich periphery of a mushy reservoir system with the Rotoiti occupying a more melt-rich central location. Uncertain is whether this was a single intrusion zoned continuously in crystallinity, or discrete adjacent intrusions, but our results illustrate and quantify complexities of magma storage across relatively short distances.

North Island, Taupo Volcanic Zone

Integrity of production wells and confining unit at the Naval Weapons Industrial Reserve Plant, Dallas, Texas, 1995

Ground water in the shallow alluvial aquifer is contaminated at the Naval Weapons Industrial Reserve Plant, Dallas, Texas. Five production wells at the site are cased through the alluvial aquifer and underlying units and are screened in either the Paluxy or Twin Mountains aquifers. Three abandoned wells, originally completed in the Twin Mountains aquifer but filled with drilling mud in 1958, also penetrate the alluvial aquifer. The Paluxy and Twin Mountains aquifers are used for drinking-water supplies in and around the Dallas-Fort Worth metroplex. Trichloroethylene and its degradation products, dichloroethylene and vinyl chloride, and the metal chromium previously have been detected in the shallow alluvial aquifer. Current (1995) analyses of water-quality samples taken from the static water column of the five production wells and one of the abandoned wells indicate no trichloroethylene, dichloroethylene, or vinyl chloride in the water column of these wells. Chromium was detected in all samples, but concentrations were less than the practical quantitation limit, which is the regulatory action level for this site. The results of borehole geophysical log analysis indicate that two of the production wells could have vertically connected intervals where cement bonding in the well annulus is poor. The other production wells have overall good bonding. Temperature logs do not indicate flow behind casing except in the screened interval of one well. Geophysical logs show the Eagle Ford Shale ranges from 147 to 185 feet thick at the site. The Eagle Ford Shale has low permeability and a high plasticity index. These physical characteristics make the Eagle Ford Shale an excellent confining unit.

Texas

Environmental monitoring of groundwater, surface water, and soil at the Ammonium Perchlorate Rocket Motor Destruction Facility at the Letterkenny Army Depot, Chambersburg, Pennsylvania, 2021

Letterkenny Army Depot in Chambersburg, Pennsylvania, built an Ammonium Perchlorate Rocket Motor Destruction (ARMD) Facility in 2016 to centralize rocket motor destruction and contain all waste during the destruction process. The U.S. Geological Survey has collected environmental samples from groundwater, surface water, and soils at ARMD since 2016. During 2021, samples were collected from four groundwater wells in September, one surface-water site in October, and five soil sites in November near the facility. Samples were analyzed for nutrients, trace metals, major ions, total volatile organic compounds, and perchlorate. Perchlorate was not detected in any 2021 samples. Groundwater results showed no constituents exceeded any U.S. Environmental Protection Agency (EPA) maximum contaminant level (MCL). Dissolved arsenic (As) was detected in one well above the reporting detection level (RDL) of 3 micrograms per liter (μg/L) at 5.4 μg/L but below its MCL of 10 μg/L. Dissolved iron (Fe) was the only inorganic constituent measured above an EPA secondary maximum contaminant level (SMCL). All groundwater samples collected in 2021 exceeded the Fe SMCL of 300 μg/L, with concentrations ranging from 390 μg/L to 3,500 μg/L. Surface-water data collected during 2021 showed no measured constituents in the surface-water sample that exceeded any EPA MCL or SMCL. Soil samples collected from 2016 through 2021 showed all concentrations of As exceeded the EPA soil screening levels of 3 milligrams per kilogram (mg/kg) but did not exceed the Pennsylvania medium-specific concentrations for As of 61 mg/kg. Arsenic concentrations in 2021 ranged from 9.1 mg/kg to 12.9 mg/kg. The 2021 results for the ARMD Facility indicate no increases in concentrations of reported compounds compared to data from 2016 to 2020. The contained burn treatment facility for demilitarization of rocket motors during 2021 appears to have operated without elevating concentrations of target compounds compared to previous years.

Pennsylvania

Characterization of selected biological, chemical, and physical conditions at fixed sites in the Upper Colorado River Basin, Colorado, 1995-98

Biological community samples were collected at 15 sites in the Upper Colorado River Basin (UCOL) in Colorado as part of the National Water-Quality Assessment (NAWQA) Program. Sites sampled in two physiographic provinces, the Southern Rocky Mountains and the Colorado Plateau, represented agriculture, mining, urban and recreation, and mixed land uses and background conditions. Nine measures of water quality, which include information on nutrients, specific conductance (a surrogate for salinity), trace elements in streambed sediment, pesticides in fish tissue, fish communities, and macroinvertebrate richness and composition and stream habitat were used for comparisons among sites within the two physiographic provinces. Sampling sites from three other NAWQA study units—the Rio Grande Valley, the South Platte River Basin, and the Upper Snake River Basin study units—were categorized on the basis of land use and stream size in order to develop a larger data set for comparison to sites in the UCOL. Three categories of land use—forested (includes mining, urban and recreation, and background), agriculture, and mixed—were used for comparison to the UCOL fixed sites. Results indicated that all sites other than the Colorado River below Baker Gulch (a background site) showed some water-quality characteristics to be significantly affected. Results indicated that the concentrations of cadmium and zinc in streambed sediment at mining land-use sites in the Southern Rocky Mountains physiographic province generally were orders of magnitude higher than streambed-sediment concentrations at the background site. Streambed-sediment concentrations at mining land-use sites in the UCOL were greater than the 75 th percentile of concentrations from sites in the three other NAWQA study units. Fish communities and habitat conditions were degraded at mining land-use sites compared to the background site. Ephemeroptera, Plecoptera, and Trichoptera (EPT) richness and the percentage of EPT were lower at mining land-use sites than at the background site and were less than the 50 th percentile of those for sites from the three other NAWQA study units. Nutrient concentrations at urban and recreation sites in the Southern Rocky Mountain physiographic province generally were greater than concentrations at the background site and generally were between the 25 th and 90 th percentile of concentrations for sites from the three other NAWQA study units. Habitat conditions and fish communities at urban and recreation sites were slightly degraded compared to the background site. EPT richness and the percentage of EPT were lower at urban and recreation sites than at the background site and were between the 25 th and 75 th percentile of those for sites from the three other NAWQA study units. The percentage of Chironomidae, which may be indicative of pollutant-tolerant organisms, was higher at urban and recreation sites than at the background site. Mixed land-use sites in the Southern Rocky Mountains physiographic province had similar nutrient concentrations and similar cadmium and zinc streambed-sediment concentrations. Fish-community degradation index values were very different among the three mixed land-use sites in the Southern Rocky Mountains physiographic province. Larger percentages of omnivores and anomalies such as lesions and deformities at two mixed land-use sites resulted in higher degradation values of the fish community. Agriculture land-use sites had higher concentrations of nutrients and selenium than the background site in the Colorado Plateau physiographic province. Concentrations of p,p '-DDE in fish tissue at agriculture sites were higher than the 75 th percentile of concentrations for sites from the three other NAWQA study units. Fish communities had degradation values near the 75 th percentile for agriculture sites. The percentage of EPT was low at agriculture sites when compared to the background site. Two mixed land-use sites in the Colorado Plateau physiographic province had similar concentrations of nutrients, selenium, and p,p '-DDE, and similar EPT richness and composition. These two sites were located downstream from agricultural and urban activities. Some water-quality measures at these two sites indicated degradation compared to a mixed land-use site upstream from most of the agriculture and urban activities in the Colorado Plateau physiographic province.

Colorado