Craig H. Jones

Active foundering of a continental arc root beneath the southern Sierra Nevada in California

Seismic data provide images of crust–mantle interactions during ongoing removal of the dense batholithic root beneath the southern Sierra Nevada mountains in California. The removal appears to have initiated between 10 and 3 Myr ago with a Rayleigh–Taylor-type instability, but with a pronounced asymmetric flow into a mantle downwelling (drip) beneath the adjacent Great Valley. A nearly horizontal shear zone accommodated the detachment of the ultramafic root from its granitoid batholith. With continuing flow into the mantle drip, viscous drag at the base of the remaining ∼35-km-thick crust has thickened the crust by ∼7 km in a narrow welt beneath the western flank of the range. Adjacent to the welt and at the top of the drip, a V-shaped cone of crust is being dragged down tens of kilometres into the core of the mantle drip, causing the disappearance of the Moho in the seismic images. Viscous coupling between the crust and mantle is therefore apparently driving present-day surface subsidence.

Origin of High Mountains in the Continents: The Southern Sierra Nevada

Active and passive seismic experiments show that the southern Sierra, despite standing 1.8 to 2.8 kilometers above its surroundings, is underlain by crust of similar seismic thickness, about 30 to 40 kilometers. Thermobarometry of xenolith suites and magnetotelluric profiles indicate that the upper mantle is eclogitic to depths of 60 kilometers beneath the western and central parts of the range, but little subcrustal lithosphere is present beneath the eastern High Sierra and adjacent Basin and Range. These and other data imply the crust of both the High Sierra and Basin and Range thinned by a factor of 2 since 20 million years ago, at odds with purported late Cenozoic regional uplift of some 2 kilometers.

Crustal thickness variations across the Colorado Rocky Mountains from teleseismic receiver functions

Variations in crustal thickness from the Great Plains of Kansas, across the Colorado Rocky Mountains, and into the eastern Colorado Plateau are determined by receiver function analysis of broadband teleseismic P waveforms recorded during the 1992 Rocky Mountain Front Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL) experiment. The receiver functions are calculated using a time domain deconvolution approach and are interpreted in terms of a single crustal layer, with thickness determined by a grid-search comparison of observed receiver functions with synthetics. The average crustal thicknesses determined by these methods are Kansas Great Plains, 43.8±0.4 km; Colorado Great Plains, 49.9±1.2 km; Colorado Rocky Mountains, 50.1±1.3 km; and northeast Colorado Plateau, 43.1±0.9 at latitudes of 38°–40°N. The main variations in crustal thickness that we observe are between the Kansas Great Plains and the Colorado Great Plains and between the Rocky Mountains and the Colorado Plateau. There is not a significant crustal thickness difference between the Colorado Great Plains and the Colorado Rocky Mountains. Together with gravity data and mass balance calculations, these results are incompatible with the hypothesis that the compensation of the Rocky Mountains relative to the Great Plains is accommodated purely by an Airy-type crustal root or any other mechanism that restricts compensation solely to the crust and requires significant support for the excess topography of the Rocky Mountains to come from the mantle. Models with a rigid elastic plate may match receiver function estimates of crustal thickness but underpredict the amplitude of the gravity low over the Rockies. Our favored model includes lateral variations in crustal velocities obtained from refraction studies and crustal thickness variations constrained by the receiver functions. These models indicate that there is a profound transition in mantle density structure near the eastern range front.

Missing roots and mantle drips : regional Pn and teleseismic arrival times in the southern Sierra Nevada and vicinity, California

Previous seismological studies have placed the source of the uplift of the Sierra either within the crust, suggesting a Mesozoic age for the source of the uplift, or in the upper mantle, consistent with late Cenozoic creation of the buoyant material producing the uplift of the range. We deployed 16 temporary seismometers in the high part of the southern Sierra Nevada to augment the permanent Southern California Seismic Network and record arrivals from regional and teleseismic earthquakes. Arrival times of P waves from 54 teleseisms recorded at these stations are advanced by over a second by a high-velocity body in the upper mantle west and northwest of Lake Isabella. Inversion of the arrival times indicates that this “Isabella anomaly” is of limited north-south extent (about 40–60 km), has compressional velocities about 4–5% higher than its surroundings, and probably extends from about 100 to 200 km depth. The limited north-south extent of the “Isabella anomaly” indicates that it is unrelated to the Sierra; we speculate that it is the downgoing part of a small scale convection system similar to that inferred beneath Southern California. This inversion does not clearly reveal either a large crustal root or a substantial low-velocity body in the upper mantle beneath the Sierra. Although the presence of either degrades the fit to the arrival times and requires high-velocity material beneath the low-velocity material of either root, Bouguer gravity anomalies require low-density material under the Sierra. Assuming that arrival times from earthquakes 150–350 km north and south from the southern Sierra come from a common refractor (the one-layer structure), the upper mantle P wave velocity (P_n) beneath the High Sierra is about 7.6–7.65 km/s; if the arrivals from north and south are from different refractors (the two-layer structure), material with a P wave velocity greater than ∼7.2 km/s (the “7 J.x” layer) would lie under a nearly flat interface from more normal crustal velocities and be separated by a north-dipping interface from underlying mantle with velocities about 7.9–8.1 km/s. The P_n velocity beneath the region immediately to the east is significantly greater (7.9–8.0 km/s) than that of material at equal depths under the Sierra. For the one-layer structure, further assuming that mean crustal velocities are uniform along north-south lines, we find little dip on the Moho in the area; using the arrival times from earthquakes to the south, we infer a depth of 33±5 km for the Moho beneath the southern High Sierra. This structure of a thin to normal crust over a low-velocity mantle can be reconciled with earlier observations that were used to infer a thick crust under the Sierra. By considering the Bouguer gravity anomaly, the surface geology, refraction profiles in this region, and our own observations, we suggest that 1/3 to 1/2 of the modern elevation in the range is supported by lateral (east-west) density contrasts in the crust; the remainder is supported by density contrasts in the uppermost mantle or lateral variations in the thickness of the “7.x” layer. Our interpretation is that the southern Sierra overlies mantle lithosphere that has been thinned and warmed in response to regional lithospheric extension in Neogene time. This part of the upper mantle might have provided the melt that migrated to the east and produced volcanics in the southwestern Great Basin; depletion of the upper mantle might have increased the seismic velocity and decreased the density of material about 60–100 km beneath the southern Sierra.

Tectonics of Pliocene removal of lithosphere of the Sierra Nevada, California

Pliocene (ca. 3.5 Ma) removal of dense eclogitic material under the Sierra Nevada has been proposed from variations in the petrology and geochemistry of Neogene volcanic rocks and their entrained xenoliths from the southern Sierra. The replacement of eclogite by buoyant, warm asthenosphere is consistent with present-day seismologic and magnetotelluric observations made in the southern Sierra. A necessary consequence of replacing eclogite with peridotite is that mean surface elevations and gravitational potential energy both increase. An increase in potential energy should increase extensional strain rates in the area. If these forces are insuffi cient to signifi cantly alter Pacifi c‐ North American plate motion, then increased extensional strain rates in the vicinity of the Sierra must be accompanied by changes in the rate and style of deformation elsewhere. Changes in deformation in California and westernmost Nevada agree well with these predictions. Existing geologic evidence indicates that a period of rapid uplift along the Sierran crest of more than ~1 km occurred between 8 and 3 Ma, most likely as a consequence of removal of lower lithosphere. About this same time, extensional deformation was initiated within ~50 km of the eastern side of the Sierra (5‐3 Ma), and regional shortening began to produce the California Coast Ranges (5‐3 Ma). We suggest that these events were induced by the >1.2 ◊ 10 12 N/m increase of gravitational potential energy generated by the Sierran uplift. Evidence for Pliocene uplift, adjoining crustal extension, and shortening in directly opposing parts of the Coast Ranges is found along nearly the entire length of the Sierra Nevada and implies that lithosphere was removed beneath all of the presentday mountain range. The uplifted area lies between two large, upper-mantle, high-Pwave-velocity bodies under the south end of the San Joaquin Valley and the north end of the Sacramento Valley. These high-velocity bodies plausibly represent the present position of material removed from the base of the crust. Lithospheric removal may also be responsible for shifting of the distribution of transform slip from the San Andreas Fault system to the Eastern California shear zone, a prediction that awaits better-defi ned slip histories on both faults. Overall, the late Cenozoic deformational history of the Sierra Nevada and vicinity illustrates that locally derived forces can infl uence deformation kinematics within plate-boundary zones.

Seismic structure of the lithosphere from teleseismic converted arrivals observed at small arrays in the southern Sierra Nevada and vicinity, California

21 well-distributed teleseisms (30° to 100° distance) were recorded by mixed broad-band (BB) and short-period (SP) seismic arrays at Mineral King (MK) and Horseshoe Meadow (HM) in the southern Sierra Nevada and at Darwin Plateau (DP) between the Inyo and Argus ranges. These arrays permit identification and separation of direct P and S arrivals, reflections from topography, scattered energy, and arrivals from different back azimuths (multipath arrivals). P-to-S conversions can be identified from beamed BB or SP seismograms. Least squares time domain processing recovered single-event receiver functions from these beamed seismograms. Converted phases attributable to the Moho are clear and uncomplicated at DP; Ps-P times of 3.9–4.2 s correspond to crustal thicknesses of 31–33 km, assuming a mean P crustal velocity of 6.2 km/s and a Poissons ratio (ν) of 0.255. Ps-P times at HM (eastern Sierra) are 3.9–4.6 s (32–37 km) excepting some unusual seismograms from events at back azimuths of 225°–239°. MK (western Sierra) times for Ps-P show a strong E-W asymmetry: 3.9–4.1 s (31–33 km) from the east, and 4.8–5.3 s (39–42 km) from the west. The MK arrivals from the east are multiple and substantially weaker by comparison with HM and DP. These results confirm the absence of a thick crust under the southern High Sierra inferred from both a refraction experiment coincident with this experiment and some earlier studies. At DP the Moho event follows an intracrustal negative polarity event defining the top of an S wave low-velocity zone. This feature dips west under the Sierran crest at HM but is absent farther west at MK. This feature appears to be a manifestation of extensional strain and thus indicates that the surficially undeformed Sierra overlies a tectonized lower crust. Sub-Moho energy is absent under the Basin and Range (DP) but is conspicuous under the High Sierra at positive arrivals ∼7.3 s (MK) and ∼9 s (HM) after the P. These arrivals might be from the base of a low-velocity, low-density upper mantle body supporting the Sierra.

Structure of the Sierra Nevada from receiver functions and implications for lithospheric foundering

Receiver functions sampling the Sierra Nevada batholith and adjacent regions exhibit significant variations in the structure of the crust and upper mantle. Crustal Vp/Vs values are lower in the core of the batholith and higher in the northern Sierra Nevada, portions of the Basin and Range, and near young volcanic fields in the eastern Sierra Nevada and Owens Valley. P- to S-wave conversions from the Moho vary from high amplitude and shallow (>25% of the direct P-arrival amplitude, 25–35 km depth) along the eastern Sierra Nevada to low amplitude and deep (

Distributed deformation in the lower crust and upper mantle beneath a continental strike-slip fault zone: Marlborough fault system, South Island, New Zealand

Converted phases from teleseisms recorded by a seismic array spanning the northern half of the Marlborough fault system, South Island, New Zealand, show a continuous unbroken Moho underlying a seismically anisotropic lower crust beneath the two north- ernmost faults of the fault system. These observations suggest that distributed deforma- tion, not slip on a narrow vertical fault, accommodates displacement in the lower crust below the 120-480 km of right-lateral slip across the Wairau fault, one splay of the Marl- borough fault system, and the northward continuation of the Alpine fault. Beneath the Wairau fault, the Moho dips 258-308 southeast from a depth of ;26 km northwest of the fault to a depth of ;34 km southeast of the fault. Farther to the southeast, Ps conversions from the Moho continue under the Awatere fault (34 6 10 km of slip) with a constant amplitude and depth of ;34 km. Across the two faults, converted energy from 16-20 km depth varies with back-azimuth in a manner suggesting the presence of anisotropy in the lower crust. These observations imply that one of the tenets of plate tectonics, that faults defining plate boundaries pass through both crust and upper mantle, does not apply to New Zealand, or to continents in general.

Variations across and along a major continental rift: An interdisciplinary study of the Basin and Range Province, western USA

Abstract Geological, geochemical, and geophysical data gathered within the central part of the Basin and Range and adjacent areas of the western USA suggests that considerable heterogeneity characterizes Cenozoic extension in this region. Good exposure and an abundance of pre-rifting markers indicate 250 km of extension of the upper crust over the past 16 m.y. Extension of several hundred percent has occurred in two distinct deformational domains, Death Valley and Lake Mead, separated by a relatively unextended block, the Spring Mountains. The limited topographic differences between extended and unextended regions imply that material with a crustal density has been added to the extended regions. Although igneous activity can provide some of this added material, kinematics of extension within the Death Valley region suggest that lateral flow of the middle and lower crust into the extended areas accounts for much of the needed material. Such flow is consistent with geochemical analysis of intermediate to silicic volcanic rocks in the Death Valley area. These volcanic rocks contain isotopic and geochemical trends similar to Mesozoic plutonic rocks from the western part of the Sierra Nevada, about 150–200 km to the west, thus suggesting that the upper crust has moved by that amount relative to deeper crustal levels. Geochemical analyses of basaltic magmas in the region indicate that two mantle reservoirs are present: an OIB-type asthenosphere, and an old, Precambrian continental lithosphere. The ancient lithospheric mantle is preserved beneath the Central Basin and Range, but to the west and north the basaltic rocks have a signature compatible with an asthenospheric origin. These differences indicate that the degree of thinning and removal of the mantle lithosphere varies considerably across the Central Basin and Range. These differences are compatible with the inference from geological and geophysical arguments that thinning of the mantle lithosphere at the latitude of the Central Basin and Range is localized beneath the Sierra Nevada. Geophysical measurements have shown that the thickness of the crust varies little from a mean of about 30 km over the entire Basin and Range; the crust under the high Sierra Nevada to the west might have about the same thickness. Estimates of the buoyancy of the crust and mantle based on P-wave crustal structures suggest that the most buoyant, and thus probably the warmest, mantle lies under the Sierra Nevada and not under areas of strongly thinned upper crust of the Death Valley and Lake Mead regions to the east. Similar analyses indicate that the extended upper crust of the Northern Basin and Range overlies an upper mantle more buoyant than that of the Southern and Central Basin and Range; this is in accord with geochemical and seismological inferences. Thus, the style of lithospheric extension varies considerably both along and across the strike of the Basin and Range.

Lithospheric gravitational potential energy and past orogenesis: Implications for conditions of initial Basin and Range and Laramide deformation

Gravitational body forces (i.e., buoyancy forces) have come to be seen as critical to the evolution of orogens. Nevertheless, constraining the role of body forces in specific geologic scenarios is made difficult by the substantial number of poorly constrained physical parameters needed to fully relate forces to deformation. By separating the calculation of buoyancy forces from the calculation of the resulting deformation, models based on relatively simple descriptions of the lithosphere can yield geologically useful constraints. Among these are the importance of paleoelevation in driving syn- and postcontractional extension and in localizing contractional strain. Although such phenomena have been considered in more complex models of continental deformation, the simpler analysis presented here clearly establishes first-order limits on lithospheric structures and paleoelevations consistent with buoyancy-driven deformation. In the early Cenozoic Great Basin of the western United States, we show that the low elevations inferred in much of the geologic literature are inconsistent with a body-force origin for observed extensional tectonism. East of the Colorado Plateau, localization of Laramide deformation coincides with pre-Laramide subsidence of the Western Interior seaway. This subsidence prestressed the lithosphere, making the Southern Rocky Mountains the weak link in responding to regional compressional stress.