Geology ReportsSearch

Geology topics

Browse USGS reports

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

1,748 records · Page 71Linked to original sources

The 17 January 1994 Northridge, California, earthquake: A retrospective analysis

The 17 January 1994 Northridge, California, earthquake was a watershed event, with far-reaching societal and scientific impacts. The earthquake, which occurred in the early days of both broadband seismic networks and the Internet, spurred advances in seismic monitoring, real-time systems, and development of data products. Motivated by the 30th anniversary of the earthquake, we present a brief retrospective of the earthquake and its impact, and reconsider both ground motions and the aftershock distribution using modern tools and best-available data. With improvements in instrumentation and analysis methodology, recent earthquakes continue to reveal increasing complexity of ground motions, fault systems, and earthquake ruptures. Even in the absence of data from state-of-the art instrumentation, a retrospective consideration of ground motion data from the Northridge earthquake reveals complexities beyond what could be characterized (and modeled) thirty years ago. Aftershock relocations for both the 1971 Sylmar and 1994 Northridge earthquakes also reveal an updated view of fault complexity. Our study does provide a cautionary tale regarding legacy data sets and research results that are not easily accessible, which can result in discrepancies between catalog data and products from best-available science. We also briefly describe outreach products produced as part of the anniversary commemoration.

California

Effects of stochastically-simulated near-fault ground motions on soil liquefaction

The scarcity of historically recorded near-fault ground motions poses a challenge to systematically understanding the influence of near-fault effects on various types of seismic demands for engineering purposes. In particular, the current state of knowledge of the influence of ground-shaking intensity on soil liquefaction and its consequences does not specifically account for the effects of near-fault ground motion characteristics. In this study, the influence of near-fault ground motions on liquefaction triggering and lateral spreading are investigated using non-linear modeling of a hypothetical liquefiable soil column in the finite-element computational platform OpenSees subjected to simulated ground motion time series that represent strong earthquake shaking in the near field. The simulated ground motion time series and resulting datasets are based on a parametric stochastic model and are developed for a range of source and path parameters to represent a realistic variability of ground motion characteristics. Dependencies between ground motion intensity measures (IMs) and liquefaction demand parameters are investigated for near-fault pulse and nonpulse-like ground motion sets. Evolutionary IMs, such as cumulative absolute velocity (CAV) and the time-varying magnitude-adjusted peak ground acceleration (PGAM), are considered in developing liquefaction triggering probability density functions. Post-liquefaction triggering responses such as lateral spreading displacements are examined in relation to PGA M and CAV. The ground motion simulations are validated by comparing their liquefaction-capacity PGA M fragilities and post-triggering CAV vulnerability relationships to historical records from the 1994 Northridge earthquake in California, USA. Finally, a path forward for future studies that includes finding systematic differences in the IM-liquefaction demand relationships between near-fault and far-field stochastic ground motion sets is outlined.

Conference Paper

Geological context and significance of the clay-sulfate transition region in Mount Sharp, Gale crater, Mars: An integrated assessment based on orbiter and rover data

On Mars, phyllosilicate (“clay”) minerals are often associated with older terrains, and sulfate minerals are associated with younger terrains, and this dichotomy is taken as evidence that Mars’ surface dried up over time. Therefore, in situ investigation of the Mount Sharp strata in Gale crater, which record a shift from dominantly clay-bearing to sulfate-bearing minerals, as seen in visible−near-infrared orbital reflectance spectra, is a key science objective for the Mars Science Laboratory (MSL) Curiosity rover mission. Here, we present regional (orbiter-based) and in situ (rover-based) evidence for a low-angle erosional unconformity that separates the lacustrine and marginal lacustrine deposits of the Carolyn Shoemaker formation from the dominantly eolian deposits of the lower Mirador formation within the orbitally defined clay-sulfate transition region. The up-section record of wetter (Carolyn Shoemaker formation) to drier (lower Mirador formation) depositional conditions is accompanied by distinct changes in diagenesis. Clay minerals occur preferentially within the Carolyn Shoemaker formation and are absent within the lower members of the Mirador formation. At and above the proposed unconformity, strata are characterized by an increase in diagenetic nodules enriched in X-ray amorphous Mg-sulfate. Early clay formation in the Carolyn Shoemaker formation may have created a hydraulic barrier such that later migrating magnesium- and sulfur-rich fluids accumulated preferentially within the lower members of the Mirador formation. The proposed unconformity may have also acted as a fluid conduit to further promote Mg-sulfate nodule formation at the Carolyn Shoemaker−Mirador formation boundary. These results confirm an association of the clay-sulfate transition with the drying of depositional environments, but they also suggest that at least some orbital sulfate signatures within the region are not time-congruent with the environmental signals extracted from primary sedimentology. Our findings highlight that complex interactions among primary depositional environment, erosion, and diagenesis contribute to the transition in clay-sulfate orbital signatures observed in the stratigraphy of Mount Sharp.

GSA Bulletin

Representation of surface-water flows using Gradient-Related Discharge in an Everglades Network

The Everglades Depth Estimation Network interpolates water-level gage data to produce daily water-level elevations for the Everglades in south Florida. These elevations were used to estimate flow vectors (gradients and directions) and volumetric flow rates using the Gradient-Related Discharge in an Everglades Network (GARDEN) application developed by the U.S. Geological Survey in cooperation with the U.S. Army Corps of Engineers. Flow rates in both the east-west and north-south directions were computed on a 400-meter square grid using modified parameters in the Manning’s equation. The frictional resistance parameter in the Manning’s equation was calibrated to measured flow rates at coastal creeks fed by Everglades Depth Estimation Network boundary flows. Levees and other features that act as barriers to flow were defined as “no-flow” grid cells where vectors were set to zero. The flow volume magnitudes were calibrated with 2020 daily values of coastal river flows, and verification was performed using 2021 data. Within a given day, the measured coastal river flows fluctuate more than the GARDEN boundary flows because of tidal and wind forcings. Because the GARDEN boundary flows were the upstream water source for the coastal rivers, calibration focused on matching average daily flow volumes rather than daily fluctuations. The Pearson’s correlation coefficient is 0.766 for the 2020 calibration period and 0.566 for the 2021 verification period. Applying GARDEN to periods with hydraulic-control-structure releases allows the propagation of structure flows to be seen in the daily flow-vector maps along with the multiday response of flows farther downgradient. Flow vectors may be overestimated near control structures because of difficulties in resolving the water gradient downstream from the structure. Flow vectors farther from the structure are more accurate than those near the structure.

Florida

Prediction of regional broadband strong ground motions using a teleseismic source model of the 18 April 2014 Mw 7.3 Papanoa, Mexico, earthquake

To estimate predicted ground motion from a teleseismic slip model, we use a low‐ and high‐frequency hybrid method to simulate the regional, strong ground motions observed following the 18 April 2014 moment magnitude ( ⁠M w ⁠ ) 7.3 Papanoa, Mexico, earthquake. To generate the regional ground motion at low frequencies (<1 Hz), a teleseismically derived, finite‐fault, kinematic model is used to define the earthquake source, taking into account slip‐model variations identified with a parameter sampling approach that considers possible errors in the fault geometry, the hypocenter depth, and the rupture velocity. A 3D crustal model is used to calculate the low‐frequency ground motions using a finite‐element calculation that includes topography and considers variations in the source model to estimate the uncertainty in the calculations. High frequencies (>1 Hz) are added using a 1D full‐wave propagation code that estimates uncertainties by considering multiple random distributions of slip with different spatial correlation lengths. The synthetic, broadband (0.05–10.0 Hz) ground motions are obtained by combining the low‐ and high‐frequency portions match filtered at 1 Hz. These synthetic ground motions are compared with the regional observations using velocity records, peak ground acceleration, and medians of the orientation‐independent response spectra of the horizontal components (RotD50) calculated at periods of 0.2, 0.3, 0.5, 1.0, 2.0, 3.0, 5.0, 7.5, and 10.0 s. The results indicate that ground motions estimated at these periods using our hybrid approach based primarily on a teleseismically derived source model are comparable to the values observed for the 2014 Papanoa earthquake at regional distances. The approach could be used to estimate strong‐motion spectral levels expected for regions with limited local and regional recordings and could also fill in magnitude or distance gaps in ground‐motion prediction relations utilized in the assessment of seismic hazard.

Papanoa

Projected sea-level rise and high tide flooding at De Soto National Memorial, Florida

Introduction National parks and preserves in the South Atlantic-Gulf Region contain valuable coastal habitats such as tidal wetlands and mangrove forests, as well as irreplaceable historic buildings and archeological sites located in low-lying areas. These natural and cultural resources are vulnerable to accelerated sea-level rise and escalating high tide flooding events. Through a Natural Resources Preservation Program-funded project during 2021–23, the U.S. Geological Survey, in collaboration with the National Park Service, estimated the probability of inundation at De Soto National Memorial, Florida, and several other parks under various sea-level rise scenarios and contemporary high tide flooding thresholds. The maps produced for this effort can be used to assess potential habitat change and explore how infrastructure and cultural resources within the park may be exposed to future flooding-related hazards.

Florida

Dissolved arsenic concentrations in surface waters within the upper portions of the Klamath River Basin, Oregon and California

Arsenic toxicity is an environmental health problem. Levels of arsenic in surface waters at some locations in the Klamath River Basin in southern Oregon and northern California can exceed the U.S. Environmental Protection Agency (EPA) standard for drinking water. There are both anthropogenic and natural sources of arsenic. The Klamath River Basin consists primarily of volcanic deposits and contains an underground geothermal system with hot springs and warm water wells, all known natural sources of arsenic. Anthropogenic sources of arsenic are related to the agricultural use of herbicides, fungicides, and insecticides. Surface water arsenic levels can also be affected by fertilizer amendments, evaporative concentration, oxygen-level depletion, and various geochemical transformations that can increase arsenic mobilization. In this study by the U.S. Geological Survey and the Bureau of Reclamation, dissolved concentrations of arsenic, copper, and lead were measured in surface waters at 39 unique sites within the upper portions of the Klamath River Basin between 2018 and 2022. In every year, except 2022, sites were sampled four times between April and November. Surface-water arsenic concentrations varied up to four-orders of magnitude among sites. Median arsenic concentration was lowest at Cherry Creek (0.03 micrograms per liter [μg/L]) and highest at Wood Kimball Spring (36.7 μg/L), two sites located north of Upper Klamath Lake. The highest arsenic concentrations (17.4±4.9 μg/L, n =3) were found in drain sites (defined here as a waterbody returning used irrigation water) while the lowest arsenic concentrations were found in an artesian well (0.8 μg/L, n =1). The elevated arsenic concentrations of the drain sites suggest that arsenic might be concentrated or mobilized by agricultural activities, water re-use practices, and (or) by geochemical processes occurring around water stored in drains (that is, in the water column and across sediment water boundaries). A source of arsenic in drain water in the Klamath Strait Drain area includes water used for irrigation originating from Ady Canal. Other potential sources include groundwater, geothermal water, and local soils and sediments. Seasonal differences in surface-water arsenic concentrations were detected at 13 sites, 10 of which had higher arsenic concentrations in summer than in either spring or fall. The sites sampled around Upper Klamath Lake, the impounded rivers, one of the two canal sites, and 5 of the 14 river sites had higher surface-water arsenic concentrations in the summer than in either spring or fall. Surface-water arsenic concentrations from groundwater sources (that is, springs and in the artesian well) did not vary significantly among seasons (p-values greater than 0.1). Median surface-water concentrations of copper and lead ranged from 0.03 to 3.7 μg/L, and from 0.013 to 0.175 μg/L ( n =2–18), respectively. Dissolved concentrations of both metals were below acute toxicity endpoints reported by the EPA for freshwater animals. Surface-water arsenic concentrations varied independently from corresponding changes in surface-water lead or copper concentrations. However, arsenic concentrations measured in bed-sediment samples collected from a subset of sites located north of Upper Klamath Lake correlated strongly and significantly with the corresponding sedimentary lead concentrations ( p =0.015). Aqueous arsenic speciation measured in a subset of sites in 2019 and 2022 showed that all the arsenic existed as arsenic (V), the most oxidized arsenic species, and presumably, the least toxic. The highest proportions of arsenite (As(III)), the presumably most toxic arsenic species, relative to total arsenic concentrations were found at drain sites. Our assessment of dissolved arsenic concentrations in various surface-water bodies in the Upper Klamath River Basin reveals geographical areas of consistently low (below 2.1 μg/L), moderate (below 10 μg/L) and high (above 10 μg/L) surface-water arsenic concentrations. South of Upper Klamath Lake, surface-water arsenic concentrations were consistently higher than 20 μg/L at two drain sites located in an area of predominant agricultural land use with extensive water re-use practices. North of Upper Klamath Lake, surface-water arsenic concentrations greater than 20 μg/L were consistently measured at sites with limited nearby agricultural activities, suggesting a geogenic source. The consistently high arsenic levels from the Wood River at Jackson F. Kimball State Park, Fort Creek, and Crooked Creek, which are sites located at or near headwater spring sources, suggest a natural background source of arsenic. Water flowing downstream from this area could be a potential source of arsenic to Upper Klamath Lake and the Upper Klamath River.

California, Oregon

Turbidite correlation for paleoseismology

Marine turbidite paleoseismology relies on the assumption of synchronous triggering of turbidity currents by earthquake shaking to infer rupture extent and recurrence. Such inference commonly depends on age dating and correlation of the physical stratigraphy of deposits carried by turbidity currents (i.e., turbidites) across great distances. Along the Cascadia subduction zone, which lies offshore the Pacific Northwest, USA, turbidite facies in core photographs, X-ray computed tomography images, and magnetic susceptibility (MS) data exhibit differences in character over relatively short distances, which implies that not all deposits can be correlated with confidence. Thus, subjective correlation based on expected similarity over great distances and weak age constraints does not independently support paleoseismic models. We present a new method for correlating turbidites along the Cascadia margin that can yield a more objective and repeatable stratigraphic framework to underpin earthquake recurrence. We use dynamic time warping to correlate MS logs and measure correlation coefficients of core pairs to evaluate correlation strength. We then compare these measures to a distribution of correlation coefficients of randomly generated turbidite sequences and find that only a small number of core pairs can be correlated more confidently than randomly stacked turbidites. This methodology promises a more robust correlation strategy for future stratigraphic studies.

Oregon, Washington

Relatively stable pressure effects and time-increasing thermal contraction control Heber geothermal field deformation

Due to geological complexities and observational gaps, it is challenging to identify the governing physical processes of geothermal field deformation including ground subsidence and earthquakes. In the west and east regions of the Heber Geothermal Field (HGF), decade-long subsidence was occurring despite injection of heat-depleted brines, along with transient reversals between uplift and subsidence. These observed phenomena contradict current knowledge that injection leads to surface uplift. Here we show that high-yield production wells at the HGF center siphon fluid from surrounding regions, which can cause subsidence at low-rate injection locations. Moreover, the thermal contraction effect by cooling increases with time and eventually overwhelms the pressure effects of pressure fluctuation and poroelastic responses, which keep relatively stable during geothermal operations. The observed subsidence anomalies result from the siphoning effect and thermal contraction. We further demonstrate that thermal contraction dominates long-term trends of surface displacement and seismicity growth, while pressure effects drive near-instantaneous changes.

California

The U.S. Geological Survey Ohio Water Microbiology Laboratory

The U.S. Geological Survey Ohio Water Microbiology Laboratory is a part of the Ohio-Kentucky-Indiana Water Science Center. The mission of the laboratory is to provide microbiological data of public health significance from surface waters, groundwaters, and sediments for a variety of study objectives. The laboratory conducts internal projects, works with external cooperators, and assists U.S. Geological Survey offices and National programs. The laboratory offers guidance, study design, and data interpretation expertise to collaborators, all following rigorous quality control and quality assurance procedures.

Fact Sheet

Groundwater, surface-water, and water-chemistry data, Black Mesa area, northeastern Arizona—2019–2021

The Navajo (N) aquifer is an extensive aquifer and the primary source of groundwater in the 5,400-square-mile Black Mesa area in northeastern Arizona. Water availability is an important issue in the Black Mesa area because of the arid climate, past industrial water use, and continued water requirements for municipal use by a growing population. Precipitation in the area typically ranges from less than 6 to more than 16 inches per year, depending on location. The U.S. Geological Survey water-monitoring program in the Black Mesa area began in 1971 and provides information about the long-term effects of groundwater withdrawals from the N aquifer for industrial and municipal uses. This report presents the results of data collected as part of the monitoring program in the Black Mesa area from calendar years 2020–2021 and, additionally, uses streamflow statistics from November and December 2019. The monitoring program includes measurements of (1) groundwater withdrawals (pumping), (2) groundwater levels, (3) spring discharge, (4) surface-water discharge, and (5) groundwater chemistry. In calendar year 2020, total groundwater withdrawals were estimated to be 2,680 acre-feet (acre-ft), and, in 2021, total withdrawals were estimated to be 2,570 acre-ft. Total withdrawals during 2021 were about 65 percent less than total withdrawals in 2005 because the Peabody Western Coal Company discontinued its use of water to transport coal in a coal slurry pipeline after 2005 and ceased mining operations in 2019. Owing to Navajo Nation and Hopi Reservation access restrictions during the Coronavirus pandemic, water levels were not collected from municipal wells in 2020 or 2021. Water levels measured in 2021 from wells completed in the unconfined areas of the N aquifer within the Black Mesa area showed a decline in 7 of 13 wells when compared with water levels from the prestress period (prior to 1965). The changes in water levels across all 13 wells ranged from +8.4 feet (ft) to −42.4 ft, and the median change was −0.4 ft. Water levels also showed decline in 11 of 12 wells measured in the confined area of the aquifer when compared to the prestress period. The median change for the confined area of the aquifer was −25.9 ft, with changes across all 12 wells ranging from +17.3 ft to −133.7 ft. Spring flow was measured at four springs between 2020 and 2021. Flow fluctuated during the period of record for Burro Spring and Pasture Canyon Spring, but a decreasing trend was statistically significant (p<0.05) at Moenkopi School Spring and Unnamed Spring near Dennehotso, Arizona. Discharge at Burro Spring has remained relatively constant since it was first measured in the 1980s, and discharge at Pasture Canyon Spring has fluctuated for the period of record. Continuous records of surface-water discharge in the Black Mesa area were collected from streamflow-gaging stations at the following sites: Moenkopi Wash at Moenkopi 09401260 (1976–2021), Dinnebito Wash near Sand Springs 09401110 (1993–2020), Polacca Wash near Second Mesa 09400568 (1994–2020), and Pasture Canyon Springs 09401265 (2004–2021). Median winter flows (November through February) of each winter were used as an estimate of the amount of groundwater discharge at the above-named sites. For the period of record, the median winter flows have generally remained constant at Polacca Wash and Pasture Canyon Springs, whereas a decreasing trend was observed at Moenkopi Wash and Dinnebito Wash. In 2020 and 2021, water samples were collected from a total of four springs in the Black Mesa area and analyzed for selected chemical constituents. Results from the four springs were compared with previous analyses from the same springs. Dissolved solids, chloride, and sulfate concentrations increased at Moenkopi School Spring during the more than 30 years of record at that site. Concentrations of dissolved solids and sulfate at Pasture Canyon Spring have not varied significantly (p>0.05) since the early 1980s, and there is no increasing or decreasing trend in those data. However, concentrations of chloride from Pasture Canyon Spring show a diminishing trend. Concentrations of dissolved solids, chloride, and sulfate at Unnamed Spring near Dennehotso have varied for the period of record, but there is no statistical trend in the data. Concentrations of dissolved solids at Burro Spring have varied for the period of record, but there is no statistical trend in the data. However, concentrations of chloride and sulfate from Burro Spring show a trend towards lower concentrations.

Arizona

Bedrock geologic map of the Woodstock quadrangle, Grafton County, New Hampshire

The bedrock geology of the Woodstock 7.5-minute quadrangle consists of highly deformed metasedimentary rocks of the Central Maine trough, including the Silurian Rangeley and Perry Mountain Formations and the Devonian Littleton Formation. The central, northern, and eastern parts of the quadrangle are underlain by the oldest rocks in the area, the Rangeley Formation. In the southwest and south-central part of the quadrangle, metaturbidites of the Perry Mountain Formation and subsequent Littleton Formation overly the Rangeley Formation in a deformed F 1 synform, herein informally called the Bagley Brook basin. The metasedimentary rocks were intruded by widespread syn- to post-tectonic granitoids of the Devonian New Hampshire Plutonic Suite and minor post-metamorphic Jurassic-Cretaceous mafic dikes of the White Mountain Plutonic-Volcanic Suite. The metasedimentary rocks were affected by at least two episodes of deformation in the Devonian Acadian orogeny. The dominant regional foliation is second-generation (S 2 /D 2 ) and formed during the development of sillimanite-muscovite mineral assemblages. Large bodies of the Early Devonian Kinsman Granodiorite intruded the metasedimentary rocks semi-concordantly during D 2 deformation. Dikes of the Late Devonian Concord Granite cut the Kinsman Granodiorite and the metasedimentary rocks and were emplaced either syn- or post-D 2 . The map pattern in the Rangeley Formation is dominated by northeast to northwest trending, moderately to steeply north-dipping F 2 and F 3 folds. Map-scale F 1 folds are defined by the Bagley Brook basin. Previous division of Rangeley Formation stratigraphy in this region into “upper” and “lower” parts was not corroborated by 1:24,000-scale mapping of lithodemic units, and rocks previously mapped as part of the Smalls Falls and Madrid Formations are here reassigned to the Rangeley Formation. Some rocks previously mapped as the lower part of the Littleton Formation are now assigned to the Perry Mountain Formation. The Littleton Formation on this map is approximately equivalent to rocks previously mapped as the upper part of the same formation. Steeply dipping fractures in the quadrangle show a preferred northeast orientation, consistent with subsurface fracture orientations in the well fields near Mirror Lake. Jurassic-Cretaceous mafic dikes and normal faults show preferred northeast orientations, similar to the fractures, suggesting that the extensional stress field that controlled dike orientation during the Mesozoic also produced the dominant brittle fabrics in the area.

New Hampshire

Examining 22 years of ambient seismic wavefield at Mount St. Helens

An increase in seismic activity precedes most volcanic eruptions. Whereas event-based forecasting approaches have been successful, some eruptions remain unanticipated, resulting in casualties and damage. Our study leverages the recent advancements in ambient field seismology. We explore features extracted from continuous ambient fields using traditional methods, for example, peak ground velocity, peak ground acceleration, root mean square, root median square, real-time seismic amplitude measurement, and novel methods (displacement seismic amplitude ratio and spectral width). In addition, we explore unsupervised learning of higher order wavelet features using scattering networks. We find that combining all the methods was necessary to disentangle the effects of seismic sources from structural changes at Mount St. Helens. Although the ambient wavefield-based approach does not yield additional or more significant precursory signals than event-based methods at Mount St. Helens, our study demonstrates that the ambient wavefield provides supplementary information, mainly about structural changes and complements traditional methods. The ambient seismic wavefield offers additional insights into long-lasting processes. We find enhanced wave attenuation correlating with geochemical measurements. We interpret this as ongoing structural changes, such as dome growth or the evolution of the volcanic conduit system. On annual and decadal timescales, we interpret seasonal seismic attenuation in the shallow subsurface as groundwater fluctuations, corroborated by observations at the nearby Spirit Lake level. This multimethod approach at Mount St. Helens sheds light on a volcanic system’s underlying dynamics and structure.

Washington

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

Colored shaded relief bathymetry and acoustic backscatter of Ozette Lake, Washington

Offshore of the Pacific Northwest of the United States is the Cascadia Subduction Zone, a 1,000-kilometer-long tectonic boundary defined by a large fault, called a megathrust, that extends from the Mendocino Junction off northern California to the Nootka Fracture Zone off Vancouver Island, Canada (U.S. Geological Survey, 2023). The Juan de Fuca and Gorda oceanic plates to the west of this boundary subduct under the North America continental plate to the east. Several other smaller faults that cut through the North America plate crust also affect the region. Although their effects upon Ozette Lake are uncertain, geological evidence for past earthquakes, such as underwater landslides, may be found in Pacific Northwest lakes. Underwater landslides caused by past earthquakes should be well preserved in these relatively undisturbed lake environments. The floor of Ozette Lake, Washington, located along the Pacific coast of the United States, west of the Puget Sound region and about 140 kilometers east of the megathrust was mapped by the U.S. Geological Survey in July of 2019 to search for evidence of past earthquakes. Mapping was completed using a SWATHplus-M 234-kHz interferometric side scan sonar system pole-mounted on the U.S. Geological Survey research vessel San Lorenzo. The system collected full-coverage bathymetric and acoustic backscatter data that were processed to 2-meter spatial resolution (Dartnell and others, 2024). This two-map series displays the results of this mapping. A colored shaded relief bathymetry map (sheet 1) and an acoustic backscatter map (sheet 2) show the lake floor morphology and backscatter intensities, respectively, that can be analyzed for evidence of past earthquakes.

Washington

Field observations and logs from the Rose Hip trench exposure across a north-facing scarp within the Seattle Fault Zone, southern Bainbridge Island, Washington

The Seattle Fault Zone is an approximately 70-km-long, east-west-trending zone of south-dipping blind reverse faults within the Puget lowland region in Washington. Because of the proximity, the Seattle Fault Zone poses a significant earthquake hazard to the Puget sound and Seattle metropolitan regions. We present preliminary mapping and trench-site information from a paleoseismic investigation across a newly identified active fault scarp located within the hanging wall of the Seattle Fault Zone on southern Bainbridge Island, Washington. The trench exposed monoclinally folded Miocene bedrock, fractured and faulted glacial-related deposits, and laminated lacustrine deposits capped by slope-derived colluvium. The observations from this investigation record late Pleistocene to Holocene north-vergent folding and faulting along this new fault scarp.

Washington

Distribution of ancient carbon in groundwater and soil gas from degradation of petroleum near the Red Hill Bulk Fuel Storage Facility, O‘ahu, Hawai‘i

The groundwater below the Red Hill Bulk Fuel Storage Facility (the facility) in Oʻahu, Hawaiʻi, contains fuel compounds from past spills. This study used carbon-14 analyses to distinguish fuel-derived carbon from background carbon, along with other biodegradation indicators, to address two goals: (1) determine the extent and migration direction of groundwater affected by residual fuel below the facility and (2) determine if residual fuel locations in the subsurface could be identified by analyzing soil gas at the surface above the facility. Groundwater from 19 wells was sampled between September 2022 and April 2023. Nonvolatile dissolved organic carbon (NVDOC) from a well presumed to be unaffected by past spills contained 38 percent ancient carbon indicating a natural source of ancient carbon in the subsurface. The NVDOC concentrations and ancient carbon percentages indicate fuel biodegradation products are likely present on the north and south of Red Hill with the greatest effects at well RHMW02 near the 2014 spill site. The NVDOC concentrations are almost three times higher than diesel range organic (DRO) concentrations in groundwater from the same sites. Major ion data indicate that iron reduction is an important biodegradation process. Soil probe samples and soil carbon traps were used to determine the carbon-14 content of soil carbon dioxide. Ancient carbon from fuel biodegradation was not detected at any soil probe or carbon trap site in contrast to a 2017 study which reported ancient carbon detections. A reanalysis of the 2017 results using a range of local values for background carbon-14 indicates that ancient carbon from fuel biodegradation was probably only detected in lower tunnel exhaust system samples and not in any soil carbon trap samples. Measurements of carbon dioxide efflux with a dynamic closed chamber were highly variable. The soil gas results indicate that soil gas measurements at land surface were not useful for detecting residual fuel at the facility.

Hawaii

Time varying crustal anisotropy at Whakaari/White Island volcano

Whakaari/White Island has been the most active New Zealand volcano in the 21st century, producing small phreatic and phreatomagmatic eruptions, which are hard to predict. The most recent eruption occurred in 2019, tragically claiming the lives of 22 individuals and causing numerous injuries. We employed shear-wave splitting analyses to investigate variations in anisotropy during quiescence, unrest, and an eruption between 2018 and 2020. We examined spatial and temporal variations in 7972 shear-wave splitting measurements and derived 1402 V p /IV s ratio measurements. Comparing shear-wave splitting parameters from similar earthquake paths across different times indicates that the observed temporal changes are unlikely to result from variations in earthquake paths through media with spatial variability. Instead, these changes may stem from variations in anisotropy over time, likely caused by changes in crack alignment due to stress or varying fluid content.

Whakaari/White Island volcano