USGS Science⌕ Search

SEARCH · USGS Science

Results for “Environmental Monitoring and Assessment”

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,662 records · Page 78Linked to original sources
USGS2024

U.S.-Mexico Borderland & vegetation community map

People on both sides of the United States-Mexico border need a high-resolution, binational vegetation community map that spans the entire United States-Mexico borderlands. Traditionally, mapping efforts in this region were impeded by complex logistics related to the international border, differing national needs and plans, and resource allocations and priorities. To address this need, scientists from the U.S. Geological Survey (USGS) Southwest Biological Science Center partnered with the Sonoran Joint Venture, the U.S. Fish and Wildlife Service (FWS) Migratory Bird Program, data engineers from the Department of Biosystems Engineering at the University of Arizona, and collaborators from the Wildlands Network, the Borderlands Program to produce the first prototype land cover map within the overlapping Mojave Desert, Sonoran Desert, and the North American Bird Conservation Initiative’s Bird Conservation Region 33 (BCR33) using Landsat satellite data . BCR33 is an area of high biodiversity, providing habitat for bird species of concern and other wildlife. The land cover map supports FWS recovery plan efforts related to conservation planning activities for many species, including Yellow-billed Cuckoo ( Coccyzus americanus ), Cactus Ferruginous Pygmy-Owl ( Glaucidium brasilianum cactorum ), Southwestern Willow Flycatcher ( Empidonax traillii extimus ), Yuma Ridgway’s Rail ( Rallus obsoletus yumanensis ), Bendire’s thrasher ( Toxostoma bendirei ), LeConte’s thrasher ( Toxostoma lecontei ), Masked Bobwhite ( Colinus virginianus ridgwayi ), jaguar ( Panthera onca ), and endangered plants such as Bartram’s stonecrop ( Graptopetalum bartramii ) and the Pima pineapple cactus ( Coryphantha robustispina ssp. robustispina ). In 2024, a Phase-II map for the full BCR33 region was completed, increasing the understanding of the binational nature of natural communities. The published map and associated paper can be found here .

Borderland↗
USGS2023

Slow slip detectability in seafloor pressure records offshore Alaska

In subduction zones worldwide, seafloor pressure data are used to observe tectonic deformation, particularly from megathrust earthquakes and slow slip events (SSEs). However, such measurements are also sensitive to oceanographic circulation-generated pressures over a range of frequencies that conflate with tectonic signals of interest. Using seafloor pressure and temperature data from the Alaska Amphibious Community Seismic Experiment, and sea surface height data from satellite altimetry, we evaluate the efficacy of various seasonal and oceanographic pressure signal proxy corrections and conduct synthetic tests to determine their impact on the timing and amplitude prediction of ramp-like signals typical of SSEs. We find that subtracting out the first mode of the complex empirical orthogonal functions of the pressure records on either the shelf or slope yields signal root-mean-square error (RMS) reductions up to 73% or 80%, respectively. Additional correction with proxies that exploit the depth-dependent spatial coherence of pressure records provides cumulative variance reductions up to 83% and 93%, respectively. Our detectability tests show that the timing and amplitude of synthetic SSE-like ramps can be well constrained for ramp amplitudes ≥4 cm on the shelf and ≥2 cm on the slope, using a fully automated detector. The principal limits on detectability are residual abrupt changes in pressure that occur as part of the transition to and from summer to winter conditions but are not adequately characterized by our seasonal corrections, as well as the inability to properly account for instrumental drift, which is not readily separated from the seasonal signal.

Alaska↗
USGS2024

Limited evidence of late Quaternary tectonic surface deformation in the eastern Tennessee seismic zone, USA

The ~300-km-long eastern Tennessee seismic zone (ETSZ), USA, is the second-most seismically active region east of the Rocky Mountains. Seismicity generally occurs below the Paleozoic fold-and-thrust belt within the Mesoproterozoic basement, at depths of 5–26 km, and earthquake magnitudes during the instrumental record have been moment magnitude ( M w )≤4.8. Evidence of surface deformation may not exist or be difficult to detect because of the vegetated and soil-mantled landscape, landslides, locally steep topography, anthropogenic landscape modification, or long, irregular recurrence intervals between surface-rupturing earthquakes. Despite the deep seismicity, analog models indicate that accumulation of strike-slip or oblique-slip displacement at depth could be expected to propagate upward through the Paleozoic section, producing a detectable surficial signal of distributed faulting. To identify potential surface deformation, we interrogated the landscape at different spatial scales. We evaluated morphotectonic and channel metrics, such as channel sinuosity and catchment-scale hypsometry. Additionally, we mapped possible fault-related topographic features on 1-m lidar. Finally, we integrated our observations with available bedrock and Quaternary surficial mapping and subsurface geophysical data. At a regional scale, most morphotectonic and channel metrics have a strong lithologic control. Within smaller regions of similar lithology, we observe changes in landscape metrics like channel sinuosity and catchment-scale hypsometry that spatially correlate with new lineaments identified in this study and previously mapped east–west Cenozoic faults. These faults have apparent left-lateral offsets, are optimally oriented to slip in the current stress field, and match kinematics from recent focal mechanisms, but do not clearly preserve evidence of late Pleistocene or Holocene tectonic surface deformation. Most newly mapped lineaments might be explained by either tectonic or non-tectonic origins, such as fluvial or karst processes. We also re-evaluated a previously described paleoseismic site and interpret that the exposure does not record evidence of late Pleistocene faulting but instead is explained by fluvial stratigraphy.

Tennessee↗
USGS2025

Broadband stochastic simulation of earthquake ground motions with multiple strong phases with an application to the 2023 Kahramanmaraş, Turkey (Türkiye), earthquake

Stochastic ground motion simulation models are often less accurate at lower frequencies than at higher frequencies when fitting recorded data unless supplemented by a deterministic forward directivity velocity pulse model. Moreover, time-modulated stochastic models, which adjust ground motion amplitudes over time, typically use functions that fail to capture multiple strong-motion phases. The February 2023 Turkey (Türkiye) earthquake exhibited diverse recordings, including near-fault and far-field motions with pulse-like and non-pulse-like characteristics, along with single and multiple strong-motion phases. To better represent such a diverse set of recordings, this study enhances a fully non-stationary site-based stochastic model without combining it with a deterministic model. Improvements include a new band-pass filter with upper- and lower-frequency limits, which refines the representation of the low-frequency content. Moreover, a time-modulating function that can represent energy arrival in multiple strong phases is introduced. The reference model’s parameters are identified by fitting to the energy content, zero-level crossings, and cumulative counts of positive-minima and negative-maxima of a target accelerogram. This fitting procedure is modified to address the increased number of parameters. These improvements broaden the reference model’s applicability while preserving its simplicity, a key aspect appealing to engineering practitioners. The improved model’s applicability is demonstrated by simulating a dataset from the February 2023 Türkiye earthquake, and the accuracy is tested using a pulse-like Next Generation Attenuation Relationships for Western United States dataset. Validations are performed based on total energy, zero-level crossings, Fourier amplitude spectrum, elastic response spectra, and peak ground motion parameters. Validations are performed schematically in the time and frequency domains and quantitatively using goodness-of-fit scores, various validation-metrics errors, and inter-period correlations. Overall, the improved stochastic model can effectively simulate a set of diverse ground motion recordings, including near-fault pulse-like records, records with multiple strong phases, and far-field motions across a broad frequency range.

Kahramanmaraş↗
USGS2026

Evaluation of stream capture related to groundwater pumping, middle Humboldt River Basin, Nevada

Historical, future, and potential stream capture from groundwater pumping in the middle Humboldt River Basin (MHRB), Nevada, is estimated using a calibrated numerical groundwater flow model. The model was developed to estimate (1) stream capture, which is the change in flux between the groundwater system and the Humboldt River and tributaries, and (2) change in streamflow, which is the change in streamflow estimated for the Imlay gage on the Humboldt River (U.S. Geological Survey streamgage 10333000). Historical stream capture for water years (WYs) 1961–2015 is estimated using recorded and estimated groundwater pumping during that period. Future (predictive) stream capture was based on historical stresses (WYs 1961–2015) using a scenario that simulated non-mine pumping from WY 2015 at a uniform rate for 100 years into the future. Potential stream capture throughout the middle Humboldt River Basin from groundwater pumping during varying durations of time are presented in a series of capture maps. Maps also are presented that show the potential to capture from groundwater evapotranspiration, as well as the storage changes for pumping duration of 100 years. Estimates of historical stream capture from the mainstem Humboldt River during the early 1960s are less than 400 acre-feet per year (acre-ft/yr) when groundwater withdrawals and pumping rates were relatively small compared to more recent times. In the late 1980s and early 1990s, groundwater withdrawals increased and estimated historical stream capture also increased from about 4,000 acre-ft/yr in the late 1980s and early 1990s to as much as 18,800 acre-feet (acre-ft) in WY 1998. In WY 2015, estimated historical stream capture declined to about 13,000 acre-ft because of decreasing groundwater withdrawals and lower streamflow during the drought of WYs 2012–15, resulting in less stream water available for capture. Stream capture was estimated for 100 years into the future based on WY 2015 non-mine pumping rates and mine-dewatering activity through WY 2015. Stream capture is forecast to increase to about 23,000 acre-ft/yr, and streamflow in the Humboldt River could decrease by as much as 19,000 acre-ft/yr. Pumping for mine-dewatering and the associated discharge of that water affects streamflow in the Humboldt River at Imlay, Nevada (U.S. Geological Survey streamgage 10333000). Historically, from WYs 1991 to 2015, streamflow was greater at Imlay gage during active mine-dewatering from mine-water discharge operations and increased by as much as 105,000 acre-ft in WY 1998. The increase was attributed mostly to the discharge of groundwater from mine-related dewatering operations directly into the mainstem Humboldt River or its tributaries, with some of this increase associated with return flows from discharge to rapid infiltration basins. Results indicate that streamflow at Imlay gage is expected to decrease by as much as 1,600 acre-ft/yr 30 years after mine-related pumping and discharge are discontinued. The streamflow reductions at the Imlay gage are expected to then decrease to around 500 acre-ft/yr, 100 years after mine-related pumping and discharge are discontinued. Potential capture maps were produced for pumping durations of 10, 25, 50, and 100 years. Capture map results indicate that areas of greater potential stream capture occur adjacent to the Humboldt River and for upstream tributaries areas north of the Humboldt River.

Nevada↗
USGS2024

Upper crustal seismic velocity structure of the Hayward fault zone, San Francisco Bay, California, USA: Results from the 2016 East Bay Seismic Experiment (EBSI-16)

We developed Vp, Vs, Vp / Vs ratio, and Poisson’s ratio models of the uppermost crust (<4 km depth) from the eastern San Francisco (SF) Bay (California, USA) to near the Calaveras fault along a 15-km-long, linear profile. Upper crustal velocities are highly variable beneath, west, and well east of the Hayward fault. We observe eight notable features, from west to east: (1) Near San Francisco Bay, there is an ~2-km-wide structure with high Vp / Vs ratios (up to 5) and Poisson’s ratios (up to 0.48) extending from the surface to the base of our model, which we suggest the structure is a near-vertical fault that lies along a straight-line projection between the Silver Creek fault to the south and the Point Richmond fault to the north. The structure may be part of an ~90-km-long fault along the eastern SF Bay. (2) The western East Bay Plain, the lower lying area between the bay and the hills, includes up to 800 m of low-velocity sediments ( Vp ~1600–3000 m/s, Vs ~500 m/s to ~1000 m/s), underlain by higher velocity basement rocks ( Vp ~3000–5800 m/s; Vs ~1000–1500 m/s). (3) Between ~1 km and 3 km east of the Bay shoreline, sediments thin in a series of steps (likely faults) toward the Hayward fault. (4) Between ~3 km west and ~1 km east of the Hayward fault (at the East Chabot fault) at depths greater than 1 km, basement Vp (up to 6000 m/s) and Vs (up to 2800 m/s) are high, and Vp / Vs ratios (<2) and Poisson’s ratios (<0.3) are low, suggesting crystalline rocks. Furthermore, a near-vertical zone of low Vp / Vs ratios and Poisson’s ratios is between near-surface traces of the Hayward and East Chabot faults, likely corresponding to the San Leandro Gabbro of Ponce et al. (2003) . (5) Eastward of the East Chabot fault in the upper 1.5 km, basement Vp (~3000 m/s to ~4200 m/s) and Vs (~1200–2000 m/s) are lower than those west of the fault. (6) In the eastern Hayward/Oakland Hills, there are zones of laterally varying, high- and low-velocity ( Vp ~2500–3000 m/s) Jurassic–Cretaceous and Tertiary sediments in the shallow subsurface that likely extend much deeper than imaged. (7) Seismic energy that propagates westward from sources east of the Hayward fault (HF) appear weaker than energy that propagates eastward from sources west of the HF, suggesting that the HF acts as a partial barrier to shallow seismic energy propagation into the more populated eastern SF Bay area. (8) Unlike many fault zones, it appears that the active trace of the Hayward fault (in our study area) is not cored by a prominent, low-velocity zone relative to rocks to the east and west of the active trace. However, the active trace does mark a prominent change from relatively higher velocities to the west and lower velocities to the east.

California↗
USGS2024

Timing and geometry of the Chemehuevi Formation reveal a late Pleistocene sediment pulse into the Lower Colorado River

The Chemehuevi Formation is a distinctive 50−150-m-thick wedge-shaped Pleistocene sedimentary unit deposited by the Colorado River. It lines the perimeters of the river’s floodplains and bedrock canyons for more than 600 km between the mouth of the Grand Canyon and the delta region in the Gulf of California. The formation is composed of a basal tan to light-yellowish-brown and pale-orange mud-dominated facies overlain and interbedded by a light-yellow-brown sand-dominated facies. The unit is one of two extensively exposed aggradational packages in the Lower Colorado River corridor, in addition to a series of other smaller alluvial terrace deposits. The Chemehuevi Formation appears to represent the response of a fully integrated Colorado River system to a significant perturbation, in contrast to the Bullhead Alluvium, which is likely a unique result of Pliocene river integration. The aggradation of the Chemehuevi Formation in the Lower Colorado River corridor may be similarly due to a unique event in the Colorado River system, or it may instead be a well-preserved sedimentary sequence recording typical behavior of the Colorado River below the Grand Canyon in the late Pleistocene. As such, multiple causal mechanisms have been proposed, but no study to date has conclusively explained the Chemehuevi Formation. To help resolve its timing, duration, and origin, we applied post-infrared infrared stimulated luminescence, carbonate U-Th series, and zircon sensitive high-resolution ion microprobe U-Th series geochronology to determine the ages of key exposures of the unit over a wide spatial area. These new data demonstrate that the Chemehuevi Formation was deposited ca. 110−90 ka. The depositional ages collectively overlap, suggesting that deposition occurred rapidly relative to the resolution of the geochronometers. The new depositional timing coincides with a shift from glacial to interglacial conditions after the marine isotope stage 5-6 transition. This observation is consistent with a climate-induced sediment pulse as a causal mechanism, yet correlations with similar deposits in the Colorado River headwaters or in neighboring catchments appear elusive. Potentially, climate transitions between glacial and interglacial periods induced a sediment pulse from hillslopes of the Colorado River system that resulted in the Chemehuevi Formation. An alternative or additional explanation is that the Chemehuevi Formation represents release of lava dam−impounded sediment in the Grand Canyon. The surface geometry of the Chemehuevi Formation projects upstream to the approximate location of lava dams, and the largest possible lava dam impoundment (the Upper Prospect dam) is comparable in volume to the formation. The lava dam hypothesis appears to be a possible explanation for the Chemehuevi Formation. However, tying deposition to a specific lava dam or series of lava dams remains challenging due to discrepancies in timing and volume. The combined effects of a series of lava dams may have led to the Chemehuevi Formation, as the last Pleistocene lava dam eruption coincides with the onset of deposition. Alternatively, the formation may result from the combined effects of both regional climate transitions and the lava dams that created a transient reservoir to compound a climate transition−driven sediment pulse. The geochronologic data presented here do not allow us to distinguish between the lava dam or climate transition hypotheses but will need to be reconciled with any future proposed depositional model.

Arizona, California, Nevada↗
USGS2024

Using stable oxygen isotope dual-inlet isotope-ratio mass spectrometry to elucidate uranium transport and mixed 230Th/U calcite formation ages at the seminal Devils Hole, Nevada, natural laboratory

Rationale Vein calcite in Devils Hole has been precipitating continuously in oxygen-isotope equilibrium at a constant temperature for over 500 000 years, providing an unmatched δ 18 O paleoclimate time series. A substantial issue is that coeval calcite (based on matching δ 18 O values) has uranium-series ages differing by 12 000 years. Methods An unparalleled high-accuracy δ 18 O chronology series from continuously submerged calcite was used to correct the published uranium-series ages of non-continuously formed calcite in two cores, cyclically exposed by water-table decline during glacial–interglacial transitions. This method relies on the premise that the δ 18 O values of coevally precipitated calcite are identical, allowing matching calcite δ 18 O values to establish formation ages. Results Exposed calcite can have apparent ages that are 12 000 years too young due to unrecognized uranium mobility and resulting mixed ages identified in over 50 mixed uranium-series ages from previous studies. Secondary uranium in fluids, sourced from the formation or dissolution of porous carbonate deposits (folia) with high uranium-238 ( 238 U) concentrations, has migrated up to 10 mm into vein calcite. Conclusions The continuously submerged Devils Hole δ 18 O chronology is not explained by orbital forcing. Rather, this chronology represents a regional climate record in the southern Great Basin of sea-surface-temperature (SST) variations off California, variations that preceded the last and penultimate deglaciations by 5000 to approximately 10 000 years. Temporal discrepancies between the continuously submerged Devils Hole chronology and other regional δ 18 O records (e.g., the Leviathan chronology) can be explained by unrecognized cryptic, pernicious uranium mobility, leading to model estimations that may be thousands of years younger than actual ages. Consequently, paleo-moisture availability, water-table, and groundwater recharge models based on these mixed uranium-series ages are too young by as much as 12 000 years. The potential for post-formation uranium addition in subaerial cores and speleothems underscores the need for caution in uranium-series dating, highlighting δ 18 O time-series comparisons as a method for identifying mixed ages.

Nevada↗
USGS2026

Revisiting the utility of regional-scale, high-quality geophysical data in mineral exploration - A case study featuring the Mammoth Magnetic Anomaly, Pinal County, Arizona

Regional aeromagnetic surveys passively measure the total magnetic intensity (TMI) and are a foundational tool used in mineral exploration (Airo, 2015). With the increased global demand and the number of critical mineral resources required for manufacturing high-tech devices, developing high-quality, regional-scale geophysical surveys could aid critical mineral exploration efforts and geologic mapping. In 2019, the U. S. Geological Survey launched the Earth Mapping Resources Initiative (Earth MRI) to modernize the geologic and geophysical mapping of regions that have the potential to contain critical mineral resources within the United States. In support of planning Earth MRI geophysical surveys, Drenth and Grauch (2019) defined five aeromagnetic data quality rankings (rank 1 through rank 5) applying them to the airborne geophysical survey inventory of the United States (Johnson et al., 2021). Rank 1 aeromagnetic surveys are of the highest quality, meeting modern standards and allowing best practices for qualitative and quantitative interpretation; whereas rank 5 aeromagnetic surveys are of the lowest quality, being useful only for qualitative interpretation of broad features. Through the Earth MRI effort, 48 high-quality, regional-scale rank 1 and 2 airborne magnetic and radiometric geophysical surveys have been planned, collected, or publicly release through May 2025 (U. S. Geological Survey, 2025). Here, a portion of a rank 1 Earth MRI aeromagnetic survey in southeast Arizona is presented and compared to a legacy rank 5 aeromagnetic survey over the Mammoth Magnetic Anomaly (MMA), demonstrating how modern, high-quality aeromagnetic data improves our view of crustal geology, aiding mineral exploration.

Arizona↗
USGS1976

Neutron activation analysis in hydrology

Neutron activation analysis in hydrogeochemical and related investigations is increasingly practical because of developments in radiation detectors and computerized data processing. While neutron activation is the most widely used kind of activation analysis other kinds (for example, charged particle activation) are available for special applications. Neutron activation analysis is promising because approximately 70 percent of the elements in the periodic table are detected with satisfactory sensitivity, contamination from reagents is eliminated, simultaneous determination of more than 20 elements is often possible, technique is versatile, is applicable to small samples and is non-destructive in the instrumental form. Developments in theU.S. Geological Survey include a special sampling technique for water; a comprehensive technique using carriers precipitated with sulfide to concentrate desired trace elements; a determination of mercury, silver, and gold; determination of aluminum, vanadium and nine other elements; determinations of two groups of rare earth elements; and adaptations of the analytical methods to sediment and biota.

Open-File Report↗
USGS2024

Thermo-hydrologic processes governing supra-permafrost talik dynamics in discontinuous permafrost near Umiujaq (Québec, Canada)

Widespread supra-permafrost talik formation is currently recognized as a critical mechanism that could accelerate permafrost thaw in the Arctic (e.g., Connon et al. 2018; Farquharson et al. 2022). However, the trajectory of permafrost dynamics following talik formation may prove difficult to predict. Physically-based cryohydrogeologic models provide a powerful tool for understanding processes and factors controlling talik dynamics and, ultimately, how permafrost will respond to climate change. Such models are typically used to represent multiple non-linear processes relevant for groundwater systems in cold regions, such as coupled heat and groundwater movement, including freeze-thaw dynamics and the effects on the surface energy balance and the subsurface thermal and hydraulic properties (Lamontagne-Hallé et al. 2020). Though cryohydrogeologic modeling advances have been made in simulating talik dynamics, few applications have been tested against robust long-term hydrometeorological and subsurface observations.

Quebec↗
USGS1976

The water quality of Livingston Reservoir on the Trinity River, Southeastern Texas

The concentrations of dissolved solids, chloride, and sulfate in Livingstone Reservoir on the Trinity River in southeastern Texas usually average less than 250 mg/1 (milligrams per litre), 40 mg/1, and 50 mg/1, respectively. The water is usually hard or moderately hard (61 to 180 mg/1 as calcium carbonate). The concentrations of principal dissolved constituents in the reservoir are usually maximum during summer and fall when evaporation is high and inflow is low. Thermal stratification of the reservoir usually begins in March and persists until September or October. Neither the seasonal variation of dissolved constituents in inflow to the reservoir nor thermal stratification has resulted in significant stratification of the principal dissolved constituents. However, thermal stratification has resulted in significant seasonal and areal variations of dissolved oxygen, which results in higher concentration of dissolved iron, dissolved manganese, total phosphorus, and total inorganic nitrogen.

Texas↗
USGS2025

Prospectivity mapping for geologic hydrogen

Geologic, or naturally occurring, hydrogen has the potential to become a new, low-carbon, primary energy resource. Often referred to as “white” or “gold” hydrogen, this gas occurs naturally in the Earth’s subsurface, similar to petroleum resources. However, unlike petroleum, which releases carbon dioxide when burned, burning hydrogen only produces water as a byproduct. Exploration for geologic hydrogen remains in an early stage and discoveries of high concentrations of subsurface hydrogen are still relatively rare. To facilitate research and exploration for this potential resource, this report presents the first publicly available prospectivity map of geologic hydrogen accumulations in the conterminous United States. Prospective regions are those regions in which all major components necessary for a hydrogen accumulation likely are present—a source of sufficient hydrogen generation, porous reservoirs for storage, and seals to prevent leakage. The midcontinent region of the United States and the central California coast are revealed as having high prospectivity. This analysis also identifies previously unrecognized prospective regions that may be favorable due to long distance lateral migration of subsurface hydrogen, such as the offshore eastern seaboard of the United States, and can provide a linkage between surface observations of hydrogen degassing and far-field source regions. The methodology developed to create this map is expandable and flexible and may be adapted to incorporate new concepts in the hydrogen system and for application to other regions of the world.

conterminous United States↗
USGS2024

Distinguishing natural sources from anthropogenic events in seismic data

As seismic data are increasingly used to investigate a diverse range of subsurface phenomena beyond regular fast-rupturing earthquakes (Peng and Gomberg, 2010; Beroza and Ide, 2011), it is important to acknowledge that human-generated ground vibrations may be mistaken for naturally generated subsurface processes (Larose et al., 2015; Li et al., 2018). Correct discrimination of natural processes from anthropogenic noise is especially pressing given the trend in seismic detection research toward automated algorithms and machine learning methods (Yoon et al., 2015; Kong et al., 2019;Mousavi and Beroza, 2022) and the growth in seismic data collection in new environments such as urban and industry settings (e.g., Díaz et al.,2017).

Seismological Research Letters↗
USGS2025

Rapid emplacement of the Keaiwa Lava Flow of 1823 from the Great Crack in the Southwest Rift Zone of Kilauea volcano

The Keaīwa Lava Flow of 1823 in the Southwest Rift Zone of Kīlauea volcano is unusual for its expansive pāhoehoe sheet flow morphology and lack of constructive vent topography, despite having a similar tholeiitic basalt composition to other lavas erupted from Kīlauea. This lava flow issued from a ∼10-km-long continuous fissure now known as the Great Crack, and has an unusually thin sheet flow morphology with margin thicknesses of ∼15–110 cm (average of 42 cm). Based on field observations of the lava flow at its fissure vent (e.g., drain-back features), we propose that the Great Crack formed, or at least significantly widened, just prior to and syn-eruptively with this 1823 eruption. The absence of pyroclastic cones or spatter ramparts indicates that the eruption consisted of a rapid outpouring of relatively degassed lava as the fissure unzipped. The rapidly moving lava flow overtopped pre-existing tumuli and scoria cones (e.g., Lava Plastered Cones) up to ∼10 m tall. Glass and whole-rock chemistry yield homogeneous compositions for the lavas erupted from the Great Crack, with glass compositions of 6.40 ± 0.10 wt% MgO and whole-rock compositions of 7.39 ± 0.07 wt% MgO. Lava pads erupted from a short western fissure system are richer in mafic minerals (e.g., olivine and clinopyroxene), and show slightly more MgO-rich whole-rock compositions (7.79 ± 0.05 wt%). MgO-in-glass thermometry on juvenile spatter yield eruption temperatures of 1153 ± 13°C that are typical of Kīlauea lavas. Thus, the extensive sheet-like lava flow morphology is not a direct consequence of unusual magmatic or rheological conditions (i.e., low viscosity). Instead, the flow morphology is associated with high effusion rates caused by sudden drainage of uprift magma as it erupted from the Great Crack. Lava flow modeling on a 2-m-resolution digital elevation model indicates that a minimum bulk effusion rate of ∼5800 m 3 /s (∼3500 m 3 /s dense rock equivalent) and a minimum flow velocity of ∼11 m/s are required for the lava flow to overcome the topography of the Lava Plastered Cones. This effusion rate is among the highest inferred for eruptions in Hawaiʻi and around the world. This study highlights a less frequent eruption style at Hawaiian volcanoes characterized by a sudden outpouring of lava from an unusual fissure system. Local eyewitness accounts indicate that the 1823 eruption was preceded by seismicity. Given the complex magmatic-volcanic-tectonic relations across Kīlauea, we speculate that the south flank could have slipped over one or more events that ultimately triggered unzipping of the Great Crack and passive release of briefly stored uprift magma. An eruption similar to 1823 at Kīlauea or Mauna Loa, with an eruptive timeframe that could be as short as an hour, with high effusion rates and rapid flow front velocities, would not easily allow for a timely response.

Hawaii↗
USGS2025

SURF: An automated method for building nonplanar 3D fault models from earthquake hypocenters

Accurately characterizing 3D fault geometry is vital for improving our understanding of earthquake behavior and informing the development of seismic hazard models. Despite their importance, subsurface fault structures tend to be poorly constrained because of limitations in observational data. Improvements to the seismic networks and earthquake detection algorithms have increased the precision and volume of earthquake catalogs, which help illuminate detailed subsurface fault structure and provide the most direct information available about fault geometries at depth. We present a Python package to automate generating 3D fault geometries directly from hypocentral seismicity patterns. This method begins with clustering events based on their spatial density, identifying coherent patterns. Nearby clusters are then merged based on the similarity of their orientations. We fit nonplanar surfaces using support vector regression to balance surface accuracy with minimal deviations from planarity. The fault models are output as quadrilateral meshes at user‐defined resolution. In the process of generating the 3D fault surfaces, we compute the spatial density of seismicity around the surface and the planarity as quantitative metrics of the model outputs. As a proof of concept, we apply this approach to the San Andreas–Calaveras fault junction region and the 2019 Ridgecrest earthquake sequence, both in California, which contain complex subparallel faults well defined at the Earth’s surface and abundant microseismicity. These case studies demonstrate the method’s ability to model complex fault structures, including long continuous fault surfaces, crossing faults, variably dipping segments, and subparallel faults. We test the method on both standard network catalogs and double‐difference relocated catalogs. We find that our seismicity‐based fault model results align with published 3D models that incorporate additional constraints and interpretations ( Plesch et al. , 2020 ; Aagaard and Hirakawa, 2021 ). This workflow provides a low‐user‐input solution for estimating fault geometries at depth from earthquake catalogs.

California↗
USGS1988

Geology of an area near Brentwood, Williamson County, Tennessee

The lithology and geologic structure of an area near Brentwood, Williamson County, Tennessee, were studied to define the potential aquifers and confining units of the ground-water flow system. Four formations were identified. They are, in descending order, the Bigby-Cannon Limestone, the Hermitage Formation, the Carters Limestone, and the Lebanon Limestone. The Bigby-Cannon Limestone and the Carters Limestone are potential aquifers. The Hermitage Formation and the Lebanon Limestone are confining units. The Bigby-Cannon Limestone and the Hermitage Formation have been affected by recent erosion. Any variation of the Carters Limestone is controlled by pre-Carters erosion of the top of the Lebanon Limestone. The thickness of the Carters Limestone ranges from 65 to 79 feet. Structurally, the area reflects the regional northwest dip. A small scale anti-cline-syncline pair also is evident. This feature is not a result of erosion and also occurs in the T-3 bentonite bed in the Carters Limestone.

Tennessee↗
USGS2026

Magmatic source of the opening phase of Kīlauea’s 2018 Lower East Rift Zone eruption

The 2018 eruption of Kīlauea volcano in its Lower East Rift Zone began with the discharge of evolved high-Ti basalt as weak lava fountains and short, slow-moving lava flows. The lavas were quickly geochemically recognized as being derived from magmas stored within the rift zone and remobilized by a new intrusion, a sequence that is common at Kīlauea. This initial phase of the 2018 eruption, referred to as phase 1a, lasted for 6 days and was followed by extrusion of mixed magma after a 3-day pause. Even though remobilization of older rift zone magmas is common within Kīlauea’s rift zones, it is difficult to determine which past intrusion(s) may have initially emplaced those stored magmas. This difficulty stems from the tendency for Kīlauea magmas to follow very similar differentiation paths without significant variations in major, minor, or even trace element chemistry. We investigate possible magma sources for the lavas erupted during phase 1a of the 2018 eruption using whole-rock, mineral, and glass major and trace element compositions from historical East Rift Zone eruptions with adjacent fissures. We consider two primary hypotheses for the phase 1a source: magmas associated with the 1955 Lower East Rift Zone eruption or the nine eruptions in the Middle and Upper East Rift Zone during the 1960s. Our results suggest that magma associated with the earliest phases of Kīlauea’s 1955 eruption was the most likely source of the 2018 phase 1a remobilized magma. We determine volatile saturation pressures from melt inclusion chemistry and find similar storage depths for the 2018 phase 1a and early 1955 magmas. The phase 1a and early 1955 lavas are nearly indistinguishable in all of the compositional criteria considered, implying that the leftover 1955 magma body barely cooled and differentiated in the 63 years between eruptions (cooling rates of ~0.1 °C/year). This study sheds light on the potential for protracted storage of eruptible magmas in rift zones at Kīlauea, and highlights some of the challenges and solutions to identifying genetic relationships between magmas at Kīlauea.

Hawaii↗