Search USGS⌕ Search

SEARCH · Search USGS

Results for “Earth Science Frontiers”

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 91 records · Page 5Linked to original sources
USGS2021

A multidisciplinary investigation into the eruptive style, processes, and duration of a Cascades back-arc tholeiitic basalt: A case study of the Brushy Butte flow field, northern California, United States

The Cascades back-arc in northern California is dominated by monogenetic tholeiitic basalts that erupted throughout the Pleistocene. Elucidating their eruptive history and processes is important for understanding potential future eruptions here. We focus on the well-exposed monogenetic volcano that emplaced the Brushy Butte flow field, which constructed a ∼150 m tall edifice, has flow lobes up to >10 km long, and in total covers ∼150 km 2 with an eruptive volume of 3.5 km 3 . We use a multidisciplinary approach of field mapping, petrography, geochemistry, paleomagnetism, geochronology, and lidar imagery to unravel the eruptive history and processes that emplaced this flow field. Tholeiitic basalts in northern California have diverse surface morphology and vegetation cover but similar petrographic appearances, which makes them hard to distinguish in the field. Geochemistry and paleomagnetism offer an independent means of distinguishing tholeiitic basalts. Brushy Butte flow field lavas are similar in major-oxide and trace-element abundances but differ from adjacent tholeiitic basalts. This is also apparent in remanent magnetic directions. Additionally, paleomagnetism indicates that the flow field was emplaced during a geologically brief time interval (10–20 years), which 36 Cl cosmogenic dating puts at 35.7 ± 1.7 ka. Lidar imagery shows that these flows erupted from at least 28 vents encompassing multiple scoria cones, spatter cones, and craters. Flows can be grouped into four pulses using stratigraphic position and volume. Pulse 1 is the most voluminous, comprising eight eruptions and ∼2.3 km 3 . Each subsequent pulse started rapidly but decayed quickly, and each successive pulse erupted less lava (i.e., 2.3 km 3 for pulse 1, 0.6 km 3 for pulse 2, 0.3 km 3 for pulse 3, and 0.2 km 3 for pulse 4). Many of these flows host well-established lava channels and levees (with channel breakouts) that lead to lava fans, with some flows hosting lava ponds. Similar flow features from tholeiitic eruptions elsewhere demonstrate that these morphologies generally occur over weeks, months, or longer (e.g., Puʻu ʻŌʻō eruption at K–llauea, Hawaiʻi). This multidisciplinary study shows the range of eruptive styles and durations of a Cascades back-arc eruption and illustrates how potential future tholeiitic eruptive activity in the western United States might progress.

California↗
USGS2021

Recency of faulting and subsurface architecture of the San Diego Bay pull-apart basin, California, USA

In southern California, plate boundary motion between the North American and Pacific plates is distributed across several sub-parallel fault systems. The offshore faults of the California Continental Borderland (CCB) are thought to accommodate ~10-15% of the total plate boundary motion, but the exact distribution of slip and the mechanics of slip partitioning remain uncertain. The Newport-Inglewood-Rose Canyon fault is the easternmost fault within the CCB whose southern segment splays out into a complex network of faults beneath San Diego Bay. A pull-apart basin model between the Rose Canyon and the offshore Descanso fault has been used to explain prominent fault orientations and subsidence beneath San Diego Bay; however this model does not account for faults in the southern portion of the bay or faulting east of the bay. To investigate the characteristics of faulting and stratigraphic architecture beneath San Diego Bay, we combined a suite of reprocessed legacy airgun multi-channel seismic profiles and high-resolution Chirp data, with age and lithology controls from geotechnical boreholes and shallow sub-surface vibracores. This combined dataset is used to create gridded horizon surfaces, fault maps, and perform a kinematic fault analysis. The structure beneath San Diego Bay is dominated by down-to-the-east motion on normal faults that can be separated into two distinct groups. The strikes of these two fault groups can be explained with a double pull-apart basin model for San Diego Bay. In our conceptual model, the western portion of San Diego Bay is controlled by a right-step between the Rose Canyon and Descanso faults, which matches both observations and predictions from laboratory models. The eastern portion of San Diego Bay appears to be controlled by an inferred step-over between the Rose Canyon and San Miguel-Vallecitos faults and displays distinct fault strike orientations, which kinematic analysis indicates should have a significant component of strike-slip partitioning that is not detectable in the seismic data. The potential of a Rose Canyon-San Miguel-Vallecitos fault connection would effectively cut the stepover distance in half and have important implications for the seismic hazard of the San Diego-Tijuana metropolitan area (population ~3 million people).

California↗
USGS2021

Monochromatic long-period seismicity prior to the 2012 earthquake swarm at Little Sitkin Volcano, Alaska

Detection of the earliest stages of unrest is one of the most challenging and yet critically needed aspects of volcano monitoring. We investigate a sequence of five unusual long-period (LP) earthquakes that occurred in the days prior to the onset of a months-long volcano-tectonic (VT) earthquake swarm beneath Little Sitkin volcano in the Aleutian Islands during late 2012. The long-period earthquakes had two distinctive characteristics: their signals were dominated by a monochromatic spectral peak at approximately 0.57 Hz and they had impulsive P and S-wave arrivals on a seismometer located on Amchitka Island 80 km to the southeast of the volcano. In each case, the monochromatic earthquakes ended with a higher-frequency event after approximately 2 min of duration. We find evidence that the five monochromatic LP earthquakes resulted from the resonance of a tabular magma body at middle crustal depths (15 km) on the western side of Little Sitkin. Based on the resonant frequency and quality factor of the monochromatic LP earthquakes, we infer the magma body to have a lateral extent of 500 m and a thickness of 9 m. We interpret that a magmatic intrusion excited the monochromatic LP earthquakes and subsequently increased the stress beneath the volcano, leading to the onset of the shallow (<10 km depth) VT swarm five days later.

Alaska↗
USGS2022

Selected crater and small caldera lakes in Alaska: Characteristics and hazards

This study addresses the characteristics, potential hazards, and both eruptive and non-eruptive role of water at selected volcanic crater lakes in Alaska. Crater lakes are an important feature of some stratovolcanoes in Alaska. Of the volcanoes in the state with known Holocene eruptive activity, about one third have summit crater lakes. Also included are two volcanoes with small caldera lakes (Katmai, Kaguyak). The lakes play an important but not well studied role in influencing eruptive behavior and pose some significant hydrologic hazards. Floods from crater lakes in Alaska are evaluated by estimating maximum potential crater lake water volumes and peak outflow discharge with a dam-break model. Some recent eruptions and hydrologic events that involved crater lakes also are reviewed. The large volumes of water potentially hosted by crater lakes in Alaska indicate that significant flowage hazards resulting from catastrophic breaching of crater rims are possible. Estimates of maximum peak flood discharge associated with breaching of lake-filled craters derived from dam-break modeling indicate that flood magnitudes could be as large as 10 3 –10 6 m 3 /s if summit crater lakes drain rapidly when at maximum volume. Many of the Alaska crater lakes discussed are situated in hydrothermally altered craters characterized by complex assemblages of stratified unconsolidated volcaniclastic deposits, in a region known for large magnitude (>M7) earthquakes. Although there are only a few historical examples of eruptions involving crater lakes in Alaska, these provide noteworthy examples of the role of external water in cooling pyroclastic deposits, acidic crater-lake drainage, and water-related hazards such as lahars and base surge.

Alaska↗
USGS2022

Submarine landslide susceptibility mapping in recently deglaciated terrain, Glacier Bay, Alaska

Submarine mass wasting events have damaged underwater structures and propagated waves that have inundated towns and affected human populations in nearby coastal areas. Susceptibility to submarine landslides can be pronounced in degrading cryospheric environments, where existing glaciers can provide high volumes of sediment, while cycles of glaciation and ice-loss can damage and destabilize slopes. Despite their contribution to potential tsunami hazard, submarine landslides can be difficult to study because of limited access and data collection in underwater environments. Here we present a method to quantify and map the submarine landslide susceptibility of sediment-covered slopes in Glacier Bay, Glacier Bay National Park and Preserve, Alaska, using multibeam-sonar bathymetric digital elevation models (DEMs) and historical maps of glacial extents over the last ∼250 years. After mapping an inventory of >7,000 landslide scarps in submarine sediments, we filtered the inventory by size to account for limitations in DEM resolution and spatial scales relevant to tsunami hazards. We then assessed landslide concentration, accounting for the age of the initial exposure of submarine slopes by deglaciation. We found a positive correlation between landslide concentration and deglaciation age, which we interpreted as a mean landslide accumulation rate over the period of record. Local deviations from this rate indicated differences in susceptibility. Additionally, we accounted for some of the effect of material and morphometric properties by estimating the submarine bedrock-sediment distribution using a morphometric model and assessing the relationship between slope angle and landslide incidence. Finally, we supplemented our susceptibility assessment with a geomorphic component based on the propensity of active submarine fans and deltas to produce landslides. Thus, our map of submarine landslide susceptibility incorporates three components: age-adjusted landslide concentration, slope angle, and geomorphology. We find that areas of mapped high susceptibility correlate broadly with areas of high sediment input and availability, locations of fans and deltas, and steep sediment-covered glacially carved fjords and troughs. Areas of high submarine landslide susceptibility in Glacier Bay moderately correspond with locations of known high-hazard subaerial slopes, but more research on submarine and subaerial landslides in degrading cryospheric environments would be beneficial to better understand landslide and tsunami hazards.

Alaska↗
USGS2022

Depositional controls on detrital zircon provenance: An example from upper Cretaceous strata, southern Patagonia

Understanding how depositional environments within a sedimentary system redistribute and sequester sediment is critical for interpreting basin-scale provenance trends. However, sedimentary source-to-sink models commonly examine temporal changes and do not consider how variation in sedimentation processes across a dispersal pathway may result in contrasting provenance signatures. In this paper, we demonstrate a down-paleoslope shift in detrital zircon provenance signatures correlated with shallow-marine lithofacies patterns from the Upper Cretaceous La Anita Formation and underlying continental slope lithofacies of the Alta Vista Formation (Magallanes-Austral Basin, southern Patagonia). New stratigraphic, sedimentologic, and lithofacies analysis results from the La Anita Formation suggest an upward shoaling succession, from a (i) storm-influenced shoreface, (ii) fluvially-dominated, wave-influenced delta, and a (iii) high-energy, gravelly foreshore. Stratigraphic sections are paired with U-Pb detrital zircon sandstone samples (N = 20; n = 5219), which provide both maximum depositional ages and provenance characteristics. While all samples contain abundant zircon derived from the Andean volcanic arc (ca. 145–75 Ma), the amount from both Jurassic distal volcanic massifs (ca. 188–162 Ma) and recycled orogenic sources exhumed during the advance of the Cretaceous fold-and-thrust belt (>200 Ma; 157–142 Ma) vary with changes in depositional environment. We argue that down-paleoslope, systematic enriching of local fold-and-thrust belt material within the La Anita Formation is reflective of progressive mixing of grains transported via shallow-marine processes, while distally enriched fluvio-deltaic transported zircons were sourced from large, regional catchments. This suggests that competition between transport processes across a shallow and marginal marine sequence of rocks affects the resulting provenance signatures recorded within a single stratigraphic succession. These data also detail the degree of sediment pathway connectivity between shallow-marine sources and deep-marine sinks. Detrital zircon results from muddy continental slope facies of the Alta Vista Formation are made up entirely locally derived material, while zircon results from deep-water, sand-rich channel facies of the Formation are indistinguishable from coeval fluvio-deltaic zircon signatures. This implies that continental shelf-to-slope connectivity in a sediment dispersal system, via submarine canyons or shelf-edge delta progradation, is necessary for detrital zircon distributions from the shallow-marine realm to propagate into the deeper marine.

Patagonia↗
USGS2022

Barkley Canyon gas hydrates: A synthesis based on two decades of seafloor observation and remote sensing

Barkley Canyon is one of the few known sites worldwide with the occurrence of thermogenic gas seepage and formation of structure-II and structure-H gas hydrate mounds on the seafloor. This site is the location of continuous seafloor monitoring as part of the Ocean Networks Canada (ONC) cabled observatory off the west coast off Vancouver Island, British Columbia, Canada. We combine repeat remotely operated vehicle (ROV) seafloor video observations, mapping with an autonomous underwater vehicle (AUV), ship-, ROV-, and AUV-based identification of gas flares, as well as seismic and Chirp data to investigate the distribution of fluid migration pathways. Geologically, the site with the prominent gas hydrate mounds and associated fluid seepage is covering an area of ∼0.15 km 2 and is situated on a remnant of a rotated fault block that had slipped off the steep flanks of the north-east facing canyon wall. The gas hydrate mounds, nearly constant in dimension over the entire observation period, are associated with gas and oil seepage and surrounded by debris of chemosynthetic communities and authigenic carbonate. The formation of gas hydrate at and near the seafloor requires additional accommodation space created by forming blisters at the seafloor that displace the regular sediments. An additional zone located centrally on the rotated fault block with more diffuse seepage (∼0.02 km 2 in extent) has been identified with no visible mounds, but with bacterial mats, small carbonate concretions, and clam beds. Gas venting is seen acoustically in the water column up to a depth of ∼300 m. However, acoustic water-column imaging during coring and ROV dives showed rising gas bubbles to much shallower depth, even <50 m, likely a result of degassing of rising oil droplets, which themselves cannot be seen acoustically. Combining all observations, the location of the gas hydrate mounds is controlled by a combination of fault-focused fluid migration from a deeper reservoir and fluid seepage along more permeable strata within the rotated slope block. Fluids must be provided continuously to allow the sustained presence of the gas hydrate mounds at the seafloor.

Barkley Canyon↗
USGS2022

HydroBench: Jupyter supported reproducible hydrological model benchmarking and diagnostic tool

Evaluating whether hydrological models are right for the right reasons demands reproducible model benchmarking and diagnostics that evaluate not just statistical predictive model performance but also internal processes. Such model benchmarking and diagnostic efforts will benefit from standardized methods and ready-to-use toolkits. Using the Jupyter platform, this work presents HydroBench, a model-agnostic benchmarking tool consisting of three sets of metrics: 1) common statistical predictive measures, 2) hydrological signature-based process metrics, including a new time-linked flow duration curve and 3) information-theoretic diagnostics that measure the flow of information among model variables. As a test case, HydroBench was applied to compare two model products (calibrated and uncalibrated) of the National Hydrologic Model - Precipitation Runoff Modeling System (NHM-PRMS) at the Cedar River watershed, WA, United States. Although the uncalibrated model has the highest predictive performance, particularly for high flows, the signature-based diagnostics showed that the model overestimates low flows and poorly represents the recession processes. Elucidating why low flows may have been overestimated, the information-theoretic diagnostics indicated a higher flow of information from precipitation to snowmelt to streamflow in the uncalibrated model compared to the calibrated model, where information flowed more directly from precipitation to streamflow. This test case demonstrated the capability of HydroBench in process diagnostics and model predictive and functional performance evaluations, along with their tradeoffs. Having such a model benchmarking tool not only provides modelers with a comprehensive model evaluation system but also provides an open-source tool that can further be developed by the hydrological community.

Frontiers in Earth Sciences↗
USGS2022

Modflow-setup: Robust automation of groundwater model construction

In an age of both big data and increasing strain on water resources, sound management decisions often rely on numerical models. Numerical models provide a physics-based framework for assimilating and making sense of information that by itself only provides a limited description of the hydrologic system. Often, numerical models are the best option for quantifying even intuitively obvious connections between human activities and water resource impacts. However, despite many recent advances in model data assimilation and uncertainty quantification, the process of constructing numerical models remains laborious, expensive, and opaque, often precluding their use in decision making. Modflow-setup aims to provide rapid and consistent construction of MODFLOW groundwater models through robust and repeatable automation. Common model construction tasks are distilled in an open-source, online code base that is tested and extensible through collaborative version control. Input to Modflow-setup consists of a single configuration file that summarizes the workflow for building a model, including source data, construction options, and output packages. Source data providing model structure and parameter information including shapefiles, rasters, NetCDF files, tables, and other (geolocated) sources to MODFLOW models are read in and mapped to the model discretization, using Flopy and other general open-source scientific Python libraries. In a few minutes, an external array-based MODFLOW model amenable to parameter estimation and uncertainty quantification is produced. This paper describes the core functionality of Modflow-setup, including a worked example of a MODFLOW 6 model for evaluating pumping impacts to a lake in central Wisconsin, United States.

Wisconsin↗
USGS2022

Stratigraphy and eruption history of maars in the Clear Lake Volcanic Field, California

The Clear Lake Volcanic Field (CLVF) is the northernmost and youngest field in a chain of volcanic provinces in the California Coast Range mountains. Effusive and explosive volcanic activity in the field has spanned at least 2.1 million years, with the youngest eruptions comprising a series of maar craters at the edges of, and within, Clear Lake itself. This work documents the first direct ages for many of these maar deposits, and builds the stratigraphic basis for interpreting eruptive processes and dynamics of the young eruptions which produced them. Detailed stratigraphy has distinguished maar eruption products from pyroclastic deposits (monolithologic falls and flows, previously mapped together with maars as a single unit), and established a set of 6 eruption facies from maar deposit lithology, grain size parameters, and depositional structures. Radiocarbon dates from carbon films found on clasts at 3 outcrops have constrained several of these maar eruptions to ~8500-13,500 years BP, coinciding with eruptive periods previously estimated based on lake core tephrachronology. Part of this period also coincides with indigenous occupation (< 12,000 years BP), which suggests that oral histories of Pomo and other local tribes may contain descriptions of volcanic phenomena experienced by local residents of the CLVF. Collaboration between volcanologists and indigenous historians may add a valuable human dimension to the youngest eruptions of the Clear Lake Volcanic Field; combined, geologic and ethnographic avenues of research will help build a richer eruption history for future volcanic hazard assessment.

California↗
USGS2022

Forecasting and communicating the dispersion and fallout of ash during volcanic eruptions: lessons from the September 20, 2020 eruptive pulse at Sangay volcano, Ecuador

Volcanic ash is a hazard to human health and activities. Accurate and timely forecasts coupled with effective communication about the dispersion and fallout of volcanic ash during explosive events are essential to reduce impacts on local communities and limit economic losses. In this study, we present the first detailed description of an eruptive pulse at Sangay volcano and its eruption source parameters. The eruptive pulse on September 20, 2020, started at about 9:20 UTC and lasted between 90 and 100 min, producing an eruptive column that rapidly separated into (1) a higher (15.2 km above sea level, asl), gas-rich cloud moving east-southeast and (2) a lower (12.2 km asl), ash-rich cloud moving west and causing ash fallout up to 280 km from the volcano. Field data collected immediately after the event allow estimating the volume of bulk tephra to be between 1.5 and 5.0 × 10 6 m 3 , corresponding to a volcanic explosivity index of 2. The eruptive pulse, identified as violent Strombolian, emitted andesitic ash that was more mafic than products ejected by Sangay volcano in recent decades. Component analysis and glass chemistry of juvenile particles support the hypothesis that this event excavated deeper into the upper conduit compared to typical Strombolian activity at Sangay volcano, while grain-size analysis allows reconstruction of the total grain-size distribution of the fallout deposit. The discrepancies between the ash fallout simulations performed with the Ash3D online tool and the actual deposit are mainly the result of inaccurate pre- and syn-eruptive configurations and highlight the importance of additional processes such as aggregation. Communication products issued during the event included (1) several standard short reports, (2) volcano observatory notices for aviation, (3) social media posts, and (4) a special report providing the results of the ash fallout simulation. Although communication was effective with the authorities and the connected population, an effort must be made to reach the most vulnerable isolated communities for future events. This study shows that forecasting ash dispersion and fallout during volcanic eruptions can guide early warnings and trigger humanitarian actions, and should become a standard in volcano observatories worldwide.

Sangay volcano↗
USGS2022

A study on the effect of site response on California seismic hazard map assessment

Prior studies have repeatedly shown that probabilistic seismic hazard maps from several different countries predict higher shaking than that observed. Previous map assessments have not, however, considered the influence of site response on hazard. Seismologists have long acknowledged the influence of near-surface geology, in particular low-impedance sediment layers, on earthquake ground-motion at frequencies of engineering concern. Although the overall effects of site response are complex, modern ground-motion models (GMMs) account for site effects using terms based on V S30 , the time-averaged shear-wave velocity in the upper 30 m of the Earth’s surface. In this study, we consider general implications of incorporating site terms from modern GMMs using site-specific V S30 as a proxy in probabilistic seismic hazard maps for California. At the long periods (1–5 s) that affect tall buildings, site terms amplify the mapped hazard by factors of 1–3 at many sites relative to maps calculated for the standard reference soft-rock site condition, V S30 = 760 m/s. However, at the short periods of ground-motion that are the main contributors to peak ground acceleration (PGA) and thus affect smaller structures, only negligible effects occur due to nonlinear deamplification of strong ground-motion at high frequencies. Nonlinear deamplification increases as the shaking level increases. For very strong shaking, deamplification can overcome the linear amplification, yielding net deamplification. We explore the implications of these results for the evaluation of hazard maps. Because site effects do not change the maps appreciably at short periods, we can exclude site response as an explanation for why the maps overpredict historically observed shaking as captured by the California Historical Intensity Mapping Project (CHIMP) dataset. The results are expected to be generalizable to regions that are comparable to California in terms of structure and seismicity rates. In low-to-moderate-seismicity regions where the hazard reflects weaker shaking, nonlinear site response is expected to be less important for the hazard.

California↗
USGS2022

Crustal permeability changes observed from seismic attenuation: Impacts on multi-mainshock sequences

We use amplitude ratios from narrowband-filtered earthquake seismograms to measure variations of seismic attenuation over time, providing unique insights into the dynamic state of stress in the Earth’s crust at depth. Our dataset from earthquakes of the 2016-2017 Central Apennines sequence allows us to obtain high-resolution time histories of seismic attenuation (frequency band: 0.5-30 Hz) characterized by strong earthquake dilatation-induced fluctuations at seismogenic depths, caused by the cumulative elastic stress drop after the sequence, as well as damage-induced ones at shallow depths caused by energetic surface waves. Cumulative stress drop causes negative dilatation, reduced permeability, and seismic attenuation, whereas strong-motion surface waves produce an increase in crack density, and so in permeability and seismic attenuation. In the aftermath of the main shocks of the sequence, we show that the M ≥ 3.5 earthquake occurrence vs. time and distance is consistent with fluid diffusion: diffusion signatures are associated with changes in seismic attenuation during the first days of the Amatrice, Visso-Norcia, and Capitignano sub-sequences. We hypothesize that coseismic permeability changes create fluid diffusion pathways that are at least partly responsible for triggering multi-mainshock seismic sequences. Here we show that anelastic seismic attenuation fluctuates coherently with our hypothesis.

Frontiers in Earth Science↗
USGS2022

Forecasting explosions at Sinabung Volcano, Indonesia, based on SO2 emission rates

Dome-building volcanic eruptions are often associated with frequent Vulcanian explosions, which constitute a substantial threat to proximal communities. One proposed mechanism driving such explosions is the sealing of the shallow volcanic system followed by pressurization due to gas accumulation beneath the seal. We investigate this hypothesis at Sinabung Volcano (Sumatra, Indonesia), which has been in a state of eruption since August 2010. In 2013, the volcano began erupting a lava dome and lava flow, and frequent explosions produced eruptive columns that rose many kilometers into the atmosphere and at times sent pyroclastic density currents down the southeast flanks. A network of scanning Differential Optical Absorption Spectrometers (DOAS) was installed on the volcano’s eastern flank in 2016 to continuously monitor SO 2 emission rates during daytime hours. Analysis of the DOAS data from October 2016 to September 2017 revealed that passive SO 2 emissions were generally lower in the 5 days leading up to explosive events (∼100 t/d) than was common in 5-day periods leading up to days on which no explosions occurred (∼200 t/d). The variability of passive SO 2 emissions, expressed as the standard deviation, also took on a slightly wider range of values before days with explosions (0–103 t/d at 1-sigma) than before days without explosions (43–117 t/d). These observations are consistent with the aforementioned seal-failure model, where the sealing of the volcanic conduit blocks gas emissions and leads to pressurization and potential Vulcanian explosions. We develop a forecasting methodology that allows calculation of a relative daily explosion probability based solely on measurements of the SO 2 emission rate in the preceding days. We then calculate forecast explosion probabilities for the remaining SO 2 emissions dataset (October 2017—September 2021). While the absolute accuracy of forecast explosion probabilities is variable, the method can inform the probability of an explosion occurring relative to that on other days in each test period. This information can be used operationally by volcano observatories to assess relative risk. The SO 2 emissions-based forecasting method is likely applicable to other open vent volcanoes experiencing dome-forming eruptions.

Frontiers in Earth Science↗
USGS2023

Estimates of volcanic mercury emissions from Redoubt Volcano, Augustine Volcano, and Mount Spurr eruption ash

Ash is a potential sink of volcanically sourced atmospheric mercury (Hg), and the concentration of particle-bound Hg may provide constraints on Hg emissions during eruptions. We analyze Hg concentrations in 227 bulk ash samples from the Mount Spurr (1992), Redoubt Volcano (2009), and Augustine Volcano (2006) volcanic eruptions to investigate large-scale spatial, temporal, and volcanic-source trends. We find no significant difference in Hg concentrations in bulk ash by distance or discrete eruptive events at each volcano, suggesting that in-plume reactions converting gaseous Hg 0 to adsorbed Hg 2+ are happening on shorter timescales than considered in this study (minutes) and any additional in-plume controls are not discernable within intra-volcanic sample variability. However, we do find a significant difference in Hg concentration of ash among volcanic sources, which indicates that volcanoes may emit comparatively high or low quantities of Hg. We combine our Hg findings with total mass estimates of ashfall deposits to calculate minimum, first-order Hg emissions of 8.23 t Hg for Mount Spurr (1992), 1.25 t Hg for Redoubt Volcano (2009), and 0.16 t Hg for Augustine Volcano (2006). In particular, we find that Mount Spurr is a high Hg emitting volcano, and that its 1992 particulate Hg emissions likely contributed substantially to the global eruptive volcanic Hg budget for that year. Based on our findings, previous approaches that use long-term Hg/SO 2 mass ratios to estimate eruptive total Hg under-account for Hg emitted in explosive events, and global volcanogenic Total Hg estimates need revisiting.

Alaska↗
USGS2023

Weak degassing from remote Alaska volcanoes characterized with a new airborne Imaging DOAS instrument and a suite of in situ sensors

Gas emissions from volcanoes occur when volatile species exsolve from magmatic and hydrothermal systems and make their way to the surface. Measurements of emitted gases therefore provide insights into volcanic processes. On 16 July 2021, we made airborne measurements of weak gas plumes emitted from four remote Alaska volcanoes: Iliamna Volcano, Mount Douglas, Mount Martin, and Mount Mageik. Integrated into a small fixed-wing aircraft, a new Imaging Differential Optical Absorption Spectroscopy (DOAS) instrument was used to map the spatial extent of SO 2 plumes as they drifted downwind. Contrary to conventional Mobile DOAS instruments, which provide only a single viewing direction, the Imaging DOAS simultaneously measures SO 2 column density along 48 individual viewing directions oriented in a swath above or below the aircraft. Each of the individual measurements have a comparable precision and sensitivity to those obtained by conventional instruments. Together, they provide high resolution 2D imagery of the volcanic plumes and allow calculation of limited emission rate time series information. Although zenith-facing DOAS measurements achieve greater accuracy and are performed here, the application of the Imaging DOAS in a nadir-facing setup is also discussed and compared to satellite observations made in similar geometries. Also onboard the aircraft, a suite of electrochemical and optical sensors measured the relative abundances of the six major volcanic volatile species H 2 O, CO 2 , SO 2 , H 2 S, HCl, and HF as the aircraft passed through the plumes. Mean SO 2 emission rates of 90 ± 10, 20 ± 3, and 13 ± 3 t/d were measured at Iliamna Volcano, Mount Douglas, and Mount Martin, respectively. SO 2 emissions were below the DOAS detection limit at Mount Mageik but CO 2 and H 2 S could be measured with the in situ sensors. The information gleaned from these measurements was used to assess and compare activity at these volcanoes, all of which were found to be in a state of background degassing but whose emissions pointed to different source conditions ranging from mixed magmatic-hydrothermal to purely hydrothermal in character. Additional measurements at Mount Spurr, Redoubt Volcano, and Augustine Volcano failed to detect the very weak gas concentrations downwind of these persistently degassing vents.

Alaska↗
USGS2023

An algorithm for correction of atmospheric scattering dilution effects in volcanic gas emission measurements using skylight differential optical absorption spectroscopy

Differential Optical Absorption Spectroscopy (DOAS) is commonly used to measure gas emissions from volcanoes. DOAS instruments measure the absorption of solar ultraviolet (UV) radiation scattered in the atmosphere by sulfur dioxide (SO 2 ) and other trace gases contained in volcanic plumes. The standard spectral retrieval methods assume that all measured light comes from behind the plume and has passed through the plume along a straight line. However, a fraction of the light that reaches the instrument may have been scattered beneath the plume and thus has passed around it. Since this component does not contain the absorption signatures of gases in the plume, it effectively “dilutes” the measurements and causes underestimation of the gas abundance in the plume. This dilution effect is small for clean-air conditions and short distances between instrument and plume. However, plume measurements made at long distance and/or in conditions with significant atmospheric aerosol, haze, or clouds may be severely affected. Thus, light dilution is regarded as a major error source in DOAS measurements of volcanic degassing. Several attempts have been made to model the phenomena and the physical mechanisms are today relatively well understood. However, these models require knowledge of the local atmospheric aerosol composition and distribution, parameters that are almost always unknown. Thus, a practical algorithm to quantitatively correct for the dilution effect is still lacking. Here, we propose such an algorithm focused specifically on SO 2 measurements. The method relies on the fact that light absorption becomes non-linear for high SO 2 loads, and that strong and weak SO 2 absorption bands are unequally affected by the diluting signal. These differences can be used to identify when dilution is occurring. Moreover, if we assume that the spectral radiance of the diluting light is identical to the spectrum of light measured away from the plume, a measured clean air spectrum can be used to represent the dilution component. A correction can then be implemented by iteratively subtracting fractions of this clean air spectrum from the measured spectrum until the respective absorption signals on strong and weak SO 2 absorption bands are consistent with a single overhead SO 2 abundance. In this manner, we can quantify the magnitude of light dilution in each individual measurement spectrum as well as obtaining a dilution-corrected value for the SO 2 column density along the line of sight of the instrument. This paper first presents the theory behind the method, then discusses validation experiments using a radiative transfer model, as well as applications to field data obtained under different measurement conditions at three different locations; Fagradalsfjall located on the Reykjanaes peninsula in south Island, Manam located off the northeast coast of mainland Papua New Guinea and Holuhraun located in the inland of north east Island.

Frontiers in Earth Science↗
USGS2023

A new permanent, low-cost, low-power SO2 camera for continuous measurement of volcanic emissions

Since its introduction to volcanology in the mid-2000 s, the SO 2 camera has become an important instrument for the acquisition of accurate and high time-resolution SO 2 emission rates, aiding in hazard assessment and volcanological research. However, with the exception of a few locations (Stromboli, Etna, Kīlauea), hitherto the majority of measurements have been made on discrete field campaigns, which provide only brief snapshots into a volcano’s activity. Here, we present the development of a new, low-cost, low-power SO 2 camera for permanent deployment on volcanoes, facilitating long-term, quasi-continuous (daylight hours only) measurements. We then discuss preliminary datasets from Lascar and Kīlauea volcanoes, where instruments are now in continuous operation. Further proliferation of such instrumentation has the potential to greatly improve our understanding of the transient nature of volcanic activity, as well as aiding volcano monitoring/eruption forecasting.

Frontiers in Earth Science↗