USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research: Solid Earth”

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,663 records · Page 63Linked to original sources
USGS2025

Accurate simulation of flow through dipping aquifers with MODFLOW 6 using enhanced cell connectivity

In simulations of groundwater flow through dipping aquifers, layers of model cells are often “deformed” to follow the top and bottom elevations of the aquifers. When this approach is used in MODFLOW, adjacent cells within the same model layer are vertically offset from one another, and the standard conductance-based (two-point) formulation for flow between cells does not rigorously account for these offsets. The XT3D multi-point flow formulation in MODFLOW 6 is designed to account for geometric irregularities in the grid, including vertical offsets, and to provide accurate results for both isotropic and anisotropic groundwater flow. A recent study evaluated the performance of the standard formulation and XT3D using a simple, synthetic benchmark model of a steeply dipping aquifer. Although XT3D generally improved the accuracy of flow simulations relative to the standard formulation as expected, neither formulation produced accurate flows in cases that involved large vertical offsets. In this paper, we explain that the inability of XT3D to produce accurate flows in the steeply dipping aquifer benchmark was not due to an inherent limitation of the flow formulation, but rather to the limited cell connectivity inherent in the most commonly used discretization packages in MODFLOW 6. Furthermore, we demonstrate that XT3D is able to produce the expected accuracy when adequate cell connectivity is introduced using MODFLOW's unstructured grid type and the aquifer is discretized vertically using at least two model layers.

Groundwater↗
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↗
USGS2024

Analyzing spatial distributions and alignments of pitted cone features in Utopia Planitia on Mars

Martian geomorphology and surface features provide links to understanding past geologic processes such as fluid movement, local and regional tectonics, and feature formation mechanisms. Pitted cones are common features in the northern plains basins of Mars. They have been proposed to have formed from upwelling volatile-rich fluids, such as magma or water-sediment slurries. In this study, we map the spatial distributions of pitted cone features across the Utopia Planitia (UP) basin. Using the Average Nearest Neighbor technique, we find that pitted cone features appear to be clustered across the basin, occurring in a narrow band around the circumferential rim of UP. Additionally, we find that pitted cone features also appear to be aligned in chains that are sub-parallel with the basin margin. Parallel bands of cones generally follow elevation contours of the UP basin, which suggests elevation or a correlated factor plays a major role in pitted cone and cone chain formation. We propose that pitted cones and cone chains may be related to vertical fractures formed around the UP basin rim from subsidence of infilling basin material.

Icarus↗
USGS2025

Conceptual and numerical groundwater flow model of the Iowa River alluvial aquifer near Tama County, Iowa, 1980 through 2022

The Iowa River alluvial aquifer is an important source of water on the Meskwaki Settlement in Tama County, Iowa, which is land owned by the Sac & Fox Tribe of the Mississippi in Iowa (commonly known as the Meskwaki Nation). The U.S. Geological Survey constructed a groundwater flow model, including a conceptual and numerical model, of the Iowa River alluvial aquifer and underlying hydrogeologic units near the Meskwaki Settlement in Tama County, Iowa, for the period of January 1980–August 2022 to estimate the fraction of water pumped from the Iowa River alluvial aquifer by Meskwaki Settlement wells that is derived from streamflow depletion in the Iowa River and its tributaries. Streamflow depletion is a reduction in streamflow caused by groundwater pumping and includes the interception by groundwater production wells of water that otherwise would have been discharged to streams (called “captured groundwater discharge”) and induced infiltration of streamflow to the production wells. Calibrated model runs were performed with no simulated pumping and simulated pumping only at Meskwaki Settlement wells, and the change in simulated flow rates between the groundwater system and streams for the two model runs represents the amount of streamflow depletion in the Iowa River and tributary streams resulting from pumping at the Meskwaki Settlement wells. Streamflow depletion in the Iowa River and its tributaries as a percentage of simulated pumping at the Meskwaki Settlement wells was calculated by dividing this difference by the total simulated pumping rate for the Meskwaki Settlement wells. The model results demonstrate that the mean monthly streamflow depletion, including induced infiltration and captured discharge, in the Iowa River and its tributary streams as a percentage of mean monthly pumping at the Meskwaki Settlement wells was 97.4 percent and ranged from 65.4 to 112 percent. Of the total streamflow depletion, mean monthly induced recharge was 20.9 percent and ranged from 4.9 to 37.2 percent. Mean monthly captured discharge was 76.5 percent and ranged from 57.1 to 97.1 percent. These results indicate that most of the water pumped from the Meskwaki Settlement wells is the result of streamflow depletion, in the form of both induced infiltration and captured discharge.

Iowa↗
USGS2026

ShakeAlert Earthquake Early Warning System performance during the Mw 7.0 offshore Cape Mendocino earthquake

The 5 December 2024 M w 7.0 Offshore Cape Mendocino earthquake was a challenging test of the U.S. West Coast ShakeAlert earthquake early warning system due to its offshore epicenter and limited near‐source station coverage. We analyzed real‐time performance of all components of the ShakeAlert system, including the seismic algorithms (earthquake point‐source integrated code [EPIC] and Finite‐fault rupture Detector [FinDer]), the geodetic algorithm (Geodetic First Approximation of Size and Time–peak ground displacement [GFAST‐PGD]), and network telemetry during the event. EPIC created the first solution for this earthquake 15 s after origin time with an initial magnitude estimate of M 5.6 and location error of 10 km from the Advanced National Seismic System epicenter. An early spurious trigger from station CE.89101 fortuitously maintained location accuracy and, correspondingly, magnitude accuracy. FinDer contributed its first solution at 18 s with a location estimate closer to the seismic network and produced two distinct rupture geometries, leading to minor fluctuations in estimated intensity contours. GFAST‐PGD did not meet alerting thresholds but otherwise performed as expected. Network latencies were <2 s for most stations, supporting the rapid detection of this earthquake by the system. Roughly five million alerts were delivered to cell phone devices in California and Oregon during this event. This was also the first instance of a school district‐wide ShakeAlert‐powered system being activated. Comparisons to recorded seismograms demonstrate that the maximum warning times before potentially damaging shaking (intensity 6+) were in the range of 5–55 s. Although the ShakeAlert system provided accurate solutions and useful alert delivery, this earthquake raised awareness of potential issues within the system, including the need for improved offshore location estimates, a combination of solutions from ShakeAlert servers, and handling of spurious triggers.

California↗
USGS2024

Conceptual hydrogeologic framework and groundwater budget near the southeastern part of Puget Sound, Washington

More than 1 million people live within the active model area (AMA) in the southeastern part of the lowlands surrounding Puget Sound, or Puget Lowland, Washington, and groundwater is the source for approximately one-half of their public, domestic, and irrigation water demands. The 887-square-mile AMA, located in King and Pierce Counties, represents the area of analysis for the conceptual hydrogeologic framework and numerical groundwater-flow models within the study area and includes the Puyallup River and Chambers-Clover Creek watersheds. To assess the potential hydrologic and anthropogenic impacts to groundwater and the connected surface-water resources, conceptual and numerical groundwater-flow models of groundwater flow were developed by the U.S. Geological Survey Washington Water Science Center in close cooperation with 18 water-resource agencies and stakeholders. This report presents information used to characterize the groundwater-flow system and the development of a numerical model in the AMA. Included are descriptions of the geology and conceptual hydrogeologic framework, groundwater levels and flow directions, groundwater recharge and discharge, numerical groundwater-flow model construction and results, and model limitations. The study area encompasses the western part of Pierce County and the southwestern part of King County, Washington. The study area extends south to the Nisqually River, southwest to Tanwax Creek, northeast to the Green River, and north through the valley near Auburn and adjacent uplands. It is bounded on the east by foothills of the Cascade Range, and on the northwest by Puget Sound.

Washington↗
USGS2026

Biome-scale spatial patterns of avian abundance reveal proactive conservation opportunities in North American grasslands

North American grassland birds have experienced steeper population declines than any other avian guild, yet conservation efforts remain largely reactive and fragmented. We used nearly four decades of North American Breeding Bird Survey data to identify biome-scale spatial patterns (clustering) of grassland bird abundance for the Great Plains. Our results reveal an ecological core in the north-central Plains where community-level abundance is either increasing by >100% or remains high and stable, providing a strategic roadmap for a “Defend the Core” conservation approach. This approach flips the script from reactive triage centered on isolated population fragments to a proactive strategy of maintaining large-scale ecosystem integrity. Conversely, we found that population losses are more spatially clustered than wins, reflecting the relentless, one-way movement of woody encroachment and agricultural conversion. This asymmetry supports prioritizing intact landscapes, as current restoration rates are often outpaced by the scale of habitat loss. Notably, we found that community-level spatial clustering is a more robust indicator of biome condition than trends of individual flagship species, suggesting that managing for ecosystem integrity provides a more effective multi-species umbrella. Given our results, there is an opportunity for operationalizing a Great Plains Conservation Design that is ecosystem-centric and rooted in the sustainability of the private-land cattle production that maintains these open spaces. By leveraging avian abundance as a biological sensor, managers and producers can deploy a shared vision that matches the spatial scale of the threats, moving from reactive triage to proactive defense of core working grasslands in North America.

Great Plains biome↗
USGS2025

Generalized Bancroft algorithm for locating earthquakes with P- and S-wave arrival times

Because of similarities between locating an earthquake with seismic stations and locating a Global Positioning System (GPS) receiver from satellites, the Bancroft algorithm developed for GPS processing can be used to locate earthquakes. Such an approach to earthquake location differs from the conventional method of choosing an initial or trial solution and then iteratively improving the solution until convergence. The Bancroft algorithm has the advantage of being a direct, noniterative solution but with the disadvantage of only being able to accommodate a homogeneous velocity model. An additional limitation of the standard Bancroft algorithm is that it considers arrival times in a medium with a single propagation velocity. This poses no problem for GPS processing because electromagnetic waves travel at the speed of light; however, for seismic waves it means the algorithm can be applied to collections of either P ‐ or S ‐wave arrival times. Here, I show how the Bancroft algorithm can be generalized to handle both P ‐ and S ‐wave arrival‐time measurements simultaneously. I also show how to accommodate depth‐varying P ‐ and S ‐wave velocity models. I apply the generalized Bancroft algorithm to microearthquakes beneath Tanaga Volcano in Alaska and compare standard locations from the widely used HYPOINVERSE location code to Bancroft locations and to the output of HYPOINVERSE when setting the trial location to the Bancroft location. I find the Bancroft locations outperform the results from the other methods for shallow earthquakes near sea level, where a quantity known as the geometric dilution of precision is large and linearized approaches such as HYPOINVERSE are expected to struggle.

Alaska↗
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↗
USGS2021

Factors influencing distribution of Coccidioides immitis in soil, Washington State, 2016

Coccidioides immitis and Coccidioides posadasii are causative agents of Valley fever, a serious fungal disease endemic to regions with hot, arid climate in the United States, Mexico, and Central and South America. The environmental niche of Coccidioide s spp. is not well defined, and it remains unknown whether these fungi are primarily associated with rodents or grow as saprotrophs in soil. To better understand the environmental reservoir of these pathogens, we used a systematic soil sampling approach, quantitative PCR (qPCR), culture, whole-genome sequencing, and soil chemical analysis to identify factors associated with the presence of C. immitis at a known colonization site in Washington State linked to a human case in 2010. We found that the same strain colonized an area of over 46,000 m 2 and persisted in soil for over 6 years. No association with rodent burrows was observed, as C. immitis DNA was as likely to be detected inside rodent holes as it was in the surrounding soil. In addition, the presence of C. immitis DNA in soil was correlated with elevated levels of boron, calcium, magnesium, sodium, and silicon in soil leachates. We also observed differences in the microbial communities between C. immitis -positive and -negative soils. Our artificial soil inoculation experiments demonstrated that C. immitis can use soil as a sole source of nutrients. Taken together, these results suggest that soil parameters need to be considered when modeling the distribution of this fungus in the environment.

Washington↗
USGS2024

SCEC/USGS Community Stress Drop Validation Study: How spectral fitting approaches influence measured source parameters

Spectral source parameters used to estimate an earthquake’s stress drop (Δσ) can vary significantly across measurement approaches. The Statewide California Earthquake Center/U.S. Geological Survey Community Stress-Drop Validation Study was initiated to compare source parameter estimates, focusing initially on a dataset from the 2019 Ridgecrest earthquake sequence. As part of that validation effort, here we focus on one potential source of uncertainty: whether spectral fitting approaches alone, applied to a common set of spectra from the 2019 Ridgecrest sequence result in different source parameter estimates. By using a common set of benchmark spectra analyzed across a consistent frequency band of 1–40 Hz, we eliminate many sources of variability. A subgroup of validation study participants volunteered to estimate the low-frequency displacement (Ω0) and corner frequency ( f c ) by fitting a smooth function to benchmark displacement spectra. Participants used linear- or log-sampled spectra, assumed a Brune or Boatwright spectral model, and applied different misfit criteria. We compare 17 approaches used to estimate Ω0, f c , and Δσ for 54 earthquake spectra. Our results reveal that 35% of events have Δσ estimates within a factor of two, whereas others exhibit variations exceeding an order of magnitude. The variability in Ω0 and f c can largely be attributed to whether a spectrum is consistent with the smooth function of an idealized simple crack model. The trade-off between Ω0 and f c may be more pronounced when using linearly sampled spectra, as higher frequency spectral bumps control the fits. As expected, methods that assumed a Boatwright model tended to have lower Ω0 and somewhat higher f c compared to those assuming a Brune model, although resulting Δσ estimates are similar. When compared to the overall validation study results, the fitting approach alone may account for between 5% and 90% (25% on average) of the total variability in spectral Δσ.

California↗
USGS2025

U.S. Geological Survey Integrated Water Availability Assessment—2010–20

This professional paper is a multichapter report that assesses water availability in the United States for water years 2010–20. This work was conducted as part of the fulfillment of the mandates of Subtitle F of the Omnibus Public Land Management Act of 2009 (Public Law 111-11), also known as the SECURE Water Act. As such, this work examines the spatial and temporal distribution of water quantity and quality in surface water and groundwater, as related to human and ecosystem needs and as affected by human and natural influences. Chapter A introduces the National Integrated Water Availability Assessment and provides important background and definitions for how the report characterizes water availability and its components. Chapter A also presents the key findings of Chapters B–F and thus acts as a summary of the entire report. Chapter B is a national assessment of water supply, which is the quantity of water supplied through climatic inputs. Chapter C is a national assessment of water quality, which is the chemical and physical characteristics of water. Chapter D assesses water use including withdrawals and consumptive use in the conterminous United States. Chapter E presents an analysis of factors affecting future water availability under changing climate conditions. The National Integrated Water Availability Assessment culminates with Chapter F, which is an integrated assessment of water availability that considers the amount and quality of water coupled with the suitability of that water for specific uses. Together, these six chapters constitute the National Integrated Water Availability Assessment for water years 2010–20.

Professional Paper↗
USGS2025

Potentiometric surface maps and groundwater-level hydrographs for confined aquifers of the New Jersey Coastal Plain, 2018

The U.S. Geological Survey, in cooperation with the New Jersey Department of Environmental Protection (NJDEP), prepared potentiometric surface maps for 10 confined aquifers of the New Jersey Coastal Plain physiographic province based on water-level measurements collected during late 2018 and early 2019 from 951 wells in New Jersey and parts of Pennsylvania and Delaware. Maps were prepared for the confined Cohansey aquifer, Rio Grande water-bearing zone, Atlantic City 800-foot sand, Piney Point aquifer, Vincentown aquifer, Wenonah-Mount Laurel aquifer, Englishtown aquifer system, and the upper, middle, and lower aquifers of the Potomac-Raritan-Magothy aquifer system. Potentiometric surface maps indicate regional cones of depression in the following aquifers and the counties in which they are centered: Atlantic City 800-foot sand in Atlantic County, the Piney Point aquifer in Cumberland County, the Wenonah-Mount Laurel aquifer and Englishtown aquifer system in Monmouth and Ocean Counties, the Wenonah-Mount Laurel aquifer in Camden and Gloucester County, the upper aquifer of the Potomac-Raritan-Magothy aquifer system in Ocean County, and the upper, middle, and lower aquifers of the Potomac-Raritan-Magothy aquifer system in Camden County. Cones of depression with smaller areal extents were in the confined Cohansey aquifer, the Rio Grande water-bearing zone, the Atlantic City 800-foot sand centered in Cape May County, the Piney Point aquifer centered in Ocean County, the Wenonah-Mount Laurel aquifer in Salem and Burlington Counties, the Englishtown aquifer system in Camden County, the upper aquifer of the Potomac-Raritan-Magothy aquifer system in Monmouth County, and the middle aquifer of the Potomac-Raritan-Magothy aquifer system in Monmouth, Ocean, and Salem Counties. No cone of depression was interpreted in the Vincentown aquifer. Long-term hydrographs are presented for 75 wells spanning each of the 10 confined aquifers, and contain a mix of discrete water-level measurements and daily mean water levels based on continuously recorded 15-minute data. Changes of water levels during 2014–19, as indicated by the hydrographs, were compared with those of previous periods to assess any departures from historical data. During 2014–19, water levels were stable and fluctuated within similar ranges as previous periods in the following aquifers and locations: all wells in the confined Cohansey aquifer, the Rio Grande water-bearing zone, the Vincentown aquifer, the Englishtown aquifer system, the Piney Point aquifer wells in Burlington and Ocean Counties, six of eight wells in the Wenonah-Mount Laurel aquifer, all wells in the upper and lower aquifers of the Potomac-Raritan-Magothy aquifer system outside NJDEP Critical Areas, and all wells in the middle aquifer of the Potomac-Raritan-Magothy aquifer system except those within NJDEP Critical Area II. Increasing water levels in 2014–19, ongoing since historical periods, were indicated in the following aquifers and locations: Atlantic County wells in the Piney Point aquifer, all wells in the upper and lower aquifers of the Potomac-Raritan-Magothy aquifer system outside NJDEP Critical Areas, and all wells in the middle aquifer of the Potomac-Raritan-Magothy aquifer system within NJDEP Critical Area II. Water levels in the Atlantic City 800-foot sand also increased during 2014–19 in wells in Atlantic County and northern Cape May County closer to the center of the cone of depression in that aquifer, which is a response unique to this period and absent from previous periods. During 2013–19, continued decreasing water levels, ongoing since previous periods, were indicated by hydrographs of Atlantic City 800-foot sand wells in southern Cape May County, Piney Point aquifer wells in Cumberland County where the regional cone of depression is located, and two wells in the Wenonah-Mount Laurel aquifer—070478, which in 2014–19 departed from previous periods, and 330020, which continued a gradual decrease throughout its period of record.

New Jersey↗
USGS2026

Decadal trends and occurrence of geogenic constituents and mixtures in groundwater across the continental United States

Worldwide, about 50% of the population is served by groundwater-sourced drinking water. Numerous groundwater quality assessments have found that geogenic constituents are among the most common contaminants in drinking-water aquifers. Documenting changing groundwater quality is a crucial aspect of water availability assessments. We assess trends and occurrence of geogenic constituent concentrations in groundwater across the continental United States using 3 decades of data from the U.S. Geological Survey’s National Water Quality Network. Thousands of groundwater wells were grouped into agricultural, urban, or domestic supply network types. Although most networks and constituents had no statistically significant change in concentration, many had increasing concentration trends, elevated concentrations, or both. Lithium, sodium, radium, sulfate, and uranium had increasing trends in more than 10% of the study networks. Urban and domestic well networks had increasing lithium and sodium trends more often than agricultural networks. Manganese most commonly increased in domestic well networks; uranium more commonly increased in agricultural and urban networks. Elevated concentration mixtures were widespread, and mixture complexities appeared to increase over time. Our results indicate that more than 2.3 million domestic-well users may be affected by elevated concentrations of one or more geogenic constituents.

continental United States↗
USGS2026

Bathymetric survey and storage capacity of Upper Lake Mary near Flagstaff, Arizona in 2024

The U.S. Geological Survey (USGS), in cooperation with the city of Flagstaff, collected bathymetric, light detection and ranging (lidar), and land-survey data of Upper Lake Mary in Arizona during the months of April and October 2024. The city of Flagstaff uses a combination of groundwater from well fields throughout the Flagstaff area and surface water, mainly from Upper Lake Mary, for its potable water supply. The purpose of the survey is to update previous surveys using new technology and compare the results to previous surveys to determine if there was a decrease in storage capacity that could affect the city’s water supply. The lakebed was mapped in April 2024 using a vessel equipped with a multibeam echosounder (MBES) and mobile lidar scanner with positioning captured using a real-time kinematic global navigation satellite system (RTK GNSS) base and receivers. In October 2024, areas of the reservoir that were too shallow for the boat and shoreline that were not captured by the vessel-based lidar were surveyed on foot using hand-held RTK GNSS receivers. At full pool (spillway elevation of 6,831.82 feet above NAVD 88 [2,082.34 meters (m)], Upper Lake Mary has a storage capacity of 16,449.80 acre-feet (20,290,611.73 cubic meters) and a surface area of 953.57 acres (3,860,926.075 square meters). The reservoir is 5.7 miles (9.7 kilometers) long and varies in width from 326 feet (99.36 m) near the central, narrow portion of the reservoir to 2,613 feet (796.44 m) in the upper portion. Comparisons between this survey and the previous two surveys from the 1950s and 2006 indicate no apparent decrease in reservoir area or storage capacity. Results of the 2024 survey indicate that Upper Lake Mary’s storage capacity increased by 0.9 percent from the 2006 survey and a 1.6 percent increase in surface area from the 2006 survey.

Arizona↗
USGS2026

Simulation of groundwater flow to evaluate hydrogeologic controls on a PFAS plume, Coakley Landfill Superfund site, Rockingham County, New Hampshire

Per- and polyfluoroalkyl substances (PFAS), including perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS), have been detected at combined concentrations above 2,000 nanograms per liter (ng/L) at groundwater seep locations near the Coakley Landfill Superfund site, in North Hampton, New Hampshire. The landfill was active from 1972 to 1985. An impermeable cap was placed on the landfill in 1998. The adjacent area to the Coakley Landfill has many water supply wells, and transport of PFAS compounds to the wells is a concern. Fracture anisotropy in the underlying bedrock aquifer complicates the understanding of PFAS transport because groundwater preferentially travels along fractures that may not align with the prevailing groundwater flow direction. In 2018, the U.S. Environmental Protection Agency and the U.S. Geological Survey began an investigation of the groundwater flow from the Coakley Landfill site. This report describes the modification of a numerical groundwater-flow model for the local area around the Coakley Landfill and summarizes findings of the investigation. In addition, this report includes a brief description of PFOA and PFOS occurrence, a discussion of model construction, evaluation of model performance through calibration, and discussion of simulation results for two periods (before and after capping). Limitations are also discussed. Results show that simulated groundwater flow moves from the Coakley Landfill to the west and north. Advective transport modeling using particle tracking shows that groundwater from the landfill discharges primarily to streams to the west and north, and a small amount is transported to distal wells. Dilution of contaminants through advection and dispersion likely plays a role in whether PFAS compounds from the landfill will be detected above laboratory reporting levels at distal wells.

New Hampshire↗
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↗
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↗