skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: The dynamic floor of Yellowstone Lake, Wyoming, USA: The last 14 k.y. of hydrothermal explosions, venting, doming, and faulting
Abstract Hydrothermal explosions are significant potential hazards in Yellowstone National Park, Wyoming, USA. The northern Yellowstone Lake area hosts the three largest hydrothermal explosion craters known on Earth empowered by the highest heat flow values in Yellowstone and active seismicity and deformation. Geological and geochemical studies of eighteen sublacustrine cores provide the first detailed synthesis of the age, sedimentary facies, and origin of multiple hydrothermal explosion deposits. New tephrochronology and radiocarbon results provide a four-dimensional view of recent geologic activity since recession at ca. 15–14.5 ka of the >1-km-thick Pinedale ice sheet. The sedimentary record in Yellowstone Lake contains multiple hydrothermal explosion deposits ranging in age from ca. 13 ka to ~1860 CE. Hydrothermal explosions require a sudden drop in pressure resulting in rapid expansion of high-temperature fluids causing fragmentation, ejection, and crater formation; explosions may be initiated by seismicity, faulting, deformation, or rapid lake-level changes. Fallout and transport of ejecta produces distinct facies of subaqueous hydrothermal explosion deposits. Yellowstone hydrothermal systems are characterized by alkaline-Cl and/or vapor-dominated fluids that, respectively, produce alteration dominated by silica-smectite-chlorite or by kaolinite. Alkaline-Cl liquids flash to steam during hydrothermal explosions, producing much more energetic events than simple vapor expansion in vapor-dominated systems. Two enormous explosion events in Yellowstone Lake were triggered quite differently: Elliott’s Crater explosion resulted from a major seismic event (8 ka) that ruptured an impervious hydrothermal dome, whereas the Mary Bay explosion (13 ka) was triggered by a sudden drop in lake level stimulated by a seismic event, tsunami, and outlet channel erosion.  more » « less
Award ID(s):
2149482
PAR ID:
10406060
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
GSA Bulletin
Volume:
135
Issue:
3-4
ISSN:
0016-7606
Page Range / eLocation ID:
547 to 574
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Chemical changes in hot springs, as recorded by thermal waters and their deposits, provide a window into the evolution of the postglacial hydrothermal system of the Yellowstone Plateau Volcanic Field. Today, most hydrothermal travertine forms to the north and south of the ca. 631 ka Yellowstone caldera where groundwater flow through subsurface sedimentary rocks leads to calcite saturation at hot springs. In contrast, low-Ca rhyolites dominate the subsurface within the Yellowstone caldera, resulting in thermal waters that rarely deposit travertine. We investigated the timing and origin of five small travertine deposits in the Upper and Lower Geyser Basins to understand the conditions that allowed for travertine deposition. New 230Th-U dating, oxygen (δ18O), carbon (δ13C), and strontium (87Sr/86Sr) isotopic ratios, and elemental concentrations indicate that travertine deposits within the Yellowstone caldera formed during three main episodes that correspond broadly with known periods of wet climate: 13.9−13.6 ka, 12.2−9.5 ka, and 5.2−2.9 ka. Travertine deposition occurred in response to the influx of large volumes of cold meteoric water, which increased the rate of chemical weathering of surficial sediments and recharge into the hydrothermal system. The small volume of intracaldera travertine does not support a massive postglacial surge of CO2 within the Yellowstone caldera, nor was magmatic CO2 the catalyst for postglacial travertine deposition. 
    more » « less
  2. Glacial and periglacial sediments and landforms record the chronology of glaciation and amount of Pleistocene erosion during colder periods that added substantially to global sediment budgets and contributed to the global CO2 cycle. The now-drained glacial Lake Devlin, dammed in a Front Range tributary valley by a glacier in the North Branch of Boulder Creek (Colorado, USA) preserves an important sedimentary archive of the ca. 32−14 ka Pinedale glaciation, recording both paleoclimate information and an integrated measure of glacial and periglacial erosion rates over a full glacial cycle. Despite rapid erosion of fine-grained deposits after the lake drained, most sediment generated during Pinedale time remains as legacy deposits in the catchment. Geomorphic evidence and dating of glaciolacustrine sediment from surface exposures demonstrate that the ca. 30 ka Pinedale glacial advance was nearly as extensive as the local Late Glacial Maximum at ca. 20 ka. Sedimentary archives dated by 14C, optically stimulated luminescence, and cosmogenic nuclides extend earlier studies (Madole et al., 1973) of pollen and magnetic susceptibility (MS) in cores from the glaciolacustrine deposits of Lake Devlin and of Pinedale climate. Records suggest short-term warming and biotic change at ca. 15 ka after ∼14 kyr of cold, dry conditions punctuated by MS peaks at ca. 26.5 ka, 20 ka, and 16.5 ka. Lake Devlin drained catastrophically after ca. 14 ka, millennia after ice had retreated upvalley from the lateral moraine that dammed the lake. Sediment production during the Pinedale was equivalent to a periglacial and glacial erosion rate of ∼70 mm kyr−1, several times higher than long-term rates in the adjacent Front Range, but much lower than rates measured where modern glaciers are eroding weak bedrock in zones of rapid rock uplift, such as SSE Alaska, USA. Data from the Lake Devlin basin contribute to contemporary discussions of how glacial erosion influences the global CO2 cycle. 
    more » « less
  3. Abstract We describe and interpret deposits associated with the final Ubehebe Crater-forming, phreatomagmatic explosive phase of the multivent, monogenetic Ubehebe volcanic center. Ubehebe volcano is located in Death Valley, California, USA. Pyroclastic deposits occur in four main facies: (1) lapilli- and blockdominated beds, (2) thinly bedded lapilli tuff, (3) laminated and cross-laminated ash, and (4) massive lapilli ash/tuff. Lapilli- and block-dominated beds are found mostly within several hundred meters of the crater and transition outward into discontinuous lenses of lapilli and blocks; they are interpreted to have been deposited by ballistic processes associated with crater-forming explosions. Thinly bedded lapilli tuff is found mainly within several hundred meters, and laminated and cross-laminated ash extends at least 9 km from the crater center. Dune forms are common within ~2 km of the crater center, while finer-grained, distal deposits tend to exhibit planar lamination. These two facies (thinly bedded lapilli tuff and laminated and cross-laminated ash) are interpreted to record multiple pyroclastic surges (dilute pyroclastic currents). Repeated couplets of coarse layers overlain by finer-grained, laminated horizons suggest that many or most of the surges were transient, likely recording individual explosions, and they traveled over complex topography in some areas. These two factors complicate the application of classical sediment-transport theory to quantify surge properties. However, dune-form data provide possible constraints on the relationships between suspended load sedimentation and bed-load transport that are consistent using two independent approaches. Massive lapilli ash/tuff beds occur in drainages below steep slopes and can extend up to ~1 km onto adjacent valley floors beneath large catchments. Although they are massive in texture, their grain-size characteristics are shared with laminated and cross-laminated ash facies, with which they are locally interbedded. These are interpreted to record concentrated granular flows sourced by remobilized pyroclastic surge deposits, either during surge transport or shortly after, while the surge deposits retained their elevated initial pore-gas pressures. Although similar surge-derived concentrated flows have been described elsewhere (e.g., Mount St. Helens, Washington, USA, and Soufriére Hills, Montserrat, West Indies), to our knowledge Ubehebe is the first case where such processes have been identified at a maar volcano. These concentrated flows followed paths that were independent of the pyroclastic surges and represent a potential hazard at similar maar volcanoes in areas with complex terrain. 
    more » « less
  4. Rare earth element (REE) deposits are commonly associated with carbonatites and (per)alkaline rocks where hydrothermal magmatic fluids can play a significant role in REE mobilization and deposition [1]. Thermodynamic modeling permits predicting the evolution of ore-forming fluids and can be used to test different controls on hydrothermal REE mobility including temperature, pressure, the solubility of REE minerals, aqueous REE speciation and pH evolution associated with fluid-rock interaction. Previous modeling studies either focused on REE fluoride/chloride complexation in acidic aqueous fluids [2] or near neutral/alkaline fluids associated with calcite vein formation [3]. Such models were also applied to interpret field observations in REE deposits Bayan Obo in China and Bear Lodge in Wyoming [3,4]. Recent hydrothermal calcite-fluid REE partitioning experiments provide new data to simulate the solubility of REE in calcite, REE carbonates/fluorocarbonates at high temperatures [5, 6]. We studied the competing effects controlling the mobility of REE in hydrothermal fluids between 100 and 400 °C at 500 bar. Speciation calculations were carried out in the Ca-F-CO2-Na-Cl-H2O system using the GEMS code package [7]. The properties of minerals and aqueous species were taken from the MINES thermodynamic database [3,5]. The Gallinas Mountains hydrothermal REE deposit in New Mexico was used as a field analogue to compare our models with the formation of calcite-fluorite veins hosting bastnäsite. Previous fluid inclusion studies hypothesized that the REE were transported as fluoride complexes [8] but more recent modeling studies have shown that fluoride essentially acts as a depositing ligand [2]. Here we show more detailed simulations predicting the stability of fluorite, calcite and REE minerals relevant to ore-forming processes in carbonatites and alkaline systems. [1] Gysi et al. (2016), Econ. Geol. 111, 1241-1276; [2] Migdisov and Williams-Jones (2014), Mineral. Deposita 49, 987-997. [3] Perry and Gysi (2018), Geofluids; [4] Liu et al. (2020), Minerals 10, 495; [5] Perry and Gysi (2020), Geochim. Cosmochim. Acta 286, 177-197; [6] Gysi and Williams-Jones (2015) Chem. Geol. 392, 87-101;[7] Kulik et al. (2013), Computat. Geosci. 17, 1-24; [8] Williams-Jones et al. (2000), Econ. Geol. 95, 327-341 
    more » « less
  5. Abstract Surface deformation and seismicity provide critical information to understand the dynamics of volcanic unrest. During 2006–2007, >80 mm/yr uplift was observed by interferometric synthetic aperture radar (InSAR) at the central Atka volcanic center, Alaska, coinciding with an increasing seismicity rate. On November 25, 2006, a phreatic eruption occurred at the Korovin volcanic vent, 5‐km north of the central Atka, following the drainage of its crater lake a month prior to the eruption. The InSAR data are assimilated into three‐dimensional finite element models using the Ensemble Kalman Filter to investigate: (1) the pressure source creating the surface deformation; (2) the triggering of the volcano‐tectonic (VT) earthquakes in the Atka volcanic center; and (3) the triggering of the phreatic eruption at Korovin. The models show that the pressure source required to create the surface deformation is a NE‐tilted, oblate ellipsoid, which rotated from steep to gentle dipping from June to November 2006 before the eruption. The modeled dilatancy in a pre‐existing weak zone, coinciding with the Amlia‐Amukta fault, driven by the pressure source has a spatial and temporal correlation with the evolution of the VT earthquakes during the unrest. The fault dilatancy may have increased the connected porosity and permeability of the fault zone allowing fluid injection which triggered the observed seismicity. In addition, the dilatated fault may have increased the fluid capacity of the fault zone by ∼105 m3, causing the discharge of the crater lake at Korovin. Consequently, the phreatic eruption of the Korovin volcano may have been triggered. 
    more » « less