SEARCH · Search USGS

Results for “Sustainability”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,297 records · Page 72Linked to original sources

USGS2018

Metamodeling and mapping of nitrate flux in the unsaturated zone and groundwater, Wisconsin, USA

Nitrate contamination of groundwater in agricultural areas poses a major challenge to the sustainability of water resources. Aquifer vulnerability models are useful tools that can help resource managers identify areas of concern, but quantifying nitrogen (N) inputs in such models is challenging, especially at large spatial scales. We sought to improve regional nitrate (NO 3 − ) input functions by characterizing unsaturated zone NO 3 − transport to groundwater through use of surrogate, machine-learning metamodels of a process-based N flux model. The metamodels used boosted regression trees (BRTs) to relate mappable landscape variables to parameters and outputs of a previous “vertical flux method” (VFM) applied at sampled wells in the Fox, Wolf, and Peshtigo (FWP) river basins in northeastern Wisconsin. In this context, the metamodels upscaled the VFM results throughout the region, and the VFM parameters and outputs are the metamodel response variables. The study area encompassed the domain of a detailed numerical model that provided additional predictor variables, including groundwater recharge, to the metamodels. We used a statistical learning framework to test a range of model complexities to identify suitable hyperparameters of the six BRT metamodels corresponding to each response variable of interest: NO 3 − source concentration factor (which determines the local NO 3 − input concentration); unsaturated zone travel time; NO 3 − concentration at the water table in 1980, 2000, and 2020 (three separate metamodels); and NO 3 − “extinction depth”, the eventual steady state depth of the NO 3 − front. The final metamodels were trained to 129 wells within the active numerical flow model area, and considered 58 mappable predictor variables compiled in a geographic information system (GIS). These metamodels had training and cross-validation testing R 2 values of 0.52 – 0.86 and 0.22 – 0.38, respectively, and predictions were compiled as maps of the above response variables. Testing performance was reasonable, considering that we limited the metamodel predictor variables to mappable factors as opposed to using all available VFM input variables. Relationships between metamodel predictor variables and mapped outputs were generally consistent with expectations, e.g. with greater source concentrations and NO 3 − at the groundwater table in areas of intensive crop use and well drained soils. Shorter unsaturated zone travel times in poorly drained areas likely indicated preferential flow through clay soils, and a tendency for fine grained deposits to collocate with areas of shallower water table. Numerical estimates of groundwater recharge were important in the metamodels and may have been a proxy for N input and redox conditions in the northern FWP, which had shallow predicted NO 3 − extinction depth. The metamodel results provide proof-of-concept for regional characterization of unsaturated zone NO 3 − transport processes in a statistical framework based on readily mappable GIS input variables.

Wisconsin↗

USGS2018

Suspended-sediment concentrations and loads in the lower Mississippi and Atchafalaya rivers decreased by half between 1980 and 2015

The Weighted Regressions on Time, Discharge, and Season (WRTDS) model was used to derive estimates of suspended-sediment concentration (SSC) and suspended-sediment load (SSL), their dependence on discharge, and their trends with confidence intervals, for one site each on the lowermost Mississippi and Atchafalaya Rivers. The WRTDS model reduces uncertainty in SSCs related to variable streamflow conditions. Flow-normalized SSCs in each river were similar, and decreased from about 260 mg/L to 130 mg/L from 1980 through 2015; combined annual SSL in the two rivers decreased from about 200 Megatons per year (MT/y) to about 100 MT/y. Declines in SSC and SSL were more gradual from 2005 through 2015 and show signs of stabilizing. Our estimates of SSL in 2015 differ markedly from several recently published estimates of current and projected future Mississippi River SSLs, which were generally around 200 MT/y. However, these values came mostly from a different site upstream on the Mississippi River. The relationship between SSC and streamflow differed in an important way between the two rivers. SSC increased as streamflow increased for the entire range of observed streamflow in the Atchafalaya River. In the Mississippi River, SSC followed the same pattern during low and moderate streamflow but decreased at the highest streamflow that tended to occur between January and July. Since much of the water flowing in the Atchafalaya originates from the Mississippi River, the difference suggests a within-basin source of suspended sediment for the Atchafalaya River that is absent in the lower Mississippi River. These findings have important implications for the restoration of deltaic wetlands in coastal Louisiana. Accurate estimates of the SSL available in each river are crucial for understanding how effective diversions of river water into adjacent estuaries will be in sustaining these wetlands. Our study demonstrates that there might be far less sediment available than previously reported. Further, the difference in the relationship between SSC and streamflow in the two rivers is highly relevant to the ongoing discussion of coastal restoration strategies because the delta building that is occurring at the mouth of the Atchafalaya River is frequently used as a model of what could be expected with controlled diversions in the lower Mississippi River delta. The differences in the SSC behavior with changes in streamflow between the two rivers needs to be considered when results from the Atchafalaya River system are projected to those of the Mississippi River.

Louisiana↗

USGS2020

Isotopic and geochemical assessment of the sensitivity of groundwater resources of Guam, Mariana Islands, to intra- and inter-annual variations in hydroclimate

Assessing the sensitivity of groundwater systems to hydroclimate variability is critical to sustainable management of the water resources of Guam, US territory. We assess spatial and temporal variability of isotopic and geochemical compositions of vadose and phreatic groundwater sampled from cave drip sites and production wells, respectively, to better understand the vulnerability of the freshwater lens on Guam to variability in hydroclimate. We independently evaluate the existing conceptual model of the Northern Guam Lens Aquifer that is largely based on physical, as opposed to geochemical, observations. Sampling was conducted from 2008 to 2015, over which rainfall gradually increased. Major ion geochemistry and Sr isotope values of groundwater show varying influence from soil, limestone bedrock, and seawater. Geochemical modeling that can explain spatial variability in groundwater Na+ and Mg2+ concentrations and Sr/Ca and 87Sr/86 Sr values indicates that groundwater compositions are dominantly controlled by mixing of freshwater with seawater and water-rock interaction. Differences between amount-weighted annual average precipitation δ18 O values and groundwater δ18 O values indicate a recharge bias toward the wet season, consistent with other tropical carbonate island aquifer settings. Intra- and inter-annual variations in Na+ concentrations and δ18 O values in groundwater reflect sensitivity of recharge to seasonal variations in rainfall amount and changes in annual rainfall amounts. Our results indicate the influence of multiple modes of recharge on groundwater compositions and spatial variability in the sensitivity of groundwater to seawater mixing. This sensitivity of the freshwater lens points to the vulnerability of groundwater resources to changes in recharge associated with climate, land-use change, and increases in population.

Journal of Hydrology↗

USGS2020

Rethinking a groundwater flow system using a multiple-tracer geochemical approach: A case study in Moab-Spanish Valley, Utah

The Glen Canyon Group Aquifer (GCGA) is the sole source of public water supply for the city of Moab, Utah, a domestic and international tourist destination. Population and tourism growth are likely to target the GCGA for future water resources, but our analysis indicates that additional withdrawals would likely be sourced from groundwater storage and not be sustained by recharge. A quantitative estimate of groundwater discharge from the GCGA is problematic because the downgradient aquifer boundary is the Colorado River, and groundwater discharge to the river is very small compared to the river flow. A water budget based on a conceptual model of GCGA discharging into an adjacent alluvial Valley-Fill Aquifer (VFA) was reported by Sumsion (1971) and numerous subsequent studies have repeated and utilized this water budget. The GCGA contains stable isotopes, tritium, 3 He/ 4 He ratios, dissolved solids, and sulfate concentrations that contrast with the VFA, indicating it is instead recharged by local streams rather than from the GCGA. Water-budget calculations, based on: (1) measured spring discharge and streamflow gains, (2) horizontal gradients in VFA groundwater age, and (3) GCGA outcrop area vadose-zone pore waters are all less than previously thought. Using a lumped parameter model and 14 C groundwater ages, we estimate recharge to the deeper GCGA (DGCGA) to be 4.2 ± 2.3 × 10 6 m 3 /yr, which is approximately equal to the measured discharge from wells and springs.

Utah↗

USGS2024

Diel temperature signals track seasonal shifts in localized groundwater contributions to headwater streamflow generation at network scale

Groundwater contributions to streamflow sustain aquatic ecosystem resilience; streams without significant groundwater inputs often have well-coupled air and water temperatures that degrade cold-water habitat during warm low flow periods. Widespread uncertainty in stream-groundwater connectivity across space and time has created disparate predictions of energy and nutrient fluxes across headwater networks, hindering predictions of cold-water habitat resilience under climate change scenarios. Recently, annual paired air and water temperature signals have been harnessed to indicate stream water thermal sensitivity and the dominance of deep versus shallow groundwater influence, although the utility of diel air–water temperature signal metrics for hydrologic inference has remained unexplored. Here we analyzed two consecutive years of locally paired, air–water temperature data from 47 headwater stream sites in the Catskill Mountains, New York, USA, and discovered characteristic seasonal patterns in diel temperature signal sinusoid metrics (amplitude ratio, phase lag, and mean ratio) driven by shifts in streamflow generation mechanisms and stream network position. Hydrologic interpretations of observed patterns were supported by stream heat budget model scenarios and additional analysis of paired air–water temperature data from two streams in Shenandoah National Park, Virginia, USA, with well characterized stream-groundwater connectivity. We found that within smaller tributaries, streamflow generation transitions from runoff to groundwater dominance were driven by hillslope drying during seasonal periods of lower precipitation. This was evidenced by significant correlations (p < 0.01) between daily water:air temperature signal amplitudes (non-linear decreases of ∼ 50 %) and derived base-flow index at 22 of the 28 sites, indicating enhanced local groundwater influence on streamflow promotes decoupling of diel air–water temperature signals. Additionally, ratios between daily water:air temperature signal means were lower in tributaries (∼0.68) when compared to main-stem (∼0.8) sites, increasing linearly throughout the observational period. In conceptual stream heat budget models, groundwater inflow had minimal effects on daily phase lags (∼0.2 hr), but increases in fractional groundwater discharge (0–50 %) depressed daily amplitude (∼20 % to 50 %) and mean ratios (∼15 %), supporting the sensitivity of daily metrics to interpreted changes in seasonal groundwater contributions to streamflow. During observational periods (i.e., April through October 2021 and 2022), significant differences (p < 0.01) between tributary and main-stem air–water metrics occurred when base-flow contributions were highest (∼0.93 vs. ∼ 0.68), as sites lower in the network had daily temperature metrics dominated by stream channel thermal inertia, rather than local groundwater connectivity, showing enhanced air–water diel signal coupling during warmer, drier periods. Divergent air temperature coupling across the network was interpreted as being driven by distance from local groundwater source zones, additional lateral groundwater inflows do not contribute a meaningful fraction to channel discharge lower in the network. Given the growing footprint of stream temperature observations, diel air–water temperature signals can provide distributed metrics sensitive to upstream groundwater discharge. Consequently, these metrics can support ongoing efforts by resource managers and researchers seeking to forecast the resilience of cold-water habitat to climate warming and changing precipitation regimes in mountain headwater streams.

New York, Virginia↗

USGS2024

Groundwater-Surface water interactions research: Past trends and future directions

Interactions between groundwater and surface water sustain groundwater-dependent ecosystems and regulate river temperature and biogeochemical cycles, amongst many other processes. These interactions occur in freshwater environments including rivers, springs, lakes, and wetlands, and in coastal environments via tidal pumping, submarine groundwater discharge, and seawater intrusion. Here, we explore groundwater-surface water interactions research using bibliometric analyses of titles, abstracts, and keywords from 20,275 journal papers published between 1970 and 2023 extracted from Scopus. Analyses show that research into groundwater-surface water interactions is highly multi-disciplinary, with growing contributions from the social and biological sciences. The number of groundwater-surface water interactions papers is rapidly increasing with over 1200 papers published per year since 2020. Drawing on our data-driven approach and expert knowledge, we synthesise current research trends and identify critical future research directions. Despite the thousands of papers on groundwater-surface water interactions, important processes are still difficult to quantify or predict at meaningful spatial scales to inform water-resources management. We see benefits in future groundwater-surface water interactions research focusing on: (1) using new technologies including internet-of-things-based sensors, uncrewed vehicles, and remote-sensing approaches for data collection to inform groundwater-surface water interactions at large scales, (2) seeking approaches to upscale site-specific findings to better inform management, and (3) continuing the movement towards multi-disciplinary investigations to better inform the understanding of groundwater-surface water interactions and processes that will enable better management outcomes.

Journal of Hydrology↗

USGS2025

Ecohydrological response of a forested headwater catchment to a flash drought in the Southeastern U.S.

Flash droughts differ from traditionally defined droughts in their rapidity of intensification and often associated high vapor-pressure deficit. These droughts can lead to declines in streamflow and water table depth and induce water stress to vegetation at a greater rate than droughts that manifest over longer periods. However, little is known regarding the response of forested environments to flash drought because most studies of impacts have been conducted in agricultural settings. In this study we investigated water-use patterns of riparian trees using sap flow methods and examined the role of groundwater as a source of moisture over three periods that were delimited by antecedent soil moisture conditions. For a longer-term perspective we also examine monthly streamflow over the 35-year record. We observed that trees at only one monitoring plot showed a decrease in water use relative to evaporative demand during a flash drought. Total reverse sap flow (flow toward the roots rather than the canopy) greatly increased during the flash drought period, suggesting the likely occurrence of hydraulic redistribution to the excessively dry soils. Over the drought period groundwater became a more dominant source of moisture for sustaining forest water use. Monthly mean streamflow during the flash drought approached levels observed in past multiyear droughts. This is the first study, to our knowledge, to specifically investigate the response of multiple water budget components to flash drought in a humid forest. As more studies are conducted, a better understanding of the range of expected responses are likely to emerge.

Georgia↗

USGS2025

Multi-lithofacies alluvial characterization via airborne electromagnetic-borehole fusion using ordinary interval kriging and geologic constraints

Airborne electromagnetic (AEM) survey provides extensive spatial coverage and detailed resolution at the near surface and can be used to develop hydrogeological models. However, utilization of AEM data is not straightforward because AEM resistivity is an indirect measurement for inferring sediment types. This study develops an ordinary interval kriging (OIK) algorithm and a resistivity-to-multi-lithofacies (R2ML) data fusion workflow for multi-lithofacies alluvial characterization. OIK utilizes irregular interval data to construct three-dimensional (3D) resistivity fields from one-dimensional inverted AEM resistivity models. The R2ML workflow maps the resistivity field generated from OIK into a multi-facies lithological model, incorporating geologic constraints derived from well logs and geological observations. The numerical and real-world cases demonstrate that OIK is computationally efficient, accounts for 3D anisotropy, and minimizes the smoothing effect, thereby preserving resistivity contrasts and reducing interpolation uncertainty. The methodology is applied to lithologic characterization of the Mississippi River Valley alluvial aquifer (MRVA) in the Shellmound area, Mississippi, U.S. A frequency-domain AEM survey was conducted to support groundwater studies for the managed aquifer recharge (MAR) to the MRVA. The resulting lithological model, including four types of lithofacies—clay, very fine sand, fine-medium sands, and graveliferous sands, illustrates the geomorphological processes of the MRVA and implies potential MAR. The alignment between the lithological model and existing geological and hydrogeological investigations demonstrates that OIK and R2ML workflow effectively capture the subsurface architecture of the MRVA. The methods have broad applicability for characterizing alluvial aquifers through AEM-borehole data fusion, supporting sustainable groundwater management.

Mississippi↗

USGS2026

Groundwater spatial variability within an atoll island: Assessing shallow aquifer heterogeneity with geophysical and physicochemical measurements

This study examines the spatial variability of shallow groundwater on Dhigelaabadhoo Island using electromagnetic induction surveys, groundwater monitoring, and sediment analyses. The research reveals how variations in island morphology—such as differences in elevation, reef flat width, and sediment composition—affect the spatial distribution of groundwater lenses and the overall aquifer dynamics. Saltwater intrusion is especially pronounced in low elevated areas, with narrow reef flat plate and areas where higher hydraulic conductivity—driven by the presence of coarser sediments—is observed, whereas regions characterized by finer sediments, higher elevation, and wider reef flat plates tend to support more symmetric and less saline groundwater lenses. The geophysical investigations reveal that tidal oscillations alter groundwater movement by markedly changing water levels and conductivity, thereby underscoring the critical need to account for temporal dynamics in atoll coastal aquifer systems and the importance of integrating tidal dynamics into the aquifer zone. The findings highlight the significant role of intrinsic morphological and external hydrodynamic factors in shaping groundwater distribution on atoll islands, offering critical insights for sustainable freshwater resource management.

Dhigelaabadhoo Island↗

USGS2015

Towards a global terrestrial species monitoring program

Introduction: The Convention for Biological Diversity’s (CBD) Strategic Plan for Biodiversity 2011-2020 envisions that “By 2050, biodiversity is valued, conserved, restored and wisely used, maintaining ecosystem services, sustaining a healthy planet and delivering benefits essential for all people.” Although 193 parties have adopted these goals, there is little infrastructure in place to monitor global biodiversity trends. Recent international conservation policy requires such data to be up-to-date, reliable, comparable among sites, relevant, and understandable; as is becoming obvious from the work plan adopted by the Intergovernmental Panel for Biodiversity and Ecosystem Services (IPBES: www.ipbes.net/; http://tinyurl.com/ohdnknq). In order to meet the five strategic goals of the Strategic Plan for Biodiversity 2011-2020 and its 20 accompanying Aichi Targets for 2020 (www.cbd.int/sp/targets/), advances need to be made in coordinating large-scale biodiversity monitoring and linking these with environmental data to develop a comprehensive Global Observation Network, as is the main idea behind GEOSS the Global Earth Observation System of Systems (Christian 2005)...Here we identify ten requirements important for the successful implementation of a global biodiversity monitoring network under the flag of GEO BON and especially a global terrestrial species monitoring program.

Journal for Nature Conservation↗

USGS2022

Invasive plant hitchhikers: Appalachian Trail thru-hiker knowledge and attitudes of invasive plants and Leave No Trace practices

Hiking and backpacking on American National Scenic Trails has increased in popularity in recent years. To encourage responsible and sustainable outdoor recreation on these much-loved trails, direct and indirect management strategies must be employed by managerial agencies. The Leave No Trace (LNT) education program aims to protect natural resources by promoting minimum-impact behaviours that lessen environmental impacts. The accidental introduction and dispersal of non-native invasive flora by hikers is little studied but can have a detrimental environmental impact on protected areas. The purpose of our study was to understand whether Appalachian Trail thru-hikers are: 1) aware of this problem, 2) adhering to LNT principles to reduce this problem, and 3) willing to learn and adopt minimum-impact behaviours to address this problem. We found that thru-hiker knowledge of invasive plants was limited and that very few thru-hikers adopted low-impact practices to minimise plant introduction and spread. Promisingly, we found that most thru-hikers, once aware of the problems, were willing to learn and apply low-impact practices to minimise plant introduction and spread. We discuss the barriers to their adoption of these behaviours and present a comprehensive list of suggested LNT practices to limit invasive plant introduction and spread. We conclude that, whilst challenging, protected area managers can help deter the spread of invasive plants along trails by improving educational messaging, signage, personal communication, and providing supporting infrastructure that encourages visitors to adopt specific practices to minimise invasive plant introduction and spread within protected areas.

Virginia, West Virginia↗

USGS2017

A discrete stage-structured model of California newt population dynamics during a period of drought

We introduce a mathematical model for studying the population dynamics under drought of the California newt ( Taricha torosa ), a species of special concern in the state of California. Since 2012, California has experienced a record-setting drought, and multiple studies predict drought conditions currently underway will persist and even increase in severity. Recent declines and local extinctions of California newt populations in Santa Monica Mountain streams motivate our study of the impact of drought on newt population sizes. Although newts are terrestrial salamanders, they migrate to streams each spring to breed and lay eggs. Since egg and larval stages occur in water, a precipitation deficit due to drought conditions reduces the space for newt egg-laying and the necessary habitat for larval development. To mathematically forecast newt population dynamics, we develop a nonlinear system of discrete equations that includes demographic parameters such as survival rates for newt life stages and egg production, which depend on habitat availability and rainfall. We estimate these demographic parameters using 15 years of stream survey data collected from Cold Creek in Los Angeles County, California, and our model captures the observed decline of the parameterized Cold Creek newt population. Based upon data analysis, we predict how the number of available newt egg-laying sites varies with annual precipitation. Our model allows us to make predictions about how the length and severity of drought can affect the likelihood of persistence and the time to critical endangerment of a local newt population. We predict that sustained severe drought will critically endanger the newt population but that the newt population can rebound if a drought is sufficiently short.

Journal of Theoretical Biology↗

USGS2005

Anatahan, Northern Mariana Islands: Reconnaissance geological observations during and after the volcanic crisis of spring 1990, and monitoring prior to the May 2003 eruption

Anatahan island is 9.5 km east–west by 3.5 km north–south and truncated by an elongate caldera 5 km east–west by 2.5 km north–south. A steep-walled pit crater ∼1 km across and ∼200 m deep occupies the eastern part of the caldera. The island is the summit region of a mostly submarine stratovolcano. The oldest subaerial rocks (stage 1) are exposed low on the outer flanks and in the caldera walls. These include thick (∼10 m) and thin (2–3 m) lava flows, well-indurated tuffs, and scoria units that make up the bulk of the island. Rock compositions range from basaltic andesite to dacite, and most are plagioclase-phyric. On the steep north and south flanks of the volcano, these rocks are cut by numerous east–west-oriented, few-hundred-m-long lineaments of undetermined origin. Indurated breccias unconformably overlie scarps cut into stage 1 units low on the south flank. Intermediate-age eruptive units (stage 2) include caldera-filling lava flows and pyroclastic deposits and, on the outer flanks, vents and valley-filling lava flows. The youngest pre-2003 volcanic unit on Anatahan (stage 3) is a hydromagmatic surge and fall deposit rich in accretionary lapilli. Prior to 2003, this unit was found over almost the entire island, and in many places original depositional surfaces and outcrops could be found in high-energy environments along the coast, indicating a young (but undetermined) age. During reconnaissance visits in 1990, 1992, 1994, and 2001, geothermal activity (fumaroles as well as pits with boiling, sediment-laden pools) was observed in the southern part of the pit crater. In March and April 1990, increased local seismicity, a large regional earthquake, and reported increased fumarolic activity in the pit crater prompted evacuation of Anatahan village, at the west end of the island. Our first field investigation took place in late April 1990 to assess the level of volcanic unrest, conduct reconnaissance geological observations, collect rock and geothermal water samples, and set up a geophysical monitoring network. Results at this time were inconclusive with respect to determining whether the activity was anomalous. Water in some of the geothermal pits within the pit crater was boiling, and pH values as low as 0.7 were recorded in the field. An electronic distance measurement (EDM) network was installed, and over a ∼1-week period, up to 9 cm of extension occurred across some lines but not others. Seismicity was characterized by intermittent local earthquakes but no sustained swarms or tremor. A brief visit in June 1990 revealed that the shallow lake near the boiling pits was gone, but activity in the pits themselves was similar to that of April 1990. Only minor extension had occurred along a single EDM line since the previous visit, and no earthquakes >M2.5 occurred during the visit. Subsequent 1- to 2-day visits occurred in October 1990, May 1992, May 1994, and June 2001. Activity within the geothermal pits was relatively constant during every visit, although during this 11-year period the level of the water in each pit decreased. In June 2001, a ∼50-m-wide region of mud pots and steaming ground in the central part of the geothermal area had developed. No geologic evidence, however, suggested that an eruption would occur < 2 years afterward. Most of the EDM lines showed slight extension between late 1990 and 1992, followed by very gradual contraction from 1992 to 2001. A more extensive seismic-monitoring system was installed on the Northern Mariana Islands during these visits, and it recorded a small seismic swarm at Anatahan from May to July 1993. The telemetry component of the seismic equipment broke prior to 2001 and had not been repaired by the time of the May 2003 eruption, so no precursory seismic data were recorded to indicate pre-eruption unrest.

Northern Mariana Islands↗

USGS2005

The 2003 phreatomagmatic eruptions of Anatahan volcano - Textural and petrologic features of deposits at an emergent island volcano

Stratigraphic and field data are used in conjunction with textural and chemical evidence (including data from scanning electron microscope, electron microprobe, X-ray fluorescence, X-ray diffraction, and instrumental neutron activation analysis) to establish that the 2003 eruption of Anatahan volcano was mainly phreatomagmatic, dominated by explosive interaction of homogeneous composition low-viscosity crystal-poor andesite magma with water. The hydromagmatic mode of eruption contributed to the significant height of initial eruptive columns and to the excavation and eruption of altered rock debris from the sub-volcanic hydrothermal system. Volatile contents of glass inclusions in equilibrium phenocrysts less abundances of these constituents in matrix glass times the estimated mass of juvenile magma indicate minimum emissions of 19 kt SO2 and 13 kt Cl. This petrologic estimate of SO2 emission is an order-of-magnitude less than an estimate from TOMS. Similarly, inferred magma volumes from the petrologic data are an order of magnitude greater than those modeled from deformation data. Both discrepancies indicate additional sources of volatiles, likely derived from a separate fluid phase in the magma. The paucity of near-source volcanic-tectonic earthquakes preceding the eruption, and the dominance of sustained long-period tremor are attributed to the ease of ascent of the hot low-viscosity andesite, followed by a shallow phreatomagmatic mode of eruption. Phreatomagmatic eruptions are probably more common at emergent tropical island volcanoes, where shallow fresh-water lenses occur at near-sea-level vents. These relations suggest that phreatomagmatic explosions contributed to the formation of many of the near-sea-level craters and possibly even to the small calderas at the other Mariana islands.

Northern Mariana Islands↗

USGS2007

Swarms of repeating long-period earthquakes at Shishaldin Volcano, Alaska, 2001-2004

During 2001–2004, a series of four periods of elevated long-period seismic activity, each lasting about 1–2 months, occurred at Shishaldin Volcano, Aleutian Islands, Alaska. The time periods are termed swarms of repeating events , reflecting an abundance of earthquakes with highly similar waveforms that indicate stable, non-destructive sources. These swarms are characterized by increased earthquake amplitudes, although the seismicity rate of one event every 0.5–5 min has remained more or less constant since Shishaldin last erupted in 1999. A method based on waveform cross-correlation is used to identify highly repetitive events, suggestive of spatially distinct source locations. The waveform analysis shows that several different families of similar events co-exist during a given swarm day, but generally only one large family dominates. A network of hydrothermal fractures may explain the events that do not belong to a dominant repeating event group, i.e. multiple sources at different locations exist next to a dominant source. The dominant waveforms exhibit systematic changes throughout each swarm, but some of these waveforms do reappear over the course of 4 years indicating repeatedly activated source locations. The choked flow model provides a plausible trigger mechanism for the repeating events observed at Shishaldin, explaining the gradual changes in waveforms over time by changes in pressure gradient across a constriction within the uppermost part of the conduit. The sustained generation of Shishaldin's long-period events may be attributed to complex dynamics of a multi-fractured hydrothermal system: the pressure gradient within the main conduit may be regulated by temporarily sealing and reopening of parallel flow pathways, by the amount of debris within the main conduit and/or by changing gas influx into the hydrothermal system. The observations suggest that Shishaldin's swarms of repeating events represent time periods during which a dominant source is activated.

Alaska↗

USGS2011

Evidence of magma intrusion at Fourpeaked volcano, Alaska in 2006-2007 from a rapid-response seismic network and volcanic gases

On September 17th, 2006, Fourpeaked volcano had a widely-observed phreatic eruption. At the time, Fourpeaked was an unmonitored volcano with no known Holocene activity, based on limited field work. Airborne gas sampling began within days of the eruption and a modest seismic network was installed in stages. Vigorous steaming continued for months; however, there were no further eruptions similar in scale to the September 17 event. This eruption was followed by several months of sustained seismicity punctuated by vigorous swarms, and SO 2 emissions exceeding a thousand tons/day. Based on observations during and after the phreatic eruption, and assuming no recent pre-historical eruptive activity at Fourpeaked, we propose that the activity was caused by a minor injection of new magma at or near 5km depth beneath Fourpeaked, which remained active over several months as this magma equilibrated into the crust. By early 2007 declining seismicity and SO 2 emission signaled the end of unrest. Because the Fourpeaked seismic network was installed in stages and the seismicity was punctuated by discrete swarms, we use Fourpeaked to illustrate quantitatively the efficacy and shortcomings of rapid response seismic networks for tracking volcanic earthquakes.

Alaska↗

USGS2013

Contrasting patterns of vesiculation in low, intermediate, and high Hawaiian fountains: A case study of the 1969 Mauna Ulu eruption

Hawaiian-style eruptions, or Hawaiian fountains, typically occur at basaltic volcanoes and are sustained, weakly explosive jets of gas and dominantly coarse, juvenile ejecta (dense spatter to delicate reticulite). Almost the entire range of styles and mass eruption rates within Hawaiian fountaining occurred during twelve fountaining episodes recorded at Mauna Ulu, Kīlauea between May and December 1969. Such diversity in intensity and style is controlled during magma ascent by many processes that can be constrained by the size and shape of vesicles in the 1969 pyroclasts. This paper describes pyroclast vesicularity from high, intermediate, and low fountaining episodes with eruption rates from 0.05 to 1.3 × 10 6 m 3 h − 1 . As each eruptive episode progressed, magma ascent slowed in and around the vent system, offering extended time for bubbles to grow and coalesce. Late ejected pyroclasts are thus characterized by populations of fewer and larger vesicles with relaxed shapes. This progression continued in the intervals between episodes after termination of fountain activity. The time scale for this process of shallow growth, coalescence and relaxation of bubbles is typically tens of hours. Rims and cores of pumiceous pyroclasts from moderate to high fountaining episodes record a second post-fragmentation form of vesicle maturation. Partially thermally insulated pyroclasts can have internal bubble populations evolve more dynamically with continued growth and coalescence, on a time scale of only minutes, during transport in the fountains. Reticulite, which formed in a short-lived fountain 540 m in height, underwent late, short-lived bubble nucleation followed by rapid growth of a uniform bubble population in a thermally insulated fountain, and quenched at the onset of permeability before significant coalescence. These contrasting patterns of shallow degassing and outgassing were the dominant controls in determining both the form and duration of fountaining episodes at Mauna Ulu, and probably for many other Hawaiian-style eruptions.

Hawaii↗

USGS2016

Quantification of the intrusion process at Kīlauea volcano, Hawai'i

The characteristic size of two types of intrusions identified beneath Kīlauea's East Rift zone are uniquely estimated by combining time constraints from fractional crystallization and the rates of magma solidification during cooling. Some intrusions were rapidly emplaced as dikes, but stalled before reaching the surface, and cooled and crystallized to feed later fractionated eruptions. More specifically, using the observed time interval between initial emplacement and eruption of fractionated lava, whose degree of fractionation is estimated from petrologic mixing calculations, the extent of solidification or cooling needed to produce this amount of fractionation can be directly inferred. And from the known erupted volumes the spatial extent or size of this fractionated volume can be analytically related to the full size of the source body itself. Two examples yield dike widths of 82 and 68 m. Other intrusions remain close to the east rift magma transport path and are observed to last for decades or longer as viable magma bodies that may participate in feeding later eruptions. The thickness of semi-permanent reservoirs near the East Rift Zone magma transport path can be estimated by assuming a resupply rate that is sufficiently frequent to restrict cooling to < 10 °C. It is inferred that both types of intrusions likely began as dike offshoots from the East Rift Zone magma transport path, but the frequently resupplied bodies may have later been converted to sills or laccoliths of heights estimated at 43–62 m. Our modeled intrusions contrast with models of rapidly emplaced thinner dikes feeding shallow intrusions, which are accompanied by intense rift earthquake swarms and are often associated with eruptions. These calculations show that long-term heating of the wallrock of the magma transport paths serves to slow conduit cooling, which may be partly responsible for sustaining long East Rift Zone eruptions. Adjacent to the vertical transport path beneath Kīlauea's summit, the combined effects of heating and ever-increasing magma supply rate may have forced a commensurate enlarging of the conduit, perhaps explaining the occurrence of a temporary burst of deep (5–15 km) long-period earthquake swarms between 1987 and 1992.

Hawaii↗