Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Soil Science”

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.

951 records · Page 53Linked to original sources

Unraveling the dynamics of magmatic CO2 degassing at Mammoth Mountain, California

The accumulation of magmatic CO 2 beneath low-permeability barriers may lead to the formation of CO 2 -rich gas reservoirs within volcanic systems. Such accumulation is often evidenced by high surface CO 2 emissions that fluctuate over time. The temporal variability in surface degassing is believed in part to reflect a complex interplay between deep magmatic degassing and the permeability of degassing pathways. A better understanding of the dynamics of CO 2 degassing is required to improve monitoring and hazards mitigation in these systems. Owing to the availability of long-term records of CO 2 emissions rates and seismicity, Mammoth Mountain in California constitutes an ideal site towards such predictive understanding. Mammoth Mountain is characterized by intense soil CO 2 degassing (up to ∼1000 t d −1 ) and tree kill areas that resulted from leakage of CO 2 from a CO 2 -rich gas reservoir located in the upper ∼4 km. The release of CO 2 -rich fluids from deeper basaltic intrusions towards the reservoir induces seismicity and potentially reactivates faults connecting the reservoir to the surface. While this conceptual model is well-accepted, there is still a debate whether temporally variable surface CO 2 fluxes directly reflect degassing of intrusions or variations in fault permeability. Here, we report the first large-scale numerical model of fluid and heat transport for Mammoth Mountain. We discuss processes (i) leading to the initial formation of the CO 2 -rich gas reservoir prior to the occurrence of high surface CO 2 degassing rates and (ii) controlling current CO 2 degassing at the surface. Although the modeling settings are site-specific, the key mechanisms discussed in this study are likely at play at other volcanic systems hosting CO 2 -rich gas reservoirs. In particular, our model results illustrate the role of convection in stripping a CO 2 -rich gas phase from a rising hydrothermal fluid and leading to an accumulation of a large mass of CO 2 (∼10 7 –10 8 t) in a shallow gas reservoir. Moreover, we show that both, short-lived (months to years) and long-lived (hundreds of years) events of magmatic fluid injection can lead to critical pressures within the reservoir and potentially trigger fault reactivation. Our sensitivity analysis suggests that observed temporal fluctuations in surface degassing are only indirectly controlled by variations in magmatic degassing and are mainly the result of temporally variable fault permeability. Finally, we suggest that long-term CO 2 emission monitoring, seismic tomography and coupled thermal–hydraulic–mechanical modeling are important for CO 2 -related hazard mitigation.

California↗

Subsidence reversal in a re-establish wetland in the Sacramento-San Joaquin Delta, California, USA

The stability of levees in the Sacramento-San Joaquin Delta is threatened by continued subsidence of Delta peat islands. Up to 6 meters of land-surface elevation has been lost in the 150 years since Delta marshes were leveed and drained, primarily from oxidation of peat soils. Flooding subsided peat islands halts peat oxidation by creating anoxic soils, but net accumulation of new material in restored wetlands is required to recover land-surface elevations. We investigated the subsidence reversal potential of two 3 hectare, permanently flooded, impounded wetlands re-established on a deeply subsided field on Twitchell Island. The shallower wetland (design water depth 25 cm) was almost completely colonized by dense emergent marsh vegetation within two years; whereas, the deeper wetland (design water depth 55 cm) which developed spatially variable depths as a result of heterogeneous colonization by emergent vegetation, still had some areas remaining as open water after nine years. Changes in land-surface elevation were quantified using repeated sedimentation-erosion table measurements. New material accumulating in the wetlands was sampled by coring. Land-surface elevations increased by an average of 4 cm/yr in both wetlands from 1997 to 2006; however, the rates at different sites in the wetlands ranged from -0.5 to +9.2 cm/yr. Open water areas of the deeper wetland without emergent vegetation had the lowest rates of land-surface elevation gain. The greatest rates occurred in areas of the deeper wetland most isolated from the river water inlets, with dense stands of emergent marsh vegetation (tules and cattails). Vegetated areas of the deeper wetland in the transition zones between open water and mature emergent stands had intermediate rates of land-surface gain, as did the entire shallower wetland. These results suggest that the dominant component contributing to land-surface elevation gain in these wetlands was accumulation of organic matter, rather than mineral sediment, and that accumulation of organic matter in emergent marshes is strongly affected by hydrologic factors. Re-established, non-tidal wetlands with managed hydrology can produce significant increases in land-surface elevations, which can help to improve levee stability and protect subsided islands from future flooding.

California↗

Emergence of a zoonotic pathogen in a coastal marine sentinel: Capillaria hepatica (syn. Calodium hepaticum)-associated hepatitis in southern sea otters (Enhydra lutris nereis)

Capillaria hepatica is a globally distributed zoonotic nematode parasite that most commonly infects feral and native rats. Soil contact, pica, and living in close proximity to rat populations are risk factors for human infection. Larval nematodes and eggs that were morphologically consistent with C. hepatica were observed microscopically in livers of stranded southern sea otters ( Enhydra lutris nereis ) from California. Large (90–100 × 45–55 μm), barrel-shaped non-embryonated parasite eggs with large polar prominences and a roughened or striated surface, or 105–120 μm diameter larval aphasmid nematode profiles with a prominent stichosome and hypodermal bands were observed in the livers of three otters. The liver of a fourth animal exhibited serpiginous tracts of necrosis, micro-cavitation and pleocellular inflammation, with intralesional linear eosinophilic material that resembled cuticle from degenerating metazoan parasites. Capillaria hepatica -associated hepatitis and capsular adhesions were the cause of death for one otter, and parasite-associated liver lesions were observed in all cases. All infected otters were adult females that stranded from 2006 through 2016 at multiple sites along the central California coast. All cases stranded from December through May; during and soon after peak seasonal precipitation and land-sea runoff in California. This same seasonal pattern has been reported for other land-based parasites infecting southern sea otters. Neither C. hepatica , nor any similar nematodes have been reported from marine mammals, and southern sea otters are not typical hosts for C. hepatica or any other nematode parasites. The most likely route of exposure was via freshwater runoff containing embryonated eggs liberated from predated or decomposing terrestrial hosts, especially rats. Similar to the land-based parasites Toxoplasma gondii and Sarcocystis neurona , C. hepatica eggs may be concentrated and transmitted through filter-feeding marine invertebrates that serve as southern sea otter prey, which may also pose an unrecognized public health risk for people who consume these species.

California↗

Geochemistry of coastal permafrost and erosion-driven organic matter fluxes to the Beaufort Sea near Drew Point, Alaska

Accelerating erosion of the Alaska Beaufort Sea coast is increasing inputs of organic matter from land to the Arctic Ocean, and improved estimates of organic matter stocks in eroding coastal permafrost are needed to assess their mobilization rates under contemporary conditions. We collected three permafrost cores (4.5–7.5 m long) along a geomorphic gradient near Drew Point, Alaska, where recent erosion rates average 17.2 m year −1 . Down-core patterns indicate that organic-rich soils and lacustrine sediments (12–45% total organic carbon; TOC) in the active layer and upper permafrost accumulated during the Holocene. Deeper permafrost (below 3 m elevation) mainly consists of Late Pleistocene marine sediments with lower organic matter content (∼1% TOC), lower C:N ratios, and higher δ 13 C values. Radiocarbon-based estimates of organic carbon accumulation rates were 11.3 ± 3.6 g TOC m −2 year −1 during the Holocene and 0.5 ± 0.1 g TOC m −2 year −1 during the Late Pleistocene (12–38 kyr BP). Within relict marine sediments, porewater salinities increased with depth. Elevated salinity near sea level (∼20–37 in thawed samples) inhibited freezing despite year-round temperatures below 0°C. We used organic matter stock estimates from the cores in combination with remote sensing time-series data to estimate carbon fluxes for a 9 km stretch of coastline near Drew Point. Erosional fluxes of TOC averaged 1,369 kg C m −1 year −1 during the 21st century (2002–2018), nearly doubling the average flux of the previous half-century (1955–2002). Our estimate of the 21st century erosional TOC flux year −1 from this 9 km coastline (12,318 metric tons C year −1 ) is similar to the annual TOC flux from the Kuparuk River, which drains a 8,107 km 2 area east of Drew Point and ranks as the third largest river on the North Slope of Alaska. Total nitrogen fluxes via coastal erosion at Drew Point were also quantified, and were similar to those from the Kuparuk River. This study emphasizes that coastal erosion represents a significant pathway for carbon and nitrogen trapped in permafrost to enter modern biogeochemical cycles, where it may fuel food webs and greenhouse gas emissions in the marine environment.

Alaska↗

Near real-time detection of winter cover crop termination using harmonized Landsat and Sentinel-2 (HLS) to support ecosystem assessment

Cover crops are planted to reduce soil erosion, increase soil fertility, and improve watershed management. In the Delmarva Peninsula of the eastern United States, winter cover crops are essential for reducing nutrient and sediment losses from farmland. Cost-share programs have been created to incentivize cover crops to achieve conservation objectives. This program required that cover crops be planted and terminated within a specified time window. Usually, farmers report cover crop termination dates for each enrolled field (∼28,000 per year), and conservation district staff confirm the report with field visits within two weeks of termination. This verification process is labor-intensive and time-consuming and became restricted in 2020–2021 due to the COVID-19 pandemic. This study used Harmonized Landsat and Sentinel-2 (HLS, version 2.0) time-series data and the within-season termination (WIST) algorithm to detect cover crop termination dates over Maryland and the Delmarva Peninsula. The estimated remote sensing termination dates were compared to roadside surveys and to farmer-reported termination dates from the Maryland Department of Agriculture database for the 2020–2021 cover crop season. The results show that the WIST algorithm using HLS detected 94% of terminations (statuses) for the enrolled fields (n = 28,190). Among the detected terminations, about 49%, 72%, 84%, and 90% of remote sensing detected termination dates were within one, two, three, and four weeks of agreement to farmer-reported dates, respectively. A real-time simulation showed that the termination dates could be detected one week after termination operation using routinely available HLS data, and termination dates detected after mid-May are more reliable than those from early spring when the Normalized Difference Vegetation Index (NDVI) was low. We conclude that HLS imagery and the WIST algorithm provide a fast and consistent approach for generating near-real-time cover crop termination maps over large areas, which can be used to support cost-share program verification.

Maryland↗

Differential rates of feldspar weathering in granitic regoliths

Differential rates of plagioclase and K-feldspar weathering commonly observed in bedrock and soil environments are examined in terms of chemical kinetic and solubility controls and hydrologic permeability. For the Panola regolith, in the Georgia Piedmont Province of southeastern United States, petrographic observations, coupled with elemental balances and 87Sr/86Sr ratios, indicate that plagioclase is being converted to kaolinite at depths > 6 m in the granitic bedrock. K-feldspar remains pristine in the bedrock but subsequently weathers to kaolinite at the overlying saprolite. In contrast, both plagioclase and K-feldspar remain stable in granitic bedrocks elsewhere in Piedmont Province, such as Davis Run, Virginia, where feldspars weather concurrently in an overlying thick saprolite sequence. Kinetic rate constants, mineral surface areas, and secondary hydraulic conductivities are fitted to feldspar losses with depth in the Panola and Davis Run regoliths using a time-depth computer spreadsheet model. The primary hydraulic conductivities, describing the rates of meteoric water penetration into the pristine granites, are assumed to be equal to the propagation rates of weathering fronts, which, based on cosmogenic isotope dating, are 7 m/106 yr for the Panola regolith and 4 m/106 yr for the Davis Run regolith. Best fits in the calculations indicate that the kinetic rate constants for plagioclase in both regoliths are factors of two to three times faster than K-feldspar, which is in agreement with experimental findings. However, the range for plagioclase and K-feldspar rates (kr = 1.5 x 10-17 to 2.8 x 10-16 mol m-2 s-1) is three to four orders of magnitude lower than for that for experimental feldspar dissolution rates and are among the slowest yet recorded for natural feldspar weathering. Such slow rates are attributed to the relatively old geomorphic ages of the Panola and Davis Run regoliths, implying that mineral surface reactivity decreases significantly with time. Differential feldspar weathering in the low-permeability Panola bedrock environment is more dependent on relative feldspar solubilities than on differences in kinetic reaction rates. Such weathering is very sensitive to primary and secondary hydraulic conductivities (qp and qs), which control both the fluid volumes passing through the regolith and the thermodynamic saturation of the feldspars. Bedrock permeability is primarily intragranular and is created by internal weathering of networks of interconnected plagioclase phenocrysts. Saprolite permeability is principally intergranular and is the result of dissolution of silicate phases during isovolumetric weathering. A secondary to primary hydraulic conductivity ratio of qs/qp = 150 in the Panola bedrock results in kinetically controlled plagioclase dissolution but thermodynamically inhibited K-feldspar reaction. This result is in accord with calculated chemical saturation states for groundwater sampled in the Panola Granite. In contrast, greater secondary conductivities in the Davis Run saprolite, qs/qp = 800, produces both kinetically controlled plagioclase and K-feldspar dissolution. Faster plagioclase reaction, leading to bedrock weathering in the Panola Granite but not at Davis Run, is attributed to a higher anorthite component of the plagioclase and a wetter and warmer climate. In addition, the Panola Granite has an abnormally high content of disseminated calcite, the dissolution of which precedes the plagioclase weathering front, thus creating additional secondary permeability. Copyright ?? 2001 Elsevier Science Ltd.

Geochimica et Cosmochimica Acta↗

Grazing and ecosystem service delivery in global drylands

Grazing represents the most extensive use of land worldwide. Yet its impacts on ecosystem services remain uncertain because pervasive interactions between grazing pressure, climate, soil properties, and biodiversity may occur but have never been addressed simultaneously. Using a standardized survey at 98 sites across six continents, we show that interactions between grazing pressure, climate, soil, and biodiversity are critical to explain the delivery of fundamental ecosystem services across drylands worldwide. Increasing grazing pressure reduced ecosystem service delivery in warmer and species-poor drylands, whereas positive effects of grazing were observed in colder and species-rich areas. Considering interactions between grazing and local abiotic and biotic factors is key for understanding the fate of dryland ecosystems under climate change and increasing human pressure.

Science↗

The distribution of uranium over Europe: Geological and environmental significance

The variation of baseline levels of uranium in soil and stream sediments over Europe is described, based on new data prepared by the Forum of European Geological Surveys (FOREGS). The samples have been collected and analysed according to the protocols established for the International Union of Geological Sciences/International Association of Geochemistry and Cosmochemistry (IUGS/IAGC) Working Group on Global Geochemical Baselines. The baseline levels of U vary between 0??21 to 53 mg kg-1 in topsoils, 0??19 to 30 mg kg-1 in subsoils and < 1 to 59 mg kg-1 in stream sediments. There is generally good agreement between the levels of U in the three sample types, and the median concentration in all three media is approximately 2 mg kg-1. The most anomalous baseline levels occur over the Variscan orogen, especially areas into which late post-orogenic radiothermal high heat production (HHP) granites were emplaced. Spiderdiagrams based on trace element levels and rare earth element (REE) plots, confirm the association between the highest U anomalies and evolved radiothermal granites. High values are also associated with parts of the Alpine terrain especially in Slovenia, where there are historical U workings, and Southern Italy, where high values of U reflect contemporary volcanism. In contrast, much of the Caledonides of North West Europe and the Precambrian of the Baltic Shield and East European craton and its overlying sedimentary cover have very low values, generally < 4 mg kg-1. The results suggest that the main concern for the environment and human health from U, and the Th and K with which it is generally associated, is the naturally occurring total gamma radiation and radon potential associated with radiothermal granites. This is likely to be especially important where the granites are mineralised and have been worked historically, for example in the North West of the Iberian Peninsula where U and its decay products are likely to be more dispersed in the surface environment. The study also indicates the value of multi-element data in distinguishing between anthropogenic and naturally occurring anomalies.

Transactions of the Institution of Mining and Meta↗

Sources and characteristics of organic matter in the Clackamas River, Oregon, related to the formation of disinfection by-products in treated drinking water

This study characterized the amount and quality of organic matter in the Clackamas River, Oregon, to gain an understanding of sources that contribute to the formation of chlorinated and brominated disinfection by-products (DBPs), focusing on regulated DBPs in treated drinking water from two direct-filtration treatment plants that together serve approximately 100,000 customers. The central hypothesis guiding this study was that natural organic matter leaching out of the forested watershed, in-stream growth of benthic algae, and phytoplankton blooms in the reservoirs contribute different and varying proportions of organic carbon to the river. Differences in the amount and composition of carbon derived from each source affects the types and concentrations of DBP precursors entering the treatment plants and, as a result, yield varying DBP concentrations and species in finished water. The two classes of DBPs analyzed in this study-trihalomethanes (THMs) and haloacetic acids (HAAs)-form from precursors within the dissolved and particulate pools of organic matter present in source water. The five principal objectives of the study were to (1) describe the seasonal quantity and character of organic matter in the Clackamas River; (2) relate the amount and composition of organic matter to the formation of DBPs; (3) evaluate sources of DBP precursors in the watershed; (4) assess the use of optical measurements, including in-situ fluorescence, for estimating dissolved organic carbon (DOC) concentrations and DBP formation; and (5) assess the removal of DBP precursors during treatment by conducting treatability "jar-test" experiments at one of the treatment plants. Data collection consisted of (1) monthly sampling of source and finished water at two drinking-water treatment plants; (2) event-based sampling in the mainstem, tributaries, and North Fork Reservoir; and (3) in-situ continuous monitoring of fluorescent dissolved organic matter (FDOM), turbidity, chlorophyll- a , and other constituents to continuously track source-water conditions in near real-time. Treatability tests were conducted during the four event-based surveys to determine the effectiveness of coagulant and powdered activated carbon (PAC) on the removal of DBP precursors. Sample analyses included DOC, total particulate carbon (TPC), total and dissolved nutrients, absorbance and fluorescence spectroscopy, and, for regulated DBPs, concentrations of THMs and HAAs in finished water and laboratory-based THM and HAA formation potentials (THMFP and HAAFP, respectively) for source water and selected locations throughout the watershed. The results of this study may not be typical given the record and near record amounts of precipitation that occurred during spring that produced streamflow much higher than average in 2010-11. Although there were algal blooms, lower concentrations of chlorophyll- a were observed in the water column during the study period compared to historical data. Concentrations of DBPs in finished (treated) water averaged 0.024 milligrams per liter (mg/L) for THMs and 0.022 mg/L for HAAs; maximum values were about 0.040 mg/L for both classes of DBPs. Although DBP concentrations were somewhat higher within the distribution system, none of the samples collected for this study or for the quarterly compliance monitoring by the water utilities exceeded levels permissible under existing U.S. Environmental Protection Agency (USEPA) regulations: 0.080 mg/L for THMs and 0.060 mg/L for HAAs. DOC concentrations were generally low in the Clackamas River, typically about 1.0-1.5 mg/L. Concentrations in the mainstem occasionally increased to nearly 2.5 mg/L during storms; DOC concentrations in tributaries were sometimes much higher (up to 7.8 mg/L). The continuous in-situ FDOM measurements indicated sharp rises in DOC concentrations in the mainstem following rainfall events; concentrations were relatively stable during summer base flow. Even though the first autumn storm mobilized appreciable quantities of carbon, higher concentrations of DBPs in finished water were observed 3-weeks later, after the ground was saturated from additional rainfall. The majority of the DOC in the lower Clackamas River appears to originate from the upper basin, suggesting terrestrial carbon was commonly the dominant source. Lower-basin tributaries typically contained the highest concentrations of DOC and DBP precursors and contributed substantially to the overall loads in the mainstem during storms. During low-flow periods, tributaries were not major sources of DOC or DBP precursors to the Clackamas River. Although the dissolved fraction of organic carbon contributed the majority of DBP precursors, at times the particulate fraction (inorganic sediment and organic particles including detritus and algal material) contributed a substantial fraction of DBP precursors. Considering just the main-stem sites, on average, 10 percent of THMFP and 32 percent of HAAFP were attributed to particulate carbon. This finding suggests water-treatment methods that remove particles prior to chlorination would reduce finished-water DBP concentrations to some degree. Overall, concentrations of THM and HAA precursors were closely linked to DOC concentrations; laboratory DBP formation potentials (DBPFPs) clearly showed that THMFP and HAAFP were greatest in the downstream tributaries that contained elevated carbon concentrations. However, carbon-normalized "specific" formation potentials for THMs and HAAs (STHMFP and SHAAFP, respectively) revealed changes in carbon character over time that affected the two types of DBP classes differently. HAA precursors were elevated in waters containing aromatic-rich soil-derived material arising from forested areas. In contrast, THM precursors were associated with carbon having a lower aromatic content; highest STHMFP occurred in autumn 2011 in the mainstem from North Fork Reservoir downstream to LO DWTP. This pattern suggests the potential for a link between THM precursors and algal-derived carbon. The highest STHMFP value was measured within North Fork Reservoir, indicating reservoir derived carbon may be important for this class of DBPs. Weak correlations between STHMFP and SHAAFP emphasize that precursor sources for these types of DBPs may be different. This highlights not only that different locations within the watershed produce carbon with different reactivity (specific DBPFP), but also that different management approaches for each class of DBP precursors could be required for control. Treatability tests conducted on source water during four basin-wide surveys demonstrated that an average of about 40 percent of DOC can be removed by coagulation. While the decrease in THMFP following coagulation was similar to DOC, the decrease in HAAFP was much greater (approximately 70 percent), indicating coagulation is particularly effective at removing HAA precursors'likely because of the aromatic nature of the carbon associated with HAA precursors. Several findings from this study have direct implications for managing drinking-water resources and for providing useful information that may help improve treatment-plant operations. For example, the use of in-situ fluorometers that measure FDOM provided an excellent proxy for DOC concentration in this system and revealed short-term, rapid changes in DOC concentration during storm events. In addition, the strong correlation between FDOM values measured in-situ and HAA5 concentrations in finished water may permit estimation of continuous HAA concentrations, as was done here. As part of this study, multiple in-situ FDOM sensors were deployed continuously and in real-time to characterize the composition of dissolved organic matter. Although the initial results were promising, additional research and engineering developments will be needed to demonstrate the full utility of these sensors for this purpose. In conclusion, although DBPFPs were strongly correlated to DOC concentration, some DBPs formed from particulate carbon, including terrestrial leaf material and algal material such as planktonic species of blue-green algae and sloughed filaments, stalks, and cells of benthic algae. Different precursor sources in the watershed were evident from the data, suggesting specific actions may be available to address some of these sources. In-situ measurements of FDOM proved to be an excellent proxy for DOC concentration as well as HAA formation during treatment, which suggests further development and refinement of these sensors have the potential to provide real-time information about complex watershed processes to operators at the drinking-water treatment plants. Follow-up studies could examine the relative roles that terrestrial and algal sources have on the DBP precursor pool to better understand how watershed-management activities may be affecting the transport of these compounds to Clackamas River drinking-water intakes. Given the low concentrations of algae in the water column during this study, additional surveys during more typical river conditions could provide a more complete understanding of how algae contribute DBP precursors. Further development of FDOM-sensor technology can improve our understanding of carbon dynamics in the river and how concentrations may be trending over time. This study was conducted in collaboration with Clackamas River Water and the City of Lake Oswego water utilities. Other research partners included Oregon Health and Science University in Hillsboro, Oregon, Alexin Laboratory in Tigard, Oregon, U.S. Geological Survey National Research Program Laboratory in Denver, Colorado, and the U.S. Geological Survey Water Science Centers in Portland, Oregon, and Sacramento, California. This project was supported with funding from Clackamas River Water, City of Lake Oswego, the U.S. Geological Survey, and the Water Research Foundation.

Oregon↗

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↗

Streams, springs, and volcanic lakes for volcano monitoring

Introduction Volcanic unrest can trigger appreciable change to surface waters such as streams, springs, and volcanic lakes. Magma degassing produces gases and soluble salts that are absorbed into groundwater that feeds streams and lakes. As magma ascends, the amount of heat and degassing will increase, and so will any related geochemical and thermal signal. Subsurface magma movement can cause pressurization that alters hydrostatic head and may induce groundwater discharge. Fluid-pressure changes have been linked to distal volcano-tectonic earthquakes (White and McCausland, 2016; Coulon and others, 2017) and phreatic eruptions (for example, Yamaoka and others, 2016). Clearly, changes in groundwater and surface waters are both indicators of unrest and clues to how and where magma is rising toward the surface. Where possible, it is prudent to incorporate real-time hydrologic data into multiparameter monitoring of restless volcanoes. Hydrologic dynamics can also be tracked by changes in groundwater levels that are commonly measured in shallow boreholes (see of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–K, 5 p., https://doi.org/10.3133/sir20245062k. ">chapter K , this volume, on boreholes; Hurwitz and Lowenstern, 2024). Although inferred to be common, relatively few volcano-hydrology anomalies are well documented, and many are essentially anecdotal (Newhall and others, 2001), reflecting the fact that high-resolution time series remain rare. Extreme examples include the 2008 eruption of Nevado del Huila, Colombia, where relatively minor phreatomagmatic eruptions were accompanied by expulsion of as much as 300 million cubic meters of groundwater from fissures high on the volcano (Worni and others, 2011), generating large lahars. Substantial decreases in flow rate from springs about 8 kilometers from the summit of Mayon Volcano, Philippines, have been noted before most eruptions in the 20th century (Newhall and others, 2001). Stream monitoring at Redoubt Volcano in 2009 allowed Werner and others (2012) to recognize that groundwater was unable to absorb (or scrub) the high flux of volcanic gas and that a high CO 2 /SO 2 precursor signal had been evident for 5 months prior to the eruption. A key to better interpreting hydrologic anomalies—or even identifying them—is therefore obtaining adequate baseline data. Most hydrologic monitoring at U.S. volcanoes has been accomplished by intermittent sampling surveys with annual or less frequent sampling (for example, https://hotspringchem.wr.usgs.gov/index.php ). More frequent sampling, however, generally is needed to establish reliable baselines. A recent hydrologic and hydrothermal monitoring experiment at 25 sites and 10 of the 12 level 4 (very high threat) volcanoes in the U.S. portion of the Cascade Range demonstrated that there is sufficient temporal variability in hydrothermal fluxes, even during quiescent periods, that one-time measurements will commonly have limited interpretive value (Crankshaw and others, 2018). Thus, surveys are best augmented with data from streamgages (for example, Evans and others, 2004; Bergfeld and others, 2008). Streamflow (water discharge) data allow measured temperature and specific conductance to be converted to heat and solute mass fluxes, which could be insightful parameters for detecting anomalous activity (McCleskey and others, 2012). At the Yellowstone Caldera, long-term monitoring of river solutes has allowed calculation of the chloride flux, a proxy for heat discharge (Hurwitz and others, 2007; McCleskey and others, 2016) from the subsurface magma. This is readily accomplished because data from streamgages are continuously recorded and archived by the U.S. Geological Survey (USGS) National Water Information System (NWIS) (USGS, 2024). Similar studies on stratovolcanoes or shield volcanoes would be scientifically useful, and yet are logistically challenging, requiring streamgages on numerous radial drainages complemented by either frequent manual sampling or numerous deployments of equipment to measure water temperature and specific conductance as a proxy for water chemistry. Another challenge is that some volcanic areas, especially shield volcanoes, are characterized by near-surface porous rocks and soils, such that surface streams are rare and replaced by distant, dilute large-volume springs with only a trace of any original volcanically sourced water (Manga, 2001; Hurwitz and others, 2021). Volcanic lakes are worthy of special attention for monitoring efforts, as their temperature and composition can provide evidence of increased flux of volatile-rich fluids from below. Quantifying changes in volatile and heat release from magma can be simpler in lakes than for volcanoes with radial drainages and no major lakes. Moreover, volcanic lakes pose a range of hazards themselves, including phreatomagmatic eruptions, debris flows, flank collapse, tsunamis, and toxic gas release (Mastin and Witter, 2000; Delmelle and others, 2015; Manville, 2015; Rouwet and others, 2015)—hazards that have historically been responsible for substantial loss of life at many volcanoes worldwide (Manville, 2015). Catastrophic CO 2 release at Lake Nyos, Cameroon, in 1986 suffocated about 1,750 people and about 3,500 livestock and was probably triggered by a large landslide into the gas-saturated lake (Kling and others, 1987; Evans and others, 1993). Gas-charged springs in Soda Bay within Clear Lake (California) have caused almost a dozen deaths to bathers in the past hundred years (ABC News, 2000). A 2005 example of lake overturn and abundant gas release was documented at Mount Chiginagak in Alaska (Schaefer and others, 2008) but did not result in any human casualties. Although thermally stratified lakes, which promote trapping of exsolved magmatic gas, tend to develop in tropical regions, the phenomenon can also arise where salinity creates meromixis (a condition in which a lake does not mix completely), as occurs in Mono Lake, California (Jellison and Melack, 1993; Jellison and others, 1998). If magma erupts or flows into a lake, the interaction between hot magma and cold water can be explosive (Mastin and others, 2004; Zimanowski and others, 2015) and substantially expand the area affected by the eruption. Another hazard is the breaching of crater rims by landslides triggered by volcanic and (or) seismic activity. Under some circumstances, substantial volumes of water can be displaced, leading to large floods and lahars. Late Holocene lake flooding from Aniakchak Crater in the Alaska Peninsula (Waythomas, 2022) and from Paulina Lake in Newberry Crater, Oregon (Chitwood and Jensen, 2000), caused by the failure of outlet sills, testify to the substantial hazards at lake-filled calderas. Several volcanic systems in the United States host lakes known to receive heat and gas from underlying magma. These lakes vary widely in area, depth, and chemical composition. Lakes are present at level 4 volcanoes, including Crater Lake and Newberry Volcano in Oregon; Yellowstone Caldera in Wyoming; Long Valley Caldera, Clear Lake volcanic field, Medicine Lake, and Salton Buttes in California; and Aniakchak Crater, Mount Katmai, Fisher Caldera, Mount Okmok, and Kaguyak Crater, among others, in Alaska. A water lake was present in Halemaʻumaʻu, the crater of Kīlauea, Hawai‘i (fig. F1), from October 2019 to December 2020. Level 3 volcanoes with lakes include Mono Lake volcanic field (Calif.), Mount Bachelor (Ore.), Ukinrek Maars and Mount Chiginagak (Alaska), and Soda Lake (Nevada). In addition, there are lakes at many levels 1 and 2 volcanoes. In the United States, there are no strongly acidic lakes that receive abundant input of magmatic gas, such as those found at Mount Ruapehu (New Zealand), Ijen and Kelud (Indonesia), and Poás (Costa Rica). Nevertheless, many contain fluids that provide clues to magmatic processes below. Since publication of a previous report on recommended instrumentation for volcano monitoring (Moran and others, 2008), continuous hydrologic monitoring has become increasingly feasible. However, changes in water pressure, temperature, and chemistry remain, in general, poorly studied phenomena at volcanoes (Sparks, 2003; National Academies of Sciences, Engineering, and Medicine, 2017). Recent efforts by the USGS have included the temporary study of Cascade Range volcanoes, which included frequent (15 minute to hourly) temporal sampling of temperature, depth, and conductivity (Crankshaw and others, 2018; Ingebritsen and Evans, 2019). At Yellowstone Caldera, many streamgages have now added thermistors and specific conductance sensors, allowing estimation of time-dependent chloride flux as a proxy for variations in subsurface heat flux (McCleskey and others, 2012, 2016). Efforts to better understand lakes have also accelerated, with bathymetric mapping and sampling carried out at several locations in the United States. Especially thorough work was done at Yellowstone Lake thanks to the Hydrothermal Dynamics of Yellowstone Lake (HD-YLAKE, https://hdylake.org ) project, funded primarily by the National Science Foundation. In addition to geophysical surveys and recovery of cores and other samples, HD-YLAKE investigations included remotely operated vehicle (ROV) investigations of hydrothermal vents on the lake floor (fig. F2). Data collected by the ROV provided a better understanding of the thermal and chemical influx from lake-bottom hydrothermal systems (Sohn and others, 2017). In this chapter, we focus on detecting changes in the chemistry, temperature, discharge, or water levels of streams, springs, and lakes that can be caused by seismicity, volumetric strains, or increases in gas flux associated with ascending magma. There is unavoidable overlap with other chapters of this report. Samples of water and gas can also be obtained in boreholes ( of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–K, 5 p., https://doi.org/10.3133/sir20245062k. ">chapter K , this volume; Hurwitz and Lowenstern, 2024), both shallow and deep. Gas monitoring ( of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–E, 11 p., https://doi.org/10.3133/sir20245062e.">chapter E , this volume; Lewicki and others, 2024) relies in part on samples from springs and wells, particularly where measurable gas plumes are absent. Water acts as a trigger and lubricant for landslides and sediment-rich floods, and so hydrology has obvious relevance for lahar monitoring, as discussed in of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–H, 6 p., https://doi.org/10.3133/sir20245062h. ">chapter H (this volume; Thelen and others, 2024). Shared situational awareness among scientists engaged in geophysical, gas, and hydrologic monitoring will improve overall understanding of the volcanic hazard.

Scientific Investigations Report↗

FLUXNET-CH4: A global, multi-ecosystem database and analysis of methane seasonality from freshwater wetlands

Methane (CH 4 ) emissions from natural landscapes constitute roughly half of global CH 4 contributions to the atmosphere, yet large uncertainties remain in the absolute magnitude and the seasonality of emission quantities and drivers. Eddy covariance (EC) measurements of CH 4 flux are ideal for constraining ecosystem-scale CH 4 emissions due to quasi-continuous and high-temporal-resolution CH 4 flux measurements, coincident carbon dioxide, water, and energy flux measurements, lack of ecosystem disturbance, and increased availability of datasets over the last decade. Here, we (1) describe the newly published dataset, FLUXNET-CH 4 Version 1.0, the first open-source global dataset of CH 4 EC measurements (available at https://fluxnet.org/data/fluxnet-ch4-community-product/ , last access: 7 April 2021). FLUXNET-CH 4 includes half-hourly and daily gap-filled and non-gap-filled aggregated CH 4 fluxes and meteorological data from 79 sites globally: 42 freshwater wetlands, 6 brackish and saline wetlands, 7 formerly drained ecosystems, 7 rice paddy sites, 2 lakes, and 15 uplands. Then, we (2) evaluate FLUXNET-CH 4 representativeness for freshwater wetland coverage globally because the majority of sites in FLUXNET-CH 4 Version 1.0 are freshwater wetlands which are a substantial source of total atmospheric CH 4 emissions; and (3) we provide the first global estimates of the seasonal variability and seasonality predictors of freshwater wetland CH 4 fluxes. Our representativeness analysis suggests that the freshwater wetland sites in the dataset cover global wetland bioclimatic attributes (encompassing energy, moisture, and vegetation-related parameters) in arctic, boreal, and temperate regions but only sparsely cover humid tropical regions. Seasonality metrics of wetland CH 4 emissions vary considerably across latitudinal bands. In freshwater wetlands (except those between 20 ∘ S to 20 ∘ N) the spring onset of elevated CH 4 emissions starts 3 d earlier, and the CH 4 emission season lasts 4 d longer, for each degree Celsius increase in mean annual air temperature. On average, the spring onset of increasing CH 4 emissions lags behind soil warming by 1 month, with very few sites experiencing increased CH 4 emissions prior to the onset of soil warming. In contrast, roughly half of these sites experience the spring onset of rising CH 4 emissions prior to the spring increase in gross primary productivity (GPP). The timing of peak summer CH 4 emissions does not correlate with the timing for either peak summer temperature or peak GPP. Our results provide seasonality parameters for CH 4 modeling and highlight seasonality metrics that cannot be predicted by temperature or GPP (i.e., seasonality of CH 4 peak). FLUXNET-CH 4 is a powerful new resource for diagnosing and understanding the role of terrestrial ecosystems and climate drivers in the global CH 4 cycle, and future additions of sites in tropical ecosystems and site years of data collection will provide added value to this database. All seasonality parameters are available at https://doi.org/10.5281/zenodo.4672601 (Delwiche et al., 2021). Additionally, raw FLUXNET-CH 4 data used to extract seasonality parameters can be downloaded from https://fluxnet.org/data/fluxnet-ch4-community-product/ (last access: 7 April 2021), and a complete list of the 79 individual site data DOIs is provided in Table 2 of this paper.

Earth System Science Data↗

Isotopic composition and origin of indigenous natural perchlorate and co-occurring nitrate in the southwestern United States

Perchlorate (ClO 4 − ) has been detected widely in groundwater and soils of the southwestern United States. Much of this ClO 4 − appears to be natural, and it may have accumulated largely through wet and dry atmospheric deposition. This study evaluates the isotopic composition of natural ClO 4 − indigenous to the southwestern U.S. Stable isotope ratios were measured in ClO 4 − (δ 18 O, Δ 17 O, δ 37 Cl) and associated NO 3 − (δ 18 O, Δ 17 O, δ 15 N) in groundwater from the southern High Plains (SHP) of Texas and New Mexico and the Middle Rio Grande Basin (MRGB) in New Mexico, from unsaturated subsoil in the SHP, and from NO 3 − -rich surface caliche deposits near Death Valley, California. The data indicate natural ClO 4 − in the southwestern U.S. has a wide range of isotopic compositions that are distinct from those reported previously for natural ClO 4 − from the Atacama Desert of Chile as well as all known synthetic ClO 4 − . ClO 4 − in Death Valley caliche has a range of high Δ 17 O values (+8.6 to +18.4 ‰), overlapping and extending the Atacama range, indicating at least partial atmospheric formation via reaction with ozone (O 3 ). However, the Death Valley δ 37 Cl values (−3.1 to −0.8 ‰) and δ 18 O values (+2.9 to +26.1‰) are higher than those of Atacama ClO 4 − . In contrast, ClO 4 − from western Texas and New Mexico has much lower Δ 17 O (+0.3 to +1.3‰), with relatively high δ 37 Cl (+3.4 to +5.1 ‰) and δ 18 O (+0.5 to +4.8 ‰), indicating either that this material was not primarily generated with O 3 as a reactant or that the ClO 4 − was affected by postdepositional O isotope exchange. High Δ 17 O values in ClO 4 − (Atacama and Death Valley) are associated with high Δ 17 O values in NO 3 − , indicating that both compounds preserve characteristics of O 3 -related atmospheric production in hyper-arid settings, whereas both compounds have low Δ 17 O values in less arid settings. Although Δ 17 O variations in terrestrial NO 3 − can be attributed to mixing of atmospheric (high Δ 17 O) and biogenic (low Δ 17 O) NO 3 − , variations in Δ 17 O of terrestrial ClO 4 − are not readily explained in the same way. This study provides important new constraints for identifying natural sources of ClO 4 − in different environments by multicomponent isotopic characteristics, while presenting the possibilities of divergent ClO 4 − formation mechanisms and(or) ClO 4 − isotopic exchange in biologically active environments.

Environmental Science & Technology↗

The Boreal-Arctic Wetland and Lake Dataset (BAWLD)

Methane emissions from boreal and arctic wetlands, lakes, and rivers are expected to increase in response to warming and associated permafrost thaw. However, the lack of appropriate land cover datasets for scaling field-measured methane emissions to circumpolar scales has contributed to a large uncertainty for our understanding of present-day and future methane emissions. Here we present the Boreal–Arctic Wetland and Lake Dataset (BAWLD), a land cover dataset based on an expert assessment, extrapolated using random forest modelling from available spatial datasets of climate, topography, soils, permafrost conditions, vegetation, wetlands, and surface water extents and dynamics. In BAWLD, we estimate the fractional coverage of five wetland, seven lake, and three river classes within 0.5 × 0.5∘ grid cells that cover the northern boreal and tundra biomes (17 % of the global land surface). Land cover classes were defined using criteria that ensured distinct methane emissions among classes, as indicated by a co-developed comprehensive dataset of methane flux observations. In BAWLD, wetlands occupied 3.2 × 106 km2 (14 % of domain) with a 95 % confidence interval between 2.8 and 3.8 × 106 km2. Bog, fen, and permafrost bog were the most abundant wetland classes, covering ∼ 28 % each of the total wetland area, while the highest-methane-emitting marsh and tundra wetland classes occupied 5 % and 12 %, respectively. Lakes, defined to include all lentic open-water ecosystems regardless of size, covered 1.4 × 106 km2 (6 % of domain). Low-methane-emitting large lakes (>10 km2) and glacial lakes jointly represented 78 % of the total lake area, while high-emitting peatland and yedoma lakes covered 18 % and 4 %, respectively. Small (<0.1 km2) glacial, peatland, and yedoma lakes combined covered 17 % of the total lake area but contributed disproportionally to the overall spatial uncertainty in lake area with a 95 % confidence interval between 0.15 and 0.38 × 106 km2. Rivers and streams were estimated to cover 0.12 × 106 km2 (0.5 % of domain), of which 8 % was associated with high-methane-emitting headwaters that drain organic-rich landscapes. Distinct combinations of spatially co-occurring wetland and lake classes were identified across the BAWLD domain, allowing for the mapping of “wetscapes” that have characteristic methane emission magnitudes and sensitivities to climate change at regional scales. With BAWLD, we provide a dataset which avoids double-accounting of wetland, lake, and river extents and which includes confidence intervals for each land cover class. As such, BAWLD will be suitable for many hydrological and biogeochemical modelling and upscaling efforts for the northern boreal and arctic region, in particular those aimed at improving assessments of current and future methane emissions. Data are freely available at https://doi.org/10.18739/A2C824F9X (Olefeldt et al., 2021).

Earth System Science Data↗

The Missoula and Bonneville floods—A review of ice-age megafloods in the Columbia River basin

The Channeled Scabland of eastern Washington State, USA, brought megafloods to the scientific forefront. A 30,000-km2 landscape of coulees and cataracts carved into the region’s loess-covered basalt attests to overwhelming volumes of energetic water. The scarred landscape, garnished by huge boulder bars and far-travelled ice-rafted erratics, spurred J Harlen Bretz’s vigorously disputed flood hypothesis in the 1920s. First known as the Spokane flood, it was rebranded the Missoula flood once understood that the water came from glacial Lake Missoula, formed when the Purcell Trench lobe of the last-glacial Cordilleran ice sheet dammed the Clark Fork valley in northwestern Idaho with ice a kilometer thick. Bretz’s flood evidence in the once-remote Channeled Scabland, widely seen and elaborated by the 1950s, eventually swayed consensus for cataclysmic flooding. Missoula flood questions then turned to some that continue today: how many? when? how big? what routes? what processes? The Missoula floods passed through eastern Washington by a multitude of valleys, coulees and scabland tracts, some contemporaneously, some sequentially. Which routings and their timing depended on the positions of various lobes of the multi-pronged Cordilleran ice sheet and the erosional development of the channels themselves. The first floods mostly followed the big bend of Columbia valley looping through north-central Washington. But the south-advancing Okanogan ice lobe soon blocked that path, forming long-lasting glacial Lake Columbia in the impounded Columbia valley. Missoula floods into this lake were diverted south out of the Columbia valley and into eastern Washington coulees and scabland tracts. At least four floods entered Moses Coulee, but then as the Okanogan lobe advanced over and blocked the head of that coulee, more eastern paths took the water, including Grand Coulee and the Telford-Crab-Creek and Cheney-Palouse scabland tracts. Flood routing also depended on the erosion of the coulees. At some point, headward erosion of upper Grand Coulee lowered the divide saddle between the west-running Columbia valley and the deep and wide Grand Coulee heading southwest. Still uncertain is when this happened and the consequences with respect to the stage and extent of glacial Lake Columbia and to flood access to the other, higher, flood routes. Downstream, all flood routes converged onto Pasco Basin, flowed through Wallula Gap and the Columbia River Gorge into the Pacific Ocean, following submarine canyons and depositing sediment layers on abyssal plains. Stratigraphic studies indicate dozens—likely more than a hundred—separate Missoula floods during the last glacial period. Over the length of the flood route, backwater areas and depositional basins preserve multiple flood beds, many of which are separated by signs of time, including volcanic ash layers and soil development in subaerial environments; and varve-like beds and pelagic mud layers in lacustrine and marine settings. Evidence also comes from the glacial Lake Missoula basin, where stratigraphy indicates dozens of filling and emptying cycles. Varve counts in conjunction of radiocarbon dating and paleomagnetic secular variation show the repeated filling-and-release cycles of glacial Lake Missoula had intervals possibly as long as 100 years early in the lake’s history but diminished to just one or two years for the last few floods. This behavior accords with jökulhlaup-style floods released by subglacial drainage from a self-dumping ice-dammed lake. But not yet clear is whether such a mechanism applies to all the floods or if some emptied more cataclysmically as hypothesized by some. Radiocarbon dating of sparse organic materials remains key to defining flood chronology but has been lately bolstered by analyses of terrestrial cosmogenic nuclides and optically stimulated luminescence. Varve counts and paleomagnetic secular variation studies help to define durations and intervals represented by sequences of flood beds. The ~16 ka Mount St. Helens Set S tephra is commonly interbedded within flood deposits, enabling correlation of deposits among sites. Tephra from the 13.7–13.4 ka eruption of Glacier Peak overlies all glacial Lake Missoula and Missoula flood deposits, defining an end time. Overall conclusions are that glacial Lake Missoula was extant and producing floods for at least 3–4 ky during 20–14 ka. At least ~75 floods preceded Mount St Helens Set S, followed by 30 or more after the tephra fall. Most floods entered glacial Lake Columbia, impounded by the Okanogan lobe, for 2–5 ky between about 18.5 and 15 ka. Glacial Lake Columbia outlived Lake Missoula by >200–400 yr but may have been born later since at least one flood came down the Columbia valley before the Okanogan ice lobe blocked the Columbia valley at 18.5–18 ka. The maximum extent of the Okanogan and Purcell Trench lobes, many Missoula floods, substantial erosion of upper Grand Coulee, and the widespread tephra falls from Mount St. Helens eruptions all happened about 17–15 ka. People, in the area since 16.6–15.3 ka, almost certainly witnessed the last of the Missoula floods and later large floods from other ice-dammed lakes in the Columbia River basin. Quantitative flow analyses give peak discharge estimates and support understanding of erosional and depositional processes. The first flow assessments were simple cross-section calculations but recent assessments employ two-dimensional hydrodynamic models. The general finding is that emplacement of the maximum stage evidence requires about 20 million m3/s near the Lake Missoula outlet and about 5–15 million m3/s through Wallula Gap and downstream in the Columbia River Gorge. These hydraulic analyses raise still-unresolved questions regarding canyon erosion and possible additional water sources. The large Pleistocene Bonneville flood entered the Columbia River system from the southeast from pluvial Lake Bonneville, the Pleistocene predecessor to Great Salt Lake in the eastern Great Basin. During the last glacial, the lake basin filled, covering >50,000 km2 with 10,400 km3 of water before reaching its maximum possible stage governed by Red Rock Pass, the lowest divide separating the basin from the Snake River basin to the north. The overtopping lake rapidly incised 108–125 m into the Red Rock Pass outlet, spilling half of its total lake volume. G.K. Gilbert described the essential sequence in the 1870s, but the flood was mostly forgotten until the late 1950s when Harold Malde linked the spectacular scabland topography and bouldery “melon gravel” on the Snake River Plain to the Lake Bonneville overflow. The Bonneville flood appears to have been a singular event at about 18 ka. No evidence of multiple or pre-last-glacial spillovers has yet been found. Its total volume was about twice that of a maximum Lake Missoula flood yet its peak discharge was ~1 million m3/s, less than a tenth of the largest Missoula floods. Its comparatively simple flow path and much steadier flow make the Bonneville flood ideal for new studies of erosional and depositional processes. At least two floods seem to have passed down the Columbia valley after the last of the Missoula floods, including a large flood about ~14 ka likely from cataclysmic demise of the thinning Okanogan ice lobe dam impounding glacial Lake Columbia. Floods from earlier glacial ages left scant yet clear evidence in the Channeled Scabland and Columbia valley. But their source, timing, and magnitudes are little understood. Some deposits are paleomagnetically reversed, thus older than ~800 ka. Last-glacial floods and perhaps older ones affected the Snake River Plain, some likely sourced in lakes dammed by alpine glaciers in central Idaho.

Idaho, Oregon, Washington↗