USGS ScienceSearch

SEARCH · USGS Science

Results for “Journal of Sedimentary Research”

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 199 records · Page 11Linked to original sources

Evidence for a shallow pluton beneath the Goat Rocks Wilderness, Washington, from gravity and magnetic data

A gravity and aeromagnetic study was conducted in and adjacent to the Goat Rocks Wilderness area, Washington. This work was done in conjunction with geologic and geochemical studies which were critically important to our interpretation. The Bouguer gravity anomaly map displays a large central positive gravity anomaly and three superimposed smaller positive anomalies. The main anomaly, about 26 by 11 km with an amplitude of about 16 mGal seems to delineate a dense large pluton or sheeted dike complex intruding less dense, older volcanic and sedimentary rocks lying directly below the volcano. The pluton or sheeted dike complex may have fed the Goat Rocks volcano throughout its 1–2 m.y. history. Three small feeders, mapped on the surface as andesite intrusions, may be the cause of three smaller gravity anomalies superimposed on the main anomaly. The magnetic anomaly maps also exhibit a large central anomaly whose source is the Goat Rocks pluton. The magnetic data suggest that the three shallow intrusive cupolas have undergone a significant amount of hydrothermal alteration. In addition, to some extent the magnetic data reflect the terrain of the area. Topographic highs produce magnetic highs, and topographic lows produce magnetic lows. The deviations from this pattern delineate subsurface geology. Magnetic lows may be of particular interest in locating zones of mineralization where magnetic minerals have been altered and have become significantly less magnetic.

Washington

Origin of two clay-mineral facies of the Potomac Group (Cretaceous) in the Middle Atlantic States

Sedimentary clay that crops out in the coastal plains of New Jersey, Delaware, the part of Maryland north and east of Washington, B.C., and the northeast half of Washington, D.C., in the nonmarine Cretaceous Potomac Group is predominantly kaolinite and illite. In contrast, in part of southeastern Maryland, the southwest half of Washington, D.C., and most of eastern Virginia, Potomac Group clay is predominantly montmorillonite. Kaolinite and illite were probably derived by intense acid weathering of metamorphic and granitic rocks to the west during the Cretaceous and were deposited in a well-drained basin. Montmorillonite was most likely produced by deep weathering of a mainly granitic source. Poor drainage and alkaline surface- and ground-water conditions probably accompanied the production, transportation, deposition, and burial of the montmorillonite. Kaolinization of the montmorillonite is taking place now at the surface in northern Virginia.

Delaware, Maryland, New Jersey, Pennsylvania, Virg

Effects of river regulation on aeolian landscapes, Colorado River, southwestern USA

Connectivity between fluvial and aeolian sedimentary systems plays an important role in the physical and biological environment of dryland regions. This study examines the coupling between fluvial sand deposits and aeolian dune fields in bedrock canyons of the arid to semiarid Colorado River corridor, southwestern USA. By quantifying significant differences between aeolian landscapes with and without modern fluvial sediment sources, this work demonstrates for the first time that the flow- and sediment-limiting effects of dam operations affect sedimentary processes and ecosystems in aeolian landscapes above the fluvial high water line. Dune fields decoupled from fluvial sand supply have more ground cover (biologic crust and vegetation) and less aeolian sand transport than do dune fields that remain coupled to modern fluvial sand supply. The proportion of active aeolian sand area also is substantially lower in a heavily regulated river reach (Marble–Grand Canyon, Arizona) than in a much less regulated reach with otherwise similar environmental conditions (Cataract Canyon, Utah). The interconnections shown here among river flow and sediment, aeolian sand transport, and biologic communities in aeolian dunes demonstrate a newly recognized means by which anthropogenic influence alters dryland environments. Because fluvial–aeolian coupling is common globally, it is likely that similar sediment-transport connectivity and interaction with upland ecosystems are important in other dryland regions to a greater degree than has been recognized previously.

Colorado

Three-dimensional basin and fault structure from a detailed seismic velocity model of Coachella Valley, Southern California

The Coachella Valley in the northern Salton Trough is known to produce destructive earthquakes, making it a high seismic hazard area. Knowledge of the seismic velocity structure and geometry of the sedimentary basins and fault zones is required to improve earthquake hazard estimates in this region. We simultaneously inverted first P wave travel times from the Southern California Seismic Network (39,998 local earthquakes) and explosions (251 land/sea shots) from the 2011 Salton Seismic Imaging Project to obtain a 3-D seismic velocity model. Earthquakes with focal depths ≤10 km were selected to focus on the upper crustal structure. Strong lateral velocity contrasts in the top ~3 km correlate well with the surface geology, including the low-velocity (<5 km/s) sedimentary basin and the high-velocity crystalline basement rocks outside the valley. Sediment thickness is ~4 km in the southeastern valley near the Salton Sea and decreases to <2 km at the northwestern end of the valley. Eastward thickening of sediments toward the San Andreas fault within the valley defines Coachella Valley basin asymmetry. In the Peninsular Ranges, zones of relatively high seismic velocities (~6.4 km/s) between 2 to 4 km depth may be related to Late Cretaceous mylonite rocks or older inherited basement structures. Other high-velocity domains exist in the model down to 9 km depth and help define crustal heterogeneity. We identify a potential fault zone in Lost Horse Valley unassociated with mapped faults in Southern California from the combined interpretation of surface geology, seismicity, and lateral velocity changes in the model.

California

Three Mw ≥ 4.7 earthquakes within the Changning (China) shale gas field ruptured shallow faults intersecting with hydraulic fracturing wells

From 2017 to 2019, three destructive earthquakes (27 January 2017 Mw 4.7, 16 December 2018 Mw 5.2, and 3 January 2019 Mw 4.8) occurred in the Changning shale gas field in the southwest Sichuan Basin, China. Previous seismological studies attributed these events to hydraulic fracturing (HF), but were unable to identify the causative seismogenic faults and their slip behaviors. Here, we use Sentinel-1 synthetic aperture radar data to measure surface deformation triggered by the three events and conduct geodetic inversions to characterize their rupture models. The resulting coseismic interferograms show prominent surface deformation with the maximum line-of-sight displacements of up to 4 cm. The inversion results show that all three earthquakes mainly ruptured sedimentary formations above the shale gas bed, in the upper 3 km of the crust, with slip magnitudes ranging from 8.5 to 15 cm, and stress drops ranging from ∼1.8 to ∼3.3 MPa. Their source faults intersect with horizontal HF wells, but do not root in the crystalline basement. Combined with the reported difficulty of increasing HF operation pressures prior to the three events, we argue that they were most likely induced by direct injection of pressurized fluids into the fault zones. Crustal deformation patterns inferred from regional topography and GPS velocities highlight that the Changning field is located within a triple junction region near the southeastern margin of the Tibetan Plateau with large deformation gradients; such conditions are not only favorable to the development of critically stressed faults, but also facilitate the occurrence of at least moderate magnitude earthquakes.

Changning shale gas field, Sichuan Basin

In situ stress, natural fracture distribution, and borehole elongation in the Auburn Geothermal Well, Auburn, New York

Hydraulic fracturing stress measurements and a borehole televiewer survey were conducted in a 1.6‐km‐deep well at Auburn, New York. This well, which was drilled at the outer margin of the Appalachian Fold and Thrust Belt in the Appalachian Plateau, penetrates approximately 1540 m of lower Paleozoic sedimentary rocks and terminates 60 m into the Precambrian marble basement. Analysis of the hydraulic fracturing tests indicates that the minimum horizontal principal stress increases in a nearly linear fashion from 9.9±0.2 MPa at 593 m to 30.6±0.4 MPa at 1482 m. The magnitude of the maximum horizontal principal stress increases in a less regular fashion from 13.8±1.2 MPa to 49.0±2.0 MPa over the same depth range. The magnitudes of the horizontal principal stresses relative to the calculated overburden stress are somewhat lower than is the norm for this region and are indicative of a strike‐slip faulting regime that, at some depths, is transitional to normal faulting. As expected from the relative aseismicity of central New York State, however, analysis of the magnitudes of the horizontal principal stresses indicates, at least to a depth of 1.5 km, that frictional failure on favorably oriented preexisting fault planes is unlikely. Orientations of the hydraulic fractures at 593 and 919 m indicate that the azimuth of the maximum horizontal principal stress at Auburn is N83°E±15°, in agreement with other stress field indicators for this region. The borehole televiewer log revealed a considerable number of planar features in the Auburn well, the great majority of which are subhorizontal (dips < 5°) and are thought to be bedding plane washouts or drill bit scour marks. In addition, a smaller number of distinct natural fractures were observed on the borehole televiewer log. Of these, the distinct steeply dipping natural fractures in the lower half of the sedimentary section at Auburn tend to strike approximately east‐west, while those in the upper part of the well and in the Precambrian basement exhibit no strong preferred orientation. The origin of this east‐west striking fracture set is uncertain, as it is parallel both to the contemporary direction of maximum horizontal compression and to a late Paleozoic fracture set that has been mapped to the south of Auburn. In addition to these planar features the borehole televiewer log indicates paired dark bands on diametrically opposite sides of the borehole throughout the Auburn well. Processing of the borehole televiewer data in the time domain revealed these features to be irregular depressions in the borehole wall. As these depressions were consistently oriented in a direction at right angles to the direction of maximum horizontal compression, we interpret them to be the result of stress‐induced spalling of the borehole wall (breakouts).

New York

The 12 November 2017 Mw 7.3 Ezgeleh–Sarpolzahab (Iran) earthquake and active tectonics of the Lurestan arc

The 12 November 2017 M w 7.3 Ezgeleh‐Sarpolzahab earthquake is the largest instrumentally recorded earthquake in the Zagros Simply Folded Belt by a factor of ∼10 in seismic moment. Exploiting local, regional, and teleseismic data and synthetic aperture radar interferometry imagery, we characterize the rupture, its aftershock sequence, background seismicity, and regional tectonics. The mainshock ruptured slowly (∼2 km/s), unilaterally southward, for ∼40 km along an oblique (dextral‐thrust) fault that dips ∼14°E beneath the northwestern Lurestan arc. Slip is confined to basement depths of ∼12–18 km, resolvably beneath the sedimentary cover which is ∼8 km thick in this area. The gentle dip angle and basement location allow for a broad slip area, explaining the large magnitude relative to earthquakes in the main Fars arc of the Zagros, where shallower, steeper faults are limited in rupture extent by weak sedimentary layers. Early aftershocks concentrate around the southern and western edges of the mainshock slip area and therefore cluster in the direction of rupture propagation, implying a contribution from dynamic triggering. A cluster of events ∼100 km to the south near Mandali (Iraq) reactivated the ∼50° dipping Zagros Foredeep Fault. The basement fault responsible for the Ezgeleh‐Sarpolzahab earthquake probably accounts for the ∼1 km elevation contrast between the Lurestan arc and the Kirkuk embayment but is distinct from sections of the Mountain Front Fault that define frontal escarpments elsewhere in the Zagros. It may be related to a seismic interface underlying the central and southern Lurestan arc, and a key concern is whether or not the more extensive regional structure is also seismogenic.

Journal of Geophysical Research: Solid Earth

The geologic history of Margaritifer basin, Mars

In this study, we investigate the fluvial, sedimentary, and volcanic history of Margaritifer basin and the Uzboi-Ladon-Morava (ULM) outflow channel system. This network of valleys and basins spans more than 8000 km in length, linking the fluvially dissected southern highlands and Argyre Basin with the northern lowlands via Ares Vallis. Compositionally, thermophysically, and morphologically distinct geologic units are identified and are used to place critical relative stratigraphic constraints on the timing of geologic processes in Margaritifer basin. Our analyses show that fluvial activity was separated in time by significant episodes of geologic activity, including the widespread volcanic resurfacing of Margaritifer basin and the formation of chaos terrain. The most recent fluvial activity within Margaritifer basin appears to terminate at a region of chaos terrain, suggesting possible communication between surface and subsurface water reservoirs. We conclude with a discussion of the implications of these observations on our current knowledge of Martian hydrologic evolution in this important region.

Journal of Geophysical Research E: Planets

Secondary hydroeruptions in pyroclastic-flow deposits: Examples from Mount St. Helens

Secondary hydroeruptions occur in pyroclastic-flow deposits when water or ice is trapped beneath hot pyroclastic debris and rapidly heated to steam. These eruptions display various styles of activity including fumarolic degassing, tephra fountaining, and explosive cratering. The deposits, which occupy the layer 3 stratigraphic position on the top of pyroclastic-flow units, can be distinguished from ash-cloud material by lateral thickness variation, clast composition, and other sedimentary features. The ejecta of secondary hydroeruptions comprise a subset of hydrovolcanic pyroclastic deposits. A small secondary hydroeruption observed on the Mount St. Helens pumice plain in 1981 produced tephra that was emplaced ballistically, by deposition from base surges, and by fallout from an eruption column. Stratigraphic descriptions and grain-size analysis of the ejecta from several secondary craters on the pumice plain demonstrate that the bedforms produced by a hydroeruption change with crater diameter. In particular, craters of small diameter are surrounded by interbedded ripple-laminated ash horizons and nonstratified, fines-depleted units; large craters have ejecta ramparts comprised of coarse dunes and antidunes. These bedform changes are related to a progressive increase in eruptive energy, which produces base surges of greater power and eruptive columns of greater height. We suggest that the style of activity displayed during a secondary hydroeruption is controlled by both the total thermal energy of the system and the permeability of the pyroclastic overburden. ?? 1987.

Journal of Volcanology and Geothermal Research

Brittle deformation and slope failure at the North Menan Butte tuff cone, Eastern Snake River Plain, Idaho

The manifestation of brittle deformation within inactive slumps along the North Menan Butte, a basaltic tuff cone in the Eastern Snake River Plain, is investigated through field and laboratory studies. Microstructural observations indicate that brittle strain is localized along deformation bands, a class of structural discontinuity that is predominant within moderate to high-porosity, clastic sedimentary rocks. Various subtypes of deformation bands are recognized in the study area based on the sense of strain they accommodate. These include dilation bands (no shear displacement), dilational shear bands, compactional shear bands and simple shear bands (no volume change). Measurements of the host rock permeability between the deformation bands indicate that the amount of brittle strain distributed throughout this part of the rock is negligible, and thus deformation bands are the primary means by which brittle strain is manifest within this tuff. Structural discontinuities that are similar in appearance to deformation bands are observed in other basaltic tuffs. Therefore deformation bands may represent a common structural feature of basaltic tuffs that have been widely misclassified as fractures. Slumping and collapse along the flanks of active volcanoes strongly influence their eruptive behavior and structural evolution. Therefore characterizing the process of deformation band and fault growth within basaltic tuff is key to achieving a more complete understanding of the evolution of basaltic volcanoes and their associated hazards.

Idaho

Crustal structure of central Lake Baikal: Insights into intracontinental rifting

The Cenozoic rift system of Baikal, located in the interior of the largest continental mass on Earth, is thought to represent a potential analog of the early stage of breakup of supercontinents. We present a detailed P wave velocity structure of the crust and sediments beneath the Central Basin, the deepest basin in the Baikal rift system. The structure is characterized by a Moho depth of 39–42.5 km; an 8-km-thick, laterally continuous high-velocity (7.05–7.4 km/s) lower crust, normal upper mantle velocity (8 km/s), a sedimentary section reaching maximum depths of 9 km, and a gradual increase of sediment velocity with depth. We interpret the high-velocity lower crust to be part of the Siberian Platform that was not thinned or altered significantly during rifting. In comparison to published results from the Siberian Platform, Moho under the basin is elevated by <3 km. On the basis of these results we propose that the basin was formed by upper crustal extension, possibly reactivating structures in an ancient fold-and-thrust belt. The extent and location of upper mantle extension are not revealed by our data, and it may be offset from the rift. We believe that the Baikal rift structure is similar in many respects to the Mesozoic Atlantic rift system, the precursor to the formation of the North Atlantic Ocean. We also propose that the Central Baikal rift evolved by episodic fault propagation and basin enlargement, rather than by two-stage rift evolution as is commonly assumed.

Lake Baikal

Surface processes recorded by rocks and soils on Meridiani Planum, Mars: Microscopic Imager observations during Opportunity's first three extended missions

The Microscopic Imager (MI) on the Mars Exploration Rover Opportunity has returned images of Mars with higher resolution than any previous camera system, allowing detailed petrographic and sedimentological studies of the rocks and soils at the Meridiani Planum landing site. Designed to simulate a geologist's hand lens, the MI is mounted on Opportunity's instrument arm and can resolve objects 0.1 mm across or larger. This paper provides an overview of MI operations, data calibration, and analysis of MI data returned during the first 900 sols (Mars days) of the Opportunity landed mission. Analyses of Opportunity MI data have helped to resolve major questions about the origin of observed textures and features. These studies support eolian sediment transport, rather than impact surge processes, as the dominant depositional mechanism for Burns formation strata. MI stereo observations of a rock outcrop near the rim of Erebus Crater support the previous interpretation of similar sedimentary structures in Eagle Crater as being formed by surficial flow of liquid water. Well‐sorted spherules dominate ripple surfaces on the Meridiani plains, and the size of spherules between ripples decreases by about 1 mm from north to south along Opportunity's traverse between Endurance and Erebus craters.

Journal of Geophysical Research E: Planets

State factor relationships of dissolved organic carbon and nitrogen losses from unpolluted temperate forest watersheds

We sampled 100 unpolluted, old-growth forested watersheds, divided among 13 separate study areas over 5 years in temperate southern Chile and Argentina, to evaluate relationships among dominant soil-forming state factors and dissolved carbon and nitrogen concentrations in watershed streams. These watersheds provide a unique opportunity to examine broad-scale controls over carbon (C) and nitrogen (N) biogeochemistry in the absence of significant human disturbance from chronic N deposition and land use change. Variations in the ratio dissolved organic carbon (DOC) to nitrogen (DON) in watershed streams differed by underlying soil parent material, with average C:N = 29 for watersheds underlain by volcanic ash and basalt versus C:N = 73 for sedimentary and metamorphic parent materials, consistent with stronger adsorption of low C:N hydrophobic materials by amorphous clays commonly associated with volcanic ash and basalt weathering. Mean annual precipitation was related positively to variations in both DOC (range: 0.2-9.7 mg C/L) and DON (range: 0.008-0.135 mg N/L) across study areas, suggesting that variations in water volume and concentration may act synergistically to influence C and N losses across dry to wet gradients in these forest ecosystems. Dominance of vegetation by broadleaf versus coniferous trees had negligible effects on organic C and N concentrations in comparison to abiotic factors. We conclude that precipitation volume and soil parent material are important controls over chemical losses of dissolved organic C and N from unpolluted temperate forest watersheds. Our results raise the possibility that biotic imprints on watershed C and N losses may be less pronounced in naturally N-poor forests than in areas impacted by land use change and chronic N deposition. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences

PDCD-DAT: A global database of pyroclastic density current deposit field data

Pyroclastic density currents represent one of the deadliest hazards posed by active volcanoes. Analysis of their deposits provides valuable insights into their internal dynamics and informs numerical simulations of pyroclastic density currents which underpin many volcanic hazard assessments. We present PDCD-DAT, a global database of pyroclastic density current deposit characteristics compiled from peer-reviewed literature. The database includes both quantitative datasets (e.g., grain size, density, bedform dimensions, thickness) and qualitative descriptors (e.g., sedimentary structures, lithofacies). PDCD-DAT includes data from 85 source publications, covering 97 eruptions or eruptive phases, and 214 individual depositional units from 55 globally distributed volcanoes. Eruptions recorded in the database range from VEI 1–8. We highlight examples of potential applications of the database, which include (i) comparison of single deposit case studies to global datasets, (ii) informing input parameters and conditions for numerical and analogue models of pyroclastic density currents, (iii) validation of numerical and analogue models against a wide variety of natural case study deposit architectures, and (iv) estimating hazard impact metrics of pyroclastic density currents from past eruptions. We show that the database represents a useful tool for improving our ability to model pyroclastic density currents, predict their associated hazards, and understand the relationships between the internal dynamics of pyroclastic density currents and the properties of their deposits. PDCD-DAT is integrated with the FlowDat Mass Flow Database, which provides a sustainable platform for the database. We aim for PDCD-DAT to be expanded in the future through addition of pyroclastic density current deposit datasets from new field studies conducted by the volcanology research community.

Journal of Applied Volcanology

Implications of the petrochemistry of palladium at Iron Canyon, Lander County, Nevada

Approximately one-half of the 270 samples from the Iron Canyon area, Nevada, analyzed for platinum-group metals, contain measurable amounts of palladium in the 0.001- to 0.02-part per million range with an average of 0.0034 ppm. The rocks include lower Paleozoic sedimentary and volcanic rocks, Tertiary granitic porphyries, and breccia, all of which exhibit various degrees of hydrothermal alteration. Fault-related iron oxides and vein quartz are also present. The area lies astride the outermost fringes of the zone of dispersed alteration visible in outcrops around the middle Tertiary porphyry-copper system at Copper Canyon, Nevada. At Iron Canyon, the palladium concentrations greater than 0.003 ppm seem to be spatially related to the Butte fault zone, a north-striking fault system active during the time of porphyry-type mineralization. Palladium contents vary directly with those of mercury, arsenic, strontium, silver, lanthanum, and boron, as well as with lead, gold, and copper in obviously metallized rocks. This relation probably reflects the mobility of palladium during porphyry-type mineralization. However, our study failed to document conclusively the overall introduction of palladium during mineralization. Among the likely sources of palladium are the lower Paleozoic volcanic rocks or fluids equilibrated with magma (s) associated with the middle Tertiary porphyries.

Nevada

Experimental abrasion of detrital gold

The physical breakdown and abrasion rates of gold were studied using a tumbler to simulate natural high-energy environments. The gold fragments were tumbled for periods ranging from 30 to 240 h with different combinations of sand, cobbles, and water at velocities of 0.5 and 2.0 mi/h (0.85 and 3.22 km/h). With sand and gravel, the common bedload of the rivers that deposited the gold-bearing Tertiary sedimentary rocks of the Sierra Nevada, gold is abraded at rates of 0.015 to 0.007 percent (by weight) per hour of travel (at 0.5 mi/h or 0.845 km/h). Cobbles, rather than sand, are responsible for most of the physical changes and abrasion of the gold. Ten gold fragments tumbled for 120 h with cobbles and water (no sand) were broken down to 68 recoverable fragments and lost about 25 percent of their weight to particles smaller than could be recovered using conventional panning techniques. Gold tumbled for 120 h with sand and water lost less than 1 percent of its weight. Gold was abraded faster by wet sand than by dry sand. Velocity appears to be more important as a factor in abrasion of gold than travel distance a fourfold increase in velocity produced a tenfold increase in hourly abrasion rates of gold. Scanning electron microscope examination of the gold fragments after the tumbling experiments revealed differences in surface texture between fragments tumbled with (1) sand, (2) sand and cobbles, and (3) cobbles only.

Journal of Research of the U.S. Geological Survey

Exploration and discovery in Yellowstone Lake: Results from high-resolution sonar imaging, seismic reflection profiling, and submersible studies

Discoveries from multi-beam sonar mapping and seismic reflection surveys of the northern, central, and West Thumb basins of Yellowstone Lake provide new insight into the extent of post-collapse volcanism and active hydrothermal processes occurring in a large lake environment above a large magma chamber. Yellowstone Lake has an irregular bottom covered with dozens of features directly related to hydrothermal, tectonic, volcanic, and sedimentary processes. Detailed bathymetric, seismic reflection, and magnetic evidence reveals that rhyolitic lava flows underlie much of Yellowstone Lake and exert fundamental control on lake bathymetry and localization of hydrothermal activity. Many previously unknown features have been identified and include over 250 hydrothermal vents, several very large (>500 m diameter) hydrothermal explosion craters, many small hydrothermal vent craters (∼1–200 m diameter), domed lacustrine sediments related to hydrothermal activity, elongate fissures cutting post-glacial sediments, siliceous hydrothermal spire structures, sublacustrine landslide deposits, submerged former shorelines, and a recently active graben. Sampling and observations with a submersible remotely operated vehicle confirm and extend our understanding of the identified features. Faults, fissures, hydrothermally inflated domal structures, hydrothermal explosion craters, and sublacustrine landslides constitute potentially significant geologic hazards. Toxic elements derived from hydrothermal processes also may significantly affect the Yellowstone ecosystem.

Wyoming

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

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.

Journal of Geophysical Research Solid Earth