Zircon geochronology
To test models for the Laramide orogeny, we compiled > 260 Pb/U ages, including 79 new zircon and titanite dates to establish the timing and duration of magmatism, metamorphism and deformation in the Southern California Batholith with the goal of resolving conflicting models. Areal addition rates were calculated from integrating igneous ages with pluton areas determined from digitized geologic maps of Southern California (see Supplementary file). Our compilation and new data encompass all major blocks in the SCB including the San Gabriel, Pine Mountain, Alamo-Frazier Mountain, Little San Bernardino, San Bernardino, and Salinian blocks (Fig. 1B). Igneous zircon ages reveal 4 discrete pulses of Late Paleozoic to Mesozoic magmatism at 260–210 Ma, 160–140 Ma, 120–118 Ma, and 90–70 Ma (Fig. 2A). The latter pulse culminated in an arc flare-up event which peaked at 85–75 Ma and was associated with widespread, voluminous plutonism throughout the SCB (Fig. 1B, 2B). These Late Cretaceous plutons have typical Cordilleran arc magma geochemistry: most are magnesian, metaluminous to weakly peraluminous, and calc-alkalic with strong crustal affinities that reflect mixtures of Proterozoic crust with juvenile (Mesozoic) sources23,24. Although the cause of the 85–75 Ma flare-up remains unclear, magmatism persisted in the frontal arc of the SCB until ca. 65 Ma and our data show that most igneous activity had ceased by ca. 70 Ma (Fig. 2C; Supplementary file).
An analysis of metamorphic zircon and titanite shows that the batholith records high temperatures during the Late Cretaceous transition to thick-skinned, Laramide shortening. The data demonstrate that the arc-flare up event in the SCB was also coeval with high-temperature metamorphism at garnet-granulite to upper amphibolite-facies metamorphism and partial melting in the lower crust of the arc. Evidence for high-temperature metamorphism is preserved in the Cucamonga terrane (eastern San Gabriel Mountains) where metamorphic zircons in gneisses, migmatites and calc-silicates give dates ranging from 86 to 76 Ma at 9 − 7 kbars25,26 (Supplementary file). Ti-in-zircon analyses on metamorphic rims indicate temperatures of 800 to 713°C and garnet-quartz oxygen isotope thermometry yields similar metamorphic temperatures of 835 to 777°C. These data are consistent with mineral exchange thermometry which give temperatures of 800 − 775°C26. In the Coast Ridge Belt (Santa Lucia Mountains, Salinia block), Kidder et al.27 also report peak metamorphic pressures of 800°C at 7.5 kbars at 81 to 76 Ma. These results are significant because they show that the root of the arc was hot and partially molten through ca. 76 Ma.
Our regional mapping of mylonitic ductile-shear zones in the SCB documents that the Late Cretaceous arc flare-up was temporally and spatially associated with syn-plutonic development of a regionally extensive, transpressional shear system. This system includes, from north to south (present-day coordinates): the Nacimiento fault28, the Tumamait shear zone (Pine Mountain block), the Alamo Mountain-Piru Creek shear zone, the Black Belt shear zone25,29, and the Eastern Peninsular Ranges shear zone30. Kinematic indicators from the shear zones generally give oblique top-to-west, sinistral reverse- to thrust-sense motion (present-day geometry). Syn-kinematic, metamorphic titanite dates and Zr-in-titanite temperatures in the Tumamait shear zone range from 77 to 74 Ma at 750°C, titanites in the Black Belt shear zone give an age of 83 Ma, and those in the eastern Peninsular Ranges shear zone give ages of 89–78 Ma at 750°C31 (Supplementary file). In the Salinian block, Kidder et al.27 document extensive thickening and high-temperature ductile deformation from 93 − 81 Ma, and Singleton and Cloos32 show that the juxtaposition of the Salinian block against the Nacimiento block occurred between ca. 75 Ma and 60 Ma along sinistral strike-slip and/or thrust faults. Thus, these dates and temperatures indicate that the frontal arc of the SCB experienced sinistral-reverse, transpressional deformation during the Late Cretaceous flare-up event (Fig. 2C). Therefore, the collective data show simultaneous magmatism and high-temperature intra-arc shortening during the Late Cretaceous transition to thick-skinned deformation.
Implications for Late Cretaceous high-temperature arc processes in the SCB
Our compilation of igneous zircon dates from the frontal arc of the SCB document a widespread surge of Late Cretaceous magmatism from 90 − 70 Ma. This surge occurred throughout all major structural blocks which make up the ca. 500-km wide segment of the Late Cretaceous arc system (Fig. 1B, 2). Outside of the frontal arc, coeval magmatism also occurred in the back-arc which is now represented by the adjacent Mojave Desert33. The widespread occurrence of plutonic rocks throughout the SCB and their arc-like geochemical features24,34 suggests that mechanisms for generating large volumes of melt were still active through the Late Cretaceous. Arc flare-ups similar in magnitude and duration also occurred in adjacent sectors of the arc slightly before the SCB event (e.g., southern SNB at 110 − 90 Ma35 and northern PRB at 99 − 91 Ma36 (Fig. 2B). Importantly, our data from the SCB are in direct conflict with existing amagmatic models that invoke underthrusting of the conjugate Shatsky plateau beneath the SCB from 88 − 75 Ma and removal of the lower crust and lithospheric mantle during the beginning of the Laramide orogeny9,14,16.
Metamorphic ages and thermometry in the Cucamonga and Salinian granulites also demonstrate that the lower crust of the SCB was hot and partially molten through 75 Ma, and these features cannot be explained by existing amagmatic models. The presence of this hot arc root provides further evidence that high-temperature arc processes were operating until 75 − 70 Ma and shut-down of the frontal arc did not occur until after 70 Ma. This observation is illustrated in Fig. 3 which shows a compilation of time-temperature profiles derived from mineral thermochronology from the major structural blocks in the SCB. These data highlight two important features of the SCB: 1) the Late Cretaceous flare-up in the SCB was coeval with intra-arc, oblique sinistral-reverse thrusting and high-temperature metamorphism, and 2) termination of arc magmatism in this region was associated with an abrupt phase of rapid regional cooling of the SCB below 350°C at ca. 75 − 70 Ma (Fig. 3). This regionally extensive and rapid cooling event signifies the onset of regional refrigeration of the SCB due to flat-slab subduction involving the cold, conjugate Shatsky plateau and tectonic underplating of trench sediments beneath the SCB after 75 − 70 Ma37.
Post 75-70 Ma underthrusting of the conjugate Shatsky plateau
One of the key results of our work is that flat-slab subduction beneath the SCB post-dates the beginning of Laramide deformation by ~ 15 Myr. Consequently, oceanic plateau underthrusting cannot be called upon as the driver for thick-skin deformation in the western US prior to ca. 75 Ma. However, several existing data sets support a link between Laramide deformation and flat-slab subduction after 75 Ma. For example, the general timing of major thick-skin, basement-cored thrusting and basin development in Utah, Colorado, Wyoming and SW Montana occurred from 70 − 50 Ma4 which agrees well with underthrusting of the conjugate Shatsky after 75 Ma. In addition, zircon ages in peridotite xenoliths from the Colorado Plateau are consistent with derivation from the lithosphere at the base of the SCB and ~ 500 km of lateral displacement inboard to the Colorado Plateau transition zone after 70 Ma9.
In the SCB, the presence of underplated schists is commonly cited as evidence for flat-slab subduction38, and our geo- and thermochronological results also support a post-75 Ma emplacement model for the schists. Underplated schists in the Transverse Ranges have zircon age distributions with maximum depositional ages ranging from 75 − 68 Ma (Pelona Schist38,39) and amphibole and muscovite 40Ar/39Ar metamorphic ages from mafic schists are no older than 60 Ma38,40,41. Field observations in the Pelona schist also show no evidence for partial melting or intrusion by Cretaceous plutons. These textural and temporal constraints indicate that underplating of schists beneath the SCB took place after the termination of the SCB flare up event after 70 Ma. Therefore, these data are also consistent with flat-slab subduction after 70 Ma.
Spatial and temporal trends in sedimentary provenance within California forearc sediments also show a pronounced and sudden influx of continent-derived detritus to the southern California margin at ca. 75 Ma. Sharman et al.42 argued that this sudden influx reflects the development of a geomorphic breach within the Cretaceous arc and an associated rapid migration of forearc drainages into the continental interior. The timing of this breach is ~ 15 Myr younger than postulated plateau underthrusting in prior models, but is compatible with our model for the arrival of the conjugate Shatsky after 75 − 70 Ma (e.g., Fig. 3). These data are also consistent with eastward migration of magmatism away from the coastal arcs after 75 Ma which has been interpreted to reflect shallowing of the subduction angle over time4,43,44. We illustrate these features in our model (Fig. 4) which shows the tectonomagmatic evolution of the SCB-arc segment from 85 − 50 Ma.
A two-stage, model for the Laramide orogeny
Data from the SCB provide critical new information that allows us to resolve the debate about the link between upper-plate deformation in the western North American Cordillera and the kinematics and geometry of the down-going plate at the beginning of the Laramide orogeny. Our results show clearly that arc magmatism was robustly active through ca. 70 Ma, and underthrusting of the conjugate Shatsky plateau must have occurred after 75 − 70 Ma. Therefore, flat-slab subduction of the conjugate Shatsky plateau cannot explain the transition from thin- to thick-skin deformation at 90 − 80 Ma. This conclusion is supported by a number of sedimentological and thermochronological studies on Mesozoic sediments in SW Montana which provide evidence for the early onset of Laramide-style deformation well before 80 Ma10,45,46. These data are problematic from the standpoint of a single flat-slab event because the timing of basin formation predates the arrival of the conjugate Shatsky plateau in all flat-slab models9,14,16. Moreover, the southwest Montana basins are well outside the commonly cited corridor of Laramide deformation caused by flat-slab underthrusting17. These relationships coupled with our data from the SCB lead us to the conclusion that the Laramide orogeny cannot have a single driving mechanism.
Therefore, we propose that the Laramide orogeny is a composite tectonic event consisting of two distinct stages of paired mountain building and basin formation: 1) an early phase at 90 − 75 Ma, which was active during flare-up activity in the SCB, and 2) a more-widespread phase of basement-involved thrusting and basin formation in the continental interior from 75 − 40 Ma. In this model, the driving mechanism for initial Laramide deformation is not related to flat-slab subduction, but is closely linked to orogen-scale dextral transpression from ca. 100 − 85 Ma along the US Cordillera19,47,48. We attribute the second stage to flat-slab subduction of the conjugate Shatsky plateau beneath the SCB following the termination of flare-up magmatism in the SCB. These new data from the SCB show that multiple driving mechanisms are required to explain the diverse and previously conflicting datasets for the development of Laramide orogeny from 90 − 40 Ma.