Geology ReportsSearch

Geology topics

Heather M. Wright

Publications and source records attributed to Heather M. Wright.

33 records · Page 2Linked to original sources

Thermal and petrologic constraints on lower crustal melt accumulation under the Salton Sea Geothermal Field

In the Salton Sea region of southern California (USA), concurrent magmatism, extension, subsidence, and sedimentation over the past 0.5 to 1.0 Ma have led to the creation of the Salton Sea Geothermal Field (SSGF)—the second largest and hottest geothermal system in the continental United States—and the small-volume rhyolite eruptions that created the Salton Buttes. In this study, we determine the flux of mantle-derived basaltic magma that would be required to produce the elevated average heat flow and sustain the magmatic roots of rhyolite volcanism observed at the surface of the Salton Sea region. We use a 2D thermal model to show that a lower-crustal, partially molten mush containing &lt; 20 &#x2013; 40 % "> <20–40% interstitial melt develops over a &#x223C; 10 5 "> ∼105 -yr timescale for basalt fluxes of 0.008 to 0.010 m 3 / m 2 / yr "> 0.010m3/m2/yr (∼0.0008 to ∼0.001 km / 3 yr "> /3yr injection rate) given extension rates at or below the current value of ∼0.01 m/yr ( Brothers et al., 2009 ). These regions of partial melt are a natural consequence of a thermal regime that scales with average surface heat flow in the Salton Trough, and are consistent with seismic observations. Our results indicate limited melting and assimilation of pre-existing rocks in the lower crust. Instead, we find that basalt fractionation in the lower crust produces derivative melts of andesitic to dacitic composition. Such melts are then expected to ascend and accumulate in the upper crust, where they further evolve to give rise to small-volume rhyolite eruptions (Salton Buttes) and fuel local spikes in surface heat flux as currently seen in the SSGF. Such upper crustal magma evolution, with limited assimilation of hydrothermally altered material, is required to explain the slight decrease in &#x3B4; 18 O "> δ18O values of zircons (and melts) that have been measured in these rhyolites.

California

Surface morphology of caldera-forming eruption deposits revealed by lidar mapping of Crater Lake National Park, Oregon- Implications for emplacement and surface modification

Large explosive eruptions of silicic magma can produce widespread pumice fall, extensive ignimbrite sheets, and collapse calderas. The surfaces of voluminous ignimbrites are rarely preserved or documented because most terrestrial examples are heavily vegetated, or severely modified by post-depositional processes. Much research addresses the internal sedimentary characteristics, flow processes, and depositional mechanisms of ignimbrites, however, surface features of ignimbrites are less well documented and understood, except for comparatively small-volume deposits of historical eruptions. The ~7,700 calendar year B.P. climactic eruption of Mount Mazama, USA vented ~50 km3 of magma, deposited first as rhyodacite pumice fall and then as a zoned rhyodacite-to-andesite ignimbrite as Crater Lake caldera collapsed. Lidar collected during summer 2010 reveals the remarkably well-preserved surface of the Mazama ignimbrite and related deposits surrounding Crater Lake caldera in unprecedented detail despite forest cover. The ±1 m lateral and ±4 cm vertical resolution lidar allows surface morphologies to be classified. Surface morphologies are created by internal depositional processes and can point to the processes at work when pyroclastic flows come to rest. We describe nine surface features including furrow-ridge sets and wedge-shaped mounds in pumice fall eroded by high-energy pyroclastic surges, flow- parallel ridges that record the passage of multiple pyroclastic flows, perched benches of marginal deposits stranded by more-mobile pyroclastic-flow cores, hummocks of dense clasts interpreted as lag deposit, transverse ridges that mark the compression and imbrication of flows as they came to rest, scarps indicating ignimbrite remobilization, fields of pit craters caused by phreatic explosions, fractures and cracks caused by extensional processes resulting from ignimbrite volume loss, and stream channels eroded in the newly formed surface. The nine morphologies presented here illustrate a dynamic depositional environment that varied spatially and with time during the eruption, and show that multiple processes modified the ignimbrite after deposition, both during and after the eruption.

Oregon

The role of crystallization-driven exsolution on the sulfur mass balance in volcanic arc magmas

The release of large amounts of sulfur to the stratosphere during explosive eruptions affects the radiative balance in the atmosphere and consequentially impacts climate for up to several years after the event. Quantitative estimations of the processes that control the mass balance of sulfur between melt, crystals, and vapor bubbles is needed to better understand the potential sulfur yield of individual eruption events and the conditions that favor large sulfur outputs to the atmosphere. The processes that control sulfur partitioning in magmas are (1) exsolution of volatiles (dominantly H 2 O) during decompression (first boiling) and during isobaric crystallization (second boiling), (2) the crystallization and breakdown of sulfide or sulfate phases in the magma, and (3) the transport of sulfur-rich vapor (gas influx) from deeper unerupted regions of the magma reservoir. Vapor exsolution and the formation/breakdown of sulfur-rich phases can all be considered as closed-system processes where mass balance arguments are generally easier to constrain, whereas the contribution of sulfur by vapor transport (open system process) is more difficult to quantify. The ubiquitous “excess sulfur” problem, which refers to the much higher sulfur mass released during eruptions than what can be accounted for by amount of sulfur originally dissolved in erupted melt, as estimated from melt inclusion sulfur concentrations (the “petrologic estimate”), reflects the challenges in closing the sulfur mass balance between crystals, melt, and vapor before and during a volcanic eruption. In this work, we try to quantify the relative importance of closed- and open-system processes for silicic arc volcanoes using kinetic models of sulfur partitioning during exsolution. Our calculations show that crystallization-induced exsolution (second boiling) can generate a significant fraction of the excess sulfur observed in crystal-rich arc magmas. This result does not negate the important role of vapor migration in sulfur mass balance but rather points out that second boiling (in situ exsolution) can provide the necessary yield to drive the excess sulfur to the levels observed for crystal-rich systems. In contrast, in crystal-poor systems, magma recharge that releases sulfur-rich bubbles is necessary and most likely the primary contributor to sulfur mass balance. Finally, we apply our model to account for the effect of sulfur partitioning during second boiling and its impact on sulfur released during the Cerro Galan supereruption in Argentina (2.08 Ma) and show the potential importance of second boiling in releasing a large amount of sulfur to the atmosphere during the eruption of large crystal-rich ignimbrites.

Journal of Geophysical Research B: Solid Earth

Juvenile magma recognition and eruptive dynamics inferred from the analysis of ash time series: The 2015 reawakening of Cotopaxi volcano

Forecasting future activity and performing hazard assessments during the reactivation of volcanoes remain great challenges for the volcanological community. On August 14, 2015 Cotopaxi volcano erupted for the first time in 73 years after approximately four months of precursory activity, which included an increase in seismicity, gas emissions, and minor ground deformation. Here we discuss the use of near real-time petrological monitoring of ash samples as a complementary aid to geophysical monitoring, in order to infer eruption dynamics and evaluate possible future eruptive activity at Cotopaxi. Twenty ash samples were collected between August 14 and November 23, 2015 from a monitoring site on the west flank of the volcano. These samples contain a range of grain types that we classified as: hydrothermal/altered, lithic, juvenile, and free crystals. The relative proportions of theses grains evolved as the eruption progressed, with increasing amounts of juvenile material and a decrease in hydrothermally altered material. In samples from the initial explosion, juvenile grains are glassy, microlite-poor and contain hydrothermal minerals (opal and alunite). The rising magma came in contact with the hydrothermal system under confinement, causing hydro-magmatic explosions that cleared the upper part of the plumbing system. Subsequently, the magmatic column produced a thermal aureole in the conduit and dried out the hydrothermal system, allowing for dry eruptions. Magma ascent rates were low enough to allow for efficient outgassing and microlite growth. Constant supply of magma from below caused quasi-continuous disruption of the uppermost magma volume through a combination of shear-deformation and gas expansion. The combination of increasing crystallinity of juvenile grains, and high measured SO 2 flux indicate decreasing integrated magma ascent rates and clearing of the hydrothermal system along transport pathways in a system open to gas loss. The near real-time monitoring of ash samples combined with traditional geophysical monitoring techniques during the reawakening of Cotopaxi allowed us to gain a much clearer understanding of events than when using traditional geophysical monitoring alone.

Cotopaxi volcano

Episodic Holocene eruption of the Salton Buttes rhyolites, California, from paleomagnetic, U-Th, and Ar/Ar dating

In the Salton Trough, CA, five rhyolite domes form the Salton Buttes: Mullet Island, Obsidian Butte, Rock Hill, North and South Red Hill, from oldest to youngest. Results presented here include 40 Ar/ 39 Ar anorthoclase ages, 238 U- 230 Th zircon crystallization ages, and comparison of remanent paleomagnetic directions with the secular variation curve, which indicate that all domes are Holocene. 238 U- 230 Th zircon crystallization ages are more precise than but within uncertainty of 40 Ar/ 39 Ar anorthoclase ages, suggesting that zircon crystallization proceeded until shortly before eruption in all cases except one. Remanent paleomagnetic directions require three eruption periods: (1) Mullet Island, (2) Obsidian Butte, and (3) Rock Hill, North Red Hill, and South Red Hill. Borehole cuttings logs document up to two shallow tephra layers. North and South Red Hills likely erupted within 100 years of each other, with a combined 238 U- 230 Th zircon isochron age of: 2.83 &plusmn; 0.60 ka (2 sigma); paleomagnetic evidence suggests this age predates eruption by hundreds of years (1800 cal BP). Rock Hill erupted closely in time to these eruptions. The Obsidian Butte 238 U- 230 Th isochron age (2.86 &plusmn; 0.96 ka) is nearly identical to the combined Red Hill age, but its Virtual Geomagnetic Pole position suggests a slightly older age. The age of aphyric Mullet Island dome is the least well constrained: zircon crystals are resorbed and the paleomagnetic direction is most distinct; possible Mullet Island ages include ca. 2300, 5900, 6900, and 7700 cal BP. Our results constrain the duration of Salton Buttes volcanism to between ca. 5900 and 500 years.

Geochemistry, Geophysics, Geosystems

The effect of pressurized magma chamber growth on melt migration and pre-caldera vent locations through time at Mount Mazama, Crater Lake, Oregon

The pattern of eruptions at long-lived volcanic centers provides a window into the co-evolution of crustal magma transport, tectonic stresses, and unsteady magma generation at depth. Mount Mazama in the Oregon Cascades has seen variable activity over the last 400 ky, including the 50 km 3 climactic eruption at ca. 7.7 ka that produced Crater Lake caldera. The physical mechanisms responsible for the assembly of silicic magma reservoirs that are the precursors to caldera-forming eruptions are poorly understood. Here we argue that the spatial and temporal distribution of geographically clustered volcanic vents near Mazama reflects the development of a centralized magma chamber that fed the climactic eruption. Time-averaged eruption rates at Mount Mazama imply an order of magnitude increase in deep magma influx prior to the caldera-forming event, suggesting that unsteady mantle melting triggered a chamber growth episode that culminated in caldera formation. We model magma chamber&ndash;dike interactions over &sim;50 ky preceding the climactic eruption to fit the observed distribution of surface eruptive vents in space and time, as well as petrologically estimated deep influx rates. Best fitting models predict an expanding zone of dike capture caused by a growing, oblate spheroidal magma chamber with 10&ndash;30 MPa of overpressure. This growing zone of chamber influence causes closest approaching regional mafic vent locations as well as more compositionally evolved Mazama eruptions to migrate away from the climactic eruptive center, returning as observed to the center after the chamber drains during the caldera-forming eruption.

Oregon

Magma fracturing and degassing associated with obsidian formation: The explosive–effusive transition

This paper explores the role of melt fracturing in degassing rhyolitic volcanic systems. The Monte Pilato-Rocche Rosse eruptions in Italy evolved from explosive to effusive in style, and H 2 O content in quenched glasses changed over time from relatively H 2 O-rich (~ 0.90 wt.%) to H 2 O-poor dense obsidian (~ 0.10&ndash;0.20 wt.%). In addition, healed fractures have been recorded in all different eruptive materials, from the glass of early-erupted tube pumice and rinds of breadcrusted obsidian pyroclasts, to the glass of late-erupted dense obsidian pyroclasts, and throughout the final effusive Rocche Rosse lava flow. These rocks show multiple fault sets, some with crenulated fault planes indicating resumption of viscous flow after faulting, complex obsidian breccias with evidence for post-brecciation folding and stretching, and centimetre- to metre-thick tuffisite preserved in pyroclasts and lava, representing collapsed foam due to fracturing of vesicle walls. These microstructural observations indicate that multiple fracturing and healing events occurred during both explosive and effusive eruptions. H 2 O content in glass decreases by as much as 0.14 wt.% towards healed fractures/faults and decreases in stretched obsidian breccias towards regions of intense brecciation. A drop in pressure and/or increase in temperature along fractures caused diffusive H 2 O migration through melt towards fracture surfaces. Repetitive and pervasive fracturing and healing thereby create conditions for diffusive H 2 O loss into fractures and subsequent escape through permeable paths. This type of progressive magma degassing provides a potential mechanism to explain the formation of dense obsidian and the evolution from explosive to effusive eruption style.

Monte Pilato-Rocche Rosse

Pyroclast textural variation as an indicator of eruption column steadiness in andesitic Plinian eruptions at Mt. Ruapehu

Between 27 and 11 cal. ka BP, a transition is observed in Plinian eruptions at Mt. Ruapehu, indicating evolution from non-collapsing (steady and oscillatory) eruption columns to partially collapsing columns (both wet and dry). To determine the causes of these variations over this eruptive interval, we examined lapilli fall deposits from four eruptions representing the climactic phases of each column type. All eruptions involve andesite to basaltic andesite magmas containing plagioclase, clinopyroxene, orthopyroxene and magnetite phenocrysts. Differences occur in the dominant pumice texture, the degree of bulk chemistry and textural variability, the average microcrystallinity and the composition of groundmass glass. In order to investigate the role of ascent and degassing processes on column stability, vesicle textures were quantified by gas volume pycnometry (porosity), X-ray synchrotron and computed microtomography (μ-CT) imagery from representative clasts from each eruption. These data were linked to groundmass crystallinity and glass geochemistry. Pumice textures were classified into six types (foamy, sheared, fibrous, microvesicular, microsheared and dense) according to the vesicle content, size and shape and microlite content. Bulk porosities vary from 19 to 95 % among all textural types. Melt-referenced vesicle number density ranges between 1.8 × 10 2 and 8.9 × 10 2 mm −3 , except in fibrous textures, where it spans from 0.3 × 10 2 to 53 × 10 2 mm −3 . Vesicle-free magnetite number density varies within an order of magnitude from 0.4 × 10 2 to 4.5 × 10 2 mm −3 in samples with dacitic groundmass glass and between 0.0 and 2.3 × 10 2 mm −3 in samples with rhyolitic groundmass. The data indicate that columns that collapsed to produce pyroclastic flows contained pumice with the greatest variation in bulk composition (which overlaps with but extends to slightly more silicic compositions than other eruptive products); textures indicating heterogeneous bubble nucleation, progressively more complex growth history and shear-localization; and the highest degrees of microlite crystallization, most evolved melt compositions and lowest relative temperatures. These findings suggest that collapsing columns in Ruapehu have been produced when strain localization is prominent, early bubble nucleation occurs and variation in decompression rate across the conduit is greatest. This study shows that examination of pumice from steady phases that precede column collapse may be used to predict subsequent column behaviour.

Mt. Ruapehu

Compaction and gas loss in welded pyroclastic deposits as revealed by porosity, permeability, and electrical conductivity measurements of the Shevlin Park Tuff

Pyroclastic flows produced by large volcanic eruptions commonly densify after emplacement. Processes of gas escape, compaction, and welding in pyroclastic-flow deposits are controlled by the physical and thermal properties of constituent material. Through measurements of matrix porosity, permeability, and electrical conductivity, we provide a framework for understanding the evolution of pore structure during these processes. Using data from the Shevlin Park Tuff in central Oregon, United States, and from the literature, we find that over a porosity range of 0%–70%, matrix permeability varies by almost 10 orders of magnitude (from 10 –20 to 10 –11 m 2 ), with over three orders of magnitude variation at any given porosity. Part of the variation at a given porosity is due to permeability anisotropy, where oriented core samples indicate higher permeabilities parallel to foliation (horizontally) than perpendicular to foliation (vertically). This suggests that pore space is flattened during compaction, creating anisotropic crack-like networks, a geometry that is supported by electrical conductivity measurements. We find that the power law equation: k 1 = 1.3 × 10 –21 × ϕ 5.2 provides the best approximation of dominant horizontal gas loss, where k 1 = permeability, and ϕ = porosity. Application of Kozeny-Carman fluid-flow approximations suggests that permeability in the Shevlin Park Tuff is controlled by crack- or disk-like pore apertures with minimum widths of 0.3 and 7.5 μm. We find that matrix permeability limits compaction over short times, but deformation is then controlled by competition among cooling, compaction, water resorption, and permeable gas escape. These competing processes control the potential for development of overpressure (and secondary explosions) and the degree of welding in the deposit, processes that are applicable to viscous densification of volcanic deposits in general. Further, the general relationships among porosity, permeability, and pore geometry are relevant for flow of any fluid through an ignimbritic host.

Oregon

Sixty thousand years of magmatic volatile history before the caldera-forming eruption of Mount Mazama, Crater Lake, Oregon

The well-documented eruptive history of Mount Mazama, Oregon, provides an excellent opportunity to use pre-eruptive volatile concentrations to study the growth of an explosive silicic magmatic system. Melt inclusions (MI) hosted in pyroxene and plagioclase crystals from eight dacitic–rhyodacitic eruptive deposits (71–7.7 ka) were analyzed to determine variations in volatile-element concentrations and changes in magma storage conditions leading up to and including the climactic eruption of Crater Lake caldera. Temperatures (Fe–Ti oxides) increased through the series of dacites, then decreased, and increased again through the rhyodacites (918–968 to ~950 to 845–895 °C). Oxygen fugacity began at nickel–nickel-oxide buffer (NNO) +0.8 (71 ka), dropped slightly to NNO +0.3, and then climbed to its highest value with the climactic eruption (7.7 ka) at NNO +1.1 log units. In parallel with oxidation state, maximum MI sulfur concentrations were high early in the eruptive sequence (~500 ppm), decreased (to ~200 ppm), and then increased again with the climactic eruption (~500 ppm). Maximum MI sulfur correlates with the Sr content (as a proxy for LREE, Ba, Rb, P 2 O 5 ) of recharge magmas, represented by basaltic andesitic to andesitic enclaves and similar-aged lavas. These results suggest that oxidized Sr-rich recharge magmas dominated early and late in the development of the pre-climactic dacite–rhyodacite system. Dissolved H 2 O concentrations in MI do not, however, correlate with these changes in dominant recharge magma, instead recording vapor solubility relations in the developing shallow magma storage and conduit region. Dissolved H 2 O concentrations form two populations through time: the first at 3–4.6 wt% (with a few extreme values up to 6.1 wt%) and the second at ≤2.4 wt%. CO 2 concentrations measured in a subset of these inclusions reach up to 240 ppm in early-erupted deposits (71 ka) and are below detection in climactic deposits (7.7 ka). Combined H 2 O and CO 2 concentrations and solubility models indicate a dominant storage region at 4–7 km (up to 12 km), with drier inclusions that diffusively re-equilibrated and/or were trapped at shallower depths. Boron and Cl (except in the climactic deposit) largely remained in the melt, suggesting vapor–melt partition coefficients and gas fractions were low. Modeled Li, F, and S vapor–melt partition coefficients are higher than those of B and Cl. The decrease in maximum MI CO 2 concentration following the earliest dacitic eruptions is interpreted to result from a broadening of the shallow storage region to greater than the diameter of subjacent feeders, so that greater proportions of reservoir magma were to the side of CO 2 -bearing vapor bubbles ascending vertically from the locus of recharge magma injection, thereby escaping recarbonation by streaming vapor bubbles. The Mazama melt inclusions provide a picture of a growing magma storage region, where chemical variations in melt and magma occur due to changes in the nature and supply rate of magma recharge, the timing of degassing, and the possible degree of equilibration with gases from below.

Oregon

Columnar jointing in vapor-phase-altered, non-welded Cerro Galán Ignimbrite, Paycuqui, Argentina

Columnar jointing is thought to occur primarily in lavas and welded pyroclastic flow deposits. However, the non-welded Cerro Gal&aacute;n Ignimbrite at Paycuqui, Argentina, contains well-developed columnar joints that are instead due to high-temperature vapor-phase alteration of the deposit, where devitrification and vapor-phase crystallization have increased the density and cohesion of the upper half of the section. Thermal remanent magnetization analyses of entrained lithic clasts indicate high emplacement temperatures, above 630&deg;C, but the lack of welding textures indicates temperatures below the glass transition temperature. In order to remain below the glass transition at 630&deg;C, the minimum cooling rate prior to deposition was 3.0&thinsp;&times;&thinsp;10 &minus;3 &ndash;8.5&thinsp;&times;&thinsp;10 &minus;2 &deg;C/min (depending on the experimental data used for comparison). Alternatively, if the deposit was emplaced above the glass transition temperature, conductive cooling alone was insufficient to prevent welding. Crack patterns (average, 4.5 sides to each polygon) and column diameters (average, 75 cm) are consistent with relatively rapid cooling, where advective heat loss due to vapor fluxing increases cooling over simple conductive heat transfer. The presence of regularly spaced, complex radiating joint patterns is consistent with fumarolic gas rise, where volatiles originated in the valley-confined drainage system below. Joint spacing is a proxy for cooling rates and is controlled by depositional thickness/valley width. We suggest that the formation of joints in high-temperature, non-welded deposits is aided by the presence of underlying external water, where vapor transfer causes crystallization in pore spaces, densifies the deposit, and helps prevent welding.

Cerro Galán

A re-appraisal of the stratigraphy and volcanology of the Cerro Galán volcanic system, NW Argentina

From detailed fieldwork and biotite 40 Ar/ 39 Ar dating correlated with paleomagnetic analyses of lithic clasts, we present a revision of the stratigraphy, areal extent and volume estimates of ignimbrites in the Cerro Gal&aacute;n volcanic complex. We find evidence for nine distinct outflow ignimbrites, including two newly identified ignimbrites in the Toconquis Group (the Pitas and Vega Ignimbrites). Toconquis Group Ignimbrites (~5.60&ndash;4.51 Ma biotite ages) have been discovered to the southwest and north of the caldera, increasing their spatial extents from previous estimates. Previously thought to be contemporaneous, we distinguish the Real Grande Ignimbrite (4.68&thinsp;&plusmn;&thinsp;0.07 Ma biotite age) from the Cueva Negra Ignimbrite (3.77&thinsp;&plusmn;&thinsp;0.08 Ma biotite age). The form and collapse processes of the Cerro Gal&aacute;n caldera are also reassessed. Based on re-interpretation of the margins of the caldera, we find evidence for a fault-bounded trapdoor collapse hinged along a regional N-S fault on the eastern side of the caldera and accommodated on a N-S fault on the western caldera margin. The collapsed area defines a roughly isosceles trapezoid shape elongated E-W and with maximum dimensions 27&thinsp;&times;&thinsp;16 km. The Cerro Gal&aacute;n Ignimbrite (CGI; 2.08&thinsp;&plusmn;&thinsp;0.02 Ma sanidine age) outflow sheet extends to 40 km in all directions from the inferred structural margins, with a maximum runout distance of ~80 km to the north of the caldera. New deposit volume estimates confirm an increase in eruptive volume through time, wherein the Toconquis Group Ignimbrites increase in volume from the ~10 km 3 Lower Merihuaca Ignimbrite to a maximum of ~390 km 3 (Dense Rock Equivalent; DRE) with the Real Grande Ignimbrite. The climactic CGI has a revised volume of ~630 km 3 (DRE), approximately two thirds of the commonly quoted value.

Cerro Galán

Melt fracturing and healing: A mechanism for degassing and origin of silicic obsidian

We present water content transects across a healed fault in pyroclastic obsidian from Lami pumice cone, Lipari, Italy, using synchrotron Fourier transform infrared spectroscopy. Results indicate that rhyolite melt degassed through the fault surface. Transects define a trough of low water content coincident with the fault trace, surrounded on either side by high-water-content plateaus. Plateaus indicate that obsidian on either side of the fault equilibrated at different pressure-temperature ( P - T ) conditions before being juxtaposed. The curves into the troughs indicate disequilibrium and water loss through diffusion. If we assume constant T, melt equilibrated at pressures differing by 0.74 MPa before juxtaposition, and the fault acted as a low- P permeable path for H 2 O that diffused from the glass within time scales of 10 and 30 min. Assuming constant P instead, melt on either side could have equilibrated at temperatures differing by as much as 100 °C, before being brought together. Water content on the fault trace is particularly sensitive to post-healing diffusion. Its preserved value indicates either higher temperature or lower pressure than the surroundings, indicative of shear heating and dynamic decompression. Our results reveal that water contents of obsidian on either side of the faults equilibrated under different P - T conditions and were out of equilibrium with each other when they were juxtaposed due to faulting immediately before the system was quenched. Degassing due to faulting could be linked to cyclical seismic activity and general degassing during silicic volcanic activity, and could be an efficient mechanism of producing low-water-content obsidian.

Geology

1:100,000-scale topographic contours derived from digital elevation models, San Francisco Bay region, California: a digital database

This report presents a consistent set of 1:100,000-scale vector topographic contours for all eleven 30x60-minute quadrangles in the San Francisco Bay region for use in visualizing the topography and preparing maps of the region. The contours were prepared by contouring an areally continuous 30-m altitude grid (National Elevation Dataset, Jan., 1999), and differ from USGS hypsographic DLG's (available for only part of the region). The report consists of 26 numbered parts, which represent text, spatial data, and 1:100,000-scale map graphics. Most of the files are provided in two or three different digital formats. All files are available for download here.

California

Preliminary maps of Quaternary deposits and liquefaction susceptibility, nine-county San Francisco Bay region, California: a digital database

This report presents a preliminary map and database of Quaternary deposits and liquefaction susceptibility for the nine-county San Francisco Bay region, together with a digital compendium of ground effects associated with past earthquakes in the region. The report consists of (1) a spatial database of fivedata layers (Quaternary deposits, quadrangle index, and three ground effects layers) and two text layers (a labels and leaders layer for Quaternary deposits and for ground effects), (2) two small-scale colored maps (Quaternary deposits and liquefaction susceptibility), (3) a text describing the Quaternary map, liquefaction interpretation, and the ground effects compendium, and (4) the databse description pamphlet. The nine counties surrounding San Francisco Bay straddle the San Andreas fault system, which exposes the region to serious earthquake hazard (Working Group on California Earthquake Probabilities, 1999). Much of the land adjacent to the Bay and the major rivers and streams is underlain by unconsolidated deposits that are particularly vulnerable to earthquake shaking and liquefaction of water-saturated granular sediment. This new map provides a modern and regionally consistent treatment of Quaternary surficial deposits that builds on the pioneering mapping of Helley and Lajoie (Helley and others, 1979) and such intervening work as Atwater (1982), Helley and others (1994), and Helley and Graymer (1997a and b). Like these earlier studies, the current mapping uses geomorphic expression, pedogenic soils, and inferred depositional environments to define and distinguish the map units. In contrast to the twelve map units of Helley and Lajoie, however, this new map uses a complex stratigraphy of some forty units, which permits a more realistic portrayal of the Quaternary depositional system. The two colored maps provide a regional summary of the new mapping at a scale of 1:275,000, a scale that is sufficient to show the general distribution and relationships of the map units but cannot distinguish the more detailed elements that are present in the database. The report is the product of years of cooperative work by the USGS National Earthquake Hazards Reduction Program (NEHRP) and National Cooperative Geologic Mapping Program, William Lettis and & Associates, Inc. (WLA) and, more recently, by the California Division of Mines and Geology as well. An earlier version was submitted to the Geological Survey by WLA as a final report for a NEHRP grant (Knudsen and others, 2000). The mapping has been carried out by WLA geologists under contract to the NEHRP Earthquake Program (Grants #14-08-0001-G2129, 1434-94-G-2499, 1434-HQ-97-GR-03121, and 99-HQ-GR-0095) and with other limited support from the County of Napa, and recently also by the California Division of Mines and Geology. The current map consists of this new mapping and revisions of previous USGS mapping.

Open-File Report