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 487 records · Page 27Linked to original sources

Radar and photoclinometric studies of wrinkle ridges on Mars

Earth-based radar altimetry and image derived photoclinometric profiles were analyzed to examine both the long- and short-wavelength topography associated with wrinkle ridges on Mars. Photoclinometrically derived elevation data across wrinkle ridges were evaluated to determine the sensitivity of profiles to two empirical photoclinometric parameters, the horizontal digital number (HDN) and the scattered light value (SLV). The photoclinometric profiles are extremely sensitive to small variations in HDN. The sense of slope of a profile can be completely reversed over a range in HDN of as little as ±1. Comparably small variations in the SLV have relatively minor effects on the photoclinometrically derived elevations. The existence of elevation offsets from one side of the ridge to the other, reported in previous photoclinometric studies of martian wrinkles, were not confirmed through photoclinometry. In addition, no evidence of elevation offsets were found in Earth-based radar altimetry profiles across wrinkle ridges. In order to more accurately model wrinkle ridge topography, we controlled photoclinometrically derived elevations with long-wavelength topography obtained from the radar altimetry. The results of this study do not support kinematic models for the origin of planetary wrinkle ridges that involve deeply rooted thrust faults which separate crustal blocks at different elevations. A kinematic model involving buckling of shallow crustal layers into concentric folds that close, leading to the development of thrust faults, is consistent with wrinkle ridge morphology and terrestrial analogs. Recent geophysical studies of terrestrial analogs and the influence of shallow subsurface structures, particularly buried craters, on the localization of many wrinkle ridges on Mars suggest that thrust faults associated with the ridges are confined to the ridged plains material and do not extend into the lithosphere.

Journal of Geophysical Research E: Planets↗

Soufrière Hills eruption, Montserrat, 1995 - 1997: volcanic earthquake locations and fault plane solutions

A total of 9242 seismic events, recorded since the start of the eruption on Montserrat in July 1995, have been uniformly relocated with station travel-time corrections. Early seismicity was generally diffuse under southern Montserrat, and mostly restricted to depths less than 7 km. However, a NE-SW alignment of epicentres beneath the NE flank of the volcano emerged in one swarm of volcano-tectonic earthquakes (VTs) and later nests of VT hypocentres developed beneath the volcano and at a separated location, under St. George's Hill. The overall spatial distribution of hypocentres suggests a minimum depth of about 5 km for any substantial magma body. Activity associated with the opening of a conduit to the surface became increasingly shallow, with foci concentrated below the crater and, after dome building started in Fall 1995, VTs diminished and repetitive swarms of ‘hybrid’ seismic events became predominant. By late-1996, as magma effusion rates escalated, most seismic events were originating within a volume about 2 km diameter which extended up to the surface from only about 3 km depth - the diminution of shear failure earthquakes suggests the pathway for magma discharge had become effectively unconstricted. Individual and composite fault plane solutions have been determined for a few larger earthquakes. We postulate that localised extensional stress conditions near the linear VT activity, due to interaction with stresses in the overriding lithospheric plate, may encourage normal fault growth and promote sector weaknesses in the volcano.

Geophysical Research Letters↗

Closing the gap between regional and global travel time tomography

Recent global travel time tomography studies by Zhou [1996] and van der Hilst et al . [1997] have been performed with cell parameterizations of the order of those frequently used in regional tomography studies (i.e., with cell sizes of 1°–2°). These new global models constitute a considerable improvement over previous results that were obtained with rather coarse parameterizations (5° cells). The inferred structures are, however, of larger scale than is usually obtained in regional models, and it is not clear where and if individual cells are actually resolved. This study aims at resolving lateral heterogeneity on scales as small as 0.6° in the upper mantle and 1.2°–3° in the lower mantle. This allows for the adequate mapping of expected small-scale structures induced by, for example, lithosphere subduction, deep mantle upwellings, and mid-ocean ridges. There are three major contributions that allow for this advancement. First, we employ an irregular grid of nonoverlapping cells adapted to the heterogeneous sampling of the Earth's mantle by seismic waves [ Spakman and Bijwaard , 1998]. Second, we exploit the global data set of Engdahl et al . [1998], which is a reprocessed version of the global data set of the International Seismological Centre. Their reprocessing included hypocenter redetermination and phase reidentification. Finally, we combine all data used ( P , pP , and pwP phases) into nearly 5 million ray bundles with a limited spatial extent such that averaging over large mantle volumes is prevented while the signal-to-noise ratio is improved. In the approximate solution of the huge inverse problem we obtain a variance reduction of 57.1%. Synthetic sensitivity tests indicate horizontal resolution on the scale of the smallest cells (0.6° or 1.2°) in the shallow parts of subduction zones decreasing to approximately 2°–3° resolution in well-sampled regions in the lower mantle. Vertical resolution can be worse (up to several hundreds of kilometers) in subduction zones with rays predominantly pointing along dip. Important features of the solution are as follows: 100–200 km thick high-velocity slabs beneath all major subduction zones, sometimes flattening in the transition zone and sometimes directly penetrating into the lower mantle; large high-velocity anomalies in the lower mantle that have been attributed to subduction of the Tethys ocean and the Farallon plate; and low-velocity anomalies continuing across the 660 km discontinuity to hotspots at the surface under Iceland, east Africa, the Canary Islands, Yellowstone, and the Society Islands. Our findings corroborate that the 660 km boundary may resist but not prevent (present day) large-scale mass transfer from upper to lower mantle or vice versa. This observation confirms the results of previous, global mantle studies that employed coarser parameterizations.

Journal of Geophysical Research B: Solid Earth↗

Morphology and stratal geometry of the Antarctic continental shelf: Insights from models

Reconstruction of past ice-sheet fluctuations from the stratigraphy of glaciated continental shelves requires understanding of the relationships among the stratal geometry, glacial and marine sedimentary processes, and ice dynamics. We investigate the formation of the morphology and the broad stratal geometry of topsets on the Antarctic continental shelf with numerical models. Our models assume that the stratal geometry and morphology are principally the results of time-integrated effects of glacial erosion and sedimentation related to the location of the seaward edge of the grounded ice. The location of the grounding line varies with time almost randomly across the shelf. With these simple assumptions, the models can successfully mimic salient features of the morphology and the stratal geometry. The models suggest that the current shelf has gradually evolved to its present geometry by many glacial advances and retreats of the grounding line to different locations across the shelf. The locations of the grounding line do not appear to be linearly correlated with either fluctuations in the 5 l s O record (which presumably represents changes in the global ice volume) or with the global sea-level curve, suggesting that either a more complex relationship exists or local effects dominate. The models suggest that erosion of preglacial sediments is confined to the inner shelf, and erosion decreases and deposition increases toward the shelf edge. Some of the deposited glacial sediments must be derived from continental erosion. The sediments probably undergo extensive transport and reworking obliterating much of the evidence for their original depositional environment. The flexural rigidity and the tectonic subsidence of the underlying lithosphere modify the bathymetry of the shelf, but probably have little effect on the stratal geometry. Our models provide several guidelines for the interpretation of unconformities, the nature of preserved topset deposits, and the significance of progradation versus aggradation of shelf sediments.

Book chapter↗

Probing the Earth's strength: Can we measure small stress at high pressure?

Simulating the conditions and processes that occur in the Earth's deep interior has been a major goal of experimental geophysics since the 1920s. In particular, pioneers such as P. W. Bridgman, David Griggs, Hugh Heard, Mervyn Paterson, William Brace, and their colleagues sought to establish the basic relations between differential stresses and rock and mineral deformation under pressure-temperature conditions that favor such inelastic processes as brittle fracture, frictional sliding, plastic deformation, and dynamic metamorphism. This work has been applied to diverse problems in the Earth sciences, including the physical nature of earth quake sources, the forces that drive and resist plate movement, deformation under partial-melting conditions at the tops of mantle plumes, the flexure of the lithosphere, post-glacial rebound, and the origin of the preferred orientation of minerals and the resultant seismic anisotropy and metamorphic textures.

Eos, Transactions, American Geophysical Union↗

Results of a reconnaissance microearthquake survey of Bucaramanga, Colombia

Six University of Wisconsin portable, continuously‐recording seismographs were operated for 3½ days in late 1976 in the Eastern Cordillera of Colombia in a 200‐km‐diameter array around Bucaramanga, where there are also three permanent stations of the Instituto Geofísico de Los Andes Colombianos. Twenty‐seven microearthquakes were recorded. Most can be well located. Only one event, located along the eastern boundary of the Eastern Cordillera, is at crustal depths. At least two events, beneath the western half of the Eastern Cordillera, are 110‐115 km in depth. Twenty earthquakes are associated with the Bucaramanga nest (6.83°N, 73.12°W), at a depth of 158 ± 4 km. For these well‐recorded events, relative arrival times between stations vary by up to ± 0.3 sec, and the epicentral area is 4‐5 km in diameter. The complex character of these seismograms, compared with those from 110‐115 km depth events outside the nest, suggests either complexity in the rupture process or, more likely, multi‐pathing which originates near the source region. For the upper mantle and crust above the nest, the average values for Q are about 450 for shear waves and probably over 1,000 for compressional waves, suggesting that the rays travel through lithospheric material for their entire path. If the two events at intermediate depths outside the nest and the nest lie at the same horizon within the same slab, they define a slab that strikes N‐S and dips to the east at 33°.

Geophysical Research Letters↗

Pleistocene high-silica rhyolites of the Coso volcanic field, Inyo County, California

The high-silica rhyolite domes and lava flows of the bimodal Pleistocene part of the Coso volcanic field provide an example of the early stages of evolution of a silicic magmatic system of substantial size and longevity. Major and trace element compositions are consistent with derivation from somewhat less silicic parental material by liquid state differentiation processes in compositionally and thermally zoned magmatic systems. Seven chemically homogeneous eruptive groups can be distinguished on the basis of trace element and K/Ar data. The oldest two groups are volumetrically minor and geochemically distinct from the younger groups, all five of which appear to have evolved from the same magmatic system. Erupted volume-time relations suggest that small amounts of magma were bled from the top of a silicic reservoir at a nearly constant long-term rate over the last 0.24Ma. The interval of repose between eruptions appears to be proportional to the volume of the preceding eruptive group. This relationship suggests that eruptions take place when some parameter which increases at a constant rate reaches a critical value; this parameter may be extensional strain accumulated in roof rocks. Extension of the lithosphere favors intrusion of basalt into the crust, attendant partial melting, and maintenance of a long-lived silicic magmatic system. The Coso silicic system may contain a few hundred cubic kilometers of magma. The Coso magmatic system may eventually have the potential for producing voluminous pyroclastic eruptions if the safety valve provided by rapid crustal extension becomes inadequate to 1) defuse the system through episodic removal of volatile-rich magma from its top and 2) prohibit migration of the reservoir to a shallow crustal level.

California↗

Extension in the Rio Grande rift

A positive gravity anomaly along the axis of the Rio Grande rift reflects a volume of anomalous mass added at the base of the crust and intruded into the crust. Part of this volume can be associated with vertical uplift of the crust. The remainder of this anomalous volume, plus the volume of surficial graben fill, can be associated with horizontal crustal extension. The volume of crustal uplift in the Rio Grande rift is unknown, but this term can be eliminated by means of an independent equation provided by assumption of generalized isostatic equilibrium. The volume and mass equations combined provide a solution for extension of the crust in terms of the following parameters: total anomalous mass deficiency in the mantle lithosphere, total anomalous mass excess in the crust and its density contrast, total anomalous mass deficiency of surficial graben fill and its density contrast, and the volume of material eroded from the uplift. Using standard density estimates and masses determined by equivalent-source modeling of gravity profiles, I obtained 1-km extension at 37°N (Colorado-New Mexico border), 13-km extension at 35°N (Albuquerque, New Mexico), and 24-km extension at 33°N in southern New Mexico. These estimates are of necessity very poorly constrained and may be as much as ±60% in error. Taking results at face value, the best fitting Euler pole occurs at about 41°N, in north central Colorado. The calculated opening angle is 1.32°. For an approximately 30-m.y. duration of the extensional system (late Oligocene to the present), calculated average angular velocity is 7.8×10 −4 rad/m.y.; calculated average spreading half-rate at 33°N (for example) is a stately 0.04 cm/yr. The ratio of extension to uplift increases southward.

Journal of Geophysical Research Solid Earth↗

Wandering terranes in southern Alaska: The Aleutia Microplate and implications for the Bering Sea

Paleomagnetic and geological data suggest that much of southern Alaska is a collage of tectonostratigraphic terranes which originated in Mesozoic time at paleolatitudes far south of their present position. The time of ‘docking’ of the terranes against cratonic Alaska is critical to defining their amalgamated size and extent during their northward motion as well as their role in the evolution of the Bering Sea. One of the largest of the tectonostratigraphic terranes, the Peninsular terrane of south central and southwestern Alaska, extends offshore along the outer Bering Sea continental margin (Beringia). Paleomagnetic data suggest that this terrane has moved northward through all of Cenozoic time, but geologic data imply that the terrane had accreted to Alaska by the end of the Mesozoic. In early Cenozoic time the eastern part of the Aleutian arc appears to have been superimposed on the Peninsular terrane, and postulated northward Cenozoic motion of the terrane would therefore have required northward motion of the arc. Two accretion models, based on docking times for terranes in Alaska, are proposed, and they illustrate that large areas of the abyssal Bering Sea, the Alaska Peninsula, the Aleutian arc, and the Beringian continental margin may be part of a superterrane or microplate called Aleutia (microplate as defined by Beck et al. (1980), i.e., a microplate is a displaced segment of lithosphere that has crustal roots, whereas a superterrane is an amalgamation of terranes which may or may not be rootless). Model A implies that the Aleutian arc developed in situ on the southern edge of Aleutia after the microplate had docked. In model B, the final docking time of the Peninsular terrane is late Cenozoic, which implies that the Aleutia microplate encompasses a mammoth area that includes parts of southern Alaska, the Alaska Peninsula, the southern Beringian margin, the abyssal Bering Sea (Kula plate), and the Aleutian arc. If model A is correct, the docking time of the Peninsular terrane is late Mesozoic or earliest Tertiary. The Aleutia microplate in this model is made up solely of the abyssal Bering Sea (Kula plate), which presumably docked at the same time or slightly after the Peninsular terrane accreted against Alaska. If model B is correct, that is, if the Aleutia collided with nuclear Alaska during the Cenozoic, then a late Cenozoic suture zone, the vestige of a large open sea that must have closed between Aleutia and Alaska, must exist in south central and southwest Alaska. Either evidence for Cenozoic closure and suturing has been obliterated in Alaska or the inferences of Cenozoic terrane motion derived from paleomagnetic data are suspect.

Alaska↗

A dislocation model of strain accumulation and release at a subduction zone

Strain accumulation and release at a subduction zone are attributed to stick slip on the main thrust zone and steady aseismic slip on the remainder of the plate interface. This process can be described as a superposition of steady state subduction and a repetitive cycle of slip on the main thrust zone, consisting of steady normal slip at the plate convergence rate plus occasional thrust events that recover the accumulated normal slip. Because steady state subduction does not contribute to the deformation at the free surface, deformation observed there is completely equivalent to that produced by the slip cycle alone. The response to that slip is simply the response of a particular earth model to embedded dislocations. For a purely elastic earth model, the deformation cycle consists of a coseismic offset followed by a linear‐in‐time recovery to the initial value during the interval between earthquakes. For an elastic‐viscoelastic earth model (elastic lithosphere over a viscoelastic asthenosphere), the postearthquake recovery is not linear in time. Records of local uplift as a function of time indicate that the long‐term postseismic recovery is approximately linear, suggesting that elastic earth models are adequate to describe the deformation cycle. However, the deformation predicted for a simple elastic half‐space earth model does not reproduce the deformation observed along the subduction zones in Japan at all well if stick slip is restricted to the main thrust zone. As recognized earlier by Shimazaki, Seno, and Kato, the uplift profiles could be explained if stick slip were postulated to extend along the plate interface beyond the main thrust zone to a depth of perhaps 100 km, but independent evidence suggests that stick slip at such depths is unlikely.

Journal of Geophysical Research B: Solid Earth↗

Geology of the Mount St. Helens area: Record of discontinuous volcanic and plutonic activity in the Cascade Arc of southern Washington

The Quaternary edifice of Mount St. Helens volcano was built upon a deeply eroded terrane of gently folded and altered volcanic and plutonic rocks that represent the core of the Tertiary Cascade magmatic arc. These rocks constitute an east dipping homoclinal sequence, several kilometers thick, of subaerially erupted mafic to silicic flows and volcaniclastic strata; K‐Ar ages from this section range from about 28 to 23 Ma (late Oligocene and earliest Miocene), which corresponds to an apparent lull in Cascade volcanism to the north of Mount Rainier. Volcanism was essentially continuous during this period of time, and neither a well‐defined base nor top of the section is exposed within the mapped area. Basalt and basaltic andesite dominate the lower part of the mid‐Tertiary section, whereas andesitic and dacitic rocks comprise most of the upper part. This section was intruded by numerous mafic to silicic dikes, sills, and irregular plutonic bodies, most no more than a few million years younger than their host rocks, and subjected to pervasive burial metamorphism and widespread hydrothermal alteration. Large areas of hornfelsed rock surrounding even relatively small intrusions indicate that the proportion of plutonic rock becomes significantly greater at shallow depth beneath the existing erosion surface. A large granitic pluton intruded the mid‐Tertiary section north of Spirit Lake at about 21 Ma. The Earl porphyry copper deposit occurs within the pluton but appears too young (17 Ma) to be genetically related to it. In contrast to the rather continuous and voluminous Oligocene to early Miocene activity, volcanism since then in the Mount St. Helens area has been localized and Volumetrically minor. Products of three younger eruptive periods have been recognized: a sequence of 15 m.y. old pyroxene andesite flows resting unconformably on mid‐Tertiary strata south of Mount St. Helens, widespread shallow dikes and sills of pyroxene andesite between 10 and 8 m.y. old, and compositionally diverse rocks erupted during the past 3 m.y. The Quaternary lavas are more potassic than the Tertiary lavas and typically contain phenocrysts of hornblende and biotite, which are absent from the older rocks. A number of Tertiary structures define a broad NNE trending zone that may reflect a deep‐seated lithospheric flaw that has controlled the locus of Cascade magmatism in southern Washington for the past 25 m.y. Mount St. Helens lies within this zone at the intersection of the NNW striking St. Helens seismic zone (SHZ) and an ENE trending alignment of Pleistocene silicic plug‐domes. No surface breakage has been detected along the SHZ, which is apparently very young. The linear zone of silicic vents is probably controlled by a fault that has been interpreted from seismic records to occur directly beneath the volcano. This zone parallels the directions of regional maximum horizontal compressive stress and North America/Juan de Fuca plate convergence. Mount St. Helens is an example of a low‐volume tectonically controlled magmatic system in an early stage of development.

Washington↗

Tectonic history of the Syria Planum province of Mars

We attribute most of the development of extensive fractures in the Tharsis region to discrete tectonic provinces within the region, rather than to Tharsis as a single entity. One of these provinces is in Syria Planum. Faults and collapse structures in the Syria Planum tectonic province on Mars are grouped into 13 sets based on relative age, areal distribution, and morphology. According to superposition and fault crosscutting relations and crater counts we designate six distinct episodes of tectonic activity in the following chronologic order: stage I is an early to late Noachian deformation forming mostly east-west structures (fault set IA); some large volcanoes also formed. Faults were produced possibly by flexural uplift. Also, arcuate, north trending grabens (set IB) indicate that faulting is transitional to the next stage. Stage II is late Noachian to early Hesperian radial faulting centered in Syria Planum, possibly due to isostatic uplift in late Noachian (set IIA) to early Hesperian (set IIB) time. Stage III is early to late Hesperian faulting tangential to Syria Planum that was related to local centers of uplift (sets III1-III3) on the periphery of Syria Planum. Stage IV is a late Hesperian graben formation that was circumferential to Syria Planum (set IVA), caused either by collapse associated with eruption of magma or by flexure of the lithosphere due to volcanic loading. In association with volcanism, minor faulting occurred, producing an oval pattern of faults in southwest Syria Planum (set IVBl) and a radial pattern south of the planum (set IVB2) that apparently rejuvenated buried stage II faults. Stage V is a late Hesperian to early Amazonian development of grabens and troughs of Noctis Labyrinthus and western Valles Marineris (set V) that was probably instigated by local uplift; exposure of groundwater or ground ice zones may have produced further collapse and trough enlargement. Stage VI is early Amazonian northwest trending faulting in Noctis Fossae (set VIl), perhaps due to Tharsis Monies-centered tectonism, and north northwest normal faulting along the eastern side of the Claritas rise (set VI2) that was due to tectonic subsidence. The duration of tectonic activity in the Syria Planum province was perhaps 2–3 b.y. Photoclinometric topographic profiles across 132 grabens and fault scarps show that Syria Planum grabens have widths (average of 2.5 km, and most range from 1 to 6 km) similar to lunar grabens, but the Martian grabens have slightly higher side walls (average about 132 m) and gentler wall slopes (average of 9° and range of 2°–25°) than lunar grabens (93 m high and 18° slopes). Scarp degradation on Mars has progressed through quakes, impact shaking, and dry slope processes; the lower slopes may be due to Mars’ higher gravity. Estimates of the amount of extension for individual grabens range from 20 to 350 m; most estimates of the thickness of the faulted layer range from 0.5 to 4.5 km (average is 1.5 km). This thickness range corresponds closely to the 0.8-to 3.6-km range in depth for pits, troughs, and canyons in Noctis Labyrinthus and along the walls of Valles Marineris. We propose that the predominant 1-to 1.5-km values obtained for both the thickness of the faulted layer and the depths of the pits, troughs, and theater heads of the canyons reflect the initial depth to the water table in this region, as governed by the depth to the base of ground ice. Maximum depths for these features may indicate lowered groundwater table depths and the base of ejecta material.

Journal of Geophysical Research Solid Earth↗

Stress origins and earthquake potentials in Cascadia

Focal mechanism solutions for shallow earthquakes throughout the Cascadia plate system indicate that the primary regional stress is northerly compression, even though the Juan de Fuca plate generally is thought to be subducting N50°E. This compressional stress is pervasive throughout the Gorda-Juan de Fuca-Explorer plate system and much of the adjoining section of North American plate. Modeling, using a discrete element code, shows that this stress primarily is due to the Pacific plate being driven into the Gorda block and Juan de Fuca plate (at the Mendocino and Blanco fracture zones). This plate collision causes compression of the offshore plate system northward through the 45°W trending Vancouver Island and leads to transfer of northerly compression into the overriding plate. The northerly compression in the Cascadia plate system has caused large, crustal earthquakes in the regions of Vancouver and northern California and may be capable of causing large earthquakes also in and offshore of Washington and Oregon. Several independent lines of evidence confirm that the Cascadia subduction interface is coupled strongly and that subduction primarily must be accompanied by significant earthquakes, rather than occurring aseismically. For numerous earthquakes within the subducted plate, focal mechanisms indicate extensional stresses that trend downdip. These earthquakes reflect plate extension that is due to sinking of the more deeply subducted plate and to plate motion being resisted at the shallow subduction interface. The slab pull force of the sinking plate also can lead to earthquakes at the shallow subduction interface and the potential exists for M 7.5–8.0 subduction earthquakes to occur at segments of the subduction boundary at Washington and Vancouver and less frequently so at segments of the subduction boundary at Oregon. The in-plate driving forces of the offshore Cascadia plate system do not appear to contribute significantly to regional stresses. Stresses and earthquakes at the Cascadia province primarily result from the action of the Pacific plate on the Gorda-Juan de Fuca-Explorer plate system, a transfer of this action into the overriding North American plate, and the independent sinking and slab pull force of the subducted oceanic lithosphere.

Journal of Geophysical Research Solid Earth↗

Seismic reflection images of the crust of the northern part of the Chugach terrane, Alaska: Results of a survey for the Trans-Alaska Crustal Transect (TACT)

Deep crustal seismic reflection data show strong reflections from the middle and lower crust of the convergent continental margin near the eastern end of the Aleutian trench. These data were collected across the Border Ranges fault system, a major suture zone that separates the Peninsular and Chugach tectonostratigraphic terranes. The Chugach terrane consists of metamorphosed, strongly deformed flysch, subduction melange as well as blueschist- and greenschist-facies rocks. Even so, highly reflective rocks lie at great depth within or beneath this terrane. The shallow part of the seismic section (0–5 s) reveals only a few reflections, which undulate and interweave, indicating that some of them originate outside of the plane of the seismic section. Other undulating events could reveal early Cenozoic granitic plutons or antiformal stacks of rock imbricated along thrust faults. In contrast to this shallow reflection pattern, three reflection bands contain most of the events on the middle part of the seismic section (5–12 s or 12–34 km). We interpret divergent reflections within the uppermost band (5–5.5 s) as gently north dipping thrust faults, and subparallel reflections truncated at the boundaries of the upper and middle bands may image a metamorphic foliation or thrust faults. If the truncated reflections represent thrust faults, then the reflection geometry suggests that subparallel roof and floor thrust faults bound imbricated rocks, forming duplex structures. The upper and middle bands are separated by an area on the seismic section that reveals few reflections. This area correlates approximately in depth with a high velocity (7.5 km/s) layer evident in refraction models. The poorly reflective, high velocity unit might be volcanic rocks like those exposed in the structurally lowest part of the Chugach terrane. If so, then the bottom of the poorly reflective rocks indicates the base of the Chugach terrane. Rocks that cause the upper and middle reflection bands and the intervening, poorly reflective volcanic rocks may form a midcrustal shear zone that is about 10 km thick. The deepest band of subparallel reflections (30–35 km) correlates closely in depth with the top of the Wadati-Benioff zone associated with the underthrusting Pacific plate. The reflections probably stem from the decollement that separates lithospheric plates within the subduction zone. This reflection band increases abruptly in thickness, possibly indicating local underplating of subducted sediment. No events from the oceanic Moho have been recognized. Neither the Border Ranges fault system, the suture between the Chugach and Peninsular terranes, nor deep crustal layers of the Peninsular terrane are evident in these seismic reflection data.

Journal of Geophysical Research Solid Earth↗

State of stress and modern deformation of the northern Basin and Range Province

Constraints on the current stress regime of the actively extending northern Basin and Range province are provided by deformation data (focal mechanisms and fault slip studies), hydraulic fracturing in situ stress measurements, borehole elongation (“breakouts”) analyses, and alignment of young volcanic vents. The integrated data indicate significant variations both in principal stress orientations and magnitudes. An approximately E-W least principal stress direction appears to characterize both the eastern and western margins of the Basin and Range province, whereas in the active interior parts of the province extension occurs in response to a least principal stress oriented NW to N60°W. The contrast in stress orientations between the province boundaries and in the interior suggests that along the margins the least principal stress direction may be locally controlled by the generally northerly trending profound lithospheric discontinuities associated with these margins. Active deformation along the southeastern and western province margins is characterized by a combination of strike-slip and normal faulting. Focal mechanisms along northeastern province margin (Wasatch front) and in central Nevada indicate a combination of normal and oblique-normal faulting. Temporal, regional, and depth-dependent variations in the relative magnitudes of the vertical and maximum horizontal stresses can explain much of the observed variations in deformation styles. However, some depth variation in faulting style inferred from focal mechanisms may be apparent and simply a function of the attitude of fault planes being reactivated. Evidence for significant temporal variation (or multiple cycles of variation) in relative stress magnitude comes from the Sierran front-Basin and Range boundary region where recent earthquakes are predominantly strike slip, whereas the profound relative vertical relief across the Sierra frontal fault zone in the last 9–10 m.y. implies a normal faulting stress regime. Using the best data on stress orientation, relative stress magnitudes are constrained from slip vectors of major earthquakes and young fault displacements. Analysis of well-constrained slip vectors in the Owens Valley, California, area indicate that large temporal variations in the magnitude of the approximately N-S oriented maximum horizontal stress are required to explain dominantly dip-slip and strike-slip offsets on subparallel faults. Similar faulting relations are observed throughout much of the boundary zone between the Basin and Range-Sierra Nevada (including the Walker Lane belt). Along the eastern province margin in the Wasatch front area in Utah, available data suggest that the maximum and minimum horizontal stresses may be approximately equal at depths of <4–5 km. Earthquake focal mechanisms in this area suggest more variability in relative magnitude of the two horizontal stresses with depth. Furthermore, superimposed sets of young fault striae along a segment of the Wasatch fault also indicate temporal variations of relative stress magnitudes. Sources of regional and temporal variations in the stress field may be linked to variable shear tractions applied to the base of the brittle crust related to intrusion, thermally induced flow, and the influence of the San Andreas plate boundary. Although difficult to date accurately, the fault slip data suggest that the temporal variations in relative magnitudes stress may occur on the time scale of both a single major earthquake cycle (1000–5000 years) and multiple earthquake cycles (10,000+ years).

Journal of Geophysical Research Solid Earth↗

Igneous history of the Koyukuk terrane, western Alaska: Constraints on the origin, evolution, and ultimate collision of an accreted island arc terrane

The Koyukuk terrane of western Alaska consists of volcanic, volcaniclastic, and plutonic rocks which range from Late Paleozoic to Early Cretaceous in age. The terrane crops out in a U-shaped belt which is roughly paralleled by outer belts of ultramafic rocks, oceanic plate basalts and cherts, and retrograded blueschist facies rocks of continental protolith. These rocks have been interpreted as components of a volcanic arc terrane that collided with the North American continental margin in Early Cretaceous time. The Koyukuk terrane consists of four time-stratigraphic units: (1) pre-Middle Jurassic basalts, (2) Middle and Late Jurassic granitic rocks, (3) lower Lower Cretaceous volcanic rocks, and (4) upper Lower Cretaceous volcanic rocks. Limited chemical data from the basalts of unit 1 indicate that they were erupted in a nonarc tectonic environment, possibly in an oceanic island or back arc setting. Units 2, 3, and 4 have the characteristics of subduction-related volcanic rocks (i.e., depleted Nb and Ta and enriched alkaline elements, relative to the light rare earth elements). Unit 3 contains tholeiitic, calc-alkaline, and alkaline rocks with chondrite-normalized rare earth element patterns that range from flat (La N /Yb N = 1) to highly light rare earth element enriched (La N /Yb N > 15). The highly alkaline or shoshonitic lavas were erupted toward the end of unit 3 time (Valanginian) during the final stages of arc-continent collision. These alkaline lavas could have been derived by very small degrees of partial melting of a similar source to that of the earlier arc lavas. Unit 4 lavas are also alkaline or shoshonitic, but their incompatible element composition indicates that they were derived from a different source than that of the earlier arc lavas. These late alkaline lavas are chemically similar to crosscutting mid-Cretaceous plutons whose isotopic compositions (Arth et al., this issue ( a )) suggest derivation by partial melting of distinctly older subcontinental lithosphere. We speculate that the parental magmas of unit 4 lavas may also have been derived by partial melting of this subcontinental mantle which was underthrust beneath the Koyukuk arc terrane during the final stage of arc-continent collision.

Journal of Geophysical Research Solid Earth↗

The nature of the crust in the Yukon-Koyukuk province as inferred from the chemical and isotopic composition of five Late Cretaceous to Early Tertiary volcanic fields in western Alaska

Late Cretaceous and early Tertiary volcanic and plutonic rocks in western Alaska comprise a vast magmatic province extending from the Alaska Range north to the Arctic Circle, south to Bristol Bay, and west to the Bering Sea Shelf. The chemical and isotopic composition of five of these Late Cretaceous to early Tertiary volcanic fields in the north central part of this province were studied to determine if Paleozoic or older continental crust underlies the Yukon-Koyukuk province. Three of the fields, the Blackburn Hills, Yukon River, and Kanuti, occur within the Yukon-Koyukuk province and two, the Sischu and Nowitna, overlie bordering Precambrian and Paleozoic metamorphic terranes to the southeast. High initial 87 Sr/ 86 Sr of 0.7075–0.7079 and moderate initial 143 Nd/ 144 Nd of 0.51244–0.51247 of rhyolite, dacite, and high-silica andesite of the Sischu volcanic field indicate that the magmas have interacted with the underlying Paleozoic or older continental crust. The relatively limited variation of isotopic (initial 87 Sr/ 86 Sr = 0.7044–0.7051; initial 143 Nd/ 144 Nd = 0.51256–0.51257) and elemental compositions of andesites from the Nowitna field can be accounted for by assimilation of small amounts of Paleozoic or older continental crust during crystal fractionation of andesite parent magmas at crustal levels. The Blackburn Hills field, which consists of medium-K basalt, andesite, and rhyolite intruded by a small granitic pluton, has a large range in initial 87 Sr/ 86 Sr and initial 143 Nd/ 144 Nd that plot in the field for 60 Ma mantle, from near mid-ocean ridge basalts to near “bulk-earth” compositions (initial 87 Sr/ 86 Sr = 0.7033–0.7052; initial 143 Nd/ 144 Nd = 0.51253–0.51290). Andesites and basalts from the Blackburn Hills are divided into two group on the basis of rare earth element (REE) and isotopic composition. Isotopic variation in the more primitive group 1 is best explained by assimilation of the lower crust of the Jurassic to Early Cretaceous Koyukuk terrane by mantle-derived basalts during crystal fractionation, though part of the isotopic variation may be due to metasomatism of an oceanic island basalt type mantle source by fluids derived from subducted sediments. Group 2 andesites from the Blackburn Hills have lower heavy REE abundances and more enriched isotopic compositions. These group 2 andesites and dacites from the Kanuti field, which have ( 87 Sr/ 86 Sr) i = 0.7043–0.7048 and ( 143 Nd/ 144 Nd) i = 0.51248–0.51267, appear to have formed by partial melting of the lower crust of the Koyukuk terrane. The Yukon River field consists of basalt, andesite, dacite, and rhyolite having ( 87 Sr/ 86 Sr) i = 0.7037–0.7051 and ( 143 Nd/ 144 Nd) i = 0.51266–0.51280; its isotopic composition does not require the presence of Paleozoic or older continental crust under the volcanic field and may have formed by interaction between mantle-derived melts and the oceanic Angayucham/Tozitna or island arc Koyukuk terrane. Most of the intrusive rocks and rhyolite domes from the Blackburn Hills volcanic field have ( 87 Sr/ 86 Sr) i = 0.7038–0.7041 and dacites from the Kanuti volcanic field have ( 87 Sr/ 86 Sr) i = 0.7043–0.7048. Thus little or no old continental crust was involved in the genesis of the Late Cretaceous and early Tertiary rocks and therefore probably does not extend beneath this part of the Yukon-Koyukuk province. However, the ultimate source of the small volumes of enriched shoshonitic andesite ( 87 Sr/ 86 Sr = 0.7075, 143 Nd/ 144 Nd = 0.5125) erupted at 118 Ma in the Yukon-Koyukuk province may be continental lithosphere, which may have been thrust under this part of the Yukon-Koyukuk province during arc-continent collision in the Early Cretaceous.

Journal of Geophysical Research Solid Earth↗

Lithology and evolution of the crust-mantle boundary region in the southwestern Basin and Range Province

Mantle and crustal xenoliths from volcanic rocks in the southwestern Basin and Range province and Colorado Plateau Transition Zone reveal histories of episodic magmatism and deformation that have profoundly influenced the crustal structure of this region. Seismic transects in this area show a strongly reflective Moho of generally low relief, which, in the area of modern transects, consists of a thin zone (<2 km thick) of short reflectors. The upper mantle is transparent and has a P n of 7.8–8.0 km/s similar to much of the western United States. A lower crustal zone, 2–13 km thick, has variable internal reflectivity and a relatively low velocity of 6.6–6.8 km/s. Upper mantle peridotite xenoliths show both ductile and brittle deformational features and have structures and compositions affected by magmatic intrusion; intrusions form complex dike systems and extensive zones of grain boundary infiltration in peridotite xenoliths. Whereas melt infiltration preceded and followed ductile deformation, brittle deformation, represented by closely spaced joint systems and faults, followed ductile deformation and is related to the youngest magmatic episodes. These structural characteristics and high uppermost mantle temperature (∼1000°C) may combine to explain the relatively low P n . Alternating layers of ductily deformed and undeformed peridotites, with or without igneous intrusions, may contribute to the reflectivity of the Moho. Lower crustal xenoliths are dominantly igneous-textured pyroxenites and mafic to intermediate gabbros identical to the dikes in peridotite xenoliths. The crustal xenoliths also commonly are jointed, and in addition many show partial melting and have abundant cavities that probably were filled with CO 2 -rich fluids. These rocks are interpreted as products of underplated magmas that were fed through the mantle dike systems and may represent the lowest crustal unit identified in the seismic records. The mafic compositions and high densities of the crustal xenoliths indicate that the low velocity of the lower crust may be caused in part by fracture systems, partial melts, and high temperatures. Garnet granulite xenoliths from a locality with no mantle peridotite xenoliths probably represent crust of the region before late Miocene extension. Felsic granulite xenoliths from two localities have velocities like those of the two lower crustal units identified seismically and could be present in the modern crust as unequilibrated remnants of old crust. The preferred model for the evolution of the lower lithosphere is one in which extension affects the upper mantle as well as the crust and is overlapped in time by multiple magmatic episodes. The earliest magmatic events preceded extension, and later events accompanied and followed extension.

Arizona, California, Nevada↗