USGS Science⌕ Search

SEARCH · USGS Science

Results for “Hydrology and Earth System Sciences (HESS)”

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,661 records · Page 21Linked to original sources
USGS2026

Fully contained laboratory earthquakes: The effect of asperity aspect ratio and free surfaces

Corner frequency ( f c ) and seismic moment ( M 0 ) are key parameters derived from seismic signals that are used to characterize earthquake stress drop, rupture area, and slip. These parameters are also affected by fault geometry and boundary conditions. However, the systematic study of these effects in laboratory settings has been challenging. This study presents laboratory earthquake experiments that examine how rupture dynamics are influenced by (a) the aspect ratio of rectangular PMMA velocity-weakening (VW) asperities surrounded by the Teflon velocity-strengthening (VS) patches, and (b) whether the sides of a VW asperity are confined with VS patches or are free surfaces. We found that increasing confinement by reducing free surfaces or increasing the VW asperity aspect ratio stabilizes fault slip, so that higher normal stress is required to transition from aseismic to seismic slip. Increased confinement and high aspect ratios also reduced M 0 and increased f c , both of which were determined from the radiated seismic waves. M 0 and f c were primarily controlled by the shorter dimension of the VW asperity. Analysis of high-frequency acoustic emission signals revealed that ruptures on high-aspect-ratio VW asperities propagated more unidirectionally, whereas ruptures on square VW asperities were more complex. Further, the high-aspect-ratio asperities were more likely to be eroded by surrounding VS regions while low-aspect-ratio asperities were more likely to rupture into the VS surroundings. These results demonstrate that both the confinement from surrounding stable areas and the geometry of the seismogenic patch can affect rupture nucleation, propagation, and seismic source characteristics.

JGR Solid Earth↗
USGS2024

Constraining mean landslide occurrence rates for non-temporal landslide inventories using high-resolution elevation data

Constraining landslide occurrence rates can help to generate landslide hazard models that predict the spatial and temporal occurrence of landslides. However, most landslide inventories do not include any temporal data due to the difficulties of dating landslide deposits. Here we introduce a method for estimating the mean landslide occurrence rate of deep-seated rotational and translational slides derived solely from high-resolution (≤3 m) elevation data and globally available estimates of the diffusion coefficient for sediment flux. The method applies a linear diffusion model to the roughest landslide deposits until they reach a representative non-landslide roughness distribution. This estimates the time for a landslide deposit to be unrecognizable in high-resolution digital elevation data, which we term the mean lifetime of the landslide. Using the mean lifetime and number of landslides within an area of interest, we can estimate the mean occurrence rate of landslides over that domain. We validate this approach using a comprehensive temporal inventory of landslides in western Oregon created using age-roughness curves that are calibrated with high-resolution elevation data and radiocarbon data. We find good agreement between our diffusion method and the existing age-roughness-derived estimates, producing mean lifetimes of 4500 and 5200 years (4% difference), respectively. Hazard maps produced using the two methodologies generally agree, with the maximum differences in landslide probability reaching 0.1. Due to the relative abundance of high-resolution elevation data compared with age-dated landslides, our method could help constrain landslide occurrence rates in areas previously considered unfeasible.

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

Three-dimensional seismic velocity models for the San Francisco Bay region, California from joint body-wave and surface-wave tomography validated by waveform simulation

A high-resolution seismic velocity model for both P and S waves is essential for the San Francisco Bay (SFB) region to accurately simulate earthquake ground motions and assess seismic hazards in this high-risk region. We present new three-dimensional P-wave and S-wave velocity (Vp and Vs) models of the SFB region developed using joint tomographic inversion with a rich set of earthquake and active-source body-wave data and ambient-noise surface-wave data. We designed a dynamic balance strategy to effectively balance the contributions of body-wave and surface-wave data during each iteration of the inversion. Owing to the complementary resolution of body-wave and surface-wave data, our new models are significantly improved compared to previous tomographic models. Major geologic features in the SFB region, including low-velocity basins, high-velocity Franciscan Complex, granite, and gabbro rocks, and across-fault velocity contrasts, are imaged in our Vp and Vs models. Compared to the geology-based U.S. Geological Survey velocity model, our new model has better defined structures for some important areas that are of great concern for seismic hazard assessment, such as the Sacramento-San Joaquin Delta, and the 2014 magnitude 6.0 South Napa earthquake area. We validated our model by simulating waveforms for a set of moderate-magnitude earthquakes. The good performance of our model in waveform simulation indicates its potential for improving ground motion simulation and seismic hazard assessment in the SFB region.

California↗
USGS2025

To heal or not to heal?: 1. The effect of pore fluid pressure on the frictional healing behavior of lithologies in Oklahoma

The competition between fault healing (i.e., re-strengthening) and fault loading determines the timing and magnitude of fault failure within the seismic cycle. Repeating earthquakes can give observational estimates of fault healing rates, however, it is difficult to link laboratory studies of frictional healing and observed healing rates from repeating earthquakes in part because of uncertainty in lithology at depth. Due to well-constrained and relatively simple geology, earthquakes in Oklahoma can be linked to the granitic basement rock and to the Arbuckle Group, which is primarily composed of dolomite at earthquake depths. Here, we conduct friction experiments to measure healing rates of the two earthquake-bearing lithologies at confining pressures representative of earthquake depths and pore pressures ranging from 0% to 80% of the confining pressure. We measure frictional healing by executing slide-hold-slide tests with hold times ranging from 3 s to 3000 s. The friction experiments on the Troy Granite indicate that pore fluid pressure does not greatly affect healing rate. On the other hand, the dolomite of the Arbuckle Group exhibits decreased healing with increased pore fluid pressure, with weakening at the highest pore pressure. We hypothesize that this is due to an increase in dissolution of dolomite at high pore pressures/low effective normal stress. These healing rates are used in the companion paper to understand the moment-recurrence time behavior of repeating earthquakes in Prague, Oklahoma. This work has implications for possible enhanced dissolution and weakening behavior of the Arbuckle Group during wastewater injection activities.

Oklahoma↗
USGS2025

Temporal and spatial comparison of coal mine ventilation methane emissions and mitigation quantified using PRISMA satellite data and on-site measurements

Emission monitoring at the facility level (bottom-up, BU) is key for accurate reporting of coal mine methane (CMM) emissions. Recent advances in aerial and satellite observations (top-down, TD) indicate that these methods have the potential to support CMM emissions monitoring and reporting of, as well as track the effectiveness of, mitigation actions. That said, studies have shown discrepancies between BU and TD estimations. Performing TD monitoring with concurrent BU measurements at the same mines may help address the observed discrepancies and improve quantification methods to narrow the gap between BU and TD data. This paper presents a comparison of concurrent BU-TD methane emission fluxes monitored from two ventilation shafts at a southwestern Virginia, USA, longwall mine to complement the existing body of studies on satellite-based monitoring of coal mines by incorporating continuous ground-based monitoring with concurrent TD monitoring of methane emission fluxes to address the gap and provide valuable insights into temporal emission patterns. The shafts were monitored on multiple dates between 2020 and 2023. BU monitoring was performed at the exhaust fans, while TD quantifications used PRISMA hyperspectral satellite data and two different wind reanalysis datasets (i.e., GEOS-FP and ERA5). This mine also offered a unique opportunity for BU-TD comparisons before and after ventilation air methane (VAM) oxidizer operation, which was installed at one of the shafts in 2022. The results showed that TD-estimated mean fluxes were generally lower than BU data, which were attributed to quantification difficulties associated with the low albedo caused by heavy vegetation and the terrain of the area. However, despite the discrepancies in mean emissions, both the interquartile range and the data range of the distributions generally overlapped, and the estimates correctly showed the emission trends.

Virginia↗
USGS2024

Deformation by pressure solution and grain boundary sliding in a retrograde shear zone in southern New England, USA

Alleghanian phyllonites in a shear zone in southern New England were formed by the retrogression and hydration of a high-grade Acadian pelitic schist. The retrogression was locally incomplete, resulting in both heterogeneous mineralogy and mineral compositions, and in many arrested reaction textures. These features, documented by backscattered electron and element map imaging, and by 40 Ar/ 39 Ar dating, provide a record of the processes active as the schist was converted to a phyllonite. On the margins of the shear zone relic biotite and garnet grains survive but are partially replaced by chlorite. Relic high-grade muscovite flakes containing up to 35 % paragonite (pg) also persist but are entrained in swarms of pg-poor muscovite flakes producing folia with an anastomosing network of muscovite grains of variable composition. Single crystals of pg-rich muscovite are truncated by muscovite in successively younger cross-cutting folia that contain decreasing pg contents, falling to <5% pg in the youngest folia. In the core of the shear zone all high-grade minerals are destroyed, but recrystallization continues with earlier greenschist facies folia replaced by lower grade and pg-poorer muscovite in younger folia. Muscovite grains in truncated folia are commonly kinked and selective recrystallization of muscovite in kink bands to pg-poor-compositions demonstrates that strain energy helped drive recrystallization. The boundaries of these new pg-poor muscovite folia are sharp, and truncations indicate that grain boundary sliding (GBS) involved dissolution. GBS also enabled folding by flexural slip along muscovite grain boundaries. Fold amplitudes increased as chlorite dissolved from fold limbs precipitated in fold hinges producing crescent-shaped saddle-reef-like structures. Together these observations of truncation and replacement by dissolution and precipitation demonstrate that pressure solution and GBS facilitated both the formation of the phyllonitic shear zone and slip along it.

Connecticut↗
USGS2025

Relationship of atmospheric nitrogen deposition to soil nitrogen cycling along an elevation gradient in the Colorado Front Range

Microbial processing of atmospheric nitrogen (N) deposition regulates the retention and mobilization of N in soils, with important implications for water quality. Understanding the links between N deposition, microbial communities, N transformations, and water quality is critical as N deposition shifts toward reduced N and remains persistently high in many regions. Here, we investigated these connections along an elevation transect in the Colorado Front Range. Although rates of N deposition and pools of extractable N increased down the elevation transect, soil microbial communities and N transformation rates did not follow clear elevational patterns. The subalpine microbial community was distinct, corresponding to a high C:N ratio and low pH, while the microbial communities at the lower elevation sites were all very similar. Net nitrification, mineralization, and nitrification potential rates were highest at the Plains (1,700 m) and Montane (2,527 m) sites, suggesting that these ecosystems mobilize N. In contrast, the net immobilization of N observed at the Foothills (1,978 m) and Subalpine (3,015 m) sites suggests that these ecosystems retain N deposition. The contrast in N transformation rates between the plains and foothills, both of which receive elevated N deposition, may be due to spatial heterogeneity not captured in this study and warrants further investigation. Stream N concentrations from the subalpine to the foothills were consistently low, indicating that these soils are currently able to process and retain N deposition, but this may be disrupted if drought, wildfire, or land-use change alter the ability of the soils to retain N.

Colorado↗
USGS1998

Peak-discharge frequency and potential extreme peak discharge for natural streams in the Brazos River basin, Texas

The 2-, 5-, 10-, 25-, 50-, and 100-year peak discharges were estimated for 186 streamflow-gaging stations with at least 8 years of data for natural streams in and near the Brazos River Basin, Texas. Multiple regression equations were developed to estimate peak-discharge frequency for the 2-, 5-, 10-, 25-, 50-, and 100-year recurrence intervals for each of three hydrologic regions that compose the Brazos River Basin. The equations for each region are a function of significant basin characteristics (explanatory variables). The significant explanatory variables among six that were tested are the contributing drainage area and stream slope for regions 1 and 2 and the contributing drainage area for region 3. For the three sets of equations, the coefficient of determination ranges from 0.59 to 0.93, and the standard error ranges from 0.184 to 0.391 log units. A larger coefficient of determination and a lower standard error generally are associated with the equations for hydrologic regions 2 and 3. Statistics from the regression analysis allow computation of the prediction interval associated with a given significance level for a peak-discharge frequency estimate. The regression equations can be used to estimate peak discharges for sites at, near, or away from sites with streamflow-gaging stations. The potential extreme peak-discharge curves as related to contributing drainage area were estimated for each of the three hydrologic regions from measured extreme peaks of record at 186 sites with streamflow-gaging stations and from measured extreme peaks at 37 sites without streamflow-gaging stations in and near the Brazos River Basin. The potential extreme peak-discharge curves generally are similar for hydrologic regions 1 and 2, and the curve for region 3 consistently is below the curves for regions 1 and 2, which indicates smaller peak discharges.

Texas↗
USGS2024

U.S. Geological Survey Karst Interest Group proceedings, Nashville, Tennessee, October 22-24, 2024

Karst hydrogeologic systems represent challenging and unique conditions to scientists studying groundwater flow and contaminant transport. Karst terrains are characterized by distinct and beautiful landscapes, caverns, and springs, and many of the exceptional karst areas are designated as national or state parks. The range and complexity of landforms and groundwater flow systems associated with karst terrains are enormous, perhaps more than any other aquifer type. The U.S. Geological Survey (USGS) Karst Interest Group (KIG), formed in 2000, is a loosely knit, grassroots organization of USGS and non-USGS scientists and researchers devoted to fostering better communication among scientists working on, or interested in, karst aquifers. The primary mission of the KIG is to encourage and support interdisciplinary collaboration and technology transfer among scientists working in karst areas. To accomplish its mission, the KIG has organized a series of workshops. To date (2024), nine KIG workshops, including the workshop documented in this report, have been held. The abstracts and extended abstracts provide a snapshot in time of past and current karst related studies. The USGS Water Availability and Use Science Program funded the workshop and proceedings. The planning committee for the ninth workshop includes Thomas D. Byl (USGS and Tennessee State University), Allan K. Clark (USGS), Laura M. DeMott (USGS), Eve L. Kuniansky (USGS, Emeritus), Benjamin V. Miller (USGS), and Lawrence E. Spangler (USGS, Emeritus). The workshop proceedings are edited by Eve L. Kuniansky and Lawrence E. Spangler. The field trip guide was produced by Benjamin V. Miller and Brian Ham (Tennessee Department of Environment and Conservation) and included in the proceedings from the KIG’s 2021 virtual workshop to be used on the optional field trip held on Thursday, October 24, 2024.

Open-File Report↗
USGS2024

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

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

New York, Virginia↗
USGS2024

Subduction zone geometry modulates the megathrust earthquake cycle: Magnitude, recurrence, and variability

Megathrust geometric properties exhibit some of the strongest correlations with maximum earthquake magnitude in global surveys of large subduction zone earthquakes, but the mechanisms through which fault geometry influences subduction earthquake cycle dynamics remain unresolved. Here, we develop 39 models of sequences of earthquakes and aseismic slip (SEAS) on variably-dipping planar and variably-curved nonplanar megathrusts using the volumetric, high-order accurate code tandem to account for fault curvature. We vary the dip, downdip curvature and width of the seismogenic zone to examine how slab geometry mechanically influences megathrust seismic cycles, including the size, variability, and interevent timing of earthquakes. Dip and curvature control characteristic slip styles primarily through their influence on seismogenic zone width: wider seismogenic zones allow shallowly-dipping megathrusts to host larger earthquakes than steeply-dipping ones. Under elevated pore pressure and less strongly velocity-weakening friction, all modeled fault geometries host uniform periodic ruptures. In contrast, shallowly-dipping and sharply-curved megathrusts host multi-period supercycles of slow-to-fast, small-to-large slip events under higher effective stresses and more strongly velocity-weakening friction. We discuss how subduction zones' maximum earthquake magnitudes may be primarily controlled by the dip and dimensions of the seismogenic zone, while second-order effects from structurally-derived mechanical heterogeneity modulate the recurrence frequency and timing of these events. Our results suggest that enhanced co- and interseismic strength and stress variability along the megathrust, such as induced near areas of high or heterogeneous fault curvature, limits how frequently large ruptures occur and may explain curved faults' tendency to host more frequent, smaller earthquakes than flat faults.

Journal of Geophysical Research: Solid Earth↗
USGS2026

Evaluating groundwater quality influences from oil field operations and other anthropogenic activities in an urban setting, Santa Fe Springs, California

Groundwater quality is often affected by anthropogenic activities in urban settings. This study examines groundwater quality in and around the Santa Fe Springs Oil Field in Los Angeles County, California, where oil and gas production commonly intersects with high density industrial, commercial and residential land uses. Utilizing a combination of new and historical data, we evaluated potential pathways that would allow for oil field formation fluids to migrate into groundwater and whether mixing may have occurred based on the distribution of groundwater and oil field formation fluid tracers in samples. Samples were analyzed for a wide array of constituents including volatile organic compounds, light hydrocarbons, major ions, and various isotopic compositions. Despite evidence of oil field infrastructure providing potential pathways of migration via uncemented annular spaces, casing breaches and historical disposal of oil field formation water in surface ponds, the distribution and occurrence of stable isotopes of water, chloride, boron, and total dissolved solids do not indicate mixing of oil field formation water and groundwater. However, methane isotopic signatures and the presence of heavier alkanes suggest gas from oil-bearing formations have migrated from depth via oil field well infrastructure. Volatile organic compound detections were mainly from manufactured compounds unrelated to oil and gas production, with a relatively limited number of petroleum hydrocarbons also detected. Volatile organic compounds were generally found in wells tapping shallow, modern aged groundwater, indicating anthropogenic activities occurring at or near land surface as the source. Study results suggest that while oil field infrastructure provides migration pathways for oil field formation fluids to be introduced into groundwater, urban land uses not related to oil and gas production are the primary drivers of groundwater quality degradation.

California↗
USGS2026

Channel morphology and large wood control postfire debris-flow erosion and deposition

Runoff-generated debris flows are a known response to wildfire, and accurately predicting the volume of these debris flows is important for estimating the magnitude of downstream hazards. Prior data collection efforts have focused on debris-flow volume measurements at catchment outlets, but few studies have considered how erosion and deposition modulate the volume of sediment arriving at catchment outlets. This study takes advantage of a high-resolution dataset to document the factors that control the total debris-flow volume reaching the catchment outlet during a fatal postfire debris flow. Using pre- and post-event airborne lidar, satellite imagery and field mapping, we found that a postfire debris flow in the Black Hollow catchment in northern Colorado eroded 136,000 ± 30,000 m 3 and redeposited 27,000 ± 7,500 m 3 in the main channel. Most of the in-channel deposition (52% by volume) occurred where a confined channel reach transitioned to an unconfined channel reach downstream, allowing the flow to widen and deposit material. Wood jams played multiple roles in the debris-flow dynamics, both nucleating deposition (25% of the deposit volume was stored behind wood jams) and exacerbating erosion (50% of the total erosion occurred downstream from a wood dam break). The remaining deposition occurred due to spatial changes in channel slope as well as deposition observed at newly formed channel bars. Using these data in this study, we identified topographic and vegetation metrics that can be used (pre-event) to anticipate where deposition may occur in channels prior to a debris flow.

Colorado↗
USGS1998

Characterization of hydrogeologic units using matrix properties, Yucca Mountain, Nevada

Determination of the suitability of Yucca Mountain, in southern Nevada, as a geologic repository for high-level radioactive waste requires the use of numerical flow and transport models. Input for these models includes parameters that describe hydrologic properties and the initial and boundary conditions for all rock materials within the unsaturated zone, as well as some of the upper rocks in the saturated zone. There are 30 hydrogeologic units in the unsaturated zone, and each unit is defined by limited ranges where a discrete volume of rock contains similar hydrogeologic properties. These hydrogeologic units can be easily located in space by using three-dimensional lithostratigraphic models based on relation- ships of the properties with the lithostratigraphy. Physical properties of bulk density, porosity, and particle density; flow properties of saturated hydraulic conductivity and moisture-retention characteristics; and the state variables (variables describing the current state of field conditions) of saturation and water potential were determined for each unit. Units were defined using (1) a data base developed from 4,892 rock samples collected from the coring of 23 shallow and 8 deep boreholes, (2) described lithostratigraphic boundaries and corresponding relations to porosity, (3) recognition of transition zones with pronounced changes in properties over short vertical distances, (4) characterization of the influence of mineral alteration on hydrologic properties such as permeability and moisture-retention characteristics, and (5) a statistical analysis to evaluate where boundaries should be adjusted to minimize the variance within layers. This study describes the correlation of hydrologic properties to porosity, a property that is well related to the lithostratigraphy and depositional and cooling history of the volcanic deposits and can, therefore, be modeled to be distributed laterally. Parameters of the hydrogeologic units developed in this study and the relation of flow properties to porosity that are described can be used to produce detailed and accurate representations of the core-scale hydrologic processes ongoing at Yucca Mountain.

Nevada↗
USGS2025

Rainfall thresholds for postfire debris-flow initiation vary with short-duration rainfall climatology

The size, frequency, and geographic scope of severe wildfires are expanding across the globe, including in the Western United States. Recently burned steeplands have an increased likelihood of debris flows, which pose hazards to downstream communities. The conditions for postfire debris-flow initiation are commonly expressed as rainfall intensity-duration thresholds, which can be estimated given sufficient observational history. However, the spread of wildfire across diverse climates poses a challenge for accurate threshold prediction in areas with limited observations. Studies of mass-movement processes in unburned areas indicate that thresholds vary with local climate, such that higher rainfall rates are required for initiation in climates characterized by frequent intense rainfall. Here, we use three independent methods to test whether initiation of postfire runoff-generated debris flows across the Western United States varies similarly with climate. Through the compilation of observed thresholds at various fires, analysis of the spatial density of observed debris flows, and quantification of feature importance at different spatial scales, we show that postfire debris-flow initiation thresholds vary systematically with short-duration rainfall-intensity climatology. The predictive power of climatological data sets that are readily available before a fire occurs offers a much-needed tool for hazard management in regions that are facing increased wildfire activity, have sparse observational history, and/or have limited resources for field-based hazard assessment. Furthermore, if the observed variation in thresholds reflects long-term adjustment of the landscape to local climate, rapid shifts in rainfall intensity related to climate change will likely induce spatially variable shifts in postfire debris-flow likelihood.

Arizona, California, Colorado, Nevada, New Mexico,↗
USGS2024

Produced water geochemistry from hydraulically stimulated Niobrara Formation petroleum wells: Origin of salinity and temporal perspectives on treatment and reuse

Produced water (i.e., a mixture of returned injection fluids and geologic formation brines) represents the largest volumetric waste stream associated with petroleum production in the United States. As such, produced water has been the focus of intense study with emphasis on understanding the geologic origin of the fluids, environmental impacts of unintended or intentional release, disposal concerns, and their commodity (e.g., lithium) potential. However, produced water geochemistry from many active petroleum plays remain poorly constrained leading to knowledge gaps associated with the origin of brine salinity and parameters (e.g., radium levels) that can impact treatment, disposal, and possible reuse. Here we evaluate the major ion geochemistry, radium concentrations, and stable water isotope composition of ~120 produced water samples collected from 17 producing unconventional petroleum wells in Weld County, Colorado from the Late Cretaceous Niobrara Formation. This sample set encompasses eight produced water time series from four new wells across production days 0 to ~365 and from four established wells across production days ~1000 to ~1700. Additionally, produced water from nine other established Niobrara Formation wells were sampled at discrete time points ranging from day 458 to day 2256, as well as hydraulic fracturing input fluids. These results expand the available Niobrara Formation produced water geochemical data, previously limited to few wells sampled within the first year of production, allowing for the heterogeneity of major ions and radium to be evaluated. Furthermore, we explore the geochemical relationships between major ion ratios and stable water isotope composition to understand the origin of salinity in Niobrara Formation brines from the Denver-Julesburg Basin. These findings are discussed with perspective toward potential treatment and reuse of Niobrara produced water prior to disposal.

Colorado, Nebraska, Wyoming↗
USGS2024

Intense alteration on early Mars revealed by high-aluminum rocks at Jezero Crater

The NASA Perseverance rover discovered light-toned float rocks scattered across the surface of Jezero crater that are particularly rich in alumina ( ~ 35 wt% Al 2 O 3 ) and depleted in other major elements (except silica). These unique float rocks have heterogeneous mineralogy ranging from kaolinite/halloysite-bearing in hydrated samples, to spinel-bearing in dehydrated samples also containing a dehydrated Al-rich phase. Here we describe SuperCam and Mastcam-Z observations of the float rocks, including the first in situ identification of kaolinite or halloysite on another planet, and dehydrated phases including spinel and apparent partially dehydroxylated kaolinite. The presence of spinel in these samples is likely detrital in origin, surviving kaolinitization, pointing to an ultramafic origin. However, the association of low hydration with increased Al 2 O 3 abundances suggests heating-induced dehydration which could have occurred during the lithification or impact excavation of these rocks. Given the orbital context of kaolinite-bearing megabreccia in the Jezero crater rim, we propose an origin for these rocks involving intense aqueous alteration of the parent material, followed by dehydration/lithification potentially through impact processes, and dispersion into Jezero crater through flood or impact-related processes.

Communications Earth & Environment↗