USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Volcanology and Geothermal Resources”

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.

10 recordsLinked to original sources
USGS2025

High-precision 40Ar/39Ar dating of Quaternary basalts from Auckland Volcanic Field, New Zealand, with implications for eruption rates and paleomagnetic correlations

The Auckland Volcanic Field (AVF), which last erupted ca. 550 years ago, is a late Quaternary monogenetic basaltic volcanic field (ca. 500 km 2 ) in the northern North Island of New Zealand. Prior to this study only 12 out of the 53 identified eruptive centres of the AVF had been reliably dated. Careful sample preparation and 40 Ar/ 39 Ar analysis has increased the number of well-dated centres in the AVF to 35. The high precision of the results is attributed to selection of fresh, non-vesicular, non-glassy samples from lava flow interiors. Sample selection was coupled with separation techniques that targeted only the groundmass of samples with < 5% glass and with groundmass feldspars > 10 μm wide, coupled with ten-increment furnace step-heating of large quantities (up to 200 mg) of material. The overall AVF age data indicate an onset at 193.2 ± 2.8 ka, an apparent six-eruption flare-up from 30 to 34 ka, and a ≤ 10 kyr hiatus between the latest and second-to-latest eruptions. Such non-uniformity shows that averaging the number of eruptions over the life-span of the AVF to yield a mean eruption rate is overly simplistic. Together with large variations in eruption volumes, and the large sizes and unusual chemistry within the latest eruptions (Rangitoto 1 and Rangitoto 2), our results illuminate a complex episodic eruption history. In particular, the rate of volcanism in AVF has increased since 60 ka, suggesting that the field is still in its infancy. Multiple centres with unusual paleomagnetic inclination and declination orientations are confirmed to fit into a number of geomagnetic excursions, with five identified in the Mono Lake, two within the Laschamp, one within the post-Blake or Blake, and two possibly within the Hilina Pali.

Auckland Volcanic Field↗
USGS2023

Elucidating the magma plumbing system of Ol Doinyo Lengai (Natron Rift, Tanzania) Using satellite geodesy and numerical modeling

Ol Doinyo Lengai, located in the southern Eastern Branch of the East African Rift had several eruptive episodes with ash falls and lava flows (VEI 3) that caused damage to the nearby communities between 2007 and 2010. The volcano is remote and access is difficult. Although this volcano has been studied for decades, its plumbing system is still poorly understood, in part, because of the lack of precise observations of surface deformation during periods of quiet and unrest. This study investigates the volcanic plumbing system of Ol Doinyo Lengai and its surroundings using data from the network of permanent Global Navigation Satellite System (GNSS) sites monitoring the volcano (the TZVOLCANO network) around the flanks of the volcano and Interferometric Synthetic Aperture Radar (InSAR) observations. We constrain surface motions using 6 GNSS sites distributed around Ol Doinyo Lengai, operating between 2016 and 2021, and InSAR data covering nearly the same time period. Because of the complex local tectonics, the interpretation of the deformation pattern is not straightforward. We first invert the GNSS deformation and InSAR observations independently to infer potential deformation sources. Then we perform a joint inversion of both GNSS and InSAR datasets to verify our findings. We compare the results from the joint inversion with the results from inverting each dataset independently. The GNSS, InSAR, and joint inversion results point to a deflating source, located east of Ol Doinyo Lengai and southwest of the dormant volcano Gelai at a depth of 3.49 ± 0.03 km (GNSS inversion), 5.2 ± 1.2 km (InSAR inversion) and 3.49 ± 0.06 km (joint inversion) relative to the summit (vent) and with a volume change ∆V of −0.04 ± 0.05 × 10 6 m 3 (GNSS inversion), −0.39 ± 0.29 × 10 6 m 3 (InSAR inversion), and − 0.04 ± 0.01 × 10 6 m 3 (joint inversion). Although this is non-unique modeling of geodetic datasets with small signals, the inversion results suggest that Ol Doinyo Lengai could be fed by an offset multi-reservoir system that includes a shallow magma reservoir (<5 km) east of Ol Doinyo Lengai, possibly connected to a deeper magma reservoir.

Natron Rift, Ol Doinyo Lengai↗
USGS2016

Multireaction equilibrium geothermometry: A sensitivity analysis using data from the Lower Geyser Basin, Yellowstone National Park, USA

A multireaction chemical equilibria geothermometry (MEG) model applicable to high-temperature geothermal systems has been developed over the past three decades. Given sufficient data, this model provides more constraint on calculated reservoir temperatures than classical chemical geothermometers that are based on either the concentration of silica (SiO 2 ), or the ratios of cation concentrations. A set of 23 chemical analyses from Ojo Caliente Spring and 22 analyses from other thermal features in the Lower Geyser Basin of Yellowstone National Park are used to examine the sensitivity of calculated reservoir temperatures using the GeoT MEG code (Spycher et al. 2013, 2014) to quantify the effects of solute concentrations, degassing, and mineral assemblages on calculated reservoir temperatures. Results of our analysis demonstrate that the MEG model can resolve reservoir temperatures within approximately ±15°C, and that natural variation in fluid compositions represents a greater source of variance in calculated reservoir temperatures than variations caused by analytical uncertainty (assuming ~5% for major elements). The analysis also suggests that MEG calculations are particularly sensitive to variations in silica concentration, the concentrations of the redox species Fe(II) and H 2 S, and that the parameters defining steam separation and CO 2 degassing from the liquid may be adequately determined by numerical optimization. Results from this study can provide guidance for future applications of MEG models, and thus provide more reliable information on geothermal energy resources during exploration.

Idaho, Montana, Wyoming↗
USGS1976

Submarine geothermal resources

Approximately 20% of the earth's heat loss (or 2 ?? 1012 cal/s) is released through 1% of the earth's surface area and takes the form of hydrothermal discharge from young (Pleistocene or younger) rocks adjacent to active seafloor-spreading centers and submarine volcanic areas. This amount is roughly equivalent to man's present gross energy consumption rate. A sub-seafloor geothermal reservoir, to be exploitable under future economic conditions, will have to be hot, porous, permeable, large, shallow, and near an energy-deficient, populated land mass. Furthermore, the energy must be recoverable using technology achievable at a competitive cost and numerous environmental, legal and institutional problems will have to be overcome. The highest-temperature reservoirs should be found adjacent to the zones of the seafloor extension or volcanism that are subject to high sedimentation rates. The relatively impermeable sediments reduce hydrothermal-discharge flow rates, forcing the heat to be either conducted away or released by high-temperature fluids, both of which lead to reservoir temperatures that can exceed 300??C. There is evidence that the oceanic crust is quite permeable and porous and that it was amenable to deep (3-5 km) penetration by seawater at least some time in the early stages of its evolution. Most of the heat escapes far from land, but there are notable exceptions. For example, in parts of the Gulf of California, thermal gradients in the bottom sediments exceed 1??C/m. In the coastal areas of the Gulf of California, where electricity and fresh water are at a premium, this potential resource lies in shallow water (< 200 m) and within sight of land. Other interesting areas include the Sea of Japan, the Sea of Okhotsk and the Andaman Sea along the margins of the western Pacific, the Tyrrhenian Sea west of Italy, and the southern California borderland and west flank of the Juan de Fuca Ridge off the west coast of the United States. Many questions remain to be answered about the physical and other characteristics of these systems before they can be considered a viable resource. Until several of the most promising areas are carefully defined and drilled, the problem will remain unresolved. ?? 1976.

Journal of Volcanology and Geothermal Research↗
USGS2015

Heat flux from magmatic hydrothermal systems related to availability of fluid recharge

Magmatic hydrothermal systems are of increasing interest as a renewable energy source. Surface heat flux indicates system resource potential, and can be inferred from soil CO 2 flux measurements and fumarole gas chemistry. Here we compile and reanalyze results from previous CO 2 flux surveys worldwide to compare heat flux from a variety of magma-hydrothermal areas. We infer that availability of water to recharge magmatic hydrothermal systems is correlated with heat flux. Recharge availability is in turn governed by permeability, structure, lithology, rainfall, topography, and perhaps unsurprisingly, proximity to a large supply of water such as the ocean. The relationship between recharge and heat flux interpreted by this study is consistent with recent numerical modeling that relates hydrothermal system heat output to rainfall catchment area. This result highlights the importance of recharge as a consideration when evaluating hydrothermal systems for electricity generation, and the utility of CO 2 flux as a resource evaluation tool.

Journal of Volcanology and Geothermal Research↗
USGS2015

Coupling geophysical investigation with hydrothermal modeling to constrain the enthalpy classification of a potential geothermal resource.

An appreciable challenge in volcanology and geothermal resource development is to understand the relationships between volcanic systems and low-enthalpy geothermal resources. The enthalpy of an undeveloped geothermal resource in the Karckar region of Armenia is investigated by coupling geophysical and hydrothermal modeling. The results of 3-dimensional inversion of gravity data provide key inputs into a hydrothermal circulation model of the system and associated hot springs, which is used to evaluate possible geothermal system configurations. Hydraulic and thermal properties are specified using maximum a priori estimates. Limited constraints provided by temperature data collected from an existing down-gradient borehole indicate that the geothermal system can most likely be classified as low-enthalpy and liquid dominated. We find the heat source for the system is likely cooling quartz monzonite intrusions in the shallow subsurface and that meteoric recharge in the pull-apart basin circulates to depth, rises along basin-bounding faults and discharges at the hot springs. While other combinations of subsurface properties and geothermal system configurations may fit the temperature distribution equally well, we demonstrate that the low-enthalpy system is reasonably explained based largely on interpretation of surface geophysical data and relatively simple models.

Karckar Region↗
USGS2025

Applying U-Th disequilbrium for dating siliceous sinters

Continental hydrothermal systems are critical avenues for the crustal transport of heat and mass captured for geothermal energy and mineral exploration. Thus, understanding their temporal evolution and longevity is important for resource characterization. Deposits of microlaminated siliceous sinter, common surface expressions of high temperature reservoirs (> 170 °C), have the potential to trace hydrothermal histories. Geothermal reservoirs are often located within uranium-bearing silicic volcanic rock where subsurface fluid-rock interactions extract U into hydrothermal fluids. U incorporated in the surface sinter deposit has the potential to provide a U— Th disequilibrium dating option. We focus on samples from El Tatio geyser field in the Altiplano of northern Chile, the largest geothermal system in the Andes. Our resulting 230 Th/U ages, along with the water and deposit elemental compositions, suggest concentrations of U and Th vary predictably along the sinter apron. While distal facies containing the highest U concentrations (> 50 μg/g) are least affected by detrital Th corrections, they can display suspected open-system behavior. In contrast, more medial facies, where bacterial mats and other porous textures are co mmonly concentrated, have only trace amounts of U (< 0.1 μg/g), which leads to unreliable or geologically improbable dates. Proximal facies tend to date most consistently. By comparing existing 14 C ages with 230 Th/U results, 230 Th/U ages tend to be younger than the 14 C ages, supporting the presence of a 14 C-dead carbon influence. New data confirm that the onset of geothermal activity at El Tatio goes back to the late Pleistocene.

Altiplano, El Tatio geyser field↗
USGS2016

The chemistry and isotopic composition of waters in the low-enthalpy geothermal system of Cimino-Vico Volcanic District, Italy

Geothermal energy exploration is based in part on interpretation of the chemistry, temperature, and discharge rate of thermal springs. Here we present the major element chemistry and the δD, δ 18 O, 87 Sr/ 86 Sr and δ 11 B isotopic ratio of groundwater from the low-enthalpy geothermal system near the city of Viterbo in the Cimino-Vico volcanic district of west-Central Italy. The geothermal system hosts many thermal springs and gas vents, but the resource is still unexploited. Water chemistry is controlled by mixing between low salinity,HCO 3 -rich fresh waters (<24.2°C) flowing in shallow volcanic rocks and SO 4 -rich thermal waters (25.3°C to 62.2°C) ascending from deep, high permeability Mesozoic limestones. The (equivalent) SO 4 /Cl (0.01–0.02), Na/Cl (2.82–5.83) and B/Cl ratios (0.02–0.38) of thermal waters differs from the ratios in other geothermal systems from Central Italy, probably implying a lack of hydraulic continuity across the region. The δ 18 O (−6.6‰ to −5.9‰) and δD (−40.60‰ to −36.30‰) isotopic composition of spring water suggest that the recharge area for the geothermal system is the summit region of Mount Cimino. The strontium isotope ratios ( 87 Sr/ 86 Sr) of thermal waters (0.70797–0.70805) are consistent with dissolution of the Mesozoic evaporite-carbonate units that constitute the reservoir, and the ratios of cold fresh waters mainly reflect shallow circulation through the volcanic cover and some minor admixture (<10%) of thermal water as well. The boron isotopic composition (δ 11 B) of fresh waters (−5.00 and 6.12‰) is similar to that of the volcanic cover, but the δ 11 B of thermal waters (−8.37‰ to −4.12‰) is a mismatch for the Mesozoic reservoir rocks and instead reflects dissolution of secondary boron minerals during fluid ascent through flysch units that overlie the reservoir. A slow and tortuous ascent enhances extraction of boron but also promotes conductive cooling, partially masking the heat present in the reservoir. Overall data from this study is consistent with previous studies that concluded that the geothermal system has a large energy potential.

Journal of Volcanology and Geothermal Research↗
USGS2010

Source and fate of inorganic solutes in the Gibbon River, Yellowstone National Park, Wyoming, USA: I. Low-flow discharge and major solute chemistry

The Gibbon River in Yellowstone National Park (YNP) is an important natural resource and habitat for fisheries and wildlife. However, the Gibbon River differs from most other mountain rivers because its chemistry is affected by several geothermal sources including Norris Geyser Basin, Chocolate Pots, Gibbon Geyser Basin, Beryl Spring, and Terrace Spring. Norris Geyser Basin is one of the most dynamic geothermal areas in YNP, and the water discharging from Norris is much more acidic (pH 3) than other geothermal basins in the upper-Madison drainage (Gibbon and Firehole Rivers). Water samples and discharge data were obtained from the Gibbon River and its major tributaries near Norris Geyser Basin under the low-flow conditions of September 2006. Surface inflows from Norris Geyser Basin were sampled to identify point sources and to quantify solute loading to the Gibbon River. The source and fate of the major solutes (Ca, Mg, Na, K, SiO 2 , Cl, F, HCO 3 , SO 4 , NO 3 , and NH 4 ) in the Gibbon River were determined in this study and these results may provide an important link in understanding the health of the ecosystem and the behavior of many trace solutes. Norris Geyser Basin is the primary source of Na, K, Cl, SO 4 , and N loads (35–58%) in the Gibbon River. The largest source of HCO 3 and F is in the lower Gibbon River reach. Most of the Ca and Mg originate in the Gibbon River upstream from Norris Geyser Basin. All the major solutes behave conservatively except for NH 4 , which decreased substantially downstream from Gibbon Geyser Basin, and SiO 2 , small amounts of which precipitated on mixing of thermal drainage with the river. As much as 9–14% of the river discharge at the gage is from thermal flows during this period.

Wyoming↗
USGS2022

Revisiting the depth distribution of seismicity before and after the 2004–2008 eruption of Mount St. Helens

Hypocenter estimation at active volcanoes improves our understanding of their magmatic systems and indicates changing conditions at depth for continuously monitored volcanoes. The most active volcano in the Cascades Range, Mount St. Helens, has a multi-decadal earthquake catalog and it shows an apparent change in the depth distribution of seismicity before and after the 2004–2008 dome-building eruption and unrest sequence. We use two new resources to evaluate the accuracy of hypocenters and consequently the change in depth distribution of seismicity before and after the 2004–2008 eruption. First, we deployed a dense array of 136 three-component nodal seismographs for one month in 2017, including sub-arrays on the newly extruded dome and crater floor. Second, for events recorded during this month, we located their hypocenters using a three-dimensional (3D) wavefront-tracking location solver and a recently developed tomography model derived from active and passive source data. The relocated hypocenters are generally shallower and more concentrated beneath the crater compared to their catalog locations. The mean hypocenter movement from catalog locations is 2.86 km, with an averaged depth shift of 2.53 km upward. Comparison between 2017 hypocenters located using all of the available phase picks and those located using only the catalog picks from the permanent network suggests the improved hypocenters mostly resulted from the location solver and the 3D velocity model, with smaller changes due to the dense three-component array. Applying the same hypocenter estimation method to phase picks for all events from 1997 to 2004 and 2008–2021, we found a concentration of seismicity between sea level and ~ 1 km above it before and after the 2004–2008 eruption. The new results suggest that the seismogenic structure in the shallow magmatic system quickly re-equilibrated to its earlier state after the dome-building eruption.

Washington↗