USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earth and Planetary Science Letters”

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 289 records · Page 16Linked to original sources
USGS2011

Microbial sulfate reduction and the sulfur budget for a complete section of altered oceanic basalts, IODP Hole 1256D (eastern Pacific)

Sulfide mineralogy and the contents and isotope compositions of sulfur were analyzed in a complete oceanic volcanic section from IODP Hole 1256D in the eastern Pacific, in order to investigate the role of microbes and their effect on the sulfur budget in altered upper oceanic crust. Basalts in the 800m thick volcanic section are affected by a pervasive low-temperature background alteration and have mean sulfur contents of 530ppm, reflecting loss of sulfur relative to fresh glass through degassing during eruption and alteration by seawater. Alteration halos along fractures average 155ppm sulfur and are more oxidized, have high SO 4 /ΣS ratios (0.43), and lost sulfur through oxidation by seawater compared to host rocks. Although sulfur was lost locally, sulfur was subsequently gained through fixation of seawater-derived sulfur in secondary pyrite and marcasite in veins and in concentrations at the boundary between alteration halos and host rocks. Negative δ34S sulfide-S values (down to -30 °) and low temperatures of alteration (down to ~40 °C) point to microbial reduction of seawater sulfate as the process resulting in local additions of sulfide-S. Mass balance calculations indicate that 15-20% of the sulfur in the volcanic section is microbially derived, with the bulk altered volcanic section containing 940ppm S, and with δ34S shifted to -6.0‰) from the mantle value (0 ‰). The bulk volcanic section may have gained or lost sulfur overall. The annual flux of microbial sulfur into oceanic basement based on Hole 1256D is 3-4 X10 10 molSyr-1, within an order of magnitude of the riverine sulfate source and the sedimentary pyrite sink. Results indicate a flux of bacterially derived sulfur that is fixed in upper ocean basement of 7-8 X 10 -8 molcm -2 yr -1 over 15m.y. This is comparable to that in open ocean sediment sites, but is one to two orders of magnitude less than for ocean margin sediments. The global annual subduction of sulfur in altered oceanic basalt lavas based on Hole 1256D is 1.5-2.0 X 10 11 moly -1 , comparable to the subduction of sulfide in sediments, and could contribute to sediment-like sulfur isotope heterogeneities in the mantle.

Earth and Planetary Science Letters↗
USGS2012

On the causes of mid-Pliocene warmth and polar amplification

The mid-Pliocene (~ 3 to 3.3 Ma ago), is a period of sustained global warmth in comparison to the late Quaternary (0 to ~ 1 Ma ago), and has potential to inform predictions of long-term future climate change. However, given that several processes potentially contributed, relatively little is understood about the reasons for the observed warmth, or the associated polar amplification. Here, using a modelling approach and a novel factorisation method, we assess the relative contributions to mid-Pliocene warmth from: elevated CO 2 , lowered orography, and vegetation and ice sheet changes. The results show that on a global scale, the largest contributor to mid-Pliocene warmth is elevated CO 2 . However, in terms of polar amplification, changes to ice sheets contribute significantly in the Southern Hemisphere, and orographic changes contribute significantly in the Northern Hemisphere. We also carry out an energy balance analysis which indicates that that on a global scale, surface albedo and atmospheric emmissivity changes dominate over cloud changes. We investigate the sensitivity of our results to uncertainties in the prescribed CO 2 and orographic changes, to derive uncertainty ranges for the various contributing processes.

Earth and Planetary Science Letters↗
USGS2012

Recycling of water, carbon, and sulfur during subduction of serpentinites: A stable isotope study of Cerro del Almirez, Spain

We use the concentrations and isotope compositions of water, carbon, and sulfur in serpentinites and their dehydration products to trace the cycling of volatiles during subduction. Antigorite serpentinites from the Cerro del Almirez complex, Spain, contain 9–12 wt.% H2O and 910 ± 730 ppm sulfur, and have bulk δ18O values of 8.6 ± 0.4‰, δD = − 54 ± 5‰, and δ34S = 5.0‰, consistent with serpentinization at temperatures of ~ 200 °C by seawater hydrothermal fluids in a seafloor setting. The serpentinites were dehydrated to chlorite–harzburgite (olivine + orthopyroxene + chlorite) at 700 °C and 1.6–1.9 GPa during subduction metamorphism, resulting in loss of water, and sulfur. The chlorite–harzburgites contain 5.7 ± 1.9 wt.% H2O, and have bulk δ18O = 8.0 ± 0.9‰, and δD = − 77 ± 11‰. The rocks contain 650 ± 620 ppm sulfur having δ34S = 1.2‰. Dehydration of serpentinite resulted in loss of 5 wt.% H2O having δ18O = 8–10‰ and δD = − 27 to − 65‰, and loss of 260 ppm sulfur as sulfate, having δ34S = 14.5‰. The contents and δ13C of total carbon in the two rock types overlap, with a broad trend of decreasing carbon contents and δ13C from ~ 1300 to 200 ppm and − 9.6 to − 20.2‰. This reflects mixing between reduced carbon in the rocks (210 ppm, δ13C ≈ − 26‰) and seawater-derived carbonate (δ13C ≈ − 1‰). Our results indicate: 1) Serpentinized oceanic peridotites carry significant amounts of isotopically fractionated water, carbon and sulfur into subduction zones; 2) Subduction of serpentinites to high P and T results in loss of water, and sulfur, which can induce melting and contribute to 18O, D, and 34S enrichments and oxidation of the sub-arc mantle wedge; and 3) Isotopically fractionated water, carbon, and sulfur in serpentinite dehydration products are recycled deeper into the mantle where they can contribute to isotope heterogeneities and may be significant for volatile budgets of the deep Earth.

Cerro Del Almirez;Betic Cordillera In Southern Spa↗
USGS2012

Sources of shaking and flooding during the Tohoku-Oki earthquake: a mixture of rupture styles

Modeling strong ground motions from great subduction zone earthquakes is one of the great challenges of computational seismology. To separate the rupture characteristics from complexities caused by 3D sub-surface geology requires an extraordinary data set such as provided by the recent Mw9.0 Tohoku-Oki earthquake. Here we combine deterministic inversion and dynamically guided forward simulation methods to model over one thousand high-rate GPS and strong motion observations from 0 to 0.25 Hz across the entire Honshu Island. Our results display distinct styles of rupture with a deeper generic interplate event (~Mw8.5) transitioning to a shallow tsunamigenic earthquake (~Mw9.0) at about 25 km depth in a process driven by a strong dynamic weakening mechanism, possibly thermal pressurization. This source model predicts many important features of the broad set of seismic, geodetic and seafloor observations providing a major advance in our understanding of such great natural hazards.

Earth and Planetary Science Letters↗
USGS2012

Morphological and chemical evidence of stromatolitic deposits in the 2.75 Ga Carajás banded iron formation, Brazil

We describe evidence of biogenicity in the morphology and carbon content of well-preserved, Neoarchean samples of banded iron formation (BIF) from Carajás, Brazil. Silica-rich BIF layers contain translucent ellipsoidal or trapezoidal structures (∼5–10 μm diameter) composed of silica, hematite, and kerogen, which are arranged in larger ring-like forms (rosettes). Stable carbon isotope analysis yields a δ 13 C value of −24.5‰ indicating that the contained carbon is likely biogenic. Raman and SEM analyses, as well as wavelength-dispersive X-ray elemental maps, show kerogen inside the rosette forms. Within the iron-rich BIF layers, tubular structures (0.5–5 μm) were observed between hematite granules and blades. Kerogen and kaolinite are present in these structures. Both the rosettes and the tubular structures resemble morphologies that are characteristic of some bacterial species. We hypothesize that the Carajás BIFs originated as biomats formed by one or more species that over time produced large stromatolitic structures. The rosettes and the tubular structures, associated with chert-rich and iron-rich BIF layers, respectively, may represent two different species, or perhaps, two phases of a bacterium life cycle. For example, some modern myxobacteria exhibit similar morphologies in their resting and vegetative stages. Fe(III) precipitation may have occurred by contact of Fe(II) with bacterial slime, leading to oxidation by chemical reactions with exposed polysaccharide hydroxyl and carboxyl groups. The Fe(III) would then have been available for use as a source of energy in a dissimilatory iron reduction type of metabolism. Organic carbon input presumably came from primary producers (not necessarily aerobic) within the local water column, perhaps in shallow-water communities. Alternatively, the carbon may have originated by Fischer–Tropsch synthesis at ocean hydrothermal vents. The observed lateral continuity of BIF layers may perhaps be explained by chemical signaling by the bacteria of favorable or unfavorable environmental conditions, leading to nearly synchronous cell morphogenesis from a vegetative to resting phase and vice versa.

Earth and Planetary Science Letters↗
USGS2012

Crustal seismicity and the earthquake catalog maximum moment magnitudes (Mcmax) in stable continental regions (SCRs): Correlation with the seismic velocity of the lithosphere

A joint analysis of global seismicity and seismic tomography indicates that the seismic potential of continental intraplate regions is correlated with the seismic properties of the lithosphere. Archean and Early Proterozoic cratons with cold, stable continental lithospheric roots have fewer crustal earthquakes and a lower maximum earthquake catalog moment magnitude ( M cmax ). The geographic distribution of thick lithospheric roots is inferred from the global seismic model S40RTS that displays shear-velocity perturbations (δ V S ) relative to the Preliminary Reference Earth Model (PREM). We compare δ V S at a depth of 175 km with the locations and moment magnitudes ( M w ) of intraplate earthquakes in the crust ( Schulte and Mooney, 2005 ). Many intraplate earthquakes concentrate around the pronounced lateral gradients in lithospheric thickness that surround the cratons and few earthquakes occur within cratonic interiors. Globally, 27% of stable continental lithosphere is underlain by δ V S ≥3.0%, yet only 6.5% of crustal earthquakes with M w >4.5 occur above these regions with thick lithosphere. No earthquakes in our catalog with M w >6 have occurred above mantle lithosphere with δ V S >3.5%, although such lithosphere comprises 19% of stable continental regions. Thus, for cratonic interiors with seismically determined thick lithosphere (1) there is a significant decrease in the number of crustal earthquakes, and (2) the maximum moment magnitude found in the earthquake catalog is M cmax =6.0. We attribute these observations to higher lithospheric strength beneath cratonic interiors due to lower temperatures and dehydration in both the lower crust and the highly depleted lithospheric root.

Earth and Planetary Science Letters↗
USGS2012

High-resolution tephrochronology of the Wilson Creek Formation (Mono Lake, California) and Laschamp event using 238 U- 230 Th SIMS dating of accessory mineral rims

Sediments of the Wilson Creek Formation surrounding Mono Lake preserve a high-resolution archive of glacial and pluvial responses along the eastern Sierra Nevada due to late Pleistocene climate change. An absolute chronology for the Wilson Creek stratigraphy is critical for correlating the paleoclimate record to other archives in the western U.S. and the North Atlantic region. However, multiple attempts to date the Wilson Creek stratigraphy using carbonates and tephras yield discordant results due to open-system effects and radiocarbon reservoir uncertainties as well as abundant xenocrysts. New ion microprobe 238 U- 230 Th dating of the final increments of crystallization recorded by allanite and zircon autocrysts from juvenile pyroclasts yield ages that effectively date eruption of key tephra beds and delimit the timing of basal Wilson Creek sedimentation to the interval between 26.8±2.1 and 61.7±1.9 ka. Tephra (Ash 15) erupted during the geomagnetic excursion originally designated the Mono Lake excursion yields an age of 40.8±1.9 ka, indicating that the event is instead the Laschamp excursion. The new ages support a depositional chronology from magnetostratigraphy that indicates quasi-synchronous glacial and hydrologic responses in the Sierra Nevada and Mono Basin to regional climate change, with intervals of lake filling and glacial-snowpack melting that are in phase with peaks in spring insolation.

California↗
USGS2012

Re–Os geochronology of the lacustrine Green River Formation: Insights into direct depositional dating of lacustrine successions, Re–Os systematics and paleocontinental weathering

Lacustrine sedimentary successions provide exceptionally high-resolution records of continental geological processes, responding to tectonic, climatic and magmatic influences. These successions are therefore essential for correlating geological and climatic phenomena across continents and furthermore the globe. Producing accurate geochronological frameworks within lacustrine strata is challenging because the stratigraphy is often bereft of biostratigraphy and directly dateable tuff horizons. The rhenium–osmium (Re–Os) geochronometer is a well-established tool for determining precise and accurate depositional ages of marine organic-rich rocks. Lake systems with stratified water columns are predisposed to the preservation of organic-rich rocks and thus should permit direct Re–Os geochronology of lacustrine strata. We present Re–Os systematics from one of the world's best documented lacustrine systems, the Eocene Green River Formation, providing accurate Re–Os depositional dates that are supported by Ar–Ar and U–Pb ages of intercalated tuff horizons. Precision of the Green River Formation Re–Os dates is controlled by the variation in initial 187 Os/ 188 Os and the range of 187 Re/ 188 Os ratios, as also documented in marine systems. Controls on uptake and fractionation of Re and Os are considered to relate mainly to depositional setting and the type of organic matter deposited, with the need to further understand the chelating precursors of Re and Os in organic matter highlighted. In addition to geochronology, the Re–Os data records the 187Os/188Os composition of lake water (1.41–1.54) at the time of deposition, giving an insight into continental runoff derived from weathering of the geological hinterland of the Green River Formation. Such insights enable us to evaluate fluctuations in continental climatic, tectonic and magmatic processes and provide the ability for chemostratigraphic correlation combined with direct depositional dates. Furthermore, initial 187Os/188Os values can be used as a diagnostic tool to distinguish between lacustrine and marine depositional settings when compared to known oceanic 187 Os/ 188 Os values.

Colorado;Utah;Wyoming↗
USGS2013

Ion-probe U–Pb dating of authigenic and detrital opal from Neogene-Quaternary alluvium

Knowing depositional ages of alluvial fans is essential for many tectonic, paleoclimatic, and geomorphic studies in arid environments. The use of U–Pb dating on secondary silica to establish the age of Neogene-Quaternary clastic sediments was tested on samples of authigenic and detrital opal and chalcedony from depths of ∼25 to 53 m in boreholes at Midway Valley, Nevada. Dating of authigenic opal present as rinds on rock clasts and in calcite/silica cements establishes minimum ages of alluvium deposition; dating of detrital opal or chalcedony derived from the source volcanic rocks gives the maximum age of sediment deposition. Materials analyzed included 12 samples of authigenic opal, one sample of fracture-coating opal from bedrock, one sample of detrital opal, and two samples of detrital chalcedony. Uranium–lead isotope data were obtained by both thermal ionization mass spectrometry and ion-microprobe. Uranium concentrations ranged from tens to hundreds of μg/g. Relatively large U/Pb allowed calculation of 206 Pb/ 238 U ages that ranged from 1.64±0.36 (2 σ ) to 6.16±0.50 Ma for authigenic opal and from 8.34±0.28 to 11.2±1.3 Ma for detrital opal/chalcedony. Three samples with the most radiogenic Pb isotope compositions also allowed calculation of 207 Pb/ 235 U ages, which were concordant with 206 Pb/ 238 U ages from the same samples. These results indicate that basin development at Midway Valley was initiated between about 8 and 6 Ma, and that the basin was filled at long-term average deposition rates of less than 1 cm/ka. Because alluvium in Midway Valley was derived from adjacent highlands at Yucca Mountain, the low rates of deposition determined in this study may imply a slow rate of erosion of Yucca Mountain. Volcanic strata underlying the basin are offset by a number of buried faults to a greater degree than the relatively smooth-sloping bedrock/alluvium contact. These geologic relations indicate that movement on most faults ceased prior to erosional planation and burial. Therefore, ages of the authigenic opal from basal alluvium indicate that the last movement on buried faults was older than about 6 Ma.

Earth and Planetary Science Letters↗
USGS2013

Zinc isotope fractionation during magmatic differentiation and the isotopic composition of the bulk Earth

he zinc stable isotope system has been successfully applied to many and varied fields in geochemistry, but to date it is still not completely clear how this isotope system is affected by igneous processes. In order to evaluate the potential application of Zn isotopes as a proxy for planetary differentiation and volatile history, it is important to constrain the magnitude of Zn isotopic fractionation induced by magmatic differentiation. In this study we present high-precision Zn isotope analyses of two sets of chemically diverse, cogenetic samples from Kilauea Iki lava lake, Hawaii, and Hekla volcano, Iceland, which both show clear evidence of having undergone variable and significant degrees of magmatic differentiation. The Kilauea Iki samples display small but resolvable variations in Zn isotope composition (0.26‰<δ 66 Zn<0.36‰; δ 66 Zn defined as the per mille deviation of a sample's 66 Zn/ 64 Zn compositional ratio from the JMC-Lyon standard), with the most differentiated lithologies exhibiting more positive δ 66 Zn values. This fractionation is likely a result of the crystallization of olivine and/or Fe–Ti oxides, which can both host Zn in their crystal structures. Samples from Hekla have a similar range of isotopic variation (0.22‰<δ 66 Zn<0.33‰), however, the degree of fractionation caused by magmatic differentiation is less significant (only 0.07‰) and no correlation between isotope composition and degree of differentiation is seen. We conclude that high temperature magmatic differentiation can cause Zn isotope fractionation that is resolvable at current levels of precision, but only in compositionally-evolved lithologies. With regards to primitive (ultramafic and basaltic) material, this signifies that the terrestrial mantle is essentially homogeneous with respect to Zn isotopes. Utilizing basaltic and ultramafic sample analyses, from different geologic settings, we estimate that the average Zn isotopic composition of Bulk Silicate Earth is δ 66 Zn=0.28±0.05‰ (2s.d.).

Earth and Planetary Science Letters↗
USGS2014

Seismic structure of the central US crust and upper mantle: Uniqueness of the Reelfoot Rift

Using seismic surface waves recorded with Earthscope's Transportable Array, we apply surface wave imaging to determine 3D seismic velocity in the crust and uppermost mantle. Our images span several Proterozoic and early Cambrian rift zones (Mid-Continent Rift, Rough Creek Graben—Rome trough, Birmingham trough, Southern Oklahoma Aulacogen, and Reelfoot Rift). While ancient rifts are generally associated with low crustal velocity because of the presence of thick sedimentary sequences, the Reelfoot Rift is unique in its association with low mantle seismic velocity. Its mantle low-velocity zone (LVZ) is exceptionally pronounced and extends down to at least 200 km depth. This LVZ is of variable width, being relatively narrow ( &#x223C; 50 km "> ∼50km wide) within the northern Reelfoot Rift, which hosts the New Madrid Seismic Zone (NMSZ). We hypothesize that this mantle volume is weaker than its surroundings and that the Reelfoot Rift consequently has relatively low elastic plate thickness, which would tend to concentrate tectonic stress within this zone. No other intraplate ancient rift zone is known to be associated with such a deep mantle low-velocity anomaly, which suggests that the NMSZ is more susceptible to external stress perturbations than other ancient rift zones.

Reelfoot Rift↗
USGS2013

Continuous gravity measurements reveal a low-density lava lake at Kīlauea Volcano, Hawai‘i

On 5 March 2011, the lava lake within the summit eruptive vent at Kīlauea Volcano, Hawai‘i, began to drain as magma withdrew to feed a dike intrusion and fissure eruption on the volcanoʼs east rift zone. The draining was monitored by a variety of continuous geological and geophysical measurements, including deformation, thermal and visual imagery, and gravity. Over the first ∼14 hours of the draining, the ground near the eruptive vent subsided by about 0.15 m, gravity dropped by more than 100 μGal, and the lava lake retreated by over 120 m. We used GPS data to correct the gravity signal for the effects of subsurface mass loss and vertical deformation in order to isolate the change in gravity due to draining of the lava lake alone. Using a model of the eruptive vent geometry based on visual observations and the lava level over time determined from thermal camera data, we calculated the best-fit lava density to the observed gravity decrease — to our knowledge, the first geophysical determination of the density of a lava lake anywhere in the world. Our result, 950 +/- 300 kg m -3 , suggests a lava density less than that of water and indicates that Kīlaueaʼs lava lake is gas-rich, which can explain why rockfalls that impact the lake trigger small explosions. Knowledge of such a fundamental material property as density is also critical to investigations of lava-lake convection and degassing and can inform calculations of pressure change in the subsurface magma plumbing system.

Hawai'i↗
USGS2013

A new model for the growth of basaltic shields based on deformation of Fernandina volcano, Galápagos Islands

Space-geodetic measurements of surface deformation produced by the most recent eruptions at Fernandina – the most frequently erupting volcano in the Galápagos Archipelago – reveal that all have initiated with the intrusion of subhorizontal sills from a shallow magma reservoir. This includes eruptions from fissures that are oriented both radially and circumferentially with respect to the summit caldera. A Synthetic Aperture Radar (SAR) image acquired 1–2 h before the start of a radial fissure eruption in 2009 captures one of these sills in the midst of its propagation toward the surface. Galápagos eruptive fissures of all orientations have previously been presumed to be fed by vertical dikes, and this assumption has guided models of the origin of the eruptive fissure geometry and overall development of the volcanoes. Our findings allow us to reinterpret the internal structure and evolution of Galápagos volcanoes and of similar basaltic shields. Furthermore, we note that stress changes generated by the emplacement of subhorizontal sills feeding one type of eruption may control the geometry of subsequent eruptive fissures. Specifically, circumferential fissures tend to open within areas uplifted by sill intrusions that initiated previous radial fissure eruptions. This mechanism provides a possible explanation for the pattern of eruptive fissures that characterizes all the western Galápagos volcanoes, as well as the alternation between radial and circumferential fissure eruptions at Fernandina. The same model suggests that the next eruption of Fernandina will be from a circumferential fissure in the area uplifted by the 2009 sill intrusion, just southwest of the caldera rim.

Galápagos Islands↗
USGS2014

Structure and tectonics of the northwestern United States from EarthScope USArray magnetotelluric data

The magnetotelluric component of the EarthScope USArray program has covered over 35% of the continental United States. Resistivity tomography models derived from these data image lithospheric structure and provide constraints on the distribution of fluids and melt within the lithosphere. We present a three-dimensional resistivity model of the northwestern United States which provides new insight into the tectonic assembly of western North America from the Archean to present. Comparison with seismic tomography models reveals regions of correlated and anti-correlated resistivity and velocity that help identify thermal and compositional variations within the lithosphere. Recent (Neogene) tectonic features reflected in the model include the subducting Juan de Fuca–Gorda plate which can be traced beneath the forearc to more than 100 km depth, high lithospheric conductivity along the Snake River Plain, and pronounced lower-crustal and upper-mantle conductivity beneath the Basin and Range. The latter is abruptly terminated to the northwest by the Klamath–Blue Mountains Lineament, which we interpret as an important structure during and since the Mesozoic assembly of the region. This boundary is interpreted to separate hot extended lithosphere from colder, less extended lithosphere. The western edge of Proterozoic North America, as indicated by the Cretaceous initial 87 Sr/ 86 Sr = 0.706 contour, is clearly reflected in the resistivity model. We further image an Archean crustal block (“Pend Oreille block”) straddling the Washington/Idaho border, which we speculate separated from the Archean Medicine Hat block in the Proterozoic. Finally, in the modern Cascades forearc, the geometry and internal structure of the Eocene Siletz terrane is reflected in the resistivity model. The apparent eastern edge of the Siletz terrane under the Cascades arc suggests that pre-Tertiary rocks fill the Washington and Oregon back-arc.

Earth and Planetary Science Letters↗
USGS2013

Crustal-scale recycling in caldera complexes and rift zones along the Yellowstone hotspot track: O and Hf isotopic evidence in diverse zircons from voluminous rhyolites of the Picabo volcanic field, Idaho

Rhyolites of the Picabo volcanic field (10.4–6.6 Ma) in eastern Idaho are preserved as thick ignimbrites and lavas along the margins of the Snake River Plain (SRP), and within a deep (>3 km) borehole near the central axis of the Yellowstone hotspot track. In this study we present new O and Hf isotope data and U–Pb geochronology for individual zircons, O isotope data for major phenocrysts (quartz, plagioclase, and pyroxene), whole rock Sr and Nd isotope ratios, and whole rock geochemistry for a suite of Picabo rhyolites. We synthesize our new datasets with published Ar–Ar geochronology to establish the eruptive framework of the Picabo volcanic field, and interpret its petrogenetic history in the context of other well-studied caldera complexes in the SRP. Caldera complex evolution at Picabo began with eruption of the 10.44±0.27 Ma (U–Pb) Tuff of Arbon Valley (TAV), a chemically zoned and normal-δ 18 O (δ 18 O magma=7.9‰) unit with high, zoned 87 Sr/ 86 Sr i (0.71488–0.72520), and low-ε Nd (0) (−18) and ε Hf (0) (−28). The TAV and an associated post caldera lava flow possess the lowest ε Nd (0) (−23), indicating ∼40–60% derivation from the Archean upper crust. Normal-δ 18 O rhyolites were followed by a series of lower-δ 18 O eruptions with more typical (lower crustal) Sr–Nd–Hf isotope ratios and whole rock chemistry. The voluminous 8.25±0.26 Ma West Pocatello rhyolite has the lowest δ 18 O value (δ 18 O melt =3.3‰), and we correlate it to a 1,000 m thick intracaldera tuff present in the INEL-1 borehole (with published zircon ages 8.04–8.35 Ma, and similarly low-δ 18 O zircon values). The significant (4–5‰) decrease in magmatic-δ 18 O values in Picabo rhyolites is accompanied by an increase in zircon δ 18 O heterogeneity from ∼1‰ variation in the TAV to >5‰ variation in the late-stage low-δ 18 O rhyolites, a trend similar to what is characteristic of Heise and Yellowstone, and which indicates remelting of variably hydrothermally altered tuffs followed by rapid batch assembly prior to eruption. However, due to the greater abundance of low-δ 18 O rhyolites at Picabo, the eruptive framework may reflect an intertwined history of caldera collapse and coeval Basin and Range rifting and hydrothermal alteration. We speculate that the source rocks with pre-existing low-δ 18 O alteration may be related to: (1) deeply buried and unexposed older deposits of Picabo-age or Twin Falls-age low-δ 18 O volcanics; and/or (2) regionally-abundant late Eocene Challis volcanics, which were hydrothermally altered near the surface prior to or during peak Picabo magmatism. Basin and Range extension, specifically the formation of metamorphic core complexes exposed in the region, could have facilitated the generation of low-δ 18 O magmas by exhuming heated rocks and creating the large water-rock ratios necessary for shallow hydrothermal alteration of tectonically (rift zones) and volcanically (calderas) buried volcanic rocks. These interpretations highlight the major processes by which supereruptive volumes of magma are generated in the SRP, mechanisms applicable to producing rhyolites worldwide that are facilitated by plume driven volcanism and extensional tectonics.

Idaho↗
USGS2013

Does calving matter? Evidence for significant submarine melt

During the summer in the northeast Pacific Ocean, the Alaska Coastal Current sweeps water with temperatures in excess of 12 °C past the mouths of glacierized fjords and bays. The extent to which these warm waters affect the mass balance of Alaskan tidewater glaciers is uncertain. Here we report hydrographic measurements made within Icy Bay, Alaska, and calculate rates of submarine melt at Yahtse Glacier, a tidewater glacier terminating in Icy Bay. We find strongly stratified water properties consistent with estuarine circulation and evidence that warm Gulf of Alaska water reaches the head of 40 km-long Icy Bay, largely unaltered. A 10–20 m layer of cold, fresh, glacially-modified water overlies warm, saline water. The saline water is observed to reach up to 10.4 °C within 1.5 km of the terminus of Yahtse Glacier. By quantifying the heat and salt deficit within the glacially-modified water, we place bounds on the rate of submarine melt. The submarine melt rate is estimated at >9 m d −1 , at least half the rate at which ice flows into the terminus region, and can plausibly account for all of the submarine terminus mass loss. Our measurements suggest that summer and fall subaerial calving is a direct response to thermal undercutting of the terminus, further demonstrating the critical role of the ocean in modulating tidewater glacier dynamics.

Alaska↗
USGS2013

Detectability of thermal signatures associated with active formation of ‘chaos terrain’ on Europa

A recent study by Schmidt et al. (2011) suggests that Thera Macula, one of the &ldquo;chaos regions&rdquo; on Europa, may be actively forming over a large liquid water lens. Such a process could conceivably produce a thermal anomaly detectable by a future Europa orbiter or flyby mission, allowing for a direct verification of this finding. Here, we present a set of models that quantitatively assess the surface and subsurface temperatures associated with an actively resurfacing chaos region using constraints from Thera Macula. The results of this numerical study suggest that the surface temperature over an active chaos region can be as high as &sim;200 K. However, low-resolution Galileo Photo-Polarimeter Radiometer (PPR) observations indicate temperatures below 120 K over Thera Macula. This suggests that Thera Macula is not currently active unless an insulating layer of at least a few centimeters in thickness is present, or activity is confined to small regions, reducing the overall intensity of the thermal signature. Alternatively, Thera may have been cooling for at least 10&ndash;100 yr and still contain a subsurface lake, which can take &sim;300,000 yr to crystallize. According to the present study, a more sensitive instrument capable of detecting anomalies &sim;5 K above ambient could detect activity at Thera Macula even if an insulating layer of &sim;50 cm is present.

Earth and Planetary Science Letters↗
USGS2014

Fine-grained linings of leveed channels facilitate runout of granular flows

Catastrophic dense granular flows, such as occur in rock avalanches, debris flows and pyroclastic flows, move as fully shearing mixtures that have approximately 60 vol.% solids and tend to segregate to form coarse-grained fronts and leveed channels. Levees restrict spreading of unconfined flows and form as coarse particles that become concentrated in the top of the flow are transported to the front and then advect to the sides in the flow head. Channels from which most material has drained away down slope are commonly lined with fine-grained deposit, widely thought to remain from the tail of the waning flow. We show how segregation in experimental dense flows of carborundum or sand (300–425 μm) mixed with spherical fine ballotini (150–250 μm), on rough slopes of 27–29°, produces fine-grained channel linings that are deposited with the levees, into which they grade laterally. Maximum runout distance is attained with mixtures containing 30–40% sand, just sufficient to segregate and form levees that are adequately robust to restrict the spreading attributable to the low-friction fines. Resin impregnation and serial sectioning of deliberately arrested experimental flows shows how fines-lined levees form from the flow head; the flows create their own stable ‘conduit’ entirely from the front, which in a geophysical context can play an important mechanistic role in facilitating runout. The flow self-organization ensures that low-friction fines at the base of the segregated channel flow shear over fine-grained substrate in the channel, thus reducing frictional energy losses. We propose that in pyroclastic flows and debris flows, which have considerable mobility attributable to pore-fluid pressures, such fine-grained flow-contact zones form similarly and not only reduce frictional energy losses but also reduce flow–substrate permeability so as to enhance pore-fluid pressure retention. Thus the granular flow self-organization that produces fine-grained channel linings can be an important factor in facilitating long runout of catastrophic geophysical flows on the low slopes (few degrees) of depositional fans and aprons around mountains and volcanoes.

Earth and Planetary Science Letters↗