The geological expression of subduction termination is poorly understood due to overprinting during the collisional stage of the Wilson Cycle. The Anatolian domain of the eastern Mediterranean represents a modern system where spatial variability can be interpreted in terms of the transition from subduction to collision. Convergence in the west is accommodated by the subduction of the last remnants of Neotethyan oceanic lithosphere, while in the east, the margin has transitioned to complete continent-continent collision. In central Anatolia, however, the expression of convergence is complicated by the underthrusting of small continental fragments and attenuated continental lithosphere. By investigating variations in the geological expression of convergence across this system, we can investigate the processes that accompany the transition from subduction to collision.Spatially variable tectonomagmatic and seismic characteristics along the Anatolian margin reflect this transition. Seismic images reveal a disjointed and disaggregating subducting slab beneath central Anatolia that interacts with, and in some cases induces, mantle flow. This spatially corresponds with Miocene-to-recent volcanism that is sourced from very shallow depths (

The North and East Anatolian Fault Zones represent plate-bounding transform faults that enable the westward tectonic escape of the Anatolian Plate away from the Arabian-Eurasian collisional zone. These fault zones are both capable of hosting large (Mw > 7) seismic events, as most recently demonstrated by the extremely damaging February 2023 Kahramanmaraş earthquake sequence. This earthquake sequence highlighted that plate boundary forces in this area are distributed over a very broad region, however what controls the location, distribution, and character of this plate-bounding strike-slip system remains enigmatic. To better understand potential contributions to deformation, we compare seismic images of the lithosphere (e.g., crustal and lithospheric mantle thickness and velocity) to deformational features and seismicity near the EAFZ, as well as further west where it joins with the Anatolia-Arabia-Africa (A3) triple junction along the southeastern margin of the Anatolian escape system. We interpret that although controls on surface deformation are commonly linked to stress in the brittle upper crust, the complex deformation and seismicity patterns in this region are likely related to variations in the location and extent of the strong lithospheric mantle of the Arabian plate, which currently underthrusts Anatolia as far north as the Sürgü-Çardak fault zone (~50 km). In addition, the Arabian lithospheric mantle extends at least as far west as at least the central Adana Basin, coincident with a zone of relatively deep (>30 km) strike-slip seismogenesis that has produced Mw > 6 earthquakes. By investigating the relationship between recent geological deformation since the inception of the East Anatolian Fault (ca. 5 Ma) and the modern record of seismic structure and seismicity, we infer that the Sürgü-Çardak fault zone and its associated near-orthogonal bend reaching into the Adana Basin will be the future southeastern boundary of the Anatolian Plate escape tectonic system.

ABSTRACT The Santa Rosalía basin (Baja California Sur, México) contains a rich record of late Cenozoic volcanism, faulting, and sedimentation that provides a crucial constraint on the timing of marine flooding from the Pacific Ocean into the nascent Gulf of California oblique rift, yet the precise age of the basin is uncertain. Previous studies used reconnaissance paleomagnetic data and a 40Ar/39Ar age of 6.76 ± 0.90 Ma on the intrabasinal Cinta Colorada tuff to estimate a depositional age of ca. 7.2–6.3 Ma for the marine Boleo Formation and initial flooding of the central Gulf of California. Here, we present a large (n = 2091) detrital zircon U-Pb geochronology data set from the Boleo Formation that indicates a maximum depositional age of 6.35 ± 0.21 Ma for pumiceous sandstone at the base (below the basal limestone), a revised age of 5.86 ± 0.06 Ma for the Cinta Colorada tuff in the middle, and a maximum depositional age of 5.70 ± 0.21 Ma for the top. Detrital zircon age spectra suggest a local provenance for the Boleo Formation involving recycling from underlying Oligocene–Miocene strata in proximal source areas. Integration of detrital zircon ages with existing paleomagnetic data suggests that the lower ~30 m of the Boleo Formation accumulated during normal-polarity subchron C3An.1n (6.27–6.02 Ma), and the middle to upper Boleo Formation was deposited entirely during reverse-polarity chron C3r (6.02–5.24 Ma). We therefore reassign the depositional age span of the Boleo Formation to ca. 6.3–5.7 Ma. Although not preferred, a minimum-duration depositional model from ca. 6.1 to 5.8 Ma is also permissible if a consistently high sedimentation rate of ~0.4– 1.0 mm/yr is inferred. This revised younger age for the Boleo Formation implies marine incursion in the central Gulf of California at ca. 6.3 Ma, ~1 m.y. younger than previously thought. We envision that regional marine flooding occurred during a very short (<100 k.y.) event that inundated a narrow tectonic trough over a distance of at least ~1000 km along the plate boundary from the central Gulf of California to the Salton Trough and reaching into the present-day Lower Colorado River Valley. This study also demonstrates the utility of large-volume and large-n detrital zircon studies in establishing the ages of sedimentary successions deposited over very short time spans (<1 m.y.) and/or during relative lulls in magmatism and geomagnetic reversals.

Explanation on the creation of the lithospheric strength profiles. Table S1. Parameters used to calculate yield strength envelopes.