Search USGSSearch

SEARCH · Search USGS

Results for “Geochemistry, Geophysics, and Geosystems”

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 37 records · Page 2Linked to original sources

Evolution of fracture permeability of ultramafic rocks undergoing serpentinization at hydrothermal conditions: An experimental study

We performed flow-through laboratory experiments on five cylindrically cored samples of ultramafic rocks, in which we generated a well-mated through-going tensile fracture, to investigate evolution of fracture permeability during serpentinization. The samples were tested in a triaxial loading machine at a confining pressure of 50 MPa, pore pressure of 20 MPa, and temperature of 260°C, simulating a depth of 2 km under hydrostatic conditions. A pore pressure difference of up to 2 MPa was imposed across the ends of the sample. Fracture permeability decreased by 1–2 orders of magnitude during the 200–330 h experiments. Electron microprobe and SEM data indicated the formation of needle-shaped crystals of serpentine composition along the walls of the fracture, and chemical analyses of sampled pore fluids were consistent with dissolution of ferro-magnesian minerals. By comparing the difference between fracture permeability and matrix permeability measured on intact samples of the same rock types, we concluded that the contribution of the low matrix permeability to flow is negligible and essentially all of the flow is focused in the tensile fracture. The experimental results suggest that the fracture network in long-lived hydrothermal circulation systems can be sealed rapidly as a result of mineral precipitation, and generation of new permeability resulting from a combination of tectonic and crystallization-induced stresses is required to maintain fluid circulation.

Geochemistry, Geophysics, Geosystems

Surface slip during large Owens Valley earthquakes

The 1872 Owens Valley earthquake is the third largest known historical earthquake in California. Relatively sparse field data and a complex rupture trace, however, inhibited attempts to fully resolve the slip distribution and reconcile the total moment release. We present a new, comprehensive record of surface slip based on lidar and field investigation, documenting 162 new measurements of laterally and vertically displaced landforms for 1872 and prehistoric Owens Valley earthquakes. Our lidar analysis uses a newly developed analytical tool to measure fault slip based on cross‐correlation of sublinear topographic features and to produce a uniquely shaped probability density function (PDF) for each measurement. Stacking PDFs along strike to form cumulative offset probability distribution plots (COPDs) highlights common values corresponding to single and multiple‐event displacements. Lateral offsets for 1872 vary systematically from ∼1.0 to 6.0 m and average 3.3 ± 1.1 m (2σ). Vertical offsets are predominantly east‐down between ∼0.1 and 2.4 m, with a mean of 0.8 ± 0.5 m. The average lateral‐to‐vertical ratio compiled at specific sites is ∼6:1. Summing displacements across subparallel, overlapping rupture traces implies a maximum of 7–11 m and net average of 4.4 ± 1.5 m, corresponding to a geologic M w ∼7.5 for the 1872 event. We attribute progressively higher‐offset lateral COPD peaks at 7.1 ± 2.0 m, 12.8 ± 1.5 m, and 16.6 ± 1.4 m to three earlier large surface ruptures. Evaluating cumulative displacements in context with previously dated landforms in Owens Valley suggests relatively modest rates of fault slip, averaging between ∼0.6 and 1.6 mm/yr (1σ) over the late Quaternary.

Book

A possible transoceanic tsunami directed toward the U.S. west coast from the Semidi segment, Alaska convergent margin

The Semidi segment of the Alaska convergent margin appears capable of generating a giant tsunami like the one produced along the nearby Unimak segment in 1946. Reprocessed legacy seismic reflection data and a compilation of multibeam bathymetric surveys reveal structures that could generate such a tsunami. A 200 km long ridge or escarpment with crests >1 km high is the surface expression of an active out-of-sequence fault zone, recently referred to as a splay fault. Such faults are potentially tsunamigenic. This type of fault zone separates the relatively rigid rock of the margin framework from the anelastic accreted sediment prism. Seafloor relief of the ridge exceeds that of similar age accretionary prism ridges indicating preferential slip along the splay fault zone. The greater slip may derive from Quaternary subduction of the Patton Murray hot spot ridge that extends 200 km toward the east across the north Pacific. Estimates of tsunami repeat times from paleotsunami studies indicate that the Semidi segment could be near the end of its current inter-seismic cycle. GPS records from Chirikof Island at the shelf edge indicate 90% locking of plate interface faults. An earthquake in the shallow Semidi subduction zone could generate a tsunami that will inundate the US west coast more than the 1946 and 1964 earthquakes because the Semidi continental slope azimuth directs a tsunami southeastward.

Geochemistry, Geophysics, Geosystems

Origin and dynamics of depositionary subduction margins

Here we propose a new framework for forearc evolution that focuses on the potential feedbacks between subduction tectonics, sedimentation, and geomorphology that take place during an extreme event of subduction erosion. These feedbacks can lead to the creation of a “depositionary forearc,” a forearc structure that extends the traditional division of forearcs into accretionary or erosive subduction margins by demonstrating a mode of rapid basin accretion during an erosive event at a subduction margin. A depositionary mode of forearc evolution occurs when terrigenous sediments are deposited directly on the forearc while it is being removed from below by subduction erosion. In the most extreme case, an entire forearc can be removed by a single subduction erosion event followed by depositionary replacement without involving transfer of sediments from the incoming plate. We need to further recognize that subduction forearcs are often shaped by interactions between slow, long-term processes, and sudden extreme events reflecting the sudden influences of large-scale morphological variations in the incoming plate. Both types of processes contribute to the large-scale architecture of the forearc, with extreme events associated with a replacive depositionary mode that rapidly creates sections of a typical forearc margin. The persistent upward diversion of the megathrust is likely to affect its geometry, frictional nature, and hydrogeology. Therefore, the stresses along the fault and individual earthquake rupture characteristics are also expected to be more variable in these erosive systems than in systems with long-lived megathrust surfaces.

Geochemistry, Geophysics, Geosystems

Determining the flux of methane into Hudson Canyon at the edge of methane clathrate hydrate stability

Methane seeps were investigated in Hudson Canyon, the largest shelf-break canyon on the northern US Atlantic Margin. The seeps investigated are located at or updip of the nominal limit of methane clathrate hydrate stability. The acoustic identification of bubble streams was used to guide water column sampling in a 32 km2 region within the canyon's thalweg. By incorporating measurements of dissolved methane concentration with methane oxidation rates and current velocity into a steady-state box model, the total emission of methane to the water column in this region was estimated to be 12 kmol methane per day (range: 6 – 24 kmol methane per day). These analyses suggest this methane is largely retained inside the canyon walls below 300 m water depth, and that it is aerobically oxidized to near completion within the larger extent of Hudson Canyon. Based on estimated methane emissions and measured oxidation rates, the oxidation of this methane to dissolved CO2 is expected to have minimal influences on seawater pH. This article is protected by copyright. All rights reserved.

Geochemistry, Geophysics, Geosystems

3-D high-speed imaging of volcanic bomb trajectory in basaltic explosive eruptions

Imaging, in general, and high speed imaging in particular are important emerging tools for the study of explosive volcanic eruptions. However, traditional 2-D video observations cannot measure volcanic ejecta motion toward and away from the camera, strongly hindering our capability to fully determine crucial hazard-related parameters such as explosion directionality and pyroclasts' absolute velocity. In this paper, we use up to three synchronized high-speed cameras to reconstruct pyroclasts trajectories in three dimensions. Classical stereographic techniques are adapted to overcome the difficult observation conditions of active volcanic vents, including the large number of overlapping pyroclasts which may change shape in flight, variable lighting and clouding conditions, and lack of direct access to the target. In particular, we use a laser rangefinder to measure the geometry of the filming setup and manually track pyroclasts on the videos. This method reduces uncertainties to 10° in azimuth and dip angle of the pyroclasts, and down to 20% in the absolute velocity estimation. We demonstrate the potential of this approach by three examples: the development of an explosion at Stromboli, a bubble burst at Halema'uma'u lava lake, and an in-flight collision between two bombs at Stromboli.

Geochemistry, Geophysics, Geosystems

Subsea ice-bearing permafrost on the U.S. Beaufort Margin: 2. Borehole constraints

Borehole logging data from legacy wells directly constrain the contemporary distribution of subsea permafrost in the sedimentary section at discrete locations on the U.S. Beaufort Margin and complement recent regional analyses of exploration seismic data to delineate the permafrost's offshore extent. Most usable borehole data were acquired on a ∼500 km stretch of the margin and within 30 km of the contemporary coastline from north of Lake Teshekpuk to nearly the U.S.-Canada border. Relying primarily on deep resistivity logs that should be largely unaffected by drilling fluids and hole conditions, the analysis reveals the persistence of several hundred vertical meters of ice-bonded permafrost in nearshore wells near Prudhoe Bay and Foggy Island Bay, with less permafrost detected to the east and west. Permafrost is inferred beneath many barrier islands and in some nearshore and lagoonal (back-barrier) wells. The analysis of borehole logs confirms the offshore pattern of ice-bearing subsea permafrost distribution determined based on regional seismic analyses and reveals that ice content generally diminishes with distance from the coastline. Lacking better well distribution, it is not possible to determine the absolute seaward extent of ice-bearing permafrost, nor the distribution of permafrost beneath the present-day continental shelf at the end of the Pleistocene. However, the recovery of gas hydrate from an outer shelf well (Belcher) and previous delineation of a log signature possibly indicating gas hydrate in an inner shelf well (Hammerhead 2) imply that permafrost may once have extended across much of the shelf offshore Camden Bay.

Geochemistry, Geophysics, Geosystems

Subsea ice-bearing permafrost on the U.S. Beaufort Margin: 1. Minimum seaward extent defined from multichannel seismic reflection data

Subsea ice-bearing permafrost (IBPF) and associated gas hydrate in the Arctic have been subject to a warming climate and saline intrusion since the last transgression at the end of the Pleistocene. The consequent degradation of IBPF is potentially associated with significant degassing of dissociating gas hydrate deposits. Previous studies interpreted the distribution of subsea permafrost on the U.S. Beaufort continental shelf based on geographically sparse data sets and modeling of expected thermal history. The most cited work projects subsea permafrost to the shelf edge (∼100 m isobath). This study uses a compilation of stacking velocity analyses from ∼100,000 line-km of industry-collected multichannel seismic reflection data acquired over 57,000 km 2 of the U.S. Beaufort shelf to delineate continuous subsea IBPF. Gridded average velocities of the uppermost 750 ms two-way travel time range from 1475 to 3110 m s −1 . The monotonic, cross-shore pattern in velocity distribution suggests that the seaward extent of continuous IBPF is within 37 km of the modern shoreline at water depths < 25 m. These interpretations corroborate recent Beaufort seismic refraction studies and provide the best, margin-scale evidence that continuous subsea IBPF does not currently extend to the northern limits of the continental shelf.

Geochemistry, Geophysics, Geosystems

Seismic imaging of the metamorphism of young sediment into new crystalline crust in the actively rifting Imperial Valley, California

Plate-boundary rifting between transform faults is opening the Imperial Valley of southern California and the rift is rapidly filling with sediment from the Colorado River. Three 65–90 km long seismic refraction profiles across and along the valley, acquired as part of the 2011 Salton Seismic Imaging Project, were analyzed to constrain upper crustal structure and the transition from sediment to underlying crystalline rock. Both first arrival travel-time tomography and frequency-domain full-waveform inversion were applied to provide P-wave velocity models down to ∼7 km depth. The valley margins are fault-bounded, beyond which thinner sediment has been deposited on preexisting crystalline rocks. Within the central basin, seismic velocity increases continuously from ∼1.8 km/s sediment at the surface to >6 km/s crystalline rock with no sharp discontinuity. Borehole data show young sediment is progressively metamorphosed into crystalline rock. The seismic velocity gradient with depth decreases approximately at the 4 km/s contour, which coincides with changes in the porosity and density gradient in borehole core samples. This change occurs at ∼3 km depth in most of the valley, but at only ∼1.5 km depth in the Salton Sea geothermal field. We interpret progressive metamorphism caused by high heat flow to be creating new crystalline crust throughout the valley at a rate comparable to the ≥2 km/Myr sedimentation rate. The newly formed crystalline crust extends to at least 7–8 km depth, and it is shallower and faster where heat flow is higher. Most of the active seismicity occurs within this new crust.

California

Fitful and protracted magma assembly leading to a giant eruption, Youngest Toba Tuff, Indonesia

The paroxysmal eruption of the 74 ka Youngest Toba Tuff (YTT) of northern Sumatra produced an extraordinary 2800 km 3 of non-welded to densely welded ignimbrite and co-ignimbrite ash-fall. We report insights into the duration of YTT magma assembly obtained from ion microprobe U-Th and U-Pb dates, including continuous age spectra over >50% of final zircon growth, for pumices and a welded tuff spanning the compositional range of the YTT. A relatively large subpopulation of zircon crystals nucleated before the penultimate caldera-related eruption at 501 ka, but most zircons yielded interior dates 100-300 ka thereafter. Zircon nucleation and growth was likely episodic and from diverse conditions over protracted time intervals of >100 to >500 ka. Final zircon growth is evident as thin rim plateaus that are in Th/U chemical equilibrium with hosts, and that give crystallization ages within tens of ka of eruption. The longevity and chemical characteristics of the YTT zircons, as well as evidence for intermittent zircon isolation and remobilization associated with magma recharge, is especially favored at the cool and wet eutectoid conditions that characterize at least half of the YTT, wherein heat fluxes could dissolve major phases but have only a minor effect on larger zircon crystals. Repeated magma recharge may have contributed to the development of compositional zoning in the YTT but, considered together with limited allanite, quartz, and other mineral dating and geospeedometry, regular perturbations to the magma reservoir over >400 ka did not lead to eruption until 74 ka ago.

Lake Toba

Post-rift magmatic evolution of the eastern North American “passive-aggressive” margin

Understanding the evolution of passive margins requires knowledge of temporal and chemical constraints on magmatism following the transition from supercontinent to rifting, to post-rifting evolution. The Eastern North American Margin (ENAM) is an ideal study location as several magmatic pulses occurred in the 200 My following rifting. In particular, the Virginia-West Virginia region of the ENAM has experienced two postrift magmatic pulses at ∼152 Ma and 47 Ma, and thus provides a unique opportunity to study the long-term magmatic evolution of passive margins. Here we present a comprehensive set of geochemical data that includes new 40 Ar/ 39 Ar ages, major and trace-element compositions, and analysis of radiogenic isotopes to further constrain their magmatic history. The Late Jurassic volcanics are bimodal, from basanites to phonolites, while the Eocene volcanics range from picrobasalt to rhyolite. Modeling suggests that the felsic volcanics from both the Late Jurassic and Eocene events are consistent with fractional crystallization. Sr-Nd-Pb systematics for the Late Jurassic event suggests HIMU and EMII components in the magma source that we interpret as upper mantle components rather than crustal interaction. Lithospheric delamination is the best hypothesis for magmatism in Virginia/West Virginia, due to tectonic instabilities that are remnant from the long-term evolution of this margin, resulting in a “passive-aggressive” margin that records multiple magmatic events long after rifting ended.

Geochemistry, Geophysics, Geosystems

Observed correlation between the depth to base and top of gas hydrate occurrence from review of global drilling data

A global inventory of data from gas hydrate drilling expeditions is used to develop relationships between the base of structure I gas hydrate stability, top of gas hydrate occurrence, sulfate-methane transition depth, pressure (water depth), and geothermal gradients. The motivation of this study is to provide first-order estimates of the top of gas hydrate occurrence and associated thickness of the gas hydrate occurrence zone for climate-change scenarios, global carbon budget analyses, or gas hydrate resource assessments. Results from publicly available drilling campaigns (21 expeditions and 52 drill sites) off Cascadia, Blake Ridge, India, Korea, South China Sea, Japan, Chile, Peru, Costa Rica, Gulf of Mexico, and Borneo reveal a first-order linear relationship between the depth to the top and base of gas hydrate occurrence. The reason for these nearly linear relationships is believed to be the strong pressure and temperature dependence of methane solubility in the absence of large difference in thermal gradients between the various sites assessed. In addition, a statistically robust relationship was defined between the thickness of the gas hydrate occurrence zone and the base of gas hydrate stability (in meters below seafloor). The relationship developed is able to predict the depth of the top of gas hydrate occurrence zone using observed depths of the base of gas hydrate stability within less than 50 m at most locations examined in this study. No clear correlation of the depth to the top and base of gas hydrate occurrences with geothermal gradient and sulfate-methane transition depth was identified.

Geochemistry, Geophysics, Geosystems

Holocene paleointensity of the Island of Hawai'i from glassy volcanics

This study presents new high‐quality paleointensity records and 14 C radiocarbon age determinations from the Island of Hawai `i during the Holocene. Previous studies on Hawai `i use experimental methods and statistical selection criteria that may produce inaccurate geomagnetic field strength estimates. Additional high‐quality paleointensity results can be used to evaluate the existing Hawaiian data set and investigate Holocene geomagnetic field behavior. New paleointensity sites from 22 lava flows were calculated using the IZZI‐Thellier laboratory technique and a strict set of selection criteria. Rapidly cooled, glassy volcanic material was collected for all sites. Isotopic age determinations range from 270 to >10, 000 years before present (nine new 14 C ages are also presented as part of this study). The median intensity for the 22 flows is 47.5 μ T, with a median absolute deviation uncertainty of 5.6 μ T; substantially greater than the present‐day field strength at Hawai `i (~36 μ T). These new results are comparable to previously published data from this location and are consistent with global paleointensity models. There is no evidence of an intensity “spike” at 3,000 years before present, as seen in the Levant and elsewhere. Previously published data vary in intensity by experimental technique relative to data using glassy material and strict selection criteria. Non‐Thellier‐type data are biased low, a result of these techniques estimating intensity from possibly nonsingle domain magnetic carriers. Thellier‐Thellier data are biased high, the reasons for which remain unclear as no cooling rate effect was demonstrated, and we were unable to reproduce the high bias with different selection criteria.

Hawai`i

Reconstructing the evolution of the submarine Monterey Canyon System from Os, Nd, and Pb isotopes in hydrogenetic Fe-Mn crusts

The sources of terrestrial material delivered to the California margin over the past 7 Myr were assessed using 187 Os/ 188 Os, Nd, and Pb isotopes in hydrogenetic ferromanganese crusts from three seamounts along the central and southern California margin. From 6.8 to 4.5 (± 0.5) Ma, all three isotope systems show more radiogenic values at Davidson Seamount, located near the base of the Monterey Canyon System, than in Fe-Mn crusts from the more remote Taney and Hoss seamounts. At the Taney seamounts, approximately 225 km farther offshore from Davidson Seamount, 187 Os/ 188 Os values, but not Pb and Nd isotope ratios, also deviate from the Cenozoic seawater curve towards more radiogenic values from 6.8 to 4.5 (± 0.5) Ma. However, none of the isotope systems in Fe-Mn crusts deviate from seawater at Hoss Seamount located approximately 450 km to the south. The regional gradients in isotope ratios indicate that substantial input of dissolved and particulate terrestrial material into the Monterey Canyon System is responsible for the local deviations in the seawater Nd, Pb, and Os isotope compositions from 6.8 to 4.5 (± 0.5) Ma. The isotope ratios recorded in Fe-Mn crusts are consistent with a southern Sierra Nevada or western Basin and Range provenance of the terrestrial material which was delivered by rivers to the canyon. The exhumation of the modern Monterey Canyon must have begun between 10 and 6.8 ± 0.5 Ma, as indicated by our data, the age of incised strata, and paleo-location of the Monterey Canyon relative to the paleo-coastline.

California

Significance of northeast-trending features in Canada Basin, Arctic Ocean

Synthesis of seismic velocity, potential field, and geological data from Canada Basin and its surrounding continental margins suggests that a northeast-trending structural fabric has influenced the origin, evolution, and current tectonics of the basin. This structural fabric has a crustal origin, based on the persistence of these trends in upward continuation of total magnetic intensity data and vertical derivative analysis of free-air gravity data. Three subparallel northeast-trending features are described. Northwind Escarpment, bounding the east side of the Chukchi Borderland, extends ∼600 km and separates continental crust of Northwind Ridge from high-velocity transitional crust in Canada Basin. A second, shorter northeast-trending zone extends ∼300 km in northern Canada Basin and separates inferred continental crust of Sever Spur from magmatically intruded crust of the High Arctic Large Igneous Province. A third northeast-trending feature, here called the Alaska-Prince Patrick magnetic lineament (APPL) is inferred from magnetic data and its larger regional geologic setting. Analysis of these three features suggests strike slip or transtensional deformation played a role in the opening of Canada Basin. These features can be explained by initial Jurassic-Early Cretaceous strike slip deformation (phase 1) followed in the Early Cretaceous (∼134 to ∼124 Ma) by rotation of Arctic Alaska with seafloor spreading orthogonal to the fossil spreading axis preserved in the central Canada Basin (phase 2). In this model, the Chukchi Borderland is part of Arctic Alaska.

Alaska, Northwest Territories, Yukon Territory

Arctic deep-water ferromanganese-oxide deposits reflect the unique characteristics of the Arctic Ocean

Little is known about marine mineral deposits in the Arctic Ocean, an ocean dominated by continental shelf and basins semi-closed to deep-water circulation. Here, we present data for ferromanganese crusts and nodules collected from the Amerasia Arctic Ocean in 2008, 2009, and 2012 (HLY0805, HLY0905, HLY1202). We determined mineral and chemical compositions of the crusts and nodules and the onset of their formation. Water column samples from the GEOTRACES program were analyzed for dissolved and particulate scandium concentrations, an element uniquely enriched in these deposits. The Arctic crusts and nodules are characterized by unique mineral and chemical compositions with atypically high growth rates, detrital contents, Fe/Mn ratios, and low Si/Al ratios, compared to deposits found elsewhere. High detritus reflects erosion of submarine outcrops and North America and Siberia cratons, transport by rivers and glaciers to the sea, and distribution by sea ice, brines, and currents. Uniquely high Fe/Mn ratios are attributed to expansive continental shelves, where diagenetic cycling releases Fe to bottom waters, and density flows transport shelf bottom water to the open Arctic Ocean. Low Mn contents reflect the lack of a mid-water oxygen minimum zone that would act as a reservoir for dissolved Mn. The potential host phases and sources for elements with uniquely high contents are discussed with an emphasis on scandium. Scandium sorption onto Fe oxyhydroxides and Sc-rich detritus account for atypically high scandium contents. The opening of Fram Strait in the Miocene and ventilation of the deep basins initiated Fe-Mn crust growth ∼15 Myr ago.

Geochemistry, Geophysics, Geosystems

Numerical models of pore pressure and stress changes along basement faults due to wastewater injection: Applications to the 2014 Milan, Kansas Earthquake

We have developed groundwater flow models to explore the possible relationship between wastewater injection and the 12 November 2014 M w 4.8 Milan, Kansas earthquake. We calculate pore pressure increases in the uppermost crust using a suite of models in which hydraulic properties of the Arbuckle Formation and the Milan earthquake fault zone, the Milan earthquake hypocenter depth, and fault zone geometry are varied. Given pre‐earthquake injection volumes and reasonable hydrogeologic properties, significantly increasing pore pressure at the Milan hypocenter requires that most flow occur through a conductive channel (i.e., the lower Arbuckle and the fault zone) rather than a conductive 3‐D volume. For a range of reasonable lower Arbuckle and fault zone hydraulic parameters, the modeled pore pressure increase at the Milan hypocenter exceeds a minimum triggering threshold of 0.01 MPa at the time of the earthquake. Critical factors include injection into the base of the Arbuckle Formation and proximity of the injection point to a narrow fault damage zone or conductive fracture in the pre‐Cambrian basement with a hydraulic diffusivity of about 3–30 m 2 /s. The maximum pore pressure increase we obtain at the Milan hypocenter before the earthquake is 0.06 MPa. This suggests that the Milan earthquake occurred on a fault segment that was critically stressed prior to significant wastewater injection in the area. Given continued wastewater injection into the upper Arbuckle in the Milan region, assessment of the middle Arbuckle as a hydraulic barrier remains an important research priority.

Kansas

Geologic history of the summit of Axial Seamount, Juan de Fuca Ridge

Multibeam (1 m resolution) and side scan data collected from an autonomous underwater vehicle, and lava samples, radiocarbon-dated sediment cores, and observations of flow contacts collected by remotely operated vehicle were combined to reconstruct the geologic history and flow emplacement processes on Axial Seamount's summit and upper rift zones. The maps show 52 post-410 CE lava flows and 20 precaldera lava flows as old as 31.2 kyr, the inferred age of the caldera. Clastic deposits 1&ndash;2 m thick accumulated on the rims postcaldera. Between 31 ka and 410 CE, there are no known lava flows near the summit. The oldest postcaldera lava (410 CE) is a pillow cone SE of the caldera. Two flows erupted on the W rim between &sim;800 and 1000 CE. From 1220 to 1300 CE, generally small eruptions of plagioclase phyric, depleted, mafic lava occurred in the central caldera and on the east rim. Larger post-1400 CE eruptions produced inflated lobate flows of aphyric, less-depleted, and less mafic lava on the upper rift zones and in the N and S caldera. All caldera floor lava flows, and most uppermost rift zone flows, postdate 1220 CE. Activity shifted from the central caldera to the upper S rift outside the caldera, to the N rift and caldera floor, and then to the S caldera and uppermost S rift, where two historical eruptions occurred in 1998 and 2011. The average recurrence interval deduced from the flows erupted over the last 800 years is statistically identical to the 13 year interval between historical eruptions.

Geochemistry, Geophysics, Geosystems