Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Volcanology and Geothermal Research”

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

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

At least 55 records · Page 3Linked to original sources

Lithic breccia and ignimbrite erupted during the collapse of Crater Lake Caldera, Oregon

The climactic eruption of Mount Mazama (6845 y.B.P.) vented a total of ∼50 km 3 of compositionally zoned rhyodacitic to basaltic magma from: (a) a single vent as a Plinian pumice fall deposit and the overlying Wineglass Welded Tuff, and (b) ring vents as ignimbrite and coignimbrite lithic breccia accompanying the collapse of Crater Lake caldera. New field and grain-size data for the ring-vent products are presented in this report. The coarse-grained, poorly bedded, clast-supported lithic breccia extends as far as 18 km from the caldera center. Like the associated ignimbrite, the breccia is compositionally zoned both radially and vertically, and silicic, mixed, and mafic types can be recognized, based on the proportion of rhyodacitic pumice. Matrix fractions in silicic breccias are depleted of fines and are lithic- and crystal-enriched relative to silicic ignimbrite due to vigorous gas sorting during emplacement. Ignimbrite occurs as a proximal veneer deposit overlying the breccia, a medial (∼ 8 to ∼ 25 km from the caldera center), compositionally zoned valley fill as much as > 110 m thick, and an unzoned distal (⪖ 20 km) facies which extends as far as 55 km from the caldera. Breccia within ∼ 9 km of the caldera center is interpreted as a coignimbrite lag breccia formed within the deflation zone of the collapsing ring-vent eruption columns. Expanded pyroclastic flows of the deflation zone were probably vertically graded in both size and concentration of blocks, as recently postulated for some turbidity currents. An inflection in the rate of falloff of lithic-clast size within the lithic breccia at ∼ 9 km may mark the outer edge of the deflation zone or may be an artifact of incomplete exposure. The onset of ring-vent activity at Mt. Mazama was accompanied by a marked increase in eruptive discharge. Pyroclastic flows were emplaced as a semicontinuous stream, as few ignimbrite flow-unit boundaries are evident. As eruption from the ring vents progressed, flow-runout distance and the extent of breccia deposition decreased due to (a) greater internal flow friction, and (b) decreasing eruption column heights. Effect (b) probably resulted from a progressive decrease in magmatic gas content and discharge rate. Waning discharge may have been promoted by the tapping of more viscous, crystal-rich magma, collapse of conduit walls, and declining caldera collapse rate.

Oregon↗

Quaternary silicic pyroclastic deposits of Atitlán Caldera, Guatemala

Atitlán caldera has been the site of several silicic eruptions within the last 150,000 years, following a period of basalt/andesite volcanism. The silicic volcanism began with 5–10 km 3 of rhyodacites, erupted as plinian fall and pyroclastic flows, about 126,000 yr. B.P. At 85,000 yr. B.P. 270–280 km 3 of compositionally distinct rhyolite was erupted in the Los Chocoyos event which produced widely dispersed, plinian fall deposits and widespread, mobile pyroclastic flows. In the latter parts of this eruption rhyodacite and minor dacite were erupted which compositionally resembled the earliest silicic magmas of the Atitlán center. As a result of this major eruption, the modern Atitlán (III) caldera formed. Following this event, rhyodacites were again erupted in smaller (5–13 km 3 ) volumes, partly through the lake, and mafic volcanism resumed, forming three composite volcanoes within the caldera. The bimodal mafic/silicic Atitlán volcanism is similar to that which has occurred elsewhere in the Guatemalan Highlands, but is significantly more voluminous. Mafic lavas are thought to originate in the mantle, but rise, intrude and underplate the lower crust and partly escape to the surface. Eventually, silicic melts form in the crust, possibly partly derived from underplated basaltic material, rise, crystallize and erupt. The renewed mafic volcanism could reflect either regional magmato-tectonic adjustment after the large silicic eruption or the onset of a new cycle.

Journal of Volcanology and Geothermal Research↗

Recent geologic history of lake Atitlán, a caldera lake in western Guatemala

Lake Atitlán, a caldera lake in western Guatemala, was investigated for evidence of recent volcanic and tectonic activity. No vents, faults, or folds are apparent on high-resolution seismic reflection profiles of lake sediment, representing at least 17,500 years and probably more than 35,000 years of deposition. Three post-caldera stratovolcanoes (San Pedro, Tolimán, and Atitlán) have grown in southern parts of the caldera and two deltas have grown from the north shore of the lake. Elsewhere, the caldera is occupied by Lake Atitlán, which is more than 300 m deep and has a relatively flat floor. Refraction profiles suggest that the original floor of the caldera lies ca. 300 m below the current lake floor, but prodigious amounts of methane gas in the lake sediment attenuated seismic signals and prevented any detailed view of the original caldera floor or faults along which the floor is presumed to have collapsed. Piston cores from deep basins of Lake Atitlán record ca. 2,000 years of unusually rapid sedimentation (ca. 0.5 cm/yr). Only one thin silicic ash layer was penetrated, and it is probably from a distant source. Cores contain evidence of Mayan disturbance of the environment around the lake and, to a lesser degree, of the lake itself; they might also record episodes of increased thermal activity, each lasting several decades. Heat-flow measurements inside and just outside the caldera are high (290 and 230 mW m −2 ), suggesting hydrothermal convection and a shallow heat source. High heat flow, a geological record of post-caldera silicic eruptions, and unexplained fluctuations of lake level (episodic tumescence ofthe lake floor?) suggest that magma remains beneath Lake Atitlán and that future eruptions are possible.

Journal of Volcanology and Geothermal Research↗

Subsidence and volcanism of the Haleakala Ridge, Hawaii

Side-looking sonar (GLORIA) mapping has revealed a series of four arcuate bands of high sonic backscatter on the crest of the Haleakala Ridge, a major rift-zone ridge extending 135 km east of the island of Maui. Dredge recovery indicates that the shallowest of these bands is a drowned coral reef, and the deeper bands are also inferred to be coral reefs. The reefs occur above a prominent submarine bench 1500–2500 m deep on the ridge (H-terrace) that marks the shoreline at the end of vigorous shield building of Haleakala volcano when lava flows ceased crossing and reworking the shoreline. Since their growth these reefs have subsided as much as 2200 m and have tilted systematically about 20 m/km southward as a result of post-reef volcanic loading on the island of Hawaii, whose center of mass is about directly south of the Haleakala Ridge. The 234 U/ 238 U age of the dredged coral is 750 ± 13 ka "> 750 ± 13ka , in reasonable agreement with an age of 850 ka for the underlying H terrace previously estimated from its relationship to other dated reefs to the southwest. Basalt glass fragments dredged from the Haleakala Ridge below the H terrace are tholeiitic and contain high sulfur indicative of eruption in water deeper than 200 m. Basalt glass fragments associated with the reefs above the H terrace are dominantly tholeiitic and contain intermediate sulfur contents, indicative of subaqueous eruption in shallow, near-shore conditions. One alkalic glass fragment was recovered above the H terrace. These relations indicate that the morphologic end of shield building as recorded by construction of the H terrace was not accompanyed by a change from tholeiitic to alkalic basalt; instead tholeiite eruptions continued for some time before the erupted lava became alkalic.

Hawaii↗

Fluid inclusions and preliminary studies of hydrothermal alteration in core hole PLTG-1, Platanares geothermal area, Honduras

The Platanares geothermal area in western Honduras consists of more than 100 hot springs that issue from numerous hot-spring groups along the banks or within the streambed of the Quebrada de Agua Caliente (brook of hot water). Evaluation of this geothermal area included drilling a 650-m deep PLTG-1 drill hole which penetrated a surface mantling of stream terrace deposits, about 550 m of Tertiary andesitic lava flows, and Cretaceous to lower Tertiary sedimentary rocks in the lower 90 m of the drill core. Fractures and cavities in the drill core are partly to completely filled by hydrothermal minerals that include quartz, kaolinite, mixed-layer illite-smectite, barite, fluorite, chlorite, calcite, laumontite, biotite, hematite, marcasite, pyrite, arsenopyrite, stibnite, and sphalerite; the most common open-space fillings are calcite and quartz. Biotite from 138.9-m depth, dated at 37.41 Ma by replicate 40Ar/39 Ar analyses using a continuous laser system, is the earliest hydrothermal mineral deposited in the PLTG-1 drill core. This mid-Tertiary age indicates that at least some of the hydrothermal alteration encountered in the PLTG-1 drill core occurred in the distant past and is unrelated to the present geothermal system. Furthermore, homogenization temperatures (Th) and melting-point temperatures (Tm) for fluid inclusions in two of the later-formed hydrothermal minerals, calcite and barite, suggest that the temperatures and concentration of dissolved solids of the fluids present at the time these fluid inclusions formed were very different from the present temperatures and fluid chemistry measured in the drill hole. Liquid-rich secondary fluid inclusions in barite and caicite from drill hole PLTG-1 have Th values that range from about 20??C less than the present measured temperature curve at 590.1-m depth to as much as 90??C higher than the temperature curve at 46.75-m depth. Many of the barite Th measurements (ranging between 114?? and 265??C) plot above the reference surface boiling-point curve for pure water assuming hydrostatic conditions; however, the absence of evidence for boiling in the fluid inclusions indicates that at the time the minerals formed, the ground surface must have been at least 80 m higher than at present and underwent stream erosion to the current elevation. Near-surface mixed-layer illite-smectite is closely associated with barite and appears to have formed at about the same temperature range (about 120?? to 200??C) as the fluid-inclusion Th values for barite. Fluid-inclusion Th values for calcite range between about 136?? and 213??C. Several of the calcite Th values are significantly lower than the present measured temperature curve. The melting-point temperatures (Tm) of fluid-inclusion ice yield calculated salinities, ranging from near zero to as much as 5.4 wt. % NaCl equivalent, which suggest that much of the barite and calcite precipitated from fluids of significantly greater salinity than the present low salinity Platanares hot-spring water or water produced from the drill hole.

Journal of Volcanology and Geothermal Research↗

Degassing and differentiation in subglacial volcanoes, Iceland

Within the neovolcanic zones of Iceland many volcanoes grew upward through icecaps that have subsequently melted. These steep-walled and flat-topped basaltic subglacial volcanoes, called tuyas, are composed of a lower sequence of subaqueously erupted, pillowed lavas overlain by breccias and hyaloclastites produced by phreatomagmatic explosions in shallow water, capped by a subaerially erupted lava plateau. Glass and whole-rock analyses of samples collected from six tuyas indicate systematic variations in major elements showing that the individual volcanoes are monogenetic, and that commonly the tholeiitic magmas differentiated and became more evolved through the course of the eruption that built the tuya. At Herdubreid, the most extensively studies tuya, the upward change in composition indicates that more than 50 wt.% of the first erupted lavas need crystallize over a range of 60°C to produce the last erupted lavas. The S content of glass commonly decreases upward in the tuyas from an average of about 0.08 wt.% at the base to < 0.02 wt.% in the subaerially erupted lava at the top, and is a measure of the depth of water (or ice) above the eruptive vent. The extensive subsurface crystallization that generates the more evolved, lower-temperature melts during the growth of the tuyas, apparently results from cooling and degassing of magma contained in shallow magma chambers and feeders beneath the volcanoes. Cooling may result from percolation of meltwater down cracks, vaporization, and cycling in a hydrothermal circulation. Degassing occurs when progressively lower pressure eruption (as the volcanic vent grows above the ice/water surface) lowers the volatile vapour pressure of subsurface melt, thus elevating the temperature of the liquidus and hastening liquid-crystal differentiation.

Journal of Volcanology and Geothermal Research↗

Deformation of the Wineglass Welded Tuff and the timing of caldera collapse at Crater Lake, Oregon

Four types of deformation occur in the Wineglass Welded Tuff on the northeast caldera rim of Crater Lake: (a) vertical tension fractures; (b) ooze-outs of fiamme: (c) squeeze-outs of fiamme; and (d) horizontal pull-apart structures. The three types of plastic deformation (b-d) developed in the lower part of the Wineglass Welded Tuff where degree of welding and density are maximum. Deformation originated from concentric normal faulting and landsliding as the caldera collapsed. The degree of deformation of the Wineglass Welded Tuff increases toward the northeast part of the caldera, where plastic deformation occurred more easily because of a higher emplacement temperature probably due to proximity to the vent. The probable glass transition temperature of the Wineglass Welded Tuff suggests that its emplacement temperature was ???750??C where the tuff is densely welded. Calculation of the conductive cooling history of the Wineglass Welded Tuff and the preclimactic Cleetwood (lava) flow under assumptions of a initially isothermal sheet and uniform properties suggests that (a) caldera collapse occurred a maximum of 9 days after emplacement of the Wineglass Welded Tuff, and that (b) the period between effusion of the Cleetwood (lava) flow and onset of the climactic eruption was <100 years. If cooling is controlled more by precipitation during quiescent periods than by conduction, these intervals must be shorter than the calculated times.

Oregon↗

Surface degassing and modifications to vesicle size distributions in active basalt flows

The character of the vesicle population in lava flows includes several measurable parameters that may provide important constraints on lava flow dynamics and rheology. Interpretation of vesicle size distributions (VSDs), however, requires an understanding of vesiculation processes in feeder conduits, and of post-eruption modifications to VSDs during transport and emplacement. To this end we collected samples from active basalt flows at Kilauea Volcano: (1) near the effusive Kupaianaha vent; (2) through skylights in the approximately isothermal Wahaula and Kamoamoa tube systems transporting lava to the coast; (3) from surface breakouts at different locations along the lava tubes; and (4) from different locations in a single breakout from a lava tube 1 km from the 51 vent at Pu'u ‘O’o. Near-vent samples are characterized by VSDs that show exponentially decreasing numbers of vesicles with increasing vesicle size. These size distributions suggest that nucleation and growth of bubbles were continuous during ascent in the conduit, with minor associated bubble coalescence resulting from differential bubble rise. The entire vesicle population can be attributed to shallow exsolution of H 2 O-dominated gases at rates consistent with those predicted by simple diffusion models. Measurements of H 2 O, CO 2 and S in the matrix glass show that the melt equilibrated rapidly at atmospheric pressure. Down-tube samples maintain similar VSD forms but show a progressive decrease in both overall vesicularity and mean vesicle size. We attribute this change to open system, “passive” rise and escape of larger bubbles to the surface. Such gas loss from the tube system results in the output of 1.2 × 10 6 g/day SO 2 , an output representing an addition of approximately 1% to overall volatile budget calculations. A steady increase in bubble number density with downstream distance is best explained by continued bubble nucleation at rates of 7–8/cm 3 s. Rates are ∼25% of those estimated from the vent samples, and thus represent volatile supersaturations considerably less than those of the conduit. We note also that the small total volume represented by this new bubble population does not: (1) measurably deplete the melt in volatiles; or (2) make up for the overall vesicularity decrease resulting from the loss of larger bubbles. Surface breakout samples have distinctive VSDs characterized by an extreme depletion in the small vesicle population. This results in samples with much lower number densities and larger mean vesicle sizes than corresponding tube samples. Similar VSD patterns have been observed in solidified lava flows and are interpreted to result from either static (wall rupture) or dynamic (bubble rise and capture) coalescence. Through comparison with vent and tube vesicle populations, we suggest that, in addition to coalescence, the observed vesicle populations in the breakout samples have experienced a rapid loss of small vesicles consistent with ‘ripening’ of the VSD resulting from interbubble diffusion of volatiles. Confinement of ripening features to surface flows suggests that the thin skin that forms on surface breakouts may play a role in the observed VSD modification.

Journal of Volcanology and Geothermal Research↗

Mechanical modeling of circumferential and radial dike intrusion on Galapagos volcanoes

A distinctive and unusual pattern of eruptive fissures is observed on the active volcanoes of the Galapagos islands, reflecting circumferential dike intrusion near the calderas and radial dike intrusion on the volcano flanks. Elastic finite-element models were used to investigate how a stress field could be produced and maintained to promote both circumferential and radial dike emplacement. Modeling results show that magma reservoirs of Galapagos volcanoes are probably diapiric, because this shape promotes both circumferential and radial intrusions, but magma pressure alone cannot create the observed pattern of dikes. Loading by volcano growth and magma reservoir pressure could produce a stress field of suitable orientation but insufficient magnitude. The intrusion of circumferential dikes could alter the stress field in a way that promotes future radial diking, and vice versa. Faulting or slumping within the calderas or on the volcano flanks in response to repeated intrusions could also create a stress field conducive to continued intrusion.

Galápagos Islands, Fernandina, Isabela↗

The Uwekahuna Ash Member of the Puna Basalt: product of violent phreatomagmatic eruptions at Kilauea volcano, Hawaii, between 2800 and 2100 14C years ago

Kilauea volcano's reputation for relatively gentle effusive eruptions belies a violent geologic past, including several large phreatic and phreatomagmatic eruptions that are recorded by Holocene pyroclastic deposits which mantle Kilauea's summit area and the southeast flank of adjacent Mauna Loa volcano. The most widespread of these deposits is the Uwekahuna Ash Member, a basaltic surge and fall deposit emplaced during two or more eruptive episodes separated by a few decades to several centuries. It is infered that the eruptions which produced the Uwekahuna were driven by water interacting with a fluctuating magma column. The volume, extent and character of the Uwekahuna deposits underscore the hazards posed by relatively infrequent but potentially devastating explosive eruptions at Kilauea, as well as at other basaltic volcanoes. -from Authors

Journal of Volcanology and Geothermal Research↗

Can magma-injection and groundwater forces cause massive landslides on Hawaiian volcanoes?

Landslides with volumes exceeding 1000 km 3 have occurred on the flanks of Hawaiian volcanoes. Because the flanks typically slope seaward no more than 12 °, the mechanics of slope failure are problematic. Limit-equilibrium analyses of wedge-shaped slices of the volcano flanks show that magma injection at prospective headscarps might trigger the landslides, but only under very restrictive conditions. If static magma weight is the sole source of magma pressure, hypothetical flank failures can have any size, but can occur only if slip-surface friction angles are less than about 16 °. If slip surfaces have friction angles more typical of fragmented or intact rocks (30–40 °), flank failures can occur only if mean magma pressures exceed static equilibrium pressures. Landslide length then scales with the excess magma pressure divided by the buoyant unit weight of the volcano flank. For typical excess magma pressures, buoyant unit weights and rock friction angles, the largest landslides that might be triggered in this manner have lengths of only several kilometers. This is at least an order of magnitude too small to explain the occurrence of giant Hawaiian landslides. The growing mass of active Hawaiian volcanoes can compress the edifice and substrate rocks and consequently produce groundwater head gradients that might destabilize larger sectors of the volcano flanks. However, calculations show that volcano growth at an estimated long-term vertical rate of 0.02 m/yr can generate significant head gradients only if an areally extensive, buried clay layer exists that has a great thickness (~ 200 m) or very low hydraulic diffusivity (~ 10 −11 m 2 /s). Additional calculations show that groundwater head gradients associated with topographically induced flow and sea-level change are less likely to be important. Thus a simple, quantitative explanation for failure of Hawaiian volcano flanks remains elusive, and more complex scenarios may merit investigation.

Journal of Volcanology and Geothermal Research↗

Inclusions in Mount St. Helens dacite erupted from 1980 through 1983

Inclusions of plutonic, metavolcanic and volcanic rocks are abundant in dacite pumice and lava from the 1980–1986 eruption sequence at Mount St. Helens. Point counts of inclusions exposed in talus blocks from the dome from 1980 through 1983 show that inclusions form approximately 3.5 vol% of the lava. Eighty-five percent of the inclusions are medium-grained gabbros with an average diameter of 6 cm. Additional rock types include quartz diorite, hornfelsic basalt, dacite, andesite and vein quartz. Disaggregated inclusions are common and define shear planes within the dome. These fragmented inclusions may significantly contaminate analyses of the dacite. The gabbroic inclusions are of four distinct types, all with mineral assemblages consistent with crystallization pressures of less than 9 kb. Textures and major-element compositions indicate that most of the gabbroic inclusions are cumulates. The most abundant inclusion type is laminated gabbronorite, which contains up to 9% interstitial glass, derived from partial melting. The presence of quartz veins and hornblende-bearing veins within sheared zones in the laminated gabbronorite indicates that the source of these inclusions was holocrystalline rock that had been penetrated by water-rich fluids. The gabbronorite contained sufficient water to be susceptible to partial melting when the magma that fed the 1980–1986 eruption sequence was emplaced nearby. Various types of gabbroic inclusions, including the laminated gabbronorite, are common in Mount St. Helens lavas of approximately the last 3000 years. This coincides with the interval in which Mount St. Helens first erupted basalt and basaltic andesite lavas. These observations, together with the fact that the gabbroic inclusions are compositionally unlike any of the Tertiary intrusive rocks in the Mount St. Helens area, strongly suggest that the inclusions are related to the introduction of basalt to the Mount St. Helens magmatic system. The source of the gabbros could be a layered mafic pluton formed through crystal accumulation from multiple batches of basaltic magma emplaced at mid-to upper-crustal depths beneath the volcano. The prevalence of explosive eruptions at Mount St. Helens may play a part in bringing the inclusions to the surface. The eruptive products of the cataclysmic eruption of May 18,1980 contain notably fewer inclusions than the pyroclastic flows and dome lavas erupted subsequently. This suggests that the May 18 eruption shattered conduit wall rock that was subsequently stoped into the magma and carried to the surface later in the eruption series.

Washington↗

Graphical and PC-software analysis of volcano eruption precursors according to the Materials Failure Forecast Method (FFM)

The Materials Failure Forecasting Method for volcanic eruptions (FFM) analyses the rate of precursory phenomena. Time of eruption onset is derived from the time of “failure” implied by accelerating rate of deformation. The approach attempts to fit data, Ω, to the differential relationship Ω &#x3A9; &#xA8; = A &#x3A9; &#x2D9; "> ¨ = AΩ˙ , where the dot superscript represents the time derivative, and the data Ω may be any of several parameters describing the accelerating deformation or energy release of the volcanic system. Rate coefficients, A and α, may be derived from appropriate data sets to provide an estimate of time to “failure”. As the method is still an experimental technique, it should be used with appropriate judgment during times of volcanic crisis. Limitations of the approach are identified and discussed. Several kinds of eruption precursory phenomena, all simulating accelerating creep during the mechanical deformation of the system, can be used with FFM. Among these are tilt data, slope-distance measurements, crater fault movements and seismicity. The use of seismic coda, seismic amplitude-derived energy release and time-integrated amplitudes or coda lengths are examined. Usage of cumulative coda length directly has some practical advantages to more rigorously derived parameters, and RSAM and SSAM technologies appear to be well suited to real-time applications. One graphical and four numerical techniques of applying FFM are discussed. The graphical technique is based on an inverse representation of rate versus time. For α = 2, the inverse rate plot is linear; it is concave upward for α < 2 and concave downward for α > 2. The eruption time is found by simple extrapolation of the data set toward the time axis. Three numerical techniques are based on linear least-squares fits to linearized data sets. The “linearized least-squares technique” is most robust and is expected to be the most practical numerical technique. This technique is based on an iterative linearization of the given rate-time series. The hindsight technique is disadvantaged by a bias favouring a too early eruption time in foresight applications. The “log rate versus log acceleration technique”, utilizing a logarithmic representation of the fundamental differential equation, is disadvantaged by large data scatter after interpolation of accelerations. One further numerical technique, a nonlinear least-squares fit to rate data, requires special and more complex software. PC-oriented computer codes were developed for data manipulation, application of the three linearizing numerical methods, and curve fitting. Separate software is required for graphing purposes. All three linearizing techniques facilitate an eruption window based on a data envelope according to the linear least-squares fit, at a specific level of confidence, and an estimated rate at time of failure.

Journal of Volcanology and Geothermal Research↗

The chemical and hydrologic structure of Poas volcano, Costa Rica

Comparison of the chemical characteristics of spring and river water draining the flanks of Poas Volcano, Costa Rica indicates that acid chloride sulfate springs of the northwestern flank of the volcano are derived by leakage and mixing of acid brines formed in the summit hydrothermal system with dilute flank groundwater. Acid chloride sulfate waters of the Rio Agrio drainage basin on the northwestern flank are the only waters on Poas that are affected by leakage of acid brines from the summit hydrothermal system. Acid sulfate waters found on the northwestern flank are produced by the interaction of surface and shallow groundwater with dry and wet acid deposition of SO2 and H2SO4 aerosols, respectively. The acid deposition is caused by a plume of acid gases that is released by a shallow magma body located beneath the active crater of Poas. -from Authors

Journal of Volcanology and Geothermal Research↗

Groundwater transport of crater-lake brine at Poa´s Volcano, Costa Rica

Poa´s Volcano is an active stratovolcano in Costa Rica that has a lake in its active crater. The crater lake has high temperatures (50–90 °C), high acidity (pH ≈ 0.0), and a high dissolved-solids content (100 g/kg). The volcano has numerous freshwater springs on its flanks, but a few on the northwestern flank are highly acidic (pH = 1.6–2.5) and have high dissolved-solids concentrations (2–22 g/kg). This study analyzes the regional groundwater system at Poa´s and demonstrates the likelihood that the water discharging from the acidic springs in the Rio Agrio watershed originates at the acidic crater lake. Both heat and solute transport are analyzed on a regional scale through numerical simulations using the HST3D finite-difference model, which solves the coupled equations for fluid flow, heat transport, and solute transport. The code allows fluid viscosity and density to be functions of both temperature and solute concentration. The simulations use estimates for recharge to the mountain and a range of values and various distributions of permeability and porosity. Several sensitivity analyses are performed to test how the uncertainty in many of the model parameters affects the simulation results. These uncertainties yield an estimated range of travel times from the crater lake to the Rio Agrio springs of 1–30 years, which is in close agreement with the results of tritium analyses of the springs. Calculated groundwater fluxes into and out of the crater lake are both about several hundred kg/s. These fluxes must be accounted for in water budgets of the crater lake.

Poa's Volcano↗

Evolution of tholeiitic diabase sheet systems in the eastern United States: examples from the Culpeper Basin, Virginia-Maryland, and the Gettysburg Basin, Pennsylvania

High-TiO2, quartz-normative (HTQ) tholeiite sheets of Early Jurassic age have intruded mainly Late Triassic sedimentary rocks in several early Mesozoic basins in the eastern US. Field observations, petrographic study, geochemical analyses and stable isotope data from three HTQ sheet systems were used to develop a general model of magmatic differentiation and magmatic-hydrothermal interaction for HTQ sheets. The three sheet systems have remarkably similar major-oxide and trace-element compositions. Cumulus and evolved diabase in comagmatic sheets separated by tens of kilometers are related by igneous differentiation. Differentiated diabase in all three sheets have petrographic and geochemical signatures and fluid inclusions indicating hydrothermal alteration beginning near magmatic temperatures and continuing to relatively low temperatures. Sulfur and oxygen isotope data are consistent with a magmatic origin for the hydrothermal fluid. -from Authors

Pennsylvania, Maryland, Virginia↗

The Hengill geothermal area, Iceland: variation of temperature gradients deduced from the maximum depth of seismogenesis

Given a uniform lithology and strain rate and a full seismic data set, the maximum depth of earthquakes may be viewed to a first order as an isotherm. These conditions are approached at the Hengill geothermal area, S. Iceland, a dominantly basaltic area. The temperature at which seismic failure ceases for the strain rates likely at the Hengill geothermal area is determined by analogy with oceanic crust, and is about 650 ?? 50??C. The topographies of the top and bottom of the seismogenic layer were mapped using 617 earthquakes. The thickness of the seismogenic layer is roughly constant and about 3 km. A shallow, aseismic, low-velocity volume within the spreading plate boundary that crosses the area occurs above the top of the seismogenic layer and is interpreted as an isolated body of partial melt. The base of the seismogenic layer has a maximum depth of about 6.5 km beneath the spreading axis and deepens to about 7 km beneath a transform zone in the south of the area. -from Author

Journal of Volcanology and Geothermal Research↗

Comparison of hydrothermal alteration of carboniferous carbonate and siliclastic rocks in the Valles caldera with outcrops from the Socorro caldera, New Mexico

Continental Scientific Drilling Program (CSDP) drill hole VC-2B [total depth 1761.7 m (5780 ft); maximum temperature 295 °C] was continuously cored through the Sulphur Springs hydrothermal system in the western ring-fracture zone of the 1.14 Ma Valles caldera. Among other units, the hole penetrated 760.2 m (2494.1 ft) of Paleozoic carbonate and siliciclastic strata underlying caldera fill and precaldera volcanic and epiclastic rocks. Comparison of the VC-2B Paleozoic rocks with corresponding lithologies within and around the 32.1 Ma Socorro caldera, 192 km ( 119 miles) to the south-southwest, provides insight into the variability of alteration responses to similar caldera-related hydrothermal regimes. The Pennsylvanian Madera Limestone and Sandia Formation from VC-2B preserve many of the sedimentological and diagenetic features observed in these units on a regional basis and where unaffected by high temperatures or hydrothermal activity. Micrites in these formations in VC-2B are generally altered and mineralized only where fractured or brecciated, that is, where hydrothermal solutions could invade carbonate rocks which were otherwise essentially impermeable. Alteration intensity (and correspondingly inferred paleopermeability) is only slightly higher in carbonate packstones and grainstones, low to intermediate in siltstones and claystones, and high in poorly cemented sandstones. Hydrothermal fracture-filling phases in these rocks comprise sericite (and phengite), chlorite, allanite, apatite, an unidentified zeolite and sphene in various combinations, locally with sphalerite, galena, pyrite and chalcopyrite. Terrigenous feldspars and clays are commonly altered to chlorite and seriate, and euhedral anhydrite “porphyroblasts” with minor chlorite occur in Sandia Formation siltstone. Fossils are typically unaltered, but the walls of some colonial bryozoans in the Madera Limestone are altered to the assemblage chlorite-sericite-epidote-allanite. La, Ce and Nd are present in an unidentified hydrothermal mineral occurring throughout much of the VC-2B Pennsylvanian sequence. Carboniferous carbonate and siliciclastic formations within and around the Socorro caldera show a similar style of alteration and mineralization to their Valles caldera counterparts, but by contrast locally host commercial, caldera-related, base-metal sulfide deposits. As in the Valles rocks, mineralization and alteration in those of the Socorro caldera were strongly controlled by porosity. Unless disrupted by fractures, breccias, or karst cavities ( not identified in Valles caldera drill holes), the rocks remained relatively unaltered. Where these features allowed ingress of mineralizing hydrothermal solutions, base-metal sulfides and rare-earth-element-bearing minerals were precipitated.

New Mexico↗