Geophysical, petrological, and geochemical data provide important clues about the composition of the deep continental crust. On the basis of seismic refraction data, we divide the crust into type sections associated with different tectonic provinces. Each shows a three-layer crust consisting of upper, middle, and lower crust, in which P wave velocities increase progressively with depth. There is large variation in average P wave velocity of the lower crust between different type sections, but in general, lower crustal velocities are high (>6.9 km s−1) and average middle crustal velocities range between 6.3 and 6.7 km s−1. Heat-producing elements decrease with depth in the crust owing to their depletion in felsic rocks caused by granulite facies metamorphism and an increase in the proportion of mafic rocks with depth. Studies of crustal cross sections show that in Archean regions, 50–85% of the heat flowing from the surface of the Earth is generated within the crust. Granulite terrains that experienced isobaric cooling are representative of middle or lower crust and have higher proportions of mafic rocks than do granulite terrains that experienced isothermal decompression. The latter are probably not representative of the deep crust but are merely upper crustal rocks that have been through an orogenic cycle. Granulite xenoliths provide some of the deepest samples of the continental crust and are composed largely of mafic rock types. Ultrasonic velocity measurements for a wide variety of deep crustal rocks provide a link between crustal velocity and lithology. Meta-igneous felsic, intermediate and mafic granulite, and amphibolite facies rocks are distinguishable on the basis of P and S wave velocities, but metamorphosed shales (metapelites) have velocities that overlap the complete velocity range displayed by the meta-igneous lithologies. The high heat production of metapelites, coupled with their generally limited volumetric extent in granulite terrains and xenoliths, suggests they constitute only a small proportion of the lower crust. Using average P wave velocities derived from the crustal type sections, the estimated areal extent of each type of crust, and the average compositions of different types of granulites, we estimate the average lower and middle crust composition. The lower crust is composed of rocks in the granulite facies and is lithologically heterogeneous. Its average composition is mafic, approaching that of a primitive mantle-derived basalt, but it may range to intermediate bulk compositions in some regions. The middle crust is composed of rocks in the amphibolite facies and is intermediate in bulk composition, containing significant K, Th, and U contents. Average continental crust is intermediate in composition and contains a significant proportion of the bulk silicate Earth's incompatible trace element budget (35–55% of Rb, Ba, K, Pb, Th, and U).

Nature and composition of the continental crust: A lower crustal perspective

A data set comprising 110 spreading rates, 78 transform fault azimuths, and 142 earthquake slip vectors has been inverted to yield a new instantaneous plate motion model, designated Relative Motion 2 (RM2). The model represents a considerable improvement over our previous estimate, RM1 [Minster et al., 1974]. The mean averaging interval for the spreading rate data has been reduced to less than 3 m.y. A detailed comparison of RM2 with angular velocity vectors which best fit the data along individual plate boundaries indicates that RM2 performs close to optimally in most regions, with several notable exceptions. The model systematically misfits data along the India-Antarctica and Pacific-India plate boundaries. We hypothesize that these discrepancies are manifestations of internal deformation within the Indian plate; the data are compatible with northwest-southeast compression across the Ninetyeast Ridge at a rate of about 1 cm/yr. RM2 also fails to satisfy the east-west trending transform fault azimuths observed in the French-American Mid-Ocean Undersea Study area, which is shown to be a consequence of closure constraints about the Azores triple junction. Slow movement between North and South America is required by the data set, although the angular velocity vector describing this motion remains poorly constrained. The existence of a Bering plate, postulated in our previous study, is not necessary if we accept the proposal of Engdahl and others that the Aleutian slip vector data are biased by slab effects. Absolute motion models are derived from several kinematical hypotheses and compared with the data from hot spot traces younger than 10 m.y. Although some of the models are inconsistent with the Wilson-Morgan hypothesis, the overall resolving power of the hot spot data is poor, and the directions of absolute motion for the several slower-moving plates are not usefully constrained.

Present-day plate motions

A data set comprising 110 spreading rates, 78 transform fault azimuths and 142 earthquake slip vectors was inverted to yield a new instantaneous plate motion model, designated RM2. The mean averaging interval for the relative motion data was reduced to less than 3 My. A detailed comparison of RM2 with angular velocity vectors which best fit the data along individual plate boundaries indicates that RM2 performs close to optimally in most regions, with several notable exceptions. On the other hand, a previous estimate (RM1) failed to satisfy an extensive set of new data collected in the South Atlantic Ocean. It is shown that RM1 incorrectly predicts the plate kinematics in the South Atlantic because the presently available data are inconsistent with the plate geometry assumed in deriving RM1. It is demonstrated that this inconsistency can be remedied by postulating the existence of internal deformation with the Indian plate, although alternate explanations are possible.

/pdf/present-day-plate-motions-1cfamgs3qp.pdf

Present‐day plate motions

Knowledge of temperature and pressure, however qualitative, has been central to our views of geology since at least the early 19th century. In 1822, for example, Charles Daubeny presented what may be the very first “Geological Thermometer,” comparing temperatures of various geologic processes (Torrens 2006). Daubeny (1835) may even have been the first to measure the temperature of a lava flow, by laying a thermometer on the top of a flow at Vesuvius—albeit several months following the eruption, after intervening rain (his estimate was 390°F). In any case, pressure ( P ) and temperature ( T ) estimation lie at the heart of fundamental questions: How hot is Earth, and at what rate has the planet cooled. Are volcanoes the products of thermally driven mantle plumes? Where are magmas stored, and how are they transported to the surface—and how do storage and transport relate to plate tectonics? Well-calibrated thermometers and barometers are essential tools if we are to fully appreciate the driving forces and inner workings of volcanic systems.

This chapter presents methods to estimate the P-T conditions of volcanic and other igneous processes. The coverage includes a review of existing geothermometers and geobarometers, and a presentation of approximately 30 new models, including a new plagioclase-liquid hygrometer. Our emphasis is on experimentally calibrated “thermobarometers,” based on analytic expressions using P or T as dependent variables. For numerical reasons (touched on below) such expressions will always provide the most accurate means of P-T estimation, and are also most easily employed. Analytical expressions also allow error to be ascertained; in the absence of estimates of error, P-T estimates are nearly meaningless. This chapter is intended to complement the chapters by Anderson et al. (2008), who cover granitic systems, and by Blundy and Cashman (2008) and Hansteen and Klugel (2008), who consider additional methods …

Thermometers and Barometers for Volcanic Systems

To date, more than 7300 in situ stress orientations have been compiled as part of the World Stress Map project. Of these, over 4400 are considered reliable tectonic stress indicators, recording horizontal stress orientations to within <±25°. Remarkably good correlation is observed between stress orientations deduced from in situ stress measurements and geologic observations made in the upper 1–2 km, well bore breakouts extending to 4–5 km depth and earthquake focal mechanisms to depths of ∼20 km. Regionally uniform stress orientations and relative magnitudes permit definition of broad-scale regional stress patterns often extending 20–200 times the approximately 20–25 km thickness of the upper brittle lithosphere. The “first-order” midplate stress fields are believed to be largely the result of compressional forces applied at plate boundaries, primarily ridge push and continental collision. The orientation of the intraplate stress field is thus largely controlled by the geometry of the plate boundaries. There is no evidence of large lateral stress gradients (as evidenced by lateral variations in stress regime) which would be expected across large plates if simple resistive or driving basal drag tractions (parallel or antiparallel to absolute motion) controlled the intraplate stress field. Intraplate areas of active extension are generally associated with regions of high topography: western U.S. Cordillera, high Andes, Tibetan plateau, western Indian Ocean plateau. Buoyancy stresses related to crustal thickening and/or lithospheric thinning in these regions dominate the intraplate compressional stress field due to plate-driving forces. These buoyancy forces are just one of several categories of “second-order” stresses, or local perturbations, that can be identified once the first-order stress patterns are recognized. These second-order stress fields can often be associated with specific geologic or tectonic features, for example, lithospheric flexure, lateral strength contrasts, as well as the lateral density contrasts which give rise to buoyancy forces. These second-order stress patterns typically have wavelengths ranging from 5 to 10+ times the thickness of the brittle upper lithosphere. A two-dimensional analysis of the amount of rotation of regional horizontal stress orientations due to a superimposed local stress constrains the ratio of the magnitude of the horizontal regional stress differences to the local uniaxial stress. For a detectable rotation of 15°, the local horizontal uniaxial stress must be at least twice the magnitude of the regional horizontal stress differences. Examples of local rotations of SHmax orientations include a 75°–85° rotation on the northeastern Canadian continental shelf possibly related to margin-normal extension derived from sediment-loading flexural stresses, a 50°–60° rotation within the East African rift relative to western Africa due to extensional buoyancy forces caused by lithospheric thinning, and an approximately 90° rotation along the northern margin of the Paleozoic Amazonas rift in central Brazil. In this final example, this rotation is hypothesized as being due to deviatoric compression oriented normal to the rift axis resulting from local lithospheric support of a dense mass in the lower crust beneath the rift (“rift pillow”). Estimates of the magnitudes of first-order (plate boundary force-derived) regional stress differences computed from modeling the source of observed local stress rotations magnitudes can be compared with regional stress differences based on the frictional strength of the crust (i.e., “Byerlee's law”) assuming hydrostatic pore pressure. The examples given here are too few to provide a definitive evaluation of the direct applicability of Byerlee's law to the upper brittle part of the lithosphere, particularly in view of uncertainties such as pore pressure and relative magnitude of the intermediate principal stresses. Nonetheless, the observed rotations all indicate that the magnitude of the local horizontal uniaxial stresses must be 1–2.5+ times the magnitude of the regional first-order horizontal stress differences and suggest that careful evaluation of such local rotations may be a powerful technique for constraining the in situ magnitude stress differences in the upper, brittle part of the lithosphere.

/pdf/first-and-second-order-patterns-of-stress-in-the-lithosphere-595hs86ljh.pdf

First‐ and second‐order patterns of stress in the lithosphere: The World Stress Map Project

Cainozoic evolution of the modern plate boundary along the western United States from subduction to a predominantly transform boundary coincided with a change from compressional to extensional deformation in the western United States. Extensions tectonism responsible for the modern Basin and Range province appears to represent a unique late-stage episode of a much longer period of extension initiated in an ‘in tra -arc ’ setting contemporaneously with calc-alkaline magmatism. Basin-range extension is distinguished from early extension on the basis of angular unconformities, differences in fault trends and spacing, and associated magmatism (basaltic). Prebasin-range extension (i.e. extension preceding the break-up of the region into ranges resembling the modern ones) was under way locally by at least 30 Ma and is now recognized by faulted and highly tilted strata exposed in uplifted range blocks, by large regions of the crust underlain by passively emplaced subvolcanic batholiths, and by the thickness and distribution of stratigraphic units. Locally, high strain rates that accompanied early extensions of as much as 50-100 % are implied. Data on preferentially orientated dyke swarms and fault slip vectors indicate a strikingly uniform WSW -ENE least principal stress orientation in the period ca. 20-10 Ma, during this early extension. The change from early extension to basin-range style faulting of the upper 15 km of crust, which resulted in broadly spaced ranges (25-35 km crest-crest spacing), was time-transgressive and probably not abrupt; locally both types occurred concurrently. Southern Basin and Range block faulting occurred largely in the period 13-10 Ma, in response to a stress field orientated similarly to that responsible for the early extension. In contrast, northern Basin and Range block faulting developed after 10 Ma and continues to the present in response to a stress field orientated approximately 45° clockwise to the earlier stress field. This modern stress field, with a WNW -ESE to E-W directed least principal stress, characterizes the entire modern Basin and Range province and Rio Grande rift region. The 45° change in least principal stress orientation is consistent with superposition of dextral shear associated with the development of the San Andreas transform fault. Inclusion of pre-basin-range extension may help resolve the discrepancy between estimates of 15-30% for basin-range block faulting and total extension estimates of 100-300 % for the Basin and Range province.

Cainozoic Evolution of the State of Stress and Style of Tectonism of the Basin and Range Province of the Western United States

Regional geophysics of the Colorado Plateau

A linear rift composed of dike swarms and graben-filling volcanic rocks marks the inception of Basin and Range rifting in northern Nevada 17 to 14 m.y. ago. This mid-Miocene rift, with its associated aeromagnetic anomaly, indicates S68/sup 0/W-N68/sup 0/E (+-5/sup 0/) extension. A small strike-slip (transform) fault at Sawtooth dike, one short element of the rift, confirms the mid-Miocene extension direction. The present extension, based on data throughout the Basin and Range province, is N65/sup 0/W-S65/sup 0/E (+-20/sup 0/), which is consistent with younger fault offsets of the rift. In the Sawtooth dike region this 45/sup 0/ change in direction of extension took place between 15 and 6 m.y. B.P. The Nevada rift is thought to be a part of a linear rift zone roughly 700 km long, which includes feeder dikes of Columbia River basalts and the mid-Miocene location of the Yellowstone hot spot. Extension in this zone was probably perpendicular to a then-active trench along the western plate margin. The clockwise change in extension direction is consistent with a superposition of right-lateral shear along the plate margin when the San Andreas fault was activated.

Basin and Range rifting in northern Nevada: Clues from a mid-Miocene rift and its subsequent offsets

Synextensional magmatism has long been recognized as a ubiquitous characteristic of highly extended terranes in the western Cordillera of the United States. Intrusive magmatism can have severe effects on the local stress field of the rocks intruded. Because a lower angle fault undergoes increased normal stress from the weight of the upper plate, it becomes more difficult for such a fault to slide. However, if the principal stress orientations are rotated away from vertical and horizontal, then a low-angle fault plane becomes more favored. We suggest that igneous midcrustal inflation occurring at rates faster than regional extension causes increased horizontal stresses in the crust that alter and rotate the principal stresses. Isostatic forces and continued magmatism can work together to create the antiformal or domed detachment surface commonly observed in the metamorphic core complexes of the western Cordillera. Thermal softening caused by magmatism may allow a more mobile mid-crustal isostatic response to normal faulting.

George A. Thompson

Papers

Cainozoic Evolution of the State of Stress and Style of Tectonism of the Basin and Range Province of the Western United States

Regional geophysics of the Colorado Plateau

Basin and Range rifting in northern Nevada: Clues from a mid-Miocene rift and its subsequent offsets

Does magmatism influence low-angle normal faulting?

Implications of selected subsurface data on the structural form and evolution of some basins in the northern Basin and Range province, Nevada and Utah