Search USGSSearch

SEARCH · Search USGS

Results for “Geochemistry Geophysics 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 109 records · Page 6Linked to original sources

Copper-nickel-rich, amalgamated ferromanganese crust-nodule deposits from Shatsky Rise, NW Pacific

A unique set of ferromanganese crusts and nodules collected from Shatsky Rise (SR), NW Pacific, were analyzed for mineralogical and chemical compositions, and dated using Be isotopes and cobalt chronometry. The composition of these midlatitude, deep-water deposits is markedly different from northwest-equatorial Pacific (PCZ) crusts, where most studies have been conducted. Crusts and nodules on SR formed in close proximity and some nodule deposits were cemented and overgrown by crusts, forming amalgamated deposits. The deep-water SR crusts are high in Cu, Li, and Th and low in Co, Te, and Tl concentrations compared to PCZ crusts. Thorium concentrations (ppm) are especially striking with a high of 152 (mean 56), compared to PCZ crusts (mean 11). The deep-water SR crusts show a diagenetic chemical signal, but not a diagenetic mineralogy, which together constrain the redox conditions to early oxic diagenesis. Diagenetic input to crusts is rare, but unequivocal in these deep-water crusts. Copper, Ni, and Li are strongly enriched in SR deep-water deposits, but only in layers older than about 3.4 Ma. Diagenetic reactions in the sediment and dissolution of biogenic calcite in the water column are the likely sources of these metals. The highest concentrations of Li are in crust layers that formed near the calcite compensation depth. The onset of Ni, Cu, and Li enrichment in the middle Miocene and cessation at about 3.4 Ma were accompanied by changes in the deep-water environment, especially composition and flow rates of water masses, and location of the carbonate compensation depth.

Geochemistry, Geophysics, Geosystems

Temporal variations of geyser water chemistry in the Upper Geyser Basin, Yellowstone National Park, USA

Geysers are rare features that reflect a delicate balance between an abundant supply of water and heat and a unique geometry of fractures and porous rocks. Between April 2007 and September 2008, we sampled Old Faithful, Daisy, Grand, Oblong, and Aurum geysers in Yellowstone National Park's Upper Geyser Basin and characterized temporal variations in major element chemistry and water isotopes (δ 18 O, δD, 3 H). We compare these temporal variations with temporal trends of Geyser Eruption Intervals (GEI). SiO 2 concentrations and geothermometry indicate that the geysers are fed by waters ascending from a reservoir with temperatures of ∼190 to 210°C. The studied geysers display small and complex chemical and isotopic seasonal variations, and geysers with smaller volume display larger seasonal variations than geysers with larger volumes. Aurum and Oblong Geysers contain detectable tritium concentrations, suggesting that erupted water contains some modern meteoric water. We propose that seasonal GEI variations result from varying degrees of evaporation, meteoric water recharge, water table fluctuations, and possible hydraulic interaction with the adjacent Firehole River. We demonstrate that the concentrations of major dissolved species in Old Faithful Geyser have remained nearly constant since 1884 despite large changes in Old Faithful's eruption intervals, suggesting that no major changes have occurred in the hydrothermal system of the Upper Geyser Basin for >120 years. Our data set provides a baseline for monitoring future changes in geyser activity that might result from varying climate, earthquakes, and changes in heat flow from the underlying magmatic system.

Wyoming

Large-scale splay faults on a strike-slip fault system: The Yakima Folds, Washington State

The Yakima Folds (YF) comprise anticlines above reverse faults cutting flows of the Miocene Columbia River Basalt Group of central Washington State. The YF are bisected by the ~1100-km-long Olympic-Wallowa Lineament (OWL), which is an alignment of topographic features including known faults. There is considerable debate about the origin and earthquake potential of both the YF and OWL, which lie near six major dams and a large nuclear waste storage site. Here I show that the trends of the faults forming the YF relative to the OWL match remarkably well the trends of the principal stress directions at the end of a vertical strike-slip fault. This comparison and the termination of some YF against the OWL are consistent with the YF initially forming as splay faults caused by an along-strike decrease in the amount of strike-slip on the OWL. The hypothesis is that the YF faults initially developed as splay faults in the early to mid Miocene under NNW-oriented principal compressive stress, but the anticlines subsequently grew with thrust motion after the principal compressive stress direction rotated to N-S or NNE after the mid-Miocene. A seismic profile across one of the YF anticlines shows folding at about 7 km depth, indicating deformation of sub-basalt strata. The seismic profile and the hypothesized relationship between the YF and the OWL suggest that the structures are connected in the middle or lower crust, and that the faults forming the YF are large-scale splay faults associated with a major strike-slip fault system.

Washington

Origin of primitive tholeiitic and calc-alkaline basalts at Newberry Volcano, Oregon

Major and trace element and radiogenic isotopic characteristics of primitive mafic Pleistocene and Holocene lavas from Newberry Volcano, Oregon, define two groups. The first consists of dry tholeiitic high-alumina olivine basalts that are slightly enriched in highly incompatible elements. The second group consists of calc-alkaline basalts that contained 2–4 wt % H 2 O prior to eruption and shows strong enrichment in the light rare earth elements, Ba, and Sr, and deficits in Nb, Ta, Hf, and Zr. The tholeiitic basalts reflect 6–11% anhydrous adiabatic decompression melting of spinel peridotite. The calc-alkaline basalts derived from compositionally distinct sources with strong LIL enrichment and relative depletion in HFSE, but with Sr, Nd, Hf, and Pb isotopic composition only slightly distinct from the sources of the tholeiitic magmas. Radiogenic Os correlates with LREE enrichment in the calc-alkaline magmas, which indicates that their source materials include a contribution from a mafic component that was melted in the garnet stability field. The calc-alkaline magmas were derived by melting of peridotite metasomatized by a fluid/melt that originated by melting of a mixture of the sediment plus MORB basalt/mantle in the underlying subducting oceanic plate. While the trace element characteristics of the calc-alkaline magmas were determined by the subduction component, their isotopic characteristics were modified during transit through the mantle by interaction with the highly magmatically processed mantle wedge beneath Newberry Volcano that, without the slab component, serves as the source of the tholeiitic magmas.

Oregon

Evolution of the 2015 Cotopaxi eruption revealed by combined geochemical & seismic observations

Through integration of multiple data streams to monitor volcanic unrest scientists are able to make more robust eruption forecast and to obtain a more holistic interpretation of volcanic systems. We examined gas emission and gas geochemistry, seismic and petrologic data recorded during the 2015 unrest of Cotopaxi (Ecuador) in order to decipher the origin and temporal evolution of this eruption. Identification of families of similar seismic events and the use of seismic amplitude ratios reveals temporal changes in volcanic processes. SO2 (300 to 24000 t/d), BrO/SO2 (5-10 x10-5), SO2/HCl (5.8 ± 4.8 and 6.6 ± 3.0) and CO2/SO2 (0.6 to 2.1) measured throughout the eruption indicate a shallow magmatic source. Bulk ash and glass chemistry indicate a homogenous andesitic (SiO2 wt%=56.94 ± 0.25) magma having undergone extensive S-exsolution and degassing during ascent. These data lead us to interpret this eruption as a magma intrusion and ascend to shallow levels. The intrusion progressively interacted with the hydrothermal system, boiled off water, and produced hydromagmatic explosions. A small volume of this intrusion continued to fragment and produced episodic ash emissions until it was sufficiently degassed and rheologically stiff. Based on the 470 kt of measured SO2 we estimate that ~ 65.3 x106 m3 of magma were required to supply the emitted gases. This volume exceeds the volume of erupted juvenile material by a factor of 50. This result emphasizes the importance of careful monitoring of Cotopaxi to identify the intrusion of a new batch of magma, which could rejuvenate the non-erupted material.

Geochemistry, Geophysics, Geosystems

Upper mantle earth structure in Africa from full-wave ambient noise tomography

Our understanding of the tectonic development of the African continent and the interplay between its geological provinces is hindered by unevenly distributed seismic instrumentation. In order to better understand the continent, we used long-period ambient noise full waveform tomography on data collected from 186 broadband seismic stations throughout Africa and surrounding regions to better image the upper mantle structure. We extracted empirical Green’s functions from ambient seismic noise using a frequency-time normalization method and retrieved coherent signal at periods of 7-340 seconds. We simulated wave propagation through a heterogeneous Earth using a spherical finite-difference approach to obtain synthetic waveforms, measured the misfit as phase delay between the data and synthetics, calculated numerical sensitivity kernels using the scattering integral approach, and iteratively inverted for structure. The resulting images of isotropic, shear wavespeed for the continent reveal segmented, low-velocity upper mantle beneath the highly magmatic northern and eastern sections of the East African Rift System (EARS). In the southern and western sections, high-velocity upper mantle dominates, and distinct, low-velocity anomalies are restricted to regions of current volcanism. At deeper depths, the southern and western EARS transitions to low-velocities. In addition to the EARS, several low-velocity anomalies are scattered through the shallow upper mantle beneath Angola and North Africa, and some of these low-velocity anomalies may be connected to a deeper feature. Distinct upper mantle high-velocity anomalies are imaged throughout the continent and suggest multiple cratonic roots within the Congo region, and possible cratonic roots within the Sahara Metacraton.

Geochemistry, Geophysics, Geosystems

Geographic attribution of soils using probabilistic modeling of GIS data for forensic search efforts

Examinations of soil traces associated with forensic evidence can be used to narrow potential source area(s) by characterizing features of the trace soil assemblage, some of which are limited to specific regions. Soil characteristics may be used to infer the likelihoods of the soil trace being derived from distinct areas within digital maps, including both maps of discrete classes such as formations on geologic maps and land cover, and continuous geospatial data, such as distance from a point source. Seldom do digital maps precisely represent the observable characteristics in a soil trace. Nevertheless, logical assigned likelihoods based on the correspondence between the mapped characteristics and the observed soil particulate assemblage permit creation of a model of the more probable sources of the soil trace. This approach is applied to a 2003 case in which forensic soil samples derived from digging tools were characterized for investigative leads and to narrow the search area of a clandestine grave. This grave site was located in 2005. The suspect traveled approximately 5,000 km before arrest, so narrowing the prioritized search area for law enforcement would be beneficial. Soil examination and case circumstances were used to assign relative likelihoods within digital maps (GIS or Geographic Information Systems data) of geology, soil mineralogy, plant distributions, power plant locations, and proximity to the known travel path. The product of these individual probability maps generates joint probability models to narrow the recommended search area. The digital model output can be easily overlaid on infrastructure maps to aid law enforcement searches.

Geochemistry, Geophysics, Geosystems

Constraining the early eruptive history of the Mono Craters rhyolites, California, based on 238U–230Th isochron dating of their explosive and effusive products

The Mono Craters are an overlapping chain of at least 28 domes and coulees located south of Mono Lake, east central California, and represent the most recent eruptions of high‐silica rhyolite magma in the Mono Lake‐Long Valley volcanic region. Regionally widespread tephra fall deposits from the Mono Craters serve as important chronostratigraphic markers for correlations of late Quaternary terrestrial archives in the western United States. A well‐resolved eruption chronology for the Mono Craters is thus not only critical for volcanic hazard considerations but also for paleoclimatic and paleomagnetic studies. Here we constrain the timing of early Mono Craters volcanism using ion microprobe 238 U‐ 230 Th dating of unpolished crystal faces of autocrystic zircon and allanite microphenocrysts, which samples the final stage of crystallization prior to eruption. The stratigraphically oldest domes (biotite‐bearing rhyolites) yield 238 U‐ 230 Th isochron dates between ca. 42 and ca. 19 ka and are unambiguously linked to dated tephras in the Wilson Creek formation of ancestral Mono Lake via coupled 238 U‐ 230 Th geochronology and titanomagnetite geochemistry. The second oldest set of domes (orthopyroxene‐bearing rhyolites) have overlapping 238 U‐ 230 Th isochron dates that are within uncertainty of published K/Ar and 40 Ar/ 39 Ar dates for the fayalite‐bearing rhyolite domes, suggesting a period of intense effusive activity for the Mono Craters near the Pleistocene‐Holocene boundary. Our new and previously published 238 U‐ 230 Th dates for tephras in the Wilson Creek formation provide robust geochronologic constraints for the controversial geomagnetic excursion recorded in the Wilson Creek formation.

California

Fault initiation in serpentinite

Serpentinite fault rheology is fundamental to tectonic and earthquake processes, yet links between deformation textures and strength evolution during fault initiation are poorly constrained. Here I present field and petrographic microstructural observations of unsheared and sheared serpentinite that demonstrate a progression of fault development. I compliment observations with a clast size distribution analysis to investigate the evolution of fault rigidity, and a numerical model to query the stress distribution of a clast-in-matrix geometry. Unsheared microstructures reveal well-aligned, elongate serpentine in the matrix and short, randomly-oriented serpentine in clasts. Sheared matrix displays cataclastic textures, discrete brittle surfaces, dissolution bands and ductile textures defined by anastomosing networks of well-aligned, fine-grained serpentine. During fault initiation, matrix serpentine anisotropy promotes slip on basal planes or fiber aggregates, and clast-on-clast interactions drive a high bulk viscosity prone to stick-slip behavior. As deformation progresses clast fracturing is focused at clast tips and smaller clasts are preferentially removed by dissolution-precipitation processes, increasing the relative abundance of matrix. Strain is continually focused in the matrix and as clast content reduces the bulk viscosity drops. This study reveals that viscosity contrasts formed by primary serpentinization textures are a major driver for the development and strength evolution of faults. On a continental-scale, similar processes may govern earthquake distributions, fault growth and segmentation patterns on young, serpentinite-hosted faults.

California

Modeling ash dispersal from future eruptions of Taupo supervolcano

Hazard analysis at caldera volcanoes is challenging due to the wide range of eruptive and environmental conditions that can plausibly occur during renewed activity. Taupo volcano, New Zealand, is a frequently active and productive rhyolitic caldera volcano that has hosted the world's youngest known supereruption and numerous smaller explosive events. To assess ashfall hazard from future eruptions, we have simulated atmospheric ash dispersal using the Ash3d model. We consider five eruption scenarios spanning magma volumes of 0.1–500 km 3 and investigate the main factors governing ash dispersal in modern atmospheric conditions. Our results are examined in the context of regional synoptic weather patterns (Kidson types) that provide a framework for assessing the variability of ashfall distribution in different wind fields. For the smallest eruptions (~0.1‐km 3 magma), ashfall thicknesses >1 cm are largely confined to the central North Island, with dispersal controlled by day‐to‐day weather and the dominance of westerly winds. With increasing eruptive volume (1–5‐km 3 magma), ashfall thicknesses >1 cm would likely reach major population centers throughout the North Island. Dispersal is less dependent on weather patterns as the formation of a radially expanding umbrella cloud forces ash upwind or crosswind, although strong stratospheric winds significantly restrict umbrella spreading. For large eruptions (50–500‐km 3 magma), powerful expansion of the umbrella cloud results in widespread ashfall at damaging thicknesses (>10 cm) across most of the North Island and top of the South Island. Synoptic climatology may prove a useful additional technique for long‐term hazard planning at caldera volcanoes.

Taupo volcano

Xenon hydrate as an analogue of methane hydrate in geologic systems out of thermodynamic equilibrium

Methane hydrate occurs naturally under pressure and temperature conditions that are not straightforward to replicate experimentally. Xenon has emerged as an attractive laboratory alternative to methane for studying hydrate formation and dissociation in multiphase systems, given that it forms hydrates under milder conditions. However, building reliable analogies between the two hydrates requires systematic comparisons, which are currently lacking. We address this gap by developing a theoretical and computational model of gas hydrates under equilibrium and nonequilibrium conditions. We first compare equilibrium phase behaviors of the Xe·H 2 O and CH 4 ·H 2 O systems by calculating their isobaric phase diagram, and then study the nonequilibrium kinetics of interfacial hydrate growth using a phase field model. Our results show that Xe·H 2 O is a good experimental analog to CH 4 ·H 2 O, but there are key differences to consider. In particular, the aqueous solubility of xenon is altered by the presence of hydrate, similar to what is observed for methane; but xenon is consistently less soluble than methane. Xenon hydrate has a wider nonstoichiometry region, which could lead to a thicker hydrate layer at the gas‐liquid interface when grown under similar kinetic forcing conditions. For both systems, our numerical calculations reveal that hydrate nonstoichiometry coupled with hydrate formation dynamics leads to a compositional gradient across the hydrate layer, where the stoichiometric ratio increases from the gas‐facing side to the liquid‐facing side. Our analysis suggests that accurate composition measurements could be used to infer the kinetic history of hydrate formation in natural settings where gas is abundant.

Geochemistry, Geophysics, Geosystems

Authigenic mineral texture in submarine 1979 basalt drill core, Surtsey volcano, Iceland

Micrometer-scale maps of authigenic microstructures in submarine basaltic tuff specimens from a 1979 Surtsey volcano, Iceland, drill core acquired 15 years after eruptions terminated provide fresh perspectives for deciphering the initial alteration of oceanic basalt in a low temperature hydrothermal system. A novel investigative approach integrates synchrotron source X-ray microdiffraction (µXRD), X-ray microfluoresence (µXRF), micro-computed tomography (µCT), and scanning transmission electron microscopy (S/TEM) coupled with Raman spectroscopy to create finely resolved spatial frameworks that record a continuum of alteration in glass and olivine. Micro-analytical maps of vesicular and fractured lapilli in specimens from 157.1, 137.9, and 102.6 m depth, and borehole temperatures of 83, 93.9 and 141.3 °C measured in 1980, respectively, record the production of nanocrystalline clay mineral, zeolites, and Al-tobermorite in diverse microenvironments. Nanocrystalline clay mineral (nontronite) and zeolite (amicite) texture in linear microstructures have concentrically-oriented crystallographic preferred orientation. Raman spectra indicating degraded organic carbonaceous matter are associated with nanocrystalline clay mineral in 10–25 nm, sub-circular nanoscale cavities in altered glass at 137.9 m depth and in a concentrically-layered, crystallographically-oriented linear microstructure in altered olivine at 102.6 m. These features have little resemblance to previously described alteration features in basalt. Irregular alteration fronts between fresh and altered glass at 157.1 depth, however, show a resemblance to microchannels in older basalts. The integrated analyses describe the complex organization of previously unrecognized mineral textures in very young basalt and provide a foundational mineralogical reference for longitudinal, time-lapse characterizations of palagonitized basalt in oceanic environments.

Geochemistry, Geophysics, Geosystems

Hydroacoustic, seismic, and bathymetric observations of the 2014 submarine eruption at Ahyi Seamount, Mariana Arc

Ahyi seamount, a shallow submarine volcano in the Northern Mariana Islands, began erupting on April 23, 2014. Hydroacoustic eruption signals were observed on the regional Mariana seismic network and on distant hydrophones, and NOAA scuba divers working in the area soon after the eruption began heard and felt underwater explosion sounds. The NOAA crew observed yellow orange bubble mats along the shore of neighboring Farallon de Pájaros island, but no other surface manifestations of the eruption were reported by the crew or observed in satellite data. Here, we detail the eruption chronology and its morphologic impacts through analysis of seismic and hydroacoustic recordings and repeat bathymetric mapping. Throughout the 2-week-long eruption, Ahyi produced several thousand short, impulsive hydroacoustic signals that we interpret as underwater explosions as well as tremor near the beginning and end of the sequence. The initial tremor, which occurred for 2 hours, is interpreted as small phreatomagmatic explosions. This tremor was followed by a 90 min pause before the characteristic impulsive signals began. Occasional tremor (lasting up to a few minutes) during the last 1.5 days of the eruption is interpreted as more sustained eruptive activity. Bathymetric changes show that a new crater, about 150 m deep, formed near the former summit and a large landslide chute formed on the southeastern flank. Comparing to other geophysically-detected submarine eruptions, we find that the signals from the 2014 Ahyi eruption were more similar to those from other shallow or at surface submarine eruptions than those at deep (>500 m) eruptions.

Geochemistry, Geophysics, Geosystems

Coexisting discrete bodies of rhyolite and punctuated volcanism characterize Yellowstone's post‐Lava Creek Tuff caldera evolution

Ion‐microprobe 206 Pb/ 238 U geochronology and trace element geochemistry of the unpolished rims and sectioned interiors of zircons from Yellowstone caldera's oldest post‐caldera lavas provide insight into the magmatic system during the prelude and aftermath of the caldera‐forming Lava Creek supereruption. The post‐caldera lavas compose the Upper Basin Member of the Plateau Rhyolite and fall into two groups based on zircon crystallization age: early lavas with zircon ages between ~750 and 550 ka and late lavas with zircon ages between ~350 and 250 ka. Zircons from the early‐erupted East Biscuit Basin flow yield U‐Pb dates and trace element compositions, which when considered with the Pb isotopic compositions of their coexisting feldspars and pyroxenes, point to an isotopically distinct parental melt present during crystallization of the Lava Creek magma but untapped by the supereruption. Distinct zircon crystallization ages and Pb‐isotope compositions of major minerals between the early and late Upper Basin Member groups suggest contrasting sources in the magma reservoir. As proxies for melt evolution, the zircons indicate that Yellowstone's post‐caldera rhyolites became more evolved between mid‐ to late‐Pleistocene time, during the same interval that melting of hydrothermally altered wall rock and recharge by new silicic magmas changed in their relative roles. The results from this study indicate that discrete and ephemeral bodies of silicic magma, at times within a mush dominated reservoir and including during the prelude to the Lava Creek eruption, have characterized Yellowstone's subvolcanic reservoir.

Wyoming

SKS splitting beneath Mount St. Helens: Constraints on subslab mantle entrainment

Observations of seismic anisotropy can provide direct constraints on the character of mantle flow in subduction zones, critical for our broader understanding of subduction dynamics. Here we present over 750 new SKS splitting measurements in the vicinity of Mount St. Helens in the Cascadia subduction zone using a combination of stations from the iMUSH broadband array and Cascades Volcano Observatory network. This provides the highest density of splitting measurements yet available in Cascadia, acting as a focused “telescope” for seismic anisotropy in the subduction zone. We retrieve spatially consistent splitting parameters (mean fast direction Φ : 74°, mean delay time ∂t : 1.0 s) with the azimuthal occurrence of nulls in agreement with the fast direction of splitting. When averaged across the array, a 90° periodicity in splitting parameters as a function of back azimuth is revealed, which has not been recovered previously with single‐station observations. The periodicity is characterized by a sawtooth pattern in Φ with a clearly defined 45° trend. We present new equations that reproduce this behavior based upon known systematic errors when calculating shear wave splitting from data with realistic seismic noise. The corrected results suggest a single layer of anisotropy with an ENE‐WSW fast axis parallel to the motion of the subducting Juan de Fuca plate; in agreement with predictions for entrained subslab mantle flow. The splitting pattern is consistent with that seen throughout Cascadia, suggesting that entrainment of the underlying asthenosphere with the subducting slab is coherent and widespread.

Washington

Protracted multipulse emplacement of a post-resurgent pluton: The case of Platoro caldera complex (Southern Rocky Mountain volcanic field, Colorado)

Many eroded calderas expose associated postcollapse plutons, but detailed fieldwork‐supported studies have rarely focused on the internal structure that can contribute to understanding of emplacement dynamics. The Alamosa River monzonite pluton is a postcollapse intrusion at the Platoro caldera complex that erupted six large ignimbrites between 30.2 and 28.8 Ma in the Southern Rocky Mountains volcanic field. Magnetic fabrics in this intrusion indicate the pulsed emplacement of a vertically extensive pluton. The magmatic pulses are documented by three concentric domains of magnetic foliations elongated in ~NE‐SW direction, corresponding to structural trends at the Platoro caldera complex and preexisting regional structures. As no evidence for deformation of wall rocks and the adjacent resurgent block has been identified, we interpret the Alamosa River pluton as a postresurgent intrusion. The space‐opening process involved magmatic stoping and small‐scale magma wedging. New SHRIMP‐RG U/Pb zircon dates (28.98 ±0.18, 27.42 ± 0.35, and 27.32 ± 0.38 Ma) suggest a magmatic lifespan of ~1.7 My for the Alamosa River pluton. Our results indicate that postcaldera magmatism includes pulsed and protracted activity from large intracaldera resurgent plutons to smaller postresurgent stocks and sheeted complexes. As demonstrated by the Alamosa River pluton, some intrusions are emplaced shortly after collapse and resurgence, but postcaldera volcano‐plutonic systems may remain active for several million years or more. We also suggest that subvolcanic magma bodies may be assembled incrementally and that the record of early composite magma lenses preserved as magma wedges are later obliterated by convective flowage and crystallization.

Colorado

Using conceptual models to relate multiparameter satellite data to subsurface volcanic processes in Latin America

Satellite data have been extensively used to identify volcanic behavior. However, the physical subsurface processes causing any individual manifestation of activity can be ambiguous. We propose a classification scheme for the cause of unrest that simultaneously considers three multiparameter satellite observations. The scheme is based on characteristics of the volcanic system (open, closed, and eruptive) and unrest mechanisms (intrusion, evolution, and withdrawal) occurring at shallow depths in the volcanic system. We applied these models to satellite observations acquired at 47 of the most active volcanoes in Latin America. Of the volcanoes studied, 44 had a robust enough dataset for classification and were clustered into 4 groups and 10 subgroups with common behavioral characteristics. By identifying that these volcanoes can be clustered into a number of groupings significantly less than the number of volcanoes, we have demonstrated that commonalities in behavior patterns exist among diverse volcanic systems. Identifying volcanoes with similar characteristics underpins the use of past observations at one volcano to forecast activity at another and diverges from typical volcanic groupings, which are focused on geologic parameters (i.e., composition, volcano type, and tectonic setting). Based on satellite data alone, we have identified preeruptive intrusion prior to 15 eruptions at 12 different volcanoes, magma evolution prior to 18 eruptions at 13 volcanoes, and magma withdrawal at 3 eruptions and 3 volcanoes. Improvements to the spatial and temporal resolution are needed to make these relations robust. This classification scheme provides a framework for future automated clustering of volcanoes.

Geochemistry, Geophysics, Geosystems

A novel method for the extraction, purification and characterization of noble gases in produced fluids

Hydrocarbon systems with declining or viscous oil production are often stimulated using enhanced oil recovery (EOR) techniques, such as the injection of water, steam and CO2, in order to increase oil and gas production. As EOR and other methods of enhancing production such as hydraulic fracturing have become more prevalent, environmental concerns about the impact of both new and historical hydrocarbon production on overlying shallow aquifers have increased. Noble gas isotopes are powerful tracers of subsurface fluid provenance and can be used to understand the impact of EOR on hydrocarbon systems and potentially overlying aquifers. In oil systems, produced fluids can consist of a mixture of oil, water and gas. Noble gases are typically measured in the gas phase; however, it is not always possible to collect gases and therefore produced fluids (which are water, oil and gas mixtures) must be analyzed. We outline a new technique to separate and analyze noble gases in multi-phase hydrocarbon-associated fluid samples. An offline double capillary method has been developed to quantitatively isolate noble gases into a transfer vessel, while effectively removing all water, oil, and less volatile hydrocarbons. The gases are then cleaned and analyzed using standard techniques. Air-saturated water reference materials (n=24) were analyzed and results show a method reproducibility of 2.9% for 4He, 3.8% for 20Ne, 4.5% for 36Ar, 5.3% for 84Kr and 5.7% for 132Xe. This new technique was used to measure the noble gas isotopic compositions in six produced fluid samples from the Fruitvale Oil Field, Bakersfield, California.

Geochemistry, Geophysics, Geosystems