Search USGSSearch

USGS · 70011892

Lead and strontium isotopes and related trace elements as genetic tracers in the Upper Cenozoic rhyolite-basalt association of the Yellowstone Plateau volcanic field

Abstract

Supported by various field geologic and petrologic data, the contents of Pb, U, Th, Rb, and Sr and the isotopic compositions of Pb and Sr for upper Cenozoic volcanic rocks of the Yellowstone Plateau volcanic field are consistent with the hypothesis of derivation of the basaltic and rhyolitic magmas by partial melting of distinct source regions in the upper mantle and lower crust, respectively. All the basalt samples analyzed but one have systematically lower values of 207 Pb/ 204 Pb and 87 Sr/ 86 Sr than the rhyolites. The values of 206 Pb/ 204 Pb are smaller, and 87 Sr/ 86 Sr are mostly larger than known values in oceanic basalts. In all but one case, the values of 207 Pb/ 204 Pb are higher than expected from an extrapolation of known values in oceanic basalts to less radiogenic values of 206 Pb/ 204 Pb. Because there are no xenoliths, phenocrysts are only moderate to sparse in abundance, REE patterns are low and flat at the radiogenic end of lead isotopic compositions, several values of Rb/Sr are low, and 80% of the basalt samples form a well-developed secondary isochron separate from the rhyolites, we favor an interpretation for basalt genesis wherein isotopic signatures of most mafic magmas were attained in a continental ‘keel’ of mantlelike character about 2.6 b.y. old or somewhat older attached to the crust, and these signatures were unaltered by magma passage through the crust. At the very least, the current data continue to cast serious doubt as to the inevitability of crustal contamination for basaltic magma intruding the continental environment and postulate that much can be learned about the mantle under continents through the study of continental basalts. One basalt unit with an unusually low value of 207 Pb/ 204 Pb and an 87 Ar/ 86 Ar less than 0.704 may represent subcontinental ‘keel’-derived magma that rose unaltered to the surface. Our data also are not consistent with formation of this rhyolite-basalt association primarily by such processes as crystal fractionation, separation of immiscible silicate liquids from a common parental magma, or fractional melting of a homogeneous source. Rather as a conceptual model, we envision large mafic intrusions to have been injected into the lower crust resulting in rhyolite generation through partial anatexis of the adjacent wall rocks which probably had a 206 Pb/ 204 Pb < 17 and 87 Sr/ 86 Sr > 0.709; a model that has much in common with that proposed by Holmes (1931). All the other hypotheses listed have the necessary added complication that either the basalt or the rhyolite or both become contaminated after the two magma types separated, have problems accounting for the lack of igneous rocks of intermediate compositions or production of such large volumes of rhyolitic material (∼5000 km 3 ), and fail to explain why rhyolitic magma is not a more common occurrence in the ocean basin. We appeal to bouyancy of rhyolites to generate a barrier for basalt magma migration and account for the great preponderance of rhyolite relative to basalt at the surface. Furthermore, the complex isotopic picture in the rhyolites indicates that many of these magmas interacted with the upper crustal geologic units that they traversed. The interactions involved diverse processes, probably including reacton with hydrothermal fluids or hydrothermally altered rocks at high levels as well as by contamination with Phanerozoic sedimentary and Precambrian crystalline rocks at deeper levels. At the very least, we feel our study adds a cautionary note to the currently increasingly popular hypothesis that differentiation of basalt or gabbro magmas to rhyolite or granite (as distinct from tonalite or dacite) is a common occurrence and is therefore an important continential building process. Models for formation of rhyolite and granite predominantly by reworking of crust (anatexis) must still be considered. The primitive Archean mantle of the region was characterized by higher Rb/Sr, U/Pb, and Th/U values than are typical of modern suboceanic mantle. The mantle residuum within the continental subcrustal lithosperic ‘keel’ that resulted from the Archean crustal differentiation event probably was depleted in Rb/Sr and U/Pb, and the crust was correspondingly enriched in these ratios. The crust probably was further differentiated by an Archean high-grade metamorphism, during or after the primary event, into a granulitic lower crust depleted in U/Pb and Rb/Sr and a lower-grade upper crust enriched in these ratios.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

B. R. Doe, W.P. Leeman, R.L. Christiansen, C. E. Hedge. 2012-09-20. Lead and strontium isotopes and related trace elements as genetic tracers in the Upper Cenozoic rhyolite-basalt association of the Yellowstone Plateau volcanic field. https://doi.org/10.1029/jb087ib06p04785

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Timescales of surface faulting preservation in low-strain intraplate regions from landscape evolution modeling and the geomorphic and historical record

Large surface-rupturing intraplate earthquakes in stable continental regions (SCRs) are uncommon globally and have recurrence intervals of thousands to hundreds of thousands of years based on the paleoseismic and geomorphic record, challenging accurate active fault identification in these regions. To constrain the timescales of preservation for scarps created by surface ruptures from dip-slip earthquakes, we use a two-dimensional scarp diffusion model for typical intraplate settings and explore which parameters influence fault scarp preservation. These parameters include the coseismic vertical surface offset, the recurrence interval of similar magnitude earthquakes, diffusivity (as a proxy for mean annual precipitation rate), and the erodibility of the surficial material. We constrain parameter ranges from a compilation of historical surface ruptures in intraplate settings in a variety of climates, including the Central and Eastern United States, Australia, Europe, Central Asia (Mongolia, China), India, and West Africa. The timescales of scarp preservation from landscape evolution modeling agree well with observations of scarp preservation in low-strain SCR and intraplate tectonic settings, with some notable exceptions for Australian scarps. We find that the erodibility of the surficial material and earthquake recurrence interval have a stronger effect on the timescales of scarp preservation than diffusivity or coseismic vertical surface offset. Our model results may aid in identifying and characterizing subtle, slow-moving active faults in low-strain SCR and intraplate tectonic settings for different tectonic, geomorphic, and climatic characteristics. Accurate fault locations and characterization from the landscape record has implications for both probabilistic seismic and fault displacement hazard analyses.

Journal of Geophysical Research Solid Earth

Exploring the dynamic interactions between the Southern San Andreas Fault and a normal fault under the Salton Sea

We investigate the dynamic interactions between the Southern San Andreas Fault (SSAF) and a proximal normal fault (NF) beneath the Salton Sea in southern California. The NF, positioned near the SSAF terminus at Bombay Beach, exhibits 11–15 displacement events across 14 stratigraphic sequences, with a range of 0.2–1.4 m of vertical offset since ∼2–3 ka. Notably, four of these events may align temporally with SSAF earthquakes, raising questions about the possible interplay between the two faults. Utilizing dynamic rupture models, we analyze the coseismic interactions between the SSAF and NF, addressing under what conditions the SSAF induces slip on the NF. Our findings reveal that a suite of SSAF ruptures, particularly those propagating from north to south, can trigger slip on the normal fault and replicate observed vertical offsets. If the SSAF extends beneath the Salton Sea, earthquakes originating south of the NF intersection are less likely to trigger normal fault slip, although we cannot exclude this possibility. Some SSAF ruptures do not trigger discernible slip on the NF, rendering such events undetectable in the stratigraphic record. Our research contributes toward discussions regarding the seismic hazard in southern California, shedding light on the interplay between the SSAF and NF.

California

Quantifying site effects and their influence on earthquake source parameter estimations using a dense array in Oklahoma

We investigate the effects of site response on source parameter estimates using earthquakes recorded by the LArge-n Seismic Survey in Oklahoma (LASSO). While it is well known that near-surface unconsolidated sediments can cause an apparent breakdown of earthquake self-similarity, the influence of laterally varying site conditions remains unclear. We analyze site conditions across the 1825-station array on a river plain within an area of 40 km by 23 km using vertical ground motions from 14 regional earthquakes. While the source radiation pattern controls P-wave ground motions below 8 Hz, the surface geology correlates with P-wave ground motions above 8 Hz and S-wave ground motions at 2–21 Hz. Stations installed in alluvial sediments have vertical ground motions that can exceed three times the array median. We use the variation of ground motion of regional earthquakes across the array as a proxy for site effects. The corner frequencies and stress drops of local earthquakes ( M L = 0.01–3) estimated using a standard single-spectra approach show negative correlations with the site-effect proxy, while the seismic moments show positive correlations. In contrast, the spectral-ratio approach effectively shows no correlation. The overall bias is small as expected for this relatively homogeneous structure; accurate estimation of site-related biases requires at least 30 stations. Correcting for site-related biases reduces the standard deviations of the source parameters by less than 13% of the total variations. Remaining variations are partially associated with source directivity and model misfits— as small earthquakes can have complex ruptures.

Oklahoma