USGS ScienceSearch

SEARCH · USGS Science

Results for “Nature Reviews Earth & Environment”

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.

635 records · Page 19Linked to original sources

Arctic speleothems reveal nearly permafrost-free Northern Hemisphere in the Late Miocene

Arctic warming is happening at nearly four times the global average rate. Long-term trends of permafrost dynamics cannot be estimated directly from monitoring of present-day thaw processes, requiring paleoclimate-proxy information. Here we use cave carbonates (speleothems) from a northern Siberian cave to determine when the Northern Hemisphere was mostly permafrost-free. At present, thick continuous permafrost in this region prevents speleothem growth. In a series of partially eroded caves, speleothems grew during the late Tortonian stage (8.68 ± 0.09 Ma), a time when the geographic position of this site was already similar to today. Paleotemperatures reconstructed from speleothems show that mean annual air temperatures (MAAT) in the region were + 6.6°C to + 11.1°C, when contemporary global MAAT were ~ 4.5 °C higher than modern. Our findings provide direct evidence that warming to Tortonian-like temperatures would leave most of the Northern Hemisphere permafrost-free. This may release up to ~ 130 petagrams of carbon, enhancing further warming.

Siberia

Modeling interconnected minerals markets with multicommodity supply curves: Examining the copper-cobalt-nickel system

Demand for many of the metals used in the energy transition is expected to grow rapidly. Many of these are by-products, often considered critical because their production responds weakly to prices and is instead tied to the economics of the host mineral. We present a model of prices and production for jointly produced commodities that accounts for interconnectivity between host and by-product markets at the mine level. We demonstrate this method using the copper–cobalt–nickel system, in which approximately 99% of cobalt is a by-product of copper or nickel mining. Our results show that the model more accurately captures the economic benefits of diversified mine outputs than previous approaches. Furthermore, changes in demand drivers for any two commodities produce non-linear effects on production and price. We challenge the prior best-practice assumption that cobalt cannot impact the copper or nickel markets. Recognizing the importance of both copper and cobalt for future electrification, we emphasize that incentivizing the copper industry to reduce cobalt supply risks could inadvertently undermine copper supply.

Nature Communications

Tectonically driven integration of the 4.8 Ma Colorado River USA tracked with detrital sanidine and fish genetics

The development of the continental-scale Colorado River system, western USA, from 8 to 4.8 Ma, is revealed using 60-40 Ma detrital sanidine tracer grains and fish phylogeny. Here we show that precursor paleoriver segments became integrated north to south as traced by 60-40 Ma sand grains that were derived from the north and sequentially appeared in the 25-8 Ma Browns Park Formation of Utah, 7-6 Ma upper Bidahochi Formation of Arizona, and 4.8 Ma Bouse Formation of the lower Colorado River and proto Gulf of California. This timing is mimicked by molecular clock estimates of divergence times among fish lineages. River integration was a response to headwater uplifts in the Yellowstone hotspot track and Rocky Mountains. 40 Ar/ 39 Ar ages refine the timing for mantle-drips that caused subsidence, then uplift, of depositional basins that influenced the integration pathway and tempo. The ~ 3-million-year timescale suggests that multiscale mantle-driven uplift, rather than lake spillover, was the primary driver for integration of the proto-Colorado River through Grand Canyon.

Colorado River

Reassessing water availability in the Nile River Basin using satellite-based irrigation water use accounting

Existing estimates of irrigation areas and water use in the Nile River Basin vary widely and are based on inconsistent data and methodologies. This study leverages advances in remote sensing to provide an up-to-date and consistent basin-wide reassessment of water use in the Nile River Basin. Land cover and actual evapotranspiration data from 2013 to 2022 were used to quantify irrigated areas, irrigation water use, and the average naturalized water yield. Across the basin, 7.1 million hectares of irrigated land consume an average of 72.5 ± 2.6 billion cubic meters (BCM) of blue water annually. Egypt and Sudan together account for 92% of the irrigated area and 97% of the irrigation water use. The basin’s average naturalized yield, evaluated at Dongola station, is determined to be between 108 and 114.4 ± 3.2 BCM. Although conservative, this updated analysis indicates the basin’s average naturalized yield is 14%–20% higher than the previous naturalized yield estimate (95 BCM) and 29%–36% higher than the commonly cited observed flow of 84 BCM, which is frequently treated as the total basin-wide available water. Importantly, the study underscores that the 84 BCM represents flow at the Aswan gauging station after upstream consumptive uses, not the total water available in the basin under naturalized conditions. This distinction is critical for accurate water accounting, planning, and governance in this water-scarce basin. This study demonstrates the value of open-source remote sensing resources in data-scarce regions while emphasizing the need for region specific validation and bias correction for improved accuracy.

Nile River Basin

Earthquake-triggered ground-failure inventory associated with the M7.1 2018 Southcentral Alaska earthquake

The 30 November 2018, magnitude (Mw) 7.1 earthquake in Southcentral Alaska triggered substantial landslides, liquefaction, and ground cracking throughout the region, resulting in widespread geotechnical damage to buildings and infrastructure. Despite a challenging reconnaissance and remote-sensing environment, we constructed a detailed digital inventory of ground failure associated with the event from several sources. Sources included information derived from remotely sensed data, and data compiled from literature, social media postings, and earthquake damage information compiled by local, state, and federal agencies. Each instance of ground failure within the inventory contains information on the location and type of observed ground failure, and the methods and data used to document the occurrence. Where high-quality data, such as LIDAR or satellite imagery, were available and showed the ground-failure instance clearly, the extent is mapped as a polygon or polyline. All other locations are mapped as points. There are a total of 886 ground-failure instances documented within the inventory (400 landslides, 286 liquefaction features, and 200 features unattributed to specific processes). A semi-quantitative confidence scheme is used to describe mapping certainty associated with each ground-failure feature. This inventory represents a relatively moderate ground-failure-triggering event that occurred in a subarctic environment. This data paper describes the content within the inventory, the inventory data collection procedures, and limitations of the data. Events of this type are not often documented in detail; thus, adding the inventory data to the US Geological Survey Open Repository of Earthquake-Triggered Ground-Failure Inventories further diversifies the datasets available to the scientific community to be used to better understand and model earthquake-triggered ground failure.

Alaska

Rapid earthquake magnitude classification via P-wave strains from borehole strainmeters and Distributed Acoustic Sensing

Distributed Acoustic Sensing (DAS) offers a promising approach for earthquake early warning (EEW) in settings where seismic networks are costly to maintain. By repurposing fiber-optic cables as dense strainmeter arrays, DAS enables real-time earthquake detection wherever those fibers are accessible. However, poor azimuthal coverage and challenges in estimating magnitude from strain measurements remain key hurdles in applying for earthquake monitoring. Here, we develop a machine learning method to distinguish large (M≥5.4) earthquakes from smaller ones within the first 4 seconds of a strain waveform after a P-wave arrival without determining location. Using ensemble decision tree models trained on borehole strainmeter data (3.5≤M≤7.1) and tested on onshore DAS waveforms (including the 2024 M7 Offshore Cape Mendocino earthquake), we find that low-frequency (0.2–0.5 Hz) continuous wavelet transform coefficients are the strongest predictors of magnitude, in addition to strain amplitude. Both DAS and borehole strainmeters effectively capture long-period strain signals, making these findings valuable for EEW systems. Our method shows high precision compared to the real-time EEW system, ShakeAlert®, supporting the position that DAS is a viable technology for earthquake monitoring and magnitude classification.

California

Woods Hole Coastal and Marine Science Center—2023 annual report

The 2023 annual report of the U.S. Geological Survey Woods Hole Coastal and Marine Science Center highlights accomplishments of 2023, includes a list of 2023 publications, and summarizes the work of the center, as well as the work of each of its science groups. This product allows readers to gain a general understanding of the focus areas of the center’s scientific research and learn more about specific projects and progress made throughout 2023, all while enjoying photographs taken in various environments and laboratories, and applicable maps and figures.

Circular

Woods Hole Coastal and Marine Science Center—2024 annual report

The 2024 annual report of the U.S. Geological Survey Woods Hole Coastal and Marine Science Center highlights accomplishments of 2024, includes a list of 2024 publications, and summarizes the work of the center, as well as the work of each of its science groups. This product allows readers to gain a general understanding of the focus areas of the center’s scientific research and learn more about specific projects and progress made throughout 2024, all while enjoying photographs taken in various environments and laboratories, and applicable maps and figures.

Circular

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

Projections of multiple climate-related coastal hazards for the US Southeast Atlantic

Faced with accelerating sea level rise and changing ocean storm conditions, coastal communities require comprehensive assessments of climate-driven hazard impacts to inform adaptation measures. Previous studies have focused on flooding but rarely on other climate-related coastal hazards, such as subsidence, beach erosion and groundwater. Here, we project societal exposure to multiple hazards along the Southeast Atlantic coast of the United States. Assuming 1 m of sea level rise, more than 70% of the coastal residents and US \ $ 1 trillion in property are in areas projected to experience shallow and emerging groundwater, 15 times higher than daily flooding. Storms increase flooding exposure by an order of magnitude over daily flooding, which could impact up to ~50% of all coastal residents and US \ $ 770 billion in property value. The loss of up to ~80% of present-day beaches and high subsidence rates that currently affect over 1 million residents will exacerbate flooding and groundwater hazard risks.

Florida, Georgia, North Carolina, South Carolina,

A 700-year rupture sequence of great eastern Aleutian earthquakes from tsunami modeling of stratigraphic records

Great Aleutian underthrusting earthquakes produced destructive tsunamis impacting Hawaiʻi in 1946 and 1957. Prior modeling of the 1957 tsunami deposit and runup records on eastern Aleutian and Hawaiian Islands jointly with tide-gauge observations across the Pacific Ocean constrained a rupture model with shallow slip up to 26 m along 600 km of the plate boundary. Here we implement this modeling approach to older deposits and show alternating deep and shallow megathrust slip up to 26, 32, and 22 m for great earthquakes along the same segment in the 18 th , 15 th , and 14 th centuries. All three modeled prehistoric Aleutian earthquakes produce tsunami inundation in Hawaiʻi with the most severe, 14 th century event having impacts exceeding the 1957 event. The along-dip variability of these four ruptures spanning seven centuries provides insights on earthquake cycles for engineering design and hazard assessment. The 15 th century and 1957 rupture models provide evidence for recurrence of tsunami earthquakes, which can produce disproportionately large tsunamis for a given moment magnitude due to reduced rigidity in the shallow megathrust. The 14 th and 18 th century events likely ruptured deeper regions that did not slip in 1957, suggesting potential for corresponding deeper failure in the next great eastern Aleutian earthquake.

Alaska, Hawaii

Biogeochemical and community ecology responses to the wetting of non-perennial streams

Transitions between dry and wet hydrologic states are the defining characteristic of non-perennial rivers and streams, which constitute the majority of the global river network. Although past work has focused on stream drying characteristics, there has been less focus on how hydrology, ecology and biogeochemistry respond and interact during stream wetting. Wetting mechanisms are highly variable and can range from dramatic floods and debris flows to gradual saturation by upwelling groundwater. This variation in wetting affects ecological and biogeochemical functions, including nutrient processing, sediment transport and the assembly of biotic communities. Here we synthesize evidence describing the hydrological mechanisms underpinning different types of wetting regimes, the associated biogeochemical and organismal responses, and the potential scientific and management implications for downstream ecosystems. This combined multidisciplinary understanding of wetting dynamics in non-perennial streams will be key to predicting and managing for the effects of climate change on non-perennial ecosystems.

Nature Water

Estimating earthquake source depth using teleseismic broadband waveform modeling at the USGS National Earthquake Information Center

The U.S. Geologic Survey National Earthquake Information Center (NEIC) monitors global seismicity, producing a catalog of earthquake source parameters in near-real-time to provide information that can help mitigate the societal impact of earthquakes. The NEIC commonly relies on teleseismic observations to constrain earthquake source parameters (e.g., location, depth, magnitude, and mechanism) due to a lack of local and regional observations. For these ‘teleseismic’ events, depth phase (i.e., pP , sP ) arrival time observations provide the best estimate on source depth. However, depth phases are often difficult to accurately identify and/or pick. Therefore, NEIC relies on waveform modeling, such as those determined from W-phase ( M ww ), body wave ( M wb ), and regional ( M wr ) moment tensor estimations, to provide constraints on source depth. While depth estimates from these approaches are informative, higher frequency observations provide more precise estimates because depth phases are more prominently observed at higher frequencies. Here, we present NEIC’s relatively high-frequency (~0.04 to 1 Hz) teleseismic waveform modeling approach, termed Synthetic Depth Phase Modeling (SynDepth), for determining source depth. SynDepth was developed to provide NEIC with a tool that enables rapid, accurate, and quantifiable estimates of earthquake source depth in cases where locator depths are not reliable. This relatively simple and fast procedure searches over 1 km-incremented source depths and an expanding triangular source-time function to find the best-fitting solution. We compare automatic SynDepth solutions for a dataset of 1,216 earthquakes (M5.5-M7.6) between 2017 and 2021 to NEIC-derived depth estimates from other methods. Our approach provides a robust depth estimate for earthquakes lacking local arrival time data, and it minimizes the need for analyst review of depth-phase picks ( pP , sP ) or using predefined ‘fixed’ depths.

Seismological Research Letters

Factors influencing landslide occurrence in low-relief formerly glaciated landscapes: Landslide inventory and susceptibility analysis in Minnesota, USA

In landscapes recently impacted by continental glaciation, landslides may occur where topographic relief has been generated by the drainage of glacial lakes and ensuing post-glacial fluvial network development into unconsolidated glacially derived sediments and exhumed bedrock. To investigate linkages among environmental variables, post-glacial landscape development, and landslides, we created a landslide inventory of nearly 10,000 landslides in five regions of the formerly glaciated low-relief state of Minnesota, USA. Multivariate logistic regression indicates the importance of slope angle, lithology, and the development of stream valleys to landslide distribution. Areas underlain by fine-grained glaciolacustrine and nearshore deposits that are incised by streams are particularly prone to shallow (<1-2 m depth) landslides. Landslides also occur in a wide range of glacial and fluvial deposits, and as rockfall in layered Paleozoic sedimentary rocks in central and southern Minnesota and Precambrian igneous and sedimentary rocks in northeastern Minnesota. Although no more than 1-2% of the studied regions are susceptible to landslides, they can pose risk to life and safety, damage infrastructure, and impact water quality. The combination of recently generated low-relief steep slopes, extensive unconsolidated sediments, and layered sedimentary bedrock make this formerly glaciated landscape more susceptible to landslides than current national-scale models indicate.

Minnesota

Gulf Coast Basin CORE-CM initiative final report

The Bureau of Economic Geology at the University of Texas at Austin (UT-BEG) is leading the Gulf Coast Carbon Ore, Rare Earth, and Critical Minerals (CORE-CM) Initiative to assess the potential to produce critical minerals (CMs), including rare earth elements (REEs) from coal, coal ash, and produced water related to oil and gas production, and related materials (alumina processing waste [red mud], heavy mineral sands, graphite, and zeolite) within the Gulf Coast Basin. This project represents the first phase in a long-term program and provides reconnaissance data that will be foundational for future work by assessing resources and suggesting plans to be conducted in future work and expanding stakeholder engagement. The project includes several tasks designed to identify, characterize, and assess several necessary aspects for development of CMs and REEs in the Gulf Coast Basin.

Gulf Coast basin

Melt generation sources and conditions in the wake of a migrating slab window: Geochemistry and petrology of the million-year history of primitive volcanism at Clear Lake volcanic field, California

Clear Lake volcanic field (CLVF) is the northernmost and youngest (~2.2 Ma to 8 ka) of the volcanic centers distributed along the San Andreas transform fault in western California. The initial phase of CLVF volcanism (interval one) occurred between ~2.2 and 1.3 Ma and extends ~35 km southeast of Clear Lake, forming a semi-continuous upland plateau capped by lava flows, with isolated volcanic remnants on the periphery. This volcanism is broadly characterized by geochemically primitive compositions that reflect three source compositions and conditions of melt generation. (1) Partial melting of upwelling asthenospheric mantle lherzolite at moderate pressures (1.2–1.4 GPa) and temperatures (1297–1329 °C) produced high-CaO (9.8–11.3 wt %) basalts with high Al 2 O 3 (16.8–17.6 wt %), Mg#s (66–70), MgO (8–10 wt %), Ni (103–262 μg/g), and Cr (284–609 μg/g). These high-CaO basalts contain olivine (Fo 87–91 ) phenocrysts with Cr-spinel inclusions ± subordinate plagioclase and crop out only in the southern part of the CLVF. (2) Partial melting of depleted sub-continental lithospheric mantle harzburgite at variable pressures (0.7–1.5 GPa) and temperatures (1097–1299 °C) produced a compositional continuum of med-K 2 O, calc-alkaline, high-MgO basalts through high-MgO andesites with high Mg#s (67–77), MgO (8–14 wt %) and high Ni and Cr abundances (154–439 and 340–1124 μg/g, respectively). Mineral assemblages are olivine (Fo 88–93 ) with Cr-spinel inclusions ± subordinate clinopyroxene, orthopyroxene and plagioclase. Small (<2.5 cm) mantle harzburgite xenoliths and mantle olivine xenocrysts are also found in several of these samples. These high-MgO basalts through andesites represent the largest volume of primitive compositions and have erupted predominantly along the main, fault-controlled northwest-southeast trending axis of volcanism with peripheral outcrops to the north, west, and east. (3) Partial melting of the Gorda eclogite slab edge produced adakitic silicic slab melts with strong depletion in the heavy rare earth elements (Yb = 0.6 μg/g). Subsequent reaction of those melts with depleted ultramafic rocks during ascent imprinted the adakitic dacites with high Mg#s (65–78) and elevated Ni (117–210 μg/g) and Cr (191–283 μg/g). Phenocrysts of orthopyroxene (En 87–94 ) with spinel inclusions (Cr# = 80–88) and extremely Ni-rich (9483 μg/g) olivine cores (Fo 84–93 ) record those reactions. Small-volume outcrops of the adakites on the eastern periphery of the CLVF track the passing slab edge. The trio of melting sources recorded by early CLVF magmatism reflect the tectonically complex environment and the hot (1097–1329 °C), shallow (0.7–1.5 GPa) melting conditions for these primitive compositions and provide estimates of the heat delivered to the crust. Over time, this flux led to maturation of the CLVF magmatic system toward the more voluminous and silicic volcanism that characterizes the balance of its subsequent volcanic history and maintains the present-day anomalously high heat flow in the region. The current interval (interval four) of volcanic activity at CLVF is characterized by low-volume, fault-controlled eruptions of basaltic andesite and andesite suggestive of mantle magma and heat delivery to the crust, similar to interval one. This analogous activity provides motivation for the current study and begs the question of whether the system is undergoing thermal priming for renewed silicic volcanism.

California

Geophysical characterization of mine influenced groundwater and surface water in the Mayflower section of the Animas River, Bonita Peak Mining District, Silverton Colorado

This report details findings from geophysical investigations to identify possible groundwater - surface water interactions near the Mayflower Section of the Animas River in Silverton, Colorado. The Mayflower Section is a mine influenced Superfund Site in Colorado. This investigation utilized electromagnetic induction (EM or EMI), magnetic, and fiber optic distributed temperature system (FODTS) geophysical methods to measure the bulk earth electrical conductivity, magnetic susceptibility, and temperature of specific surveyed volumes of the earth near the Mayflower Section of the Animas River. These physical parameters are used to understand the groundwater – surface water interactions, which can guide decision makers in their assessment of mine-impacted surface water. This report details the results from characterization and monitoring technologies to provide high data density and continuous monitoring of bulk earth electrical conductivity and temperature in the Mayflower section of the Animas River to identify zones of groundwater – surface water interactions and potential metal loading from mine influenced water. The investigation separated right and left bank characterization for each method and indicates more groundwater is entering from the right bank than the left bank and these predominantly right bank discharges potentially contain metal-rich water compared to the left bank. Results also indicate mineral veins facilitate preferential groundwater discharge to the river due to possible jointing, fractures, and permeability differences sometimes occurring along veins relative to host rock. For example, Boulder Gulch is likely groundwater dominated and may be receive water impacted by Mayflower Impoundments #1 and #2. Additionally, the beaver ponds near Blair Gulch may influence groundwater discharge and Mayflower Impoundment #4 is possibly impacting groundwater and surface water and may be connected to the wetlands to the west of the impoundment. These data could be further analyzed for smaller spatial scale analysis within areas of interest. The identification of these locations along sections of the river likely impacted by mine influenced groundwater potentially entering the Animas River and can be used by site investigators, decision makers, and stakeholders in mitigation decisions and strategies.

Colorado

Forecasting volcanic activity in Germany—A multi-criteria approach

Igneous activity, including shallow intrusions and volcanism, has the potential to disrupt underground critical infrastructure. Notably, future underground infrastructure projects like high-level radioactive waste repositories must be sited in areas of extremely low disruption probability by igneous activity. In Germany, according to the Repository Site Selection Act of 2017 (Standortauswahlgesetz, or StandAG), areas in which Quaternary volcanism is either present or future volcanic activity is expected within the next 1 million years (m.y.) must be excluded from the site selection process. Although the locations of regions with Quaternary volcanism are reasonably well known in Germany, forecasting potential igneous activity at intraplate volcanic fields is challenging, as many processes and their interactions control the spatial distribution of volcanic centers. Here, a semi-quantitative, multi-criteria approach is proposed for a regional evaluation of the relative potential of future igneous activity in Germany. A variety of geoscientific indicators are used, including seismic anomalies in Earth’s mantle, gravity data, tectonic activity, sutures, ground motion, earthquakes, mantle degassing centers, and geochronological data of volcanic rocks. The indicators describe the sequence of processes from potential melt generation in Earth’s mantle, through ascent and accumulation of melt within the lithosphere, to eruption at Earth’s surface. In total, 15 out of 30 proposed geoscientific indicators are selected and quantified using 20 total assigned parameters. Defined threshold values are used to spatially delimit relevant parameter properties to focus on areas with higher potential of future magmatic activity. To consider uncertainties of parameters and their underlying processes, which are usually more spatially extensive below ground, buffer zones are defined in which values of relevance decrease with increasing distance from the initial lateral shape of a parameter. Normalized parameters are combined into an index, whose spatial value distribution is used to differentiate the relative potential of future igneous activity (within the next 1 m.y.). The sensitivity of the results is shown by varying the weighting factors for the relevant parameters in country-wide index maps. Thereby, profiles illustrate the distribution of the resulting index values and respective index fractions of various parameters. Different index maps for the relative potential of future igneous activity are presented and can be used for hazard assessments.

Professional Paper