USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earth 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.

At least 577 records · Page 32Linked to original sources

Wildland fire ash: Production, composition and eco-hydro-geomorphic effects

Fire transforms fuels (i.e. biomass, necromass, soil organic matter) into materials with different chemical and physical properties. One of these materials is ash, which is the particulate residue remaining or deposited on the ground that consists of mineral materials and charred organic components. The quantity and characteristics of ash produced during a wildland fire depend mainly on (1) the total burned fuel (i.e. fuel load), (2) fuel type and (3) its combustion completeness. For a given fuel load and type, a higher combustion completeness will reduce the ash organic carbon content, increasing the relative mineral content, and hence reducing total mass of ash produced. The homogeneity and thickness of the ash layer can vary substantially in space and time and reported average thicknesses range from close to 0 to 50 mm. Ash is a highly mobile material that, after its deposition, may be incorporated into the soil profile, redistributed or removed from a burned site within days or weeks by wind and water erosion to surface depressions, footslopes, streams, lakes, reservoirs and, potentially, into marine deposits. Research on the composition, properties and effects of ash on the burned ecosystem has been conducted on material collected in the field after wildland and prescribed fires as well as on material produced in the laboratory. At low combustion completeness (typically T < 450 °C), ash is organic-rich, with organic carbon as the main component. At high combustion completeness (T > 450 °C), most organic carbon is volatized and the remaining mineral ash has elevated pH when in solution. It is composed mainly of calcium, magnesium, sodium, potassium, silicon and phosphorous in the form of inorganic carbonates, whereas at T > 580 °C the most common forms are oxides. Ash produced under lower combustion completeness is usually darker, coarser, and less dense and has a higher saturated hydraulic conductivity than ash with higher combustion completeness, although physical reactions with CO 2 and when moistened produce further changes in ash characteristics. As a new material present after a wildland fire, ash can have profound effects on ecosystems. It affects biogeochemical cycles, including the C cycle, not only within the burned area, but also globally. Ash incorporated into the soil increases temporarily soil pH and nutrient pools and changes physical properties such as albedo, soil texture and hydraulic properties including water repellency. Ash modifies soil hydrologic behavior by creating a two-layer system: the soil and the ash layer, which can function in different ways depending on (1) ash depth and type, (2) soil type and (3) rainfall characteristics. Key parameters are the ash's water holding capacity, hydraulic conductivity and its potential to clog soil pores. Runoff from burned areas carries soluble nutrients contained in ash, which can lead to problems for potable water supplies. Ash deposition also stimulates soil microbial activity and vegetation growth. Further work is needed to (1) standardize methods for investigating ash and its effects on the ecosystem, (2) characterize ash properties for specific ecosystems and wildland fire types, (3) determine the effects of ash on human and ecosystem health, especially when transported by wind or water, (4) investigate ash's controls on water and soil losses at slope and catchment scales, (5) examine its role in the C cycle, and (6) study its redistribution and fate in the environment.

Earth-Science Reviews↗

Desert wetlands in the geologic record

Desert wetlands support flora and fauna in a variety of hydrologic settings, including seeps, springs, marshes, wet meadows, ponds, and spring pools. Over time, eolian, alluvial, and fluvial sediments become trapped in these settings by a combination of wet ground conditions and dense plant cover. The result is a unique combination of clastic sediments, chemical precipitates, and organic matter that is preserved in the geologic record as ground-water discharge (GWD) deposits. GWD deposits contain information on the timing and magnitude of past changes in water-table levels and, therefore, are a potential source of paleohydrologic and paleoclimatic information. In addition, they can be important archeological and paleontological archives because desert wetlands provide reliable sources of fresh water, and thus act as focal points for human and faunal activities, in some of the world's harshest and driest lands. Here, we review some of the physical, sedimentological, and geochemical characteristics common to GWD deposits, and provide a contextual framework that researchers can use to identify and interpret geologic deposits associated with desert wetlands. We discuss several lines of evidence used to differentiate GWD deposits from lake deposits (they are commonly confused), and examine how various types of microbiota and depositional facies aid in reconstructing past environmental and hydrologic conditions. We also review how late Quaternary GWD deposits are dated, as well as methods used to investigate desert wetlands deeper in geologic time. We end by evaluating the strengths and limitations of hydrologic and climatic records derived from GWD deposits, and suggest several avenues of potential future research to further develop and utilize these unique and complex systems.

Earth-Science Reviews↗

Analyzing high resolution topography for advancing the understanding of mass and energy transfer through landscapes: A review

The study of mass and energy transfer across landscapes has recently evolved to comprehensive considerations acknowledging the role of biota and humans as geomorphic agents, as well as the importance of small-scale landscape features. A contributing and supporting factor to this evolution is the emergence over the last two decades of technologies able to acquire high resolution topography (HRT) (meter and sub-meter resolution) data. Landscape features can now be captured at an appropriately fine spatial resolution at which surface processes operate; this has revolutionized the way we study Earth-surface processes. The wealth of information contained in HRT also presents considerable challenges. For example, selection of the most appropriate type of HRT data for a given application is not trivial. No definitive approach exists for identifying and filtering erroneous or unwanted data, yet inappropriate filtering can create artifacts or eliminate/distort critical features. Estimates of errors and uncertainty are often poorly defined and typically fail to represent the spatial heterogeneity of the dataset, which may introduce bias or error for many analyses. For ease of use, gridded products are typically preferred rather than the more information-rich point cloud representations. Thus many users take advantage of only a fraction of the available data, which has furthermore been subjected to a series of operations often not known or investigated by the user. Lastly, standard HRT analysis work-flows are yet to be established for many popular HRT operations, which has contributed to the limited use of point cloud data. In this review, we identify key research questions relevant to the Earth-surface processes community within the theme of mass and energy transfer across landscapes and offer guidance on how to identify the most appropriate topographic data type for the analysis of interest. We describe the operations commonly performed from raw data to raster products and we identify key considerations and suggest appropriate work-flows for each, pointing to useful resources and available tools. Future research directions should stimulate further development of tools that take advantage of the wealth of information contained in the HRT data and address the present and upcoming research needs such as the ability to filter out unwanted data, compute spatially variable estimates of uncertainty and perform multi-scale analyses. While we focus primarily on HRT applications for mass and energy transfer, we envision this review to be relevant beyond the Earth-surface processes community for a much broader range of applications involving the analysis of HRT.

Earth-Science Reviews↗

A review of single-sample-based models and other approaches for radiocarbon dating of dissolved inorganic carbon in groundwater

Numerous methods have been proposed to estimate the pre-nuclear-detonation 14 C content of dissolved inorganic carbon (DIC) recharged to groundwater that has been corrected/adjusted for geochemical processes in the absence of radioactive decay ( 14 C 0 ) - a quantity that is essential for estimation of radiocarbon age of DIC in groundwater. The models/approaches most commonly used are grouped as follows: (1) single-sample-based models, (2) a statistical approach based on the observed (curved) relationship between 14 C and δ 13 C data for the aquifer, and (3) the geochemical mass-balance approach that constructs adjustment models accounting for all the geochemical reactions known to occur along a groundwater flow path. This review discusses first the geochemical processes behind each of the single-sample-based models, followed by discussions of the statistical approach and the geochemical mass-balance approach. Finally, the applications, advantages and limitations of the three groups of models/approaches are discussed. The single-sample-based models constitute the prevailing use of 14 C data in hydrogeology and hydrological studies. This is in part because the models are applied to an individual water sample to estimate the 14 C age, therefore the measurement data are easily available. These models have been shown to provide realistic radiocarbon ages in many studies. However, they usually are limited to simple carbonate aquifers and selection of model may have significant effects on 14 C 0 often resulting in a wide range of estimates of 14 C ages. Of the single-sample-based models, four are recommended for the estimation of 14 C 0 of DIC in groundwater: Pearson's model, (Ingerson and Pearson, 1964; Pearson and White, 1967), Han & Plummer's model (Han and Plummer, 2013), the IAEA model (Gonfiantini, 1972; Salem et al., 1980), and Oeschger's model (Geyh, 2000). These four models include all processes considered in single-sample-based models, and can be used in different ranges of 13 C values. In contrast to the single-sample-based models, the extended Gonfiantini & Zuppi model (Gonfiantini and Zuppi, 2003; Han et al., 2014) is a statistical approach. This approach can be used to estimate 14 C ages when a curved relationship between the 14 C and 13 C values of the DIC data is observed. In addition to estimation of groundwater ages, the relationship between 14 C and δ 13 C data can be used to interpret hydrogeological characteristics of the aquifer, e.g. estimating apparent rates of geochemical reactions and revealing the complexity of the geochemical environment, and identify samples that are not affected by the same set of reactions/processes as the rest of the dataset. The investigated water samples may have a wide range of ages, and for waters with very low values of 14 C, the model based on statistics may give more reliable age estimates than those obtained from single-sample-based models. In the extended Gonfiantini & Zuppi model, a representative system-wide value of the initial 14 C content is derived from the 14 C and δ 13 C data of DIC and can differ from that used in single-sample-based models. Therefore, the extended Gonfiantini & Zuppi model usually avoids the effect of modern water components which might retain ‘bomb’ pulse signatures. The geochemical mass-balance approach constructs an adjustment model that accounts for all the geochemical reactions known to occur along an aquifer flow path (Plummer et al., 1983; Wigley et al., 1978; Plummer et al., 1994; Plummer and Glynn, 2013), and includes, in addition to DIC, dissolved organic carbon (DOC) and methane (CH 4 ). If sufficient chemical, mineralogical and isotopic data are available, the geochemical mass-balance method can yield the most accurate estimates of the adjusted radiocarbon age. The main limitation of this approach is that complete information is necessary on chemical, mineralogical and isotopic data and these data are often limited. Failure to recognize the limitations and underlying assumptions on which the various models and approaches are based can result in a wide range of estimates of 14 C 0 and limit the usefulness of radiocarbon as a dating tool for groundwater. In each of the three generalized approaches (single-sample-based models, statistical approach, and geochemical mass-balance approach), successful application depends on scrutiny of the isotopic ( 14 C and 13 C) and chemical data to conceptualize the reactions and processes that affect the 14 C content of DIC in aquifers. The recently developed graphical analysis method is shown to aid in determining which approach is most appropriate for the isotopic and chemical data from a groundwater system.

Earth-Science Reviews↗

Trace elements at the intersection of marine biological and geochemical evolution

Life requires a wide variety of bioessential trace elements to act as structural components and reactive centers in metalloenzymes. These requirements differ between organisms and have evolved over geological time, likely guided in some part by environmental conditions. Until recently, most of what was understood regarding trace element concentrations in the Precambrian oceans was inferred by extrapolation, geochemical modeling, and/or genomic studies. However, in the past decade, the increasing availability of trace element and isotopic data for sedimentary rocks of all ages has yielded new, and potentially more direct, insights into secular changes in seawater composition – and ultimately the evolution of the marine biosphere. Compiled records of many bioessential trace elements (including Ni, Mo, P, Zn, Co, Cr, Se, and I) provide new insight into how trace element abundance in Earth's ancient oceans may have been linked to biological evolution. Several of these trace elements display redox-sensitive behavior, while others are redox-sensitive but not bioessential (e.g., Cr, U). Their temporal trends in sedimentary archives provide useful constraints on changes in atmosphere-ocean redox conditions that are linked to biological evolution, for example, the activity of oxygen-producing, photosynthetic cyanobacteria. In this review, we summarize available Precambrian trace element proxy data, and discuss how temporal trends in the seawater concentrations of specific trace elements may be linked to the evolution of both simple and complex life. We also examine several biologically relevant and/or redox-sensitive trace elements that have yet to be fully examined in the sedimentary rock record (e.g., Cu, Cd, W) and suggest several directions for future studies.

Earth-Science Reviews↗

Is “morphodynamic equilibrium” an oxymoron?

Morphodynamic equilibrium is a widely adopted yet elusive concept in the field of geomorphology of coasts, rivers and estuaries. Based on the Exner equation, an expression of mass conservation of sediment, we distinguish three types of equilibrium defined as static and dynamic, of which two different types exist. Other expressions such as statistical and quasi-equilibrium which do not strictly satisfy the Exner conditions are also acknowledged for their practical use. The choice of a temporal scale is imperative to analyse the type of equilibrium. We discuss the difference between morphodynamic equilibrium in the “real world” (nature) and the “virtual world” (model). Modelling studies rely on simplifications of the real world and lead to understanding of process interactions. A variety of factors affect the use of virtual-world predictions in the real world (e.g., variability in environmental drivers and variability in the setting) so that the concept of morphodynamic equilibrium should be mathematically unequivocal in the virtual world and interpreted over the appropriate spatial and temporal scale in the real world. We draw examples from estuarine settings which are subject to various governing factors which broadly include hydrodynamics, sedimentology and landscape setting. Following the traditional “tide-wave-river” ternary diagram, we summarize studies to date that explore the “virtual world”, discuss the type of equilibrium reached and how it relates to the real world.

Earth-Science Reviews↗

The added value of time-variable microgravimetry to the understanding of how volcanoes work

During the past few decades, time-variable volcano gravimetry has shown great potential for imaging subsurface processes at active volcanoes (including some processes that might otherwise remain “hidden”), especially when combined with other methods (e.g., ground deformation, seismicity, and gas emissions). By supplying information on changes in the distribution of bulk mass over time, gravimetry can provide information regarding processes such as magma accumulation in void space, gas segregation at shallow depths, and mechanisms driving volcanic uplift and subsidence. Despite its potential, time-variable volcano gravimetry is an underexploited method, not widely adopted by volcano researchers or observatories. The cost of instrumentation and the difficulty in using it under harsh environmental conditions is a significant impediment to the exploitation of gravimetry at many volcanoes. In addition, retrieving useful information from gravity changes in noisy volcanic environments is a major challenge. While these difficulties are not trivial, neither are they insurmountable; indeed, creative efforts in a variety of volcanic settings highlight the value of time-variable gravimetry for understanding hazards as well as revealing fundamental insights into how volcanoes work. Building on previous work, we provide a comprehensive review of time-variable volcano gravimetry, including discussions of instrumentation, modeling and analysis techniques, and case studies that emphasize what can be learned from campaign, continuous, and hybrid gravity observations. We are hopeful that this exploration of time-variable volcano gravimetry will excite more scientists about the potential of the method, spurring further application, development, and innovation.

Earth-Science Reviews↗

Distinguishing brackish lacustrine from brackish marine deposits in the stratigraphic record: A case study from the late Miocene and early Pliocene Bouse Formation, Arizona and California, USA

Brackish marine and brackish continental environments are fundamentally different from a compositional perspective. Brackish water is often defined as having salinity lower than that of standard seawater but higher than that of freshwater, but less regard is given to the origin of the salts involved. The simple dilution of standard seawater by freshwater in a coastal or estuarine setting constitutes a brackish environment, but so do lakes where continental fresh water is impounded and becomes more saline through a variety of solute evolution pathways. The range of potential compositions of brackish lake water is diverse and includes water with “seawater-like” compositions. Isolated brackish lake environments located hundreds of kilometers inland can evolve towards sodium chloride-dominated, low alkalinityenvironments that mimic the composition of brackish seawater environments. These types of lakes can harbor a variety of continentally invasive but typically marine organisms, including but not limited to algae, foraminifers, mollusks, diatoms, and crustaceans . Distinguishing brackish marine from brackish lake environments in the geologic record can be difficult. In this paper, the enigmatic late Miocene and early Pliocene southern Bouse Formation of southern Arizona and California, USA, considered by many to represent a marine transgression along the lower Colorado River corridor, is discussed within a broad framework that incorporates hydrochemical, biogeographical, and species niche concepts. A brackish lake interpretation provides a powerful platform that can comprehensively account for the enigmatic mixed marine and continental fossil assemblage and possible tidal rhythmites that feature prominently in the southern Bouse Formation controversy. A review of the broader regional (paleo)environmental context for the southern Bouse supports a sodium chloride-dominated, low alkalinity, mildly brackish (10-5 ppt) Colorado River-fed lake depositional environment that was populated by an intriguing but predictable array of euryhaline, opportunistic, and continentally invasive marginal marine organisms.

Arizona, California↗

Igneous rocks in the Fish Creek Mountains and environs, Battle Mountain area, north-central Nevada: A microcosm of Cenozoic igneous activity in the northern Great Basin, Basin and Range Province, USA

The Great Basin of the western United States, the northern component of the Basin and Range Province, is a region of Cenozoic lithospheric extension with multiple periods and types of igneous activity. The composition and volume of Cenozoic magmas reflect a complex interaction between mantle-derived magmas and highly diverse crust, where both mantle sources and magmatic processes were modulated by tectonic environment. The Fish Creek Mountains in north-central Nevada underwent multiple igneous events ranging from ca. 40 Ma to 1 Ma that span all of the complex Cenozoic tectono-magmatic episodes of the Great Basin. The Fish Creek Mountains, therefore, is an ideal location to evaluate the different sources and processes involved in magma generation. Many plutons were emplaced in the region between about 40 and 38 Ma, several of which host base and precious metal deposits. Between 36 and 33 Ma, lava fields and calderas of the 37–19 Ma Ignimbrite Flare-up were emplaced. Both these and the preceding plutons resulted from southwestward rollback of the Farallon plate beneath North America during by far the most voluminous phase of Cenozoic magmatism. The lavas range from rare basalt and basaltic andesite to andesite, dacite, and rhyolite, have continental arc-like incompatible element patterns, and high initial 87 Sr/ 86 Sr and low εNd that require a metasomatized lithospheric mantle source combined with minor crustal component. Ignimbrites of the 34.4 Ma Cove Mine (trachydacite to rhyolite) and 34.0 Ma Caetano calderas (rhyolite to high-silica rhyolite) are abundantly porphyritic, include hydrous phases, were largely derived from partial melts of crustal rocks, but likely include 20–30% of a mantle-derived component. Igneous activity ceased in the region as the rollback-arc migrated to the southwest, but at 24.9 Ma a new caldera formed in the southern Fish Creek Mountains that was filled by ignimbrites of the Fish Creek Mountains Tuff. Intracaldera rhyolite ignimbrites range from aphyric, pumice-rich deposits at the base to progressively more quartz-feldspar phyric ignimbrites at higher levels; all flow units lack hydrous phases. No contemporaneous mafic or intermediate igneous activity accompanied caldera formation, but initial 87 Sr/ 86 Sr values in the Fish Creek Mountains tuffs are lower than in the Caetano Tuff, suggesting a greater mantle contribution to the 24.9 Ma ignimbrites. After another hiatus in igneous activity, the region was intruded and overlain by basalt to rhyolite dykes and lavas of the northern Nevada rift between 16.8 and 15.1 Ma. The primarily tholeiitic igneous suite is of the same age, chemistry, and isotopic composition as the Grande Ronde Formation of the Columbia River flood basalts, and evolved members (trachydacite and rhyolite) are crustally contaminated. The youngest northern Nevada rift lava is an alkali olivine basalt with isotopic affinity to basalts of the eastern Snake River Plain. After 10 Ma of quiescence, the region was locally covered by mafic lava flows with high-alumina olivine tholeiite compositions, represented by the 5.4 Ma Pumpernickel Valley flows. Their mid-ocean ridge-like incompatible element compositions indicate a depleted mantle source for the lavas, but radiogenic isotopic compositions indicate that the lavas of this region include a significant contribution from a mafic to ultramafic, high- 87 Sr/ 86 Sr source. The final igneous event in the Fish Creek Mountains region, the 4.0 to 1.0 Ma Buffalo Valley volcanic field, includes flows and spatter cones of transitional to alkalic basalt that are divided into two geochemical groups with identical isotopic compositions. They represent variable, low percent partial melts of the asthenosphere at different depths, yielding different rare earth element characteristics. Similar to the Lunar Crater volcanic field, the Buffalo Valley rocks represent a rare case where the lithosphere in the central Great Basin is now thin enough to allow melting of the underlying asthenosphere. Cenozoic magmatism in the northern Great Basin exhibits several transitions in magma sources and tectonic setting with time. Magmatism began as pre-extension, subduction-related, primarily lithospherically-derived magmas emplaced on/in tectonically-thickened crust. The onset of extension was partially driven by impingement of the Yellowstone plume that resulted in emplacement of rift-related volcanic and intrusive rocks in the northern Nevada rift, followed by the eruption of extension-related HAOT lavas along the northwest margin of the Great Basin. Finally, lithospheric thinning allowed for partial melting of the asthenosphere and eruption of alkaline basaltic lavas.

Nevada↗

Debris-flow monitoring and warning: Review and examples

Debris flows represent one of the most dangerous types of mass movements, because of their high velocities, large impact forces and long runout distances. This review describes the available debris-flow monitoring techniques and proposes recommendations to inform the design of future monitoring and warning/alarm systems. The selection and application of these techniques is highly dependent on site and hazard characterization, which is illustrated through detailed descriptions of nine monitoring sites: five in Europe, three in Asia and one in the USA. Most of these monitored catchments cover less than ∼10 km 2 and are topographically rugged with Melton Indices greater than 0.5. Hourly rainfall intensities between 5 and 15 mm/h are sufficient to trigger debris flows at many of the sites, and observed debris-flow volumes range from a few hundred up to almost one million cubic meters. The sensors found in these monitoring systems can be separated into two classes: a class measuring the initiation mechanisms, and another class measuring the flow dynamics. The first class principally includes rain gauges, but also contains of soil moisture and pore-water pressure sensors. The second class involves a large variety of sensors focusing on flow stage or ground vibrations and commonly includes video cameras to validate and aid in the data interpretation. Given the sporadic nature of debris flows, an essential characteristic of the monitoring systems is the differentiation between a continuous mode that samples at low frequency (“non-event mode”) and another mode that records the measurements at high frequency (“event mode”). The event detection algorithm, used to switch into the “event mode” depends on a threshold that is typically based on rainfall or ground vibration. Identifying the correct definition of these thresholds is a fundamental task not only for monitoring purposes, but also for the implementation of warning and alarm systems.

Earth-Science Reviews↗

History and evolution of seepage meters for quantifying flow between groundwater and surface water: Part 1 – Freshwater settings

More than 75 years after its introduction, the seepage meter remains the only device for directly quantifying exchange across the sediment-water interface between groundwater and surface water. This device, first presented in the literature in the 1940s, has been in a state of near-constant improvement and design change, necessitating a review of the history and evolution of the device and a description of current best-measurement practices. Part 1 of this two-part review documents the evolution of seepage meters deployed in freshwater settings, including a listing of suggestions for best-measurement and deployment practices. Part 2 covers the same scope for seepage meters deployed in marine settings. Traditional seepage meters isolate a portion of the sediment bed; seepage commonly is determined by routing the volume of flow across that isolated interface to or from a submerged measurement bag over a known time interval. The time-integrated volume is then divided by the bed area covered by the meter to obtain a seepage flux expressed in distance per time. Both the instrument and the measurement are deceptively simple, leading some early users to question the viability of the measurement. Numerous sources of error have been identified and addressed over the decades, resulting in large improvements in measurement consistency and accuracy. Duration of each measurement depends on the seepage rate and can vary from minutes to days, leading to the erroneous and yet common assumption that seepage is relatively stable over time. Designs that replace the measurement bag with a flowmeter eliminate bag-related errors and provide much finer temporal resolution. Resulting data indicate seepage is highly variable in many settings and responds to numerous sub-daily processes, including evapotranspiration, rainfall, seiches and waves. Combining direct measurements from seepage meters with other measurements, such as vertical hydraulic gradients and vertical temperature profiles, provides far better understanding of the processes that control exchange between groundwater and surface water.

Earth-Science Reviews↗

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↗

Palaeotsunamis in the Sino-Pacific region

Palaeotsunami research in the Sino-Pacific region has increased markedly following the 2011 Tōhoku-oki tsunami. Recent studies encompass a variety of potential sources and cover a full range of research activities from detailed studies at individual sites through to region-wide data collation for the purposes of database development. We synthesise palaeotsunami data from around the region drawing on key examples to highlight the progress made since 2011. We focus on a wide range of spatial and temporal scales, from region-wide to local events, from multi-millennial site records to estimates of magnitude and frequency along national coastlines. The review is based on sub-regions but in reviewing the combined records highlights common events and anomalies. In doing so we identify future research opportunities and notable findings arising from our review.

Earth-Science Reviews↗

Detrital zircon record of magmatism and sediment dispersal across the North American Cordilleran arc system (28-48°N)

As zircon U-Pb geochronology has become a leading method in sediment provenance studies and basin analysis over the past 20 years, the volume of detrital zircon data made available in published literature has enabled researchers to go beyond source-to-sink provenance studies to explore increasingly complex geologic problems. In this review, we utilize the growing body of detrital zircon data acquired from Jurassic-Paleocene forearc and foreland basin strata of the North American Cordillera to investigate the Mesozoic to earliest Cenozoic evolution of the arc and its associated basins between 28°N and 48°N. Our compilation includes 830 detrital zircon samples (101,898 individual ages) from 70 studies published between 2000 and 2020. For comparative purposes, we also compile 1307 igneous zircon U-Pb ages that characterize the magmatic history of the arc. We place primary emphasis on detrital zircon ages between 251 and 56 Ma that we infer to be uniquely derived from magmatic sources in the arc. Informed by existing knowledge of magmatic, structural, and sedimentological processes that acted on the orogen, we investigate spatial and temporal trends in these “arc-derived zircon” to establish a detrital record of arc magmatism, investigate source-to-sink relationships between the arc and adjacent basins, and discuss controls on sediment dispersal across the orogen. Our review shows that compilations of detrital zircon data from the Cordilleran forearc and foreland basin systems are excellent proxies for arc magmatism because the basins are enriched in arc-derived zircon and compilations provide space- and time-integrated records of crystallization ages. The compiled detrital zircon data support a history of continuous arc magmatism throughout Mesozoic and earliest Cenozoic time, characterized by low-volume magmatism from Triassic-Early Jurassic time (~251–174 Ma) and episodic higher-volume magmatism from Middle Jurassic-Late Cretaceous time (~174–66 Ma). These trends elucidate the initiation and timing of magmatic events at the orogen-scale and corroborate our understanding of cyclic arc behavior. Detrital zircon distributions are spatially and temporally variable both within and across basins, which we discuss relative to topographic development of the orogen and attendant responses of sediment dispersal systems. Detrital zircon distributions in the forearc signal rapid transfer of sediment from the arc to basins dominantly via fluvial processes. In contrast, detrital zircon distributions across the foreland reflect the presence of topographic barriers in the hinterland region of the arc that effectively isolated parts of the foreland. The presence of hinterland topography in turn highlights the important role of ash-fall events in delivering arc-derived zircon to the foreland, underscoring the need to consider ash-fall processes in paleodrainage reconstructions. These broad regional trends, and in general the close linkage between orogenic process and sediment dispersal, emerge from our compilation because it averages out much of the local variability observed in studies of more limited geographic or temporal extent.

North American Cordilleran arc system↗

Climate change hotspots and implications for the global subsea telecommunications network

A global network of subsea telecommunications cables underpins our daily lives, enabling >95% of global digital data transfer, $trillions/day in financial trading, and providing critical communications links, particularly to remote, low-income countries. Despite their importance, subsea cables and their landing stations are vulnerable to damage by natural hazards, including storm surges, waves, cyclones, earthquakes, floods, volcanic eruptions , submarine landslides and ice scour. However, the likelihood or recurrence interval of many of these types of events will likely change under future projected climate change scenarios, compounded by sea-level rise; potentially increasing hazard severity, creating previously unanticipated hazards, or hazards may shift to new locations during the 20–30-year operational life of cable systems. To date, no study has assessed the wide-reaching impacts of future climate change on subsea cables and landing stations on a global scale. Here, for the first time we synthesize the current evidence base, based on published peer-reviewed datasets, to fill this crucial knowledge gap, specifically to assess how and where future climate change is likely to impact subsea cables and their shore-based infrastructure. We find that ocean conditions are highly likely to change on a global basis as a result of climate change, but the feedbacks and links between climate change, natural processes and human activities are often complicated, resulting in a high degree of geographic variability. We identify climate change ‘hotspots’ (regions and locations likely to experience the greatest impacts) but find that not all areas will be affected in the same manner, nor synchronously by the same processes. We conclude that cable routes should carefully consider locally-variable drivers of hazard frequency and magnitude. Consideration should be given both to instantaneous events (e.g. landslides, tropical cyclones) as well as longer-term, sustained impacts (e.g. seabed currents that circulate even in deep water). Multiple factors can combine to increase the risk posed to subsea cables, hence a holistic approach is essential to assess the compounded effects of both natural processes and human activities in the future.

Earth Science Reviews↗

Influence of inherited structure on flexural extension in foreland basin systems: Evidence from the northern Arkoma basin and southern Ozark dome, USA

Extensional faults are key components of foreland basin systems. They form within the upper crust in response to flexure of the lithosphere and accommodate subsidence within the foredeep and forebulge depozones. Such faults are excellent proxies for orogenic system evolution and control the distribution of natural resources and hazards. However, the spatiotemporal evolution of flexural extension has not been documented previously at a regional scale, thereby limiting our understanding of underlying geodynamic controls. Here, we resolve late Paleozoic flexural extension in the northern Arkoma basin and southern Ozark dome, USA. We synthesize a large database of previous mapping, existing research, subsurface data, and geophysical data into 3D geologic and 2D kinematic models. Mesh surfaces representing several key horizons from the Carboniferous Period (ca. 335-306 Ma) were constructed. These surfaces were built from oil and gas well tops (n = ∼10,000) and surface geologic map contacts using an advanced kriging method. The mesh surfaces are offset by a complex 3D fault network, allowing detailed analysis of along-strike and down-dip variations in fault displacement. Analysis of the 3D model reveals a regular and repeated fault segmentation pattern wherein E -W striking, left- and foreland-stepping en échelon normal faults are segmented by inherited NE striking basement faults. Maximum vertical separation along the E -W normal faults is generally focused between the inherited NE-trending faults. This suggests that the inherited basement faults delocalized extensional strain during late Paleozoic normal faulting. Maximum vertical separation and fault localization may correlate to areas with high-amplitude positive magnetic anomalies interpreted as Mesoproterozoic granitic rocks. Speculative covariance of magnetic anomalies and fault displacements implies that the relatively strong basement granite concentrated stress, leading to localized faulting within the relatively thin sedimentary cover. Lastly, we show that flexural extension migrated southeast to northwest from the Chesterian-Morrowan (ca. 335-319 Ma) to the Desmoinesian (ca. 306 Ma). The migratory flexural extension may be explained by diachronous loading during Pangean assembly, or by synchronous loading but variable load compensation due to inherent factors.

Arkansas, Oklahoma↗

Palaeontological signatures of the Anthropocene are distinct from those of previous epochs

The “Great Acceleration” beginning in the mid-20th century provides the causal mechanism of the Anthropocene, which has been proposed as a new epoch of geological time beginning in 1952 CE. Here we identify key parameters and their diagnostic palaeontological signals of the Anthropocene, including the rapid breakdown of discrete biogeographical ranges for marine and terrestrial species, rapid changes to ecologies resulting from climate change and ecological degradation, the spread of exotic foodstuffs beyond their ecological range, and the accumulation of reconfigured forest materials such as medium density fibreboard (MDF) all being symptoms of the Great Acceleration. We show: 1) how Anthropocene successions in North America, South America, Africa, Oceania, Europe, and Asia can be correlated using palaeontological signatures of highly invasive species and changes to ecologies that demonstrate the growing interconnectivity of human systems; 2) how the unique depositional settings of landfills may concentrate the remains of organisms far beyond their geographical range of environmental tolerance; and 3) how a range of settings may preserve a long-lived, unique palaeontological record within post-mid-20th century deposits. Collectively these changes provide a global palaeontological signature that is distinct from all past records of deep-time biotic change, including those of the Holocene.

Earth Science Reviews↗

Dust in the Critical Zone: North American case studies

The dust cycle facilitates the exchange of particles among Earth's major systems, enabling dust to traverse ecosystems, cross geographic boundaries, and even move uphill against the natural flow of gravity. Dust in the atmosphere is composed of a complex and ever-changing mixture that reflects the evolving human footprint on the landscape. The emission, transport, and deposition of dust interacts with and connects Critical Zone processes at all spatial and temporal scales. Landscape properties, land use, and climatic factors influence the wind erosion of soil and nutrient loss, which alters the long-term ecological dynamics at erosional locations. Once in the atmosphere, dust particles influence the amount of solar radiation reaching Earth, and interact with longwave (terrestrial) radiation, with cascading effects on the climate system. Finally, the wet and dry deposition of particles influences ecosystem structure, composition, and function over both short and long-term scales. Tracking dust particles from source to sink relies on monitoring and measurement of the geochemical composition and size distribution of the particles, space-borne and ground-based remote sensing, and dust modeling. Dust is linked to human systems via land use and policies that contribute to dust emissions and the health-related consequences of particulate loads and composition. Despite the significant influence dust has in shaping coupled natural-human systems, it has not been considered a key component of the Critical Zone. Here, we demonstrate that dust particles should be included as a key component of the Critical Zone by outlining how dust interacts with and shapes Earth System processes from generation, through transport, to deposition. We synthesize current understanding from global research and identify critical data and knowledge gaps while showcasing case studies from North America.

Earth-Science Reviews↗