USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Ecology”

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 61Linked to original sources
USGS2024

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↗
USGS2025

A high-resolution 3-D P-wave velocity structure of the south-central Cascadia subduction zone from wide-angle shore-crossing seismic refraction data

This study addresses a significant gap in understanding the features of the south-central Cascadia subduction zone, a region characterized by complex geologic, tectonic, and seismic transitions both offshore and onshore. Unlike other segments along this margin, this area lacks a 3-D velocity model to delineate its structural and geological features on a fine scale. To address this void, we developed a high-resolution 3-D P-wave velocity model using active source seismic data from ship-borne seismic shots recorded on temporary and permanent onshore seismic stations and ocean-bottom seismometers. Our model shows velocity variations across the region with distinct velocity-depth profiles for the Siletz, Franciscan, and Klamath terranes in the overlying plate. We identified seaward dipping high-velocity static backstops associated with the Siletz and Klamath terranes, situated near the shoreline and further inland, respectively. Regions of reduced crustal velocity are associated with crustal faults. Moreover, there is significant along-strike depth variation in the subducting slab, which is about 4 km deeper near the thick, dense Siletz terrane and becomes shallower near the predominantly less-dense Franciscan terrane. This highlights a sudden tectonic and geologic transition at the southern boundary of the Siletz terrane. Our velocity model also indicates slightly increased hydration, though still minimal, in both the oceanic crust and the upper mantle of the subducting plate compared to other parts of the margin.

British Columbia, California, Oregon, Washington↗
USGS2025

Groundwater age estimates for the Mississippi River Valley alluvial aquifer based on tracer data collected during 2018–20

This study characterized groundwater age across the Mississippi River Valley alluvial aquifer (MRVA). Groundwater samples from 69 MRVA wells and 19 wells in Tertiary units of the Mississippi embayment aquifer system (MEAS) were analyzed for sulfur hexafluoride (SF 6 ), tritium ( 3 H), helium (He), and (or) carbon-14 of dissolved inorganic carbon ( 14 C). The age distributions of 89 samples were estimated by fitting lumped parameter models to processed tracer concentrations with the U.S. Geological Survey software TracerLPM. Mean ages of MRVA groundwater samples ranged from 12 to 22,000 years, with a median of 140 years. Mean ages of MEAS groundwater samples ranged from 230 to 52,000 or more years, with a median of 13,500 years. The spatial distribution of MRVA groundwater ages was found to be influenced by depth, inflow of groundwater from deeper units, and soil saturated hydraulic conductivity. In parts of the MRVA, the spatial distribution of MRVA groundwater ages was found to be influenced by annual recharge and (or) annual groundwater pumpage.

Alabama, Arkansas, Illinois, Kentucky, Louisiana, ↗
USGS2024

Wave propagation in layered soil deposits

Recent advances in the general theory of viscoelastic waves and rays in layered media provide a rigorous mathematical framework for site-specific, soil-response models used for earthquake resistant design. The advances provide general closed-form anelastic solutions for the classic problems of the response of a stack of soil layers to S and P waves, ray theory for reflected and refracted waves, Rayleigh- and Love-Type surface waves, and head waves. These general solutions valid for anelastic media regardless of the amount of material damping yield new insights regarding the characteristics of seismic waves and their ray paths that are not provided by conventional models. They provide corresponding numerical ground-response models and ray-tracing computation algorithms that account for changes in velocity and attenuation of anelastic waves associated with changes in inhomogeneity of the waves induced by anelastic soil and soil-rock boundaries. These anelastic effects manifest themselves as variations in amplitude response, amplitude attenuation, raypath location, and travel time as observed at the Earth’s surface. Implications of these anelastic effects for soil-response models used for earthquake resistant design are provided herein.

Conference Paper↗
USGS2025

Evaluating Laramide orogenesis via flexural basin response in the San Juan basin, New Mexico and Colorado

A challenge in interpreting the location, timing, and magnitude of ancient orogenic events is that ongoing uplift and erosion in the hinterlands often destroys much of the primary record of these events. However, basin-thickness patterns in the sedimentary record can provide complimentary evidence of uplift via flexural effects. Here, we deploy well-log correlation, isochores, basin modeling, flexural modeling, and subcrop mapping to evaluate the Late Cretaceous to Paleogene basin response to Laramide tectonism in the San Juan basin. A wedge of upper Campanian to Maastrichtian sedimentary rock thickens from 200 to 800 meters from southeast to northwest in the basin. This pattern can be successfully simulated via flexural modeling if we infer early Laramide uplift along the northwest basin flank that produced a 0.8 km high topographic load. The Laramide unconformity bounds the top of this Upper Cretaceous sedimentary wedge and truncates progressively older strata to the east, further supporting a westward tilt of the basin. The onset of Campanian Laramide flexure may have also contributed to the profound transgression from the upper Menefee Formation to the Lewis Shale. The Paleocene isochore map displays an approximately symmetrical pattern, with thickening towards the center of the basin. This suggests the possibility of competing flexural loads. The base Eocene structure indicates an asymmetric deep on the northeast flank of the basin, providing flexural evidence of contemporaneous uplift/loading of the Nacimiento uplift and Archuleta arch; this has been modeled as ~2.1 km load height. Both Cretaceous and Paleocene sedimentary wedges are narrow, suggesting low flexural rigidity; modeled effective elastic thicknesses (EET) are 20-30 km, comparable to estimates of modern EET for the region.

Colorado, New Mexico↗
USGS2026

Geochemical, mineralogical, and isotopic evidence for multi-stage genesis of the Hicks Dome REE + Y-HFSE-fluorite deposit, Illinois, USA

Hicks Dome hosts breccias enriched in rare earth elements (REE), Y, Th, F, Ba, Ti, Nb, and Be, alongside spatially associated lamprophyre dikes (ca. 271 Ma). Hicks Dome is located within the Illinois–Kentucky Fluorspar District, which hosts fluorite, Pb–Zn, and barite resources. This study investigates the genetic relationships between Hicks Dome mineralization in breccias, alkaline magmatism, and Illinois–Kentucky Fluorspar District mineralization. Lamprophyre dikes are light REE–enriched with chondrite-normalized abundances decreasing from La to Lu. The Host Breccia exhibits middle and heavy REE–enriched patterns that mirror those of the principal REE–Th host minerals, including fluorapatite, xenotime, and thorite. Textural evidence suggests recrystallization of phosphates, sulfates, and Ti–Nb oxides in the Host Breccia. U–Pb geochronology constrains multiple mineralizing events, with ages of 277 ± 18 Ma from low-Th apatite interpreted as main-stage mineralization, and 121.6 ± 9.7 Ma from high-Th apatite indicating later overprinting. O–H–C stable isotope data provide evidence for multiple stages of fluid-rock interaction and fluid mixing: (1) early magmatic fluids dissolved limestone country rock, (2) mixing between magmatic fluids and basinal brines led to main-stage mineralization in the Host Breccia, and (3) late-stage mineralization occurred following mixing of meteoric water and basinal brine. These results indicate that heavy REEs, high field strength elements, and fluorine precipitated proximal to its alkaline magmatic source because of fluid–rock interactions and fluid mixing. Subsequent fluid mixing drove late-stage recrystallization and additional fluorite formation, a process that may be similar to mineralization in the Illinois-Kentucky Fluorspar District.

Illinois↗
USGS2026

Flood-inundation maps of the Blue River and major tributaries in and near Kansas City, Missouri, 2023–25

Digital flood-inundation maps for 35.5 miles of the Blue River, in and near Kansas City, Missouri, were created by the U.S. Geological Survey (USGS) in cooperation with the City of Kansas City. Maps were also developed for a combined 7.3 miles of Brush Creek (4.4 miles) and Indian Creek (2.9 miles), two primary tributaries to the Blue River in the study reach. The flood-inundation maps, available through the USGS Flood Inundation Mapping Program website at https://water.usgs.gov/osw/flood_inundation/ , depict estimates of the areal extent and depth of flooding, depth of water above first-floor elevations in selected buildings, water velocity, water-surface elevations, and time of travel corresponding to selected water levels or gage heights (hereafter referred to as “stages”) at nine reference USGS streamgages—five on the Blue River (USGS station numbers 06893590, 06893578, 06893500, 06893150, and 06893100), two on Brush Creek (USGS station numbers 06893562 and 06893557), one on Indian Creek (USGS station number 06893390), and one on the Missouri River (USGS station number 06893000). Near-real-time stages at these streamgages may be obtained from USGS Water Data for the Nation at https://doi.org/10.5066/F7P55KJN or the National Weather Service National Water Prediction Service at http://water.noaa.gov/ , which also forecasts flood hydrographs at seven of these sites (USGS station numbers 06893000, 06893590, 06893578, 06893500, 06893150, 06893557, and 06893390). Flood profiles were computed for eight map reaches by means of two-dimensional hydraulic models. The models were calibrated using the USGS-developed stage–streamflow relations at each reference streamgage. Two map reaches on the lower Blue River also include the effects of backwater from the Missouri River on flood inundation extent and hydraulic characteristics. The hydraulic models were used to compute water-surface profiles for flood stages at 1-foot intervals referenced to the streamgage datums. The profile stages ranged from the National Weather Service “Action Stage” or near bankfull to a stage exceeding the highest recorded water level at each streamgage. The simulated water-surface profiles were then combined with a digital elevation model (derived from light detection and ranging data having a nonvegetated vertical accuracy of a maximum 10-centimeter root-mean-square error) to delineate the area flooded at each water level and the associated hydraulic characteristics. The availability of these maps, along with information regarding current stage from the USGS streamgage and forecasted high-flow stages from the National Weather Service, will provide emergency management personnel, resource managers, and residents with information that could be critical for flood-response activities such as evacuations and road closures, as well as for postflood recovery efforts.

Kansas, Missouri↗
USGS2026

Climate-adaptive urban planning: Quantitative assessment of drought impacts and practical strategies for climate-resilient urban green spaces

Urban green spaces (UGSs) are vital for enhancing a city’s resilience and livability; however, their functionality is increasingly jeopardized by drought, particularly in water-scarce regions. This study evaluates drought impact on UGSs in Metropolitan Adelaide, Australia, a representative semi-arid urban system, using satellite-derived Normalized Difference Vegetation Index (NDVI) time-series data spanning 2000–2020. Vegetation dynamics were analyzed through Seasonal-Trend decomposition using Loess (STL), standardized anomaly assessment, lagged Pearson correlation, Ordinary Least Squares (OLS) regression, and Mann–Kendall trend analysis. To isolate climatically sensitive signals, 29 urban lawn patches were examined separately from mixed urban canopy, given their shallow root systems and direct dependence on surface moisture. NDVI declined by approximately 0.09 units during the Millennium Drought (2001–2009), with summer greenness deficits reaching 24% below the 20-year benchmark. Temperature was the dominant driver of lawn NDVI variability (r = −0.863, R 2 = 74.5%), substantially exceeding the effect of rainfall (r = 0.156, R 2 = 2.4%). El Niño–Southern Oscillation (ENSO) cycles modulated vegetation responses, with La Niña years supporting recovery and El Niño years amplifying decline. Post-drought recovery remained incomplete, with NDVI deficits of 8–20% persisting through 2020; full recovery was observed only in 2017, coinciding with the highest recorded summer rainfall. No significant directional trend was detected over the full study period (Mann–Kendall τ = 0.005, p = 0.908). These findings demonstrate that heat, rather than water limitation alone, is the primary driver of vegetation stress in urban systems, highlighting the benefits of integrated management strategies that address both warming and moisture deficits to sustain urban green infrastructure under future climate conditions. We introduce the concept of “urban greenery drought,” referring to a form of vegetation stress in managed urban landscapes where greenness is reduced primarily by elevated temperature and atmospheric demand despite water availability.

Adelaide↗
USGS2025

Petrogenesis and mineralization potential of spinifex komatiitic basalts in the Bradley Peak greenstone terrane, Wyoming Province

Komatiitic volcanic rocks are important hosts of Ni sulfide mineralization and record early Earth evolution; however, those in the well-studied Archean Wyoming Province have received little attention. Here, we elucidate the timing and petrogenesis of the Bradley Peak komatiitic volcanic rocks using field and textural observations, geochronology, and geochemistry. Detrital and igneous zircon U-Pb ages for two samples from previously undated units support published age determinations, placing the eruption age at 2.72 Ga. Stratigraphy of the volcanic flows was mapped and 36 samples including cumulates, greenschists, and spinifex-textured rocks were collected. Whole-rock geochemistry was used to classify the spinifex-textured samples as Al-undepleted komatiitic basalts (11–17 wt% MgO). Platinum-group element concentrations (n = 25) are like those in global Al-undepleted komatiitic basalts, and PGE/Ti ratios do not indicate the volcanic flows likely host sulfide mineralization. Initial ε Nd values of −0.5 to +4.7 (n = 16), indicate that these lavas were derived from a depleted mantle source and have negligible evolved crust contamination. The primary magma to the komatiitic basalt flows is estimated to have had 19 wt% MgO and be derived from ∼15 to 25 % mantle partial melting at 3–4 GPa. Trace element chemistry and thermodynamic modeling suggest the primary melt assimilated local banded iron formation. Although the Bradley Peak komatiitic basalts do not contain positive evidence of magmatic sulfide deposits, depleted Au in the flows suggests they could be source rocks for nearby orogenic gold deposits.

Wyoming↗
USGS2026

Microtextural characteristics of adularia in banded quartz veins from the Midas low-sulfidation epithermal deposit, Nevada

High-grade ores at the Miocene Midas low-sulfidation epithermal deposit in northern Nevada are confined to crustiform quartz veins containing abundant adularia. Micro-X-ray fluorescence elemental mapping reveals that adularia is a common gangue mineral occurring in colloform bands, bands showing bladed textures, and bands with dendritic terminations. The adularia aggregates have delicate shapes and are comprised of stacked, submillimeter crystals hosted by fine-grained quartz. The textural evidence suggests that the adularia aggregates originally formed within a gel-like, noncrystalline silica matrix, which subsequently transformed into quartz. This indicates that the adularia did not precipitate in open space along the vein walls. Correlative microscopy, involving scanning electron microscopy-based automated mineralogy and optical petrography, demonstrates that bands containing abundant adularia are not the primary host to ore minerals. The ore minerals occur in different bands within the crustiform veins, implying that adularia and ore mineral precipitation did not always occur simultaneously. It is hypothesized here that fluid flow at Midas involved intermittent short-lived events of fluid flashing, causing rapid solute supersaturation in the liquid. During each flashing event, different amounts of vapor were produced along a given vein. Compositional differences between adjacent bands in the crustiform quartz veins may, therefore, be linked to variations in the amount of vapor formed during each flash event.

Nevada↗
USGS2025

The GorDAS Distributed Acoustic Sensing experiment above the Cascadia locked zone and subducted Gorda Slab

The southernmost portion of the Cascadia Subduction zone in Northern California produces high rates of moderate and large earthquakes owing to subduction of the Gorda slab and deformation associated with the Mendocino Triple Junction. Distributed Acoustic Sensing (DAS) is rapidly advancing as a method for detecting earthquakes and imaging crustal structure. We have begun a long-term DAS monitoring experiment on buried telecom fiber in Arcata, California, with the goal of increasing the available recordings of moderate to large earthquakes as well as imaging seismogenic structures. We have recorded over a year's worth of data, including most aftershocks of the 2022 M w 6.4 Ferndale earthquake, though not the mainshock itself. The dataset includes numerous magnitude 3.5 and larger earthquakes including the 2023/01/01 M w 5.4 Rio Dell earthquake. Here we present initial results comparing an earthquake detection algorithm, run in real-time on the processing unit of the interrogator system, with both the ShakeAlert earthquake early warning system as well as a post-processed earthquake catalog developed with deep-learning phase-picker algorithms. The rapid onboard processing of the detector demonstrates the potential utility of DAS-based edge computing for earthquake early warning. We also verify the quality of the strain waveforms both in terms of peak amplitudes and waveform similarity using about five months of nodal seismometer data. These instruments were deployed roughly every 300 m along the ~15km long cable and validate large variations in peak strain over short distances that are seen in the DAS data. All data from time windows surrounding both the local and teleseismic earthquakes are publicly available, which will improve our understanding of both the performance of DAS systems in moderate earthquakes and earthquake hazards associated with the Gorda subduction zone.

California↗
USGS2026

Neutron scattering reveals fractionation of natural gas mixtures in unconventional petroleum reservoir pores: Perspectives on energy resource recovery and storage

In unconventional petroleum reservoirs hydrocarbon fluids are hosted by both mineral and organic matter pores. These pores can have diameters that range from microns to less than a single nanometer and, for unconventional reservoirs, there is evidence that small pores ( <20 nm diameter) may constitute a large proportion of the available space. Understanding subsurface volumes and how fluids behave in them can be helpful for predicting hydrocarbon production and storage in the subsurface. One area with knowledge gaps regarding hydrocarbon behavior in small pores is the possibility for mixtures to fractionate (i.e., unmix) based on pore size or pore type. Mixture fractionation as a function of pore size could impact recovery of hydrocarbons, drive compositional shifts during production, and limit fluid storage within candidate reservoirs. To investigate natural gas fractionation in small geologic pores, we applied total neutron scattering to probe methane-ethane mixtures at reservoir pressures (up to ≈30 MPa) and temperature (60°C) within a sample from the Upper Cretaceous Niobrara Formation. Neutron scattering data reveal only minor fractionation occurs between methane and ethane in 20-nm diameter sample mesopores. Increased fractionation is observed for sample micropores, with up to 72% (±1% at 1-sigma) methane found in 2 nm diameter pores following injection of a 50%-50% methane-ethane mixture. These data provide rarely available direct experimental observations of hydrocarbon mixture behavior under nanoconfinement in a sample from an important unconventional petroleum reservoir. Our results are discussed in the context of evaluating hydrocarbon resources in unconventional reservoir meso- and micropores, reconciling observed gas composition changes during production, and more broadly, understanding subsurface pore volumes within an energy storage framework.

Fuel↗
USGS2026

2025 USGS National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands: Overview of model and hazard results

The U.S. Geological Survey recently updated the National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). The first version of the PRVI NSHM was released in 2003, and therefore this 2025 update includes over 20 years of new geologic, geophysical, and engineering data, methods, and models. Updates follow similar efforts performed in the recent 2023 50-state NSHM. However, this is the first NSHM in which we: (1) apply an inversion methodology to subduction interface fault sources in the earthquake rupture forecast (ERF) model; (2) develop scaled backbone median ground-motion models and independent aleatory variability models that are applied in the ground-motion characterization (GMC) model; and (3) calculate epistemic uncertainty related to alternative scenarios in the ERF and GMC models for all grid points in the study region. Long-term time-independent mean hazard calculations were performed for peak ground acceleration and 5%-damped pseudospectral acceleration at 21 spectral periods from 0.01- to 10.0-s, for eight National Earthquake Hazards Reduction Program site conditions ranging from V S30 = 150 to 1500 m/s, and for 2%, 5%, and 10% in 50-year probabilities of exceedance (return periods of 2475, 975, and 475 years, respectively). Epistemic uncertainty, in the form of selected percentiles, is also provided for a suite of test sites and all grid points in the study region for limited periods, site conditions, and probabilities of exceedance. Selected results, including comparisons with the 2003 PRVI NSHM, are shown and discussed for selected periods, site conditions, and probabilities of exceedance. When comparing the 2025 PRVI NSHM with the 2003 PRVI NSHM, hazard is generally higher at shorter periods and lower at longer periods, as a result of updates in both ERF and GMC models. The 2025 PRVI NSHM is applicable for return periods greater than ∼475 or less than ∼10,000 years.

Puerto Rico, U.S. Virgin Islands↗
USGS1998

Water-quality assessment of the Ozark Plateaus study unit, Arkansas, Kansas, Missouri, and Oklahoma — Nutrients, bacteria, organic carbon, and suspended sediment in surface water, 1993-95

Nutrient, bacteria, organic carbon, and suspended- sediment samples were collected from 1993-95 at 43 surface-water-quality sampling sites within the Ozark Plateaus National Water- Quality Assessment Program study unit. Most surface-water-quality sites have small or medium drainage basins, near-homogenous land uses (primarily agricultural or forest), and are located predominantly in the Springfield and Salem Plateaus. The water-quality data were analyzed using selected descriptive and statistical methods to determine factors affecting occurrence in streams in the study unit. Nitrogen and phosphorus fertilizer use increased in the Ozark Plateaus study unit for the period 1965-85, but the application rates are well below the national median. Fertilizer use differed substantially among the major river basins and physiographic areas in the study unit. Livestock and poultry waste is a major source of nutrient loading in parts of the study unit. The quantity of nitrogen and phosphorus from livestock and poultry wastes differed substantially among the river basins of the study unit's sampling network. Eighty six municipal sewage-treatment plants in the study unit have effluents of 0.5 million gallons per day or more (for the years 1985-91). Statistically significant differences existed in surface-water quality that can be attributed to land use, physiography, and drainage basin size. Dissolved nitrite plus nitrate, total phosphorus, fecal coliform bacteria, and dissolved organic carbon concentrations generally were larger at sites associated with agricultural basins than at sites associated with forested basins. A large difference in dissolved nitrite plus nitrate concentrations occurred between streams draining basins with agricultural land use in the Springfield and Salem Plateaus. Streams draining both small and medium agricultural basins in the Springfield Plateau had much larger concentrations than their counterparts in the Salem Plateau. Drainage basin size was not a significant factor in affecting total phosphorus, fecal coliform bacteria, or dissolved organic carbon concentrations. Suspended-sediment concentrations generally were small and indicative of the clear water in streams in the Ozark Plateaus. A comparison of the dissolved nitrite plus nitrate, total phosphorus, and fecal coliform data collected at the fixed and synoptic sites indicates that generally the data for streams draining basins of similar physiography, land-use setting, and drainage basin size group together. Many of the variations are most likely the result of differences in percent agricultural land use between the sites being compared or are discharge related. The relation of dissolved nitrite plus nitrate, total phosphorus, and fecal coliform concentration to percent agricultural land use has a strong positive 2 Water-Quality Assessment-Nutrients, Bacteria, Organic Carbon, and Suspended Sediment in Surface Water, 1993-95 correlation, with percent agricultural land use accounting for between 42 and 60 percent of the variation in the observed concentrations.

Arkansas, Kansas, Missouri, Oklahoma↗
USGS2025

Nodal seismic deployment on Mauna Loa volcano, Hawaii: Dataset and preliminary insights

Mauna Loa is the largest active volcano on Earth, comprising ∼51% of the Island of Hawai‘i’s landmass and posing significant risks to the island’s communities, infrastructure, and natural environment. Historical eruptions have produced lava flows that have reached the ocean in as little as 3 hr. The timing and location of such lava flows in the past 200 yr underscore how critical determination of the location and geometry of magma storage and structure is for volcanic hazard assessment and eruption forecasting. Now, after nearly 38 yr of volcanic quiescence, Mauna Loa has erupted again. On 27 November 2022, fissures initiated within the summit caldera and then migrated to the northeast rift zone, where they generated a large lava flow that threatened a major highway. To improve our understanding of the geometry of this magma system, we deployed a temporary nodal array on Mauna Loa in the summer of 2024. This increased our seismic coverage sufficiently to image this magmatic system. This nodal array consists of 33 seismometers distributed on and around the volcano and was deployed for over three months to record seismic signals. The primary objective of this project is to resolve the high‐resolution seismic velocity structure and characterize seismic features associated with magma storage and ascent pathways. In this article, we present an overview of the deployment, evaluate the quality of the data, and show example recordings to evaluate the suitability of the data set for future seismic investigations, including earthquake relocation, seismic tomography, and receiver function analysis. Comparisons with nearby permanent broadband and short‐period seismic stations demonstrate that the nodal array recorded high‐quality waveforms, making it a valuable resource for constraining the magmatic system beneath Mauna Loa at multiple scales.

Hawaii↗
USGS2024

Testing megathrust rupture models using tsunami deposits

The 26 January 1700 CE Cascadia subduction zone earthquake ruptured much of the plate boundary and generated a tsunami that deposited sand in coastal marshes from northern California to Vancouver Island. Although the depositional record of tsunami inundation is extensive in some of these marshes, few sites have been investigated in enough detail to map the inland extent of sand deposition and depict variability in tsunami deposit thickness and grain size. We collected 129 cores in marshes of the Salmon River estuary in Oregon and reanalyzed 114 core logs from a 1987–88 study that mapped the inland extent of circa 1700 CE sandy tsunami deposits. The ca. 1700 CE tsunami deposit in the Salmon River estuary is easily recognized in cores ≤1 m deep in which a buried marsh peat is overlain by a well sorted sand bed with a sharp lower contact that thins and fines inland. We use tsunami deposit data and models of sandy tsunami sediment transport (using Delft3D-FLOW) to test 15 rupture models that could represent a ca. 1700 CE earthquake. At least 12–16 m of slip offshore of the Salmon River, which results in 0.8–1.0 m of coastal coseismic subsidence, is required to match the ca. 1700 CE sand deposit's inland extent, which is consistent with models of heterogeneous megathrust slip in ca. 1700 CE. Our methods of detailed tsunami deposit mapping, combined with sediment transport modeling, can be used to test models of megathrust ruptures and their tsunamis to potentially improve earthquake and tsunami hazard assessments.

JGR Earth Surface↗
USGS1983

Base flow and ground water in upper Sweetwater Valley, Tennessee

The upper Sweetwater Valley area has a flow system with complex interaction between surface and ground water. A water budget study indicated that during dry years approximately three-fourths of the annual flow to Sweetwater Creek may be derived from ground-water sources. Hydrograph analysis showed seasonal variation of recharge to the ground-water flow system. Streamflow records were analyzed to estimate the frequency of low flow of Sweetwater Creek at river mile 16.7, and indicated the lowest average flow for 1 day in 20 years to be about 5.1 cubic feet per second. Two periods of base-flow measurements of Sweetwater Creek identified channel reaches with significant gains and losses of streamflow. Base flow measurements also showed interbasin transfer of water among sub-basins of the valley. Major flow surpluses were associated with areas in which the majority of flow originated at a spring. Topographically low areas adjacent to the main stem of Sweetwater Creek generally have surplus flow. Topographically higher areas generally have deficient surface outflow unless significant spring flow occurs in the basin. Ground-water recharge occurs by water draining into sinkholes, faults, and fractures. Ground-water flow is regionally diffused across formation strikes from the topographically low areas unless the water is exposed to highly permeable formations or impervious formations. Ground water infiltrates the highly permeable formations and flows along strike. Ground water encountering impervious formations may discharge at small springs at the contact, or may reroute along the contact if the upgradient rock is sufficiently permeable or has well developed secondary porosity. Groundwater discharges to streams at innumerable seeps and at a few large springs. Areas of ground-water flow up-gradient of large springs are hypothesized as likely areas of significant ground-water reservoirs.

Tennessee↗
USGS2025

Structural setting and geothermal potential of northeastern Reese River Valley, north-central Nevada: Highly prospective detailed study site for the INGENIOUS project

The northeastern part of the Reese River basin situated ~15 km southeast of Battle Mountain, Nevada, scored highly in the Nevada geothermal play fairway analysis (PFA) for hosting potential hidden geothermal systems. This site (also referred to as Argenta Rise) was therefore chosen for detailed study in the INGENIOUS project (INnovative Geothermal Exploration through Novel Investigations Of Undiscovered Systems). The high PFA scores resulted primarily from favorable structural settings (e.g., fault intersections and pull aparts) with relatively high slip rates on Quaternary faults. The INGENIOUS project is utilizing additional parameters and more rigorous analytical techniques to further advance exploration at this site. This includes integration of geological (e.g., Quaternary fault mapping) and new geophysical datasets (e.g., gravity, magnetics, MT data, and five reprocessed seismic reflection profiles) to build a structural model and to identify specific favorable sites for potential geothermal upwellings. Two-meter temperature surveys were also conducted in the area (139 measurements). This part of north-central Nevada is characterized by systems of intersecting northerly and ENE-striking faults within the broader Humboldt structural zone, a poorly understood belt of ENE-striking faults and relatively high heat flow extending across northern Nevada. Kinematic analysis of exposed fault surfaces shows that ENE-striking faults have accommodated sinistral-normal slip, and normal slip characterizes N- to NNE-striking faults. Northeastern Reese River Valley lies within a broad left step between major ENE-striking fault zones on the northern flanks of the Argenta Rim and Shoshone Range and thus corresponds to a broad pull-apart in the ENE-striking sinistral-normal fault system. Notably, the nearby Beowawe geothermal system in Whirlwind Valley (with abundant sinter, hot springs, and a geothermal power plant) occupies a fault intersection in a relatively small left step in a major ENE-striking sinistral-normal fault and may serve as an analogue for a potential hidden system in northeastern Reese River Valley. Existing geological maps, high-resolution lidar, and seismic reflection data demonstrate that northeastern Reese River Valley is structurally complex with multiple intersections between the ENE- and N- to NNE-striking fault systems. Some of these fault intersections correspond to low resistivity anomalies, magnetic lows, and/or very subtle 2-m temperature anomalies, which may indicate hidden geothermal upwellings. Three-dimensional modeling and temperature-gradient drilling are planned to further evaluate these sites for geothermal activity.

Nevada↗