Mid-Cretaceous thrusting in the southern Coast Belt, British Columbia and Washington, after strike-slip fault reconstruction

A major thrust system of mid-Cretaceous age is present along much of the Coast Belt of northwestern North America. Thrusting was concurrent, and spatially coincided, with emplacement of a great volume of are intrusives and minor local strike-slip faulting. In the southern Coast Belt (52 degrees to 47 degrees N), thrusting was followed by major dextral-slip faulting, which resulted in significant translational shuffling of the thrust system. In this paper, we restore the displacements on major dextral-slip faults of the southern Coast Belt and then analyze the mid-Cretaceous thrust system. Two reconstructions were made that use dextral faulting on the Yalakom fault (115 km), Castle Pass and Ross Lake faults (10 km), and Fraser fault (100 km). The reconstructions differ in the amount of dextral offset on the Straight Creek fault (160 and 100 km) and how much the NE part of the Cascades crystalline core expanded (30 km and 0 km) during Eocene extension. Reconstruction A produces the best match of lithotectonic units and thrust systems. Our synthesis shows that the southern Coast Belt thrust system was greater than or equal to 250 - 180 km wide after thrusting. The thrust system was mainly southwest vergent but had a belt of northeast vergent back thrusts on the northeast side associated with the Tyaughton-Methow basin, which may indicate large-scale tectonic wedging. Thrust faults are commonly low to moderate angle, but high angle faults also occur, especially as late stage, out-of-sequence, structures involving plutons. The amount of thrust displacement across the system is unknown but must be at least 100 km and may be many hundreds of kilometers. Most thrusting occurred from similar to 100 to similar to 80 Ma and did not migrate systematically until after similar to 90 Ma, when thrusting and magmatism shifted to the northeast for a few million years. Widespread thrusting occurred both near plutons and where there are no (or small) plutons, which strongly suggests that thrust faulting was caused by regional- to plate-scale forces such as rapid plate convergence and/or are-continent collision

Upper-plate response to ridge subduction and oceanic plateau accretion, Washington Cascades and surrounding region: Implications for plate tectonic evolution of the Pacific Northwest (USA and southwestern Canada) in the Paleogene

The interaction between subduction zones and oceanic spreading centers is a common tectonic process, and yet our understanding of how it is manifested in the geologic record is limited to a few well-constrained modern and ancient examples. In the Paleogene, at least one oceanic spreading center interacted with the northwestern margin of North America. Several lines of evidence place this triple junction near Washington (USA) and southern British Columbia (Canada) in the early to middle Eocene, and we summarize a variety of new data sets that permit us to track the plate tectonic setting and geologic evolution of this region from 65 to 40 Ma. The North Cascades segment of the voluminous Coast Mountains continental magmatic arc experienced a magmatic lull between ca. 60 and 50 Ma interpreted to reflect low-angle subduction. During this period of time, the Swauk Basin began to subside inboard of the paleo-trench in Washington, and the Siletzia oceanic plateau began to develop along the Farallon plate– Kula plate or Farallon plate– Resurrection plate spreading center. Farther east, peraluminous magmatism occurred in the Omineca belt and Idaho batholith. Accretion of Siletzia and ridge-trench interaction occurred between ca. 53 and 49 Ma, as indicated by: (1) near-trench magmatism from central Vancouver Island to northwestern Washington, (2) disruption and inversion of the Swauk Basin during a short-lived contractional event, (3) voluminous magmatism in the Kamloops-Challis belt accompanied by major E-W extension east of the North Cascades in metamorphic core complexes and supra-detachment basins and grabens, and (4) southwestward migration of magmatism across northeastern Washington. These events suggest that flat-slab subduction from ca. 60 to 52 Ma was followed by slab rollback and breakoff during accretion of Siletzia. A dramatic magmatic flare-up was associated with rollback and breakoff between ca. 49.4 and 45 Ma and included bimodal volcanism near the eastern edge of Siletzia, intrusion of granodioritic to granitic plutons in the crystalline core of the North Cascades, and extensive dike swarms in the North Cascades. Transtension during and shortly before the flare-up led to \u3e300 km of total offset on dextral strike-slip faults, formation of the Chumstick strike-slip basin, and subhorizontal ductile stretching and rapid exhumation of rocks metamorphosed to 8– 10 kbar in the North Cascades crystalline core. By ca. 45 Ma, the Farallon– Kula (or Resurrection)–North American triple junction was likely located in Oregon (USA), subduction of the Kula or Resurrection plate was established outboard of Siletzia, and strikeslip faulting was localized on the north-striking Straight Creek– Fraser River fault. Motion of this structure terminated by 35 Ma. These events culminated in the establishment of the modern Cascadia convergent margin

Stratigraphy, age, and provenance of the Eocene Chumstick basin, Washington Cascades; implications for paleogeography, regional tectonics, and development of strike-slip basins: Reply

We welcome the comment by Evans (2022) and the opportunity to further discuss our study of the Chumstick Formation. The correlation of fault-bound nonmarine sedimentary units in central and western Washington has been a topic of interest, and debate, for many years (Frizzell, 1979; Taylor et al., 1988; Gresens et al., 1981; Gresens, 1983; Evans and Johnson, 1989; Evans, 1994; Cheney and Hayman, 2009). However, many questions about the regional correlation of these units were resolved with the publication of a suite of internally consistent high-precision 206Pb/238U zircon dates from volcanic interbeds throughout the early to middle Eocene stratigraphy (Eddy et al., 2016). This data set confirmed the timing of sediment deposition of the different members within the Chumstick Formation. Donaghy et al. (2021) provides a detailed study of the Chumstick Formation, which builds on earlier research by Gresens et al. (1981, 1983), McClincy (1986), and Evans (1994) by incorporating new geochronologic information and additional clast counts, detrital zircon geochronology, and facies mapping. We interpret large parts of the Chumstick Formation to represent a spatially and temporally distinct sedimentary system between the Leavenworth and Entiat fault zones that likely formed as a pull-apart basin. Evans (2022) objects to several of the interpretations presented in Donaghy et al. (2021) regarding the relationship between different members of the Chumstick Formation and surrounding sedimentary units, the timing of strike-slip faulting, and the regional tectonic setting of these rocks. We discuss each of these points in the following sections

Provenance and Paleogeography of the 25-17 Ma Rainbow Gardens Formation: Evidence for Tectonic Activity at Ca. 19 Ma and Internal Drainage rather than Throughgoing Paleorivers on the Southwestern Colorado Plateau

The paleogeographic evolution of the Lake Mead region of southern Nevada and northwest Arizona is crucial to understanding the geologic history of the U.S. Southwest, including the evolution of the Colorado Plateau and formation of the Grand Canyon. The ca. 25–17 Ma Rainbow Gardens Formation in the Lake Mead region, the informally named, roughly coeval Jean Conglomerate, and the ca. 24–19 Ma Buck and Doe Conglomerate southeast of Lake Mead hold the only stratigraphic evidence for the Cenozoic pre-extensional geology and paleogeography of this area. Building on prior work, we present new sedimentologic and stratigraphic data, including sandstone provenance and detrital zircon data, to create a more detailed paleogeographic picture of the Lake Mead, Grand Wash Trough, and Hualapai Plateau region from 25 to 18 Ma. These data confirm that sediment was sourced primarily from Paleozoic strata exposed in surrounding Sevier and Laramide uplifts and active volcanic fields to the north. In addition, a distinctive signal of coarse sediment derived from Proterozoic crystalline basement first appeared in the southwestern corner of the basin ca. 25 Ma at the beginning of Rainbow Gardens Formation deposition and then prograded north and east ca. 19 Ma across the southern half of the basin. Regional thermochronologic data suggest that Cretaceous deposits likely blanketed the Lake Mead region by the end of Sevier thrusting. Post-Laramide northward cliff retreat off the Kingman/Mogollon uplifts left a stepped erosion surface with progressively younger strata preserved northward, on which Rainbow Gardens Formation strata were deposited. Deposition of the Rainbow Gardens Formation in general and the 19 Ma progradational pulse in particular may reflect tectonic uplift events just prior to onset of rapid extension at 17 Ma, as supported by both thermochronology and sedimentary data. Data presented here negate the California and Arizona River hypotheses for an “old” Grand Canyon and also negate models wherein the Rainbow Gardens Formation was the depocenter for a 25–18 Ma Little Colorado paleoriver flowing west through East Kaibab paleocanyons. Instead, provenance and paleocurrent data suggest local to regional sources for deposition of the Rainbow Gardens Formation atop a stripped low-relief western Colorado Plateau surface and preclude any significant input from a regional throughgoing paleoriver entering the basin from the east or northeast

Modeling Flow and Transport at Pore Scale with Obstructions

In this thesis we study mathematical and computational models for phenomena of flow and transport in porous media in the presence of changing pore scale geometries. The differential equations for the flow and transport models at Darcy scale involve the coefficients of permeability, porosity, and tortuosity which depend on the pore scale geometry. The models we propose help to understand how the presence of obstruc- tions impacts the Darcy scale models. The particular changes in pore scale geometry we consider are due to the formation of obstructions to the flow, and come from two important applications of interest, biofilm clogging and gas hydrate crystal plugging up the pores. The direct simulations or experiments of these processes at pore scale is generally unfeasible or impractical. We propose two computationally efficient mathematical and computational models to simulate the formation of the obstructions. The first method extends the phase separation model based on the Allen-Cahn equation; in our variant we add volume constraints and additional localization functions. The second method we propose is a Markov Chain Monte Carlo method inspired by the Ising model; here we use heuristics to choose the particular coefficients which guide the formation of obstructions of a particular type. After we generate independent realizations of the obstructed geometries, we solve flow and transport problems at pore scale. Next we use the technique called upscaling which carries the information to larger scale by averaging, and we are able to derive the ensemble of Darcy scale properties for a collection of generated pore scale geometries with obstructions. We show how these techniques can be used in synthetic geometries as well as in geometries obtained from imaging. In addition, we see that the permeability coefficient is not merely a function of porosity, but is rather highly dependent on the type of obstruction growing at the pore scale

Eocene dike orientations across the Washington Cascades in response to a major strike-slip faulting episode and ridge-trench interaction

The northern Cascade Mountains in Washington (USA) preserve an exceptional shallow to mid-crustal record of Eocene transtension marked by dextral strike-slip faulting, intrusion of dike swarms and plutons, rapid non-marine sedimentation, and ductile flow and rapid cooling in parts of the North Cascades crystalline core. Transtension occurred during ridge-trench interaction with the formation of a slab window, and slab rollback and break-off occurred shortly after collision of the Siletzia oceanic plateau at ca. 50 Ma. Dike swarms intruded a \u3e1250 km2 region between ca. 49.3 Ma and 44.9 Ma, and orientations of more than 1500 measured dikes coupled with geochronologic data provide important snapshots of the regional strain field. The mafic Teanaway dikes are the southernmost and most voluminous of the swarms. They strike NE (mean = 036°) and average ~15 m in thickness. To the north, rhyolitic to basaltic dikes overlap spatially with 49.3-46.5 Ma, mainly granodioritic plutons, but they typically predate the nearby plutons by ca. 500 k.y. The average orientations of five of the six dike domains range from 010° to 058°; W-NW- to NW-striking dikes characterize one domain and are found in lesser amounts in a few other domains. Overall, the mean strike for all Eocene dikes is 035°, and the average extension direction (305°-125°) is oblique to the strike (~320°) of the North Cascades orogen. Extension by diking reached ~45% in one \u3e7-km-long transect through the Teanaway swarm and ranged from ~5% to locally ~79% in shorter transects across other swarms, which corresponds to a minimum of ~12 km of extension. The dominant NE-striking dikes are compatible with the dextral motion on the N- to NW-striking (~355-320°) regional strike-slip faults. Some of the W-NW- to NW-striking dikes were arguably influenced by pre-existing faults, shear fractures, and foliations, and potentially in one swarm where both NE-and lesser W-NW-striking dikes are present, by a switch in principal stress axes induced by dike emplacement. Alternatively, the W-NW- to NW-striking dikes may reflect a younger regional strain field, as ca. 49.3-47.5 Ma U-Pb zircon ages of the NE-striking dikes are older than those of the few dated W-NW- to NW-trending dikes. In one scenario, NE-striking dikes intruded during an interval when strain mainly reflected dextral strike-slip faulting, and the younger dikes record a switch to more arc-normal extension. Diking ended as magmatism migrated into a N-S-trending belt west of the North Cascades core that marks the initiation of the ancestral Cascade arc

Farallon slab detachment and deformation of the Magdalena Shelf, southern Baja California

Subduction of the Farallon plate beneath northwestern Mexico stalled by ~12 Ma when the Pacific-Farallon spreading-ridge approached the subduction zone. Coupling between remnant slab and the overriding North American plate played an important role in the capture of the Baja California (BC) microplate by the Pacific Plate. Active-source seismic reflection and wide-angle seismic refraction profiles across southwestern BC (~24.5 degrees N) are used to image the extent of remnant slab and study its impact on the overriding plate. We infer that the hot, buoyant slab detached ~40 km landward of the fossil trench. Isostatic rebound following slab detachment uplifted the margin and exposed the Magdalena Shelf to wave-base erosion. Subsequent cooling, subsidence and transtensional opening along the shelf (starting ~8 Ma) starved the fossil trench of terrigenous sediment input. Slab detachment and the resultant rebound of the margin provide a mechanism for rapid uplift and exhumation of forearc subduction complexes

Seismic structure of the southern Gulf of California from Los Cabos block to the East Pacific Rise

Multichannel reflection and coincident wide-angle seismic data collected during the 2002 Premier Experiment, Sea of Cortez, Addressing the Development of Oblique Rifting (PESCADOR) experiment provide the most detailed seismic structure to date of the southern Gulf of California. Multichannel seismic (MCS) data were recorded with a 6-km-long streamer, 480-channel, aboard the R/V Maurice Ewing, and wide-angle data was recorded by 19 instruments spaced every similar to 12 km along the transect. The MCS and wide-angle data reveal the seismic structure across the continent-ocean transition of the rifted margin. Typical continental and oceanic crust are separated by a similar to 75-km-wide zone of extended continental crust dominated by block-faulted basement. Little lateral variation in crustal thicknesses and seismic velocities is observed in the oceanic crust, suggesting a constant rate of magmatic productivity since seafloor spreading began. Oceanic crustal thickness and mean crustal velocities suggest normal mantle temperature (1300 degrees C) and passive mantle upwelling at the early stages of seafloor spreading. The crustal thickness, width of extended continental crust, and predicted temperature conditions all indicate a narrow rift mode of extension. On the basis of upper and lower crust stretching factors, an excess of lower crust was found in the extended continental crust. Total extension along transect 5W is estimated to be similar to 35 km. Following crustal extension, new oceanic crust similar to 6.4-km-thick was formed at a rate of similar to 48 mm a(-1) to accommodate plate separation

Late Quaternary Faulting History of the Carrizal and Related Faults, La Paz Region, Baja California Sur, Mexico

The southwest margin of the Gulf of California has an array of active normal faults despite this being an oblique-divergent plate boundary with spreading centers that localized deformation along the plate boundary 2–3 million years ago. The Carrizal and Centenario faults form the western border fault of the Gulf of California marginal fault system within and south of La Paz Bay, and ∼20–30 km west of the capital city of La Paz, Baja California Sur, Mexico. Geologic and geomorphic mapping, optically stimulated luminescence (OSL) geochronology, and paleoseismic investigations onshore, compressed high-intensity radar pulse (CHIRP) profiling offshore, and analysis of uplifted marine terraces in the footwall of the offshore Carrizal fault provide some of the first numerical and geometrical constraints on late Pleistocene–Holocene faulting along the Carrizal fault. The onshore Carrizal fault has ruptured with up to ∼1–2 m of vertical displacement per event, likely producing ∼M 6.3–6.9 earthquakes, and at least two to three surface rupturing earthquakes have occurred since 22 ka. Onshore paleoseismic excavations and uplifted marine terraces on the western side of La Paz Bay both suggest offset rates of 0.1–0.2 mm/yr, with a footwall uplift rate of 0.13 mm/yr since 128 ka, and an approximately constant rate since marine oxygen-isotope stage (MIS) 11 terraces (420 ka). A CHIRP survey identified underwater fault scarps with heights ranging from 21 to 86 m on the Carrizal fault in La Paz Bay and from 3 to 5 m along the Centenario fault. The offshore Carrizal fault lies 8–10 km east of the western edge of La Paz Bay, forming a right step from the onshore Carrizal fault. The offshore Carrizal fault is the oldest fault of the fault system, and the fault likely grew in the latest Miocene to Pliocene in a complex way to the south toward the onshore Centenario and Carrizal faults. When the Alarcon spreading center started its modern rates at 2.4 Ma, the Carrizal fault likely slowed to the 0.1–0.2 mm/yr rates of the late Quaternary determined in this study