SEARCH · USGS Science

Results for “Geological Society of America Bulletin”

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

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,666 records · Page 62Linked to original sources

USGS2025

Ultralong, supershear rupture of the 2025 Mw 7.7 Mandalay earthquake reveals unaccounted risk

The 28 March 2025 moment magnitude ( M w ) 7.7 earthquake in Mandalay, Burma (Myanmar), ruptured 475 kilometers of the Sagaing Fault, which was more than twice the length predicted by magnitude scaling relationships. Kinematic slip models and observation of a Rayleigh Mach wave that passed through parts of Thailand confirmed that rupture occurred at supershear velocities of greater than 5 kilometers per second. The anomalous length exposed a vast population to violent near-fault shaking. The Mandalay earthquake is a modern analog for the M w 7.9 1906 San Francisco earthquake, another atypically long and fast rupture. Probabilistic seismic hazard analyses use scaling relations that do not account for such long ruptures at moderate magnitudes. This limitation, in conjunction with a likely increased population and infrastructure exposure for atypically long ruptures, contributes to a potential mischaracterization of seismic risk.

Mandalay↗

USGS2025

Ubiquity and causes of soil water preferential flow across 17 ecoregions

Preferential flow (PF) in soil causes the rapid transport of water, nutrients, and contaminants into the subsurface, influencing groundwater recharge and streamflow. Data scarcity has hindered the quantification of PF occurrence and the identification of its drivers across diverse ecoregions. We address this gap by analyzing high-frequency, multi-depth soil moisture data across 17 ecoregions in the USA, using ∼1,500 sensors at 40 sites. We discovered that PF is widespread, with sites experiencing PF in up to 60% of rainfall events ≥2 mm. Multiple approaches consistently show that PF is more likely to occur with increased peak rainfall intensity, finer textured material, low soil moisture variability, humid climate, and higher net primary productivity. This suggests that PF patterns could shift with projected climate changes, increasing uncertainty in predictions of groundwater recharge, water quality, and streamflow generation.

conterminous United States and Puerto Rico↗

USGS2026

Assessing environmental drivers of denitrification in restored riverine floodplains

Restoration of impaired floodplains is an increasingly prevalent strategy for alleviating water quality concerns and reducing downstream flooding at watershed scales. Floodplains temporarily store water and slow flow velocity to promote sedimentation during overbank flooding and remove inorganic nitrogen from floodwater and groundwater via denitrification. Evaluating the impacts of different restoration strategies on denitrification can inform more strategic investments into floodplain modifications that improve water quality outcomes. Our research investigates how denitrification rates in floodplains respond to environmental factors that are actionable from an engineering perspective through design and water resources management. We seasonally measured soil denitrification enzyme activity and various environmental characteristics in 4 floodplains with different restoration design and management approaches at the confluence of the Wabash and Tippecanoe Rivers in Indiana, United States. Our results showed that denitrification rates in an agricultural floodplain were significantly lower than in restored floodplains with native vegetation. Certain soil conditions characteristic of floodplain wetlands were associated with higher denitrification, particularly elevated total nitrogen, moisture, silt, and organic matter contents. Vegetation species composition was correlated with denitrification rates. This link may reflect the direct effects of vegetation on soil conditions, such as supplying labile organic carbon, or indirect effects, such as vegetation acting as an indicator of hydrologic regime and land use. Denitrification seasonally varied, peaking in winter when nitrate supply from rivers draining agricultural watersheds in the region is also high. Substrate limitation of soil denitrification enzyme activity was most significant during the summer when overbank flooding, which replenishes soil nitrogen stocks, rarely occurs. Our findings indicate that denitrification capacity will likely be maximized in riverine floodplains that are restored as wetlands with diverse native vegetation and enhanced hydrologic connectivity. Such restoration activities promote higher denitrification rates via elevated moisture, fine sediment deposition, and soil organic matter.

Tippecanoe River, Wabash River↗

USGS2024

Afterslip and creep in the rate-dependent framework: Joint inversion of borehole strain and GNSS displacements for the Mw 7.1 Ridgecrest earthquake

The elusive transition toward afterslip following an earthquake is challenging to capture with typical data resolution limits. A dense geodetic network recorded the Mw 7.1 Ridgecrest earthquake, including 16 Global Navigation Satellite System (GNSS) stations and 3 borehole strainmeters (BSM). The sub-nanostrain precision and sub-second sampling rate of BSMs bridges a gap between conventional seismologic and geodetic methods, exemplified by atypical postseismic shear strain reversals observed at nearfield (<2 km) station B921 that remain unexplained. We jointly invert GNSS displacements and BSM strains for coseismic and postseismic slip spanning hours to months over 7 independent periods. Cosiesmically, our model resolves the largest slip magnitudes of up to 6.6 m on the mainshock rupture plane, with similar patterns to other inferred slip distributions. The foreshock fault appears to slip coincidently with mainshock, revealing potential asperities activated during the preceding Mw 6.4 event. Postseismically, the best-fitting models adhere to mechanical rate-and-state expectations of logarithmically decaying slip adjacent to the coseismic rupture terminus, and where deep rheologic conditions favor creep. Most spatial variation occurs in the early postseismic timeframe (<1–2 weeks), with evidence for regional rheologic control and static stress dependence. Triggered creep on the neighboring Garlock Fault unexpectedly persists for >178 days—further highlighting the importance of fault networks in postseismic stress redistribution, critical to assessing future hazard.

Journal of Geophysics Research↗

USGS2025

The feasibility of using lidar-derived digital elevation models for gravity data reduction

Gravity data require submeter elevation accuracy for data processing, and differential global navigation satellite system (dGNSS) equipment is commonly used to acquire three-dimensional positional data to achieve such accuracy. However, lidar (light detection and ranging) data are commonly used to develop digital elevation models (DEMs) of Earth’s surface. Therefore, using elevations from lidar-derived DEMs for gravity-data acquisition and reduction may improve field efficiency and reduce cost. This study examines the feasibility of using DEMs for gravity-data reduction by comparing dGNSS elevation data from 435 gravity stations in Michigan, Wyoming, and Colorado with their respective DEM elevations. The results show that the average difference between DEM and dGNSS elevations is 13 centimeters (cm) and that 93 percent of those differences are less than 50 cm, even in areas with steep terrain. Because an elevation discrepancy of 50 cm corresponds to an error of roughly 0.1 milligals (mGal) in the simple Bouguer gravity anomaly, the results suggest that lidar-derived DEMs are a viable source for acquiring the elevation data needed to process gravity data, thus improving both the cost and efficiency of data collection for regional surveys where an accuracy of less than 1.0 mGal is desired.

Open-File Report↗

USGS2025

Metal fingerprints of Eocene rhyolite magmas coincident with Carlin-type gold deposition in Nevada USA

Eocene magmatic systems contemporaneous with world-class Carlin-type Au deposits in Nevada (USA) have been proposed by some researchers as a key ingredient for Au mineralization, though evidence conclusively demonstrating their genetic relationship remains tenuous. This study provides the first direct evidence of the pre-eruptive metal budget of volatile- and metal-charged silicic magmas coincident in time (~41 to 34 Ma) and space (within 5 km) with Carlin-type Au deposits. We characterize the pre-eruptive metal fingerprints of these diverse magmatic systems to assess their potential as sources of metals for Carlin-type Au mineralization. Metal abundances from quartz-hosted melt inclusions (Au, Te, Ag, Sb, Tl, Mo, W, Sn, As, Pb, Co, Cu, Ni, and Zn) characterized in situ by SHRIMP-RG and LA-ICP-MS represent our best (and only) estimates for the pre-eruptive metal budget in these systems. Median metal concentrations are generally within one order of magnitude of average upper crust and average continental rhyolite values. But there are two notable exceptions, with median Au contents extending >1 order of magnitude higher than average upper crust and median Cu contents ranging >1 order of magnitude lower than upper crust. Despite this, melts contain lower Au/Cu (<0.1), Au/Ag (<5), and Au/Tl (<0.3) than most ore-grade Carlin-type rock samples and quartz-hosted fluid inclusions, regardless of their age and timing relative to nearby Carlin-type Au mineralization. The metal fingerprints of these magmatic systems, de-fined both by traditional and multivariate compositional data analysis techniques, are distinct from one another. Yet none are particularly specialized, e.g., high Au/Cu, in terms of being ideal ingredients as postulated by magmatic models for Carlin-type Au mineralization. Magmatic Au contents do not appear to be correlated with rhyolite “flavors” in the way that Cu, Sn, and Nb contents are. Fluid/melt partitioning modeling and magma volume estimates support the idea that a diverse array of non-specialized silicic magmas could feasibly contribute some or potentially all of the Au, Ag, and Cu in Carlin-type systems. The compositional diversity among contemporaneous magmatic systems could possibly contribute to some of the diversity observed across Carlin-type Au districts in Nevada.

Nevada↗

USGS2026

Hydrothermal manganese-oxide mineralization of a carbonate ooze, Samoan hotspot region, South Pacific Ocean

Low-temperature hydrothermal manganese oxides occur throughout the global oceans. However, the hydrothermal replacement of a carbonate ooze by manganese oxides is described here for the first time. The 24 samples dredged from three locations around the Territory of American Samoa in the South Pacific Ocean include six samples with remnant carbonate and volcaniclastic sediments, which we refer to as “low Mn,” and 18 samples in which the mineralization is pervasive and has replaced most original sediment, termed “high Mn.” The 18 high-Mn samples exhibit a mean Mn content of 51 wt.%. Higher Li and Mn contents and lower Fe contents in the high-Mn samples indicate a hydrothermal origin and distinguish these samples from hydrogenetic ferromanganese crusts. Mn-oxide layers are up to 90 mm thick, some with columns to 44 mm long and 10 mm wide, the magnitude of which has not been described previously. Samples are composed of birnessite and 10 Å phyllomanganate minerals. Textures of the thickest Mn layers indicate mineralization below the seabed from ascending fluids during multiple hydrothermal pulses. Mineralization took place by complete to partial replacement and cementation of foraminiferal sediments intermixed with volcaniclastic sediments in varying amounts. Our results highlight the production of carbonate sediment-hosted hydrothermal Mn oxides from multiple hydrothermal sources within the Samoan volcanic chain. The potential extensive distribution on the regional scale of this newly described mineralization process and unique element enrichments raise questions about its broader distribution globally and potential importance to hydrothermal processes, element mass balance, and seabed mineral resources.

Geochemistry, Geophysics, Geosystems↗

USGS2026

Hazard potential of compound flooding from rainfall, storm surge, and groundwater in coastal New York and Connecticut

Compound flood events, the co-occurrence of multiple flood drivers, can result in flood hazard potential exceeding that of any single driver alone. To evaluate compound flooding in a semi-urbanized coastal area, historical records dating back to 1970 are used to study the co-occurrences of high precipitation, storm surge, and shallow groundwater conditions along the coastlines of New York and Connecticut. Joint return periods for coincident precipitation-surge events were computed using statistical dependence models and compared to the assumption of independence as a ratio, referred to here as a return period adjustment. Results indicate distinct seasonality where compound events in the area disproportionately occur in the cold season between October and April. Return period adjustments range from a factor of 1 to almost 9, demonstrating the range in precipitation-storm surge dependence across the study area. Across all 24 station triad locations, groundwater levels were elevated during times of precipitation- surge co-occurrence, reflecting the tendency for coastal storms and shallow groundwater conditions to co-occur seasonally. The result is a pseudo-trivariate compound flood hazard score and corresponding hazard map that integrates dependence between daily precipitation-surge events and overall monthly groundwater levels (as a precondition) into a relative compound hazard score. The location with the highest compound flood hazard score is on the south shore of Long Island, as well as locations across coastal Connecticut where groundwater levels compound the co-occurrence of heavy precipitation and storm surge.

Connecticut, New York↗

USGS2024

Indications of preferential groundwater seepage feeding northern peatland pools

Groundwater seepage from underlying permeable glacial sedimentary structures, such as eskers, has been hypothesized to directly feed pools in northern peat bogs. These hypotheses directly contradict classical peat bog models for ombrogenous systems, wherein meteoric water is the sole water input to these systems. Variations in the underlying mineral sediment in contact with the peat imply that unrecognized hydrogeologic connectivity may exist with pools in northern peat bogs, particularly where high permeability materials are in contact with the peat. Seepage dynamics originating from these structural variations were investigated using a suite of thermal and hydrogeophysical methods deployed around pools in a peat bog of northeastern Maine, USA. Thermal characterization methods mapped anomalies that were confirmed as matrix seepage or preferential flow pathways (PFPs). Geochemical methods were employed at identified thermal anomalies to confirm upwelling of solute-rich groundwater. Conduits around pools were associated with surficial terminations of suspected peat pipes, based on the inference of pathways extending down into the peat, that focus flow through PFPs in the peat matrix. Discharge also occurred through the peat matrix adjacent to suspected pipe structures and matrix seepage rates were quantified using analysis of diurnal temperature signals recorded at multiple depths. Seepage rates, with a maximum of nearly 0.4 m/d, were measured at localized points around pools. Periods of synchronized temperatures paired with highly muted diurnal temperature signals, recorded in diurnal temperature with depth data, were interpreted qualitatively as activation of strong upward discharge rates through suspected peat pipes. These time periods correlated strongly with local precipitation events around the peatland. Ground-penetrating radar surveys revealed discontinuities in the low permeability glacio-marine clay at the mineral sediment-peat interface, interpreted to be regional glacial esker deposits, which were located beneath and around pools. Heat tracing, specific conductance contrasts, seepage rates, and trace metal concentrations all imply groundwater seepage originating from underlying permeable glacial esker deposits and directly sourcing pools. Preferential groundwater inputs into northern peat bogs may play a key role in developing and maintaining pool systems, with enhanced solute transport impacting peatland ecology, water resources, and carbon cycling.

Maine↗

USGS2025

Examining the compositional selectivity of hydrocarbon oxidation products using liquid–liquid extraction and solid-phase extraction techniques

The effect of extraction methods on detecting hydrocarbon oxidation products (HOPs) in groundwater remains unclear. HOPs are polar, water-soluble byproducts of petroleum biodegradation. Our previous work showed that liquid–liquid extraction (LLE), a method commonly used in regulatory monitoring, has a significantly lower extraction efficiency for HOPs compared to solid-phase extraction (SPE). In this study, we evaluate the analytical limitations and compositional selectivity of LLE and SPE using groundwater samples from the Bemidji, MN, crude oil spill site. Optical properties were characterized using excitation–emission matrix spectroscopy (EEMs), and a three-component PARAFAC model was validated, showing consistent trends across both extracts and whole water samples. Ultrahigh-resolution mass spectrometry (UHR-MS) revealed that LLE selectively recovered aliphatic-like compounds but underrepresented more polar oxygenated HOPs. In contrast, SPE methods were more effective at isolating highly oxidized compound classes. These differences were consistent across a gradient of contamination. Overall, the LLE was less precise and less representative of polar HOPs, introducing bias in the characterization of HOPs. This study is the first to quantitatively demonstrate the compositional selectivity and analytical bias of LLE versus SPE for HOPs using combined EEM-PARAFAC and UHR-MS techniques, with implications for long-term monitoring and site assessment protocols.

Environmental Science and Technology↗

USGS2026

Estimation of parental abundance using hierarchical Bayesian modeling with data augmentation

Pedigree-based estimation methods leverage the fact that each offspring in a cohort is genotypically “marked” by its parents and represent a recent and promising toolset for estimating population dynamics. This includes pedigree accumulation estimators that model the “accumulation” of inferred unique parents within a given cohort to estimate parental abundance. Unlike close-kin mark-recapture approaches, which rely on intercohort comparisons, pedigree accumulation modeling can be completed solely using intracohort samples. This is particularly advantageous for semelparous species, where intercohort pairs are impossible and adult life stages can be difficult to sample without affecting their likelihood of successfully reproducing. Previous work has evaluated a range of estimators for such datasets, concluding that the non-parametric Chao estimator provides the most accurate and precise estimates for feasible levels of sampling effort. We used simulated data to evaluate an alternative estimator based on hierarchical modeling and data augmentation in a Bayesian framework. Results indicate that estimates from the hierarchical Bayesian estimator had comparable accuracy and better precision than both the previously tested Chao1 estimator and the improved iChao formulation across a range of sample sizes and sex ratios. Furthermore, the Bayesian estimator was far more robust to simulated errors in pedigree reconstruction, especially the presence of false negatives. Hierarchical Bayesian pedigree accumulation models can also provide additional insight into underlying reproductive ecology through their use of an explicit observation process, allowing for the incorporation or estimation of species- and population-specific reproductive dynamics. More broadly, the parametric nature of these models offers opportunities to efficiently pool information among datasets as well as to propagate uncertainty within more complex models.

Ecology and Evolution↗

USGS2024

Probabilistic assessment of postfire debris-flow inundation in response to forecast rainfall

Communities downstream of burned steep lands face increases in debris-flow hazards due to fire effects on soil and vegetation. Rapid postfire hazard assessments have traditionally focused on quantifying spatial variations in debris-flow likelihood and volume in response to design rainstorms. However, a methodology that provides estimates of debris-flow inundation downstream of burned areas based on forecast rainfall would provide decision-makers with information that directly addresses the potential for downstream impacts. We introduce a framework that integrates a 24 h lead-time ensemble precipitation forecast with debris-flow likelihood, volume, and runout models to produce probabilistic maps of debris-flow inundation. We applied this framework to simulate debris-flow inundation associated with the 9 January 2018 debris-flow event in Montecito, California, USA. When the observed debris-flow volumes were used to drive the probabilistic forecast model, analysis of the simulated inundation probabilities demonstrates that the model is both reliable and sharp. In the fully predictive model, however, in which debris-flow likelihood and volume were computed from the atmospheric model ensemble's predictions of peak 15 min rainfall intensity, I 15 , the model generally under-forecasted the inundation area. The observed peak I 15 lies in the upper tail of the atmospheric model ensemble spread; thus a large fraction of ensemble members forecast lower I 15 than observed. Using these I 15 values as input to the inundation model resulted in lower-than-observed flow volumes which translated into under-forecasting of the inundation area. Even so, approximately 94 % of the observed inundated area was forecast to have an inundation probability greater than 1 %, demonstrating that the observed extent of inundation was generally captured within the range of outcomes predicted by the model. Sensitivity analyses indicate that debris-flow volume and two parameters associated with debris-flow mobility exert significant influence on inundation predictions, but reducing uncertainty in postfire debris-flow volume predictions will have the largest impact on reducing inundation outcome uncertainty. This study represents a first step toward a near-real-time hazard assessment product that includes probabilistic estimates of debris-flow inundation and provides guidance for future improvements to this and similar model frameworks by identifying key sources of uncertainty.

California↗

USGS2026

Earthquake rupture arrest from depth‐dependent frictional stability on the Pütürge segment of the East Anatolian Fault Zone, Türkiye (Turkey)

Determining why earthquake ruptures stop is a central challenge in earthquake science and seismic hazard assessment. The Pütürge segment of the East Anatolian Fault Zone, Türkiye, exhibits shallow creep (<∼3 km depth) yet hosts M > 6.5 earthquake ruptures at greater depth. Here, we evaluate whether variations in frictional stability along this segment aided arrest of the 2020 M 6.7 Elaziğ and 2023 M 7.8 Pazarcık earthquake ruptures. Analysis of Sentinel‐1 Synthetic Aperture Radar imagery indicates the 2023 M 7.8 rupture propagated laterally into a metamorphic massif within the Pütürge segment, where slip rapidly decayed below detection limits. Creepmeters along this segment recorded no significant surface afterslip (<3 mm) in the following year. To investigate this fault‐slip behavior, we conducted triaxial friction experiments on Pütürge fault gouge sampled from an outcrop exposure. The gouge, composed primarily of muscovite, quartz, and calcite, is velocity strengthening at conditions approximating 0–2.5 km depth and velocity weakening at 4–5 km depth. This transition to velocity‐weakening friction is associated with enhanced comminution and shear localization observed microstructurally. Our results suggest that depth‐dependent frictional stability of the Pütürge fault segment facilitates rupture nucleation and propagation at depth while maintaining rupture‐arresting behavior near Earth’s surface.

Pütürge segment of the East Anatolian Fault Zone↗

USGS2025

Foundational uncertainties in terminal Ediacaran chronostratigraphy revealed by high-precision zircon U-Pb geochronology of the Nama Group, Namibia

The Nama Group of southern Namibia and northwestern South Africa hosts the best-dated mixed carbonate-siliciclastic foreland basin succession of the terminal Ediacaran [ca. 551 million years (Ma) ago to <538 Ma] and is key for resolving the chronology of early metazoan evolution. Numerous silicified volcanic tuff interbeds are present, but differing interpretations regarding the fidelity of their ages lead to different regional stratigraphic correlations, especially for the Urusis Formation of the Schwarzrand Subgroup. An expanded record of the Urusis Formation is found in the Swartpunt area of southern Namibia, which has yielded an important metazoan biota. But the succession in this area is preserved as a series of thrusts at the leading edge of the Gariep orogenic belt and zircon U-Pb data show systematic age repetition. We use regional stratigraphic and structural mapping, integrated with carbonate carbon isotope (δ 13 C carb ) chemostratigraphy and high-precision radioisotope U-Pb zircon geochronology from outcrop and recently acquired drill core to develop a temporally calibrated basin-wide depositional model. This integrated dataset either reflects complex zircon reworking, inheritance, or potential analytical biases (Scenario 1) or the presence of a Gariep-related cryptic décollement within the Spitskop Member that has resulted in stratigraphic repetition (Scenario 2). We investigate the evidence for and against both scenarios and consider their implications for stratigraphic and δ 13 C carb correlations between the Swartpunt area and coeval exposures along the Orange River border with South Africa. Given that these issues are in an area that hosts numerous silicified ash beds and extensive exposure, an inability to confidently discount either scenario highlights a level of compounding uncertainty in zircon U-Pb geochronology that must be considered when attempting to build global chronostratigraphic frameworks. Scenario 1 implies that some of the weighted mean ages and Bayesian eruption ages from the Swartpunt area may be >1 Myr older than the depositional age of their respective ash beds when assuming existing stratigraphic correlations. If this scenario is preferred, then a cautious approach would be to consider all weighted mean zircon U-Pb ages from ash beds to reflect maximum depositional ages. Both scenarios support deposition of the Huns Member >540 Ma in the Swartpunt area if the oldest weighted mean age reported here represents a near-depositional age, which has significant implications for the temporal calibration of important terminal Ediacaran ichnofossil assemblages and future cyclostratigraphic studies. Stratigraphic correlations common to both scenarios allow us to temporally calibrate a basin evolution model for the Nama Group. Temporal trends in initial hafnium isotope (εHf) compositions of zircon grains from ash beds throughout the succession may support progressive crustal thickening associated with underplating of the Damara orogenic belt along the northern periphery of the Kalahari craton from ca. 547 Ma to ca. 538 Ma. The compilation of new and published zircon U-Pb ages may also imply that the locus of carbonate platform development migrated from north to south (present co-ordinates), tracking the migration of foredeep subsidence.

Earth-Science Reviews↗

USGS2026

Structural evolution and slip rate variations through time of the Puente-Hills blind-thrust fault beneath Los Angeles: Implications for seismic hazard and folding kinematics

Using seismic reflection profiles, historical well logging data, and luminescence and radiocarbon ages, we determine a Pleistocene-Holocene slip history for the central, Santa Fe Springs segment of the Puente Hills blind-thrust fault (PHT), a major seismogenic fault situated beneath the urbanized Los Angeles metropolitan region. We analyze the geometry of correlative stratigraphic units in the forelimb and backlimb of the overlying growth-fold, the Santa Fe Springs anticline, and determine the uplift of seven age-correlative markers. Uplift measurements are converted to thrust displacements on the underlying PHT using a structural method laid out by Don et al. (2022, https://doi.org/10.1785/0120220048 ) that accounts for the geometry of the fault. These data indicate that deep thrust displacement on the PHT is partially consumed updip in the creation of a hanging-wall fault-bend fold, with forelimb growth strata recording <80% of the slip documented within the backlimb. Chronological data from growth strata yield age constraints for folding and faulting on the underlying PHT, providing a detailed incremental slip history derived from both the forelimb and backlimb folding for seven discrete growth horizons, spanning the past 1.4 million years. The resulting six incremental slip rates demonstrate that fault slip has varied through time from the middle Pleistocene to the Holocene. Moreover, these results reveal synchronous acceleration of both the central, Santa Fe Springs and western, Los Angeles segments of the fault system since late Pleistocene time (after 200 ka) and slip rates of greater than 2 mm/yr on the downdip, backlimb, fault-ramp below the anticline.

California↗

USGS2026

Effects of human activity on Khumbu Glacier: Towards a sustainable Everest Base Camp

The Everest Base Camp (EBC), a critical staging area for mountaineers attempting to summit Mount Everest, has seen a significant increase in human activity over the past few decades. The increasing number of climbers, trekkers, and support staff, particularly at the campsite area, has intensified environmental impacts on the Khumbu Glacier. This study investigates the effects of human activity at EBC, focusing on three contributors to glacier melt: global climate change, local fossil fuel consumption, and human urinary discharge. Through field data collected during the spring 2023 climbing season, this study examines the energy released through the use of liquefied petroleum gas (LPG), kerosene, and petrol for cooking, heating, and electricity, as well as the heat generated by human urinary discharge. The total heat energy released from these activities at EBC in the spring of 2023 was 849,174 ± 179,774 MJ, which would be sufficient to melt 2492 ± 528 tons of glacier ice and snow. Trends in land surface temperature at EBC and adjacent environments during the 1991–2023 period were calculated from Landsat 5, 7, and 9 satellite data. The 32-year Landsat record reveals that EBC surface temperature increased by 0.28 °C year −1 , which was roughly twice the rate of warming at the surface of Khumbu Glacier adjacent to EBC and about 15% higher than the temperature increase of the debris-covered terrain immediately to the north of EBC. The findings suggest that addressing these anthropogenic influences could help to preserve the Khumbu Glacier and support the sustainability of mountaineering in the region. This paper also takes a transformative approach and explores stable and safer locations in case there is a necessity to relocate the current EBC. Two sites have been identified southwest of the existing EBC, which, unlike the current supraglacial site, are situated on stable ground.

Everest Base Camp, Khumbu Glacier↗

USGS1989

An investigation of shallow ground-water quality near East Fork Poplar Creek, Oak Ridge, Tennessee

Alluvial soils and fill materials in and near the floodplain of East Fork Poplar Creek in Oak Ridge, Tennessee, are contaminated with various trace metals (primarily mercury), organic compounds, and radionuclides that were lost to the stream as a result of past operations at the Y-12 Plant, a nuclear-processing facility located within the U.S. Department of Energy’s Oak Ridge Reservation. Observation wells were installed in the shallow (above bedrock) aquifer in and near the floodplain of East Fork Poplar Creek and water-quality samples were collected to determine if contaminants found in the soils and fill are also present in the shallow ground water. The shallow aquifer in the East Fork Poplar Creek flood plain consists primarily of alluvial silt and clay with lesser amounts of sand and gravel. Thickness of the shallow aquifer ranges from essentially zero to as much as 20 feet. A silty-clay glei horizon is present between the base of the alluvium and the top of bedrock at most flood-plain locations and, where present, likely impedes downward ground-water movement. Water in the shallow aquifer near East Fork Poplar Creek occurs under water-table conditions. Recharge to the shallow aquifer is principally from precipitation, and discharge is through springs and seeps to East Fork Poplar Creek and its tributaries. During spring, summer; and fall, evapotranspiration also accounts for the removal of water in storage in the shallow aquifer. Water levels in the shallow aquifer fluctuate seasonally in response to variations in recharge and evapotranspiration. Generally, the depth to water in the observation wells ranged from about 1 to 4 feet below land surface in late winter, and from about 2 to 7 feet below land surface in late fall. During extremely dry periods, the water table recedes below the top of bedrock in some flood-plain areas, possibly causing East Fork Poplar Creek to lose water to the shallow aquifer along some reaches. Contaminants found in water samples collected from several of the observation wells in concentrations (total and (or) total-recoverable) which equaled or exceeded drinking-water standards established by the U.S. Environmental Protection Agency are antimony, chromium, lead, mercury, selenium, total phenols, and strontium-90. Total and dissolved uranium concentrations exceeded 1.0 microgram per liter in samples from nearly 70 percent of the wells in the East Fork Poplar Creek flood plain. Organic compounds that were identified in low concentrations in samples from a few wells in the flood plain are: Arochlor 1260, benzo(a)anthracene, benzo(b)fluoranthene, chrysene, 3,3-dichlorobenzidine, di-n-butylphthalate, N-nitrosodiphenylamine, and pyrene. Water from one well in the East Fork Poplar Creek flood plain at a contaminated fill site contained 37and 8 micrograms per liter of trichloroethene and trans-1,2-dichloroethene, respectively. Comparison of the results of total and (or) total-recoverable trace-metal determinations with those from dissolved determinations demonstrates that elevated concentrations of these substances in water collected from several of the wells in the East Fork Poplar Creek flood plain resulted from sorption of trace metals (and possibly organic compounds and radionuclides) by fine sediment suspended in the samples. The occurrence of contaminated sediment in these samples is suspected to be the result of borehole contamination during well installation.

Tennessee↗

USGS2025

A glimpse into the future of tectonic tremor monitoring

Tectonic tremor is a weak, long-duration seismic signal often observed in subduction zones and on some other plate-bounding faults. Because of tremor's characteristically low amplitude (and low signal-to-noise) and lack of clear phase arrivals, detecting and locating tremor usually requires techniques distinct from those applied to typical earthquakes. Major advances in detection and understanding of tremor have derived in the past from a powerful combination of new data and new analysis techniques. In a recent study, Sagae et al. (2025, https://doi.org/10.1029/2025jb031348 ) exploit that combination again, developing a new machine-learning based workflow and applying it to the S-net cabled seismic network in the Japan trench offshore northern Honshu. Their approach, although complex, succeeds in detecting several times more tremor activity than earlier studies, resulting in new insights and providing a blueprint for similar approaches that could be applied elsewhere. As real-time earthquake monitoring adopts similar tools, it may present an opportunity to bring tremor monitoring into operational workflows. In turn, this could solidify tremor monitoring as a component of future operational earthquake forecasting.

JGR Solid Earth↗