Search USGS⌕ Search

SEARCH · Search USGS

Results for “Frontiers in Earth Science”

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

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

At least 37 records · Page 2Linked to original sources

USGS “Did You Feel It?” — Science and lessons from twenty years of citizen science-based macroseismology

The U.S. Geological Survey (USGS) “Did You Feel It?” (DYFI) system is an automatic method for rapidly collecting macroseismic intensity data from Internet users’ shaking and damage reports and for generating intensity maps immediately following felt earthquakes. DYFI has been in operation for nearly two decades (1999-2019) in the United States, and for nearly 15 years globally. During that period the amount of data collected is astounding: Over 5 million individual DYFI intensity reports—spanning all magnitude and distance ranges—have been amassed and archived. DYFI allows for macroseismic data collection at rates and quantities never before imagined, and thus high-quality macroseismic intensity maps can be made almost immediately, and with more complete coverage at higher resolution than in the past. DYFI also allows for valuable positive interactions of the citizenry with a Federal science agency. In essence, the widespread adoption of DYFI—along with ShakeMap—has facilitated the general acceptance of the very concept of shaking intensity, fundamentally improving our agency’s ability to communicate both hazard and risk to the population. DYFI effectively confirms the importance of reporting and inculcating the public’s understanding of intensity—in addition to magnitude—for a proper perspective of earthquake risk-related decision-making. Furthermore, the vast amount of DYFI data allow for data-rich analyses of otherwise intractable seismological, sociological, and earthquake impact studies that range from quantifying the shaking due to induced earthquakes, to human response and risk perception, to relating recorded shaking metrics to macroseismic effects, and the attenuation of intensity with magnitude and distance. Naturally, web-based data collection also poses challenges. After two decades of experience acquiring data with the DYFI system, we address some of these hallenges by documenting refinements to our algorithmic and operational procedures that have evolved over that time. Lastly, we outline new opportune research and development directions for our DYFI approach to citizen seismology.

Frontiers in Earth Science↗

Aeolian remobilisation of volcanic ash: Outcomes of a workshop in the Argentinian Patagonia

During explosive volcanic eruptions, large quantities of tephra can be dispersed and deposited over wide areas. Following deposition, subsequent aeolian remobilisation of ash can potentially exacerbate primary impacts on timescales of months to millennia. Recent ash remobilisation events (e.g., following eruptions of Cordón Caulle 2011; Chile, and Eyjafjallajökull 2010, Iceland) have highlighted this to be a recurring phenomenon with consequences for human health, economic sectors, and critical infrastructure. Consequently, scientists from observatories and Volcanic Ash Advisory Centers (VAACs), as well as researchers from fields including volcanology, aeolian processes and soil sciences, convened at the San Carlos de Bariloche headquarters of the Argentinian National Institute of Agricultural Technology to discuss the “state of the art” for field studies of remobilised deposits as well as monitoring, modeling and understanding ash remobilisation. In this article, we identify practices for field characterisation of deposits and active processes, including mapping, particle characterisation and sediment traps. Furthermore, since forecast models currently rely on poorly-constrained dust emission schemes, we call for laboratory and field measurements to better parameterise the flux of volcanic ash as a function of friction velocity. While source area location and extent are currently the primary inputs for dispersion models, once emission schemes become more sophisticated and better constrained, other parameters will also become important (e.g., source material volume and properties, effective precipitation, type and distribution of vegetation cover, friction velocity). Thus, aeolian ash remobilisation hazard and associated impact assessment require systematic monitoring, including the development of a regularly-updated spatial database of resuspension source areas.

Frontiers in Earth Science↗

Editorial: Enabling people-centered risk communication for geohazards

In the field of natural hazards, communicating science with the public and stakeholders (i.e., interested parties) involves entering the challenging and complex world of hazard and risk communication, the ultimate purpose of which is to reduce the impact of impending hazards on people and property at risk. Hazard and risk communication are adequate if they reach people with the information that they need, at the right time, and in a form that they can use. This task appears to be particularly difficult when decisions by the public and stakeholders have to be made in the presence of uncertainty about what could happen, as is often the case with geohazards. Moreover, decision-making is complex when there are time pressures, human and economic resources are limited, and multiple sources of information need to be considered. This poses several challenges for the development of two-way and people-centered risk communication for geohazards. The “Enabling People-Centered Risk Communication for Geohazards” Research Topic analyses these challenges and identifies innovative pathways to address them. More precisely, it draws together 13 state-of-the-art articles from around the world on improving communication practices, strategies, and understandings relating to a range of various geohazards and weather-related hazards.

Frontiers in Earth Science↗

Dominant factors in controlling marine gas pools in South China

In marine strata from Sinian to Middle Triassic in South China, there develop four sets of regional and six sets of local source rocks, and ten sets of reservoir rocks. The occurrence of four main formation periods in association with five main reconstruction periods, results in a secondary origin for the most marine gas pools in South China. To improve the understanding of marine gas pools in South China with severely deformed geological background, the dominant control factors are discussed in this paper. The fluid sources, including the gas cracked from crude oil, the gas dissolved in water, the gas of inorganic origin, hydrocarbons generated during the second phase, and the mixed pool fluid source, were the most significant control factors of the types and the development stage of pools. The period of the pool formation and the reconstruction controlled the pool evolution and the distribution on a regional scale. Owing to the multiple periods of the pool formation and the reconstruction, the distribution of marine gas pools was complex both in space and in time, and the gas in the pools is heterogeneous. Pool elements, such as preservation conditions, traps and migration paths, and reservoir rocks and facies, also served as important control factors to marine gas pools in South China. Especially, the preservation conditions played a key role in maintaining marine oil and gas accumulations on a regional or local scale. According to several dominant control factors of a pool, the pool-controlling model can be constructed. As an example, the pool-controlling model of Sinian gas pool in Weiyuan gas field in Sichuan basin was summed up. ?? Higher Education Press and Springer-Verlag 2007.

Frontiers of Earth Science in China↗

Federated or cached searches: providing expected performance from multiple invasive species databases

Invasive species are a universal global problem, but the information to identify them, manage them, and prevent invasions is stored around the globe in a variety of formats. The Global Invasive Species Information Network is a consortium of organizations working toward providing seamless access to these disparate databases via the Internet. A distributed network of databases can be created using the Internet and a standard web service protocol. There are two options to provide this integration. First, federated searches are being proposed to allow users to search “deep” web documents such as databases for invasive species. A second method is to create a cache of data from the databases for searching. We compare these two methods, and show that federated searches will not provide the performance and flexibility required from users and a central cache of the datum are required to improve performance.

Frontiers of Earth Science↗

A habitat overlap analysis derived from maxent for tamarisk and the south-western willow flycatcher

Biologic control of the introduced and invasive, woody plant tamarisk ( Tamarix spp, saltcedar) in south-western states is controversial because it affects habitat of the federally endangered South-western Willow Flycatcher ( Empidonax traillii extimus ). These songbirds sometimes nest in tamarisk where floodplain-level invasion replaces native habitats. Biologic control, with the saltcedar leaf beetle ( Diorhabda elongate ), began along the Virgin River, Utah, in 2006, enhancing the need for comprehensive understanding of the tamarisk-flycatcher relationship. We used maximum entropy (Maxent) modeling to separately quantify the current extent of dense tamarisk habitat (>50% cover) and the potential extent of habitat available for E. traillii extimus within the studied watersheds. We used transformations of 2008 Landsat Thematic Mapper images and a digital elevation model as environmental input variables. Maxent models performed well for the flycatcher and tamarisk with Area Under the ROC Curve (AUC) values of 0.960 and 0.982, respectively. Classification of thresholds and comparison of the two Maxent outputs indicated moderate spatial overlap between predicted suitable habitat for E. traillii extimus and predicted locations with dense tamarisk stands, where flycatcher habitat will potentially change flycatcher habitats. Dense tamarisk habitat comprised 500 km 2 within the study area, of which 11.4% was also modeled as potential habitat for E. traillii extimus . Potential habitat modeled for the flycatcher constituted 190 km 2 , of which 30.7% also contained dense tamarisk habitat. Results showed that both native vegetation and dense tamarisk habitats exist in the study area and that most tamarisk infestations do not contain characteristics that satisfy the habitat requirements of E. traillii extimus . Based on this study, effective biologic control of Tamarix spp. may, in the short term, reduce suitable habitat available to E. traillii extimus , but also has the potential in the long term to increase suitable habitat if appropriate mixes of native woody vegetation replace tamarisk in biocontrol areas.

Frontiers of Earth Science↗

Habitat suitability of patch types: a case study of the Yosemite toad

Understanding patch variability is crucial in understanding the spatial population structure of wildlife species, especially for rare or threatened species. We used a well-tested maximum entropy species distribution model (Maxent) to map the Yosemite toad (Anaxyrus (= Bufo) canorus) in the Sierra Nevada mountains of California. Twenty-six environmental variables were included in the model representing climate, topography, land cover type, and disturbance factors (e.g., distances to agricultural lands, fire perimeters, and timber harvest areas) throughout the historic range of the toad. We then took a novel approach to the study of spatially structured populations by applying the species-environmental matching model separately for 49 consistently occupied sites of the Yosemite toad compared to 27 intermittently occupied sites. We found that the distribution of the entire population was highly predictable (AUC = 0.95±0.03 SD), and associated with low slopes, specific vegetation types (wet meadow, alpine-dwarf shrub, montane chaparral, red fir, and subalpine conifer), and warm temperatures. The consistently occupied sites were also associated with these same factors, and they were also highly predictable (AUC = 0.95±0.05 SD). However, the intermittently occupied sites were associated with distance to fire perimeter, a slightly different response to vegetation types, distance to timber harvests, and a much broader set of aspect classes (AUC = 0.90±0.11 SD). We conclude that many studies of species distributions may benefit by modeling spatially structured populations separately. Modeling and monitoring consistently-occupied sites may provide a realistic snapshot of current species-environment relationships, important climatic and topographic patterns associated with species persistence patterns, and an understanding of the plasticity of the species to respond to varying climate regimes across its range. Meanwhile, modeling and monitoring of widely dispersing individuals and intermittently occupied sites may uncover environmental thresholds and human-related threats to population persistence.

Frontiers of Earth Science↗

Moisture and temperature controls on nitrification differ among ammonia oxidizer communities from three alpine soil habitats

Climate change is altering the timing and magnitude of biogeochemical fluxes in many high elevation ecosystems. The consequent changes in alpine nitrification rates have the potential to influence ecosystem scale responses. In order to better understand how changing temperature and moisture conditions may influence ammonia oxidizers and nitrification activity, we conducted laboratory incubations on soils collected in a Colorado watershed from three alpine habitats (glacial outwash, talus, and meadow). We found that bacteria, not archaea, dominated all ammonia oxidizer communities. Nitrification increased with moisture in all soils and under all temperature treatments. However, temperature was not correlated with nitrification rates in all soils. Site-specific temperature trends suggest the development of generalist ammonia oxidizer communities in soils with greater in situ temperature fluctuations and specialists in soils with more steady temperature regimes. Rapidly increasing temperatures and changing soil moisture conditions could explain recent observations of increased nitrate production in some alpine soils.

Colorado↗

Pulverization provides a mechanism for the nucleation of earthquakes at low stress on strong faults

An earthquake occurs when rock that has been deformed under stress rebounds elastically along a fault plane ( Gilbert, 1884 ; Reid, 1911 ), radiating seismic waves through the surrounding earth. Rupture along the entire fault surface does not spontaneously occur at the same time, however. Rather the rupture starts in one tiny area, the rupture nucleation zone, and spreads sequentially along the fault. Like a row of dominoes, one bit of rebounding fault triggers the next. This triggering is understood to occur because of the large dynamic stresses at the tip of an active seismic rupture. The importance of these crack tip stresses is a central question in earthquake physics. The crack tip stresses are minimally important, for example, in the time predictable earthquake model ( Shimazaki and Nakata, 1980 ), which holds that prior to rupture stresses are comparable to fault strength in many locations on the future rupture plane, with bits of variation. The stress/strength ratio is highest at some point, which is where the earthquake nucleates. This model does not require any special conditions or processes at the nucleation site; the whole fault is essentially ready for rupture at the same time. The fault tip stresses ensure that the rupture occurs as a single rapid earthquake, but the fact that fault tip stresses are high is not particularly relevant since the stress at most points does not need to be raised by much. Under this model it should technically be possible to forecast earthquakes based on the stress-renewaql concept, or estimates of when the fault as a whole will reach the critical stress level, a practice used in official hazard mapping ( Field, 2008 ). This model also indicates that physical precursors may be present and detectable, since stresses are unusually high over a significant area before a large earthquake.

Frontiers in Earth Science↗

Pan-arctic trends in terrestrial dissolved organic matter from optical measurements

Climate change is causing extensive warming across Arctic regions resulting in permafrost degradation, alterations to regional hydrology and shifting amounts and composition of dissolved organic matter (DOM) transported by streams and rivers. Here, we characterize the DOM composition and optical properties of the six largest Arctic rivers draining into the Arctic Ocean to examine the ability of optical measurements to provide meaningful insights into terrigenous carbon export patterns and biogeochemical cycling. The chemical composition of aquatic DOM varied with season, spring months were typified by highest lignin phenol and dissolved organic carbon (DOC) concentrations with greater hydrophobic acid content, and lower proportions of hydrophilic compounds, relative to summer and winter months. Chromophoric DOM (CDOM) spectral slope ( S 275–295 ) tracked seasonal shifts in DOM composition across river basins. Fluorescence and parallel factor analysis identified seven components across the six Arctic rivers. The ratios of “terrestrial humic-like” vs. “marine humic-like” fluorescent components co-varied with lignin monomer ratios over summer and winter months, suggesting fluorescence may provide information on the age and degradation state of riverine DOM. CDOM absorbance ( a 350 ) proved a sensitive proxy for lignin phenol concentrations across all six river basins and over the hydrograph, enabling for the first time the development of a single pan-arctic relationship between a 350 and terrigenous DOC ( R 2 = 0.93). Combining this lignin proxy with high-resolution monitoring of a 350 , pan-arctic estimates of annual lignin flux were calculated to range from 156 to 185 Gg, resulting in shorter and more constrained estimates of terrigenous DOM residence times in the Arctic Ocean (spanning 7 months to 2½ years). Furthermore, multiple linear regression models incorporating both absorbance and fluorescence variables proved capable of explaining much of the variability in lignin composition across rivers and seasons. Our findings suggest that synoptic, high-resolution optical measurements can provide improved understanding of northern high-latitude organic matter cycling and flux, and prove an important technique for capturing future climate-driven changes.

Frontiers in Earth Science↗

What do we really know about the role of microorganisms in iron sulfide mineral formation?

Iron sulfide mineralization in low-temperature systems is a result of biotic and abiotic processes, though the delineation between these two modes of formation is not always straightforward. Here we review the role of microorganisms in the precipitation of extracellular iron sulfide minerals. We summarize the evidence that links sulfur-metabolizing microorganisms and sulfide minerals in nature and we present a critical overview of laboratory-based studies of the nucleation and growth of iron sulfide minerals in microbial cultures. We discuss whether biologically derived minerals are distinguishable from abiotic minerals, possessing attributes that are uniquely diagnostic of biomineralization. These inquiries have revealed the need for additional thorough, mechanistic and high-resolution studies to understand microbially mediated formation of a variety of sulfide minerals across a range of natural environments.

Frontiers in Earth Science↗

An empirical approach for estimating stress-coupling lengths for marine-terminating glaciers

Variability in the dynamic behavior of marine-terminating glaciers is poorly understood, despite an increase in the abundance and resolution of observations. When paired with ice thicknesses, surface velocities can be used to quantify the dynamic redistribution of stresses in response to environmental perturbations through computation of the glacier force balance. However, because the force balance is not purely local, force balance calculations must be performed at the spatial scale over which stresses are transferred within glacier ice, or the stress-coupling length (SCL). Here we present a new empirical method to estimate the SCL for marine-terminating glaciers using high-resolution observations. We use the empirically-determined periodicity in resistive stress oscillations as a proxy for the SCL. Application of our empirical method to two well-studied tidewater glaciers (Helheim Glacier, SE Greenland, and Columbia Glacier, Alaska, USA) demonstrates that SCL estimates obtained using this approach are consistent with theory (i.e., can be parameterized as a function of the ice thickness) and with prior, independent SCL estimates. In order to accurately resolve stress variations, we suggest that similar empirical stress-coupling parameterizations be employed in future analyses of glacier dynamics.

Frontiers in Earth Science↗

The relative effectiveness of empirical and physical models for simulating the dense undercurrent of pyroclastic flows under different emplacement conditions

High concentration pyroclastic density currents (PDCs) are hot avalanches of volcanic rock and gas and are among the most destructive volcanic hazards due to their speed and mobility. Mitigating the risk associated with these flows depends upon accurate forecasting of possible impacted areas, often using empirical or physical models. TITAN2D, VolcFlow, LAHARZ, and Δ H/L or energy cone models each employ different rheologies or empirical relationships and therefore differ in appropriateness of application for different types of mass flows and topographic environments. This work seeks to test different statistically- and physically-based models against a range of PDCs of different volumes, emplaced under different conditions, over different topography in order to test the relative effectiveness, operational aspects, and ultimately, the utility of each model for use in hazard assessments. The purpose of this work is not to rank models, but rather to understand the extent to which the different modeling approaches can replicate reality in certain conditions, and to explore the dynamics of PDCs themselves. In this work, these models are used to recreate the inundation areas of the dense-basal undercurrent of all 13 mapped, land-confined, Soufrière Hills Volcano dome-collapse PDCs emplaced from 1996 to 2010 to test the relative effectiveness of different computational models. Best-fit model results and their input parameters are compared with results using observation- and deposit-derived input parameters. Additional comparison is made between best-fit model results and those using empirically-derived input parameters from the FlowDat global database, which represent “forward” modeling simulations as would be completed for hazard assessment purposes. Results indicate that TITAN2D is able to reproduce inundated areas well using flux sources, although velocities are often unrealistically high. VolcFlow is also able to replicate flow runout well, but does not capture the lateral spreading in distal regions of larger-volume flows. Both models are better at reproducing the inundated area of single-pulse, valley-confined, smaller-volume flows than sustained, highly unsteady, larger-volume flows, which are often partially unchannelized. The simple rheological models of TITAN2D and VolcFlow are not able to recreate all features of these more complex flows. LAHARZ is fast to run and can give a rough approximation of inundation, but may not be appropriate for all PDCs and the designation of starting locations is difficult. The Δ H/L cone model is also very quick to run and gives reasonable approximations of runout distance, but does not inherently model flow channelization or directionality and thus unrealistically covers all interfluves. Empirically-based models like LAHARZ and Δ H/L cones can be quick, first-approximations of flow runout, provided a database of similar flows, e.g., FlowDat, is available to properly calculate coefficients or Δ H/L . For hazard assessment purposes, geophysical models like TITAN2D and VolcFlow can be useful for producing both scenario-based or probabilistic hazard maps, but must be run many times with varying input parameters. LAHARZ and Δ H/L cones can be used to produce simple modeling-based hazard maps when run with a variety of input volumes, but do not explicitly consider the probability of occurrence of different volumes. For forward modeling purposes, the ability to derive potential input parameters from global or local databases is crucial, though important input parameters for VolcFlow cannot be empirically estimated. Not only does this work provide a useful comparison of the operational aspects and behavior of various models for hazard assessment, but it also enriches conceptual understanding of the dynamics of the PDCs themselves.

Frontiers in Earth Science↗

The Polar WRF downscaled historical and projected twenty-first century climate for the coast and foothills of Arctic Alaska

Climate change is most pronounced in the northern high latitude region. Yet, climate observations are unable to fully capture regional-scale dynamics due to the sparse weather station coverage, which limits our ability to make reliable climate-based assessments. A set of simulated data products was therefore developed for the North Slope of Alaska through a dynamical downscaling approach. The polar-optimized Weather Research and Forecast (Polar WRF) model was forced by three sources: The ERA-interim reanalysis data (for 1979–2014), the Community Earth System Model 1.0 (CESM1.0) historical simulation (for 1950–2005), and the CESM1.0 projected (for 2006–2100) simulations in two Representative Concentration Pathways (RCP4.5 and RCP8.5) scenarios. Climatic variables were produced in a 10-km grid spacing and a 3-h interval. The ERA-interim forced WRF (ERA-WRF) proves the value of dynamical downscaling, which yields more realistic topographical-induced precipitation and air temperature, as well as corrects underestimations in observed precipitation. In summary, dry and cold biases to the north of the Brooks Range are presented in ERA-WRF, while CESM forced WRF (CESM-WRF) holds wet and warm biases in its historical period. A linear scaling method allowed for an adjustment of the biases, while keeping the majority of the variability and extreme values of modeled precipitation and air temperature. CESM-WRF under RCP 4.5 scenario projects smaller increase in precipitation and air temperature than observed in the historical CESM-WRF product, while the CESM-WRF under RCP 8.5 scenario shows larger changes. The fine spatial and temporal resolution, long temporal coverage, and multi-scenario projections jointly make the dataset appropriate to address a myriad of physical and biological changes occurring on the North Slope of Alaska.

Alaska↗

Opal-A in glassy pumice, acid alteration, and the 1817 phreatomagmatic eruption at Kawah Ijen (Java), Indonesia

At Kawah Ijen (Indonesia), vigorous SO 2 and HCl degassing sustains a hyperacid lake (pH ~0) and intensely alters the subsurface, producing widespread residual silica and advanced argillic alteration products. In 1817, a VEI 2 phreatomagmatic eruption evacuated the lake, depositing a widespread layer of muddy ash fall, and sending lahars down river drainages. We discovered multiple types of opaline silica in juvenile low-silica dacite pumice and in particles within co-erupted laharic sediments. Most spectacular are opal-replaced phenocrysts of plagioclase and pyroxene adjacent to pristine matrix glass and melt inclusions. Opal-bearing pumice has been found at numerous sites, including where post-eruption infiltration of acid water is unlikely. Through detailed analyses of an initial sampling of 1817 eruption products, we find evidence for multiple origins of opaline materials in pumice and laharic sediments. Evidently, magma encountered acid-altered materials in the subsurface and triggered phreatomagmatic eruptions. Syn-eruptive incorporation of opal-alunite clasts, layered opal, and fragment-filled vesicles of opal and glass, all suggest magma-rock interactions in concert with vesiculation, followed by cooling within minutes. Our experiments at magmatic temperature confirm that the opaline materials would show noticeable degradation in time periods longer than a few tens of minutes. Some glassy laharic sedimentary grains are more andesitic than the main pumice type and may represent older volcanic materials that were altered beneath the lake bottom and were forcefully ejected during the 1817 eruption. A post-eruptive origin remains likely for most of the opal-replaced phenocrysts in pumice. Experiments at 25°C and 100°C reveal that when fresh pumice is bathed in Kawah Ijen hyperacid fluid for 6 weeks, plagioclase is replaced without altering either matrix glass or melt inclusions. Moreover, lack of evidence for high-temperature annealing of the opal suggests that post-eruption alteration of pumice is more likely than pre-eruption envelopment of euhedral opal-replaced phenocrysts in dacitic melt. At Ijen and elsewhere, the ascent of magma into hydrous acid-altered mineral assemblages (e.g., opal, kaolinite, alunite) could induce rapid dehydration of hydrous minerals and amorphous materials, generating considerable steam and contributing to magmatic-hydrothermal and phreatomagmatic explosions.

Kawah Ijen↗

A snow density dataset for improving surface boundary conditions in Greenland ice sheet firn modeling

The surface snow density of glaciers and ice sheets is of fundamental importance in converting volume to mass in both altimetry and surface mass balance studies, yet it is often poorly constrained. Site-specific surface snow densities are typically derived from empirical relations based on temperature and wind speed. These parameterizations commonly calculate the average density of the top meter of snow, thereby systematically overestimating snow density at the actual surface. Therefore, constraining surface snow density to the top 0.1 m can improve boundary conditions in high-resolution firn-evolution modeling. We have compiled an extensive dataset of 200 point measurements of surface snow density from firn cores and snow pits on the Greenland ice sheet. We find that surface snow density within 0.1 m of the surface has an average value of 315 kg m −3 with a standard deviation of 44 kg m −3 , and has an insignificant annual air temperature dependency. We demonstrate that two widely-used surface snow density parameterizations dependent on temperature systematically overestimate surface snow density over the Greenland ice sheet by 17–19%, and that using a constant density of 315 kg m −3 may give superior results when applied in surface mass budget modeling.

Frontiers in Earth Science↗

Analysis of the Alaska Volcano Observatory’s response time to volcanic explosions-1989 to 2016

A major goal of volcano monitoring is the rapid identification of volcanic explosions and subsequent warning of associated hazards. Between 1988 and 2016 the Alaska Volcano Observatory (AVO) responded to at least 54 separate volcanic eruptions. During this period, AVO's monitoring program relied principally on seismic and satellite remote sensing data, supplemented with geodetic, gas, and visual observations to track volcanic unrest. In this study we focus on AVO's response time, or the time required for AVO to (1) identify seismic signals associated with large ash-producing volcanic explosions and (2) initiate public warnings. We restrict this analysis to volcanoes monitored by a local seismic network and explosive in character. We focus on the 1989–90 eruption of Redoubt Volcano (VEI 3), the 1992 eruption of Mount Spurr (VEI 4), the 1999 eruption of Shishaldin Volcano (VEI 3), the 2006 eruption of Augustine Volcano (VEI 3) and the 2016 eruption of Pavlof Volcano (VEI 2) as detailed records of the timing of formal warnings are preserved. These eruption sequences allow us to evaluate AVO's response time under a number of monitoring scenarios, including both expected (those with recognized precursory unrest) and surprise eruptions (those without identified precursory unrest) as well as individual and repetitive sequences of explosive events. Recorded response time ranges from ~1 to 86 min. The shorter response times (~1–13 min) were achieved during sequences of explosive events at Redoubt (1989–90), Spurr (1992) and Augustine (2006). The longer response times (31– 86 min) are recorded for unexpected or surprise explosions such as Spurr (August 18, 1992) and Pavlof (2016) and the only or first explosions in an eruptive sequence such as Shishaldin (1999) and Augustine (2006).

Alaska↗

Alaska Volcano Observatory alert and forecasting timeliness: 1989–2017

The Alaska Volcano Observatory (AVO) monitors volcanoes in Alaska and issues notifications and warnings of volcanic unrest and eruption. We evaluate the timeliness and accuracy of eruption forecasts for 53 eruptions at 20 volcanoes, beginning with Mount Redoubt's 1989–1990 eruption. Successful forecasts are defined as those where AVO issued a formal warning before eruption onset. These warning notifications are now part of AVO's Aviation Color Code and Volcanic Alert Level. This analysis considers only the start of an eruption, although many eruptions have multiple phases of activity. For the 21 eruptions at volcanoes with functioning local seismic networks, AVO has high forecasting success at volcanoes with: >15 years repose intervals and magmatic eruptions (4 out of 4, 100%); or larger eruptions (Volcanic Explosivity Index (VEI) 3 or greater; 6 out of 10, 60%). Therefore, AVO successfully forecast all four monitored, longer-repose period, VEI 3+ eruptions: Redoubt 1989–1990 and 2009, Spurr 1992, and Augustine 2005–2006. For volcanoes with functioning seismic monitoring networks, success rates are lower for: volcanoes with shorter repose periods (3 out of 16, 19%); more mafic compositions (3 out of 18, 17%); or smaller eruption size (VEI 2 or less, 1 out of 11, 9%). These eruptions (Okmok, Pavlof, Veniaminof, and Shishaldin) often lack detectable precursory signals. For 32 eruptions at volcanoes without functioning local seismic networks, the forecasting success rate is much lower (2, 6%; Kasatochi 2008 and Shishaldin 2014). For remote volcanoes where the main hazard is to aviation, rapid detection is a goal in the absence of in situ monitoring. Eruption detection has improved in recent years, shown by a decrease in the time between eruption onset and notification. Even limited seismic monitoring can detect precursory activity at volcanoes with certain characteristics (intermediate composition, longer repose times, larger eruptions), but difficulty persists in detecting subtle precursory activity at frequently active volcanoes with more mafic compositions. This suggests that volcano-specific characteristics should be considered when designing monitoring programs and evaluating forecasting success. More proximally-located sensors and data types are likely needed to forecast eruptive activity at frequently-active, more mafic volcanoes that generally produce smaller eruptions.

Alaska↗