Provenance signature of a forearc basin modified by spreading ridge subduction: Detrital zircon geochronology and detrital modes from the Paleogene Arkose Ridge Formation, southern Alaska

Upper Paleocene–Eocene boulder conglomerate, cross-stratified sandstone, and laminated carbonaceous mudstone of the Arkose Ridge Formation exposed in the southern Talkeetna Mountains record fluvial-lacustrine deposition proximal to the volcanic arc in a forearc basin modified by Paleogene spreading ridge subduction beneath southern Alaska. U-Pb ages of detrital zircon grains and modal analyses were obtained from stratigraphic sections spanning the 2,000 m thick Arkose Ridge Formation in order to constrain the lithology, age, and location of sediment sources that provided detritus. Detrital modes from 24 conglomerate beds and 54 sandstone thin sections aredominated by plutonic and volcanic clasts and plagioclase feldspar with minor quartz, schist, hornblende, argillite, and metabasalt. Westernmost sandstone and conglomerate strata contain <5% volcanic clasts whereas easternmost sandstone and conglomerate strata contain 40 to >80% volcanic clasts. Temporally, eastern sandstones andconglomerates exhibit an upsection increase in volcanic detritus from <40 to >80% volcanic clasts. U-Pb ages from >1400 detrital zircons in 15 sandstone samples reveal three main populations: late Paleocene–Eocene (60-48 Ma; 16% of all grains), Late Cretaceous–early Paleocene (85–60 Ma; 62%) and Jurassic–Early Cretaceous (200–100 Ma; 12%). A plot of U/Th vs U-Pb ages shows that >97% of zircons are <200 Ma and>99% of zircons have <10 U/Th ratios, consistent with mainly igneous source terranes. Strata show increased enrichment in late Paleocene–Eocene detrital zircons from <2% in the west to >25% in the east. In eastern sections, this younger age population increases temporally from 0% in the lower 50 m of the section to >40% in samples collected >740 m above the base. Integration of the compositional and detrital geochronologic data suggests: (1) Detritus was eroded mainly from igneous sources exposed directly north of the Arkose Ridge Formation strata, mainly Jurassic–Paleocene plutons and Paleocene–Eocenevolcanic centers. Subordinate metamorphic detritus was eroded from western Mesozoic low-grade metamorphic sources. Subordinate sedimentary detritus was eroded from eastern Mesozoic sedimentary sources. (2) Eastern deposystems received higher proportions of juvenile volcanic detritus through time, consistent with construction of adjacent slab-window volcanic centers during Arkose Ridge Formation deposition. (3)Western deposystems transported detritus from Jurassic–Paleocene arc plutons that flank the northwestern basin margin. (4) Metasedimentary strata of the Chugach accretionaryprism, exposed 20-50 km south of the Arkose Ridge Formation, did not contribute abundant detritus. Conventional provenance models predict reduced input of volcanic detritus to forearc basins during exhumation of the volcanic edifice and increasing exposure ofsubvolcanic plutons (Dickinson, 1995; Ingersoll and Eastmond, 2007). In the forearc strata of these conventional models, sandstone modal analyses record progressive increases upsection in quartz and feldspar concomitant with decreases in lithic grains, mainly volcanic lithics. Additionally, as the arc massif denudes through time, theyoungest detrital U-Pb zircon age populations become significantly older than the age of forearc deposition as the arc migrates inboard or ceases magmatism. Westernmost strata of the Arkose Ridge Formation are consistent with this conventional model. However, easternmost strata of the Arkose Ridge Formation contain sandstone modes that record an upsection increase in lithic grains accompanied by a decrease in quartz and feldspar, and detrital zircon age populations that closely match the age of deposition. This deviation from the conventional model is due to the proximity of the easternmost strata to adjacent juvenile volcanic rocks emplaced by slab-window volcanic processes. Provenance data from the Arkose Ridge Formation show that forearc basins modified by spreading ridge subduction may record upsection increases in non-arc, syndepositional volcanic detritusdue to contemporaneous accumulation of thick volcanic sequences at slab-window volcanic centers. This change may occur locally at the same time that other regions of the forearc continue to receive increasing amounts of plutonic detritus as the remnant arc denudes, resulting in complex lateral variations in forearc basin petrofacies and chronofacies.

Sedimentology, Depositional Age, and Provenance of Sedimentary and Volcanic Rocks Exposed Along Willow Creek, Eastern Susitna Basin, South-Central Alaska: Implications for Modification of a Forearc Basin by Spreading Ridge Subduction

This is the first detailed study of the westernmost portion of the outcrop belt, which extends along the western flank of the Talkeetna Mountains and includes thick, well-exposed outcrops along Willow Creek in the eastern Susitna basin. New sedimentologic, compositional, and geochronologic data were obtained from stratigraphic sections within Arkose Ridge Formation strata at Willow Creek. This data combined with new geologic mapping and geochronologic data from Willow Bench and Kashwitna River Bluff (north of Willow Creek), and from the Government Peak area (east of Willow Creek), help constrain depositional processes and source terranes that provided detritus to the westernmost Arkose Ridge Formation strata.

Benthic foraminiferal record of the Pleistocene uplift of the sedimentary deposits of the Burica Peninsula (Costa Rica-Panama) as a result of Cocos Ridge subduction beneath the Central American Arc

This study used paleobathymetric interpretations from benthic foraminiferal assemblages to determine the timing of the initial subduction of the Cocos Ridge beneath the Costa Rica-Panama Arc and the associated amount of uplift. Forty-seven samples from four stratigraphic sections of the Burica and the Armuelles formations yielded 217 identified species. Foraminiferal paleoecology and cluster analyses indicated paleobathymetric differences within and between the sections. Maximum and minimum uplift rates were calculated from the biochronology, elevation, paleobathymetry and estimated eustatic changes. Large decreases in paleobathymetry were mainly due to the initial middle Pleistocene subduction of the Cocos Ridge in less than 0.5 Ma. Uplift rates were uneven across the Burica Peninsula, as follows: La Vaca 4.5–0.8 m/ky, Rabo de Puerco 2.8–0.7 m/ky, San Bartolo–Chiquito 2.7–0.8 m/ky, and eastern coast 8.0–1.5 m/ky. These differences probably resulted from Cocos Ridge asymmetry and differential uplift of tectonic blocks.

Elastic thickness structure of the Andaman subduction zone: Implications for convergence of the Ninetyeast Ridge

We use the Bouguer coherence (Morlet isostatic response function) technique to compute the spatial variation of effective elastic thickness (T-e) of the Andaman subduction zone. The recovered T-e map resolves regional-scale features that correlate well with known surface structures of the subducting Indian plate and the overriding Burma plate. The major structure on the India plate, the Ninetyeast Ridge (NER), exhibits a weak mechanical strength, which is consistent with the expected signature of an oceanic ridge of hotspot origin. However, a markedly low strength (0< T-e <3 km) in that region, where the NER is close to the Andaman trench (north of 10 N), receives our main attention in this study. The subduction geometry derived from the Bouguer gravity forward modeling suggests that the NER has indented beneath the Andaman arc. We infer that the bending stresses of the viscous plate, which were reinforced within the subducting oceanic plate as a result of the partial subduction of the NER buoyant load, have reduced the lithospheric strength. The correlation, T-e < T-s (seismogenic thickness) reveals that the upper crust is actively deforming beneath the frontal arc Andaman region. The occurrence of normal-fault earthquakes in the frontal arc, low Te zone, is indicative of structural heterogeneities within the subducting plate. The fact that the NER along with its buoyant root is subducting under the Andaman region is inhibiting the subduction processes, as suggested by the changes in trench line, interrupted back-arc volcanism, variation in seismicity mechanism, slow subduction, etc. The low T-e and thinned crustal structure of the Andaman back-arc basin are attributed to a thermomechanically weakened lithosphere. The present study reveals that the ongoing back-arc spreading and strike-slip motion along the West Andaman Fault coupled with the ridge subduction exerts an important control on the frequency and magnitude of seismicity in the Andaman region. (C) 2013 Elsevier Ltd. All rights reserved.

Geomorphology, sedimentary processes and development of the Zenisu deep-sea channel, northern Philippine Sea

The Zenisu deep-sea channel originates on the Izu-Ogasawara island arc, and disappears in the Shikoku Basin of the Philippine Sea. The geomorphology, sedimentary processes, and the development of the Zenisu deep-sea channel were investigated on the basis of swath bathymetry, side-scan sonar imagery, submersible observations, and seismic data. The deep-sea channel can be divided into three segments according to the downslope gradient and channel orientation. They are the Zenisu Canyon, the E-W fan channel, and the trough-axis channel. The sediment fill is characterized by turbidite and debrite deposition and blocky-hummocky avalanche deposits on the flanks of the Zenisu Ridge. In the Zenisu Canyon and the Zenisu deep-sea channel, sediment transport by turbidity currents generates sediment waves (dunes) observed during the Shinkai 6500 dive 371. The development of the Zenisu Canyon is controlled by a N-S shear fault, whereas the trough-axis channel is controlled by basin subsidence associated with the Zenisu Ridge. The E-W fan channel was probably affected by the E-W fault and the basement morphology.

Volcanic ash from the Peru continental margin

During Leg 112 off Peru, volcanic material was recorded from middle Eocene to Holocene time. The petrographical and chemical composition of tephra is consistent with an origin from the Andean volcanic arc. The amount and thickness of ash layers provide valuable evidence for explosive volcanic episodicity. The first indication of volcanism was found in mid-Eocene sediments. Three volcanic pulses date from Miocene time. Two intense episodes took place in upper Pliocene and from Pleistocene to Holocene time. Pliocene-Pleistocene tephra are restricted to the southern upper-slope and shelf sites, indicating a removal of the volcanic arc and the extinction of the northern Peru volcanoes. The Cenozoic tectonic phases of the Andean margin may be correlated with the Leg 112 volcanic records. The explosive supply of evolved magmatic products succeeded the Incaic and Quechua tectonic phases. Acidic glasses are related to both andesitic and shoshonitic series. The calc-alkaline factor (CAF) of these glasses exhibited moderate magmatic variations during middle and late Miocene time. A dramatic change occurred during the Pliocene-Pleistocene, reflected in a strong CAF increase and the appearance of potassium-rich evolved shoshonitic glasses. This took place when the Nazca Ridge subduction began. This change in the magma genesis and/or differentiation conditions is probably related to thickening of the upper continental plate and to a new configuration of the Benioff Zone.

The hazard potential of the western segment of the Makran subduction zone, northern Arabian Sea

Evaluating the hazard potential of the Makran subduction zone requires understanding the previous records of the large earthquakes and tsunamis. We address this problem by searching for earthquake and tectonic proxies along the Makran Coast and linking those observations with the available constraints on historical seismicity and the tell-tale characteristics of sea floor morphology. The earthquake of Mw 8.1 of 1945 and the consequent tsunami that originated on the eastern part of the Makran are the only historically known hazardous events in this region. The seismic status of the western part of the subduction zone outside the rupture area of the 1945 earthquake remains an enigma. The near-shore shallow stratigraphy of the central part of Makran near Chabahar shows evidence of seismically induced liquefaction that we attribute to the distant effects of the 1945 earthquake. The coastal sites further westward around Jask are remarkable for the absence of liquefaction features, at least at the shallow level. Although a negative evidence, this possibly implies that the western part of Makran Coast region may not have been impacted by near-field large earthquakes in the recent past-a fact also supported by the analysis of historical data. On the other hand, the elevated marine terraces on the western Makran and their uplift rates are indicative of comparable degree of long-term tectonic activity, at least around Chabahar. The offshore data suggest occurrences of recently active submarine slumps on the eastern part of the Makran, reflective of shaking events, owing to the great 1945 earthquake. The ocean floor morphologic features on the western segment, on the contrary, are much subdued and the prograding delta lobes on the shelf edge also remain intact. The coast on the western Makran, in general, shows indications of progradation and uplift. The various lines of evidence thus suggest that although the western segment is potentially seismogenic, large earthquakes have not occurred there in the recent past, at least during the last 600 years. The recurrence period of earthquakes may range up to 1,000 years or more, an assessment based on the age of the youngest dated coastal ridge. The long elapsed time points to the fact that the western segment may have accumulated sufficient slip to produce a major earthquake.

Plate tectonic constraints on flat subduction and paleomagnetic constraints on rifting

Plate tectonics shapes our dynamic planet through the creation and destruction of lithosphere. This work focuses on increasing our understanding of the processes at convergent and divergent boundaries through geologic and geophysical observations at modern plate boundaries. Recent work had shown that the subducting slab in central Mexico is most likely the flattest on Earth, yet there was no consensus about what caused it to originate. The first chapter of this thesis sets out to systematically test all previously proposed mechanisms for slab flattening on the Mexican case. What we have discovered is that there is only one model for which we can find no contradictory evidence. The lack of applicability of the standard mechanisms used to explain flat subduction in the Mexican example led us to question their applications globally. The second chapter expands the search for a cause of flat subduction, in both space and time. We focus on the historical record of flat slabs in South America and look for a correlation between the shallowing and steepening of slab segments with relation to the inferred thickness of the subducting oceanic crust. Using plate reconstructions and the assumption that a crustal anomaly formed on a spreading ridge will produce two conjugate features, we recreate the history of subduction along the South American margin and find that there is no correlation between the subduction of a bathymetric highs and shallow subduction. These studies have proven that a subducting crustal anomaly is neither a sufficient or necessary condition of flat slab subduction. The final chapter in this thesis looks at the divergent plate boundary in the Gulf of California. Through geologic reconnaissance mapping and an intensive paleomagnetic sampling campaign, we try to constrain the location and orientation of a widespread volcanic marker unit, the Tuff of San Felipe. Although the resolution of the applied magnetic susceptibility technique proved inadequate to contain the direction of the pyroclastic flow with high precision, we have been able to detect the tectonic rotation of coherent blocks as well as rotation within blocks.

Seismic structure along transitions from flat to normal subduction: central Mexico, southern Peru, and southwest Japan

The fine-scale seismic structure of the central Mexico, southern Peru, and southwest Japan subduction zones is studied using intraslab earthquakes recorded by temporary and permanent regional seismic arrays. The morphology of the transition from flat to normal subduction is explored in central Mexico and southern Peru, while in southwest Japan the spatial coincidence of a thin ultra-slow velocity layer (USL) atop the flat slab with locations of slow slip events (SSEs) is explored. This USL is also observed in central Mexico and southern Peru, where its lateral extent is used as one constraint on the nature of the flat-to-normal transitions.

In western central Mexico, I find an edge to this USL which is coincident with the western boundary of the projected Orozco Fracture Zone (OFZ) region. Forward modeling of the 2D structure of the subducted Cocos plate using a finite-difference algorithm provides constraints on the velocity and geometry of the slab’s seismic structure in this region and confirms the location of the USL edge. I propose that the Cocos slab is currently fragmenting into a North Cocos plate and a South Cocos plate along the projection of the OFZ, by a process analogous to that which occurred when the Rivera plate separated from the proto-Cocos plate 10 Ma.

In eastern central Mexico, observations of a sharp transition in slab dip near the abrupt end of the Trans Mexican Volcanic Belt (TMVB) suggest a possible slab tear located within the subducted South Cocos plate. The eastern lateral extent of the USL is found to be coincident with these features and with the western boundary of a zone of decreased seismicity, indicating a change in structure which I interpret as evidence of a possible tear. Analysis of intraslab seismicity patterns and focal mechanism orientations and faulting types provides further support for a possible tear in the South Cocos slab. This potential tear, together with the tear along the projection of the OFZ to the northwest, indicates a slab rollback mechanism in which separate slab segments move independently, allowing for mantle flow between the segments.

In southern Peru, observations of a gradual increase in slab dip coupled with a lack of any gaps or vertical offsets in the intraslab seismicity suggest a smooth contortion of the slab. Concentrations of focal mechanisms at orientations which are indicative of slab bending are also observed along the change in slab geometry. The lateral extent of the USL atop the horizontal Nazca slab is found to be coincident with the margin of the projected linear continuation of the subducting Nazca Ridge, implying a causal relationship, but not a slab tear. Waveform modeling of the 2D structure in southern Peru provides constraints on the velocity and geometry of the slab’s seismic structure and confirms the absence of any tears in the slab.

In southwest Japan, I estimate the location of a possible USL along the Philippine Sea slab surface and find this region of low velocity to be coincident with locations of SSEs that have occurred in this region. I interpret the source of the possible USL in this region as fluids dehydrated from the subducting plate, forming a high pore-fluid pressure layer, which would be expected to decrease the coupling on the plate interface and promote SSEs.

I. A seismotectonic study of the Middle America Subduction Zone. II. Lithosphere and upper mantle structure of the Canadian Shield and Eastern North America

This thesis consists of two separate parts. Part I (Chapter 1) is concerned with seismotectonics of the Middle America subduction zone. In this chapter, stress distribution and Benioff zone geometry are investigated along almost 2000 km of this subduction zone, from the Rivera Fracture Zone in the north to Guatemala in the south. Particular emphasis is placed on the effects on stress distribution of two aseismic ridges, the Tehuantepec Ridge and the Orozco Fracture Zone, which subduct at seismic gaps. Stress distribution is determined by studying seismicity distribution, and by analysis of 190 focal mechanisms, both new and previously published, which are collected here. In addition, two recent large earthquakes that have occurred near the Tehuantepec Ridge and the Orozco Fracture Zone are discussed in more detail. A consistent stress release pattern is found along most of the Middle America subduction zone: thrust events at shallow depths, followed down-dip by an area of low seismic activity, followed by a zone of normal events at over 175 km from the trench and 60 km depth. The zone of low activity is interpreted as showing decoupling of the plates, and the zone of normal activity as showing the breakup of the descending plate. The portion of subducted lithosphere containing the Orozco Fracture Zone does not differ significantly, in Benioff zone geometry or in stress distribution, from adjoining segments. The Playa Azul earthquake of October 25, 1981, M_s=7.3, occurred in this area. Body and surface wave analysis of this event shows a simple source with a shallow thrust mechanism and gives M_o=1.3x10²⁷ dyne-cm. A stress drop of about 45 bars is calculated; this is slightly higher than that of other thrust events in this subduction zone. In the Tehuantepec Ridge area, only minor differences in stress distribution are seen relative to adjoining segments. For both ridges, the only major difference from adjoining areas is the infrequency or lack of occurrence of large interplate thrust events.

Part II involves upper mantle P wave structure studies, for the Canadian shield and eastern North America. In Chapter 2, the P wave structure of the Canadian shield is determined through forward waveform modeling of the phases P_nl, P, and PP. Effects of lateral heterogeneity are kept to a minimum by using earthquakes just outside the shield as sources, with propagation paths largely within the shield. Previous mantle structure studies have used recordings of P waves in the upper mantle triplication range of 15-30°; however, the lack of large earthquakes in the shield region makes compilation of a complete P wave dataset difficult. By using the phase PP, which undergoes triplications at 30-60°, much more information becomes available. The WKBJ technique is used to calculate synthetic seismograms for PP, and these records are modeled almost as well as the P. A new velocity model, designated S25, is proposed for the Canadian shield. This model contains a thick, high-Q, high-velocity lid to 165 km and a deep low-velocity zone. These features combine to produce seismograms that are markedly different from those generated by other shield structure models. The upper mantle discontinuities in S25 are placed at 405 and 660 km, with a simple linear gradient in velocity between them. Details of the shape of the discontinuities are not well constrained. Below 405 km, this model is not very different from many proposed P wave models for both shield and tectonic regions.

Chapter 3 looks in more detail at recordings of P_nl in eastern North America. First, seismograms from four eastern North American earthquakes are analyzed, and seismic moments for the events are calculated. These earthquakes are important in that they are among the largest to have occurred in eastern North America in the last thirty years, yet in some cases were not large enough to produce many good long-period teleseismic records. A simple layer-over-a-halfspace model is used for the initial modeling, and is found to provide an excellent fit for many features of the observed waveforms. The effects on P_nl of varying lid structure are then investigated. A thick lid with a positive gradient in velocity, such as that proposed for the Canadian shield in Chapter 2, will have a pronounced effect on the waveforms, beginning at distances of 800 or 900 km. P_nl records from the same eastern North American events are recalculated for several lid structure models, to survey what kinds of variations might be seen. For several records it is possible to see likely effects of lid structure in the data. However, the dataset is too sparse to make any general observations about variations in lid structure. This type of modeling is expected to be important in the future, as the analysis is extended to more recent eastern North American events, and as broadband instruments make more high-quality regional recordings available.

Plutonism in three orogenic pulses, Eastern Blue Ridge Province, southern Appalachians

The Eastern Blue Ridge Province of the southern Appalachians contains, in part, remnants of an Ordovician accretionary wedge complex formed during subduction of an oceanic tract before mid-Ordovician accretion with Laurentia. The Eastern Blue Ridge Province consists of metapelite and amphibolite intruded by low-K plutons, high-temperature (T > 750 degrees C) Ordovician eclogite, and other high-pressure metamafic and meta-ultramatic rocks. Felsic plutons in the Eastern Blue Ridge Province are important time markers for regional-scale tectonics, deformation, and metamorphism. Plutons were thought to be related to either Taconian (Ordovician) or Acadian (Devonian-Silurian) tectonothermal events.We dated five plutonic or metaplutonic rocks to constrain pluton crystallization ages better and thus the timing of tectonism. The Persimmon Creek gneiss yielded a protolith crystallization age of 455.7 +/- 2.1 Ma, Chalk Mountain 377.7 +/- 2.5 Ma, Mt. Airy 334 +/- 3 Ma, Stone Mountain 335.6 +/- 1.0 Ma, and Rabun 335.1 +/- 2.8 Ma. The latter four plutons were thought to be part of the Acadian Spruce Pine Suite, but instead our new ages indicate that Alleghanian (Carboniferous-Permian) plutonism is widespread and voluminous in the Eastern Blue Ridge Province. The Chattahoochee fault, which was considered an Acadian structure, cuts the Rabun pluton and thus must have been active during the Alleghanian orogeny. The new ages indicate that Persimmon Creek crystallized less than 3 m.y. after zircon crystallization in Eastern Blue Ridge eclogite and is nearly synchronous with nearby high-grade metamorphism and migmatization. The three phases of plutonism in the Eastern Blue Ridge Province correspond with established metamorphic ages for each of the three major orogenic pulses along the western flank of the southern Appalachians.

Plutonismn in three orogenic pulses, Eastern Blue Ridge Province, southern Appalachians

The Eastern Blue Ridge Province of the southern Appalachians contains, in part, remnants of an Ordovician accretionary wedge complex formed during subduction of an oceanic tract before mid-Ordovician accretion with Laurentia. The Eastern Blue Ridge Province consists of metapelite and amphibolite intruded by low-K plutons, high-temperature (T >750 °C) Ordovician eclogite, and other high-pressure metamafic and meta-ultramafic rocks. Felsic plutons in the Eastern Blue Ridge Province are important time markers for regional-scale tectonics, deformation, and metamorphism. Plutons were thought to be related to either Taconian (Ordovician) or Acadian (Devonian-Silurian) tectonothermal events. We dated five plutonic or metaplutonic rocks to constrain pluton crystallization ages better and thus the timing of tectonism. The Persimmon Creek gneiss yielded a protolith crystallization age of 455.7 ± 2.1 Ma, Chalk Mountain 377.7 ± 2.5 Ma, Mt. Airy 334 ± 3Ma, Stone Mountain 335.6 ± 1.0 Ma, and Rabun 335.1 ± 2.8 Ma. The latter four plutons were thought to be part of the Acadian Spruce Pine Suite, but instead our new ages indicate that Alleghanian (Carboniferous-Permian) plutonism is widespread and voluminous in the Eastern Blue Ridge Province. The Chattahoochee fault, which was considered an Acadian structure, cuts the Rabun pluton and thus must have been active during the Alleghanian orogeny. The new ages indicate that Persimmon Creek crystallized less than 3 m.y. after zircon crystallization in Eastern Blue Ridge eclogite and is nearly synchronous with nearby high-grade metamorphism and migmatization. The three phases of plutonism in the Eastern Blue Ridge Province correspond with established metamorphic ages for each of the three major orogenic pulses along the western flank of the southern Appalachians. © 2006 Geological Society of America.

Testing the influence of far-field topographic forcing on subduction initiation at a passive margin

Despite favourable gravitational instability and ridge-push, elastic and frictional forces prevent subduction initiation fromarising spontaneously at passive margins. Here,we argue that forces arising fromlarge continental topographic gradients are required to initiate subduction at passivemargins. In order to test this hypothesis,we use 2Dnumerical models to assess the influence of the Andean Plateau on stressmagnitudes and deformation patterns at the Brazilian passive margin. The numerical results indicate that “plateau-push” in this region is a necessary additional force to initiate subduction. As the SE Brazilianmargin currently shows no signs of self-sustained subduction, we examined geological and geophysical data to determine if themargin is in the preliminary stages of subduction initiation. The compiled data indicate that the margin is presently undergoing tectonic inversion, which we infer as part of the continental–oceanic overthrusting stage of subduction initiation. We refer to this early subduction stage as the “Brazilian Stage”, which is characterized by N10 kmdeep reverse fault seismicity at themargin, recent topographic uplift on the continental side, thick continental crust at themargin, and bulging on the oceanic side due to loading by the overthrusting continent. The combined results of the numerical simulations and passivemargin analysis indicate that the SE Brazilian margin is a prototype candidate for subduction initiation.