USGS Science⌕ Search

SEARCH · USGS Science

Results for “Lithosphere”

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

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18Linked to original sources

Growth and degradation of Hawaiian volcanoes

The 19 known shield volcanoes of the main Hawaiian Islands—15 now emergent, 3 submerged, and 1 newly born and still submarine—lie at the southeast end of a long-lived hot spot chain. As the Pacific Plate of the Earth’s lithosphere moves slowly northwestward over the Hawaiian hot spot, volcanoes are successively born above it, evolve as they drift away from it, and eventually die and subside beneath the ocean surface. The massive outpouring of lava flows from Hawaiian volcanoes weighs upon the oceanic crust, depressing it by as much as 5 km along an axial Hawaiian Moat. The periphery of subsidence is marked by the surrounding Hawaiian Arch. Subsidence is ongoing throughout almost all of a volcano’s life. During its active life, an idealized Hawaiian volcano passes through four eruptive stages: preshield, shield, postshield, and rejuvenated. Though imperfectly named, these stages match our understanding of the growth history and compositional variation of the Hawaiian volcanoes; the stages reflect variations in the amount and rate of heat supplied to the lithosphere as it overrides the hot spot. Principal growth occurs in the first 1–2 million years as each volcano rises from the sea floor or submarine flank of an adjacent volcano. Volcanic extinction ensues as a volcano moves away from the hot spot. Eruptive-stage boundaries are drawn somewhat arbitrarily because of their transitional nature. Preshield-stage lava is alkalic as a consequence of a nascent magma-transport system and less extensive melting at the periphery of the mantle plume fed by the hot spot. The shield stage is the most productive volcanically, and each Hawaiian volcano erupts an estimated 80–95 percent of its ultimate volume in tholeiitic lavas during this stage. Shield-stage volcanism marks the time when a volcano is near or above the hot spot and its magma supply system is robust. This most active stage may also be the peak time when giant landslides modify the flanks of the volcanoes, although such processes begin earlier and extend later in the life of the volcanoes. Late-shield strata extend the silica range as alkali basalt and even hawaiite lava flows are sparsely interlayered with tholeiite at some volcanoes. Rare are more highly fractionated shield-stage lava flows, which may reach 68 weight percent SiO 2 . Intervolcano compositional differences result mainly from variations in the part of the mantle plume sampled by magmatism and the distribution of magma sources within it. Volcanism wanes gradually as Hawaiian volcanoes move away from the hot spot, passing from the shield stage into the postshield stage. Shallow magma reservoirs (1–7-km depth) of the shield-stage volcanoes cannot be sustained as magma supply lessens, but smaller reservoirs at 20–30-km depth persist. The rate of extrusion diminishes by a factor of 10 late in the shield stage, and the composition of erupted lava becomes more alkalic—albeit erratically—as the degree of melting diminishes. The variation makes this transition, from late shield to postshield, difficult to define rigorously. Of the volcanoes old enough to have seen this transition, eight have postshield strata sufficiently distinct and widespread to map separately. Only two, Ko‘olau and Lāna‘i, lack rocks of postshield composition. Five Hawaiian volcanoes have seen rejuvenated-stage volcanism following quiescent periods that ranged from 2.0 to less than 0.5 million years. The rejuvenated stage can be brief—only one or two eruptive episodes—or notably durable. That on Ni‘ihau lasted from 2.2 to 0.4 million years ago; on Kaua‘i, the stage has been ongoing since 3.5 million years ago. As transitions go, the rejuvenated stage may be thought of as the long tail of alkalic volcanism that begins in late-shield time and persists through the postshield (+rejuvenated-stage) era. Because successive Hawaiian volcanoes erupt over long and overlapping spans of time, there is a wide range in the age of volcanism along the island chain, even though the age of Hawaiian shields is progressively younger to the southeast. For example, almost every island from Ni‘ihau to Hawai‘i had an eruption in the time between 0.3 and 0.4 million years ago, even though only the Island of Hawai‘i had active volcanoes in their shield stage during that time. Once they have formed, Hawaiian volcanoes become subject to a spectrum of processes of degradation. Primary among these are subaerial erosion, landslides, and subsidence. The islands, especially those that grow high above sea level, experience mean annual precipitation that locally exceeds 9 m, leading to rapid erosion that can carve deep canyons in less than 1 million years. Hawaiian volcanoes have also been modified by giant landslides. Seventeen discrete slides that formed in the past 5 m.y. have been identified around the main Hawaiian Islands, and fully 70 are known along the Hawaiian Ridge between Midway Islands and the Island of Hawai‘i. These giant landslides displace large amounts of seawater to generate catastrophic giant waves (megatsunami). The geologic evidence for megatsunami in the Hawaiian Islands includes chaotic coral and lava-clast breccia preserved as high as 155 m above sea level on Lāna‘i and Moloka‘i. Large Hawaiian volcanoes can persist as islands through the rapid subsidence by building upward rapidly enough. But in the long run, subsidence, coupled with surface erosion, erases any volcanic remnant above sea level in about 15 m.y. One consequence of subsidence, in concert with eustatic changes in sea level, is the drowning of coral reefs that drape the submarine flanks of the actively subsiding volcanoes. At least six reefs northwest of the Island of Hawai‘i form a stairstep configuration, the oldest being deepest.

Hawaii↗

Tectonic evolution and Cretaceous gold metallogenesis of southwestern Alaska

Cretaceous gold metallogenesis in southwestern Alaska comprises three distinct episodes related to the accretionary evolution of northwestern North America. The oldest mineralizing event is characterized by 112 Ma Cu-Au-Bi-Te porphyry-type(?) veining in the zoned Bonanza and adjacent plutons that intruded rocks of the Nyac terrane. Tectonic reconstructions and limited geological and geochemical data suggest that Cu-Au mineralization in the Nyac district may be related to terminal subduction during accretion of the Togiak-Koyukuk arc. The subsequent 100 to 89 Ma metallogenic event is a product of subduction-related magmatism immediately following accretion of the Peninsular-Alexander-Wrangellia superterrane and includes formation of the giant Pebble porphyry Cu-Au-Mo deposit. Pebble underwent a complex history of highly oxidized magmatism that is isotopically linked to enriched lithosphere or metasomatized mantle sources. Pebble and other porphyryrelated plutons were emplaced as a consequence of changes in plate motion and onset of dextral transpression along the continental margin to the southeast of their present-day locations. The final 75 to 65 Ma metallogenic event is regionally the most extensive. It followed a ~15-m.y.-long magmatic lull and is related to enigmatic subduction-related magmatism in the western part of the Alaska Range and the Kuskokwim basin. This event resulted in the formation of porphyry, reduced pluton-related, orogenic, and possible epithermal Au deposits. In the better-studied relatively low-relief areas of the Kuskokwim basin, many of the mineralized systems are spatially associated with ilmenite-series monzonite to quartz monzonite composite plutons that have isotopic signatures consistent with derivation from crustally contaminated mantle sources. These pluton-hosted Au deposits comprise low-sulfide stockworks, sheeted veins, and/or breccias in the cupolas of moderately differentiated intrusions that contain high Au; anomalous As, Bi, Sb, and/or Te; and have variable, but not uneconomic, Cu concentrations. In addition, Au, Sb, and/or Hg deposits hosted by competent NE-trending granite porphyry dike-sill complexes or along faults in flysch are widespread in the Kuskokwim basin. These deposits, including the giant Donlin Creek Au deposit, formed at shallow crustal levels and are classified here as epizonal orogenic Au deposits and related Hg and Sb lodes. Mesozonal orogenic Au deposits were formed along the margins of the uplifting Willow Creek pluton that is now exposed along the southern side of the Talkeetna Mountains batholith. Recent discoveries in the rugged Alaska Range comprise less well-studied porphyry Au-Cu (Whistler, Island Mountain), epithermal(?) Au (Terra), and reduced pluton-hosted Au deposits (Estelle). The cause of the broad latest Cretaceous magmatism associated with the final metallogenic event is enigmatic, but indicates widespread high heat flow possibly related to flat-slab subduction, lithospheric mantle delamination, or escape tectonics. Plutonism postdated regional folding and coincided with periods of movement along regional faults and formation of the NE-trending structural fabric. Magmatism resulted from initiation of transpressional faulting in more landward positions during oblique subduction of the Kula plate. Crustal contamination of mantle melts, either near their source or during ascent through the thick flysch of the Kuskokwim basin, produced low oxidation state magmas and controlled much of the metallogeny, particularly for those deposits that display similarities to both reduced porphyry Au-Cu deposits and, to a lesser degree, reduced intrusion-related Au systems. High heat flow also induced crustal melting and metamorphic devolatilization to form orogenic Au deposits. A transition to extension at ~60 Ma recorded elsewhere in Alaska temporally corresponds to termination of gold deposit formation in southwestern Alaska. The multiple periods of gold metallogenesis in southwestern Alaska offer a wide variety of targets for exploration within discrete parts of the region.

Society of Economic Geologists Special Publication↗

I. Thermal evolution of Ganymede and implications for surface features. II. Magnetohydrodynamic constraints on deep zonal flow in the giant planets. III. A fast finite-element algorithm for two-dimensional photoclinometry

The work is divided into three independent papers: PAPER I: Thermal evolution models are presented for Ganymede, assuming a mostly differentiated initial state of a water ocean overlying a rock layer. The only heat sources are assumed to be primordial heat (provided by accretion) and the long-lived radiogenic heat sources in the rock component. As Ganymede cools, the ocean thins, and two ice layers develop, one above composed of ice I, and the other below composed of high-pressure polymorphs of ice. Subsolidus convection proceeds separately in each ice layer, its transport of heat calculated using a simple parameterized convection scheme and the most recent data on ice rheology. The model requires that the average entropy of the deep ice layer exceed that of the ice I layer. If the residual ocean separating these layers becomes thin enough, then a Rayleigh-Taylor-like ("diapiric") instability may ensue, driven by the greater entropy of the deeper ice and merging the two ice mantles into a single convective layer. This instability is not predicted by linear analysis but occurs for plausible finite amplitude perturbations associated with large Rayleigh number convection. The resulting warm ice diapirs may lead to a dramatic "heat pulse" at the surface and to fracturing of the lithosphere, and may be directly or indirectly responsible for resurfacing and grooved terrain formation on Ganymede. The timing of this event depends rather sensitively on poorly known rheological parameters but could be consistent with chronologies deduced from estimated cratering rates. Irrespective of the occurrence or importance of the heat pulse, we find that lithospheric fracturing requires rapid stress loading (on a timescale ≾ 10 4 ) years). Such a timescale can be realized by warm ice diapirism, but not directly by gradual global expansion. In the absence of any quantitative and self-consistent model for the resurfacing of Ganymede by liquid water, we favor resurfacing by warm ice flows,which we demonstrate to be physically possible, a plausible consequence of our models, compatible with existing observations, and a hypothesis testable by Galileo. We discuss core formation as an alternative driver for resurfacing, and conclude that it is less attractive. We also consider anew the puzzle of why Callisto differs so greatly from Ganymede, offering several possible explanations. The models presented do not provide a compelling explanation for all aspects of Ganymedean geological evolution, since we have identified several potential problems, most notably the apparently extended period of grooved terrain formation (several hundred million years), which is difficult to reconcile with the heat pulse phenomenon. PAPER II: The observed zonal flows of the giant planets will, if they penetrate below the visible atmosphere, interact significantly with the planetary magnetic field outside the metalized core. The appropriate measure of this interaction is the Chandrasekhar number Q = ( H 2 )/(4πρνα 2 λ) (where H = radial component of the magnetic field, ν = eddy viscosity, λ = magnetic diffusivity, α -1 = lengthscale on which λ varies); at depths where Q ≳ 1 the velocity will be forced to oscillate on a small lengthscale or decay to zero. We estimate the conductivity due to semiconduction in H 2 (Jupiter, Saturn) and ionization in H 2 O (Uranus, Neptune) as a function of depth; the value λ ≃ 10 10 cm 2 s -1 needed for Q = 1 is readily obtained well outside the metallic core (where λ ≃ 10 2 cm 2 s -1 ). These assertions are quantified by a simple model of the equatorial zonal jet in which the flow is assumed uniform on cylinders concentric with the spin axis, and the viscous and magnetic torques on each cylinder are balanced. We solve this "Taylor constraint" simultaneously with the dynamo equation to obtain the velocity and magnetic field in the equatorial plane. With this model we reproduce the widely differing jet widths of Jupiter and Saturn (though not the flow at very high or low latitudes) using ν = 2500 cm 2 s -1 , consistent with the requirement that viscous dissipation not exceed the specific luminosity. A model Uranian jet consistent with the limited Voyager data can also be constructed, with appropriately smaller ν, but only if one assumes a two-layer interior. We tentatively predict a wide Neptunian jet. For Saturn (but not Jupiter or Uranus) the model has a large magnetic Reynolds number where Q = 1 and hence exhibits substantial axisymmetrization of the field in the equatorial plane . This effect may or may not persist at higher latitudes. The one-dimensional model presented is only a first step. Variation of the velocity and magnetic field parallel to the spin axis must be modeled in order to answer several important questions, including: 1) What is the behavior of flows at high latitudes, whose Taylor cylinders are interrupted by the region with Q ≳ 1? 2) To what extent is differential rotation in the envelope responsible for the spin-axisymmetry of Saturn's magnetic field? PAPER III: It is shown that the problem of two-dimensional photoclinometry (PC) -- the reconstruction of a surface z ( x , y ) from a brightness image B( x , y ) -- may be formulated in a natural way in terms of finite elements. The resulting system of equations is underdetermined as a consequence of the lack of boundary conditions for z , but a unique solution may be chosen by minimizing a function S expressing the "roughness" of the surface. An efficient PC algorithm based on this formulation is presented, requiring ~ 10.66 (four-byte) memory locations and ~10 4 floating multiplications/additions per pixel, and incorporating: 1) Minimization of the roughness by the penalty method, which yields the smallest set of equations. 2) Iterative solution of the nonlinear equations by Newton's method. 3) Solution of the linearized equations by an inner iterative cycle of successive over-relaxation, which takes advantage of the extreme sparseness of the system. 4) Multigridding, in which the solutions to the smaller problems obtained by reducing the resolution are used recursively to greatly speed convergence at the higher resolutions, and 5) A rapid noniterative initial estimate of z obtained by exploiting the special symmetry of the equations obtained in the first linearization. The algorithm is extensively demonstrated on 200 by 200 pixel synthetic "images" generated from digital topographic data for northern Utah over a range of phase angles. Rms error in the solution is ~ 22 m, out of ~ 660 m total relief. The error is dominated by "stripes" with the same azimuth as the light source, resulting from use of the roughness criterion in lieu of boundary conditions; the rms error along profiles parallel to the stripes is only ~ 2-8 m, depending on the phase angle. Satisfactory solutions are obtained even in the presence of quantization error, noise, and moderate blur in the image. Applications of the PC algorithm to both remote sensing and photomicrography are sketched; a photoclinometric map of a low-relief Precambrian era fossil is presented as an example of the latter. Prospects for dealing with photometrically inhomogeneous surfaces, and an extension of the method to the analysis of side-looking radar data ("radarclinometry") are also discussed.

Thesis↗

Plate-tectonic evolution of the western U.S.A.

Changing interactions of lithospheric plates provide the framework for this review of the 3100 m.y. geological history of some 3 million km2 of mountains, deserts, plateaux and plains. The Precambrian to Neogene development of the western U.S.A. is outlined in terms of plate collisions, subduction events and deformation of lithospheric slabs, with some interpretations based on SE Asia and other regions of complex tectonics.-R.A.H.

Episodes↗

GPS constraints on the kinematics of continental deformation

Recent GPS observations from the western United States, New Zealand, central Greece, and Japan indicate that present-day continental deformation is typically focused in narrow deforming zones whose extent is much smaller than the intervening largely inactive regions. However, these narrow zones are heterogeneously distributed, reflecting the inherent heterogeneity of continental lithospheric strength and internal buoyancy. Plate driving and resisting forces stress plate boundary zones and plate interiors and drive deformation. These forces change continuously and discontinuously, leading to continental deformation that typically evolves and migrates with time. Magmatic and tectonic processes alter lithospheric rheology and internal buoyancy and also contribute to the time-varying character of continental deformation.

International Geology Review↗

Earth's magnetic field complex: U.S. National activities during the Decade of Geopotential Field Research

The US geomagnetism community is supported by NASA, NOAA, USGS, NSF, DOD, and US universities. During the Decade of Geopotential Field Research, inaugurated in 1999 with the launch of the Danish satellite Ørsted on a US rocket, the US community has been involved in satellite mission development and analysis, instrument development, model development, and in the discovery and understanding of new processes with satellite magnetic signatures. The ESA Swarm mission has been a primary focus of the US community, with three US scientists on Swarm's Mission Advisory Group. Swarm will measure, for the first time, the E-W gradient of the magnetic field. One of us (T. Sabaka) is involved with the development of a Comprehensive inversion scheme as part of the SMART consortium. This effort is an outgrowth of the Comprehensive Model [1]. Swarm will also provide valuable observations for ionospheric specification and forecast. The geomagnetism group at NOAA (S. Maus, P. Alken and C. Manoj) has developed algorithms to estimate the strength of the eastward electric field (EEF). As the driver of the equatorial plasma fountain, the EEF is an important space weather parameter. ESA is considering the implementation of the EEF as a dedicated inversion chain in the Level-2 Facility. In 2006, NASA launched a minisatellite magnetometer constellation mission (ST-5) to test technologies and software. The ST-5 constellation featured the first along-track gradient measurements. NASA has also initiated efforts to study geomagnetism mission concepts after Swarm. One of the ideas under consideration is the systematic measurement of radial field gradients. Instrument development, and geomagnetic observatories, are also an integral part of the US effort. The past decade has seen significant advances in the development of a self-calibrating vector helium magnetometer, and in the automation of the US observatory network. Working in coordination with Intermagnet, the USGS Geomagnetism Program has made operational 1-second data acquisition at 13 of its magnetic observatories. The Program is also developing a realtime 1-minute and 1-hour Dst service. Within the past decade, US scientists have been leaders in the development of models that describe the global geomagnetic environment, including comprehensive models (the CM series), maps of the lithospheric field from satellite (MF-series), near surface maps of the lithospheric field (WDMAM-series), models of the thickness of the magnetic crust, the IGRF and World Magnetic Model series, ionospheric models such as the EEJM1, JVDM1, and the IRI, and data assimilation-based models (MoSST-series) that predict the future state of the geomagneic field.

Conference Paper↗

An overview of the global variability in radiated energy and apparent stress

a global study of radiated seismic energies ER and apparent stresses τ a reveals systematic patterns. earthquakes with the highest apparent stress occur in regions of intense deformation and rupture strong lithosphere. in oceanic settings, these are strike-slip earthquakes (τ a up to 27 Mpa) occurring intraplate or at evolving ends of transform faults. at subduction zones and intracontinental settings, these are strike slip earthquakes with τ a up to 7 Mpa. normal-fault earthquakes exhibit a more complex pattern. higher τ a ’s (up to5Mpa) are found for intraslab events at depths from 35 to 70 km that occur near zones of intense deformation such as a sharp slab bend or the juncture of colliding slabs. lower τ a ’s (< 1 Mpa) are found for normal-fault earthquakes at the outer rise and outer trench wall or deep in flat warm slabs. the lowest average τ a (0.3 Mpa) is found for thrust-fault earthquakes at subduction zones. the variation of average apparent stress with tectonics suggests a relationship with lithospheric strength and fault maturity. Mature faults, such as plate boundaries that have experienced large cumulative slip, appear to have low strength and tend to yield earthquakes with low apparent stresses. immature faults, in contrast, are stronger and yield high apparent stresses because either they are the result of fresh-rock fracture or at least their cumulative fault slip is quite small. these results have implications of use to the seismic engineering community because ER and its magnitude counterpart Me are reliable indicators of the potential for damaging ground motion.

Book chapter↗

Importance of understanding bottom-up control when characterizing geothermal systems

Methods designed to identify favorable areas for geothermal resources have traditionally been focused on near-surface information, namely data that can be compiled into a 2D map. However, these methods fail to account for the third dimension: depth. As a result, they do not incorporate deep crustal and mantle features like heat sources. Geophysical methods with multi-scale capabilities, such as magnetotellurics (MT), provide tools to image deeper structures and bottom-up control on location of near surface hydrothermal systems in 3D. This study briefly demonstrates the advantage of understanding bottom-up control on hydrothermal systems to aid assessment and characterization. A regional 3D electrical resistivity model of the Great Basin is developed from MT data that images the near surface down to the mantle. From the 3D model, electrical conductance (depth integrated electrical conductivity) maps are created at logical depth intervals to identify anomalies. The conductance maps image discrete zones of high-conductance between 15-20 km depth indicative of fluid collection at the brittle-ductile transition; high-conductance zones around Moho depths (30-50 km) suggestive of partial melt; high-conductance zones in the upper mantle indicative of higher temperature and larger melt fraction, and low-conductance zones indicative of lithospheric material descending in the mantle. One anomalous zone of low-conductance in the mantle is under the north central Nevada, suggesting vertical mantle flow transport of heat to the crust related to sinking lithospheric material.

Conference Paper↗

Ice sheet load cycling and fluid underpressures in the Eastern Michigan Basin, Ontario, Canada

Strong fluid underpressures have been detected in Paleozoic strata in the eastern Michigan Basin, with hydraulic heads reaching ~400&thinsp;m below land surface (~4&thinsp;MPa underpressure) and ~200&thinsp;m below sea level in strata where unusually low permeabilities (~10 &minus;20 &ndash;10 &minus;23 &thinsp;m 2 ) were measured in situ. Multiple glaciations, including three with as much as 3&thinsp;km of ice cover at the site in the last 120&thinsp;ka, suggest a causal link with the underpressures. We examined this possibility using a one-dimensional groundwater flow model incorporating mechanical loading from both ice weight and lithospheric flexure. Because hydrologic and mechanical changes during glaciation are not well characterized and subsurface properties are imperfectly known, the model was used inversely to estimate flexural loads and loosely constrained permeabilities by matching observed pressures. Acceptable matches were obtained for a surprisingly wide range of scenarios with permeabilities close to measured values and plausible flexural loads. Matches were not obtained when too many parameters were preselected, or when permeabilities were constrained to be significantly larger than measured values. In successful model runs groundwater expulsion under glacial-mechanical loads caused the underpressuring, and flexural loads were important if aquifer and sub-glacial pressures were significantly elevated during glaciation. Simulated fluid pressures in the low-permeability strata fluctuated by 30&ndash;40&thinsp;MPa during glacial cycles but resulted in advective transport of only tens of meters or less. Although other mechanisms cannot be ruled out, we conclude that glacial-mechanical forcing of a water-saturated system can explain the observed underpressures.

Ontario↗

Eocene and Miocene extension, meteoric fluid infiltration, and core complex formation in the Great Basin (Raft River Mountains, Utah)

Metamorphic core complexes (MCCs) in the North American Cordillera reflect the effects of lithospheric extension and contribute to crustal adjustments both during and after a protracted subduction history along the Pacific plate margin. While the Miocene-to-recent history of most MCCs in the Great Basin, including the Raft River-Albion-Grouse Creek MCC, is well documented, early Cenozoic tectonic fabrics are commonly severely overprinted. We present stable isotope, geochronological ( 40 Ar/ 39 Ar), and microstructural data from the Raft River detachment shear zone. Hydrogen isotope ratios of syntectonic white mica ( δ 2 H ms ) from mylonitic quartzite within the shear zone are very low (−90‰ to −154‰, Vienna SMOW) and result from multiphase synkinematic interaction with surface-derived fluids. 40 Ar/ 39 Ar geochronology reveals Eocene (re)crystallization of white mica with δ 2 H ms ≥ −154‰ in quartzite mylonite of the western segment of the detachment system. These δ 2 H ms values are distinctively lower than in localities farther east ( δ 2 H ms ≥ −125‰), where 40 Ar/ 39 Ar geochronological data indicate Miocene (18–15 Ma) extensional shearing and mylonitic fabric formation. These data indicate that very low δ 2 H surface-derived fluids penetrated the brittle-ductile transition as early as the mid-Eocene during a first phase of exhumation along a detachment rooted to the east. In the eastern part of the core complex, prominent top-to-the-east ductile shearing, mid-Miocene 40 Ar/ 39 Ar ages, and higher δ 2 H values of recrystallized white mica, indicate Miocene structural and isotopic overprinting of Eocene fabrics.

Utah↗

Seismic velocity structure of the crust and shallow mantle of the Central and Eastern United States by seismic surface wave imaging

Seismic surface waves from the Transportable Array of EarthScope's USArray are used to estimate phase velocity structure of 18 to 125 s Rayleigh waves, then inverted to obtain three-dimensional crust and upper mantle structure of the Central and Eastern United States (CEUS) down to &sim;200 km. The obtained lithosphere structure confirms previously imaged CEUS features, e.g., the low seismic-velocity signature of the Cambrian Reelfoot Rift and the very low velocity at >150 km depth below an Eocene volcanic center in northwestern Virginia. New features include high-velocity mantle stretching from the Archean Superior Craton well into the Proterozoic terranes and deep low-velocity zones in central Texas (associated with the late Cretaceous Travis and Uvalde volcanic fields) and beneath the South Georgia Rift (which contains Jurassic basalts). Hot spot tracks may be associated with several imaged low-velocity zones, particularly those close to the former rifted Laurentia margin.

Geophysical Research Letters↗

Post-rift magmatic evolution of the eastern North American “passive-aggressive” margin

Understanding the evolution of passive margins requires knowledge of temporal and chemical constraints on magmatism following the transition from supercontinent to rifting, to post-rifting evolution. The Eastern North American Margin (ENAM) is an ideal study location as several magmatic pulses occurred in the 200 My following rifting. In particular, the Virginia-West Virginia region of the ENAM has experienced two postrift magmatic pulses at ∼152 Ma and 47 Ma, and thus provides a unique opportunity to study the long-term magmatic evolution of passive margins. Here we present a comprehensive set of geochemical data that includes new 40 Ar/ 39 Ar ages, major and trace-element compositions, and analysis of radiogenic isotopes to further constrain their magmatic history. The Late Jurassic volcanics are bimodal, from basanites to phonolites, while the Eocene volcanics range from picrobasalt to rhyolite. Modeling suggests that the felsic volcanics from both the Late Jurassic and Eocene events are consistent with fractional crystallization. Sr-Nd-Pb systematics for the Late Jurassic event suggests HIMU and EMII components in the magma source that we interpret as upper mantle components rather than crustal interaction. Lithospheric delamination is the best hypothesis for magmatism in Virginia/West Virginia, due to tectonic instabilities that are remnant from the long-term evolution of this margin, resulting in a “passive-aggressive” margin that records multiple magmatic events long after rifting ended.

Geochemistry, Geophysics, Geosystems↗

Continental rupture and the creation of new crust in the Salton Trough rift, southern California and northern Mexico: Results from the Salton Seismic Imaging Project

A refraction and wide-angle reflection seismic profile along the axis of the Salton Trough, California and Mexico, was analyzed to constrain crustal and upper mantle seismic velocity structure during active continental rifting. From the northern Salton Sea to the southern Imperial Valley, the crust is 17-18 km thick and approximately one-dimensional. The transition at depth from Colorado River sediment to underlying crystalline rock is gradual and is not a depositional surface. The crystalline rock from ~3 to ~8 km depth is interpreted as sediment metamorphosed by high heat flow. Deeper felsic crystalline rock could be stretched pre-existing crust or higher grade metamorphosed sediment. The lower crust below ~12 km depth is interpreted to be gabbro emplaced by rift-related magmatic intrusion by underplating. Low upper-mantle velocity indicates high temperature and partial melting. Under the Coachella Valley, sediment thins to the north and the underlying crystalline rock is interpreted as granitic basement. Mafic rock does not exist at 12-18 depth as it does to the south, and a weak reflection suggests Moho at ~28 km depth. Structure in adjacent Mexico has slower mid-crustal velocity and rocks with mantle velocity must be much deeper than in the Imperial Valley. Slower velocity and thicker crust in the Coachella and Mexicali valleys define the rift zone between them to be >100 km wide in the direction of plate motion. North American lithosphere in the central Salton Trough has been rifted apart and is being replaced by new crust created by magmatism, sedimentation, and metamorphism.

Journal of Geophysical Research↗

Refining the formation and early evolution of the Eastern North American Margin: New insights from multiscale magnetic anomaly analyses

To investigate the oceanic lithosphere formation and early seafloor spreading history of the North Atlantic Ocean, we examine multiscale magnetic anomaly data from the Jurassic/Early Cretaceous age Eastern North American Margin (ENAM) between 31 and 40°N. We integrate newly acquired sea surface magnetic anomaly and seismic reflection data with publicly available aeromagnetic and composite magnetic anomaly grids, satellite-derived gravity anomaly, and satellite-derived and shipboard bathymetry data. We evaluate these data sets to (1) refine magnetic anomaly correlations throughout the ENAM and assign updated ages and chron numbers to M0–M25 and eight pre-M25 anomalies; (2) identify five correlatable magnetic anomalies between the East Coast Magnetic Anomaly (ECMA) and Blake Spur Magnetic Anomaly (BSMA), which may document the earliest Atlantic seafloor spreading or synrift magmatism; (3) suggest preexisting margin structure and rifting segmentation may have influenced the seafloor spreading regimes in the Atlantic Jurassic Quiet Zone (JQZ); (4) suggest that, if the BSMA source is oceanic crust, the BSMA may be M series magnetic anomaly M42 (~168.5 Ma); (5) examine the along and across margin variation in seafloor spreading rates and spreading center orientations from the BSMA to M25, suggesting asymmetric crustal accretion accommodated the straightening of the ridge from the bend in the ECMA to the more linear M25; and (6) observe anomalously high-amplitude magnetic anomalies near the Hudson Fan, which may be related to a short-lived propagating rift segment that could have helped accommodate the crustal alignment during the early Atlantic opening.

Journal of Geophysical Research B: Solid Earth↗

Temporal stress changes caused by earthquakes: A review

Earthquakes can change the stress field in the Earth’s lithosphere as they relieve and redistribute stress. Earthquake-induced stress changes have been observed as temporal rotations of the principal stress axes following major earthquakes in a variety of tectonic settings. The stress changes due to the 2011 Mw9.0 Tohoku-Oki, Japan, earthquake were particularly well documented. Earthquake stress rotations can inform our understanding of earthquake physics, most notably addressing the long-standing problem of whether the Earth’s crust at plate boundaries is “strong” or “weak.” Many of the observed stress rotations, including that due to the Tohoku-Oki earthquake, indicate near-complete stress drop in the mainshock. This implies low background differential stress, on the order of earthquake stress drop, supporting the weak crust model. Earthquake stress rotations can also be used to address other important geophysical questions, such as the level of crustal stress heterogeneity and the mechanisms of postseismic stress reloading. The quantitative interpretation of stress rotations is evolving from those based on simple analytical methods to those based on more sophisticated numerical modeling that can capture the spatial-temporal complexity of the earthquake stress changes.

Journal of Geophysical Research B: Solid Earth↗

Deformation of the Pacific/North America plate boundary at Queen Charlotte Fault: The possible role of rheology

The Pacific/North America (PA/NA) plate boundary between Vancouver Island and Alaska is similar to the PA/NA boundary in California in its kinematic history and the rate and azimuth of current relative motion, yet their deformation styles are distinct. The California plate boundary shows a broad zone of parallel strike slip and thrust faults and folds, whereas the 49‐mm/yr PA/NA relative plate motion in Canada and Alaska is centered on a single, narrow, continuous ~900‐km‐long fault, the Queen Charlotte Fault (QCF). Using gravity analysis, we propose that this plate boundary is centered on the continent/ocean boundary (COB), an unusual location for continental transform faults because plate boundaries typically localize within the continental lithosphere, which is weaker. Because the COB is a boundary between materials of contrasting elastic properties, once a fault is established there, it will probably remain stable. We propose that deformation progressively shifted to the COB in the wake of Yakutat terrane's northward motion along the margin. Minor convergence across the plate boundary is probably accommodated by fault reactivation on Pacific crust and by an eastward dipping QCF. Underthrusting of Pacific slab under Haida Gwaii occurs at convergence angles >14°–15° and may have been responsible for the emergence of the archipelago. The calculated slab entry dip (5°–8°) suggests that the slab probably does not extend into the asthenosphere. The PA/NA plate boundary at the QCF can serve as a structurally simple site to investigate the impact of rheology and composition on crustal deformation and the initiation of slab underthrusting.

Journal of Geophysical Research B: Solid Earth↗

Contrasting volcanism in Hawaiʻi and the Galápagos

The archipelagos of Hawai&lsquo;i and the Gal&aacute;pagos originated at mantle hotspots, yet the volcanoes that make up the island chains differ in most respects. Some of the most important differences include the dynamics of magma supply, characteristics of magma storage and transport, morphology, and compositional and structural evolution. Of particular significance in the Gal&aacute;pagos is the lack of well-developed rift zones, which may be related to higher rates of pre-eruptive inflation compared to Hawai&lsquo;i, and the absence of widespread flank instability&mdash;a common feature of Hawai&lsquo;i's volcanoes. The close proximity of the Gal&aacute;pagos to a mid-ocean-ridge system may account for many of the differences between Hawaiian and Gal&aacute;pagos volcanoes. The Gal&aacute;pagos archipelago is built on young, thin oceanic crust, which might allow for contemporaneous growth of numerous volcanoes, and its volcanoes are fed by a mix of plume and asthenospheric melt sources. Hawaiian volcanoes, in contrast, grew in the middle of the Pacific Plate on older, thicker crust, where localized changes in mantle and lithosphere structure and composition did not exert dominant control over volcano evolution.

Book chapter↗

Basins in ARC-continental collisions

Arc-continent collisions occur commonly in the plate-tectonic cycle and result in rapidly formed and rapidly collapsing orogens, often spanning just 5-15 My. Growth of continental masses through arc-continent collision is widely thought to be a major process governing the structural and geochemical evolution of the continental crust over geologic time. Collisions of intra-oceanic arcs with passive continental margins (a situation in which the arc, on the upper plate, faces the continent) involve a substantially different geometry than collisions of intra-oceanic arcs with active continental margins (a situation requiring more than one convergence zone and in which the arc, on the lower plate, backs into the continent), with variable preservation potential for basins in each case. Substantial differences also occur between trench and forearc evolution in tectonically erosive versus tectonically accreting margins, both before and after collision. We examine the evolution of trenches, trench-slope basins, forearc basins, intra-arc basins, and backarc basins during arc-continent collision. The preservation potential of trench-slope basins is low; in collision they are rapidly uplifted and eroded, and at erosive margins they are progressively destroyed by subduction erosion. Post-collisional preservation of trench sediment and trench-slope basins is biased toward margins that were tectonically accreting for a substantial length of time before collision. Forearc basins in erosive margins are usually floored by strong lithosphere and may survive collision with a passive margin, sometimes continuing sedimentation throughout collision and orogeny. The low flexural rigidity of intra-arc basins makes them deep and, if preserved, potentially long records of arc and collisional tectonism. Backarc basins, in contrast, are typically subducted and their sediment either lost or preserved only as fragments in melange sequences. A substantial proportion of the sediment derived from collisional orogenesis ends up in the foreland basin that forms as a result of collision, and may be preserved largely undeformed. Compared to continent-continent collisional foreland basins, arc-continent collisional foreland basins are short-lived and may undergo partial inversion after collision as a new, active continental margin forms outboard of the collision zone and the orogen whose load forms the basin collapses in extension.

Book chapter↗