USGS ScienceSearch

SEARCH · USGS Science

Results for “International Journal for Disaster Risk Reduction”

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

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

1,443 records · Page 39Linked to original sources

Conceptualization and simulation of groundwater flow and groundwater availability in the Boone and Roubidoux aquifers in northeastern Oklahoma, 1980–2017

Oklahoma Groundwater Law (Oklahoma Statute § 82-1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations to determine the maximum annual yield for the State’s groundwater basins. The Boone and Roubidoux aquifers (also known as the Springfield Plateau aquifer and Ozark aquifer, respectively) are bedrock aquifers that extend from northeastern Oklahoma into Kansas, Arkansas, and Missouri. At present (2024), the Oklahoma Water Resources Board has yet to legally issue orders for the final determination of maximum annual yields for the Boone and Roubidoux aquifers. To support determination of a maxi­mum annual yield, the U.S. Geological Survey, in coopera­tion with the Oklahoma Water Resources Board, developed a hydrogeologic framework, a conceptual groundwater-flow model, and a calibrated numerical groundwater-flow model for the Boone and Roubidoux aquifers. Three types of groundwater-availability scenarios were simulated by using the calibrated numerical model. These scenarios were used to (1) estimate equal-proportionate-share groundwater withdrawal rates (groundwater withdrawal applied equally over the aquifer), (2) quantify the potential effects of projected groundwater withdrawals on groundwater storage over a 50-year period, and (3) simulate the poten­tial effects of a hypothetical 10-year drought. For the Boone aquifer, equal-proportionate-share groundwater withdrawal rates were 1.10, 0.98, and 0.96 acre-feet per acre per year for the 20-, 40-, and 50-year scenarios, respectively. For the Roubidoux aquifer, equal-proportionate-share groundwater withdrawal rates were 1.76, 1.34, and 1.25 acre-feet per acre per year for the 20-, 40-, and 50-year simulations, respectively. For the 50-year scenarios, stream seepage was minimally affected. Over the 10-year drought scenario, groundwater storage in the Boone and Roubidoux aquifers decreased by 660,451 acre-feet (6.7 percent) and 508,472 acre-feet (1.0 per­cent), respectively.

Oklahoma

Sediment thickness map of United States Atlantic and Gulf Coastal Plain Strata, and their influence on earthquake ground motions

With the recent successful accounting of basin depth ground-motion adjustments in seismic hazard analyses for select areas of the western United States, we move toward implementing similar adjustments in the Atlantic and Gulf Coastal Plains by constructing a sediment thickness model and evaluating multiple relevant site amplification models for central and eastern United States seismic hazard analyses. We digitize and combine existing sediment thickness data sets into a composite surface that delineates the base of Cretaceous sediments under the Atlantic Coastal Plain and the base of Mesozoic sediments under the Gulf Coastal Plain. Amplification models dependent on sediment thickness, site natural period, and source-to-site path length are compared with data sets of observed ground motions to evaluate the ability of the new models to improve ground motion estimates. We find that the amplification models can account for observed trends in sediment-thickness and period-dependent residuals, but some tuning is required. For example, the model of Chapman and Guo requires a reference V S 30 , the time-averaged shear-wave velocity within 30 m of the Earth’s surface, for non-Coastal Plain sites, which we estimate to be between about 1 and 2 km/s. Along with our sediment thickness model, we estimate a velocity profile for application to the Harmon et al. site-natural-period-based model in order to best match the Chapman and Guo period dependence for a broad range of sediment thicknesses. The Next Generation of Attenuation models for the eastern United States Gulf Coast path-based adjustment models can also account for seismic attenuation in the Coastal Plain sediments and reduce the standard deviation of total residuals. If enacted in the U.S. Geological Survey National Seismic Hazard Model, these amplification models will reduce predicted short-period (<1 s) and increase predicted long-period (>1 s) ground motions in the Coastal Plains appreciably.

Earthquake Spectra

The Sedimentary Geochemistry and Paleoenvironments Project Phase 2 data release: An open data resource for the study of Earth's environmental history

Geochemical data from sedimentary rocks are the primary source of information regarding Earth's surface evolution through time, including its air and water envelopes and interactions with life and deep Earth processes. The Sedimentary Geochemistry and Paleoenvironments Project (SGP) is a scientific consortium centered around open data and community-driven development of cyberinfrastructure tools and resources for sedimentary geochemistry and Earth history. Here we describe the SGP Phase 2 data release, which focused on incorporating Paleoproterozoic and Mesoproterozoic (2500–1000 million years ago) data and better accommodating carbonate data. This data release was built through the involvement of >200 researchers worldwide in academia, government, and industry, and provides the largest available public data resource for our user community in the academic fields of geochemistry, sedimentology, tectonics, paleontology, Earth history, and paleoclimate, as well as the petroleum and minerals industries. The dataset now encompasses 126,006 samples and 4,132,371 geochemical analyses. In addition to direct entry by SGP Team Members, we have ingested and incorporated datasets from the Geoscience Australia OZCHEM database, the Alberta Geological Survey, and the Deep-Time Marine Sedimentary Element Database (DM-SED) compilation. This paper details sampling in the Phase 2 dataset with respect to age, geography, lithology, and other geological characteristics, documents access via our search website and API, discusses possible issues and/or biases in the dataset that could impact analyses, describes plans for governance and stewardship of data from Indigenous lands, and serves as the citable reference paper for the data release.

Chemical Geology

Integrating Sr isotopes, microchemistry, and genetics to reconstruct Salmonidae species and life history

Recent approaches to fisheries research emphasize the importance of the coproduction of knowledge in building resilient and culturally mindful fisheries management frameworks. Despite widespread recognition of the need for Indigenous knowledge and historical reference points as baseline data, archaeological data are rarely included in conservation biology research designs. Here we propose a novel multiproxy method to learn from former fisheries stewards by generating archaeological data on past salmonid population parameters. We used a newly developed, high throughput qPCR (HT-qPCR) chip, originally designed for environmental DNA (eDNA), for species identification of archaeological salmonid vertebrae. We combine this with the laser ablation split-stream (LASS) approach to identify ocean-migration versus freshwater residency. We test this multidisciplinary approach using both contemporary and archaeological salmonid samples and new radiocarbon dates from the Tronsdal Site on the Skagit River, Washington State, USA. This is a useful approach for extracting information about Salmonidae species and life history diversity from archaeological remains to reconstruct historic baselines for several population parameters in anadromous species with long periods of freshwater residency. The approach outlined in this paper may be particularly useful for research investigating past fisheries dynamics, offering hundreds to thousands of years of temporal depth for modern fisheries management, harvest policies, restoration ecology, and conservation biology.

Idaho, Oregon, Washington

Timescales of cumulate mobilization and mixing for the 1868 A.D. eruption of Mauna Loa, Island of Hawai‘i

The deadly 1868 A.D. eruption of Mauna Loa’s lower Southwest Rift Zone (Island of Hawai‘i) included a M7.9 earthquake and associated tsunami and landslides, demonstrating the severe hazards posed by Earth’s largest active subaerial volcano. To better understand the relationship between intense seismic activity, dike emplacement, magma storage, transport histories, and mobilization of olivine cumulates at Mauna Loa, we examine compositional zoning of olivine in the 1868 lava flows. Samples range from basalt (< 10% olivine) to picrite (30–40% olivine). The olivine cargo is heterogeneous (Fo 78.2–89.2 ; forsterite = [Mg/(Mg + Fe) × 100]) but dominated by ~ Fo 89 cores that lie above the Fe-Mg equilibrium field of host glasses. Crystal rims < Fo 80 are due to post-eruptive modification in slow cooling lava flows. Minor element compositions fall within the range of other Mauna Loa olivine erupted in the past 200 years. Olivine crystals exhibit both normal and complex Fo zoning patterns that yield timescales of diffusive re-equilibration that range from 3 to 258 days, with 72% of crystals recording 71 days or less. These timescales correspond to magmatic priming of the summit reservoir system ~ 2 months prior to the eruption and the M7.9 earthquake likely facilitated the transport of the crystal-rich summit-derived magmas downrift shortly prior to eruption. If the recently proposed faster Fe-Mg diffusion coefficient is used, timescales instead range from < 1 day to 25 days, with most recording 1 week or less. In this scenario, most of the olivine zoning would have to have been generated after the M7.9 earthquake perturbed the system.

Hawaii

Permafrost–wildfire interactions: active layer thickness estimates for paired burned and unburned sites in northern high latitudes

As the northern high-latitude permafrost zone experiences accelerated warming, permafrost has become vulnerable to widespread thaw. Simultaneously, wildfire activity across northern boreal forest and Arctic/subarctic tundra regions impacts permafrost stability through the combustion of insulating organic matter, vegetation, and post-fire changes in albedo. Efforts to synthesis the impacts of wildfire on permafrost are limited and are typically reliant on antecedent pre-fire conditions. To address this, we created the FireALT dataset by soliciting data contributions that included thaw depth measurements, site conditions, and fire event details with paired measurements at environmentally comparable burned and unburned sites. The solicitation resulted in 52 466 thaw depth measurements from 18 contributors across North America and Russia. Because thaw depths were taken at various times throughout the thawing season, we also estimated end-of-season active layer thickness (ALT) for each measurement using a modified version of the Stefan equation. Here, we describe our methods for collecting and quality-checking the data, estimating ALT, the data structure, strengths and limitations, and future research opportunities. The final dataset includes 48 669 ALT estimates with 32 attributes across 9446 plots and 157 burned–unburned pairs spanning Canada, Russia, and the United States. The data span fire events from 1900 to 2022 with measurements collected from 2001 to 2023. The time since fire ranges from 0 to 114 years. The FireALT dataset addresses a key challenge: the ability to assess impacts of wildfire on ALT when measurements are taken at various times throughout the thaw season depending on the time of field campaigns (typically June through August) by estimating ALT at the end-of-season maximum. This dataset can be used to address understudied research areas, particularly algorithm development, calibration, and validation for evolving process-based models as well as extrapolating across space and time, which could elucidate permafrost–wildfire interactions under accelerated warming across the high-northern-latitude permafrost zone. The FireALT dataset is available through the Arctic Data Center ( https://doi.org/10.18739/A2RN3092P , Talucci et al., 2024).

Arctic

From critical minerals to food security, the benefits of data collaboration

The volume of data in the public geoscience sphere is rapidly and continually expanding. At Geoscience Australia (GA) we saw an over 500% increase in data points within our relational databases between 2018 and 2024, over the life of the Exploring for the Future (EFTF) program. With the Resourcing Australia’s Prosperity initiative, a continued increase in data quantity will be seen for the next 10 to 35 years. At the same time, a broadening audience for geoscience data is increasing the desire to enhance the diversity of delivery streams. This ranges from data-dense highly technical outputs for geoscience specialists to curated interpretive products for people who are non-geoscientists. Development of these curated outputs has contributed to our awareness of the need for data to be collected and compiled in a way that ensures its reuse, with a focus on quality metadata and data provenance.

Conference Paper

Non-destructive sampling of larval amphibians shows tail tissues reflect whole-body methylmercury concentrations and trophic-related differences in bioaccumulation

Mercury (Hg), and specifically methylmercury (MeHg), is a contaminant of global concern to humans and wildlife, but there is still limited understanding of its effects on many taxa. Amphibians are closely associated with water, where Hg is converted to MeHg, and high concentrations of MeHg can reduce survival of amphibians in the wild. However, limited non-lethal proxies exist for estimating MeHg bioaccumulation in amphibian larvae, which could be used to reduce destructive sampling and allow sampling of imperiled species. We evaluated whether tail clips were indicative of whole-body MeHg concentrations of larval amphibians, whether dragonflies were an effective bioindicator of Hg concentrations of larval amphibians, and compared MeHg bioaccumulation in predatory caudate (salamanders, newts) and non-predatory anuran (frogs, toads) larvae. Tail-clip and whole-body MeHg concentrations were strongly correlated, especially for caudates. There was high variability among sites and species, but caudates had higher MeHg than anurans, and MeHg concentrations of caudates and anurans were strongly correlated. Dragonflies were a weak indicator of MeHg for anurans but strong indicators for larval caudates, likely because caudates and dragonflies are carnivorous and occupy similar trophic positions. Our study revealed that sampling tails can be an effective non-destructive index of whole-body MeHg for larval amphibians and demonstrates differences in MeHg bioaccumulation related to trophic position in larval amphibians and dragonflies. All post-embryo life stages for a broad suite of amphibian species can now be sampled non-lethally for MeHg, allowing for a broader examination of population-level effects and focus on imperiled species that cannot be sampled destructively.

conterminous United States

Quality of groundwater used for domestic supply in the eastern Sacramento Valley and adjacent foothills, California

Summary More than 2 million Californians rely on groundwater from privately owned domestic wells for drinking-water supply. This report summarizes a water-quality survey of domestic and small-system drinking-water supply wells in the eastern Sacramento Valley and adjacent foothills where more than 25,000 residents are estimated to use privately owned domestic wells. Study results show that inorganic and organic constituents in groundwater were present above regulatory (maximum contaminant level, MCL) benchmarks for public drinking-water quality in 8 and 3 percent, respectively, of the aquifer area used for domestic drinking-water supply (herein, “domestic groundwater resources”; fig. 1). The only inorganic constituent detected above regulatory benchmarks was arsenic. The only organic constituent exceeding regulatory benchmarks was the fumigant 1,2,3-trichloropropane (1,2,3-TCP). Three additional organic constituents—the disinfection by-product chloroform, the gasoline oxygenate methyl tert -butyl ether (MTBE), and the solvent tetrachloroethene (PCE)—were detected at low concentrations below one-tenth of regulatory benchmarks in 34, 10, and 10 percent of domestic groundwater resources, respectively. Total dissolved solids (TDS), iron, and manganese exceeded non-regulatory aesthetic guidelines for drinking water in 5, 10, and 26 percent of domestic groundwater resources, respectively. Per- and polyfluoroalkyl substances (PFASs) were detected in 29 percent of domestic groundwater resources,with 5 percent exceeding the recently enacted (April 2024) U.S. Environmental Protection Agency MCLs. Total coliform and enterococci bacteria were detected in 13 and 8 percent of domestic groundwater resources, respectively. Redox sensitive constituents in this study included arsenic, manganese, nitrate, and iron. In the lower elevation portions of the eastern Sacramento Valley study area, reducing conditions in groundwater aquifers promote elevated arsenic, iron, and manganese, and conversely lower concentrations of nitrate. The presence of the volatile organic compound (VOC) 1,2,3-TCP was related to its past history in select agricultural land uses (on orchards or vineyards) in the Sacramento Valley; however, unlike in the San Joaquin Valley where orchards and vineyards are more common, its detection frequency was low (only detected in one well in this study). Chloroform was frequently detected in this study at low levels. Chloroform is a disinfection byproduct commonly found in domestic wells treated by shock chlorination. The solvent PCE is among the most frequently detected VOCs in groundwater, which is primarily related to its long history of use and its persistence in groundwater in oxic conditions. The gasoline oxygenate MTBE was a contaminant introduced to groundwater through atmospheric exchange when it was used as a fuel additive to decrease smog inducing emissions from vehicles. Its occurrence in groundwater at low levels is expected and makes it a potentially useful tracer of relatively recent recharge water being withdrawn from wells. The PFASs are anthropogenic chemicals with hundreds of uses, and they have been incorporated into many different products, processes, and applications worldwide. Like MTBE, the occurrence of PFASs in groundwater may be in part due to atmospheric exchange, but there are several other pathways that contribute PFASs to the environment.

California

Ground deformation and gravity for volcano monitoring

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

Scientific Investigations Report

Evaluation of nutrient, alkalinity, and acid-neutralizing capacity stabilities in water samples analyzed by the U.S. Geological Survey National Water Quality Laboratory, 2023–24

The U.S. Geological Survey evaluated the stability of water-sample chemical analysis of nutrient, alkalinity, and acid-neutralizing capacity constituents with respect to the duration between sample collection and laboratory analysis, also known as the sample holding time. A study began in the spring of 2023 to evaluate the sample stability, between 2 and 180 days after sample collection, of the chemical properties and chemical constituents of alkalinity as calcium carbonate, filtered; acid-neutralizing capacity as calcium carbonate, unfiltered; total ammonia as nitrogen, filtered; total ammonia plus organic nitrogen as nitrogen, filtered and unfiltered; nitrite as nitrogen, filtered; nitrate plus nitrite as nitrogen, filtered; total nitrogen, filtered and unfiltered; orthophosphate as phosphorous, filtered; and total phosphorus as phosphorus (filtered and unfiltered) in water. Both surface water and groundwater matrices were represented. Sample instability varied by observed property and matrix; therefore, providing general guidance for sample holding time is not possible based on matrices alone. No correlations between field measurements of sample characteristics and sample instability were observed. Although observations for some properties indicate sample stability that exceeds the recognized U.S. Geological Survey National Water Quality Laboratory method holding times, this is not necessarily the case for matrices and seasonal characteristics that were not investigated. Based on the limited number of six sample sources used in this study, some patterns emerge for the 12 observed properties studied. Five observed properties generally indicate stability for as many as 180 days after sampling (total nitrogen as nitrogen, both filtered and unfiltered; orthophosphate as phosphorus, filtered; and phosphorus as phosphorus, both filtered and unfiltered). Other observed properties indicate stability for as many as 180 days for some matrices, but not for others. Finally, some observed properties indicate instability well before 180 days.

Scientific Investigations Report

Use of a numerical groundwater-flow model and projected climate scenarios to simulate the effects of future climate conditions on base flow for reach 1 of the Washita River alluvial aquifer and Foss Reservoir storage, western Oklahoma

To better understand the relation between climate variability and future groundwater resources in reach 1 of the Washita River alluvial aquifer and Foss Reservoir in western Oklahoma, the U.S. Geological Survey, in cooperation with the Bureau of Reclamation, used a previously published numerical groundwater-flow model and climate-model data to investigate changes in base flow and reservoir storage by evaluating three scenarios. The three projected climate scenarios were (1) a central-tendency scenario, (2) a warmer/drier scenario, and (3) a less-warm/wetter scenario. To estimate future base flow and groundwater availability in western Oklahoma, specifically in reach 1 of the Washita River alluvial aquifer, downscaled climate-model data from 231 Coupled Model Intercomparison Project phase 5 (CMIP5) projections coupled with a previously published numerical groundwater-flow model were used to compare the effects of different climate scenarios on the aquifer. Changes in base flow and groundwater-level elevations during a 30-year baseline scenario (1985–2014) and the three 30-year projected climate scenarios (2050–79) under central-tendency, warmer/drier, and less-warm/wetter climatic conditions were assessed by using the calibrated model. In the simulations, the amount of base flow and reservoir storage declined in the central-tendency and warmer/drier scenarios compared to the amount of base flow and reservoir storage under historical climatic conditions (baseline scenario). Mean annual change in reservoir storage decreased from the baseline scenario the most in the warmer/drier scenario, followed by the central-tendency scenario, but increased in the less-warm/wetter scenario compared to the baseline scenario. At the end of the simulation period (2079), the largest magnitude differences in groundwater-level elevations in all three projected climate scenarios relative to the baseline scenario occurred upstream from Foss Reservoir. Results from incorporating downscaled climate projections into localized numerical groundwater-flow models can highlight potential future changes in and implications for groundwater resources and availability.

Oklahoma

Shallow lake, strong shake: Record of seismically triggered lacustrine sedimentation from the 1959 M7.3 Hebgen Lake earthquake within Henrys Lake, Idaho

We investigate a shallow lake basin for evidence of a large historic intraplate earthquake in western North America. Henrys Lake, Idaho is an atypical candidate for a lacustrine paleoseismic study given its shallow depth (~7 m) and low relief (≤2° slopes ). Here, we test the earthquake-recording capacity of this basin type by showing sedimentological evidence of the 1959 M7.3 Hebgen Lake earthquake within sediment cores, using anthropogenically produced 137 Cs activity to constrain timing. In addition to expanding the morphologic range of basins targeted for lacustrine paleoseismic studies, this work has implications for sediment response in dam-enhanced basins. Lack of sedimentological evidence for other earthquakes coupled with radiocarbon chronology reveals that the 1959 event is the only clearly recorded earthquake within Henrys Lake since the mid-Holocene. Henrys Lake offers a proxy for paleo-earthquake signatures within similar lacustrine environments and underscores the importance of further paleoseismic studies in the region.

Idaho

Estimated hydrogeologic, spatial, and temporal distribution of self-supplied domestic groundwater withdrawals for aquifers of the Virginia Coastal Plain

Water use from private-domestic wells accounts for nearly 40 percent of total groundwater withdrawals in the Virginia Coastal Plain Physiographic Province (henceforth called the Virginia Coastal Plain). However, because self-supplied domestic water use generally falls below the Virginia Department of Environmental Quality (VDEQ) reporting and management threshold of 300,000 gallons per month, quantifying these withdrawals is challenging. This report builds upon the foundation of previous U.S. Geological Survey investigations by providing revised techniques to improve estimates of the aquifer source, spatial distribution, and monthly magnitude of these groundwater withdrawals. The aquifer sources of private-domestic wells in the Virginia Coastal Plain were estimated by cross-referencing 8,264 well records from the VDEQ and the Virginia Department of Health to a digital model of the Virginia Coastal Plain hydrogeologic framework. This analysis highlights the regional importance of the Yorktown-Eastover, Potomac, and surficial aquifers. Collectively, these three aquifers account for 80 percent of self-supplied domestic groundwater withdrawals. The population using self-supplied domestic water was estimated using census blocks, well-use ratios, building footprints, and land-use and land-cover data to produce a high-resolution, disaggregated, raster-based dataset. This approach improves upon previous models at the census-block or road-network scale by reducing the low-density spread of the self-supplied domestic population across undeveloped areas and concentrating the population and its corresponding water use in the areas where it is most likely to occur. Results show that an estimated 475,332 people comprise the 2020 self-supplied domestic population of the Virginia Coastal Plain, an increase of 5.7 percent since 2010, and the greatest concentrations of self-supplied domestic population surround large cities. Estimates could be further refined with the addition of current and complete spatial data on public water-system service areas. The quantity of water used by the self-supplied domestic population was estimated by modifying published state per-capita water-use coefficients with the corresponding monthly variability assessed from Virginia Coastal Plain public water-system withdrawal data. This analysis estimates an average increase of 12 percent from June through August and an average decrease of 8 percent from December through March from the baseline annual average of 80 gallons per day per capita, which generally matches similar studies in the eastern United States. The application of these revised methodologies for the estimation of private-domestic wells and the self-supplied domestic population improves understanding of domestic groundwater use in the Virginia Coastal Plain across hydrogeologic, spatial, and temporal scales. These revisions help better inform water-resource managers and decision makers and support higher resolution groundwater modeling. Furthermore, these methods are transferrable to other areas where self-supplied domestic water withdrawals are important to the overall water budget.

Virginia

Regionalization of low-flow characteristics of Tennessee streams

Procedures for estimating 3-day 2-year, 3-day 10-year, 3-day 20-year, and 7-day 10-year low flows at ungaged stream sites in Tennessee are based on surface geology and drainage area size. One set of equations applies to west Tennessee streams, and another set applies to central and east Tennessee streams. The equations do not apply to streams where flow is significantly altered by activities of man. Standard errors of estimate of equations for west Tennessee are 24 to 32% and for central and east Tennessee 31 to 35%. Streamflow recession indexes, in days/log cycle, are used to account for effects of geology of the drainage basin on low flow of streams. The indexes in Tennessee range from 32 days/log cycle for clay and shale to 350 days/log cycle for gravel and sand, indicating different aquifer characteristics of the geologic units that sustain streamflows during periods of no surface runoff. Streamflow recession rate depends primarily on transmissivity and storage characteristics of the aquifers, and the average distance from stream channels to basin divides. Geology and drainage basin size are the most significant variables affecting low flow in Tennessee streams according to regression analyses. (Author 's abstract)

Tennessee

Geophysical constraints on continental rejuvenation in central China: Implications for outward growth of the Tibetan Plateau

Continental rejuvenation results from the tectonic reactivation of crustal structures and lithospheric reworking by mantle flow. Geochemical observations and field mapping have traditionally provided the primary evidence for the secular evolution of crustal composition and tectonic processes during continental rejuvenation. Nonetheless, the impact of continental rejuvenation on the observed present-day strain rate and orogenic-scale lithospheric structure has not been well constrained. The pre-existing E-W–trending Central China Orogenic Belt has been overprinted by the N-S–trending Central Longitudinal Seismic Belt and constitutes the intracontinental West Qinling Syntaxis in central China, where the tectonic setting changes eastward from contraction to extension. Combining updated global positioning system data and high-resolution crustal seismic tomography, we reveal a modern continental rejuvenation process within the West Qinling Syntaxis in central China. The northward extrusion of the Tibetan Plateau's weak lithospheric layer (middle-lower crust and lithospheric mantle) of southwestern China relative to the rigid Sichuan Basin/Ordos Block of the eastern West Qinling Syntaxis results in regional dextral shearing that shapes the Central Longitudinal Seismic Belt and defines the eastern Tibetan Plateau margin. The pre-existing E-W–trending Central China Orogenic Belt has been preserved above the brittle-ductile transition zone, and the northward movement of the deep lithospheric layer drives the deformation of the upper crust in the West Qinling Syntaxis. Our results, along with previous studies, suggest the presence of an intracontinental lithospheric interchange structure in central China. The continental rejuvenation of the West Qinling Syntaxis results from a combination of fault reactivation in the upper crust (Stage I, Eocene–Oligocene) and reworking of the deep lithosphere (Stage II, middle–late Miocene) related to the plateau-wide shift in stress accommodation ultimately driven by the redistribution of mass outward from the central Tibetan Plateau. At present, the transition zone between the high- and low-velocity anomalies along the Central Longitudinal Seismic Belt not only shapes the landscape boundary but controls the size and recurrence interval of earthquakes within the West Qinling Syntaxis in central China.

GSA Bulletin

A deep dive into subduction zones and the Japan Trench

Deep beneath the ocean, giant slabs of rock slide beneath one another in areas called subduction zones. These powerful movements shape the planet by building mountains, creating fiery volcanoes, and even causing huge earthquakes. Scientists study subduction zones to understand why these natural events happen and where and how they might occur so we can be better prepared. One important subduction zone is the Japan Trench, where a big earthquake happened in 2011. To learn more about subduction zones, scientists go on ocean drilling expeditions. They sail on research ships, drill into the seafloor, and collect rock samples. These samples help scientists understand what happened in the past and how subduction zones work. In turn, this information helps communities better understand and prepare for future earthquakes.

Frontiers for Young Minds

Combining multisite tsunami and deformation modeling to constrain slip distributions for the 1700 C.E. Cascadia earthquake

A major earthquake ruptured the Cascadia subduction zone (CSZ) on 26 January 1700. Key paleoseismic evidence associated with this event include tsunami deposits, stratigraphic evidence of coastal coseismic subsidence, written Japanese records of a tsunami unaccompanied by earthquake shaking, and margin‐wide turbidites found offshore and in lacustrine environments. Despite this wealth of independent clues, important details about this event remain unresolved. Dating uncertainties do not conclusively establish whether the proxies are from one earthquake or a sequence of them, and we have limited knowledge of the likely slip distributions of the event or events. Here, we use a catalog of 37,500 candidate synthetic ruptures between M w 7.8 and 9.2 and simulate their resulting coseismic deformation and tsunami inundation. Each model is then compared against estimated Japan tsunami arrivals, regional coastal subsidence records, and local paleotsunami deposits mapped at six different coastal marshes and one coastal lake along the CSZ. We find that seven full‐margin ruptures with a median magnitude of M w 9.1 satisfy all three constraints. We favor one M w 9.11 model that best matches all site paleoseismic observations and suggests that the Cascadia megathrust slipped up to ∼30 m and must have shallow geodetic coupling. We also find that some sequences composed of three or four ruptures can still satisfy the observations, yet no sequences of two ruptures can. Sequences are differentiated into three groups based on whether they contain a mainshock rupture located in the south (>44° N) or further north. All sequences contain unruptured portions of the megathrust and most contain mainshocks with peak slip above 40 m. The fit of the geologic evidence from sequences is poor in comparison to single‐event models. Therefore, sequences are generally less favored compared to full‐margin events.

British Columbia, California, Oregon, Washington