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 163 records · Page 9Linked to original sources

USGS2004

Shear wave velocity, seismic attenuation, and thermal structure of the continental upper mantle

Seismic velocity and attenuation anomalies in the mantle are commonly interpreted in terms of temperature variations on the basis of laboratory studies of elastic and anelastic properties of rocks. In order to evaluate the relative contributions of thermal and non-thermal effects on anomalies of attenuation of seismic shear waves, Q −1 s , and seismic velocity, V s , we compare global maps of the thermal structure of the continental upper mantle with global Q −1 s and V s maps as determined from Rayleigh waves at periods between 40 and 150 s. We limit the comparison to three continental mantle depths (50, 100 and 150 km), where model resolution is relatively high. The available data set does not indicate that, at a global scale, seismic anomalies in the upper mantle are controlled solely by temperature variations. Continental maps have correlation coefficients of <0.56 between V s and T and of <0.47 between Q s and T at any depth. Such low correlation coefficients can partially be attributed to modelling artefacts; however, they also suggest that not all of the V s and Q s anomalies in the continental upper mantle can be explained by T variations. Global maps show that, by the sign of the anomaly, V s and Q s usually inversely correlate with lithospheric temperatures: most cratonic regions show high V s and Q s and low T , while most active regions have seismic and thermal anomalies of the opposite sign. The strongest inverse correlation is found at a depth of 100 km, where the attenuation model is best resolved. Significantly, at this depth, the contours of near-zero Q s anomalies approximately correspond to the 1000 °C isotherm, in agreement with laboratory measurements that show a pronounced increase in seismic attenuation in upper mantle rocks at 1000–1100 °C. East–west profiles of V s , Q s and T where continental data coverage is best (50°N latitude for North America and 60°N latitude for Eurasia) further demonstrate that temperature plays a dominant, but non-unique, role in determining the value of lithospheric V s and Q s . At 100 km depth, where the resolution of seismic models is the highest, we compare observed seismic V s and Q s with theoretical V T s and Q T s values, respectively, that are calculated solely from temperature anomalies and constrained by experimental data on temperature dependencies of velocity and attenuation. This comparison shows that temperature variations alone are sufficient to explain seismic V s and Q s in ca 50 per cent of continental regions. We hypothesize that compositional anomalies resulting from Fe depletion can explain the misfit between seismic and theoretical V s in cratonic lithosphere. In regions of active tectonics, temperature effects alone cannot explain seismic V s and Q s in the lithosphere. It is likely that partial melts and/or fluids may affect seismic parameters in these regions. This study demonstrates that lithospheric temperature plays the dominant role in controlling V s and Q s anomalies, but other physical parameters, such as compositional variations, fluids, partial melting and scattering, may also play a significant role in determining V s and Q s variations in the continental mantle.

Geophysical Journal International↗

USGS2007

Propagation of the Hawaiian-Emperor volcano chain by Pacific plate cooling stress

The lithosphere crack model, the main alternative to the mantle plume model for age-progressive magma emplacement along the Hawaiian-Emperor volcano chain, requires the maximum horizontal tensile stress to be normal to the volcano chain. However, published stress fields calculated from Pacific lithosphere tractions and body forces (e.g., subduction pull, basal drag, lithosphere density) are not optimal for southeast propagation of a stress-free, vertical tensile crack coincident with the Hawaiian segment of the Hawaiian-Emperor chain. Here we calculate the thermoelastic stress rate for present-day cooling of the Pacific plate using a spherical shell finite element representation of the plate geometry. We use observed seafloor isochrons and a standard model for lithosphere cooling to specify the time dependence of vertical temperature profiles. The calculated stress rate multiplied by a time increment (e.g., 1 m.y.) then gives a thermoelastic stress increment for the evolving Pacific plate. Near the Hawaiian chain position, the calculated stress increment in the lower part of the shell is tensional, with maximum tension normal to the chain direction. Near the projection of the chain trend to the southeast beyond Hawaii, the stress increment is compressive. This incremental stress field has the form necessary to maintain and propagate a tensile crack or similar lithosphere flaw and is thus consistent with the crack model for the Hawaiian volcano chain.

Special Paper of the Geological Society of America↗

USGS2023

Breaking plates: Creation of the East Anatolian fault, the Anatolian plate, and a tectonic escape system

Lateral movement of lithospheric fragments along strike-slip faults in response to collision (escape tectonics) has characterized convergent settings since the onset of plate tectonics and is a mechanism for the formation of new plates. The Anatolian plate was created by the sequential connection of strike-slip faults following ≥10 m.y. of distributed deformation that ultimately localized into plate-bounding faults. Thermochronology data and seismic images of lithosphere structure near the East Anatolian fault zone (EAFZ) provide insights into the development of the new plate and escape system. Low-temperature thermochronology ages of rocks in and near the EAFZ are significantly younger than in other fault zones in the region, e.g., apatite (U-Th)/He: 11–1 Ma versus 27–13 Ma. Young apatite (U-Th)/He ages and thermal history modeling record thermal resetting along the EAFZ over the past ~5 m.y. and are interpreted to indicate thermal activity triggered by strike-slip faulting in the EAFZ as it formed as a through-going, lithosphere-scale structure. The mechanism for EAFZ formation may be discerned from S-wave velocity images from the Continental Dynamics–Central Anatolian Tectonics (CD-CAT) seismic experiment. These images indicate that thin but strong Arabian lithospheric mantle extends ~50–150 km north beneath Anatolian crust and would have been located near the present surficial location of the Bitlis-Zagros suture zone (co-located with the EAFZ in our study area) at ca. 5 Ma. Underthrusting of strong Arabian lithosphere facilitated localization of the EAFZ and thus was a fundamental control on the formation of the Anatolian plate and escape system.

East Anatolian Fault↗

USGS2023

Mantle origin and crustal differentiation of basalts and hawaiites of northern Harrat Rahat, Kingdom of Saudi Arabia

Quaternary volcanic rocks of northern Harrat Rahat, Kingdom of Saudi Arabia, are chiefly alkali basalts with subordinate transitional basalts, hawaiites, mugearites, benmoreites, and trachytes. Geochemical and isotopic results indicate that crystallization-differentiation, mixing, and cumulate reassimilation within the magmatic system produced most of its compositional diversity, with only minor involvement of Neoproterozoic crust. With increasing evolution, crystal assemblages that separated from and drove basalt-to-hawaiite differentiation passed from (1) dunitic or troctolitic to (2) olivine gabbroic to (3) titanomagnetite-bearing olivine gabbroic, with typical hawaiites representing about 20 weight percent residual liquids from an estimated primary alkali basaltic parent. Crystallization-differentiation simulations for midcrustal pressures yield the closest compositional match to the basalt-hawaiite suite, and stagnation in the midcrustal area (near 20 kilometers [km] depth) may result from density trapping beneath the lower density upper continental crust. Least differentiated alkali basalts have fractionation-adjusted major-oxide compositions that are consistent with separation from the deeper parts (60–70 km) of the spinel-lherzolite stability field at pressures that are close to the local lithosphere-asthenosphere boundary (60–80 km). Mantle potential temperature estimates are strongly sensitive to modeling approach, but potential temperatures (1,345–1,390 degrees Celsius [°C]) are not discernably greater than for midocean ridge basalts (MORB; 1,350–1,410 °C) if adherence to spinel-lherzolite melting relations is required. Inversion of the trace-element concentrations of the lesser fractionated basalts indicates a depleted mantle source, similar to MORB-source estimates, but one that is enriched in Sr and includes a greater relative proportion of melting in the garnet-peridotite stability field. These geochemical and thermal relations, as well as radiogenic isotopes, point to a dominantly depleted mantle asthenospheric source for Harrat Rahat basalts, admixed with subordinate materials either from the Afar mantle plume or an enriched MORB component in the ambient asthenosphere. The lithosphere-asthenosphere boundary is shallower beneath the belt of major volcanic fields on the Arabia Plate, and restoration of rifting across the Red Sea and Gulf of Aden places the south end of this belt adjacent to the northern part of the Afar region, suggesting a once-continuous structure, possibly an arch, a weakness, or a discontinuity along the base of the lithosphere. Magma generation can be ascribed to focused upwelling and decompression melting, perhaps driven by a magmatic-feedback mechanism whereby basaltic intrusions into the deep lithosphere solidify as eclogites, causing lithospheric foundering and further asthenospheric upwelling and decompression melting in a restricted region.

northern Harrat Rahat↗

USGS2026

Toward a four-dimensional petrogenetic model of a distributed volcanic field on the southern edge of the Colorado Plateau

A detailed characterization of the >3,000 square kilometer (km 2 ) Springerville volcanic field, located on the southern tip of the Colorado Plateau in Arizona, United States, with its more than 501 volcanic units and widely distributed >420 cinder cones and lava flows, provides constraints toward an integrated petrogenetic model for the field. Large-volume effusive tholeiitic eruptions at 2–1.5 mega-annum (Ma) transitioned to more numerous, smaller volume alkali olivine basalt (AOB) events at 1.5–1.0 Ma, with increasing abundances of evolved alkalic rocks (EARs), and a final 1.0–0.3 Ma period dominated by smaller volume, more explosive alkalic eruptions. Early large-scale melting in a relatively enriched lithospheric mantle (EM) source generated large-volume effusive tholeiitic magmas. Depths of tholeiite magma generation average about 90 kilometers (km) across the field, but depths for individual units decreased southward, consistent with lithospheric thinning toward the Colorado Plateau margin. Early and middle-stage transitional basalts, alkali olivine basalts (AOBs), and basanites originate from a progressively deeper (>100 km) region in a prevalent mantle (PREMA)-like asthenospheric source produced by increasingly smaller degrees of melting, as low as about 2 percent. The chemical signature of the basanites is consistent with small degrees of melting in a carbonated, asthenospheric source to depths of about 140 km. As heat waned, the last phase of volcanism was dominated by more explosive EARs derived at shallower lithospheric pressures but that have isotopic and trace element similarities to the deeper asthenospheric magmas. This suggests mixing between deeper basanitic and shallower tholeiitic magmas. With waning heat, eruptions became more localized along alignments, likely related to boundaries between blocks of Proterozoic crust with differing properties that affected magma ascent. The petrogenetic patterns are consistent with a variety of processes. Basin and Range Province extension, melting, and heat-induced weakening progressively eroded the Colorado Plateau’s thicker lithosphere, giving rise to relatively high degrees of partial melting from shallower (<90 km) sources that produced the early large volumes of effusive tholeiitic magma. Possible lithosphere delamination and removal, and a resulting steep boundary step with the asthenosphere, may have given rise to mantle edge convection, resulting in the mixing of basanitic and tholeiitic magmas. With plate motion, shear-driven upwelling likely gave an eastward component to convection, affecting relative amounts of melting at the field’s western boundary and corresponding with volcanism moving eastward during times of peak vent production by eruption of alkalic lavas, giving rise to more explosive and lower volume eruptions.

Arizona, New Mexico↗

USGS1992

ALVIN investigation of an active propagating rift system, Galapagos 95.5° W

ALVIN investigations have defined the fine-scale structural and volcanic patterns produced by active rift and spreading center propagation and failure near 95.5° W on the Galapagos spreading center. Behind the initial lithospheric rifting, which is propagating nearly due west at about 50 km m.y. −1 , a triangular block of preexisting lithosphere is being stretched and fractured, with some recent volcanism along curving fissures. A well-organized seafloor spreading center, an extensively faulted and fissured volcanic ridge, develops ~ 10 km (~ 200,000 years) behind the tectonic rift tip. Regional variations in the chemical compositions of the youngest lavas collected during this program contrast with those encompassing the entire 3 m.y. of propagation history for this region. A maximum in degree of magmatic differentiation occurs about 9 km behind the propagating rift tip, in a region of diffuse rifting. The propagating spreading center shows a gentle gradient in magmatic differentiation culminating at the SW-curving spreading center tip. Except for the doomed rift, which is in a constructional phase, tectonic activity also dominates over volcanic activity along the failing spreading system. In contrast to the propagating rift, failing rift lavas show a highly restricted range of compositions consistent with derivation from a declining upwelling zone accompanying rift failure. The lithosphere transferred from the Cocos to the Nazca plate by this propagator is extensively faulted and characterized by ubiquitous talus in one of the most tectonically disrupted areas of seafloor known. The pseudofault scarps, where the preexisting lithosphere was rifted apart, appear to include both normal and propagator lavas and are thus more lithologically complex than previously thought. Biological communities, probably vestimentiferan tubeworms, occur near the top of the outer pseudofault scarp, although no hydrothermal venting was observed.

Marine Geophysical Research↗

USGS1993

Isotopic and trace-element constraints on mantle and crustal contributions to Siberian continental flood basalts, Noril'sk area, Siberia

We present a tightly controlled and comprehensive set of analytical data for the 250-Ma Siberian flood-basalt province. Consideration of major- and trace-element compositions, along with strontium, lead and neodymium isotopic compositions, strongly supports earlier Russian subdivision of this magmatism into three magmatic cycles, giving rise to three assemblages of eleven basalt suites in the ascending order Ivakinsky-Gudchikhinsky, Khakanchansky-Nadezhdinsky and Morongovsky-Samoedsky. Geochemical and isotopic discontinuities of varying magnitude characterize most of the boundaries between the eleven recognized basalt suites in the Noril'sk area. Although we conclude that the dominant volume of erupted magma originated from an asthenospheric mantle plume, none of the lavas is interpreted to directly represent asthenospheric melts, which would have been far more magnesian. On the basis of thermal considerations, we consider it unlikely that vast volumes of basaltic melt were produced directly from the continental lithospheric mantle beneath the Siberian craton. Moreover, there is little evidence from mantle xenoliths that the geochemical signatures of such melts would correspond to those of the Siberian flood basalts. Studies of melt migration lead us to conclude that transport of asthenospheric melt through the lithospheric mantle would be rapid, by fracture propagation. Lavas from the Gudchikhinsky suite have negligible Ta-Nb anomalies and positive ϵ Nd values and their parental magmas presumably interacted little with the continental lithospheric mantle or crust. All other lavas have negative Ta-Nb anomalies and lower ϵ Nd values that we attribute to interaction with continental crust. The model that we have developed requires discrete contributions from the plume and complex processing of all erupted magmas in the continental crust. The earliest magmas represent small percentages of melt formed in equilibrium with garnet. Over time, the percentage of melting in the source region and the volume of magma produced increased, and garnet was no longer stable in the plume source. All of the plume-derived melts initially contained more than 20 wt% MgO and became less Mg rich by fractionation of olivine as they traversed the lithospheric mantle. We conclude, however, that the most significant control on the geochemical and isotopic compositions of all the erupted lavas was processing of mantle-derived magma in crustal reservoirs during periodic replenishment, periodic tapping, continuous crystal fractionation and wallrock assimilation. Rapid eruption of an extremely large volume of processed magma that varied little in chemical and isotopic composition produced the sequence of relatively monotonous tholeiitic basalts that constitute the 2,300-m-thick third assemblage of the Siberian flood-basalt province near Noril'sk.

Geochimica et Cosmochimica Acta↗

USGS2022

Seismic evidence for magmatic underplating along the Kodiak-Bowie Seamount Chain, Gulf of Alaska

Oceanic crust formed at mid-ocean ridges may be later modified by off-ridge magmatism forming seamounts, guyots, and islands. We investigate processes associated with seamount formation in the Gulf of Alaska Seamount Province using two coincident seismic reflection/wide-angle profiles. A north-south profile crosses the Kodiak-Bowie Seamount Chain and Aja fracture zone (FZ), and an orthogonal east-west profile is located about 90 km south of the seamount chain over Pacific plate oceanic crust. Structure along the profile away from the seamount chain is consistent with typical oceanic crust. Crust in our study region is thinnest (about 5.6 km) at the Aja FZ. Unlike observations from active transform faults, no low-velocity anomaly is observed at the Aja FZ suggesting that the crustal velocities have recovered to normal values through crack closure and crack healing. Higher lower crustal velocities (∼7.3 and > 7.5 km/s) and thicker crust (∼8.5 and ∼7.0 km) are observed near the Pratt and Durgin Seamounts and at the intersection of the Kodiak-Bowie Seamount Chain linear trend, respectively. These observations are attributed to magmatic underplating associated with seamount province magmatism. Lithospheric thickness variations across the Aja FZ may form a barrier or impediment to magmatic flow. The thickest crust (8.5 km) along our two profiles is located on the younger side of the FZ, and we suggest that the majority of magmatism jumped south of the Aja FZ when thinner lithosphere was encountered by the Bowie hot spot. The crustal structure near the Kodiak-Bowie Seamount Chain is most similar to that of other seamounts and guyots that formed on similarly young lithosphere (8–12 Ma). Our results suggest that lithospheric thickness at the time of hot spot interaction has a large control on magmatic underplating at seamounts and seamount provinces.

Alaska↗

USGS2003

Density of the continental roots: Compositional and thermal contributions

The origin and evolution of cratonic roots has been debated for many years. Precambrian cratons are underlain by cold lithospheric roots that are chemically depleted. Thermal and petrologic data indicate that Archean roots are colder and more chemically depleted than Proterozoic roots. This observation has led to the hypothesis that the degree of depletion in a lithospheric root depends mostly on its age. Here we test this hypothesis using gravity, thermal, petrologic, and seismic data to quantify differences in the density of cratonic roots globally. In the first step in our analysis we use a global crustal model to remove the crustal contribution to the observed gravity. The result is the mantle gravity anomaly field, which varies over cratonic areas from -100 to +100 mGal. Positive mantle gravity anomalies are observed for cratons in the northern hemisphere: the Baltic shield, East European Platform, and the Siberian Platform. Negative anomalies are observed over cratons in the southern hemisphere: Western Australia, South America, the Indian shield, and Southern Africa. This indicates that there are significant differences in the density of cratonic roots, even for those of similar age. Root density depends on temperature and chemical depletion. In order to separate these effects we apply a lithospheric temperature correction using thermal estimates from a combination of geothermal modeling and global seismic tomography models. Gravity anomalies induced by temperature variations in the uppermost mantle range from -200 to +300 mGal, with the strongest negative anomalies associated with mid-ocean ridges and the strongest positive anomalies associated with cratons. After correcting for thermal effects, we obtain a map of density variations due to lithospheric compositional variations. These maps indicate that the average density decrease due to the chemical depletion within cratonic roots varies from 1.1% to 1.5%, assuming the chemical boundary layer has the same thickness as the thermal boundary layer. The maximal values of the density drop are in the range 1.7-2.5%, and correspond to the Archean portion of each craton. Temperatures within cratonic roots vary strongly, and our analysis indicates that density variations in the roots due to temperature are larger than the variations due to chemical differences.

Earth and Planetary Science Letters↗

USGS2000

Quantifying precambrian crustal extraction: The root is the answer

We use two different methods to estimate the total amount of continental crust that was extracted by the end of the Archean and the Proterozoic. The first method uses the sum of the seismic thickness of the crust, the eroded thickness of the crust, and the trapped melt within the lithospheric root to estimate the total crustal volume. This summation method yields an average equivalent thickness of Archean crust of 49±6 km and an average equivalent thickness of Proterozoic crust of 48± 9 km. Between 7 and 9% of this crust never reached the surface, but remained within the continental root as congealed, iron-rich komatiitic melt. The second method uses experimental models of melting, mantle xenolith compositions, and corrected lithospheric thickness to estimate the amount of crust extracted through time. This melt column method reveals that the average equivalent thickness of Archean crust was 65±6 km, and the average equivalent thickness of Early Proterozoic crust was 60±7 km. It is likely that some of this crust remained trapped within the lithospheric root. The discrepancy between the two estimates is attributed to uncertainties in estimates of the amount of trapped, congealed melt, overall crustal erosion, and crustal recycling. Overall, we find that between 29 and 45% of continental crust was extracted by the end of the Archean, most likely by 2.7 Ga. Between 51 and 79% of continental crust was extracted by the end of the Early Proterozoic, most likely by 1.8–2.0 Ga. Our results are most consistent with geochemical models that call upon moderate amounts of recycling of early extracted continental crust coupled with continuing crustal growth (e.g. McLennan, S.M., Taylor, S.R., 1982. Geochemical constraints on the growth of the continental crust. Journal of Geology, 90, 347–361; Veizer, J., Jansen, S.L., 1985. Basement and sedimentary recycling — 2: time dimension to global tectonics. Journal of Geology 93(6), 625–643). Trapped, congealed, iron-rich melt within the lithospheric root may represent some of the iron that is ‘missing’ from the lower crust. The lower crust within Archean cratons may also have an unexpectedly low iron content because it was extracted from more primitive, undepleted mantle.

Tectonophysics↗

USGS2019

Azimuthal seismic anisotropy of 70 Ma Pacific‐plate upper mantle

Plate formation and evolution processes are predicted to generate upper mantle seismic anisotropy and negative vertical velocity gradients in oceanic lithosphere. However, predictions for upper mantle seismic velocity structure do not fully agree with the results of seismic experiments. The strength of anisotropy observed in the upper mantle varies widely. Further, many refraction studies observe a fast direction of anisotropy rotated several degrees with respect to the paleospreading direction, suggesting that upper mantle anisotropy records processes other than 2D corner flow and plate‐driven shear near mid‐ocean ridges. We measure 6.0 ± 0.3% anisotropy at the Moho in 70 Ma lithosphere in the central Pacific with a fast direction parallel to paleospreading, consistent with mineral alignment by 2D mantle flow near a mid‐ocean ridge. We also find an increase in the strength of anisotropy with depth, with vertical velocity gradients estimated at 0.02 km/s/km in the fast direction and 0 km/s/km in the slow direction. The increase in anisotropy with depth can be explained by mechanisms for producing anisotropy other than intrinsic effects from mineral fabric, such as aligned cracks or other structures. This measurement of seismic anisotropy and gradients reflects the effects of both plate formation and evolution processes on seismic velocity structure in mature oceanic lithosphere, and can serve as a reference for future studies to investigate the processes involved in lithospheric formation and evolution.

Journal of Geophysical Research B: Solid Earth↗

USGS2022

Silicate volcanism on Europa’s seafloor and implications for habitability

Habitable ocean environments on Europa require an influx of reactants to maintain chemical disequilibrium. One possible source of reactants is seafloor volcanism. Modeling has shown that dissipation of tidal energy in Europa's asthenosphere can generate melt, but melt formation cannot be equated with volcanism. Melt must also be transported through Europa's cold lithosphere to erupt at the seafloor. Here, we use two models of dike propagation to show that dikes can only traverse the lithosphere if either the fracture toughness of the lithosphere or the flux into the dike is large (>500 MPa m 1/2 or ∼1 m 2 s −1 , respectively). We conclude that cyclic volcanic episodes might provide reactants to Europa's ocean if magma accumulates at the base of the lithosphere for several thousand years. However, if dikes form too frequently, or are too numerous, the magma flux into each will be insufficient, and volcanism cannot support a habitable ocean environment.

Geophysical Research Letters↗

USGS2012

Rift flank uplifts and Hinterland Basins: Comparison of the Transantarctic Mountains with the Great Escarpment of southern Africa

Uplifted rift margins are a common feature of continents and oceans. Two variants of rift flank morphologies have been recognized: One in which the topography warps down from an inland high toward the continental margin, and one where the tropographic peak lies close to the continental margin. The Great Escarpment of southern Africa and the Transantarctic Mountains are examples of the first and the second variants of rift flanks, respectively. Both rift flanks are bordered on their landward side by broad continental basins: the Kalahari and the Wilkes hinterland basins. If these basins are interpreted as flexural “outer lows” that deepen in unison with the uplift of the rift flanks, the lithosphere on the uplifted side is very rigid in both cases (elastic thickness T e of 100 ± 20 km for southern Africa and 110 ± 20 km for East Antarctica). We suggest that the variation in rift flank morphology is caused by the isostatic response to uplift forces of elastic plates sharing different boundary conditions. We model the uplift of the Transantarctic Mountains as an upward deflection of an elastic plate which is broken at the front of the Transantarctic Mountains, and we model the uplift of the Great Escarpment as an upward deflection of a continuous elastic plate that is modified by the downward load of sediments on the continental margin. Although the Transantarctic Mountain uplift is young (60–0 Ma) and the southern African uplift is old (<100 Ma), the different isostatic responses of the two margins are not a function of age, because most loading (sedimentation) and unloading (erosion) took place shortly after rifting. Detailed modeling of topography, gravity, geological markers, and the locations of depocenters suggests that lithospheric rigidity decreases under the Transantarctic Mountains, whereas in southern Africa the decrease occurs not under the Great Escarpment but far seaward under the continental shelf and slope. If the distribution of lithospheric rigidity is indicative of the thermal regime of the lithosphere, then uplifted rift flanks are not always underlain by a thermal anomaly. This and other geological evidence indicate that a single mechanism cannot explain the uplift of both the Antarctic and the African margins.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Nonlinear strain buildup and the earthquake cycle on the San Andreas fault

Two contrasting models of the earthquake deformation cycle on strike slip faults predict significant temporal declines in shear strain rate near the fault, accompanied by a progressive broadening of the zone of deformation adjacent to it. In the thin lithosphere model, transient deformation results from flow in the asthenosphere due to stress relaxation following faulting through most or all of the lithosphere. For an earth model with a thick elastic lithosphere (plate thickness » depth of seismic slip), transient motions are due to postearthquake aseisrnic slip below the coseismic fault plane. Data from the San Andreas fault indicate a long-term temporal decrease in strain rate that persists for at least 30 years and may extend through the entire earthquake cycle. Observations support a cycle-long rate decrease and a temporal spreading of the deformation profile only if movement cycles on the northern and southern locked sections of the fault are basically similar. If so, the usually lower strain rates and broader deformation zone currently observed on the southern San Andreas represent a later evolutionary stage of the northern locked section, where a great earthquake is a more recent occurrence. Although the data allow some extreme models to be discarded, no sufficiently strong constraints exist to decide between the thin and thick lithosphere models. Regardless of the appropriate model the geodetic observations themselves indicate that strain buildup is sufficiently nonlinear to cause significant departures from recurrence estimates based on linear strain accumulation and the time-predictable model.

Journal of Geophysical Research Solid Earth↗

USGS2012

A plate flexure approximation to postseismic and interseismic deformation

The rather large postseismic deformation that is associated with two‐dimensional dip‐slip faulting in the lithosphere is related to the bending of a free plate generated by dip‐slip faulting. In the absence of gravity, asthenosphere relaxation eventually permits the faulted lithosphere to assume the dihedral configuration of a faulted free plate. For thrust faulting, the faulted area is depressed into the asthenosphere, and the flanks of the plate slope uniformly upward. In the presence of gravity, buoyancy forces act upon the plate, and the ultimate ( t →∞) postseismic uplift is approximated ( y /α > 0.2) by w 1 exp (− y /α) (cos y /α ‐ sin y /α), where w 1 is an explicit function of the fault parameters and lithosphere properties, y is the horizontal distance from the downdip end of the fault, and α is the flexural parameter for the lithosphere plate. The relaxed‐asthenosphere response is the sum of the coseismic deformation (unrelaxed‐astheriosphere response) and this ultimate postseismic deformation. The annual, steady deformation associated with strain accumulation at a subduction zone is simply the relaxed‐asthenosphere response to virtual, normal slip on the main thrust zone of an amount equal to the annual plate convergence. From these relations we have estimated the deformation (strain and uplift) expected along the Pacific coast of northern Honshu. The estimated deformation is about twice the observed deformation, but the predicted geographical distribution of the deformation is similar to the observed deformation.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Midcontinent rift volcanism in the Lake Superior region: Sr, Nd, and Pb isotopic evidence for a mantle plume origin

Between 1091 and 1098 Ma, most of a 15- to 20-km thickness of dominantly tholeiitic basalt erupted in the Midcontinent Rift System of the Lake Superior region, North America. The Portage Lake Volcanics in Michigan, which are the youngest MRS flood basalts, fall into distinctly high- and low-TiO 2 types having different liquid lines of descent. Incompatible trace elements in both types of tholeiites are enriched compared to depleted or primitive mantle (La/Yb = 4.3–5.3; Th/Ta = 2.12–2.16; Zr/Y = 4.3–4.4), and both basalt types are isotopically indistinguishable. Sr, Nd, and Pb isotopic compositions of the Portage Lake tholeiites have 87 Sr/ 86 Sr i ≈ 0.7038, ε Nd (1095 Ma) ≈ 0±2, and μ 1 ≈ 8.2. Model ages with respect to a depleted mantle source ( T DM ) average about 1950–2100 Ma. Portage Lake rhyolites fall into two groups. Type I rhyolites have Nd and Pb isotopic characteristics (ε Nd (1095 Ma) ≈ 0 to −4.7; μ 1 ≈ 8.2–7.8) consistent with contamination of tholeiitic rocks by 5–10% Archean crust. The one type II rhyolite analyzed has Nd and Pb isotopic compositions (ε Nd (1095 Ma) ≈ −13 to −16; μ 1 ≈ 7.6–7.7) which are consistent with partial melting of Archean crust. Early Proterozoic crust was not a major contaminant of MRS rocks in the Lake Superior region. Most reported Nd and Pb isotopic compositions of MRS tholeiites from the main stage of volcanism in the Lake Superior region and of the Duluth Complex are comparable to the Nd and Pb isotopic data for Portage Lake tholeiites. The isotopic enrichment of the MRS source compared to depleted mantle is striking and must have occurred at least 700 m.y. before 1100 Ma. There are two likely sources for such enriched MRS tholeiitic magmatism: subcontinental lithospheric mantle enriched during the early Proterozoic or enriched mantle derived from an upwelling plume. Subcontinental lithospheric mantle alone as a source cannot be ruled out strictly on isotopic grounds, but melting of enriched lithosphere potentially ranging in age from 3.6 to 1.8 Ga would have been less likely to generate either the homogeneity or the volume of MRS tholeiites. Decompression melting of an upwelling enriched mantle plume in a region of lithosphere thinned by extension could have successfully generated the enormous volume (850×10 3 km 3 ) of relatively homogeneous magma in a restricted time interval.

Journal of Geophysical Research Solid Earth↗

USGS2022

Geophysical extent of the Wyoming Province, western USA: Insights into ancient subduction and craton stability

A new 3-D resistivity model, estimated from inversion of magnetotelluric data, images crustal and upper-mantle structure of the Wyoming Province and adjacent areas. The Archean province is imaged as a coherent resistive domain, in sharp contrast to active tectonic domains of the western U.S. Prominent high-conductivity belts define the northern, eastern, and southern margins of the Wyoming Province and are interpreted as sutures marking the remnants of Paleoproterozoic orogens. The model results suggest the northern boundary of the Wyoming Province is located 150 km south of its traditional placement and adjacent to a composite orogen separating the Wyoming Province and Medicine Hat block. The eastern province boundary is clearly imaged along the Black Hills, whereas the western margin is obscured by Cenozoic extension and magmatism. An internal boundary within the Wyoming Province is interpreted to represent a Neoarchean suture; in stark contrast to Proterozoic sutures, though, it is not marked by a high-conductivity belt. This difference in conductivity is speculated to reflect changes in the subduction process through time. The absence of high-conductivity along Archean sutures appears to be global in nature and related to reduced continental freeboard in the Archean which limited continental weathering and the delivery of carbon-rich sediments to the seafloor. Although the entire Wyoming Province has been proposed to have undergone lithospheric modification that lessened its stability, the resistivity model suggests a thick lithospheric root remains in place except along its western margin. These results suggest that Archean cratons may be more resistant to lithospheric modification by influx of heat and fluids associated with extension and plumes than previously thought, and that metasomatism does not necessarily weaken the lithosphere and set a craton on the path to destruction.

Colorado, Montana, Utah, Wyoming↗

USGS2007

Thermal structure of oceanic transform faults

We use three-dimensional finite element simulations to investigate the temperature structure beneath oceanic transform faults. We show that using a rheology that incorporates brittle weakening of the lithosphere generates a region of enhanced mantle upwelling and elevated temperatures along the transform; the warmest temperatures and thinnest lithosphere are predicted to be near the center of the transform. Previous studies predicted that the mantle beneath oceanic transform faults is anomalously cold relative to adjacent intraplate regions, with the thickest lithosphere located at the center of the transform. These earlier studies used simplified rheologic laws to simulate the behavior of the lithosphere and underlying asthenosphere. We show that the warmer thermal structure predicted by our calculations is directly attributed to the inclusion of a more realistic brittle rheology. This temperature structure is consistent with a wide range of observations from ridge-transform environments, including the depth of seismicity, geochemical anomalies along adjacent ridge segments, and the tendency for long transforms to break into small intratransform spreading centers during changes in plate motion.

Geology↗