USGS Science⌕ Search

SEARCH · USGS Science

Results for “Frontiers in Earth Science - Volcanology”

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 109 records · Page 6Linked to original sources
USGS2024

Twenty years of explosive-effusive activity at El Reventador volcano (Ecuador) recorded in its geomorphology

Shifts in activity at long-active, open-vent volcanoes are difficult to forecast because precursory signals are enigmatic and can be lost in and amongst daily activity. Here, we propose that crater and vent morphologies, along with summit height, can help us bring some insights into future activity at one of Ecuador’s most active volcanoes El Reventador. On 3 November 2002, El Reventador volcano experienced the largest eruption in Ecuador in the last 140 years and has been continuously active ever since with transitions between and coexistence of explosive and effusive activity, characterized by Strombolian and Vulcanian behavior. Based on the analysis of a large dataset of thermal and visual images, we determined that in the last 20 years of activity, the volcano faced three destructive events: A. Destruction of the upper part of the summit leaving a north-south breached crater (3 November 2002), B. NE border crater collapse (2017), and C. NW flank collapse (2018), with two periods of reconstruction of the edifice: Period 1. Refill of the crater (2002-early 2018) and Period 2. Refill of the 2018 scar (April 2018–December 2022). Through photogrammetric analysis of visual and thermal images acquired in 11 overflights of the volcano, we created a time-series of digital elevation models (DEMs) to determine the maximum height of the volcano at each date, quantify the volume changes between successive dates, and characterize the morphological changes in the summit region. We estimate that approximately 34.1x10 6 m 3 of volcanic material was removed from the volcano due to destructive events, whereas 64.1x10 6 m 3 was added by constructive processes. The pre-2002 summit height was 3,560 m and due to the 2002 eruption it decreased to 3,527 m; it regained its previous height between 2014 and 2015 and the summit crater was completely filled by early April 2018. Event A resulted from an intrusion of magma that erupted violently; we proposed that Events B and C could be a result of an intrusion as well but may also be due to a lack of stability of the volcano summit which occurs when it reaches its maximum height of approximately 3,590 and 3,600 m.

El Reventador↗
USGS2023

Diving deeper into seep distribution along the Cascadia Convergent Margin, USA

Previous margin-wide studies of methane seep distribution along the Cascadia Subduction Zone indicate peaks in seep density within the landward limit of the of gas hydrate stability zone (GHSZ; ≤500 m depth), suggesting a link between current ocean warming, acceleration of hydrate dissociated, and methane emissions. This inferred connection, however, may not account for regional geologic and/or structural complexities driving methane seepage. Expanding upon an existing seep database by adding new seeps data, we conducted statistical and spatial analyses to determine margin-wide distribution trends and offer a tectonic framework for understanding the tendency toward non-normality and spatial clustering. We then highlight the role of local-scale drivers of seep formation in addition to the first-order tectonic framework, using systematic geologic/geomorphic characterization of seep emission sites in southern Cascadia and case studies using meta-attribute analysis of seismic reflection data.. Seep distribution along the margin is non-random, but instead of clustering along the 500-m isobath, regions of high seep density occur in canyons and topographic highs. New findings from this study conclude that co-location of the outer arc high (OAH) and the landward limit of the GHSZ may explain high concentrations of seeps where deformation is the greatest and hydrates are unstable. Detailed analysis of the spatial relationships between seep sites and geologic-geomorphic features in southern Cascadia reveal a link between seeps and anticlines, with 52% of the seeps found in association with anticlines, 36% found at faults, 16% associated with canyons, and 11% at seafloor failure scarps. Given that a majority of anticlines are located along or seaward of the OAH in the actively deforming outer wedge, we suggest that the location of the OAH is a primary structural control on seep distribution. This scenario is supported by neural network analysis of multichannel seismic data revealing zones of probable fluid migration along vertical pipes, faults, and chimneys in the vicinity of active seep sites on anticlines. Determining linkages between seeps and submarine tectonic geomorphology is a crucial first step for understanding and forecasting the distribution of methane seepage, but also a necessity for evaluating causal relationships between ocean warming and gas-hydrate stability.

California, Oregon, Washington↗
USGS2023

Current and future sinkhole susceptibility in karst and pseudokarst areas of the conterminous United States

Sinkholes in karst and pseudokarst regions threaten infrastructure, property, and lives. We mapped closed depressions in karst and pseudokarst regions of the conterminous United States (U.S.) from 10-m-resolution elevation data using high-performance computing, and then created a heuristic additive model of sinkhole susceptibility that also included nationally consistent data for factors related to geology, soils, precipitation extremes, and development. Maps identify potential sinkhole hotspots based on current conditions and projections for 50 years into the future (the years 2070–2079) based on climate change and urban development scenarios. Areas characterized as having either high or very high sinkhole susceptibility contain 94%–99% of known or probable sinkhole locations from three U.S. state databases. States and counties with the highest amounts and percentages of land in zones of highest sinkhole susceptibility are identified. Projected changes in extreme precipitation and development did not substantially change current hotspots of highest sinkhole susceptibility. Results provide a uniform index of sinkhole potential that can support national planning, instead of existing assessments produced through various methods within individual states or smaller areas.

conterminous United States↗
USGS2023

Editorial: Remote sensing of volcanic gas emissions from the ground, air, and space

When magma rises in volcanic systems, volatile species exsolve from the melt and are outgassed to the atmosphere. The melt composition and temperature, depth at which degassing occurs, extent of gas-water-rock interactions, and volume of ascending magma are all factors that determine the composition and rate of gas emissions at the surface. Interpreted in a petrological framework, gas measurements thus provide information on these fundamental parameters of volcanic systems. Volcanic gases have traditionally been sampled in the field and later analyzed with standard laboratory methods, but remote sensing measurements are playing an increasingly central role in characterizing emissions and the volcanoes from which they originate. The 17 contributions in this Research Topic summarize the state-of-the-art in volcanic gas remote sensing and identify key areas in which the field could further improve our understanding of global volcanism and its impact on Earth’s environment in the next decade.

Frontiers in Earth Science↗
USGS2023

A novel non-destructive workflow for examining germanium and co-substituents in ZnS

A suite of complementary techniques was used to examine germanium (Ge), a byproduct critical element, and co-substituent trace elements in ZnS and mine wastes from four mineral districts where germanium is, or has been, produced within the United States. This contribution establishes a comprehensive workflow for characterizing Ge and other trace elements, which captures the full heterogeneity of samples through extensive pre-characterization. This process proceeded from optical microscopy, to scanning electron microscopy and cathodoluminescence (CL) imaging, to electron microprobe analysis, prior to synchrotron-based investigations. Utilizing non-destructive techniques enabled reanalysis, which proved essential for verifying observations and validating unexpected results. In cases where the Fe content was <0.3 wt% in ZnS, cathodoluminescence imaging proved to be an efficient means to qualitatively identify trace element zonation that could then be further explored by other micro-focused techniques. Micro-focused X-ray diffraction was used to map the distribution of the non-cubic ZnS polymorph, whereas micro-focused X-ray fluorescence (μ-XRF) phase mapping distinguished between Ge 4+ hosted in primary ZnS and a weathering product, hemimorphite [Zn 4 Si 2 O 7 (OH) 2 ·H 2 O]. Microprobe data and μ-XRF maps identified spatial relationships among trace elements in ZnS and implied substitutional mechanisms, which were further explored using Ge and copper (Cu) X-ray absorption near-edge spectroscopy (XANES). Both oxidation states of Ge (4+ and 2+) were identified in ZnS along with, almost exclusively, monovalent Cu. However, the relative abundance of Ge oxidation states varied among mineral districts and, sometimes, within samples. Further, bulk XANES measurements typically agreed with micro-focused XANES (μ-XANES) spectra, but unique micro-environments were detected, highlighting the importance of complementary bulk and micro-focused measurements. Some Ge μ-XANES utilized a high energy resolution fluorescence detector, which improved spectral resolution and spectral signal-to-noise ratio. This detector opens new opportunities for exploring byproduct critical elements in complex matrices. Overall, the non-destructive workflow employed here can be extended to other byproduct critical elements to more fully understand fundamental ore enrichment processes, which have practical implications for critical element exploration, resource quantification, and extraction.

Frontiers in Earth Science↗
USGS2024

Automatic identification and quantification of volcanic hotspots in Alaska using HotLINK: The hotspot learning and identification network

An increase in volcanic thermal emissions can indicate subsurface and surface processes that precede, or coincide with, volcanic eruptions. Space-borne infrared sensors can detect hotspots—defined here as localized volcanic thermal emissions—in near-real-time. However, automatic hotspot detection systems are needed to efficiently analyze the large quantities of data produced. While hotspots have been automatically detected for over 20 years with simple thresholding algorithms, new computer vision technologies, such as convolutional neural networks (CNNs), can enable improved detection capabilities. Here we introduce HotLINK: the Hotspot Learning and Identification Network, a CNN trained to detect hotspots with a dataset of −3,800 satellite-based, Visible Infrared Imaging Radiometer Suite (VIIRS) images from Mount Veniaminof and Mount Cleveland volcanoes, Alaska. We find that our model achieves an accuracy of 96% (F1-score 0.92) when evaluated on −1,700 unseen images from the same volcanoes, and 95% (F1-score 0.67) when evaluated on −3,000 images from six additional Alaska volcanoes (Augustine Volcano, Bogoslof Island, Okmok Caldera, Pavlof Volcano, Redoubt Volcano, Shishaldin Volcano). In comparison with an existing threshold-based hotspot detection algorithm, MIROVA (Coppola et al., Geological Society, London, Special Publications, 2016, 426, 181–205), our model detects 22% more hotspots and produces 12% fewer false positives. Additional testing on −700 labeled Moderate Resolution Imaging Spectroradiometer (MODIS) images from Mount Veniaminof demonstrates that our model is applicable to this sensor’s data as well, achieving an accuracy of 98% (F1-score 0.95). We apply HotLINK to 10 years of VIIRS data and 22 years of MODIS data for the eight aforementioned Alaska volcanoes and calculate the radiative power of detected hotspots. From these time series we find that HotLINK accurately characterizes background and eruptive periods, similar to MIROVA, but also detects more subtle warming signals, potentially related to volcanic unrest. We identify three advantages to our model over its predecessors: 1) the ability to detect more subtle volcanic hotspots and produce fewer false positives, especially in daytime images; 2) probabilistic predictions provide a measure of detection confidence; and 3) its transferability, i.e., the successful application to multiple sensors and multiple volcanoes without the need for threshold tuning, suggesting the potential for global application.

Alaska↗
USGS2024

Seismic attenuation and stress on the San Andreas Fault at Parkfield: Are we critical yet?

The Parkfield transitional segment of the San Andreas Fault (SAF) is characterized by the production of frequent quasi-periodical M6 events that break the very same asperity. The last Parkfield mainshock occurred on 28 September 2004, 38 years after the 1966 earthquake, and after the segment showed a ∼22 years average recurrence time. The main reason for the much longer interevent period between the last two earthquakes is thought to be the reduction of the Coulomb stress from the M6.5 Coalinga earthquake of 2 May 1983, and the M6 Nuñez events of June 11th and 22 July 1983. Plausibly, the transitional segment of the SAF at Parkfield is now in the late part of its seismic cycle and current observations may all be relative to a state of stress close to criticality. However, the behavior of the attenuation parameter in the last few years seems substantially different from the one that characterized the years prior to the 2004 mainshock. A few questions arise: (i) Does a detectable preparation phase for the Parkfield mainshocks exist, and is it the same for all events? (ii) How dynamically/kinematically similar are the quasi-periodic occurrences of the Parkfield mainshocks? (iii) Are some dynamic/kinematic characteristics of the next mainshock predictable from the analysis of current data? (e.g., do we expect the epicenter of the next failure to be co-located to that of 2004?) (iv) Should we expect the duration of the current interseismic period to be close to the 22-year “undisturbed” average value? We respond to the questions listed above by analyzing the non-geometric attenuation of direct S-waves along the transitional segment of the SAF at Parkfield, in the close vicinity of the fault plane, between January 2001 and November 2023. Of particular interest is the preparatory behavior of the attenuation parameter as the 2004 mainshock approached, on both sides of the SAF. We also show that the non-volcanic tremor activity modulates the seismic attenuation in the area, and possibly the seismicity along the Parkfield fault segment, including the occurrence of the mainshocks.

California↗
USGS2024

Editorial: From cold seeps to hydrothermal vents: Geology, chemistry, microbiology, and ecology in marine and coastal environments

This Research Topic compiles contemporary studies on cold seeps, hydrothermal vents, mud volcanoes, and related seafloor features that are associated with focused fluid emissions and the transfer of carbon, other chemical species, and sometimes heat from the geosphere to the ocean. Because these features sometimes tap fluids and gas originating kilometers below the seafloor, they provide an important window into deep processes that are otherwise inaccessible to scientists. At the shallow portion of their journey, migrating fluids nearing the seafloor contribute to a range of unique biological, physical, and chemical processes within the sediments themselves and at the sediment-water interface. Seafloor fluid emissions play a critical role in global biogeochemical cycles, ocean chemistry, and possibly even climate change. Seafloor leakage points often emit hydrocarbon gases (especially methane and CO 2 ) and are sometimes the loci for deposition of seafloor minerals that have economic value. A burgeoning area of research focuses on natural products generated at these features, seeking compounds with potential pharmaceutical or other applications. Multidisciplinary studies have become routine for characterization of seafloor fluid emission sites, attesting to the inseparability of geologic, physical, chemical, and biological processes in these settings. It is increasingly common for researchers to combine in a single research cruise: subbottom imaging and seafloor mapping; porewater and water column geochemistry and gas sampling; sediment retrieval for lithologic, biostratigraphic, and solid phase analyses; and studies of benthic and subseafloor communities at the microbial to macrofaunal scales. This multidisciplinary approach has the advantage of ensuring the spatial and temporal coincidence of surveys and samples, an important factor at highly dynamic seafloor fluid emission sites. In addition, researchers often use remotely operated vehicles (ROVs), autonomous underwater vehicles (AUVs), or human-occupied vehicles (HOVs) to record video of the seafloor, compile photomosaics, collect targeted samples, and survey with high-resolution geophysical near-seafloor systems, providing a degree of detail about seafloor fluid emission sites that is unprecedented compared to most areas of the deep ocean. While rarer, long-term cabled observatories or shorter-term deployments of portable observatories are also used at some loci for seafloor fluid flux and are particularly helpful for capturing temporal variations at these dynamic features. Here we summarize the Research Topic’s contribution to multidisciplinary seafloor emission studies in the categories of cold seeps, mud volcanoes, and hydrothermal vents. Figure 1 shows the geographic distribution of the studies in this Research Topic and key features referred to in this Introduction.

Frontiers in Earth Science↗
USGS2025

An exploration of the relative influence of physical models for Omori’s law

Omori’s law states that the rate of aftershocks decays as a function of inverse time. There are multiple physical explanations that we reduce into a nonlinear mixed effects relation of three terms: (1) a Rate/State expression that can account for static/dynamic and viscoelastic triggering caused directly by the mainshock, (2) a fluid diffusion triggering term, and (3) a randomized secondary triggering (cascade) term. We fit free physical-model parameters to an observed aftershock sequence through two nonlinear regression methods to find the relative contributions of physics-based models in an observed aftershock sequence. Results from both methods show that Rate/State models overpredict aftershock rates by ∼0–30%. Secondary aftershocks cause a net negative contribution (seismicity rate reduction that corrects overprediction by other terms) ranging between ∼0 and 30%. All regression solutions yield negative secondary triggering contributions without being guided to do so. A physical explanation for this is that aftershock occurrence relieves stress from the crust, ultimately causing the sequence to extinguish itself. Fluid diffusion triggering contributions range from ∼0 to 20%. Diffusion processes are observed to be shorter in time than the full duration of an aftershock sequence and they are also spatially limited, diminishing their influence. Our results apply to an aftershock decay curve from the 2016 Central Apennines earthquake sequence, meaning that our specific results may not be general. Our primary conclusion is that any one physical model cannot alone fit the observed sequence as well as the combination of three we investigated.

Frontiers in Earth Science↗
USGS2025

Estimation of dynamic geologic CO2 storage resources in the Illinois Basin, including effects of brine extraction, anisotropy, and hydrogeologic heterogeneity

Since the vast majority of carbon dioxide (CO 2 ) storage resources in the United States are in deep saline aquifers, optimizing the use of these saline storage resources could be crucial for efficient development of geologic CO 2 storage (GCS) resources and basin- or larger-scale deployment of GCS in the country. Maximum CO 2 injection rates can be enhanced by extracting brine from the CO 2 storage unit. However, disposal of the extracted brine is both a technological and economic challenge. The lowest-cost option would likely be reinjection of the extracted brine into another formation above or below the CO 2 storage unit. Therefore, it is important to estimate brine injectivity as it will constrain the potential to increase CO 2 injectivity at an injection site that has access to multiple geologic storage units where either CO 2 or brine can be injected. Using a simulation-optimization framework, coupled with a non-isothermal, multiphase CO 2 -water-salt equation-of-state module, we developed a computationally efficient method for evaluating optimization of simultaneous CO 2 injection, brine extraction, and brine (re)injection at hypothetical injection sites deployed across a geologic basin. The Illinois basin is ideal for testing our methodology because it contains multiple geologic storage units with seals in between them to isolate injection of CO 2 in one unit from interfering with the injection of either brine or CO 2 in another unit above or below it. In addition, we investigated the relative effects of variation in key geologic parameters as well as two reservoir structures (hydrogeologic heterogeneity/anisotropy and homogeneity/isotropy) on CO 2 injectivities and enhancement of CO 2 injectivity through extracting brine. Results suggest that permeability, depth, and especially thickness of the storage unit could be the most influential parameters determining CO 2 injectivity. They also suggest that only injecting CO 2 into the storage unit with the greatest injectivity, enhancing that unit’s injectivity by extracting brine, and disposing of the produced brine in other suitable units could maximize total CO 2 injectivity in limited regions of the basin. At the majority of simulated injection sites, however, we found that injecting CO 2 into all of the accessible and suitable storage units was more likely to maximize the CO 2 storage resource.

Illinois, Kentucky, Indiana↗
USGS2026

Pliocene (Piacenzian) planktic foraminiferal biogeography: Insights into how planktic biogeography is linked to climate and socioeconomic impacts on marine ecosystems

The Piacenzian Age (3.60–2.58 Ma) of the Pliocene Epoch was characterized by globally warmer climates, higher sea levels, and atmospheric CO 2 similar to present. Utilizing a robust dataset of 2,101 samples and over 637,000 foraminifer specimens from 77 deep-sea core sites worldwide, we document planktic foraminifer biogeography and biodiversity during the mid-Piacenzian Warm Period (mPWP). Cluster analysis and multidimensional scaling reveal five major bioregions: tropical, warm subtropical, transitional, polar, and a distinct North Atlantic polar bioregion. Each bioregion was dominated by characteristic species with well-established temperature preferences. Analyses demonstrate higher species richness and evenness in low and mid-latitudes, with increased diversity associated with periods of climatic warming and poleward expansion of warm water assemblages. The long-term stability of biogeographic patterns underscores ecological conservatism but also highlights potential vulnerability to rapid anthropogenic climate change. Our findings emphasize the critical role of planktic foraminifers in reconstructing past ocean conditions and offer valuable insights into links between planktic biogeography, climate, and socioeconomic impacts on marine ecosystems. This work advances our understanding of marine ecosystem responses to climate extremes and provides a foundation for future regional and temporal analyses of planktic foraminifer biogeography under global change scenarios.

Frontiers in Earth Science↗
USGS2018

Science at the frontier: Multimethod research to evaluate ecosystem change across multiple scales

Changes in the Earth system occur across the full spectrum of spatial and temporal scales, yet our research approaches to understanding and predicting those changes are typically restricted to a pre-defined window of space and time. For this reason, there is substantial power in integrating different approaches, particularly for research associated with the multifaceted nature of ecosystem responses to global change. Within a given research approach – for example, remote sensing, field experimentation, modeling – science promotes the continued advancement of tools and techniques. As technical advancements continue at an unprecedented rate, new opportunities for integrated, multi-approach research emerge, which could more effectively capture the mechanisms and patterns that drive ecosystem structure and function. A capacity to move beyond comparison and into the realm of integration shows promise for promoting significant advances in Earth system science, as evidenced by the Organized Session Science at the Frontier: Using Multimethod Research to Create New Knowledge and Assess Tools Across Spatial and Temporal Scales on December 12, 2017 at the Fall meeting of the American Geophysical Union in New Orleans, Louisiana.

New Phytologist↗
USGS2016

Seasonal flows of international British Columbia-Alaska rivers: The nonlinear influence of ocean-atmosphere circulation patterns

The northern portion of the Pacific coastal temperate rainforest (PCTR) is one of the least anthropogenically modified regions on earth and remains in many respects a frontier area to science. Rivers crossing the northern PCTR, which is also an international boundary region between British Columbia, Canada and Alaska, USA, deliver large freshwater and biogeochemical fluxes to the Gulf of Alaska and establish linkages between coastal and continental ecosystems. We evaluate interannual flow variability in three transboundary PCTR watersheds in response to El Ni&ntilde;o-Southern Oscillation (ENSO), Pacific Decadal Oscillation (PDO), Arctic Oscillation (AO), and North Pacific Gyre Oscillation (NPGO). Historical hydroclimatic datasets from both Canada and the USA are analyzed using an up-to-date methodological suite accommodating both seasonally transient and highly nonlinear teleconnections. We find that streamflow teleconnections occur over particular seasonal windows reflecting the intersection of specific atmospheric and terrestrial hydrologic processes. The strongest signal is a snowmelt-driven flow timing shift resulting from ENSO- and PDO-associated temperature anomalies. Autumn rainfall runoff is also modulated by these climate modes, and a glacier-mediated teleconnection contributes to a late-summer ENSO-flow association. Teleconnections between AO and freshet flows reflect corresponding temperature and precipitation anomalies. A coherent NPGO signal is not clearly evident in streamflow. Linear and monotonically nonlinear teleconnections were widely identified, with less evidence for the parabolic effects that can play an important role elsewhere. The streamflow teleconnections did not vary greatly between hydrometric stations, presumably reflecting broad similarities in watershed characteristics. These results establish a regional foundation for both transboundary water management and studies of long-term hydroclimatic and environmental change.

Alaska, British Columbia↗