Gary S. Fuis

Community Fault Model (CFM) for Southern California

We present a new three-dimensional model of the major fault systems in southern California. The model describes the San Andreas fault and associated strike- slip fault systems in the eastern California shear zone and Peninsular Ranges, as well as active blind-thrust and reverse faults in the Los Angeles basin and Transverse Ranges. The model consists of triangulated surface representations (t-surfs) of more than 140 active faults that are defined based on surfaces traces, seismicity, seismic reflection profiles, wells, and geologic cross sections and models. The majority of earthquakes, and more than 95% of the regional seismic moment release, occur along faults represented in the model. This suggests that the model describes a comprehen- sive set of major earthquake sources in the region. The model serves the Southern California Earthquake Center (SCEC) as a unified resource for physics-based fault systems modeling, strong ground-motion prediction, and probabilistic seismic hazards assessment.

Mapping the megathrust beneath the northern Gulf of Alaska using wide‐angle seismic data

In the northern Gulf of Alaska and Prince William Sound, we have used wide-angle seismic reflection/refraction profiles, earthquake studies, and laboratory measurements of physical properties to determine the geometry of the Prince William and Yakutat terranes, the Aleutian megathrust, and the subducting Pacific plate. In this complex region, the Yakutat terrane is underthrust beneath the Prince William terrane, and both terranes are interpreted to be underlain by the Pacific plate. Wide-angle seismic reflection/refraction profiles recorded along five seismic lines are used to unravel this terrane geometry. Modeled velocities in the upper crust of the Prince William terrane (to 18 km depth) agree closely with laboratory velocity measurements of Orca Group phyllites and quartzofeldspathic graywackes (the chief components of the Prince William terrane) to hydrostatic pressures as high as 600 MPa (6 kbar). A landward dipping reflector at depths of 16–24 km is interpreted as the base of the Prince William terrane. This reflector corresponds to the top of the Wadati-Benioff zone seismicity and is interpreted as the megathrust. Immediately beneath the megathrust is a 4-km-thick 6.9-km/s refractor, which we infer to be the source of a prominent magnetic anomaly and which is interpreted by us and previous workers to be gabbro in Eocene age oceanic crust of the underthrust Yakutat terrane. Wide-angle seismic data, magnetic anomaly data, and tectonic reconstructions indicate that the Yakutat terrane has been underthrust beneath the Prince William terrane for at least a few hundred kilometers. Wide-angle seismic data are consistent with a 9° to 10° landward dip of the subducting Pacific plate beneath the outer shelf and slope, distinctly different from the inferred average 3° to 4° dip of the overlying 6.9-km/s refractor and the Wadati-Benioff seismic zone beneath the inner shelf. Our preferred interpretation of the geophysical data is that one composite plate, composed of the Pacific plate of a fairly uniform thickness and the Yakutat plate of varying thickness, is subducting beneath southern Alaska.

Trans-Alaska Crustal Transect and continental evolution involving subduction underplating and synchronous foreland thrusting

We investigate the crustal structure and tectonic evolution of the North American continent in Alaska, where the continent has grown through magmatism, accretion, and tectonic under-plating. In the 1980s and early 1990s, we conducted a geological and geophysical investigation, known as the Trans-Alaska Crustal Transect (TACT), along a 1350-km-long corridor from the Aleutian Trench to the Arctic coast. The most distinctive crustal structures and the deepest Moho along the transect are located near the Pacific and Arctic margins. Near the Pacific margin, we infer a stack of tectonically underplated oceanic layers interpreted as remnants of the extinct Kula (or Resurrection) plate. Continental Moho just north of this underplated stack is more than 55 km deep. Near the Arctic margin, the Brooks Range is underlain by large-scale duplex structures that overlie a tectonic wedge of North Slope crust and mantle. There, the Moho has been depressed to nearly 50 km depth. In contrast, the Moho of central Alaska is on average 32 km deep. In the Paleogene, tectonic underplating of Kula (or Resurrection) plate fragments overlapped in time with duplexing in the Brooks Range. Possible tectonic models linking these two regions include flat-slab subduction and an orogenic-float model. In the Neogene, the tectonics of the accreting Yakutat terrane have differed across a newly interpreted tear in the subducting Pacific oceanic lithosphere. East of the tear, Pacific oceanic lithosphere subducts steeply and alone beneath the Wrangell volcanoes, because the overlying Yakutat terrane has been left behind as underplated rocks beneath the rising St. Elias Range, in the coastal region. West of the tear, the Yakutat terrane and Pacific oceanic lithosphere subduct together at a gentle angle, and this thickened package inhibits volcanism.

Crustal structure and tectonics from the Los Angeles basin to the Mojave Desert, southern California

A seismic refraction and low-fold reflection survey, known as the Los Angeles Region Seismic Experiment (LARSE), was conducted along a transect (line 1) extending from Seal Beach, California, to the Mojave Desert, crossing the Los Angeles and San Gabriel Valley basins and San Gabriel Mountains. The chief result of this survey is an interpreted cross section that addresses a number of questions regarding the crustal structure and tectonics of southern California that have been debated for decades and have important implications for earthquake hazard assessment. The results (or constraints) are as follows. (1) The maximum depth of the Los Angeles basin along line 1 is 8–9 km. (2) The deep structure of the Sierra Madre fault zone in the northern San Gabriel Valley is as follows. The Duarte branch of the Sierra Madre fault zone forms a buried, 2.5-km-high, moderately north dipping buttress between the sedimentary and volcanic rocks of the San Gabriel Valley and the igneous and metamorphic rocks of the San Gabriel Mountains. (For deeper structure, see following.) (3) There are active crustal decollements in southern California. At middle-crustal depths, the Sierra Madre fault zone appears to sole into a master decollement that terminates northward at the San Andreas fault and projects southward beneath the San Gabriel Valley to the Puente Hills blind thrust fault. (4) The dip and depth extent of the San Andreas fault along line 1 dips steeply (∼83°) northward and extends to at least the Moho. (5) The subsurface lateral extent of the Pelona Schist in southern California is as follows. Along line 1, the Pelona Schist underlies much, if not all of the San Gabriel Mountains south of the San Andreas fault to middle-crustal depths. North of the San Andreas fault, it is apparently not present along the transect.

Anatomy of a metamorphic core complex: Seismic refraction/wide‐angle reflection profiling in southeastern California and western Arizona

The metamorphic core complex belt in southeastern California and western Arizona is a NW-SE trending zone of unusually large Tertiary extension and uplift. Midcrustal rocks exposed in this belt raise questions about the crustal thickness, crustal structure, and the tectonic evolution of the region. Three seismic refraction/wide-angle reflection profiles, acquired and analyzed as a part of the U.S. Geological Surveys Pacific to Arizona Crustal Experiment, were collected to address these issues. The results presented here, which focus on the Whipple and Buckskin-Rawhide mountains, yield a consistent three-dimensional image of this part of the metamorphic core complex belt. The seismic refraction/wide-angle reflection data are of excellent quality and are characterized by six principal phases that can be observed on all three profiles. These phases include refractions from the near-surface and crystalline basement, reflections from boundaries in the middle and lower crust, and reflections and refractions from the upper mantle. The final model consists of a thin veneer (<2 km) of upper plate and fractured lower plate rocks (1.5–5.5 kms−1) overlying a fairly homogeneous basement (∼6.0 km s−1) and a localized high-velocity (6.4 km s−1) body situated beneath the western Whipple Mountains. A prominent midcrustal reflection is identified beneath the Whipple and Buckskin-Rawhide mountains between 10 and 20 km depth. This reflector has an arch-like shape and is centered beneath, or just west of, the metamorphic core complex belt. This event is underlain by a weaker, approximately subhorizontal reflection at 24 km depth. Together, these two discontinuities define a lens-shaped midcrustal layer with a velocity of 6.35–6.5 km s−1. The apex of this midcrustal layer corresponds roughly to a region of major tectonic denudation and uplift (∼10 km) defined by surface geologic mapping and petrologic barometry studies. The layer thins to the northeast and is absent in the Transition Zone. The 6.35–6.5 km s−1 velocities are compatible with a diorite composition or a mixture of mafic and silicic rocks. This midcrustal layer is underlain by a higher-velocity lower crustal layer that is modeled as only 3–6 km thick beneath the metamorphic core complex belt and regions to the southwest. To the northeast, however, this layer thickens to 8–10 km as the midcrustal layer pinches out above it. The velocity of the lower crust is constrained by traveltime modeling and is 6.6±0.15 kms−1 beneath the western Transition Zone and the metamorphic core complex belt; higher velocities may be present farther to the southwest where the layer is thin. The velocity of the lower crust is too low to accommodate significant amount of mafic underplating at the base of the crust. Instead, we interpret the velocities to indicate that the lower crust is passively thinned beneath these regions without significant addition of mafic mantle-derived intrusions. The crust-mantle boundary does not dome up beneath the core complexes but remains approximately subhorizontal at a depth of 26–28 km or, in the case of the Whipple Mountains, actually deepens; a 3-km crustal root is modeled. This lack of upward doming of the Moho, together with the vertical alignment of the metamorphic core complex belt over what are believed to be extension-related structures in the middle and lower crust, suggest that there is no lateral offset of upper crustal deformation from deeper zones of extension, as one would expect if extension occurred along crust-penetrating shear zones (Wernicke, 1981; Wernicke et al., 1985). Instead, domed and inflated middle crust and thinned lower crust directly underlie the region of greatest thinning of the upper crust.

Crustal structure of accreted terranes in southern Alaska, Chugach Mountains and Copper River Basin, from seismic refraction results

Seismic refraction data were collected along a 320-km-long “Transect” line in southern Alaska, crossing the Prince William, Chugach, Peninsular, and Wrangellia terranes, and along several shorter lines within individual terranes. Velocity structure in the upper crust (less than 9-km depth) differs among the four terranes. In contrast, layers in the middle crust (9- to 25-km depth) in some cases extend across projected terrane boundaries. The following observations can be made: (1) An intermediate-velocity layer (6.4 km/s) at 9-km depth extends across the deep projection of the suture between the Chugach and Peninsular terranes, suggesting that the northern Chugach and southern Peninsular terranes are detached and rest on a deeper terrane of unknown origin. (2) The top of a gently north dipping sequence of low- and high-velocity layers (5.7–7.8 km/s), more than 10 km thick, extends from near the surface in the southern Chugach terrane to more than 20-km depth beneath the southern Peninsular terrane. This sequence, truncated by the suture between the Prince William and Chugach terranes, is interpreted to be an underplated “terrane” made up of fragments of the Kula plate and its sedimentary overburden that were accreted during subduction in the late Mesozoic and/or early Tertiary, during or between times of accretion of the Prince William and Chugach terranes. (3) A thick crustal “root”, with a laminated sequence at its top, extends from a depth of 19 km to as much as 57 km beneath the northern Peninsular and Wrangellia terranes. This root extends across the deep projection of the suture between the Peninsular and Wrangellia terranes, although resolution of this apparent crosscutting relationship is relatively poor. This root may represent tectonically or, possibly, magmatically emplaced rocks. The lower crust beneath the Prince William, Chugach, and southern Peninsular terranes includes a north dipping, 3- to 8-km-thick section of subducting oceanic crust.

Fault systems of the 1971 San Fernando and 1994 Northridge earthquakes, southern California: Relocated aftershocks and seismic images from LARSE II

We have constructed a composite image of the fault systems of the M 6.7 San Fernando (1971) and Northridge (1994), California, earthquakes, using industry reflection and oil test well data in the upper few kilometers of the crust, relocated aftershocks in the seismogenic crust, and LARSE II (Los Angeles Region Seismic Experiment, Phase II) reflection data in the middle and lower crust. In this image, the San Fernando fault system appears to consist of a decollement that extends 50 km northward at a dip of ∼25° from near the surface at the Northridge Hills fault, in the northern San Fernando Valley, to the San Andreas fault in the middle to lower crust. It follows a prominent aseismic reflective zone below and northward of the main-shock hypocenter. Interpreted upward splays off this decollement include the Mission Hills and San Gabriel faults and the two main rupture planes of the San Fernando earthquake, which appear to divide the hanging wall into shingle- or wedge-like blocks. In contrast, the fault system for the Northridge earthquake appears simple, at least east of the LARSE II transect, consisting of a fault that extends 20 km southward at a dip of ∼33° from ∼7 km depth beneath the Santa Susana Mountains, where it abuts the interpreted San Fernando decollement, to ∼20 km depth beneath the Santa Monica Mountains. It follows a weak aseismic reflective zone below and southward of the main-shock hypocenter. The middle crustal reflective zone along the interpreted San Fernando decollement appears similar to a reflective zone imaged beneath the San Gabriel Mountains along the LARSE I transect, to the east, in that it appears to connect major reverse or thrust faults in the Los Angeles region to the San Andreas fault. However, it differs in having a moderate versus a gentle dip and in containing no mid-crustal bright reflections.

Accretion and subduction tectonics in the Chugach Mountains and Copper River Basin, Alaska: initial results of the Trans-Alaska Crustal Transect

Geologic, seismic, gravity, and magnetic data from the northern Chugach Mountains and southern Copper River Basin, Alaska, indicate that the Chugach terrane (CGT) and the composite Peninsular/Wrangellia terrane (PET/WRT) are thin (< 10 km), rootless sheets bounded on the south by north-dipping thrust faults that sole into a shallow, horizontal, low-velocity zone. The CGT has been thrust at least 40 km beneath the PET/WRT along the Border Ranges fault system (BRFS). Adjacent to the BRFS, uplift and erosion of 30-40 km since Jurassic time have exposed blueschist-facies rocks in the CGT and mafic and ultramafic cumulate rocks in the PET/WRT. Four paired north-dipping layers of low and high seismic velocities extend beneath the northern CGT and southern PET/WRT and may be slices of subducted oceanic crust and upper mantle; the upper two pairs may now be joined to the continental plate. 15 references, 5 figures.

A geologic interpretation of seismic-refraction results in northeastern California

In 1981, the U.S. Geological Survey conducted a seismic-refraction experiment in northeastern California designed to study the Klamath Mountains, Cascade Range, Modoc Plateau, and Basin and Range provinces. Key profiles include 135-km-long, north-south lines in the Klamath Mountains and Modoc Plateau provinces and a 260-km-long, east-west line crossing all of the provinces. The seismic-velocity models for the Klamath and Modoc lines are comparatively homogeneous laterally but are quite different from each other. The Klamath model is finely layered from the surface to at least 14-km depth, consisting of a series of high-velocity layers (6.1–6.7 km/s), ranging in thickness from 1 to 4 km, with alternating positive and negative velocity gradients. A layer with an unreversed velocity of 7.0 km/s extends from 14 km to an unknown depth. The Modoc model, in contrast, is relatively thickly layered and has lower velocities than does the Klamath model at all depths down to 25 km. An upper layer, 4.5 km thick, of low-velocity material (2.1–4.4 km/s) overlies a basement with a considerably higher velocity (6.2 km/s). Velocity increases slowly with depth, with a small velocity step (to 6.4 km/s) at 11 km and a 7.0-km/s layer beginning at 25-km depth. Moho is probably 38–45 km deep under the Modoc Plateau, but its depth is unknown under the Klamath Mountains. A combined velocity-density model for the east-west line consists of a western part similar in configuration to the Klamath velocity model, an eastern part similar to the Modoc velocity model, and laterally changing velocity-density structure in between, in the Cascade Range. Beneath its upper layer, the velocity model for the Modoc Plateau is similar to that determined by other researchers for the adjacent Sierra Nevada. The velocity model is unlike those for rift areas, to which the Modoc Plateau has been compared by some authors. We theorize that beneath a veneer of volcanic and sedimentary rocks (the upper layer), the Modoc Plateau is underlain by a basement of granitic and metamorphic rocks that, like rocks in the Sierra Nevada, are the roots of one or more magmatic arcs. The fine layering in the Klamath seismic-velocity model is consistent with the geologic structure of the Klamath Mountains, characterized by imbricate thrusting of oceanic rock layers of various compositions and ages. Independent modeling of aeromagnetic data indicates that the base of the Trinity ultramafic sheet, the second major rock layer down in the structural sequence, corresponds to a velocity step to 6.7 km/s at 7-km depth in our model. The 6.7-km/s layer beneath the Trinity ultramafic sheet apparently corresponds to rocks of the central metamorphic belt, which are mafic schists. Rock units structurally deeper than rocks of the central metamorphic belt can be correlated with velocity layers below the 6.7-km/s layer, but with less certainty. In the model for the east-west line, the region of laterally changing velocity structure beneath the Cascade Range includes a 10-km step down to the east in the top of the 7.0-km/s layer. This region of lateral velocity change we interpret to be a fault, fold, or intrusive contact (or some combination of the three) between the stack of oceanic rock layers that underlie the Klamath Mountains and the buried roots of magmatic arcs inferred to underlie the Modoc Plateau. Magmas forming the modern Cascade Range arc apparently rise through this region.

Geophysical evidence for the evolution of the California Inner Continental Borderland as a metamorphic core complex

We use new seismic and gravity data collected during the 1994 Los Angeles Region Seismic Experiment (LARSE) to discuss the origin of the California Inner Continental Borderland (ICB) as an extended terrain possibly in a metamorphic core complex mode. The data provide detailed crustal structure of the Borderland and its transition to mainland southern California. Using tomographic inversion as well as traditional forward ray tracing to model the wide-angle seismic data, we find little or no sediments, low (≤6.6 km/s) P wave velocity extending down to the crust-mantle boundary, and a thin crust (19 to 23 km thick). Coincident multichannel seismic reflection data show a reflective lower crust under Catalina Ridge. Contrary to other parts of coastal California, we do not find evidence for an underplated fossil oceanic layer at the base of the crust. Coincident gravity data suggest an abrupt increase in crustal thickness under the shelf edge, which represents the transition to the western Transverse Ranges. On the shelf the Palos Verdes Fault merges downward into a landward dipping surface which separates “basement” from low-velocity sediments, but interpretation of this surface as a detachment fault is inconclusive. The seismic velocity structure is interpreted to represent Catalina Schist rocks extending from top to bottom of the crust. This interpretation is compatible with a model for the origin of the ICB as an autochthonous formerly hot highly extended region that was filled with the exhumed metamorphic rocks. The basin and ridge topography and the protracted volcanism probably represent continued extension as a wide rift until ∼13 m.y. ago. Subduction of the young and hot Monterey and Arguello microplates under the Continental Borderland, followed by rotation and translation of the western Transverse Ranges, may have provided the necessary thermomechanical conditions for this extension and crustal inflow.