SEARCH · USGS Science

Results for “Earth Surface Processes and Landforms”

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,665 records · Page 25Linked to original sources

USGS2024

Ground deformation and gravity for volcano monitoring

Introduction When magma accumulates or migrates, it can cause pressurization and related ground deformation. Characterization of surface deformation provides important constraints on the potential for future volcanic activity, especially in combination with seismic activity, gas emissions, and other indicators. A wide variety of techniques and instrument types have been applied to the study of ground deformation at volcanoes (sidebar, p. 2; Dzurisin, 2000, 2003, 2007). Geodetic instruments include continuously recording Global Navigation Satellite System (GNSS; of which the United States’ Global Positioning System is one example) stations (fig. D1), borehole tiltmeters, and interferometric synthetic aperture radar (InSAR) measurements (from satellites, occupied and unoccupied aircraft systems, and ground-based sensors). Additional geodetic measurements like continuous- and survey-mode gravity (fig. D2) can contribute substantially to interpreting these data. Borehole strainmeters (see chapter K , this volume, by Hurwitz and Lowenstern, 2024) also have outstanding utility for monitoring deformation, although because of cost and permitting challenges, we do not include them as part of standard volcano monitoring networks for U.S. volcanoes. Still other techniques like light detection and ranging (lidar), structure from motion, and optical satellite data can be used to derive gross topographic changes, which can be used to map volcanic deposits, infer eruption rates, and gain insights into the source processes associated with eruptive activity (see chapter G , this volume, on tracking surface changes caused by volcanic activity; Orr and others, 2024). Experience has shown that no single geodetic monitoring technique is adequate to detect and track the entire range of ground-motion patterns that occur at volcanoes, primarily because of the temporal and spatial diversity of volcano deformation (fig. D3). Similarly, the magnitude of surface deformation varies widely. Geodetic monitoring strategies should therefore include multiple techniques and instrument types to cover a wide range of spatial and temporal scales. In identifying recommendations for geodetic instrumentation for volcano monitoring networks, we attempted to maximize the diversity of instrument types to measure the full range of deformation signals and minimize their expense and number; thus, we do not include several well-known deformation-monitoring techniques in our recommendations. Extensometers, for example, measure strains over distances of a few meters and have an excellent record of success in detecting changes in preeruptive localized ground motion across existing cracks, including at Mount St. Helens, Washington (Iwatsubo and others, 1992), and Piton de la Fournaise, Réunion Island (Peltier and others, 2006). Despite being relatively inexpensive, extensometers are best used primarily when localized ground displacements (for example, ground cracks) need to be tracked, and are not necessary at all volcanoes. In considering volcano deformation monitoring strategies, two complicating factors are deserving of special attention. First, not all deformation is driven by subsurface magmatic activity—for example, at many large stratovolcanoes (for example, Mount Rainier), flank collapses and landslides are significant geologic hazards (Reid and others, 2001) that may occur even in the absence of magmatic activity. Monitoring the stability of volcanoes is thus another critical application of geodetic monitoring networks to inform hazard assessment. One of the most famous examples of edifice instability is the large flank collapse that initiated the May 18, 1980, eruption of Mount St. Helens. Deformation monitoring had detected a bulge on the north flank of the mountain in April 1980 that was expanding by several meters per day (Lipman and others, 1981). Given that flank collapses can happen at any time during a period of volcanic unrest (or even outside a period of unrest), the capability to assess edifice stability is critical. Second, although volcanoes are commonly treated as idealized structures that erupt from single points, like centralvent stratovolcanoes, many are characterized by long rift zones from which eruptions may originate, and distributed volcanic fields are characterized by broadly spaced vents. For example, linear dikes are common at Kīlauea, Mauna Loa, and between Mount Shasta and Medicine Lake in California. At Kīlauea, one of these linear dikes emerged more than 40 kilometers (km) away from the summit of the volcano during the lower East Rift Zone eruption in 2018. Other volcanic fields, like Lassen volcanic center, California, or the San Francisco Volcanic Field, Arizona, have many small vents spread over a wide area. Although the instrumentation guidelines presented in this chapter remain phrased for central-vent volcanoes, they should be modified as needed in the context of the eruptive characteristics of each individual volcanic system. Spatial analysis of geodetic network coverage could help to ensure adequate instrumentation in areas where volcanism can occur over a broad area as opposed to a central vent. As an example, consider the adjacent volcanoes Mount Shasta and Medicine Lake. If station locations are chosen based only on the distance from the centers of the volcanoes, then any geodetic anomalies between the two volcanoes—an area of potential volcanism as indicated by the presence of volcanic features—may remain undetected by ground-based instrumentation. The spatial analysis is accomplished via a grid of pressure point sources (Mogi, 1958) evenly distributed across the map area, at a depth of 5 km in this example (fig. D4). Each source is inflated until predicted deformations exceed the GNSS white noise uncertainty estimates at one site (Langbein, 2017; Murray and Svarc, 2017). This volume of detectable magma provides a measure of the quality of the coverage (fig. D4). The results indicate that, as of 2022, there is a large area between Mount Shasta and Medicine Lake volcano with existing mapped dikes in which a substantial amount of magma could intrude without being detected geodetically. Applying this style of analysis to individual volcanic systems can provide a guide for designing network geometry given the expected locations of future eruptions.

Scientific Investigations Report↗

USGS2025

Preferential groundwater discharges along stream corridors are disregarded sources of greenhouse gases

Groundwater delivery of greenhouse gases (GHGs) to stream banks and riparian areas, before mixing with surface waters, has not been well quantified. We measured preferential groundwater delivery of GHGs to stream banks within three stream reaches, and found that stream banks with discharging groundwater emitted more CO 2 and were sources of N 2 O compared to stream banks without actively discharging groundwater, which emitted less CO 2 and were N 2 O sinks. At one of our stream reaches, groundwater CO 2 and N 2 O concentrations were 1.4–19.2 and 1.1–40.6 times higher than those in surface water, respectively, and groundwater delivery rates of CO 2 and N 2 O were 1.5 and 1.6 times higher than surface water emissions per unit area. On average, 21% (range 0%–100%) of CO 2 and N 2 O were emitted at the stream bank before mixing with surface waters. Preferential groundwater GHG emissions may contribute substantially to stream corridor emissions and may be underestimated when using a channel-centric approach to estimate riverine GHG budgets.

JGR Biogeosciences↗

USGS2025

Estimated average annualized tsunami losses for the United States

Tsunami hazards are substantial threats to coastal communities across the United States (U.S.) and its territories. U.S. states and territories collaborate through the National Tsunami Hazard Mitigation Program (NTHMP) to develop their own tsunami-hazard information for outreach and evacuation planning. An effort to curate this tsunami-hazard information to support comprehensive risk analysis at the national level has not yet been completed. In support of this effort, the Federal Emergency Management Agency (FEMA) collaborated with the NTHMP, the National Oceanic and Atmospheric Administration (NOAA) and the U.S. Geological Survey (USGS) starting in 2023. This collaboration included the collection and analysis of existing tsunami hazard data and methods in the U.S. Tsunami subject matter experts identified and selected scientifically defensible methods for estimating the risks to buildings and populations in coastal communities. These efforts may support decision making regarding resilience policies, priorities, strategies and funding levels. Tsunamis can be triggered by earthquakes, subaerial or submarine landslides, volcanic eruptions, glacial calving, near-earth objects, weather or other events. These events can cause severe destruction, injuries, and loss of life due to powerful currents and flooding. Tsunamis pose a substantial threat to the western United States and all U.S. territories, as described below. ■ Hawaii is threatened by distant tsunamis due to its central location in the Pacific Ocean basin and has a history of local events. ■ Alaska, particularly the Aleutian Islands, faces local tsunami threats due to proximity to the Alaska-Aleutian Subduction Zone, as well as distant tsunamis from around the Pacific Ocean basin. ■ The western coast of the U.S. is threatened by distant tsunamis from around the Pacific Ocean basin and local source tsunamis from earthquakes generated within the Cascadia Subduction Zone in the Pacific Northwest. ■ American Samoa faces local tsunami threats from earthquakes generated in the nearby Tonga Trench, as well as distant tsunami threats. ■ Guam and the Commonwealth of the Northern Mariana Islands are threatened by local tsunamis from the nearby Mariana Subduction Zone, as well as distant sources from around the Pacific Ocean Basin. ■ Puerto Rico and the United States Virgin Islands are threatened by multiple local and distant tsunami sources, such as the Puerto Rico Trench (PRT), given their location in the complex seismic region of the Caribbean Sea. Several historical events stand out because of their catastrophic impacts. ■ In the Pacific Northwest, the 1700 Cascadia earthquake caused a tsunami that affected coastal Native American communities, though the extent of the damage is not fully documented (Ludwin, et al., 2005). ■ In Puerto Rico, the 1918 earthquake triggered a tsunami that caused $77 million in damage in 2022 dollars and 116 fatalities, primarily along the western coast (Coffman et al., 1982). ■ The 1946 Aleutian Islands earthquake triggered a massive tsunami that devastated Hilo, Hawaii, killing 158 people and resulting in approximately $375 million in damage (adjusted to 2022 dollars) (Fisher et al., 2023). ■ The 1964 Alaska earthquake (M 9.2) generated tsunamis that caused severe destruction in some communities across Alaska, Oregon, and California. This disaster led to a total of 124 fatalities and approximately $2.9 billion in property damage (adjusted to 2022 dollars) (Brocher et al., 2014) (Alaska Science Center, 2024). ■ In American Samoa, a tsunami generated by the 2009 Samoa earthquake (Mw 8.1) caused widespread devastation, resulting in 34 confirmed fatalities (Apatu et al., 2013) and economic losses exceeding $160 million (adjusted to 2022 dollars) (DHS, 2011). More recent events, including the 2010 Chile earthquake, the 2011 Japan earthquake, and the 2022 Tonga volcanic eruption, resulted in millions of dollars in damage to numerous ports and harbors in the U.S. South Pacific territories, Hawaii, and along the west coast of the U.S. (Lynett, et al., 2022) (Wilson, et al., 2013). Since these events, the expansion of the built environment in lowlying areas along the coast has increased the exposure of buildings and people, thereby further escalating community risk from tsunamis. This report provides a comprehensive national assessment of earthquake-generated tsunami risk. It does not include impacts from tsunamis generated by landslides, volcanic eruptions, glacial calving, near-earth objects, weather, or other events. This study is based on the best available hazard data from the U.S. Pacific Coast (California, Oregon and Washington), Alaska, Hawaii, U.S. Pacific Territories (American Samoa, Guam and Commonwealth of the Northern Mariana Islands) and Caribbean Territories (Puerto Rico and United States Virgin Islands). Tsunami risks associated with states along the East Coast, Gulf Coast, and Great Lakes are not included in this study because Hazus 6.1 software (FEMA 2024a) does not currently include the ability to analyze tsunami risk in those states. Once modeling capabilities and tsunami hazard data become available for additional states, FEMA may incorporate these data into future editions of this study.

Alaska, California, Hawaii Oregon, Washington↗

USGS1999

Environmental setting and water-quality issues in the lower Tennessee River basin

The goals of the National Water-Quality Assessment Program are to describe current water-quality conditions for a large part of the Nation's water resources, identify water-quality changes over time, and identify the primary natural and human factors that affect water quality. The lower Tennessee River Basin is one of 59 river basins selected for study. The water-quality assessment of the lower Tennessee River Basin study unit began in 1997. The lower Tennessee River Basin study unit encompasses an area of about 19,500 square miles and extends from Chattanooga, Tennessee, to Paducah, Kentucky. The study unit had a population of about 1.5 million people in 1995. The study unit was subdivided into subunits with relatively homogeneous geology and physiography. Subdivision of the study unit creates a framework to assess the effects of natural and cultural settings on water quality. Nine subunits were delineated in the study unit; their boundaries generally coincide with level III and level IV ecoregion boundaries. The nine subunits are the Coastal Plain, Transition, Western Highland Rim, Outer Nashville Basin, Inner Nashville Basin, Eastern Highland Rim, Plateau Escarpment and Valleys, Cumberland Plateau, and Valley and Ridge.The lower Tennessee River Basin consists of predominantly forest (51 percent) and agricultural land (40 percent). Activities related to agricultural land use, therefore, are the primary cultural factors likely to have a widespread effect on surface- and ground-water quality in the study unit. Inputs of total nitrogen and phosphorus from agricultural activities in 1992 were about 161,000 and 37,900 tons, respectively. About 3.7 million pounds (active ingredient) of pesticides was applied to crops in the lower Tennessee River Basin in 1992. State water-quality agencies identified nutrient enrichment and pathogens as water-quality issues affecting both surface and ground water in the lower Tennessee River Basin. Water-quality data collected by State and Federal agencies between 1980 and 1996 were summarized to characterize surface- and ground-water quality of the subunits with respect to these issues. Median concentrations of nitrogen species generally were less than 1 milligram per liter in surface and ground water in all subunits, and were highest throughout the subunits that had the largest percentages of agricultural land use. Median phosphorus concentrations also were less than 1 milligram per liter in all subunits. Phosphatic limestones present in two subunits had a larger effect on phosphorus concentrations in surface and ground water than did the amount of agricultural land use in these subunits. Median counts of fecal coliform were higher in surface water than in ground water in all subunits. The highest median counts in surface water were in the Valley and Ridge (7,500 colonies per 100 milliliters) and the Outer Nashville Basin subunits (5,000 colonies per 100 milliliters). Highest median counts in ground water were in the Inner and Outer Nashville Basin subunit. Natural setting likely has an important effect with respect to fecal contamination of surface and ground water in the lower Tennessee River Basin.

Alabama, Georgia, Kentucky, Mississippi, Tennessee↗

USGS2024

Post-wildfire debris flows

Post-wildfire debris flows pose severe hazards to communities and infrastructure near and within recently burned mountainous terrain. Intense heat of wildfires changes the runoff characteristics of a watershed by combusting the vegetative canopy, litter, and duff, introducing ash into the soil and creating water repellant soils. Following wildfire, rainfall on bare ground is less able to infiltrate into the fire-altered soils and overland flow is less impeded by vegetation. Rainfall runoff in recently burned areas can erode hillslopes owing to the removal of soil binding organic matter near the soil surface by fire. In channels, loose, dry-ravel deposits composed of sand and gravel are readily entrained by concentrated runoff in channels. Entrainment of soil on hillslopes and in channels bulks up the sediment concentration of the rainfall runoff to generate debris flows capable of transporting boulders and large woody debris. Post-wildfire debris flows can be triggered by rainfall conditions that would typically produce little runoff during unburned conditions. The primary rainfall trigger for post-wildfire debris flows is high intensity rainfall during short duration convective rainstorms or periods of high rainfall intensity embedded within a long-duration frontal storm. Numerous observations of debris flows triggered by storms lasting less than an hour following periods of little to no rainfall indicate that antecedent rainfall is not a requirement for initiation of post-wildfire debris flows. Post-wildfire debris-flow hazard assessment entails estimating probability and magnitude of debris flows in the burned area, estimating debris-flow runout and intensity, and defining rainfall intensity-duration thresholds for debris-flow initiation. In the United States, probability and magnitude is estimated using empirically derived models largely based on data collected in southern California. The models provide maps to identify watersheds and drainage paths where post-wildfire hazards are most pronounced. Rainfall intensity-duration thresholds can be incorporated into flood hazard forecasting tools. Currently, work is underway to identify how to best implement debris-flow runout models in burned areas with efficiency and accuracy. Post-wildfire debris flows have been a long-recognized process in the Transverse Ranges of southern California; however, climate change is driving more frequent wildfires to burn more mountainous terrain throughout the western United States and worldwide. As a result, post-wildfire debris flows are becoming a more common threat in areas where they were once infrequent. As the threat of post-wildfire debris flow expands into new areas, evaluating the hazard becomes challenging because the degree to which wildfire increases debris-flow susceptibility varies from region to region. This chapter summarizes the knowledge to date for evaluating post-wildfire debris-flow susceptibility and hazard assessment. We summarize the characteristics of wildfire burn severity, topography, underlying soil and geology, and rainfall conditions that contribute to making a watershed most likely to produce post-wildfire debris flows. Methods for hazard assessment in the United States and other countries are summarized. We highlight knowledge gaps for how post-wildfire debris-flow susceptibility varies throughout the western United States and worldwide and identify research needs to improve hazard assessment methods in different geographies.

Book chapter↗

USGS2026

Insights into Mountain Pass carbonatite formation from in-situ sulfur isotopes and geochemistry of sulfate and sulfide minerals

The Mountain Pass carbonatite stock hosts a world-class rare earth element deposit and may be classified as a carbonate-sulfate igneous rock, as it contains on average > 50 volume percent carbonate minerals and 20 to 30 volume percent sulfate minerals. The sulfates range in composition from barite to celestine and locally occur with sparse sulfide minerals. We investigate the origin of sulfur enrichment and the occurrence of sulfur-bearing minerals in the Mountain Pass carbonatite with in-situ sulfur isotope and mineral chemistry. Barite cores with δ 34 S of 1 to 3‰ do not coexist with sulfides, whereas celestine rims with δ 34 S of > 3‰ are associated with sulfides with δ 34 S < -10‰. We propose a model in which sulfur-bearing sediments were subducted during episodes of plate convergence in the Mojave Province that preceded Mountain Pass magmatism. Metasomatism of the overlying mantle by melts derived from the subducted sediments generated an unusually carbon- and sulfur-rich source to yield carbonatite magmas. Sulfur from primary carbonatite magmas and ~ 1 to 7% sulfur from subducted sediment melts yielded a slightly enriched δ 34 S composition (relative to depleted mantle δ 34 S of -1‰) for early crystallizing barite. Celestine rims on magmatic barite cores formed at low, hydrothermal temperatures (< 350 °C) based on S isotope thermometry for equilibrium celestine-galena and celestine-pyrite pairs. The sparse sulfides in the carbonatite stock are not in equilibrium with the primary barite cores and therefore do not permit S isotope thermometry estimates of magmatic temperatures. The S/Se ratios of sulfide minerals (> 3,400) typically exceed primitive mantle values (S/Se of 3,340), also consistent with their derivation from hydrothermal fluids. Trace occurrences of sulfide and sulfate minerals in alkaline silicate stocks related to the carbonatite stock have similar δ 34 S compositions and yield similarly low formation temperatures, suggesting regionally extensive and chemically similar sulfur-bearing hydrothermal fluids that imparted lithologically diverse rocks with a consistent sulfur isotope fingerprint.

California, Nevada↗

USGS2025

Pre-Acadian tectonics of the eastern Orange-Milford Belt, south-central Connecticut

This excursion presents a reinterpretation of mapping and new analytical data from the eastern Orange-Milford belt (OMB) in south-central Connecticut. The OMB is a fault-bound terrane of argillites and mafic rocks of anomalously low metamorphic grade—and of poorly constrained ages and tectonic affinity—wedged between kyanite/sillimanite-grade peri- Laurentian rocks to the west and anatectic peri-Gondwanan rocks to the east (Fig. 1A). Our data demonstrate that Ordovician(?) igneous and sedimentary rocks of the OMB were variably metamorphosed in the Ordovician and Silurian but escaped regionally pervasive, high-grade Devonian and later metamorphism. Previous interpretations (Fritts 1963a, 1965a, 1965b; Burger, 1967; Burger and others, 1968; Rodgers, 1985) described these rocks as a conformable sequence of low-grade, Ordovician to Devonian metasediments and metavolcanics. Our results reveal that the “metavolcanics” are not extrusive rocks but rather slivers of lower oceanic crust with complicated high- and low-grade metamorphic fabrics, intruded by a swarm of Silurian sheeted basalt dikes, and in fault contact with the surrounding metasediments. The purpose of this trip is to show evidence of early Paleozoic (Taconic) deformation and metamorphism preserved in rocks of the eastern OMB. These rocks remained shallow, cool, and sufficiently dry during the regionally dominant Acadian and Alleghanian orogenies to have avoided significant overprinting. As such, these rocks serve as windows into a geologic history otherwise unavailable between anatectic rocks of the peri-Gondwanan Bronson Hill, Avalon, and Gander terranes east of the Hartford basin and sillimanite-grade rocks of the peri-Laurentian Hartland and gneiss dome belts west of the OMB. We present major and trace element geochemistry including rare-earth element patterns for all mafic units in the eastern OMB as well as 40 Ar/ 39 Ar age spectra of amphibole, muscovite, and K-feldspar from rocks of the Maltby Lakes complex (of Deasy and others, 2017), Savin Schist, and Wepawaug Schist. Our evidence demonstrates that the units of the OMB have been assembled by faulting or intrusion, and that no stratigraphic relationships exist between the argillaceous schists and the metaigneous rocks.

Connecticut↗

USGS2024

Field evidence and indicators of rockfall fragmentation and implications for mobility

Rockfall fragmentation can play an important role in hazard studies and the design of protective measures. However, the current lack of modeling tools that incorporate rock fragmentation mechanics is a limitation to enhancing studies and design. This research investigates the fragmentation patterns of rockfalls and analyzes the resulting distribution of fragment sizes within corresponding rockfall deposits. We focus on small rock fragments, which provide insights into the dynamics of the rockfall event and can be used as input for numerical modeling. We analyzed multiple rockfall events from locations worldwide, each exhibiting different degrees of fragmentation. Using image analysis techniques, we mapped all visible blocks, determined their volumes, and measured the distances they travelled from the initial point of impact. A key finding is the identification of three indicators of fragmentation. First, in cases where fragmentation was largely absent, we observed a trend of increasing block size with distance from the impact point or source area, which aligns with previously published findings. However, for energetic rockfall events characterized by intense fragmentation, we observed that small fragments exhibited longer travel distances compared to larger fragments. This distinction allowed us to differentiate blocks primarily resulting from the disaggregation process from those primarily resulting from dynamic fragmentation, with implications for rockfall mobility. Second, although the size distribution of rockfall deposits exhibits a power-law scaling for volumes larger than a minimum size threshold corresponding to a rollover of the distribution, in some case studies a deviation from power-law scaling is observed, indicating a process of larger block comminution due to fragmentation. Third, we found that rockfalls with fragmentation experience reduced mobility, indicated by higher reach angles, and higher lateral dispersion showing a wider distribution of trajectories. We interpret these findings as being directly related to the energy-consuming nature of fragmentation, which prevents farther deposition of fragmented rock blocks.

Albacete province, Lombardy and Aosta Valley, Yose↗

USGS2025

Geochemistry and radiogenic isotopes constrain the mantle source region of the Mountain Pass Intrusive Suite, California

The Mountain Pass carbonatite stock is the largest rare earth element (REE) deposit and only active REE mine in the United States. The carbonatite intrusion and spatially associated alkaline silicate intrusions constitute the Mountain Pass Intrusive Suite, which is located within the Mojave Province in California. Both the carbonatite and the alkaline silicate rocks are enriched in large ion lithophile elements and light REEs, and less enriched to depleted in high field strength elements, indicating the mantle source region was metasomatically enriched in incompatible trace elements. The cause of this metasomatic mantle enrichment and the genetic relationship between the carbonatite and the alkaline silicate stocks are poorly understood. In this study, major and trace element geochemical data and isotopic (Rb-Sr, Sm-Nd, and Lu-Hf) data are presented to constrain genesis of the Mountain Pass Intrusive Suite, from mantle source region to the intrusion of the stocks. Our geochemical data are consistent with derivation of the alkaline silicate and carbonatite melts through partial melting from a shared mantle source region rather than through liquid immiscibility or fractional crystallization and separation of a carbothermal fluid. Although the Rb-Sr isotopic system in the Mountain Pass Intrusive Suite is disturbed at the whole-rock scale, the isotopic systems for whole-rock Sm-Nd (εNd i = ‐2.2 ± 0.8) and zircon Lu-Hf (εHf i = 0.1 ± 1.1) are robust and support mantle derivation of the magmas. Geochemical modeling using experimentally derived partition coefficients was used to identify possible causes of enrichment in incompatible elements through metasomatism in the mantle source region. Modeling of metasomatism by melts derived by partial melting of deeply subducted carbonated sediments approximates observed Mountain Pass Intrusive Suite trace element chemistry. Scattered εHf i in inherited zircon (2.8 ± 2.6) is consistent with derivation from an arc-related environment with substantial crustal contamination. Paleotectonic studies in the Mojave Province indicate that regional subduction preceded emplacement of the Mountain Pass Intrusive Suite by ∼300 Ma. Melting of the Mountain Pass source region may have been caused by post-collisional thermal relaxation and extension.

California↗

USGS2024

Groundwater, surface-water, and water-chemistry data, Black Mesa area, northeastern Arizona—2019–2021

The Navajo (N) aquifer is an extensive aquifer and the primary source of groundwater in the 5,400-square-mile Black Mesa area in northeastern Arizona. Water availability is an important issue in the Black Mesa area because of the arid climate, past industrial water use, and continued water requirements for municipal use by a growing population. Precipitation in the area typically ranges from less than 6 to more than 16 inches per year, depending on location. The U.S. Geological Survey water-monitoring program in the Black Mesa area began in 1971 and provides information about the long-term effects of groundwater withdrawals from the N aquifer for industrial and municipal uses. This report presents the results of data collected as part of the monitoring program in the Black Mesa area from calendar years 2020–2021 and, additionally, uses streamflow statistics from November and December 2019. The monitoring program includes measurements of (1) groundwater withdrawals (pumping), (2) groundwater levels, (3) spring discharge, (4) surface-water discharge, and (5) groundwater chemistry. In calendar year 2020, total groundwater withdrawals were estimated to be 2,680 acre-feet (acre-ft), and, in 2021, total withdrawals were estimated to be 2,570 acre-ft. Total withdrawals during 2021 were about 65 percent less than total withdrawals in 2005 because the Peabody Western Coal Company discontinued its use of water to transport coal in a coal slurry pipeline after 2005 and ceased mining operations in 2019. Owing to Navajo Nation and Hopi Reservation access restrictions during the Coronavirus pandemic, water levels were not collected from municipal wells in 2020 or 2021. Water levels measured in 2021 from wells completed in the unconfined areas of the N aquifer within the Black Mesa area showed a decline in 7 of 13 wells when compared with water levels from the prestress period (prior to 1965). The changes in water levels across all 13 wells ranged from +8.4 feet (ft) to −42.4 ft, and the median change was −0.4 ft. Water levels also showed decline in 11 of 12 wells measured in the confined area of the aquifer when compared to the prestress period. The median change for the confined area of the aquifer was −25.9 ft, with changes across all 12 wells ranging from +17.3 ft to −133.7 ft. Spring flow was measured at four springs between 2020 and 2021. Flow fluctuated during the period of record for Burro Spring and Pasture Canyon Spring, but a decreasing trend was statistically significant (p<0.05) at Moenkopi School Spring and Unnamed Spring near Dennehotso, Arizona. Discharge at Burro Spring has remained relatively constant since it was first measured in the 1980s, and discharge at Pasture Canyon Spring has fluctuated for the period of record. Continuous records of surface-water discharge in the Black Mesa area were collected from streamflow-gaging stations at the following sites: Moenkopi Wash at Moenkopi 09401260 (1976–2021), Dinnebito Wash near Sand Springs 09401110 (1993–2020), Polacca Wash near Second Mesa 09400568 (1994–2020), and Pasture Canyon Springs 09401265 (2004–2021). Median winter flows (November through February) of each winter were used as an estimate of the amount of groundwater discharge at the above-named sites. For the period of record, the median winter flows have generally remained constant at Polacca Wash and Pasture Canyon Springs, whereas a decreasing trend was observed at Moenkopi Wash and Dinnebito Wash. In 2020 and 2021, water samples were collected from a total of four springs in the Black Mesa area and analyzed for selected chemical constituents. Results from the four springs were compared with previous analyses from the same springs. Dissolved solids, chloride, and sulfate concentrations increased at Moenkopi School Spring during the more than 30 years of record at that site. Concentrations of dissolved solids and sulfate at Pasture Canyon Spring have not varied significantly (p>0.05) since the early 1980s, and there is no increasing or decreasing trend in those data. However, concentrations of chloride from Pasture Canyon Spring show a diminishing trend. Concentrations of dissolved solids, chloride, and sulfate at Unnamed Spring near Dennehotso have varied for the period of record, but there is no statistical trend in the data. Concentrations of dissolved solids at Burro Spring have varied for the period of record, but there is no statistical trend in the data. However, concentrations of chloride and sulfate from Burro Spring show a trend towards lower concentrations.

Arizona↗

USGS2026

Simulation of the impacts of spring diversions on streamflow in the Strawberry Creek watershed, San Bernardino County, California, using an integrated hydrological model

The Strawberry Creek watershed, in the San Bernardino Mountains of southern California, features a group of natural springs known as Arrowhead Springs that have been augmented with diversions in the form of sub-horizontal borings and tunnels. Understanding the effects of these structures on streamflow through groundwater capture is crucial for managing surface-water resources in this watershed. In this study, we constructed the Strawberry Creek integrated hydrological model to increase this understanding. The Strawberry Creek integrated hydrological model is an integrated surface runoff and groundwater model, GSFLOW to simulate surface runoff and infiltration and groundwater flow. The model has 263 rows and 176 columns in each of its three layers. The model area includes the Strawberry Creek and four adjacent watersheds. The precipitation-runoff modeling system was calibrated using two streamflow gaging stations, and the groundwater and surface-water flow model was calibrated to the reported spring diversion discharge and a sparse number of groundwater-level measurements. The Strawberry Creek integrated hydrological model was run with and without diversions active, and simulated streamflow was compared, indicating that in the headwaters of Strawberry Creek, about 35 percent of the diversion flow was captured from base flow.

California↗

USGS2025

Simulation of the impacts of spring diversions on streamflow in the Strawberry Creek watershed, San Bernardino County, California, using an integrated hydrological model

The Strawberry Creek watershed, situated in the San Bernardino Mountains of southern California, features a group of natural springs known as Arrowhead Springs that have been augmented with diversions in the form of sub-horizontal borings and tunnels. Understanding the impact of these structures on streamflow through groundwater capture is crucial for managing surface-water resources in this watershed. In this study we constructed the Strawberry Creek integrated hydrological model (SCIHM) to increase this understanding. The SCIHM is an integrated surface runoff and groundwater model that uses the coupled groundwater and surface-water flow model (GSFLOW), which is based on the integration of the precipitation-runoff modeling system (PRMS) and the modular groundwater flow model commonly called MODFLOW, version MODFLOW-2005 software to simulate surface runoff and infiltration and groundwater flow. The model has three layers, 263 rows, and 176 columns. The model area includes the Strawberry Creek and four adjacent watersheds. The PRMS model was calibrated using two streamflow gaging stations and the GSFLOW model was calibrated to reported spring diversion discharge and a sparse number of groundwater-level measurements. The SCIHM was run with and without diversions active and simulated streamflow was compared, finding that in the headwaters of Strawberry Creek about 35 percent of the diversion flow was captured from base flow.

EarthArXiv↗

USGS2026

Ten new insights in climate science 2025

Interdisciplinary understanding is vital for delivering sound climate policy advice. However, navigating the ever-growing and increasingly diverse scholarly literature on climate change is challenging for any individual researcher. This annual synthesis highlights and explains recent advances across a variety of fields of climate change research. This year, the 10 insights focus on: (1) the record-warmth of 2023/2024 and the elevated Earth energy imbalance; (2) acceleration of ocean warming and intensifying marine heatwaves; (3) northern land carbon sinks under strain; (4) reinforcing feedback between biodiversity loss and climate change; (5) accelerated depletion of groundwater; (6) global dengue incidence; (7) global income losses and labour productivity declines; (8) strategic scaling of CDR; (9) integrity challenges in carbon credit markets and emerging responses; and (10) effective policy mixes for emissions reductions. The insights have been written to be accessible to researchers from different fields, serving as entry-points to specific topics, as well as providing an overview of the evolving landscape of climate change research. In the final section, the insights are used to develop overarching policy-relevant messages. This paper provides the basis for a science-policy report that was shared with all Party delegations ahead of COP30 in Belém, Brazil.

Global Sustainability↗

USGS1977

Acid-rain weathering of a metasedimentary rock basin, Herkimer County, New York

A geochemical balance of major chemical constituents of precipitation inflow with those of surface-water outflow was used to evaluate the effect of acid rainfall (average pH = 4.3) on weathering rates of metasedimentary rocks underlying a surface-water basin in New York State. Precipitation samples and surface-water samples show that approximately 45 percent of the net cation outflow is due to weathering by precipitation acidity. The remaining 55 percent of the cation load in the outflow is attributed to interactions between atmosphere, rock, and soil water in the unsaturated zone. Comparison of the acidity of acid precipitation and air-equilibrated distilled water suggests that since preindustrial time, precipitation acidity has probably increased weathering rates for this basin approximately 75 percent. The increased acidity of the precipitation may also increase the rate of dissolution of heavy metals from water-piping systems.

New York↗

USGS2026

Aspergillosis (Avian) case definition for wildlife

Diagnostic laboratories receive carcasses and samples for diagnostic evaluation and pathogen/toxin detection. Case definitions bring clarity and consistency to the evaluation process. Their use within and between organizations allows more uniform reporting of diseases and etiologic agents. The intent of a case definition is to provide scientifically based criteria for determining: (a) if an individual carcass has a specific disease and degree of confidence in that diagnosis and (b) if there is evidence of a pathogen or toxin in a carcass or sample (for example, swab, tissue sample, skin scraping, blood/serum sample, environmental sample, or other). This case definition is specific to aspergillosis and applies to all avian species.

Techniques and Methods↗

USGS2026

Water-withdrawal and groundwater-level fluctuations, physical and infiltration characteristics of Ohio River sediment, and isotope-based estimates of water sources relative to a riverbank filtration system in an alluvial aquifer, Louisville, Kentucky, 2000–20

The U.S. Geological Survey and the Louisville Water Company, from 2000 to 2020, cooperatively investigated physical and hydrologic conditions in the Ohio River alluvial aquifer and adjacent riverbed sediments in the northeast part of Louisville, Kentucky, to increase knowledge of processes affecting water sources to collector wells during operation of a riverbank filtration system (RBF). The RBF, in 2020, included an initial collector well (CW1) completed in 1999 and a tunnel and collector well system (TCWS) composed of four collector wells (CW2–CW5) completed in late 2010 that were connected by a tunnel that routed water to a pumping station and a treatment plant. Daily mean water-withdrawal rates typically ranged from about 10.4 to 20.7 million gallons per day (Mgal/d) from 2005 to 2007, then were cyclic from 2008 to 2010 and ranged from as low as about 9 to 13 Mgal/d in winter-spring months up to about 19 to 20 Mgal/d in summer-autumn months. From 2011 through 2017, peak daily mean water-withdrawal rates from CW2–CW5 fluctuated between winter-spring lows of about 20–30 Mgal/d to summer-autumn peak rates of about 41.6 Mgal/d in 2013 to about 48.2 Mgal/d in 2012. From 2017 to 2020, TCWS peak daily mean water-withdrawal rates decreased from about 45.7 Mgal/d in 2017 to about 37.5 Mgal/d in 2020 and peak daily mean rates from CW1 ranged from about 20.7 Mgal/d in 2020 to about 21.6 Mgal/d in 2018. Periods in 2019–20 with tunnel pumping station daily mean water-withdrawal rates smaller than about 30 Mgal/d had generally stable water levels in CW2 and indicated the sustainable yield of the TCWS collector wells under the prevailing hydrologic conditions. Drawdown during July–August of 2008 and 2009 was associated with CW1, the sole operating collector well. Groundwater-flow directions outside the CW1 cone of depression during July–August of 2008 and 2009 were from areas near the alluvial valley margin toward the Ohio River and in 2009 toward Harrods Creek. Collector wells CW2 through CW5 derived water during 2011–19 surveys from surface-water infiltration, as indicated by drawdown that extended beneath the Ohio River and Harrods Creek. Increases in drawdown toward the alluvial valley margins in July–August of 2013 and 2015, compared with July 2011, indicated less available recharge from those flow directions and decreased yield of groundwater from the aquifer after progressively longer TCWS water withdrawals. Drawdown at CW2–CW5 and at observation wells between them lessened in October 2019 and indicated that onshore groundwater levels can recover when collector well withdrawals decrease. Overlapping cones of depression in 2011–15 potentiometric surfaces centered around CW2, CW3, CW4, and CW5, and drawdown along divides between those wells indicated well interference between collector wells. Well-interference effects increased during 2013 and 2015 survey periods, despite decreased water-withdrawal rates compared with July 2011, and expanded to include CW1 when its use resumed. Groundwater-level fluctuations in observation wells near the Ohio River were consistent with gradual depletion of groundwater yield during tunnel pumping station operation at river elevations of about 420 feet above the National Geodetic Vertical Datum of 1929 (NGVD 29) and enhanced recharge of the aquifer when river elevations were above flood stage of about 430 feet above NGVD 29. Groundwater-level declines in observation wells near the Ohio River were most apparent in summer and some autumn months between 2012 and 2017 during periods when tunnel pumping station mean daily water-withdrawal rates were about 35 Mgal/d or more. Periods of relatively stable groundwater levels during 2019–20 in observation wells between CW2–CW5, when tunnel pumping station water-withdrawal rates were about 30–33 Mgal/d or less, also indicated that the prevailing sustainable yield was about 30 Mgal/d or less. Wells associated with CW3 and CW4 had the largest increases in slopes of best-fit lines relating water-level difference to water-withdrawal rate between 2012 and 2020, an indication of decreased efficiencies of those wells to supply water or of changes in groundwater yield around the well. Increased slopes of best-fit lines relating water-level difference to water-withdrawal rate from 2019 to 2020, when groundwater levels were otherwise relatively stable, also indicated that the specific capacity of one or more TCWS collector wells may have declined during that period. Riverbed surface elevation, roughness, and hardness changes between hydroacoustic surveys indicated active accretion and erosion of bed sediment offshore from the riverbank filtration system in 2002–4. Riverbed hydraulic conductivity computed from 2013 to 2015 in tests at 11 sites ranged from 0.2 to 230 feet per day. Ten temporary piezometer sites with adjacent operating collector wells had drawdowns that ranged from about 4 to 11.7 feet and hammer blows greater than or equal to 49 counts, which indicated a relatively great density of alluvial sediments and the potential for diminished riverbed infiltration rates. Contributions of Ohio River infiltration to collector wells estimated from deuterium and oxygen-18 isotopic ratios in representative Ohio River and groundwater samples were considerably larger in warmer months than cooler months and ranged from 80–90 percent in CW1–CW5 in May 2020 to 100 percent in May 2019 at CW2, CW4, and CW5. Groundwater contributions to water withdrawals by collector wells increased in cooler autumn months and ranged from 40 percent at CW3 and CW5 in November 2020 to 70 percent at CW2, CW4, and CW5 in November–December 2019. November 2019 and November 2020 isotopic compositions of water from the collector wells had enriched deuterium and oxygen-18 isotopic ratios consistent with mixtures of evaporation-modified water from the Ohio River and meteoric-origin groundwater.

Indiana, Kentucky↗

USGS2000

Quality assurance and analysis of water levels in wells on Pahute Mesa and vicinity, Nevada Test Site, Nye County, Nevada

Periodic and continual water-level data from 1963 to 1998 were compiled and quality assured for 65 observation wells on Pahute Mesa and vicinity, Nye County, Nevada. As part of the quality assurance of all water levels, ancillary data pertinent to computing hydraulic heads in wells were compiled and analyzed. Quality-assured water levels that were not necessarily in error but which did not represent static heads in the regional aquifer system, or required some other qualification, were flagged. Water levels flagged include those recovering from recent pumping or well construction, water levels affected by nuclear tests, and measurements affected by borehole deviations. A cursory examination of about 30 wells with available water-level and down-hole temperature data indicate that water levels in most wells on Pahute Mesa would not be significantly affected by temperature if corrected to 95 degrees Fahrenheit. Wells with large corrections (greater than 10 feet) are those with long water columns (greater than 1,500 feet of water above the assumed point of inflow) in combination with mean water-column temperatures exceeding 105 degrees Fahrenheit. Water-level fluctuations in wells on Pahute Mesa are caused by several factors including infiltration of precipitation, barometric pressure, Earth tides, ground-water pumpage, and seismic events caused by tectonic activity and underground nuclear testing. No observed water-level fluctuations were attributed to a naturally occurring earthquake. The magnitude and duration of changes in water levels caused by nuclear tests are affected by the test size and the distance from a well to the test. Identifying water levels that might be affected by past nuclear tests is difficult because pre-testing water-level data are sparse. Hydrologically significant trends were found in 13 of 25 wells with multiple years of water-level record. The largest change in water levels (1,029 feet in 25 years) occurred in well U-19v PS 1D as a result of the Almendro nuclear test. Likely explanations for trends in most of the wells are either changes in precipitation patterns that affect recharge rates to the ground-water system, pumping effects from water-supply well U-20 WW, or a combination of these two factors.

Nevada↗

USGS2024

Using the horizontal-to-vertical spectral ratio method to estimate thickness of the Barry Arm landslide, Prince William Sound, Alaska

Conducting detailed investigations of large landslides is difficult, especially in the subsurface, largely due to environmental factors such as steep slopes, difficult access, and numerous objective hazards. These factors have made it challenging to accurately estimate the depth to the failure surface of the Barry Arm landslide, a large (roughly 10 8 cubic meters), deep-seated bedrock landslide in Prince William Sound, Alaska, recognized in 2019. The landslide has exhibited accelerated movement in recent years and poses a potential tsunamigenic hazard if rapid failure occurs. Failure surface depth, equivalent to landslide thickness, is a necessary metric for landslide-volume calculations and associated tsunami wave models. In this report, we used seismic noise recorded by a seismometer located on the Barry Arm landslide in Alaska to calculate the horizontal-to-vertical spectral ratio (HVSR) to investigate the site fundamental frequency ( f 0 ) and depth of the failure surface. To ensure that observed peak frequencies in the spectral ratio were related to the underlying stratigraphy (and not caused by other noise sources like nearby glaciers, topographic resonance, weather, or human activities), we also calculated HVSRs using earthquake signals, HVSRs at other seismic stations within a 2.5-kilometer radius, and a standard spectral ratio between the landslide station and other sites. We observed multiple peaks in the landslide HVSR curves at 1.5 hertz (Hz), 4–5 Hz, and 7–11 Hz. The frequencies of these peaks were consistent at the landslide site through time and across methods and were dissimilar to those identified at other seismic stations in the area, making it unlikely the peaks were caused by local noise. Directional HVSRs calculated at 15-degree intervals showed amplification of the higher frequency peaks in the direction parallel to slip, indicating two-dimensional site effects. We used the distinct frequency peaks in the seismic record to develop a 4-layer conceptual model of the landslide wherein the top of the deepest layer represents the primary failure surface, or the boundary between damaged (mobile) and undamaged material. We inverted Rayleigh wave ellipticity curves within this 4-layer configuration with constraints on S-wave velocity and layer thickness based on analogous material properties identified in the literature. This was necessary absent any site-specific subsurface S-wave velocity data. The best-fitting models indicate a mean slope-normal depth to the failure surface of 188 (±9) meters (m), with additional stratigraphic boundaries at 4 and 20 m below ground surface, potentially representing layered motion. These results agree with and improve upon ranges estimated by previous studies and can support future modeling and assessment efforts at Barry Arm.

Alaska↗