SEARCH · Search USGS

Results for “Lithosphere”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

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

USGS2012

A tomographic glimpse of the upper mantle source of magmas of the Jemez lineament, New Mexico

The 800-km-long Jemez lineament is the most active volcanic feature in the southwestern United States. It is the southeastern tectonic boundary of the Colorado Plateau and crosses the Rio Grande rift at the Jemez Mountains. The primary volcanism of the lineament is basaltic and has occurred in the last 4.5 m.y. To infer spatial distributions of partial melt in the upper mantle source zones for the Rio Grande rift and the Jemez lineament, we investigated the lateral variations of P wave velocity in the upper mantle beneath these features. We used teleseismic P wave delays recorded at a 22-station network to perform a damped least squares, three-dimensional inversion for these lateral variations. Our technique employed velocity interpolation within a three-dimensional grid of points, rather than using blocks of constant P wave velocity. This method allows highly realistic computation of seismic ray paths as well as accurate computation of the matrix elements in our system of equations. Determinations of resolution of results were done in two independent ways, both of which gave consistent estimates of resolution. In our best resolved volume the inversion showed no significant concentration of relative low velocity for P waves beneath the Rio Grande rift. However, directly beneath the Jemez lineament there is a ∼100-km-wide, 1–2% low-velocity feature in the depth range of 50–160 km. Because of the association of the low P wave velocity with the Jemez volcanic lineament but not with the Rio Grande rift, because lowered P wave velocity can be associated with increased partial melt, and because the volume of recent volcanism at the lineament greatly exceeds that at the rift, we infer that a large magmatic source zone exists beneath the Jemez lineament but not beneath the Rio Grande rift. This implies that the volcanic potential of the Jemez lineaments continues to greatly exceed that of the Rio Grande rift. The mantle source zones for volcanics of the Jemez lineament are not overridden by, but rather track, the motion of the North American plate; this implies that these sources are within the lithospheric plate, as is clarified in the discussion. The magmatic source zones of the Jemez lineament are modeled as due to clockwise rotation of the Colorado Plateau about a pole in northeastern Colorado. This rotation caused extension of the lithosphere beneath the Jemez lineament, permitting concentration there of partially melted rock in the upper mantle.

Journal of Geophysical Research Solid Earth↗

USGS2012

Late Cenozoic volcanism, subduction, and extension in the Lassen region of California, southern Cascade Range

Some 537 volcanic vents younger than 7 Ma are identified and these are classified into five age intervals and five compositional categories based on SiO2 content. Maps of vents by age and composition illustrate regionally representative volcanic trends. By 2 Ma, the eastern limit of volcanism had contracted westward toward the late Quaternary arc. Late Quaternary volcanism is concentrated around and north of the silicic Lassen volcanic center. The belt of most recent volcanism (25-0 ka) has been active since at least 2 Ma. Most mafic volcanism is calcalkaline basalt and basaltic andesite. Normal faults and linear groups of vents are evidence of widespread crustal extension throughout most of the Lassen region. NNW orientation of maximum horizontal stress is indicated. The Lassen volcanic region is thought to occur above the subducting Gorda North plate but also within a broad zone of distributed extension that occurs in the North American lithosphere east and southeast of the present Cascadia subduction zone. The scarcity of volcanic rocks older than 7 Ma suggests that a more compressive lithospheric stress regime prior to the late Miocene extensional episode may have suppressed volcanism. -from Authors

Journal of Geophysical Research Solid Earth↗

USGS1984

Upwarp of anomalous asthenosphere beneath the Rio Grande rift

Continental rifts are possible analogues of mid-ocean ridges, although major plate tectonic features are less clearly observed1. Current thermal models of mid-ocean ridges2-4 consist of solid lithospheric plates overlying the hotter, less viscous asthenosphere, with plate thickness increasing away from the ridge axis. The lithospheric lower boundary lies at or near the melting point isotherm, so that at greater depths higher temperatures account for lower viscosity, lower seismic velocities and possibly partial melting. Upwarp of this boundary at the ridge axis concentrates heat there, thus lowering densities by expansion and raising the sea floor to the level of thermal isostatic equilibrium. At slow spreading ridges, a major central graben forms owing to the mechanics of magma injection into the crust5. Topography, heat flow, gravity and seismic studies support these models. On the continents, a low-velocity channel has been observed, although it is poorly developed beneath ancient cratons6-9. Plate tectonic models have been applied to continental basins and margins10-12, but further similarities to the oceanic models remain elusive. Topographic uplift is often ascribed to Airy type isostatic compensation caused by crustal thickening, rather than thermal compensation in the asthenosphere. Here we discuss the Rio Grande rift, in southwestern United States. Teleseismic P-wave residuals show that regional uplift is explained by asthenosphere uplift rather than crustal thickening. ?? 1984 Nature Publishing Group.

Nature↗

USGS1992

Formation of harzburgite by pervasive melt/rock reaction in the upper mantle

Many mantle peridotite samples are too rich in SiO2 (in the form of orthopyroxene) and have ratios of light to heavy rare earth elements that are too high to be consistent with an origin as the residuum of partial melting of the primitive mantle. Trace element studies of melt/rock reaction zones in the Trinity peridotite provide evidence for reaction of the mantle lithosphere with ascending melts, which dissolved calcium-pyroxene and precipitated orthopyroxene as magma mass decreased. This process can account for the observed major and trace element compositions of lithospheric mantle samples, and may accordingly be prevalent in the upper mantle.

Nature↗

USGS2005

Crustal rheology of the Himalaya and Southern Tibet inferred from magnetotelluric data

The Cenozoic collision between the Indian and Asian continents formed the Tibetan plateau, beginning about 70 million years ago. Since this time, at least 1,400 km of convergence has been accommodated by a combination of underthrusting of Indian and Asian lithosphere, crustal shortening, horizontal extrusion and lithospheric delamination. Rocks exposed in the Himalaya show evidence of crustal melting and are thought to have been exhumed by rapid erosion and climatically forced crustal flow. Magnetotelluric data can be used to image subsurface electrical resistivity, a parameter sensitive to the presence of interconnected fluids in the host rock matrix, even at low volume fractions. Here we present magnetotelluric data from the Tibetan-Himalayan orogen from 77??E to 92??E, which show that low resistivity, interpreted as a partially molten layer, is present along at least 1,000 km of the southern margin of the Tibetan plateau. The inferred low viscosity of this layer is consistent with the development of climatically forced crustal flow in Southern Tibet. ?? 2005 Nature Publishing Group.

Nature↗

USGS2004

Global land ice measurements from space (GLIMS): remote sensing and GIS investigations of the Earth's cryosphere

Concerns over greenhouse‐gas forcing and global temperatures have initiated research into understanding climate forcing and associated Earth‐system responses. A significant component is the Earth's cryosphere, as glacier‐related, feedback mechanisms govern atmospheric, hydrospheric and lithospheric response. Predicting the human and natural dimensions of climate‐induced environmental change requires global, regional and local information about ice‐mass distribution, volumes, and fluctuations. The Global Land‐Ice Measurements from Space (GLIMS) project is specifically designed to produce and augment baseline information to facilitate glacier‐change studies. This requires addressing numerous issues, including the generation of topographic information, anisotropic‐reflectance correction of satellite imagery, data fusion and spatial analysis, and GIS‐based modeling. Field and satellite investigations indicate that many small glaciers and glaciers in temperate regions are downwasting and retreating, although detailed mapping and assessment are still required to ascertain regional and global patterns of ice‐mass variations. Such remote sensing/GIS studies, coupled with field investigations, are vital for producing baseline information on glacier changes, and improving our understanding of the complex linkages between atmospheric, lithospheric, and glaciological processes.

Geocarto International↗

USGS2023

Episodic evolution of a protracted convergent margin revealed by detrital zircon geochronology in the Greater Caucasus

Convergent margins play a fundamental role in the construction and modification of Earth's lithosphere and are characterized by poorly understood episodic processes that occur during the progression from subduction to terminal collision. On the northern margin of the active Arabia-Eurasia collision zone, the Greater Caucasus Mountains provide an opportunity to study a protracted convergent margin that spanned most of the Phanerozoic and culminated in Cenozoic continental collision. However, the main episodes of lithosphere formation and deformation along this margin remain enigmatic. Here, we use detrital zircon U–Pb geochronology from Paleozoic and Mesozoic (meta)sedimentary rocks in the Greater Caucasus, along with select zircon U–Pb and Hf isotopic data from coeval igneous rocks, to link key magmatic and depositional episodes along the Caucasus convergent margin. Devonian to Early Carboniferous rocks were deposited prior to Late Carboniferous accretion of the Greater Caucasus crystalline core onto the Laurussian margin. Permian to Triassic rocks document a period of northward subduction and forearc deposition south of a continental margin volcanic arc in the Northern Caucasus and Scythian Platform. Jurassic rocks record the opening of the Caucasus Basin as a back-arc rift during southward migration of the arc front into the Lesser Caucasus. Cretaceous rocks have few Jurassic-Cretaceous zircons, indicating a period of relative magmatic quiescence and minimal exhumation within this basin. Late Cenozoic closure of the Caucasus Basin juxtaposed the Lesser Caucasus arc to the south against the crystalline core of the Greater Caucasus to the north and led to the formation of a hypothesized terminal suture. We expect this suture to be within ~20 km of the southern range front of the Greater Caucasus because all analysed rocks to the north exhibit a provenance affinity with the crystalline core of the Greater Caucasus.

Basin Research↗

USGS2004

The global short-period wavefield modelled with a Monte Carlo seismic phonon method

At high frequencies (∼1 Hz), much of the seismic energy arriving at teleseismic distances is not found in the main phases (e.g. P , PP , S , etc.) but is contained in the extended coda that follows these arrivals. This coda results from scattering off small-scale velocity and density perturbations within the crust and mantle and contains valuable information regarding the depth dependence and strength of this heterogeneity as well as the relative importance of intrinsic versus scattering attenuation. Most analyses of seismic coda to date have concentrated on S -wave coda generated from lithospheric scattering for events recorded at local and regional distances. Here, we examine the globally averaged vertical-component, 1-Hz wavefield (>10° range) for earthquakes recorded in the IRIS FARM archive from 1990 to 1999. We apply an envelope-function stacking technique to image the average time–distance behavior of the wavefield for both shallow (≤50 km) and deep (≥500 km) earthquakes. Unlike regional records, our images are dominated by P and P coda owing to the large effect of attenuation on PP and S at high frequencies. Modelling our results is complicated by the need to include a variety of ray paths, the likely contributions of multiple scattering and the possible importance of P -to- S and S -to- P scattering. We adopt a stochastic, particle-based approach in which millions of seismic phonons are randomly sprayed from the source and tracked through the Earth. Each phonon represents an energy packet that travels along the appropriate ray path until it is affected by a discontinuity or a scatterer. Discontinuities are modelled by treating the energy normalized reflection and transmission coefficients as probabilities. Scattering probabilities and scattering angles are computed in a similar fashion, assuming random velocity and density perturbations characterized by an exponential autocorrelation function. Intrinsic attenuation is included by reducing the energy contained in each particle as an appropriate function of traveltime. We find that most scattering occurs in the lithosphere and upper mantle, as previous results have indicated, but that some lower-mantle scattering is likely also required. A model with 3 to 4 per cent rms velocity heterogeneity at 4-km scale length in the upper mantle and 0.5 per cent rms velocity heterogeneity at 8-km scale length in the lower mantle (with intrinsic attenuation of Q α = 450 above 200 km depth and Q α = 2500 below 200 km) provides a reasonable fit to both the shallow- and deep-earthquake observations, although many trade-offs exist between the scale length, depth extent and strength of the heterogeneity.

Geophysical Journal International↗

USGS2004

Apparent stress, fault maturity and seismic hazard for normal-fault earthquakes at subduction zones

The behavior of apparent stress for normal-fault earthquakes at subduction zones is derived by examining the apparent stress (?? a = ??Es/Mo, where E s is radiated energy and Mo is seismic moment) of all globally distributed shallow (depth, ?? < 70 km) earthquakes with normal-fault mechanisms that occurred in or near subduction zones between 1987 and 2001 for which Es and M o are available. Accurately determined hypocentres from the Engdahl-Hilst-Buland (EHB) catalogue establish the fine detail of the Wadati-Benioff zone. In many cases, we can relate trends in apparent stress to specific features within the subduction zone and compare these with trends for interplate-thrust earthquakes in the same subduction zones. There are two depth ranges over which ??a, for normal-fault earthquakes shows maxima. The highest and most anomalous values of ??a are found in the deeper depth range from 35 to 70 km. The high-??a, events (up to 5 MPa) are characteristically intraslab and in proximity to zones of intense deformation such as a sharp slab bend or where opposing slabs collide. High-??a events in the region of the shallower maximum (hypocentres between 10-35 km and ??a > 1 MPa) are also generally intraslab, but occur where the lithosphere has just begun subduction beneath the overriding plate. They usually occur in cold slabs near trenches where the direction of plate motion across the trench is oblique to the trench axis, or where there are local contortions or geometrical complexities of the plate boundary. Lower ??a (< 1 MPa) is associated with events occurring at the outer rise (OR) complex (between the OR and the trench axis), as well as with intracrustal events occurring just landward of the trench. The average apparent stress of intraslab-normal-fault earthquakes is considerably higher than the average apparent stress of interplate-thrust-fault earthquakes. In turn, the average ?? a of strike-slip earthquakes in intraoceanic environments is considerably higher than that of intraslab-normal-fault earthquakes. The variation of average ??a with focal mechanism and tectonic regime suggests that the level of ?? a is related to fault maturity. Lower stress drops are needed to rupture mature faults such as those found at plate interfaces that have been smoothed by large cumulative displacements (from hundreds to thousands of kilometres). In contrast, immature faults, such as those on which intraslab-normal-fault earthquakes generally occur, are found in cold and intact lithosphere in which total fault displacement has been much less (from hundreds of metres to a few kilometres). Also, faults on which high ??a oceanic strike-slip earthquakes occur are predominantly intraplate or at evolving ends of transforms. At subduction zones, earthquakes occurring on immature faults are likely to be more hazardous as they tend to generate higher amounts of radiated energy per unit of moment than earthquakes occurring on mature faults. We have identified earthquake pairs in which an interplate-thrust and an intraslab-normal earthquake occurred remarkably close in space and time. The intraslab-normal member of each pair radiated anomalously high amounts of energy compared to its thrust-fault counterpart. These intraslab earthquakes probably ruptured intact slab mantle and are dramatic examples in which Mc (an energy magnitude) is shown to be a far better estimate of the potential for earthquake damage than Mw. This discovery may help explain why loss of life as a result of intraslab earthquakes was greater in the 20th century in Latin America than the fatalities associated with interplate-thrust events that represented much higher total moment release. ?? 2004 RAS.

Geophysical Journal International↗

USGS2004

A physical model for strain accumulation in the San Francisco Bay Region

Strain accumulation in tectonically active regions is generally a superposition of the effects of background tectonic loading, steady-state dislocation processes, such as creep, and transient deformation. In the San Francisco Bay region (SFBR), the most uncertain of these processes is transient deformation, which arises primarily in association with large earthquakes. As such, it depends upon the history of faulting and the rheology of the crust and mantle, which together determine the pattern of longer term (decade-scale) post-seismic response to earthquakes. We utilize a set of 102 GPS velocity vectors in the SFBR in order to characterize the strain rate field and construct a physical model of its present deformation. We first perform an inversion for the continuous velocity gradient field from the discrete GPS velocity field, from which both tensor strain rate and rotation rate may be extracted. The present strain rate pattern is well described as a nearly uniform shear strain rate oriented approximately N34??W (140 nanostrain yr-1) plus a N56??E uniaxial compression rate averaging 20 nanostrain yr-1 across the shear zone. We fit the velocity and strain rate fields to a model of time-dependent deformation within a 135-kin-wide, arcuate shear zone bounded by strong Pacific Plate and Sierra Nevada block lithosphere to the SW and NE, respectively. Driving forces are purely lateral, consisting of shear zone deformation imposed by the relative motions between the thick Pacific Plate and Sierra Nevada block lithospheres. Assuming a depth-dependent viscoelastic structure within the shear zone, we account for the effects of steady creep on faults and viscoelastic relaxation following the 1906 San Francisco and 1989 Loma Prieta earthquakes, subject to constant velocity boundary conditions on the edges of the shear zone. Fault creep is realized by evaluating dislocations on the creeping portions of faults in the fluid limit of the viscoelastic model. A priori plate-boundary(PB)-parallel motion is set to 38 mm yr -1. A grid search based on fitting the observed strain rate pattern yields a mantle viscosity of 1.2 ?? 1019 Pa s and a PB-perpendicular convergence rate of ???3 mm yr-1. Most of this convergence appears to be uniformly distributed in the Pacific-Sierra Nevada plate boundary zone. ?? 2005 RAS.

Geophysical Journal International↗

USGS2006

Mechanical deformation model of the western United States instantaneous strain-rate field

We present a relationship between the long-term fault slip rates and instantaneous velocities as measured by Global Positioning System (GPS) or other geodetic measurements over a short time span. The main elements are the secularly increasing forces imposed by the bounding Pacific and Juan de Fuca (JdF) plates on the North American plate, viscoelastic relaxation following selected large earthquakes occurring on faults that are locked during their respective interseismic periods, and steady slip along creeping portions of faults in the context of a thin-plate system. In detail, the physical model allows separate treatments of faults with known geometry and slip history, faults with incomplete characterization (i.e. fault geometry but not necessarily slip history is available), creeping faults, and dislocation sources distributed between the faults. We model the western United States strain-rate field, derived from 746 GPS velocity vectors, in order to test the importance of the relaxation from historic events and characterize the tectonic forces imposed by the bounding Pacific and JdF plates. Relaxation following major earthquakes (M ??? 8.0) strongly shapes the present strain-rate field over most of the plate boundary zone. Equally important are lateral shear transmitted across the Pacific-North America plate boundary along ???1000 km of the continental shelf, downdip forces distributed along the Cascadia subduction interface, and distributed slip in the lower lithosphere. Post-earthquake relaxation and tectonic forcing, combined with distributed deep slip, constructively interfere near the western margin of the plate boundary zone, producing locally large strain accumulation along the San Andreas fault (SAF) system. However, they destructively interfere further into the plate interior, resulting in smaller and more variable strain accumulation patterns in the eastern part of the plate boundary zone. Much of the right-lateral strain accumulation along the SAF system is systematically underpredicted by models which account only for relaxation from known large earthquakes. This strongly suggests that in addition to viscoelastic-cycle effects, steady deep slip in the lower lithosphere is needed to explain the observed strain-rate field. ?? 2006 The Authors Journal compilation ?? 2006 RAS.

Geophysical Journal International↗

USGS2003

Mantle fault zone beneath Kilauea Volcano, Hawaii

Relocations and focal mechanism analyses of deep earthquakes (???13 kilometers) at Kilauea volcano demonstrate that seismicity is focused on an active fault zone at 30-kilometer depth, with seaward slip on a low-angle plane, and other smaller, distinct fault zones. The earthquakes we have analyzed predominantly reflect tectonic faulting in the brittle lithosphere rather than magma movement associated with volcanic activity. The tectonic earthquakes may be induced on preexisting faults by stresses of magmatic origin, although background stresses from volcano loading and lithospheric flexure may also contribute.

Science↗

USGS1985

Genetic interpretation of lead-isotopic data from the Columbia River basalt group, Oregon, Washington, and Idaho

Lead-isotopic signatures of basalts from the Columbia River Basalt Group of central Oregon, Washington, and western Idaho, and of the high-alumina olivine basalts from southeastern Oregon plot mostly in the field of Pb-isotopic data from volcanic rocks of the Cascade Range. In contrast to the mid-ocean ridge basalts (MORB), the Pb-isotopic data for several of the formations form small, tight clusters, and the Nd- and Sr-isotopic data show discrete variation between these basalt groups. Pseudoisochrons are not observed. A single model will account for the observed isotopic and trace-element data from samples from most of the Columbia River Basalt Group. This model calls for partial melting of the convecting oceanic-type mantle and mixing with, or contamination by, fluids derived from continental-derived sediments. It requires that crustal components are not only subducted along the trench, but that they are also transported in the low-velocity zone of the mantle at least 400 km behind the active volcanic arc into the back-arc environment represented by the Columbia Plateaus province. Lead isotopes in the Saddle Mountains Basalt, the youngest formation, suggest that these magmas were formed by partial melting of a 2.6-b.y.-old subcontinental lithosphere that extends as far west as south-central Washington near Ice Harbor Dam. The source material for these flows is characterized by high Th/U, Th/Pb, Rb/Sr, and Nd/Sm ratios and light REE enrichment. An additional component, presumably granitic crust, is necessary to explain the Pb-isotopic data of the Umatilla and perhaps the Pomona Members of the Saddle Mountains Basalt. This contamination model is supported by published δ 18 O data and petrological interpretation. Published studies of basalts from the Deccan Traps, the Tasmanian dolerites, the Snake River Group, and the volcanics of the Yellowstone Plateau province also suggest that melting of old subcontinental lithosphere beneath the continental crust may be an important process in basalt genesis on a global scale.

Idaho, Oregon, Washington↗

USGS1991

Neogene rotations and quasicontinuous deformation of the Pacific Northwest continental margin

Paleomagnetically determined rotations about vertical axes of 15 to 12 Ma flows of the Miocene Columbia River Basalt Group of Oregon and Washington decrease smoothly with distance from the plate margin, consistent with a simple physical model for continental deformation that assumes the lithosphere behaves as a thin layer of fluid. The average rate of northward translation of the continental margin since 15 Ma calculated from the rotations, using this model, is about 15 mm/yr, which suggests that much of the tangential motion between the Juan de Fuca and North American plates since middle Miocene time has been taken up by deformation of North America. The fluid-like character of the large-scale deformation implies that the brittle upper crust follows the motions of the deeper parts of the lithosphere.

Geology↗

USGS1994

Closing of the Midcontinent-Rift - a far-field effect on Grenvillian compression

The Midcontinent rift formed in the Laurentian supercontinent between 1109 and 1094 Ma. Soon after rifting, stresses changed from extensional to compressional, and the central graben of the rift was partly inverted by thrusting on original extensional faults. Thrusting culminated at about 1060 Ma but may have begun as early as 1080 Ma. On the southwest-trending arm of the rift, the crust was shortened about 30 km; on the southeast-trending arm, strike-slip motion was dominant. The rift developed adjacent to the tectonically active Grenville province, and its rapid evolution from an extensional to a compressional feature at ca. 1080 Ma was coincident with renewal of northwest-directed thrusting in the Grenvllle, probably caused by continent-continent collision. A zone of weak lithosphere created by rifting became the locus for deformation within the otherwise strong continental lithosphere. Stresses transmitted from the Grenville province utilized this weak zone to close and invert the rift.

Geology↗

USGS1996

Three-dimensional crustal structure of the southern Sierra Nevada from seismic fan profiles and gravity modeling

Traveltime data from the 1993 Southern Sierra Nevada Continental Dynamics seismic refraction experiment reveal low crustal velocities in the southern Sierra Nevada and Basin and Range province of California (6.0 to 6.6 km/s), as well as low upper mantle velocities (7.6 to 7.8 km/s). The crust thickens from southeast to northwest along the axis of the Sierra Nevada from 27 km in the Mojave Desert to 43 km near Fresno, California. A crustal welt is present beneath the Sierra Nevada, but the deepest Moho is found under the western slopes, not beneath the highest topography. A density model directly derived from the crustal velocity model but with constant mantle density satisfies the pronounced negative Bouguer anomaly associated with the Sierra Nevada, but shows large discrepancies of >50 mgal in the Great Valley and in the Basin and Range province. Matching the observed gravity with anomalies in the crust alone is not possible with geologically reasonable densities; we require a contribution from the upper mantle, either by lateral density variations or by a thinning of the lithosphere under the Sierra Nevada and the Basin and Range province. Such a model is consistent with the interpretation that the uplift of the present Sierra Nevada is caused and dynamically supported by asthenospheric upwelling or lithospheric thinning under the Basin and Range province and eastern Sierra Nevada.

California↗

USGS2008

A mantle plume beneath California? The mid-Miocene Lovejoy Flood Basalt, northern California

The Lovejoy basalt represents the largest eruptive unit identified in California, and its age, volume, and chemistry indicate a genetic affinity with the Columbia River Basalt Group and its associated mantle-plume activity. Recent field mapping, geochemical analyses, and radiometric dating suggest that the Lovejoy basalt erupted during the mid-Miocene from a fissure at Thompson Peak, south of Susanville, California. The Lovejoy flowed through a paleovalley across the northern end of the Sierra Nevada to the Sacramento Valley, a distance of 240 km. Approximately 150 km 3 of basalt were erupted over a span of only a few centuries. Our age dates for the Lovejoy basalt cluster are near 15.4 Ma and suggest that it is coeval with the 16.1-15.0 Ma Imnaha and Grande Ronde flows of the Columbia River Basalt Group. Our new mapping and age dating support the interpretation that the Lovejoy basalt erupted in a forearc position relative to the ancestral Cascades arc, in contrast with the Columbia River Basalt Group, which erupted in a backarc position. The arc front shifted trenchward into the Sierran block after 15.4 Ma. However, the Lovejoy basalt appears to be unrelated to volcanism of the predominantly calc-alkaline Cascade arc; instead, the Lovejoy is broadly tholeiitic, with trace-element characteristics similar to the Columbia River Basalt Group. Association of the Lovejoy basalt with mid-Miocene flood basalt volcanism has considerable implications for North American plume dynamics and strengthens the thermal "point source" explanation, as provided by the mantle-plume hypothesis. Alternatives to the plume hypothesis usually call upon lithosphere-scale cracks to control magmatic migrations in the Yellowstone-Columbia River basalt region. However, it is difficult to imagine a lithosphere-scale flaw that crosses Precambrian basement and accreted terranes to reach the Sierra microplate, where the Lovejoy is located. Therefore, we propose that the Lovejoy represents a rapid migration of plume-head material, at ??20 cm/yr to the southwest, a direction not previously recognized. ?? 2008 The Geological Society of America.

Special Paper of the Geological Society of America↗

USGS2017

Linkages and feedbacks in orogenic systems: An introduction

Orogenic processes operate at scales ranging from the lithosphere to grain-scale, and are inexorably linked. For example, in many orogens, fault and shear zone architecture controls distribution of heat advection along faults and also acts as the primary mechanism for redistribution of heat-producing material. This sets up the thermal structure of the orogen, which in turn controls lithospheric rheology, the nature and distribution of deformation and strain localization, and ultimately, through localized mechanical strengthening and weakening, the fundamental shape of the developing orogenic wedge ( Fig. 1 ). Strain localization establishes shear zone and fault geometry, and it is the motion on these structures, in conjunction with climate, that often focuses erosional and exhumational processes. This climatic focusing effect can even drive development of asymmetry at the scale of the entire wedge ( Willett et al., 1993 ).

Book chapter↗