Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Interpretation”

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 1,297 records · Page 72Linked to original sources
USGS2010

Seismicity and fluid geochemistry at Lassen Volcanic National Park, California: Evidence for two circulation cells in the hydrothermal system

Seismic analysis and geochemical interpretations provide evidence that two separate hydrothermal cells circulate within the greater Lassen hydrothermal system. One cell originates south to SW of Lassen Peak and within the Brokeoff Volcano depression where it forms a reservoir of hot fluid (235–270°C) that boils to feed steam to the high-temperature fumarolic areas, and has a plume of degassed reservoir liquid that flows southward to emerge at Growler and Morgan Hot Springs. The second cell originates SSE to SE of Lassen Peak and flows southeastward along inferred faults of the Walker Lane belt (WLB) where it forms a reservoir of hot fluid (220–240°C) that boils beneath Devils Kitchen and Boiling Springs Lake, and has an outflow plume of degassed liquid that boils again beneath Terminal Geyser. Three distinct seismogenic zones (identified as the West , Middle , and East seismic clusters) occur at shallow depths (<6 km) in Lassen Volcanic National Park, SW to SSE of Lassen Peak and adjacent to areas of high-temperature (≤161°C) fumarolic activity (Sulphur Works, Pilot Pinnacle, Little Hot Springs Valley, and Bumpass Hell) and an area of cold, weak gas emissions (Cold Boiling Lake). The three zones are located within the inferred Rockland caldera in response to interactions between deeply circulating meteoric water and hot brittle rock that overlies residual magma associated with the Lassen Volcanic Center. Earthquake focal mechanisms and stress inversions indicate primarily N–S oriented normal faulting and E–W extension, with some oblique faulting and right lateral shear in the East cluster. The different focal mechanisms as well as spatial and temporal earthquake patterns for the East cluster indicate a greater influence by regional tectonics and inferred faults within the WLB. A fourth, deeper (5–10 km) seismogenic zone (the Devils Kitchen seismic cluster) occurs SE of the East cluster and trends NNW from Sifford Mountain toward the Devils Kitchen thermal area where fumarolic temperatures are ≤123°C. Lassen fumaroles discharge geothermal gases that indicate mixing between a N 2 -rich, arc-type component and gases derived from air-saturated meteoric recharge water. Most gases have relatively weak isotopic indicators of upper mantle or volcanic components, except for gas from Sulphur Works where δ 13 C–CO 2 , δ 34 S–H 2 S, and δ 15 N–N 2 values indicate a contribution from the mantle and a subducted sediment source in an arc volcanic setting.

California↗
USGS2000

Sea-floor environments within Long Island Sound: A regional overview

Modern sea-floor sedimentary environments within the glaciated, topographically complex Long Island Sound estuary have been interpreted and mapped from an extensive collection of sidescan sonographs, bottom samples, and video-camera observations together with supplemental bathymetric, marine-geologic, and bottom-current data. Four categories of environments are present that reflect the dominant long-term processes of erosion or nondeposition; coarsegrained bedload transport; sediment sorting and reworking; and fine-grained deposition. (1) Environments of erosion or nondeposition contain exposures of glacial drift, coarse lag deposits, and possibly bedrock and include sediments which range from boulder fields to gravelly coarse-to-medium sands. (2) Environments of coarse-grained bedload transport are mantled by sand ribbons and sand waves and contain mostly coarse-to-fine sands with only small amounts of mud. (3) Environments of sediment sorting and reworking comprise both uniform and heterogeneous sediment types and contain variable amounts of fine sand and mud. (4) Environments of fine-grained deposition are blanketed by muds and sandy muds. The patchy distribution of sedimentary environments within Long Island Sound reflects both regional and local changes in bottom processes. Regional changes are primarily the result of a strong, east-to-west decreasing gradient of bottom tidal-current speeds, coupled with the net (westward) estuarine bottom drift. The regional current regime has produced a westward succession of environments along the basin floor beginning with erosion or nondeposition at the narrow eastern entrance to the Sound, changing to an extensive area of coarse-grained bedload transport, passing into a contiguous band of sediment sorting, and ending with broad areas of fine-grained deposition in the central and western Sound. However, local changes in processes are superimposed on the regional conditions within the central and western parts of the basin and along the nearshore margins. Within the central and western basin, localized sedimentary environments are produced where the bottom flow is enhanced by, and interacts with, the bottom topography, whereas along the nearshore margins, they variously reflect wave-produced currents, the irregular bathymetry, the indented shoreline, and the proximal supply of sediments. Results from this study (1) confirm the high trapping efficiency of fine-grained sediments in the Sound, (2) suggest that fine-grained sediments accumulate at an average (regional) rate of 0.08 g/cm 2 /y , and (3) indicate that the postglacial delta in the eastern Sound was a significant source of fine-grained sediments now buried beneath depositional areas.

Connecticut, New York↗
USGS2008

Methanogenic pathways of coal-bed gas in the Powder River Basin, United States: The geologic factor

Coal-bed gas of the Tertiary Fort Union and Wasatch Formations in the Powder River Basin in Wyoming and Montana, U.S. was interpreted as microbial in origin by previous studies based on limited data on the gas and water composition and isotopes associated with the coal beds. To fully evaluate the microbial origin of the gas and mechanisms of methane generation, additional data for 165 gas and water samples from 7 different coal-bed methane-bearing coal-bed reservoirs were collected basinwide and correlated to the coal geology and stratigraphy. The C1/(C2 + C3) ratio and vitrinite reflectance of coal and organic shale permitted differentiation between microbial gas and transitional thermogenic gas in the central part of the basin. Analyses of methane ??13C and ??D, carbon dioxide ??13C, and water ??D values indicate gas was generated primarily from microbial CO2 reduction, but with significant gas generated by microbial methyl-type fermentation (aceticlastic) in some areas of the basin. Microbial CO2 reduction occurs basinwide, but is generally dominant in Paleocene Fort Union Formation coals in the central part of the basin, whereas microbial methyl-type fermentation is common along the northwest and east margins. Isotopically light methane ??13C is distributed along the basin margins where ??D is also depleted, indicating that both CO2-reduction and methyl-type fermentation pathways played major roles in gas generation, but gas from the latter pathway overprinted gas from the former pathway. More specifically, along the northwest basin margin gas generation by methyl-type fermentation may have been stimulated by late-stage infiltration of groundwater recharge from clinker areas, which flowed through highly fractured and faulted coal aquifers. Also, groundwater recharge controlled a change in gas composition in the shallow Eocene Wasatch Formation with the increase of nitrogen and decrease of methane composition of the coal-bed gas. Other geologic factors, such as burial, thermal and maturation history, lateral and vertical continuity, and coalification of the coal beds, also played a significant role in controlling methanogenic pathways and provided new perspectives on gas evolution and emplacement. The early-stage gas produced by CO2 reduction has mixed with transitional thermogenic gas in the deeper, central parts of the Powder River Basin to form 'old' gas, whereas along the basin margins the overprint of gas from methyl-type fermentation represents 'new' gas. Thus, a clear understanding of these geologic factors is necessary to relate the microbiological, biogeochemical, and hydrological processes involved in the generation of coal-bed gas.

International Journal of Coal Geology↗
USGS2017

Magmatic degassing, lava dome extrusion, and explosions from Mount Cleveland volcano, Alaska, 2011–2015: Insight into the continuous nature of volcanic activity over multi-year timescales

Mount Cleveland volcano (1730 m) is one of the most active volcanoes in the Aleutian arc, Alaska, but heightened activity is rarely accompanied by geophysical signals, which makes interpretation of the activity difficult. In this study, we combine volcanic gas emissions measured for the first time in August 2015 with longer-term measurements of thermal output and lava extrusion rates between 2011 and 2015 calculated from MODIS satellite data with the aim to develop a better understanding of the nature of volcanic activity at Mount Cleveland. Degassing measurements were made in the month following two explosive events (21 July and 7 August 2015) and during a period of new dome growth in the summit crater. SO 2 emission rates ranged from 400 to 860 t d − 1 and CO 2 /SO 2 ratios were < 3, consistent with the presence of shallow magma in the conduit and the observed growth of a new lava dome. Thermal anomalies derived from MODIS data from 2011 to 2015 had an average repose time of only 4 days, pointing to the continuous nature of volcanic activity at this volcano. Rapid increases in the cumulative thermal output were often coincident with visual confirmation of dome growth or accumulations of tephra in the crater. The average rate of lava extrusion calculated for 9 periods of rapid increase in thermal output was 0.28 m 3 s − 1 , and the total volume extruded from 2011 to 2015 was 1.9–5.8 Mm 3 . The thermal output from the lava extrusion events only accounts for roughly half of the thermal budget, suggesting a continued presence of shallow magma in the upper conduit, likely driven by convection. Axisymmetric dome morphology and occasional drain back of lava into the conduit suggests low-viscosity magmas drive volcanism at Mount Cleveland. It follows also that only small overpressures can be maintained given the small domes and fluid magmas, which is consistent with the low explosivity of most of Mount Cleveland's eruptions. Changes between phases of dome growth and explosive activity are somewhat unpredictable and likely result from plugs that are related to the dome obtaining a critical dimension, or from small variations in the magma ascent rate that lead to crystallization-induced blockages in the upper conduit, thereby reducing the ability of magma to degas. We suggest the small magma volumes, slow ascent rates, and low magma viscosity lead to the overall lack of anomalous geophysical signals prior to eruptions, and that more continuous volcanic degassing measurements might lead to more successful eruption forecasting at this continuously-active open-vent volcano.

Alaska↗
USGS2010

Depositional setting, petrology and chemistry of Permian coals from the Paraná Basin: 2. South Santa Catarina Coalfield, Brazil

In Brazil economically important coal deposits occur in the southern part of the Paraná Basin, where coal seams occur in the Permian Rio Bonito Formation, with major coal development in the states of Rio Grande de Sul and Santa Catarina. The current paper presents results on sequence stratigraphic interpretation of the coal-bearing strata , and petrological and geochemical coal seam characterization from the South Santa Catarina Coalfield, Paraná Basin. In terms of sequence stratigraphic interpretation the precursor mires of the Santa Catarina coal seams formed in an estuarine-barrier shoreface depositional environment , with major peat accumulation in a high stand systems tract (Pre-Bonito and Bonito seams), a lowstand systems tract (Ponta Alta seam, seam A, seam B) and a transgressive systems tract (Irapuá, Barro Branco and Treviso seams). Seam thicknesses range from 1.70 to 2.39 m, but high proportions of impure coal (coaly shale and shaley coal), carbonaceous shale and partings reduce the net coal thickness significantly. Coal lithoypes are variable, with banded coal predominant in the Barro Branco seam, and banded dull and dull coal predominantly in Bonito and Irapuá seams, respectively. Results from petrographic analyses indicate a vitrinite reflectance range from 0.76 to 1.63 %Rrandom (HVB A to LVB coal). Maceral group distribution varies significantly, with the Barro Branco seam having the highest vitrinite content (mean 67.5 vol%), whereas the Irapuá seam has the highest inertinite content (33.8 vol%). Liptinite mean values range from 7.8 vol% (Barro Branco seam) to 22.5 vol% (Irapuá seam). Results from proximate analyses indicate for the three seams high ash yields (50.2 – 64.2 wt.%). Considering the International Classification of in-Seam Coals, all samples are in fact classified as carbonaceous rocks (> 50 wt.% ash). Sulfur contents range from 3.4 to 7.7 wt.%, of which the major part occurs as pyritic sulfur. Results of X-ray diffraction indicate the predominance of quartz and kaolinite (also pyrite). Gypsum , gibbsite , jarosite and calcite were also identified in some samples. Feldspar was noted but is rare. The major element distribution in the three seams (coal basis) is dominated by SiO 2 (31.3 wt.%, mean value), Al 2 O 3 (14.5 wt.%, mean value) and Fe 2 O 3 (6.9 wt.%, mean value). Considering the concentrations of trace elements that are of potential environmental hazards the Barro Branco, Bonito and Irapuá seams (coal base) are significantly enriched in Co (15.7 ppm), Cr (54.5 ppm), Li (59.3 ppm), Mn (150.4 ppm), Pb (58.0 ppm) and V (99.6 ppm), when compared to average trace elements contents reported for U. S. coals. Hierarchical cluster analysis identified, based on similarity levels, three groups of major elements and seven groups of trace elements. Applying discriminant analyses using trace and major element distribution, it could be demonstrated that the three seams from Santa Catarina show distinct populations in the discriminant analyses plots, and also differ from the coals of Rio Grande do Sul analyzed in a previous study.

Santa Catarina↗
USGS2016

RMT focal plane sensitivity to seismic network geometry and faulting style

Modern tectonic studies often use regional moment tensors (RMTs) to interpret the seismotectonic framework of an earthquake or earthquake sequence; however, despite extensive use, little existing work addresses RMT parameter uncertainty. Here, we quantify how network geometry and faulting style affect RMT sensitivity. We examine how data-model fits change with fault plane geometry (strike and dip) for varying station configurations. We calculate the relative data fit for incrementally varying geometries about a best-fitting solution, applying our workflow to real and synthetic seismograms for both real and hypothetical station distributions and earthquakes. Initially, we conduct purely observational tests, computing RMTs from synthetic seismograms for hypothetical earthquakes and a series of well-behaved network geometries. We then incorporate real data and station distributions from the International Maule Aftershock Deployment (IMAD), which recorded aftershocks of the 2010 M W 8.8 Maule earthquake, and a set of regional stations capturing the ongoing earthquake sequence in Oklahoma and southern Kansas. We consider RMTs computed under three scenarios: (1) real seismic records selected for high data quality; (2) synthetic seismic records with noise computed for the observed source-station pairings and (3) synthetic seismic records with noise computed for all possible station-source pairings. To assess RMT sensitivity for each test, we observe the ‘fit falloff’, which portrays how relative fit changes when strike or dip varies incrementally; we then derive the ranges of acceptable strikes and dips by identifying the span of solutions with relative fits larger than 90 per cent of the best fit. For the azimuthally incomplete IMAD network, Scenario 3 best constrains fault geometry, with average ranges of 45° and 31° for strike and dip, respectively. In Oklahoma, Scenario 3 best constrains fault dip with an average range of 46°; however, strike is best constrained by Scenario 1, with a range of 26°. We draw two main conclusions from this study. (1) Station distribution impacts our ability to constrain RMTs using waveform time-series; however, in some tectonic settings, faulting style also plays a significant role and (2) increasing station density and data quantity (both the number of stations and the number of individual channels) does not necessarily improve RMT constraint. These results may be useful when organizing future seismic deployments (e.g. by concentrating stations in alignment with anticipated nodal planes), and in computing RMTs, either by guiding a more rigorous data selection process for input data or informing variable weighting among the selected data (e.g. by eliminating the transverse component when strike-slip mechanisms are expected).

Geophysical Journal International↗
USGS1989

Multiple hydrothermal and metamorphic events in the Kidd Creek volcanogenic massive sulphide deposit, Timmins, Ontario: evidence from tourmalines and chlorites

Tourmaline and chlorite are the principal ferromagnesian silicate minerals in the Kidd Creek massive sulphide deposit. Tourmaline is most common in sphalerite-rich peripheral margins of the chalcopyrite stringer zone. Within the north orebody, samples typically contain <1% tourmaline, but small areas (hand-specimen scale) may have 10–20%. Chlorite is more widely distributed and in places constitutes 30–50% of rock volumes. Associated assemblages may include quartz, sulphides (principally chalcopyrite, sphalerite, and (or) pyrite), carbonate, albite, sericite, and rare fluorite, allanite, or zoisite(?).The tourmalines and chlorites record a series of multiple hydrothermal and metamorphic events. Paragenetic studies suggest that tourmaline was deposited during several discrete stages of mineralization, as evidenced by brecciation and cross-cutting relationships. Most of the tourmalines have two concentric growth zones defined by different colours (green, brown, blue, yellow). Some tourmalines also display pale discordant rims that cross-cut and embay the inner growth zones and polycrystalline, multiple-extinction domains. Late sulphide veinlets (chalcopyrite, pyrrhotite) transect the inner growth zones and pale discordant rims of many crystals. The concentric growth zones are interpreted as primary features developed by the main ore-forming hydrothermal system, whereas the discordant rims, polycrystalline domains, and cross-cutting sulphide veinlets reflect post-ore metamorphic processes.Detailed electron microprobe analyses of tourmalines show a wide compositional range, from Fe-rich dravite nearly to end-member schorl, with Fe/(Fe + Mg) ratios varying from 0.33 to 0.92; only minor amounts of Ca are present, yielding uniformly high Na/(Na + Ca) ratios of 0.84–0.99. Two sets of chemical zoning trends are identified in the tourmalines, involving systematic changes in Fe/(Fe + Mg), Na/(Na + Ca), Al, and Ti that are believed to reflect internal coupled substitutions (e.g.,  + Ti = Na + Al) and local mineral equilibria (e.g., tourmaline–chlorite). Analyses of the pale discordant reaction rims show consistent depletion of Fe, Ca, and Ti, presumably by fluid–solid reactions during post-ore metamorphism.Chlorites also show an extensive range in composition, from ripidolite nearly to end-member daphnite, with Fe/(Fe + Mg) ratios of 0.43–0.98 and Si cation values of 5.00–5.39. Chlorites from the fringes of the footwall stringer zone have narrow compositional ranges, whereas chlorites near footwall rhyolite sills in the core of the stringer zone display major variations in Fe/(Fe + Mg) ratios, including one sample with a range of 0.68–0.95. The former group of chlorites has Fe/(Fe + Mg) ratios that correlate well with those of coexisting tourmalines (exclusive of late reaction rims). Data for the latter group, in contrast, fall off equilibrium K D curves, indicating that the tourmalines and chlorites within these samples are not in chemical equilibrium. The chlorites are believed to have been altered (overprinted) by Fe-rich hydrothermal fluids apparently generated during intrusion of the rhyolite sills. The tourmalines, however, are unaffected and retain primary chemical signatures.Variations in mineral proportions and mineral chemistry within the deposit mainly depend on fluctuations in temperature, pH, water/rock ratios, and amounts of entrained seawater. The major proposed control is mixing between high-temperature, Fe-rich end-member hydrothermal fluids and cold, Mg-rich entrained seawater. Fe/(Fe + Mg) variations in footwall tourmalines (and equilibrium chlorites) are believed to largely reflect the progressive infiltration of Mg-rich seawater into the margins and top of the hydrothermal system. The more Fe-rich compositions of Kidd Creek tourmalines relative to those from sediment-hosted massive sulphide deposits (e.g., Sullivan, British Columbia) may be related to the preferential generation of end-member hydrothermal fluids in proximal volcanic environments like that at Kidd Creek.

Ontario↗
USGS2013

Merapi 2010 eruption—Chronology and extrusion rates monitored with satellite radar and used in eruption forecasting

Despite dense cloud cover, satellite-borne commercial Synthetic Aperture Radar (SAR) enabled frequent monitoring of Merapi volcano's 2010 eruption. Near-real-time interpretation of images derived from the amplitude of the SAR signals and timely delivery of these interpretations to those responsible for warnings, allowed satellite remote sensing for the first time to play an equal role with in situ seismic, geodetic and gas monitoring in guiding life-saving decisions during a major volcanic crisis. Our remotely sensed data provide an observational chronology for the main phase of the 2010 eruption, which lasted 12 days (26 October–7 November, 2010). Unlike the prolonged low-rate and relatively low explosivity dome-forming and collapse eruptions of recent decades at Merapi, the eruption began with an explosive eruption that produced a new summit crater on 26 October and was accompanied by an ash column and pyroclastic flows that extended 8 km down the flanks. This initial explosive event was followed by smaller explosive eruptions on 29 October–1 November, then by a period of rapid dome growth on 1–4 November, which produced a summit lava dome with a volume of ~ 5 × 10 6 m 3 . A paroxysmal VEI 4 magmatic eruption (with ash column to 17 km altitude) destroyed this dome, greatly enlarged the new summit crater and produced extensive pyroclastic flows (to ~ 16 km radial distance in the Gendol drainage) and surges during the night of 4–5 November. The paroxysmal eruption was followed by a period of jetting of gas and tephra and by a second short period (12 h) of rapid dome growth on 6 November. The eruption ended with low-level ash and steam emissions that buried the 6 November dome with tephra and continued at low levels until seismicity decreased to background levels by about 23 November. Our near-real-time commercial SAR documented the explosive events on 26 October and 4–5 November and high rates of dome growth (> 25 m 3 s − 1 ). An event tree analysis for the previous 2006 Merapi eruption indicated that for lava dome extrusion rates > 1.2 m 3 s − 1 , the probability of a large (1872-scale) eruption was ~ 10%. Consequently, the order-of-magnitude greater rates in 2010, along with the explosive start of the eruption on 26 October, the large volume of lava accumulating at the summit by 4 November, and the rapid and large increases in seismic energy release, deformation and gas emissions were the basis for warnings of an unusually large eruption by the Indonesian Geological Agency's Center for Volcanology and Geologic Hazard Mitigation (CVGHM) and their Volcano Research and Technology Development Center (BPPTK) in Yogyakarta — warnings that saved thousands of lives.

Merapi Volcano↗
USGS2007

Fault locking, block rotation and crustal deformation in the Pacific Northwest

We interpret Global Positioning System (GPS) measurements in the northwestern United States and adjacent parts of western Canada to describe relative motions of crustal blocks, locking on faults and permanent deformation associated with convergence between the Juan de Fuca and North American plates. To estimate angular velocities of the oceanic Juan de Fuca and Explorer plates and several continental crustal blocks, we invert the GPS velocities together with seafloor spreading rates, earthquake slip vector azimuths and fault slip azimuths and rates. We also determine the degree to which faults are either creeping aseismically or, alternatively, locked on the block-bounding faults. The Cascadia subduction thrust is locked mainly offshore, except in central Oregon, where locking extends inland. Most of Oregon and southwest Washington rotate clockwise relative to North America at rates of 0.4–1.0 ° Myr −1 . No shear or extension along the Cascades volcanic arc has occurred at the mm/yr level during the past decade, suggesting that the shear deformation extending northward from the Walker Lane and eastern California shear zone south of Oregon is largely accommodated by block rotation in Oregon. The general agreement of vertical axis rotation rates derived from GPS velocities with those estimated from palaeomagnetic declination anomalies suggests that the rotations have been relatively steady for 10–15 Ma. Additional permanent dextral shear is indicated within the Oregon Coast Range near the coast. Block rotations in the Pacific Northwest do not result in net westward flux of crustal material—the crust is simply spinning and not escaping. On Vancouver Island, where the convergence obliquity is less than in Oregon and Washington, the contractional strain at the coast is more aligned with Juan de Fuca—North America motion. GPS velocities are fit significantly better when Vancouver Island and the southern Coast Mountains move relative to North America in a block-like fashion. The relative motions of the Oregon, western Washington and Vancouver Island crustal blocks indicate that the rate of permanent shortening, the type that causes upper plate earthquakes, across the Puget Sound region is 4.4 ± 0.3 mm yr −1 . This shortening is likely distributed over several faults but GPS data alone cannot determine the partitioning of slip on them. The transition from predominantly shear deformation within the continent south of the Mendocino Triple Junction to predominantly block rotations north of it is similar to changes in tectonic style at other transitions from shear to subduction. This similarity suggests that crustal block rotations are enhanced in the vicinity of subduction zones possibly due to lower resisting stress.

Geophysical Journal International↗
USGS2001

Anomalously high b-values in the South Flank of Kilauea volcano, Hawaii: Evidence for the distribution of magma below Kilauea's East rift zone

The pattern of b-value of the frequency-magnitude relation, or mean magnitude, varies little in the Kaoiki-Hilea area of Hawaii, and the b-values are normal, with b = 0.8 in the top 10 km and somewhat lower values below that depth. We interpret the Kaoiki-Hilea area as relatively stable, normal Hawaiian crust. In contrast, the b-values beneath Kilauea's South Flank are anomalously high (b = 1.3-1.7) at depths between 4 and 8 km, with the highest values near the East Rift zone, but extending 5-8 km away from the rift. Also, the anomalously high b-values vary along strike, parallel to the rift zone. The highest b-values are observed near Hiiaka and Pauahi craters at the bend in the rift, the next highest are near Makaopuhi and also near Puu Kaliu. The mildest anomalies occur adjacent to the central section of the rift. The locations of the three major and two minor b-value anomalies correspond to places where shallow magma reservoirs have been proposed based on analyses of seismicity, geodetic data and differentiated lava chemistry. The existence of the magma reservoirs is also supported by magnetic anomalies, which may be areas of dike concentration, and self-potential anomalies, which are areas of thermal upwelling above a hot source. The simplest explanation of these anomalously high b-values is that they are due to the presence of active magma bodies beneath the East Rift zone at depths down to 8 km. In other volcanoes, anomalously high b-values correlate with volumes adjacent to active magma chambers. This supports a model of a magma body beneath the East Rift zone, which may widen and thin along strike, and which may reach 8 km depth and extend from Kilauea's summit to a distance of at least 40 km down rift. The anomalously high b-values at the center of the South Flank, several kilometers away from the rift, may be explained by unusually high pore pressure throughout the South Flank, or by anomalously strong heterogeneity due to extensive cracking, or by both phenomena. The major b-value anomalies are located SSE of their parent reservoirs, in the direction of motion of the flank, suggesting that magma reservoirs leave an imprint in the mobile flank. We hypothesize that the extensive cracking may have been acquired when the anomalous parts of the South Flank, now several kilometers distant from the rift zone, were generated at the rift zone near persistent reservoirs. Since their generation, these volumes may have moved seaward, away from the rift, but earthquakes occurring in them still use the preexisting complex crack distribution. Along the decollement plane at 10 km depth, the b-values are exceptionally low (b = 0.5), suggesting faulting in a more homogeneous medium. ?? 2001 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗
USGS1999

Hydrated salt minerals on Europa's surface from the Galileo near-infrared mapping spectrometer (NIMS) investigation

We reported evidence of heavily hydrated salt minerals present over large areas of Europa's surface from analysis of reflectance spectra returned by the Galileo mission near infrared mapping spectrometer (NIMS) [ McCord et al. , 1997a, b, 1998a, b]. Here we elaborate on this earlier evidence, present spatial distributions of these minerals, examine alternate water-ice interpretations, expand on our hydrated-salts interpretation, consider salt mineral stability on Europa, and discuss the implications. Extensive well-defined areas on Europa show distinct, asymmetric water-related absorption bands in the 1 to 2.5-μm region. Radiative transfer modeling of water ice involving different particle sizes and layers at Europa temperatures does not reproduce the distinctive Europa water bands. However, ice near its melting temperature, such as in terrestrial environments, does have some characteristics of the Europa spectrum. Alternatively, some classes of heavily hydrated minerals do exhibit such water bands. Among plausible materials, heavily hydrated salt minerals, such as magnesium and sodium sulfates, sodium carbonate and their mixtures, are preferred. All Europa spectral features are present in some salt minerals and a very good match to the Europa spectrum can be achieved by mixing several salt spectra. However, no single or mix of salt mineral spectra from the limited library available has so far been found to perfectly match the Europa spectrum in every detail. The material is concentrated at the lineaments and in chaotic terrain, which are technically disrupted areas on the trailing side. Since the spectrum of the material on Europa is nearly the same everywhere so-far studied, the salt or salt-mixture composition may be nearly uniform. This suggests similar sources and processes over at least a near-hemispheric scale. This would suggest that an extensive subsurface ocean containing dissolved salts is the source, and several possible mechanisms for deposit emplacement are considered. The hydrogen bonds associated with hydration of these salts are similar or greater in strength and energy to those in pure water ice. Thus, once on the surface, the salt minerals should be as stable to disruption as water ice at the Europa temperatures, and mechanisms are suggested to enhance the stability of both materials. Spectra obtained of MgSO 4 ·6H 2 O at 77 K show only small differences from room temperature spectra. The main difference is the appearance of the individual absorptions composing the broad, composite water features and associated with the several different H 2 O sites in the salt hydrate molecule. This suggests that the Europa absorption bands are also composites. Thus higher spectral resolution may reveal these diagnostic features in Europa's spectrum. The specific salts present and their relative abundances would be indicators of the chemistry and conditions of an ocean environment, and areas of fresh, heavy concentration of these minerals should make ideal lander mission sampling sites.

Journal of Geophysical Research E: Planets↗
USGS2008

Trace fossils and sedimentary facies from a Late Cambrian‐Early Ordovician tide‐dominated shelf (Santa Rosita Formation, northwest Argentina): Implications for ichnofacies models of shallow marine successions

The Santa Rosita Formation is one the most widely distributed lower Paleozoic units of northwest Argentina. At the Quebrada del Salto Alto section, east of Purmamarca, Jujuy Province, it is represented by four sedimentary facies: thick‐bedded planar cross‐stratified quartzose sandstones (A), thin‐bedded planar cross‐stratified quartzose sandstones and mudstones (B), wave‐rippled sandstones and bioturbated mudstones (C), and black and greenish gray shales (D). Paleocurrent data, sandstone architecture, and sedimentary structures from faciès A and B indicate bipolar/bimodal paleoflows, suggesting the action of tidal currents. The succession is interpreted as that of a tide‐dominated shelf, with only secondary influence of wave processes. Trace fossils are restricted to facies B and C. The Cruziana ichnocoenosis is preserved on the soles of thin‐bedded planar cross‐stratified quartzose sandstones (faciès B). This ichnocoenosis consists of Conostichus isp., Cruziana omanica, C. semiplicata, C. cf. tortworthi, Cruziana isp. Helminthopsis abeli, Monomorphichnus bilinearis, M. multilineatus, Palaeophycus tubularis, Rusophycus carbonarias, R. latus , and R. isp. The occurrence of Cruziana semiplicata , C. omanica, C. cf. tortworthi , and Rusophycus latus supports a Late Cambrian‐Tremadoc age. Slabbing of Cruziana shows complex interactions between biologic and sedimentologic processes, and suggests a predominance of exhumed traces, washed out and recast by tractive sand deposition. Sandstone soles are densely packed with biogenic structures and exhibit distinctive clusters of Rusophycus isp. that most likely represent trilobite nesting burrows. The Cruziana ichnocoenosis records the resident fauna of a protected, lower intertidal to subtidal interbar setting. The Skolithos ichnocoenosis is represented by high to low density vertical burrows of Skolithos linearis , which extend downwards to the quartzose sandstone soles of faciès B and cross the Cruziana ichnocoenosis. The Skolithos ichnocoenosis represents colonization by suspension‐feeding organisms following a major change in environmental conditions, related to the migration of lower intertidal to subtidal sandwaves. The Planolites ichnocoenosis consists exclusively of Planolites montanus within mudstones overlying wave‐rippled sandstones (facies C). The Planolites ichnocoenosis records opportunistic colonization by inf aunal deposit feeders that mined the organic‐rich fine‐grained sediment during the waning phase of storms that scoured organic detritus from the sea bottom. The section records, from base to top, a Cruziana‐Skolithos ichnofacies zone, a Skolithos ichnofacies zone and an unbioturbated zone typified by the thick‐bedded cross‐stratified quartzose sandstone (fades A). This trend reflects progressively higher energy conditions linked to the establishment of a large sand wave complex. The presence of a mixed Cruziana‐Skolithos ichnofacies in the lower interval reflects changes in substrate and energy levels, rather than water depth. Accordingly, contrasting ichnocoenoses from interbars (Cruziana) and sandwaves (Skolithos) must be considered an example of ichnofacies controlled by local parameters instead of general bathymétrie trends. Conversely, the vertical replacement of the Cruziana ichnofacies by the Skolithos ichnofacies towards the middle interval of the section reflects the environmental changes associated with the transition between the intertidal and subtidal zones. As overall tidal energy increases from supratidal to subtidal settings, the Skolithos ichnofacies tends to occur seaward of the Cruziana ichnofacies in tide‐dominated shallow marine environments. Therefore, onshore‐offshore ichnofacies replacement in tide‐dominated shallow seas is opposite to that in wave‐dominated marine settings.

Ichnos: An International Journal for Plant and Ani↗
USGS1990

Holocene paleoclimatic evidence and sedimentation rates from a core in southwestern Lake Michigan

Preliminary results of a multidisciplinary study of cores in southwestern Lake Michigan suggest that the materials in these cores can be interpreted in terms of both isostatically and climatically induced changes in lake level. Ostracodes and mollusks are well preserved in the Holocene sediments, and they provide paleolimnologic and paleoclimatic data, as well as biogenic carbonate for stable-isotope studies and radiocarbon dating. Pollen and diatom preservation in the cores is poor, which prevents comparison with regional vegetation records. New accelerator-mass spectrometer 14C ages, from both carbon and carbonate fractions, provide basin-wide correlations and appear to resolve the longstanding problem of anomalously old ages that result from detrital organic matter in Great Lakes sediments. Several cores contain a distinct unconformity associated with the abrupt fall in lake level that occurred about 10.3 ka when the isostatically depressed North Bay outlet was uncovered by the retreating Laurentide Ice Sheet. Below the unconformity, ostracode assemblages imply deep, cold water with very low total dissolved solids (TDS), and bivalves have ?? 18O (PDB) values as light as - 10 per mil. Samples from just above the unconformity contain littoral to sublittoral ostracode species that imply warmer, higher-TDS (though still dilute) water than that inferred below the unconformity. Above this zone, another interval with ?? 18O values more negative than - 10 occurs. The isotopic data suggest that two influxes of cold, isotopically light meltwater from Laurentide ice entered the lake, one shortly before 10.3 ka and the other about 9 ka. These influxes were separated by a period during which the lake was warmer, shallower, but still very low in dissolved solids. One or both of the meltwater influxes may be related to discharge from Lake Agassiz into the Great Lakes. Sedimentation rates appear to have been constant from about 10 ka to 5 ka. Bivalve shells formed between about 8 and 5 ka have ?? 18O values that range from-2.3 to-3.3 per mil and appear to decrease toward the end of the interval. The ostracode assemblages and the stable isotopes suggest changes that are climatically controlled, including fluctuating water levels and increasing dissolved solids, although the water remained relatively dilute (TDS < 300 mg/l). A dramatic decrease in sedimentation rates occurred at about 5 ka, about the time of the peak of the Nippissing high lake stage. This decrease in sedimentation rate may be associated with a large increase in effective wave base as the lake approached its present size and fetch. A dramatic reduction in ostracode and mollusk abundances during the late Holocene is probably due to this decrease in sedimentation rates, which would result in increased carbonate dissolution. Ostracode productivity may also have declined due to a reduction in bottom-water oxygen caused by increased epilimnion algal productivity.

Lake Michigan↗
USGS2012

Basin and range crustal and upper mantle structure, northwest to central Nevada

We present an interpretation of the crustal and uppermost mantle structure of the Basin and Range of northwestern Nevada based on seismic refraction/wide‐angle reflection, near‐vertical reflection, and gravity data. In comparison to most previous estimates, we find that the crust is somewhat thicker (32–36 km versus 22–30 km), and the uppermost mantle velocity is somewhat higher (8.0 km/s versus 7.3–7.9 km/s). Along our transects, the crust is thinnest (32 km) in the Carson Sink‐Buena Vista Valley region and increases by 2–4 km to the west and east, respectively. There is considerable complexity throughout the crust where velocities range from of 2.5 km/s at the surface to 7.4 km/s in the lowermost crust. Variations in velocity and structure of the upper crustal layers reveal apparent basement velocity depressions (areas of lower velocities extending up to 10 km in depth) that underlie some surface ranges as well as the basins. The middle crust rises from about 20 km beneath central Nevada to within 12 km of the surface beneath the area of thinnest crust and is characterized by a modest (∼0.1 km/s) change in velocity and low‐velocity gradients. These midcrustal layers mark the onset of high crustal reflectivity and the apparent limiting depth to which Basin and Range faults can be traced in near‐vertical reflection profiles, suggesting that these midcrustal layers represent the transition between the brittle and ductile zones of the crust. The lower crust is more structurally complex, with layers thickening and thinning in a systematic manner with the upper crustal layers; generally, where there are velocity depressions in the upper crust, the lower crust is thickest and shallowest. The geometry of these lower crustal layers (derived from refraction modeling) coincides with changes in the crustal reflectivity, determined from the Consortium of Continental Reflection Profiling reflection data. The lower crustal layer is unusually high in velocity (7.4 km/s) and is likely the layer identified as mantle in some previous studies. We do not identify the 7.4 km/s layer as mantle because (1) there is an underlying layer with a velocity (8.0 km/s) that is more consistent with the worldwide average velocity for the upper mantle, and (2) the 7.4 km/s layer does not correspond to the “reflection” Moho. Gravity modeling and comparison to existing seismic models show a general consensus in many aspects with respect to crustal structure. This new model forms the basis for speculation on some of the processes associated with rifting of the Basin and Range Province. One such process, lithospheric magmatism, is inferred from the strong attenuation of transmitted seismic waves, which occurs at the same interface at which high‐amplitude, bright spot reflections originate. Unlike previous models, the overall structure and velocity of the crust and uppermost mantle of our new model are similar to other regions worldwide which have undergone high degrees of extension.

Nevada↗
USGS2012

Rupture process of a multiple main shock sequence: analysis of teleseismic, local and field observations of the Tennant Creek, Australia, earthquakes of January 22, 1988

On January 22, 1988, three large intraplate earthquakes (with M S 6.3, 6.4, and 6.7) occurred within a 12-hour period near Tennant Creek, Australia. These earthquakes, which occurred over a small interval of time and within a small volume of space, present a unique opportunity to study the rupture process of the class of intraplate earthquakes that occur as multiple main shocks. Broadband displacement and velocity records of body waves from teleseismically recorded data are analyzed to determine source mechanisms, depths, and complexity of rupture of each of the three main shocks. Hypocenters of an additional 150 foreshocks and aftershocks constrained by local arrival time data and field observations of surface rupture are used to complement the source characteristics of the main shocks in order to derive as complete a description of the rupture process as possible. The interpretation of the combined data sets suggests that the overall rupture process involved unusually complicated stress release. As locations of the main shock hypocenters progressively moved from west to east, we infer that the first and third main shocks, denoted as MS 1 and MS 3 , produced the southeast-northwest trending scarps observed at the western end (the Kunayungku fault) and at the eastern end (the east end of the Lake Surprise fault), respectively, of the rupture zone. The epicenter of the only immediate foreshock was located in the gap between these two fault scarps. MS 1 nucleated near this epicenter and ruptured upward and to the northwest from a depth of 6.5 km. MS 3 ruptured predominantly to the SE at a depth of 4.5 km. The second main shock, MS 2 , is inferred to have produced the deformation of the southwest trending central scarp segment (the western end of the Lake Surprise fault). From the sense of thrusting seen at the surface and from the distribution of aftershock hypocenters, the south dipping nodal planes derived from waveform modeling are identified as the fault planes for earthquakes MS1 and MS 3 . In contrast, the dip of the central fault scarp is reversed relative to the dips of the western and eastern fault scarps. The rupture process Of MS 2 turns out to be commensurately complex and sufficiently explains the geological complexity. MS 2 consisted of three subevents. The southeast dipping nodal plane of the first two subevents is coplanar with a southeast dipping plane implied by locations of aftershocks which did not break the surface. Choice of the north dipping plane as the rupture plane of the third subevent, consistent with the surface deformation and coplanar with a second plane delineated by aftershocks, would imply conjugate faulting. The majority of the aftershocks are concentrated near the edges of the fault planes, and there is an absence of activity in the center of the planes. The areas of absent activity may represent the failed asperities of the main shocks in which substantial stress relief occurred. The rupture process of each main shock is characterized by the rapid release of energy followed by a much slower release of moment and by aftershock zones whose dimensions exceed the inferred dimensions of the rupture. These characteristics suggest that substantial slow slip occurred on each of the three fault interfaces that was not accompanied by major energy release. The first main shock nucleated at the deepest part of its aftershock zone and ruptured upward. In contrast, MS 2 and MS 3 nucleated at depths that were one third to one half of the maximum depth of their aftershock zones. This variation of focal depth and the strong increase of moment and radiated energy with each main shock imply that lateral variations of strength were more important than vertical gradients of shear stress in controlling the progression of rupture.

Journal of Geophysical Research Solid Earth↗
USGS2016

Erosional and depositional history of the Atlantic passive margin as recorded in detrital zircon fission-track ages and lithic detritus in Atlantic Coastal plain sediments

Comparison of fission-track (FT) ages of detrital zircons recovered from Atlantic Coastal Plain sediments to FT ages of zircons from bedrock in source terranes in the Appalachians provides a key to understanding the provenance of the sediments and, in turn, the erosional and depositional history of the Atlantic passive margin. In Appalachian source terranes, the oldest zircon fission-track (ZFT) ages from bedrock in the western Appalachians (defined for this paper as the Appalachian Plateau, Valley and Ridge, and far western Blue Ridge) are notably older than the oldest ages from bedrock in the eastern Appalachians (Piedmont and main part of the Blue Ridge). The age difference is seen both in ZFT sample ages and in individual zircon grain ages and reflects differences in the thermotectonic history of the rocks. In the east, ZFT data indicate that the rocks cooled from temperatures high enough to partially or totally reset ZFT ages during the Paleozoic and (or) Mesozoic. The majority of the rocks are interpreted to have cooled through the ZFT closure temperature (∼235 °C) at various times during the late Paleozoic Alleghanian orogeny. In contrast, most of the rocks sampled in the western Appalachians have never been heated to temperatures high enough to totally reset their ZFT ages. Reflecting their contrasting thermotectonic histories, nearly 80 percent of the sampled western rocks yield one or more zircon grains with very old FT ages, in excess of 800 Ma; zircon grains yielding FT ages this old have not been found in rocks in the Piedmont and main part of the Blue Ridge. The ZFT data suggest that the asymmetry of zircon ages of exposed bedrock in the eastern and western Appalachians was in evidence by no later than the Early Cretaceous and probably by the Late Triassic. Detrital zircon suites from sands collected in the Atlantic Coastal Plain provide a record of detritus eroded from source terranes in the Appalachians during the Mesozoic and Cenozoic. In Virginia and Maryland, sands of Early Cretaceous through late early Oligocene age do not yield any old zircons comparable in age to the old zircons found in bedrock in the western Appalachians. Very old zircons yielding FT ages >800 Ma are only encountered in Coastal Plain sands of middle early Miocene and younger age. Miocene and younger fluvial-deltaic deposits associated with the major mid-Atlantic Coastal Plain rivers that now head in the western Appalachians (the Hudson, Delaware, Susquehanna, Potomac, James, and Roanoke) contain abundant clasts of fossiliferous chert and quartzite and other distinctive rock types derived from Paleozoic rocks of the western Appalachians. These distinctive clasts have not been reported in older Coastal Plain sediments. The ZFT and lithic detritus data indicate that the drainage divide for one or more east-flowing mid-Atlantic rivers migrated west into the western Appalachians, and the river(s) began transporting western Appalachian detritus to the Atlantic Coastal Plain, sometime between the late early Oligocene and middle early Miocene. By no later than late middle Miocene most if not all of the major rivers that now head west of the Blue Ridge were transporting western Appalachian detritus to the Coastal Plain. Prior to the drainage divide migrating into the western Appalachians, the ZFT data are consistent with the dominant source of Atlantic Coastal Plain sediments being detritus from the Piedmont and main part of the Blue Ridge, with possible input from distant volcanic sources. The ZFT data suggest that the rapid increase in the rate of siliciclastic sediment accumulation in middle Atlantic margin offshore basins that peaked in the middle Miocene and produced almost 30 percent of the total volume of post-rift siliciclastic sediments in the offshore basins began in the early Miocene when Atlantic river(s) gained access to the relatively easily eroded Paleozoic sedimentary rocks of the western Appalachians.

American Journal of Science↗
USGS1982

Application of automated image analysis to coal petrography

The coal petrologist seeks to determine the petrographic characteristics of organic and inorganic coal constituents and their lateral and vertical variations within a single coal bed or different coal beds of a particular coal field. Definitive descriptions of coal characteristics and coal facies provide the basis for interpretation of depositional environments, diagenetic changes, and burial history and determination of the degree of coalification or metamorphism. Numerous coal core or columnar samples must be studied in detail in order to adequately describe and define coal microlithotypes, lithotypes, and lithologic facies and their variations. The large amount of petrographic information required can be obtained rapidly and quantitatively by use of an automated image-analysis system (AIAS). An AIAS can be used to generate quantitative megascopic and microscopic modal analyses for the lithologic units of an entire columnar section of a coal bed. In our scheme for megascopic analysis, distinctive bands 2 mm or more thick are first demarcated by visual inspection. These bands consist of either nearly pure microlithotypes or lithotypes such as vitrite/vitrain or fusite/fusain, or assemblages of microlithotypes. Megascopic analysis with the aid of the AIAS is next performed to determine volume percentages of vitrite, inertite, minerals, and microlithotype mixtures in bands 0.5 to 2 mm thick. The microlithotype mixtures are analyzed microscopically by use of the AIAS to determine their modal composition in terms of maceral and optically observable mineral components. Megascopic and microscopic data are combined to describe the coal unit quantitatively in terms of (V) for vitrite, (E) for liptite, (I) for inertite or fusite, (M) for mineral components other than iron sulfide, (S) for iron sulfide, and (VEIM) for the composition of the mixed phases (X i ) i = 1,2, etc. in terms of the maceral groups vitrinite V, exinite E, inertinite I, and optically observable mineral content M. The volume percentage of each component present is indicated by a subscript. For example, a lithologic unit was determined megascopically to have the composition (V) 13 (I) 1 (S) 1 (X 1 ) 83 (X 2 ) 2 . After microscopic analysis of the mixed phases, this composition was expressed as (V) 13 (I) 1 (S) 1 (V 63 E 19 I 14 M 4 ) 83 (V 67 E 11 I 13 M 9 ) 2 . Finally, these data were combined in a description of the bulk composition as V 67 E 16 I 13 M 3 S 1 . An AIAS can also analyze textural characteristics and can be used for quick and reliable determination of rank (reflectance). Our AIAS is completely software based and incorporates a television (TV) camera that has optimum response characteristics in the range of reflectance less than 5%, making it particularly suitable for coal studies. Analysis of the digitized signal from the TV camera is controlled by a microprocessor having a resolution of 64 gray levels between full illumination and dark current. The processed image is reconverted for display on a TV monitor screen, on which selection of phases or features to be analyzed is readily controlled and edited by the operator through use of a lightpen. We expect that automated image analysis, because it can rapidly provide a large amount of pertinent information, will play a major role in the advancement of coal petrography.

International Journal of Coal Geology↗
USGS2008

The evolution of volcano-hosted geothermal systems based on deep wells from Karaha-Telaga Bodas, Indonesia

Temperature and pressure surveys, fluid samples, and petrologic analyses of rock samples from deep drill holes at the Karaha - Telaga Bodas geothermal field on the volcanic ridge extending northward from Galunggung Volcano, West Java, have provided a unique opportunity to characterize the evolution of an active volcano-hosted geothermal system. Wells up to 3 km in depth have encountered temperatures as high as 353??C and a weakly altered granodiorite that intruded to within 2 to 3 km of the surface. The intrusion is shallowest beneath the southern end of the field where an acid lake overlies a nearly vertical low resistivity structure (<10 ohm-m) defined by magnetotelluric measurements. This structure is interpreted to represent a vapor-dominated chimney that provides a pathway to the surface for magmatic gases. Four distinct hydrothermal mineral assemblages document the evolution of the geothermal system and the transition from liquid- to vapor-dominated conditions. The earliest assemblage represents the initial liquid-dominated system generated during emplacement of the granodiorite between 5910 ?? 76 and 4200 ?? 150 y BP. Tourmaline, biotite, actinolite, epidote and clay minerals were deposited contemporaneously at progressively greater distances from the intrusive contact (assemblage 1). At 4200 ?? 150 y BP, flank collapse and the formation of the volcano's crater, Kawah Galunggung, resulted in catastrophic decompression and boiling of the hydrothermal fluids. This event initiated development of the modern vapor-dominated regime. Chalcedony and then quartz were deposited as the early low salinity liquids boiled (assemblage 2). Both vapor- and liquid-rich fluid inclusions were trapped in the quartz crystals. Liquid-rich fluid inclusions from the southern part of the field record salinities ranging from 0 to 26 weight percent NaCl- CaCl2 equivalent and locally contain fluorite daughter crystals. We suggest, based on temperature-salinity relationships and evidence of boiling, that these fluids were progressively concentrated as steam was lost from the system. However, mixing with fluids derived from the underlying intrusion or generated during the formation of acid SO4 water on the vapor-dominated chimney margins could have contributed to the observed salinities. As pressures declined, CO2- and SO4-rich steam-heated water drained downward, depositing anhydrite and calcite (assemblage 3) in the fractures, limiting further recharge. Fluid inclusions with salinities up to 31 weight percent NaCl equivalent were trapped in these minerals as the descending water vaporized. The final assemblage is represented by precipitates of NaCl, KCl and FeClx deposited on rock surfaces in portions of the vapor-dominated zone that boiled dry. Vapor-dominated conditions extend over a distance of at least 10 km and to depths below sea level. Deep wells drilled into the underlying liquid-dominated reservoir in the northern and central part of the volcanic ridge produce low salinity fluids representing recent recharge of meteoric and steam-heated water. The evolution of volcanic-hosted vapor-dominated geothermal systems can be described by a five stage model. Stage 1 involves the formation of an over-pressured liquid-dominated geothermal system soon after magmatic intrusion. In Stages 2 and 3, pressures progressively decrease, and a curtain of steam-heated water surrounding a magmatic vapor-dominated chimney at 350??C and 14 ?? 2 MPa develops. The relatively low pressure near the base of the chimney causes liquid inflow adjacent to the intrusion and the development of a secondary marginal vapor-dominated zone. In Stage 4, the magmatic vapor discharge from the intrusion becomes small, vapor pressure declines, and the secondary vapor-dominated zone expands above the intrusion. In Stage 5, the vapor-dominated zone floods because heat from the intrusion is insufficient to boil all liquid inflow. A more common, liquid-dominated volcanic-hosted system the

American Journal of Science↗