SEARCH · Geology Reports

Results for “World & I”

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,711 records · Page 95Linked to original sources

USGS2008

Assessment of crop growth and soil water modules in SWAT2000 using extensive field experiment data in an irrigation district of the Yellow River Basin

SWAT, a physically-based, hydrological model simulates crop growth, soil water and groundwater movement, and transport of sediment and nutrients at both the process and watershed scales. While the different versions of SWAT have been widely used throughout the world for agricultural and water resources applications, little has been done to test the performance, variability, and transferability of the parameters in the crop growth, soil water, and groundwater modules in an integrated way with multiple sets of field experimental data at the process scale. Using an multiple years of field experimental data of winter wheat (Triticum aestivum L.) in the irrigation district of the Yellow River Basin, this paper assesses the performance of the plant-soil-groundwater modules and the variability and transferability of SWAT2000. Comparison of the simulated results by SWAT to the observations showed that SWAT performed quite unsatisfactorily in LAI predictions during the senescence stage, in yield predictions, and in soil-water estimation under dry soil-profile conditions. The unsatisfactory performance in LAI prediction might be attributed to over-simplified senescence modeling; in yield prediction to the improper computation of the harvest index; and in soil water under dry conditions to the exclusion of groundwater evaporation from the soil water balance in SWAT. In this paper, improvements in crop growth, soil water, and groundwater modules in SWAT were implemented. The saturated soil profile was coupled to the oscillating groundwater table. A variable evaporation coefficient taking into account soil water deficit index, groundwater depth, and crop root depth was used to replace the fixed coefficient in computing groundwater evaporation. The soil water balance included the groundwater evaporation. The modifications improved simulations of crop evapotranspiration and biomass as well as soil water dynamics under dry soil-profile conditions. The evaluation shows that the crop growth and soil water components of SWAT could be further refined to better simulate the hydrology of agricultural watersheds. ?? 2008 Elsevier B.V. All rights reserved.

Journal of Hydrology↗

USGS2013

Fens as whole-ecosystem gauges of groundwater recharge under climate change

Currently, little is known about the impact of climate change on groundwater recharge in the Sierra Nevada and southern Cascade Range of California or other mountainous regions of the world. The purpose of this study was to determine whether small alpine peat lands called fens can be used as whole-ecosystem gauges of groundwater recharge through time. Fens are sustained by groundwater discharge and are highly sensitive to changes in groundwater flow due to hydrologic disturbance including climate change. Seven fens in the Sierra Nevada and southern Cascade Range were studied over a 50-80 year period using historic aerial photography. In each aerial photograph, fen areas were identified as open lawn and partially treed areas that exhibited (1) dark brownish-green coloring or various shades of gray and black in black and white imagery and (2) mottling of colors and clustering of vegetation, which signified a distinct moss canopy with overlying clumped sedge vegetation. In addition to the aerial photography study, a climate analysis for the study sites was carried out using both measured data (U.S. Department of Agriculture Natural Resources Conservation Service SNOwpack TELemetry system) and modeled data (a downscaled version of the Parameter-elevation Regressions on Independent Slopes Model) for the period from 1951 to 2010. Over the study period, the five fens in the Sierra Nevada were found to be decreasing between 10% and 16% in delineated area. The climate analysis revealed significant increases through time in annual mean minimum temperature (T min ) between 1951-1980 and 1981-2010. In addition, April 1 snow water equivalent and snowpack longevity also decreased between 1951-1980 and 1981-2010. For the fens in the Cascade Range, there were no discernible changes in delineated area. At these sites, increases in T min occurred only within the past 20-25 years and decreases in snowpack longevity were more subtle. A conceptual model is presented, which illustrates that basic differences in hydrogeology of the Sierra Nevada vs. the Cascade Range may control the threshold at which changes in delineated fen areas are discernible. Overall, the results from this study show that fens in the Sierra Nevada have strong potential as whole ecosystem gauges for determining long-term changes in groundwater recharge under climate change. Due to either more moderate climate change and/or hydrogeological differences, fens in the southern Cascade Range currently do not appear to have the same utility. A greater sample size of fens in the Sierra Nevada is needed to confirm the general applicability of this method. In addition, future work needs to focus on integrating fen monitoring with geochemical and/or isotopic process-level studies in order to quantify changes in groundwater recharge identified using this new approach.

Journal of Hydrology↗

USGS2013

Appraising options to reduce shallow groundwater tables and enhance flow conditions over regional scales in an irrigated alluvial aquifer system

Some of the world’s key agricultural production systems face big challenges to both water quantity and quality due to shallow groundwater that results from long-term intensive irrigation, namely waterlogging and salinity, water losses, and environmental problems. This paper focuses on water quantity issues, presenting finite-difference groundwater models developed to describe shallow water table levels, non-beneficial groundwater consumptive use, and return flows to streams across two regions within an irrigated alluvial river valley in southeastern Colorado, USA. The models are calibrated and applied to simulate current baseline conditions in the alluvial aquifer system and to examine actions for potentially improving these conditions. The models provide a detailed description of regional-scale subsurface unsaturated and saturated flow processes, thereby enabling detailed spatiotemporal description of groundwater levels, recharge to infiltration ratios, partitioning of ET originating from the unsaturated and saturated zones, and groundwater flows, among other variables. Hybrid automated and manual calibration of the models is achieved using extensive observations of groundwater hydraulic head, groundwater return flow to streams, aquifer stratigraphy, canal seepage, total evapotranspiration, the portion of evapotranspiration supplied by upflux from the shallow water table, and irrigation flows. Baseline results from the two regional-scale models are compared to model predictions under variations of four alternative management schemes: (1) reduced seepage from earthen canals, (2) reduced irrigation applications, (3) rotational lease fallowing (irrigation water leased to municipalities, resulting in temporary dry-up of fields), and (4) combinations of these. The potential for increasing the average water table depth by up to 1.1 and 0.7 m in the two respective modeled regions, thereby reducing the threat of waterlogging and lowering non-beneficial consumptive use from adjacent fallow and naturally-vegetated lands, is demonstrated for the alternative management intervention scenarios considered. Net annual average savings of up to about 9.9 million m 3 (8000 ac ft) and 2.3 million m 3 (1900 ac ft) of non-beneficial groundwater consumptive use is demonstrated for the study periods in each of the two respective study regions. Alternative water management interventions achieve varying degrees of benefits in each of the two regions, suggesting a need to adopt region-specific interventions and avoid a ‘one-size-fits-all’ approach. Impacts of the considered interventions on return flows to the river were predicted to be significant, highlighting the need for flow augmentation to comply with an interstate river compact and portending beneficial impacts on solute loading.

Colorado↗

USGS2014

Surface-groundwater interactions in hard rocks in Sardon Catchment of western Spain: an integrated modeling approach

The structural and hydrological complexity of hard rock systems (HRSs) affects dynamics of surface–groundwater interactions. These complexities are not well described or understood by hydrogeologists because simplified analyses typically are used to study HRSs. A transient, integrated hydrologic model (IHM) GSFLOW (Groundwater and Surface water FLOW) was calibrated and post-audited using 18 years of daily groundwater head and stream discharge data to evaluate the surface–groundwater interactions in semi-arid, ∼80 km 2 granitic Sardon hilly catchment in Spain characterized by shallow water table conditions, relatively low storage, dense drainage networks and frequent, high intensity rainfall. The following hydrological observations for the Sardon Catchment, and more generally for HRSs were made: (i) significant bi-directional vertical flows occur between surface water and groundwater throughout the HRSs; (ii) relatively large groundwater recharge represents 16% of precipitation (P, 562 mm.y −1 ) and large groundwater exfiltration (∼11% of P) results in short groundwater flow paths due to a dense network of streams, low permeability and hilly topographic relief; deep, long groundwater flow paths constitute a smaller component of the water budget (∼1% of P); quite high groundwater evapotranspiration (∼5% of P and ∼7% of total evapotranspiration); low permeability and shallow soils are the main reasons for relatively large components of Hortonian flow and interflow (15% and 11% of P, respectively); (iii) the majority of drainage from the catchment leaves as surface water; (iv) declining 18 years trend (4.44 mm.y −1 ) of groundwater storage; and (v) large spatio-temporal variability of water fluxes. This IHM study of HRSs provides greater understanding of these relatively unknown hydrologic systems that are widespread throughout the world and are important for water resources in many regions.

Sardon Catchment↗

USGS2014

Identifying dominant controls on hydrologic parameter transfer from gauged to ungauged catchments: a comparative hydrology approach

Daily streamflow information is critical for solving various hydrologic problems, though observations of continuous streamflow for model calibration are available at only a small fraction of the world’s rivers. One approach to estimate daily streamflow at an ungauged location is to transfer rainfall–runoff model parameters calibrated at a gauged (donor) catchment to an ungauged (receiver) catchment of interest. Central to this approach is the selection of a hydrologically similar donor. No single metric or set of metrics of hydrologic similarity have been demonstrated to consistently select a suitable donor catchment. We design an experiment to diagnose the dominant controls on successful hydrologic model parameter transfer. We calibrate a lumped rainfall–runoff model to 83 stream gauges across the United States. All locations are USGS reference gauges with minimal human influence. Parameter sets from the calibrated models are then transferred to each of the other catchments and the performance of the transferred parameters is assessed. This transfer experiment is carried out both at the scale of the entire US and then for six geographic regions. We use classification and regression tree (CART) analysis to determine the relationship between catchment similarity and performance of transferred parameters. Similarity is defined using physical/climatic catchment characteristics, as well as streamflow response characteristics (signatures such as baseflow index and runoff ratio). Across the entire US, successful parameter transfer is governed by similarity in elevation and climate, and high similarity in streamflow signatures. Controls vary for different geographic regions though. Geology followed by drainage, topography and climate constitute the dominant similarity metrics in forested eastern mountains and plateaus, whereas agricultural land use relates most strongly with successful parameter transfer in the humid plains.

Journal of Hydrology↗

USGS2014

Developing and testing temperature models for regulated systems: a case study on the Upper Delaware River

Water temperature is an important driver of many processes in riverine ecosystems. If reservoirs are present, their releases can greatly influence downstream water temperatures. Models are important tools in understanding the influence these releases may have on the thermal regimes of downstream rivers. In this study, we developed and tested a suite of models to predict river temperature at a location downstream of two reservoirs in the Upper Delaware River (USA), a section of river that is managed to support a world-class coldwater fishery. Three empirical models were tested, including a Generalized Least Squares Model with a cosine trend (GLScos), AutoRegressive Integrated Moving Average (ARIMA), and Artificial Neural Network (ANN). We also tested one mechanistic Heat Flux Model (HFM) that was based on energy gain and loss. Predictor variables used in model development included climate data (e.g., solar radiation, wind speed, etc.) collected from a nearby weather station and temperature and hydrologic data from upstream U.S. Geological Survey gages. Models were developed with a training dataset that consisted of data from 2008 to 2011; they were then independently validated with a test dataset from 2012. Model accuracy was evaluated using root mean square error (RMSE), Nash Sutcliffe efficiency (NSE), percent bias (PBIAS), and index of agreement (d) statistics. Model forecast success was evaluated using baseline-modified prime index of agreement (md) at the one, three, and five day predictions. All five models accurately predicted daily mean river temperature across the entire training dataset (RMSE = 0.58–1.311, NSE = 0.99–0.97, d = 0.98–0.99); ARIMA was most accurate (RMSE = 0.57, NSE = 0.99), but each model, other than ARIMA, showed short periods of under- or over-predicting observed warmer temperatures. For the training dataset, all models besides ARIMA had overestimation bias (PBIAS = −0.10 to −1.30). Validation analyses showed all models performed well; the HFM model was the most accurate compared other models (RMSE = 0.92, both NSE = 0.98, d = 0.99) and the ARIMA model was least accurate (RMSE = 2.06, NSE = 0.92, d = 0.98); however, all models had an overestimation bias (PBIAS = −4.1 to −10.20). Aside from the one day forecast ARIMA model (md = 0.53), all models forecasted fairly well at the one, three, and five day forecasts (md = 0.77–0.96). Overall, we were successful in developing models predicting daily mean temperature across a broad range of temperatures. These models, specifically the GLScos, ANN, and HFM, may serve as important tools for predicting conditions and managing thermal releases in regulated river systems such as the Delaware River. Further model development may be important in customizing predictions for particular biological or ecological needs, or for particular temporal or spatial scales.

Delaware, New York, Pennsylvania↗

USGS2021

Uncertainty in remote sensing of streams using noncontact radars

Accounting for freshwater resources and monitoring floods are vital functions for societies throughout the world. Remote-sensing methods offer great prospects to expand stream monitoring in developing countries and to smaller, headwater streams that are largely ungauged worldwide. This study evaluates the potential to estimate discharge using eight radar units that have been installed over streams in diverse hydrologic and hydraulic settings across the United States. The research highlights error characteristics associated with the measurements of stage using pulsed wave radars, mean channel velocity from continuous wave Doppler radars, and their combined use to estimate discharge at sites that were collocated with conventional streamgauges. Potential stage biases caused by the thermal expansion and contraction of supporting structures due to diurnal temperature changes were examined. A dry concrete, flume showed the temperature-dependent stage variations were no more than 2 cm. Surface velocity retrievals needed to be adjusted to represent the mean channel velocity when estimating discharge. Different approaches were evaluated and application of two different, depth-dependent adjustment factors was found to yield the most accurate estimates. This study found that it is possible to get accurate discharge estimates from noncontact radar measurements, providing cost-effective solutions for remote sensing of ungauged streams. Lastly, radar measurements of the raw variables (i.e., stage and surface velocity) can be used in an early alerting context to detect flash floods in ungauged streams.

Arizona, Colorado, Oklahoma, Texas↗

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↗

USGS2013

Updated methodology for nuclear magnetic resonance characterization of shales

Unconventional petroleum resources, particularly in shales, are expected to play an increasingly important role in the world’s energy portfolio in the coming years. Nuclear magnetic resonance (NMR), particularly at low-field, provides important information in the evaluation of shale resources. Most of the low-field NMR analyses performed on shale samples rely heavily on standard T 1 and T 2 measurements. We present a new approach using solid echoes in the measurement of T 1 and T 1 –T 2 correlations that addresses some of the challenges encountered when making NMR measurements on shale samples compared to conventional reservoir rocks. Combining these techniques with standard T 1 and T 2 measurements provides a more complete assessment of the hydrogen-bearing constituents (e.g., bitumen, kerogen, clay-bound water) in shale samples. These methods are applied to immature and pyrolyzed oil shale samples to examine the solid and highly viscous organic phases present during the petroleum generation process. The solid echo measurements produce additional signal in the oil shale samples compared to the standard methodologies, indicating the presence of components undergoing homonuclear dipolar coupling. The results presented here include the first low-field NMR measurements performed on kerogen as well as detailed NMR analysis of highly viscous thermally generated bitumen present in pyrolyzed oil shale.

Journal of Magnetic Resonance↗

USGS2016

Value-focused framework for defining landscape-scale conservation targets

Conservation of natural resources can be challenging in a rapidly changing world and require collaborative efforts for success. Conservation planning is the process of deciding how to protect, conserve, and enhance or minimize loss of natural and cultural resources. Establishing conservation targets (also called indicators or endpoints), the measurable expressions of desired resource conditions, can help with site-specific up to landscape-scale conservation planning. Using conservation targets and tracking them through time can deliver benefits such as insight into ecosystem health and providing early warnings about undesirable trends. We describe an approach using value-focused thinking to develop statewide conservation targets for Florida. Using such an approach allowed us to first identify stakeholder objectives and then define conservation targets to meet those objectives. Stakeholders were able to see how their shared efforts fit into the broader conservation context, and also anticipate the benefits of multi-agency and -organization collaboration. We developed an iterative process for large-scale conservation planning that included defining a shared framework for the process, defining the conservation targets themselves, as well as developing management and monitoring strategies for evaluation of their effectiveness. The process we describe is applicable to other geographies where multiple parties are seeking to implement collaborative, large-scale biological planning.

Journal for Nature Conservation↗

USGS2024

Grand challenges in anticipating and responding to critical materials supply risks

Critical materials are resources that are vulnerable to supply disruptions, where those disruptions can have significant adverse impacts on society. In the coming years, materials supply risks associated with the energy transition and geopolitics are likely to intensify and new risks are expected to emerge. This perspective identifies three “Grand Challenges” that represent frontier areas for critical materials research and highlights some promising new directions for each area: (1) extending visibility downstream to value-added materials beyond elemental forms; (2) quantifying the risks associated with market dynamics; and (3) developing tools to inform policy interventions. Emerging digital capabilities have the potential to play a significant role addressing long-standing limitations in data quality and access to unlock progress on these challenges. Progress in these areas can equip decision-makers across industry, government, and finance with tools to understand the complexity and uncertainty introduced by these real-world challenges.

Joule↗

USGS2007

Geochemistry of rare earth elements in Permian coals from the Huaibei Coalfield, China

The rare earth elements (REEs) in coals are important because of: (a) REE patterns can be an indicator of the nature of source rocks of the mineral matter as well as sedimentary environments; (b) REEs abundance in coal may have industrial-significance. In this study, a total of thirty-four samples of Permian coal, partings, roof, and floor were collected from the Huaibei Coalfield, Anhui Province, China. Abundances of rare earth elements (REEs) and other elements in the samples were determined by inductively coupled-plasma mass spectrometry (ICP-MS) and inductively coupled-plasma atomic emission spectrometry (ICP-AES). The results show that the REEs are enriched in coals in the Huaibei Coalfield as compared with Chinese and U.S. coals and the world coal average. Coals in the Lower Shihezi Formation (No. 7, 5, and 4 Coals) and Upper Shihezi Formation (No. 3) have higher REE abundances than the coals in Shanxi Formation (No. 10). Magmatic intrusion resulted in high enrichment of REEs concentrations in No. 5 and 7 Coals. The REE abundances are positively correlated with the ash content. The mineral matter in these coals is mainly made up of clay minerals and carbonates. The REEs are positively correlated with lithophile elements including Si, Al, Ti, Fe, and Na, which are mainly distributed in clay minerals, indicating that REEs are contained mainly in clay minerals. The REE abundances in coals normalized by the ash are higher than that in partings. REEs abundances of coals cannot be accounted for by the REE content in the mineral matter, and some REEs associated with organic matter in coals. ?? 2007 Elsevier Ltd. All rights reserved.

Journal of Asian Earth Sciences↗

USGS2009

A multidisciplinary effort to assign realistic source parameters to models of volcanic ash-cloud transport and dispersion during eruptions

During volcanic eruptions, volcanic ash transport and dispersion models (VATDs) are used to forecast the location and movement of ash clouds over hours to days in order to define hazards to aircraft and to communities downwind. Those models use input parameters, called “eruption source parameters”, such as plume height H , mass eruption rate Ṁ , duration D , and the mass fraction m 63 of erupted debris finer than about 4ϕ or 63 μm, which can remain in the cloud for many hours or days. Observational constraints on the value of such parameters are frequently unavailable in the first minutes or hours after an eruption is detected. Moreover, observed plume height may change during an eruption, requiring rapid assignment of new parameters. This paper reports on a group effort to improve the accuracy of source parameters used by VATDs in the early hours of an eruption. We do so by first compiling a list of eruptions for which these parameters are well constrained, and then using these data to review and update previously studied parameter relationships. We find that the existing scatter in plots of H versus Ṁ yields an uncertainty within the 50% confidence interval of plus or minus a factor of four in eruption rate for a given plume height. This scatter is not clearly attributable to biases in measurement techniques or to well-recognized processes such as elutriation from pyroclastic flows. Sparse data on total grain-size distribution suggest that the mass fraction of fine debris m 63 could vary by nearly two orders of magnitude between small basaltic eruptions (∼ 0.01) and large silicic ones (> 0.5). We classify eleven eruption types; four types each for different sizes of silicic and mafic eruptions; submarine eruptions; “brief” or Vulcanian eruptions; and eruptions that generate co-ignimbrite or co-pyroclastic flow plumes. For each eruption type we assign source parameters. We then assign a characteristic eruption type to each of the world's ∼ 1500 Holocene volcanoes. These eruption types and associated parameters can be used for ash-cloud modeling in the event of an eruption, when no observational constraints on these parameters are available.

Journal of Volcanology and Geothermal Research↗

USGS2012

The 2010 explosive eruption of Java's Merapi volcano—A ‘100-year’ event

Merapi volcano (Indonesia) is one of the most active and hazardous volcanoes in the world. It is known for frequent small to moderate eruptions, pyroclastic flows produced by lava dome collapse, and the large population settled on and around the flanks of the volcano that is at risk. Its usual behavior for the last decades abruptly changed in late October and early November 2010, when the volcano produced its largest and most explosive eruptions in more than a century, displacing at least a third of a million people, and claiming nearly 400 lives. Despite the challenges involved in forecasting this ‘hundred year eruption’, we show that the magnitude of precursory signals (seismicity, ground deformation, gas emissions) was proportional to the large size and intensity of the eruption. In addition and for the first time, near-real-time satellite radar imagery played an equal role with seismic, geodetic, and gas observations in monitoring eruptive activity during a major volcanic crisis. The Indonesian Center of Volcanology and Geological Hazard Mitigation (CVGHM) issued timely forecasts of the magnitude of the eruption phases, saving 10,000–20,000 lives. In addition to reporting on aspects of the crisis management, we report the first synthesis of scientific observations of the eruption. Our monitoring and petrologic data show that the 2010 eruption was fed by rapid ascent of magma from depths ranging from 5 to 30 km. Magma reached the surface with variable gas content resulting in alternating explosive and rapid effusive eruptions, and released a total of ~ 0.44 Tg of SO 2 . The eruptive behavior seems also related to the seismicity along a tectonic fault more than 40 km from the volcano, highlighting both the complex stress pattern of the Merapi region of Java and the role of magmatic pressurization in activating regional faults. We suggest a dynamic triggering of the main explosions on 3 and 4 November by the passing seismic waves generated by regional earthquakes on these days.

Merapi volcano↗

USGS2015

An automated SO2 camera system for continuous, real-time monitoring of gas emissions from Kīlauea Volcano's summit Overlook Crater

SO 2 camera systems allow rapid two-dimensional imaging of sulfur dioxide (SO 2 ) emitted from volcanic vents. Here, we describe the development of an SO 2 camera system specifically designed for semi-permanent field installation and continuous use. The integration of innovative but largely “off-the-shelf” components allowed us to assemble a robust and highly customizable instrument capable of continuous, long-term deployment at Kīlauea Volcano's summit Overlook Crater. Recorded imagery is telemetered to the USGS Hawaiian Volcano Observatory (HVO) where a novel automatic retrieval algorithm derives SO 2 column densities and emission rates in real-time. Imagery and corresponding emission rates displayed in the HVO operations center and on the internal observatory website provide HVO staff with useful information for assessing the volcano's current activity. The ever-growing archive of continuous imagery and high-resolution emission rates in combination with continuous data from other monitoring techniques provides insight into shallow volcanic processes occurring at the Overlook Crater. An exemplary dataset from September 2013 is discussed in which a variation in the efficiency of shallow circulation and convection, the processes that transport volatile-rich magma to the surface of the summit lava lake, appears to have caused two distinctly different phases of lake activity and degassing. This first successful deployment of an SO 2 camera for continuous, real-time volcano monitoring shows how this versatile technique might soon be adapted and applied to monitor SO 2 degassing at other volcanoes around the world.

Hawai'i↗

USGS2019

Gas and ash emissions associated with the 2010–present activity of Sinabung Volcano, Indonesia

Sinabung Volcano (Sumatra, Indonesia) awoke from over 1200 years of dormancy with multiple phreatic explosions in 2010. After a period of quiescence, Sinabung activity resumed in 2013, producing frequent explosions, lava dome extrusion, and pyroclastic flows from dome collapses, becoming one of the world's most active volcanoes and displacing over 20,000 citizens. This study presents a compilation of the geochemical datasets collected by the Indonesian Center for Volcanology and Geological Hazard Mitigation (CVGHM) from 2010 - current (2016), which provides insights into the evolution of the eruption. Based on observations of SO2 emissions, ash componentry, leachate chemistry, and bulk ash geochemistry, the eruption can be split into five distinct phases. The initial stage of phreatic summit explosions occurred from August - October 2010, during which background SO2 emissions averaged ~550 ± 180 t/d (1 s.d.). An eruptive pause (phase two) starting in October 2010 abruptly ended in September 2013 with a resumption of conduit-clearing eruptions. This third phase had a relatively modest background SO2 emission rate (avg. ~410 ± 275 t/d) and produced ash consisting entirely of accidental ejecta with high S/Cl leachate ratios (up to 30), suggestive of deep-sourced magma and the incorporation of hydrothermal sulfur-bearing phases. The most intense phase of the eruption (phase four) occurred from December 2013 to February 2014, when juvenile magma first reached the surface. This period included dozens of large eruptions per day, high SO2 emission rates (average: 1,120 ± 1,030 t/d, peak: ~3,800 t/d), the onset of lava dome extrusion, and a dramatic drop in S/Cl ash leachates to ratios < 5, all reflecting increased degassing from shallow magma and the clearing out of sulfurous phases from the old hydrothermal system. From late February 2014 through the time of writing (September 2016), Sinabung settled into a relatively steady state of lower activity (phase five). Ash emissions now consist of dominantly juvenile material, and background SO2 emission rates have been progressively decreasing to an average of ~250 - 300 t/d. Starting August 2016, SO2 emissions started being measured in a continuous manner using a network of permanent scanning DOAS instruments. We find that long-term SO2 emission rates have been gradually declining at Sinabung since early 2014, consistent with an apparent decrease in magma supply. Our degassing model suggests that large explosions and pyroclastic flows could continue in the near-term owing to conduit plugging and dome collapses, remaining a major threat until the magma supply rate decreases further and the eruption ends.

Sumatra↗

USGS2019

Globally detected volcanic lightning and umbrella dynamics during the 2014 eruption of Kelud, Indonesia

Volcanic lightning shows considerable promise as a monitoring and research tool to characterize explosive eruptions. Its key strengths are rapid and remote detection, because the radio signals produced by lightning can propagate thousands of km at the speed of light. Despite these tantalizing properties, the scientific work on volcanic lightning has only recently started gaining momentum. Much more is needed to understand what lightning reveals about the evolution of an eruption in near-real time. Here we examine the timing and energy release of lightning generated by the eruption of Kelud volcano in Indonesia on 13 February 2014, as detected by the World Wide Lightning Location Network (WWLLN). The eruption column reached at least 26 km above sea level, representing the highest plume since the advent of global lightning networks in the last decade. Therefore, it provides valuable constraints on the electrification of end-member, sustained Plinian columns. We investigate the lightning in context with satellite images, photographs, and other published studies. Results show that the earliest satellite-detected activity was a thermal anomaly at ~15:46 UTC, corresponding to a directed blast at the onset of eruption (and only a few lightning strokes). Following a brief pause, the eruption produced a sustained column and umbrella cloud that spread outward into the tropical stratosphere. Rates of umbrella expansion provide an average mass eruption rate (MER) in the range of 8 × 10 7 –1 × 10 8 kg s −1 . A more nuanced picture emerges from the time-varying MERs (determined between each satellite pass), which show rapid intensification during the first hour of eruption, followed by constant MER for about an hour, and waning toward the end (after ~17:50 UTC). At this stage, decreasing flux into the umbrella cloud coincides with column instability and formation of pyroclastic density currents, as recorded by photos from the ground ~17:45 UTC. We infer that some of the erupted mass partitioned into ground-hugging currents, leading to a lower apparent MER. Interestingly, there is not a 1:1 correlation between lightning intensity and MER over the course of eruption. Stroke rates increase sharply within the first 30–40 min (during rapid intensification of the plume), and then drop below 2 strokes per min once the MER remains constant. This suggests that electrification was controlled by the rate of increase in MER—in other words, the acceleration of particles out of the vent. We also show that lightning stroke-rates and energies are greatest within 50 km of the vent, even when the ash cloud extends >200 km downwind, indicating that lightning was focused in the regions of highest particle concentration and turbulence. Overall, we conclude that abrupt changes in lightning rates are clearly linked to changes in eruption behavior, and that rapid detection could aid monitoring efforts to characterize eruption rates or styles.

Journal of Volcanology and Geothermal Research↗

USGS2020

Geophysical characterization of the Northwest Geysers geothermal field, California

The Clear Lake Volcanic Field in northern California is the youngest and northern-most part of a long-lived volcanic system that has produced recent (~10 ka) eruptions. Adjacent to the Clear Lake Volcanic Field is the worlds largest energy producing geothermal field, The Geysers. The hottest part of The Geysers geothermal field is in the northwest where temperatures reach ~400 C at 3 km depth. Low permeability, high thermal gradients, and low steam saturation prescribed development of an enhanced geothermal system (EGS) in the Northwest Geysers to increase energy producing capacity. Though the Northwest Geysers is known to be the hottest part of the field, geophysical methods have failed to adequately image any inferred heat source. This project aims to image the heat source of the Northwest Geysers using newly collected gravity and magnetotelluric (MT) measurements. Gravity data were jointly modeled with existing magnetic data along a two-dimensional profile aligned with an existing geologic cross-section. The key feature of the potential field model is a low-density, low-susceptibility body at 5 km depth (bmsl) under the EGS. MT data were modeled in three-dimensions to characterize subsurface resistivity structure, where the upper 3 km of the resistivity model agrees well with existing data. Lithologic and steam saturation are estimated from modeled resistivity values using existing geophysical data. Below 3 km depth (bmsl), the resistivity model images a possible young intrusion under the EGS. A possible zone of partial melt (<5%) below 7 km depth (bmsl) in the northwestern part of the field is also imaged which extends northeast towards the main part of the Clear Lake Volcanic Field.

California↗