SEARCH · Geology Reports

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 37 records · Page 2Linked to original sources

USGS2012

Equivalent strike-slip earthquake cycles in half-space and lithosphere-asthenosphere earth models

By virtue of the images used in the dislocation solution, the deformation at the free surface produced throughout the earthquake cycle by slippage on a long strike-slip fault in an Earth model consisting of an elastic plate (lithosphere) overlying a viscoelastic half-space (asthenosphere) can be duplicated by prescribed slip on a vertical fault embedded in an elastic half-space. For the case in which each earthquake ruptures the entire lithosphere (thickness H ), the half-space equivalent slip rate is as follows: Depth interval 0- H , slip identical to that in lithosphere-asthenosphere model (i.e., abrupt coseismic slip and no subsequent slip); depth interval (2 n −1) H to (2 n +1)H ( n = 1,2,…), slip rate uniform in space and dependent upon time as F n ( t ) exp (− t /τ a ) where F n is a ( n - 1) degree polynomial in t , τ a is twice the asthenosphere relaxation time (η/μ), and t is measured from the instant after the preceding earthquake. The slip rate averaged over the seismic cycle in each depth interval equals the secular rate of relative plate motion. For reasonable values of τ a , slip rates below 5 H do not vary much from that mean value and can be treated as constant. Thus the surface deformation due to the earthquake cycle in the lithosphere-asthenosphere model can be calculated very simply from the half-space model with time-dependent slip in the two depth intervals H −3 H and 3 H −5 H , and uniform slip at a rate equal to the secular relative plate velocity below depth 5 H . Inversion of 1973–1988 geodetic measurements of deformation across the segment of the San Andreas fault in the Transverse Ranges north of Los Angeles for the half-space equivalent slip distribution suggests no significant slip on the fault above 30 km and a uniform slip rate of 36 mm/yr below 30 km. One equivalent lithosphere-asthenosphere model would have a 30-km thick lithosphere and an asthenosphere relaxation time greater than 33 years, but other models are possible.

Journal of Geophysical Research↗

USGS2002

Strike-slip earthquakes in the oceanic lithosphere: Observations of exceptionally high apparent stress

The radiated energies, E S , and seismic moments, M 0 , for 942 globally distributed earthquakes that occurred between 1987 to 1998 are examined to find the earthquakes with the highest apparent stresses (τ a =μ E S / M 0 , where μ is the modulus of rigidity). The globally averaged τ a for shallow earthquakes in all tectonic environments and seismic regions is 0.3 MPa. However, the subset of 49 earthquakes with the highest apparent stresses (τ a greater than about 5.0 MPa) is dominated almost exclusively by strike-slip earthquakes that occur in oceanic environments. These earthquakes are all located in the depth range 7–29 km in the upper mantle of the young oceanic lithosphere. Many of these events occur near plate-boundary triple junctions where there appear to be high rates of intraplate deformation. Indeed, the small rapidly deforming Gorda Plate accounts for 10 of the 49 high-τ a events. The depth distribution of τ a , which shows peak values somewhat greater than 25 MPa in the depth range 20–25 km, suggests that upper bounds on this parameter are a result of the strength of the oceanic lithosphere. A recently proposed envelope for apparent stress, derived by taking 6 per cent of the strength inferred from laboratory experiments for young (less than 30 Ma) deforming oceanic lithosphere, agrees well with the upper-bound envelope of apparent stresses over the depth range 5–30 km. The corresponding depth-dependent shear strength for young oceanic lithosphere attains a peak value of about 575 MPa at a depth of 21 km and then diminishes rapidly as the depth increases. In addition to their high apparent stresses, which suggest that the strength of the young oceanic lithosphere is highest in the depth range 10–30 km, our set of high-τ a earthquakes show other features that constrain the nature of the forces that cause interplate motion. First, our set of events is divided roughly equally between intraplate and transform faulting with similar depth distributions of τ a for the two types. Secondly, many of the intraplate events have focal mechanisms with the T -axes that are normal to the nearest ridge crest or subduction zone and P -axes that are normal to the proximate transform fault. These observations suggest that forces associated with the reorganization of plate boundaries play an important role in causing high-τ a earthquakes inside oceanic plates. Extant transform boundaries may be misaligned with current plate motion. To accommodate current plate motion, the pre-existing plate boundaries would have to be subjected to large horizontal transform push forces. A notable example of this is the triple junction near which the second large aftershock of the 1992 April Cape Mendocino, California, sequence occurred. Alternatively, subduction zone resistance may be enhanced by the collision of a buoyant lithosphere, a process that also markedly increases the horizontal stress. A notable example of this is the Aleutian Trench near which large events occurred in the Gulf of Alaska in late 1987 and the 1998 March Balleny Sea M = 8.2 earthquake within the Antarctic Plate.

Geophysical Journal International↗

USGS2003

The relationship between the instantaneous velocity field and the rate of moment release in the lithosphere

Instantaneous velocity gradients within the continental lithosphere are often related to the tectonic driving forces. This relationship is direct if the forces are secular, as for the case of loading of a locked section of a subduction interface by the downgoing plate. If the forces are static, as for the case of lateral variations in gravitational potential energy, then velocity gradients can be produced only if the lithosphere has, on average, zero strength. The static force model may be related to the long-term velocity field but not the instantaneous velocity field (typically measured geodetically over a period of several years) because over short time intervals the upper lithosphere behaves elastically. In order to describe both the short- and long-term behaviour of an (elastic) lithosphere-(viscoelastic) asthenosphere system in a self-consistent manner, I construct a deformation model termed the expected interseismic velocity (EIV) model. Assuming that the lithosphere is populated with faults that rupture continually, each with a definite mean recurrence time, and that the Earth is well approximated as a linear elastic-viscoelastic coupled system, I derive a simple relationship between the instantaneous velocity field and the average rate of moment release in the lithosphere. Examples with synthetic fault networks demonstrate that velocity gradients in actively deforming regions may to a large extent be the product of compounded viscoelastic relaxation from past earthquakes on hundreds of faults distributed over large ( ≥106 km2) areas.

Geophysical Journal International↗

USGS2019

Melt equilibration depths as sensors of lithospheric thickness during Eurasia-Arabia collision and the uplift of the Anatolian Plateau

A co-investigation of mantle melting conditions and seismic structure revealed an evolutionary record of mantle dynamics accompanying the transition from subduction to collision along the Africa-Eurasia margin and the >1 km uplift of the Anatolian Plateau. New 40 Ar/ 39 Ar dates of volcanic rocks from the Eastern Taurides (southeast Turkey) considerably expand the known spatial extent of Miocene-aged mafic volcanism following a magmatic lull over much of Anatolia that ended at ca. 20 Ma. Mantle equilibration depths for these chemically diverse basalts are interpreted to indicate that early to middle Miocene lithospheric thickness in the region varied from ∼50 km or less near the Bitlis suture zone to ∼80 km near the Inner Tauride suture zone. This southward-tapering lithospheric base could be a vestige of the former interface between the subducted (and now detached) portion of the Arabian plate and the overriding Eurasian plate, and/or a reflection of mantle weakening associated with greater mantle hydration trenchward prior to collision. Asthenospheric upwelling driven by slab tearing and foundering along this former interface, possibly accompanied by convective removal of the lithosphere, could have led to renewed volcanic activity after 20 Ma. Melt equilibration depths for late Miocene and Pliocene basalts together with seismic imaging of the present lithosphere indicate that relatively invariant lithospheric thicknesses of 60–70 km have persisted since the middle Miocene. Thus, no evidence is found for large-scale (tens of kilometers) Miocene delamination of the lower lithosphere from the overriding plate, which has been proposed elsewhere to account for late Miocene and younger uplift of Anatolia.

Anatolian Plate↗

USGS2017

Magnetotelluric imaging of lower crustal melt and lithospheric hydration in the Rocky Mountain Front transition zone, Colorado, USA

We present an electrical resistivity model of the crust and upper mantle from two‐dimensional (2‐D) anisotropic inversion of magnetotelluric data collected along a 450 km transect of the Rio Grande rift, southern Rocky Mountains, and High Plains in Colorado, USA. Our model provides a window into the modern‐day lithosphere beneath the Rocky Mountain Front to depths in excess of 150 km. Two key features of the 2‐D resistivity model are (1) a broad zone (~200 km wide) of enhanced electrical conductivity (<20 Ωm) in the midcrust to lower crust that is centered beneath the highest elevations of the southern Rocky Mountains and (2) hydrated lithospheric mantle beneath the Great Plains with water content in excess of 100 ppm. We interpret the high conductivity region of the lower crust as a zone of partially molten basalt and associated deep‐crustal fluids that is the result of recent (less than 10 Ma) tectonic activity in the region. The recent supply of volatiles and/or heat to the base of the crust in the late Cenozoic implies that modern‐day tectonic activity in the western United States extends to at least the western margin of the Great Plains. The transition from conductive to resistive upper mantle is caused by a gradient in lithospheric modification, likely including hydration of nominally anhydrous minerals, with maximum hydration occurring beneath the Rocky Mountain Front. This lithospheric “hydration front” has implications for the tectonic evolution of the continental interior and the mechanisms by which water infiltrates the lithosphere.

Colorado↗

USGS1998

Mountains and Calderas on Io: Possible Implications for Lithosphere Structure and Magma Generation

The combination of Voyager images and newly acquired Galileo images with low illumination and resolutions ranging from 2 to 6 km/pixel now allows determination of the global distribution of mountains and volcanic centers on Io. The mountains generally do not have characteristics typical of terrestrial volcanic landforms, they are evenly distributed across the surface and show no obvious correlation with known hot spots or plumes. Relative elevations, determined by shadow measurements and stereoscopy, indicate that mountains in the newly imaged area range in elevation up to at least 7.6 km. The origin of the mountains remains uncertain. Some appear to be multitiered volcanic constructs; others enclosing the partial remains of large circular depressions appear to be remnants of old volcanoes; yet others show extensive tectonic disruption. Volcanic centers also appear to be distributed evenly across the surface except for an apparently somewhat lower density at high latitudes. The low latitudes have one volcanic center per 7 ?? 104km2, and, on average, the centers are spaced roughly 250 km apart. The global distribution of high mountains suggests that the lithosphere over most of Io is thick. Although the thickness cannot be calculated, the previously suggested 30 km appears reasonable as a lower limit. The high rates of resurfacing combined with the likely dissipation of most of the tidal energy in the asthenosphere and underlying mantle implies a very low temperature gradient in the upper part of the lithosphere and steep gradients in the lower lithosphere. The slow rate of separation of melt from host rock in the magma source regions as a consequence of the low gravity on Io, coupled with the high rate of magma production, will likely result in larger melt fractions than is typical for source regions on Earth. The variety of volcanic landforms suggests that volcanic products with a range of compositions are deposited on the surface. This mixture will be carried downward through the lithosphere as a consequence of the 0.5-1.5 cm/yr resurfacing rates. During descent, the more volatile components will tend to be driven off early, but complete or near-complete melting at the base of the lithosphere may result in rehomogenization of the silicate mixture that remains. ?? 1998 Academic Press.

Icarus↗

USGS2012

Thermal stresses due to cooling of a viscoelastic oceanic lithosphere

Theories based upon thermal contraction of cooling oceanic lithosphere provide a successful basis for correlating seafloor bathymetry and heat flow. The horizontal components of the contraction of the lithosphere as it cools potentially give rise to large thermal stresses. Current methods to calculate these stresses assume that on the time scales of cooling, the lithosphere initially behaves as an inviscid fluid and instantly freezes into an elastic solid at some critical temperature. These instant-freezing methods inaccurately predict transient thermal stresses in rapidly cooling silicate glass plates because of the temperature dependent rheology of the material. The temperature dependent rheology of the lithosphere may affect the transient thermal stress distribution in a similar way, and for this reason we use a thermoviscoelastic model to estimate thermal stresses in young oceanic lithosphere. This theory is formulated here for linear creep processes that have an Arrhenius rate dependence on temperature. Our results show that the stress differences between instant freezing and linear thermoviscoelastic theory are most pronounced at early times (0–20 m.y.) when the instant freezing stresses may be twice as large. The solutions for the two methods asymptotically approach the same solution with time. A comparison with intraplate seismicity shows that both methods underestimate the depth of compressional stresses inferred from the seismicity in a systematic way.

Journal of Geophysical Research Solid Earth↗

USGS2003

Lithospheric roots beneath western Laurentia: The geochemical signal in mantle garnets

This study presents major and trace element data for 243 mantle garnet xenocrysts from six kimberlites in parts of western North America. The geochemical data for the garnet xenocrysts are used to infer the composition, thickness, and tectonothermal affinity of the mantle lithosphere beneath western Laurentia at the time of kimberlite eruption. The garnets record temperatures between 800 and 1450??C using Ni-in-garnet thermometry and represent mainly lherzolitic mantle lithosphere sampled over an interval from about 110-260 km depth. Garnets with sinuous rare-earth element patterns, high Sr, and high Sc/V occur mainly at shallow depths and occur almost exclusively in kimberlites interpreted to have sampled Archean mantle lithosphere beneath the Wyoming Province in Laurentia, and are notably absent in garnets from kimberlites erupting through the Proterozoic Yavapai Mazatzal and Trans-Hudson provinces. The similarities in depths of equilibration, but differing geochemical patterns in garnets from the Cross kimberlite (southeastern British Columbia) compared to kimberlites in the Wyoming Province argue for post-Archean replacement and (or) modification of mantle beneath the Archean Hearne Province. Convective removal of mantle lithosphere beneath the Archean Hearne Province in a "tEctonic vise" during the Proterozoic terminal collisions that formed Laurentia either did not occur, or was followed by replacement of thick mantle lithosphere that was sampled by kimberlite in the Triassic, and is still observed there seismically today.

Canadian Journal of Earth Sciences↗

USGS2014

Time causal operational estimation of electric fields induced in the Earth's lithosphere during magnetic storms

In support of projects for monitoring geomagnetic hazards for electric power grids, we develop a simple mathematical formalism, consistent with the time causality of deterministic physics, for estimating electric fields that are induced in the Earth's lithosphere during magnetic storms. For an idealized model of the lithosphere, an infinite half‐space having uniform electrical conductivity properties described by a galvanic tensor, we work in the Laplace‐transformed frequency domain to obtain a transfer function which, when convolved with measured magnetic field time series, gives an estimated electric field time series. Using data collected at the Kakioka, Japan observatory, we optimize lithospheric conductivity parameters by minimizing the discrepancy between model‐estimated electric field variation and that actually measured. With our simple model, we can estimate 87% of the variance in storm time Kakioka electric field data; a more complicated model of lithospheric conductivity would be required to estimate the remaining 13% of the variance. We discuss how our estimation formalism might be implemented for geographically coordinated real‐time monitoring of geoelectric fields.

Geophysical Research Letters↗

USGS2015

A random-walk algorithm for modeling lithospheric density and the role of body forces in the evolution of the Midcontinent Rift

This paper develops a Monte Carlo algorithm for extracting three-dimensional lithospheric density models from geophysical data. Empirical scaling relationships between velocity and density create a 3D starting density model, which is then iteratively refined until it reproduces observed gravity and topography. This approach permits deviations from uniform crustal velocity-density scaling, which provide insight into crustal lithology and prevent spurious mapping of crustal anomalies into the mantle. We test this algorithm on the Proterozoic Midcontinent Rift (MCR), north-central U.S. The MCR provides a challenge because it hosts a gravity high overlying low shear-wave velocity crust in a generally flat region. Our initial density estimates are derived from a seismic velocity/crustal thickness model based on joint inversion of surface-wave dispersion and receiver functions. By adjusting these estimates to reproduce gravity and topography, we generate a lithospheric-scale model that reveals dense middle crust and eclogitized lowermost crust within the rift. Mantle lithospheric density beneath the MCR is not anomalous, consistent with geochemical evidence that lithospheric mantle was not the primary source of rift-related magmas and suggesting that extension occurred in response to far-field stress rather than a hot mantle plume. Similarly, the subsequent inversion of normal faults resulted from changing far-field stress that exploited not only warm, recently faulted crust but also a gravitational potential energy low in the MCR. The success of this density modeling algorithm in the face of such apparently contradictory geophysical properties suggests that it may be applicable to a variety of tectonic and geodynamic problems.

Iowa, Kansas, Michigan, Minnesota, Nebraska, Wisco↗

USGS1985

Rock mechanics observations pertinent to the rheology of the continental lithosphere and the localization of strain along shear zones

Emphasized in this paper are the deformation processes and rheologies of rocks at high temperatures and high effective pressures, conditions that are presumably appropriate to the lower crust and upper mantle in continental collision zones. Much recent progress has been made in understanding the flexure of the oceanic lithosphere using rock-mechanics-based yield criteria for the inelastic deformations at the top and base. At mid-plate depths, stresses are likely to be supported elastically because bending strains and elastic stresses are low. The collisional tectonic regime, however, is far more complex because very large permanent strains are sustained at mid-plate depths and this requires us to include the broad transition between brittle and ductile flow. Moreover, important changes in the ductile flow mechanisms occur at the intermediate temperatures found at mid-plate depths. Two specific contributions of laboratory rock rheology research are considered in this paper. First, the high-temperature steady-state flow mechanisms and rheology of mafic and ultramafic rocks are reviewed with special emphasis on olivine and crystalline rocks. Rock strength decreases very markedly with increases in temperature and it is the onset of flow by high temperature ductile mechanisms that defines the base of the lithosphere. The thickness of the continental lithosphere can therefore be defined by the depth to a particular isotherm Tc above which (at geologic strain rates) the high-temperature ductile strength falls below some arbitrary strength isobar (e.g., 100 MPa). For olivine Tc is about 700??-800??C but for other crustal silicates, Tc may be as low as 400??-600??C, suggesting that substantial decoupling may take place within thick continental crust and that strength may increase with depth at the Moho, as suggested by a number of workers on independent grounds. Put another way, the Moho is a rheological discontinuity. A second class of laboratory observations pertains to the general phenomenon of ductile faulting in which ductile strains are localized into shear zones. Ductile faults have been produced in experiments of five different rock types and is generally expressed as strain softening in constant-strain-rate tests or as an accelerating-creep-rate stage at constant differential stress. A number of physical mechanisms have been identified that may be responsible for ductile faulting, including the onset of dynamic recrystallization, phase changes, hydrothermal alteration and hydrolytic weakening. Microscopic evidence for these processes as well as larger-scale geological and geophysical observations suggest that ductile faulting in the middle to lower crust and upper mantle may greatly influence the distribution and magnitudes of differential stresses and the style of deformation in the overlying upper continental lithosphere.

Tectonophysics↗

USGS2007

Regional uplift associated with continental large igneous provinces: The roles of mantle plumes and the lithosphere

The timing and duration of surface uplift associated with large igneous provinces provide important constraints on mantle convection processes. Here we review geological indicators of surface uplift associated with five continent-based magmatic provinces: Emeishan Traps (260??million years ago: Ma), Siberian Traps (251??Ma), Deccan Traps (65??Ma), North Atlantic (Phase 1, 61??Ma and Phase 2, 55??Ma), and Yellowstone (16??Ma to recent). All five magmatic provinces were associated with surface uplift. Surface uplift can be measured directly from sedimentary indicators of sea-level in the North Atlantic and from geomorpholocial indicators of relative uplift and tilting in Yellowstone. In the other provinces, surface uplift is inferred from the record of erosion. In the Deccan, North Atlantic and Emeishan provinces, transient uplift that results from variations in thermal structure of the lithosphere and underlying mantle can be distinguished from permanent uplift that results from the extraction and emplacement of magma. Transient surface uplift is more useful in constraining mantle convection since models of melt generation and emplacement are not required for its interpretation. Observations of the spatial and temporal relationships between surface uplift, rifting and magmatism are also important in constraining models of LIP formation. Onset of surface uplift preceded magmatism in all five of the provinces. Biostratigraphic constraints on timing of uplift and erosion are best for the North Atlantic and Emeishan Provinces, where the time interval between significant uplift and first magmatism is less than 1??million years and 2.5??million years respectively. Rifting post-dates the earliest magmatism in the case of the North Atlantic Phase 1 and possibly in the case of Siberia. The relative age of onset of offshore rifting is not well constrained for the Deccan and the importance of rifting in controlling magmatism is disputed in the Emeishan and Yellowstone Provinces. In these examples, rifting is not a requirement for onset of LIP magmatism but melting rates are significantly increased when rifting occurs. Models that attempt to explain emplacement of these five LIPs without hot mantle supplied by mantle plumes often have difficulties in explaining the observations of surface uplift, rifting and magmatism. For example, small-scale convection related to craton or rift boundaries (edge-driven convection) cannot easily explain widespread (1000??km scale) transient surface uplift (Emeishan, Deccan, North Atlantic), and upper mantle convection initiated by differential incubation beneath cratons (the hotcell model) is at odds with rapid onset of surface uplift (Emeishan, North Atlantic). The start-up plume concept is still the most parsimonious way of explaining the observations presented here. However, observations of surface uplift cannot directly constrain the depth of origin of the hot mantle in a plume head. The short time interval between onset of transient surface uplift and magmatism in the North Atlantic and Emeishan means that the associated starting plume heads were probably not large (??? 1000??km diameter) roughly spherical diapirs and are likely to have formed narrow (??? 100??km radius) upwelling jets, with hot mantle then spreading rapidly outward within the asthenosphere. In cases where rifting post-dates magmatism (N Atlantic Phase 1) or where the degree of lithospheric extension may not have been great (Siberia), a secondary mechanism of lithospheric thinning, such as gravitational instability or delamination of the lower lithosphere, may be required to allow hot mantle to decompress sufficiently to explain the observed volume of magma with a shallow melting geochemical signature. Any such additional thinning mechanisms are probably a direct consequence of plume head emplacement. ?? 2007 Elsevier B.V. All rights reserved.

Chemical Geology↗

USGS2017

Lithospheric density structure beneath the Tarim basin and surroundings, northwestern China, from the joint inversion of gravity and topography

Intraplate strain generally focuses in discrete zones, but despite the profound impact of this partitioning on global tectonics, geodynamics, and seismic hazard, the processes by which deformation becomes localized are not well understood. Such heterogeneous intraplate strain is exemplified in central Asia, where the Indo-Eurasian collision has caused widespread deformation while the Tarim block has experienced minimal Cenozoic shortening. The apparent stability of Tarim may arise either because strain is dominantly accommodated by pre-existing faults in the continental suture zones that bound it—essentially discretizing Eurasia into microplates—or because the lithospheric-scale strength (i.e., viscosity) of the Tarim block is greater than its surroundings. Here, we jointly analyze seismic velocity, gravity, topography, and temperature to develop a 3-D density model of the crust and upper mantle in this region. The Tarim crust is characterized by high density, v s , v p , and v p /v s , consistent with a dominantly mafic composition and with the presence of an oceanic plateau beneath Tarim. Low-density but high-velocity mantle lithosphere beneath southern (southwestern) Tarim underlies a suite of Permian plume-related mafic intrusions and A-type granites sourced in previously depleted mantle lithosphere; we posit that this region was further depleted, dehydrated, and strengthened by Permian plume magmatism. The actively deforming western and southern margins of Tarim—the Tien Shan, Kunlun Shan, and Altyn Tagh fault—are underlain by buoyant upper mantle with low velocity; we hypothesize that this material has been hydrated by mantle-derived fluids that have preferentially migrated along Paleozoic continental sutures. Such hydrous material should be weak, and herein strain focuses there because of lithospheric-scale variations in rheology rather than the pre-existence of faults in the brittle crust. Thus this world-class example of strain partitioning arises not simply from the pre-existence of brittle faults but from the thermo-chemical and therefore rheological variations inherited from prior tectonism.

Earth and Planetary Science Letters↗

USGS2025

Exploring the disparate tectonic manifestations of lithospheric shortening on ocean worlds: Enceladus and Europa

The moons Europa and Enceladus both have water-ice dominated shells tens of kilometers thick overlying global oceans. Both moons are tidally heated and have a high rock mass fractions relative to other moons in their system, and both have youthful, tectonically deformed surfaces. Despite their many commonalities, the surfaces of the two moons bear little resemblance to one another. Europa’s tectonics are dominated by extensional features, including wide bands that imply extremely large surface strains. Features formed by the requisite lithospheric shortening accompanying such extension are largely absent, but long-wavelength folding may have occurred in some locations. Enceladus’ tectonics include features formed in both extension and contraction. Unlike Europa, folds, fold and thrust belts, and isolated thrust faults are plentiful there. In this investigation we use numerical simulations of the shortening of ice lithospheres to bring new insight into why two moons with so much in common appear so different from each other. We find that the temperature structure of the lithosphere controls the style of deformation that results from shortening. Warm conditions resulting from higher surface temperatures and heat fluxes produce folding, which is essentially a viscous response to shortening. Cold temperatures resulting from lower surface temperatures and heat fluxes results in faulting, which is a brittle response. Intermediate conditions can lead to complex combinations of folding and faulting. We hypothesize that Europa’s surface is generally warm enough for shortening to be accommodated by subtle folds that are difficult to detect in the sparse topography data currently available, whereas Enceladus’ colder temperatures favor a combination of folding and faulting. We argue that the dependence of tectonic style on lithospheric temperature self-consistently explains the pattern of tectonics observed across Enceladus, and we emphasize that non-rigid deformation (folding) is not inconsistent with the emerging paradigm of plate tectonics on Europa.

Icarus↗

USGS2004

Timing of Precambrian melt depletion and Phanerozoic refertilization events in the lithospheric mantle of the Wyoming Craton and adjacent Central Plains Orogen

Garnet peridotite xenoliths from the Sloan kimberlite (Colorado) are variably depleted in their major magmaphile (Ca, Al) element compositions with whole rock Re-depletion model ages generally consistent with this depletion occurring in the mid-Proterozoic. Unlike many lithospheric peridotites, the Sloan samples are also depleted in incompatible trace elements, as shown by the composition of separated garnet and clinopyroxene. Most of the Sloan peridotites have intermineral Sm-Nd and Lu-Hf isotope systematics consistent with this depletion occurring in the mid-Proterozoic, though the precise age of this event is poorly defined. Thus, when sampled by the Devonian Sloan kimberlite, the compositional characteristics of the lithospheric mantle in this area primarily reflected the initial melt extraction event that presumably is associated with crust formation in the Proterozoic-a relatively simple history that may also explain the cold geotherm measured for the Sloan xenoliths. The Williams and Homestead kimberlites erupted through the Wyoming Craton in the Eocene, near the end of the Laramide Orogeny, the major tectonomagmatic event responsible for the formation of the Rocky Mountains in the late Cretaceous-early Tertiary. Rhenium-depletion model ages for the Homestead peridotites are mostly Archean, consistent with their origin in the Archean lithospheric mantle of the Wyoming Craton. Both the Williams and Homestead peridotites, however, clearly show the consequences of metasomatism by incompatible-element-rich melts. Intermineral isotope systematics in both the Homestead and Williams peridotites are highly disturbed with the Sr and Nd isotopic compositions of the minerals being dominated by the metasomatic component. Some Homestead samples preserve an incompatible element depleted signature in their radiogenic Hf isotopic compositions. Sm-Nd tie lines for garnet and clinopyroxene separates from most Homestead samples provide Mesozoic or younger "ages" suggesting that the metasomatism occurred during the Laramide. Highly variable Rb-Sr and Lu-Hf mineral "ages" for these same samples suggest that the Homestead peridotites did not achieve intermineral equilibrium during this metasomatism. This indicates that the metasomatic overprint likely was introduced shortly before kimberlite eruption through interaction of the peridotites with the host kimberlite, or petrogenetically similar magmas, in the Wyoming Craton lithosphere. ?? 2004 Elsevier B.V. All rights reserved.

LITHOS↗

USGS2016

Asthenosphere–lithosphere interactions in Western Saudi Arabia: Inferences from 3 He/ 4 He in xenoliths and lava flows from Harrat Hutaymah

Extensive volcanic fields on the western Arabian Plate have erupted intermittently over the last 30 Ma following emplacement of the Afar flood basalts in Ethiopia. In an effort to better understand the origin of this volcanism in western Saudi Arabia, we analyzed 3 He/ 4 He, and He, CO 2 and trace element concentrations in minerals separated from xenoliths and lava flows from Harrat Hutaymah, supplemented with reconnaissance He isotope data from several other volcanic fields (Harrat Al Birk, Harrat Al Kishb and Harrat Ithnayn). Harrat Hutaymah is young (< 850 ka) and the northeasternmost of the volcanic fields. There is a remarkable homogeneity of 3 He/ 4 He trapped within most xenoliths, with a weighted mean of 7.54 ± 0.03 R A (2σ, n = 20). This homogeneity occurs over at least eight different xenolith types (including spinel lherzolite, amphibole clinopyroxenite, olivine websterite, clinopyroxenite and garnet websterite), and encompasses ten different volcanic centers within an area of ~ 2500 km 2 . The homogeneity is caused by volatile equilibration between the xenoliths and fluids derived from their host magma, as fluid inclusions are annealed during the infiltration of vapor-saturated magmas along crystalline grain boundaries. The notable exceptions are the anhydrous spinel lherzolites, which have a lower weighted mean 3 He/ 4 He of 6.8 ± 0.3 R A (2σ, n = 2), contain lower concentrations of trapped He, and have a distinctly depleted light rare earth element signature. 3 He/ 4 He values of ~ 6.8 R A are also commonly found in spinel lherzolites from harrats Ithnayn, Al Birk, and from Zabargad Island in the Red Sea. Olivine from non-xenolith-bearing lava flows at Hutaymah spans the He isotope range of the xenoliths. The lower 3 He/ 4 He in the anhydrous spinel lherzolites appears to be tied to remnant Proterozoic lithosphere prior to metasomatic fluid overprinting. Elevated 3 He/ 4 He in the western harrats has been observed only at Rahat (up to 11.8 R A ; Murcia et al., 2013), a volcanic field situated above thinned lithosphere beneath the Makkah-Medinah-Nafud volcanic lineament. Previous work established that spinel lherzolites at Hutaymah are sourced near the lithosphere-asthenosphere boundary (LAB), while other xenolith types there are derived from shallower depths within the lithosphere itself (Thornber, 1992). Helium isotopes are consistent with melts originating near the LAB beneath many of the Arabian harrats, and any magma derived from the Afar mantle plume currently appears to be of minor importance.

LITHOS↗

USGS2019

Lithospheric signature of late Cenozoic extension in electrical resistivity structure of the Rio Grande rift, New Mexico, USA

We present electrical resistivity models of the crust and upper mantle from two‐dimensional (2‐D) inversion of magnetotelluric (MT) data collected in the Rio Grande rift, New Mexico, USA. Previous geophysical studies of the lithosphere beneath the rift identified a low‐velocity zone several hundred kilometers wide, suggesting that the upper mantle is characterized by a very broad zone of modified lithosphere. In contrast, the surface expression of the rift (e.g., high‐angle normal faults and synrift sedimentary units) is confined to a narrow region a few tens of kilometers wide about the rift axis. MT data are uniquely suited to probing the depths of the lithosphere that fill the gap between surface geology and body wave seismic tomography, namely the middle to lower crust and uppermost mantle. We model the electrical resistivity structure of the lithosphere along two east‐west trending profiles straddling the rift axis at the latitudes of 36.2 and 32.0°N. We present results from both isotropic and anisotropic 2‐D inversions of MT data along these profiles, with a strong preference for the latter in our interpretation. A key feature of the anisotropic resistivity modeling is a broad (~200‐km wide) zone of enhanced conductivity (<20 Ωm) in the middle to lower crust imaged beneath both profiles. We attribute this lower crustal conductor to the accumulation of free saline fluids and partial melt, a direct result of magmatic activity along the rift. High‐conductivity anomalies in the midcrust and upper mantle are interpreted as fault zone alteration and partial melt, respectively.

New Mexico↗

USGS2010

Rheology of the lithosphere

During the quadrennial term 1979–1982, major advances have been made in our knowledge of the rheology of the oceanic lithosphere by the skillful combination of experimental and theoretical rock mechanics, seismology and marine geophysics in increasingly sophisticated models for the flexure of the oceanic lithosphere at seamounts and island chains, along transform faults, and at subduction zones. The relative simplicity of plate bending geometry, thermal history, and mineralogical and chemical compositions of the oceanic plates in these settings make the geophysical observations very powerful constraints on the in situ rheology of the oceanic lithosphere. In the first part of this paper, I review the laboratory work on materials appropriate to the oceanic lithosphere with emphasis on contributions during the quadrennial period and the need for future work. The important results of flexure models incorporating realistic material properties are then summarized.

Reviews of Geophysics and Space Physics↗