USGS Science⌕ Search

SEARCH · USGS Science

Results for “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 739 records · Page 41Linked to original sources

Born of fire: In search of volcanoes in U.S. national parks, four striking examples

Geologic features, particularly volcanic features, have been protected by the National Park Service since its inception. Some volcanic areas were nationally protected even before the National Park Service was established. The first national park, Yellowstone National Park, is one of the most widely known geothermal and volcanic areas in the world. It contains the largest volcanic complex in North America and has experienced three eruptions which rate among the largest eruptions known to have occurred on Earth. Half of the twelve areas established as national parks before the 1916 Organic Act which created the National Park Service are centered on volcanic features. The National Park Service now manages lands that contain nearly every conceivable volcanic resource, with at least seventy-six managed lands that contain volcanoes or volcanic rocks. Given that so many lands managed by the National Park Service contain volcanoes and volcanic rocks, we cannot give an overview of the history of each one; rather we highlight four notable examples of parks that were established on account of their volcanic landscapes. These parks all helped to encourage the creation and success of the National Park Service by inspiring the imagination of the public. In addition to preserving and providing access to the nation's volcanic heritage, volcanic national parks are magnificent places to study and understand volcanoes and volcanic landscapes in general. Scientists from around the world study volcanic hazards, volcanic history, and the inner working of the Earth within U.S. national parks. Volcanic landscapes and associated biomes that have been relatively unchanged by human and economic activities provide unique natural laboratories for understanding how volcanoes work, how we might predict eruptions and hazards, and how these volcanoes affect surrounding watersheds, flora, fauna, atmosphere, and populated areas.

Earth Sciences History↗

The John Wesley Powell Center for Analysis and Synthesis

The Powell Center provides an environment for cross-disciplinary scientific collaboration. The Center expands U.S. Geological Survey earth system science synthesis research activities by fostering the innovation that results from accumulated knowledge, constructive errors, and the "information spillover" that emerges from collaborative settings. Working Groups at the Powell Center use existing data to produce new knowledge..

Fact Sheet↗

An interdisciplinary assessment of regional-scale nonpoint source ground-water vulnerability: Theory and application

An integrated earth science-economics model, developed within a geographic information system (GIS), combines a regional-scale nonpoint source vulnerability assessment with a specific remediation measure to avoid unnecessary agricultural production costs associated with the use of agrochemicals in the Pearl Harbor basin on the island of Oahu, Hawaii. This approach forms the core of a risk-based regulation for the application of agrochemicals and estimates the benefits of an information-based approach to decisionmaking.

Hawaii↗

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↗

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↗

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↗

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↗

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↗

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↗