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

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↗
USGS2011

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↗
USGS2011

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↗
USGS2011

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↗
USGS2015

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↗
USGS2014

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↗
USGS2016

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↗
USGS2016

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↗
USGS2017

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↗
USGS2018

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↗
USGS2018

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↗
USGS2018

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↗
USGS2018

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↗
USGS2018

Prevalence of seismic rate anomalies preceding volcanic eruptions in Alaska

Seismic rate increases often precede eruptions at volcanoes worldwide. However, many eruptions occur without such precursors. Additionally, identifying seismic rate increases near volcanoes with high levels of background seismicity is non-trivial and many periods of elevated seismicity occur without ensuing eruptions, limiting their usefulness for forecasting in some cases. Although these issues are commonly known, efforts to quantify them are limited. In this study, we consistently apply a common statistical tool, the β-statistic, to seismically monitored eruptions in Alaska of various styles to determine the overall prevalence of seismic rate anomalies immediately preceding eruptions. We find that 6 out of 20 (30%) eruptions have statistically significant precursory seismic rate increases. Of these 6 eruptions, 3 of them occur at volcanoes with relatively felsic compositions, repose periods >15 years, and VEI ≥ 3. Overall, our results confirm that seismic rate increases are common prior to larger eruptions at long dormant, “closed-system” volcanoes, but uncommon preceding smaller eruptions at more frequently active, “open-system” volcanoes with more mafic magmas. We also explore the rate of other anomalies not precursory to eruptions and investigate their origins. Some of these non-eruptive anomalies can be explained by aftershocks of regional seismic events, magmatic activity that did not lead to eruption, or unrest at other nearby volcanoes. Some open-system volcanoes have high non-eruptive anomaly rates and low pre-eruptive anomaly rates and are thus not amenable to forecasting based on earthquake catalogs. In this study, we find that 31% of anomalies lead to eruption. With continued calibration at more volcanoes, the β-statistic that we apply may be used more broadly to analyze future periods of seismic unrest at other volcanoes, properly placing such episodes into the context of the long-term background rate. These results may be useful for informing future eruption forecasts around the world, and the statistical tool may aid volcano observatories in identifying future seismic rate anomalies under changing network conditions.

Alaska↗
USGS2018

Short-term forecasting and detection of explosions during the 2016–2017 eruption of Bogoslof volcano, Alaska

We describe a multidisciplinary approach to forecast, rapidly detect, and characterize explosive events during the 2016–2017 eruption of Bogoslof volcano, a back-arc shallow submarine volcano in Alaska’s Aleutian arc. The eruptive sequence began in December 2016 and included about 70 discrete explosive events. Because the volcano has no local monitoring stations, we used distant stations on the nearest volcanoes, Okmok (54 km) and Makushin (72 km), combined with regional infrasound sensors and lightning detection from the Worldwide Lightning Location Network (WWLLN). Pre-eruptive seismicity was detected for 12 events during the first half of the eruption; for all other events co-eruptive signals allowed for detection only. Monitoring of activity used a combination of scheduled checks combined with automated alarms. Alarms triggered on real-time data included real-time seismic amplitude measurement (RSAM); infrasound from several arrays, the closest being on Okmok; and lightning strokes detected from WWLLN within a 20-km radius of the volcano. During periods of unrest, a multidisciplinary response team of four people fulfilled specific roles to evaluate geophysical and remote-sensing data, run event-specific ash-cloud dispersion models, ensure interagency coordination, and develop and distribute of formalized warning products. Using this approach, for events that produced ash clouds ≥7.5 km above sea level, Alaska Volcano Observatory (AVO) called emergency response partners 15 min, and issued written notices 30 min, after event onset (mean times). Factors that affect timeliness of written warnings include event size and number of data streams available; bigger events and more data both decrease uncertainty and allow for faster warnings. In remote areas where airborne ash is the primary hazard, the approach used at Bogoslof is an effective strategy for hazard mitigation.

Alaska↗
USGS2018

Mount St. Helens retrospective: Lessons learned since 1980 and remaining challenges

Since awakening from a 123-year repose in 1980, Mount St. Helens has provided an opportunity to study changes in crustal magma storage at an active arc volcano—a process of fundamental importance to eruption forecasting and hazards mitigation. There has been considerable progress, but important questions remain unanswered. Was the 1980 eruption triggered by an injection of magma into an upper crustal reservoir? If so, when? How did magma rise into the edifice without producing detectable seismicity deeper than ∼2.5 km or measurable surface deformation beyond the volcano’s north flank? Would precursory activity have been recognized earlier if current monitoring techniques had been available? Despite substantial improvements in monitoring capability, similar questions remain after the dome-forming eruption of 2004–2008. Did additional magma accumulate in the reservoir between the end of the 1980–1986 eruption and the start of the 2004–2008 eruption? If so, when? What is the significance of a relative lull in seismicity and surface deformation for several years prior to the 2004–2008 eruption onset? How did magma reach the surface without producing seismicity deeper than ∼2 km or measurable deformation more than a few hundred meters from the vent? Has the reservoir been replenished since the eruption ended, and is it now primed for the next eruption? What additional precursors, if any, should be expected? This paper addresses these questions, explores possible answers, and identifies unresolved issues in need of additional study. The 1980–1986 and 2004–2008 eruptions could have resulted from second boiling during crystallization of magma long-resident in an upper crustal reservoir, rather than from injection of fresh magma from below. If reservoir pressurization and magma ascent were slow enough, resulting strain might have been accommodated by viscoelastic deformation, without appreciable seismicity or surface deformation, until rising magma entered a brittle regime within 2–2.5 km of the surface. Given the remarkably gas-poor nature of the 2004–2008 dome lava, future eruptive activity might require a relatively long period of quiescence and reservoir pressurization or a large injection of fresh magma—an event that arguably has not occurred since the Kalama eruptive period (C.E. 1479–1720).

Washington↗
USGS2018

Measuring SO2 emission rates at Kīlauea Volcano, Hawaii, using an array of upward-looking UV spectrometers, 2014-2017

Retrieving accurate volcanic sulfur dioxide (SO 2 ) gas emission rates is important for a variety of purposes. It is an indicator of shallow subsurface magma, and thus may signal impending eruption or unrest. SO 2 emission rates are significant for accurately assessing climate impact, and providing context for assessing environmental, agricultural, and human health effects during volcanic eruptions. The U.S. Geological Survey Hawaiian Volcano Observatory uses an array of ten fixed, upward-looking ultraviolet spectrometer systems to measure SO 2 emission rates at 10-s sample intervals from the Kīlauea summit. We present Kīlauea SO 2 emission rates from the volcano’s summit and middle East Rift Zone during 2014–2017 and discuss the major sources of error for these measurements. Due to the wide range of SO 2 emissions encountered at the summit vent, we used a variable wavelength spectral analysis range to accurately quantify both high and low SO 2 column densities. We compare measured emission rates from the fixed spectrometer array to independent road and helicopter-based traverse measurements and evaluate the magnitudes and sources of uncertainties for each method. To address the challenge of obtaining accurate plume speed measurements, we examine ground-based wind-speed, plume speed tracking via spectrometer, and SO 2 camera derived plume speeds. Our analysis shows that: (1) the summit array column densities calculated using a dual fit window, are within -6 to +22% of results obtained with a variety of other conventional and experimental retrieval methods; (2) emission rates calculated from the summit array located ∼3 km downwind provide the best, practical estimate of summit SO 2 release under normal trade wind conditions; (3) ground-based anemometer wind speeds are 22% less than plume speeds determined by cross-correlation of plume features; (4) our best estimate of average Kīlauea SO 2 release for 2014–2017 is 5100 t/d, which is comparable to the space-based OMI emissions of 5518 t/d; and (5) short-term variability of SO 2 emissions reflects Kīlauea lava lake dynamics.

Hawaii↗
USGS2018

Best practices for elevation-based assessments of sea-level rise and coastal flooding exposure

Elevation data are critical for assessments of sea-level rise (SLR) and coastal flooding exposure. Previous research has demonstrated that the quality of data used in elevation-based assessments must be well understood and applied to properly model potential impacts. The cumulative vertical uncertainty of the input elevation data substantially controls the minimum increments of SLR and the minimum planning horizons that can be effectively used in assessments. For regional, continental, or global assessments, several digital elevation models (DEMs) are available for the required topographic information to project potential impacts of increased coastal water levels, whether a simple inundation model is used or a more complex process-based or probabilistic model is employed. When properly characterized, the vertical accuracy of the DEM can be used to report assessment results with the uncertainty stated in terms of a specific confidence level or likelihood category. An accuracy evaluation has been conducted of global DEMs to quantify their inherent vertical uncertainty to demonstrate how accuracy information should be considered when planning and implementing a SLR or coastal flooding assessment. The evaluation approach includes comparison of the DEMs with high-accuracy geodetic control points as the independent reference data over a variety of coastal relief settings. The global DEMs evaluated include SRTM, ASTER GDEM, ALOS World 3D, TanDEM-X, NASADEM, and MERIT. High-resolution, high-accuracy DEM sources, such as airborne lidar and stereo imagery, are also included to give context to the results from the global DEMs. The accuracy characterization results show that current global DEMs are not adequate for high confidence mapping of exposure to fine increments (<1 m) of SLR or with shorter planning horizons (<100 years) and thus they should not be used for such mapping, but they are suitable for general delineation of low elevation coastal zones. In addition to the best practice of rigorous accounting for vertical uncertainty, other recommended procedures are presented for delineation of different types of impact areas (marine and groundwater inundation) and use of regional relative SLR scenarios. The requirement remains for a freely available, high-accuracy, high-resolution global elevation model that supports quantitative SLR and coastal inundation assessments at high confidence levels.

Frontiers in Earth Science↗