Geochemical indications for the Paleocene-Eocene Thermal Maximum (PETM) and Eocene Thermal Maximum 2 (ETM-2) hyperthermals in terrestrial sediments of the Canadian Arctic

During the late Paleocene to early Eocene, clastic fluvial sediments and coals were deposited in northern high latitudes as part of the Marga­ret Formation at Stenkul Fiord (Ellesmere Island, Nunavut, Canada). Syn-sedimentary tectonic movements of the Eurekan deformation continu­ously affected these terrestrial sediments. Different volcanic ash layers occur, and unconformities subdivide the deposits into four sedimentary units. Rare vertebrate fossils indicate an early Eocene (Graybullian) age for the upper part of the Stenkul Fiord outcrop. Here, we present carbon isotope data of bulk coal, related organic-rich mud and siltstones, a plant leaf wax-derived alkane, and additional plant remains. These data provide a complete carbon isotope record of one stratigraphic section with defined unconformity positions and in relation to other Eurekan deformation features. A previously dated ash layer MA-1 provided a U-Pb zircon age of 53.7 Ma and is used as a stratigraphic tie point, together with a discrete negative carbon isotope excursion found above MA-1 in a closely sampled coal seam. The excursion is identified as the likely expression of the I-1 hyperthermal event. Based on our isotope data that reflect the early Eocene dynamics of the carbon cycle, this tie point, and previous paleontological constraints from vertebrate fossils, the locations of the Paleocene-Eocene Thermal Maximum (PETM) and Eocene Thermal Maximum 2 (ETM-2) hyperthermals and their extent along the complete section are herein identified. Within the intervals of the PETM and ETM-2 hyperthermal events, increasing amounts of clastic sediments reached the site toward the respective end of the event. This is interpreted as a response of the fluvial depositional system to an intensified hydrological system during the hyperthermal events. Our study establishes an enhanced stratigraphic framework allowing for the calcula­tion of average sedimentation rates of different intervals and considerations on the completeness of the stratigraphic record. As one of the few high-latitude outcrops of early Eocene terrestrial sediments, the Stenkul Fiord location offers further possibilities to study the effects of extreme warming events in the Paleogene.

Tectonostratigraphy and major structures of the Georgian Greater Caucasus: Implications for structural architecture, along-strike continuity, and orogen evolution

Although the Greater Caucasus Mountains have played a central role in absorbing late Cenozoic convergence between the Arabian and Eurasian plates, the orogenic architecture and the ways in which it accommodates modern shortening remain debated. Here, we addressed this problem using geologic mapping along two transects across the southern half of the western Greater Caucasus to reveal a suite of regionally coherent stratigraphic packages that are juxtaposed across a series of thrust faults, which we call the North Georgia fault system. From south to north within this system, stratigraphically repeated ~5–10-km-thick thrust sheets show systematically increasing bedding dip angles (<30° in the south to subvertical in the core of the range). Likewise, exhumation depth increases toward the core of the range, based on low-temperature thermochronologic data and metamorphic grade of exposed rocks. In contrast, active shortening in the modern system is accommodated, at least in part, by thrust faults along the southern margin of the orogen. Facilitated by the North Georgia fault system, the western Greater Caucasus Mountains broadly behave as an in-sequence, southward-propagating imbricate thrust fan, with older faults within the range progressively abandoned and new structures forming to accommodate shortening as the thrust propagates southward. We suggest that the single-fault-centric “Main Caucasus thrust” paradigm is no longer appropriate, as it is a system of faults, the North Georgia fault system, that dominates the architecture of the western Greater Caucasus Mountains.

Cooling of the continental plate during flat-slab subduction

Most flat-slab subduction regions are marked by an absence of arc volcanism, which is consistent with closure of the hot mantle wedge as the subducting plate flattens below the continent. Farther inland, low surface heat flow is observed, which is generally attributed to cooling of the continent by the underlying flat slab. However, modern flat slabs have only been in place for <20 Ma, and it is unclear whether there has been sufficient time for cooling to occur. We use numerical models to assess temporal variations in continental thermal structure during flat-slab subduction. Our models show that the flat slab leads to continental cooling on timescales of tens of millions of years. Cool slab temperatures must diffuse through the continental lithosphere, resulting in a delay between slab emplacement and surface cooling. Therefore, the timescales primarily depend on the flat-slab depth with shallower slabs resulting in shorter timescales. The magnitude of cooling increases for a shallow or long-lived flat slab, old subducting plate, and fast convergence rates. For regions with flat slabs at 45–70 km depth (e.g., Mexico and Peru), shallow continental cooling initiates 5–10 Ma after slab emplacement, and low surface heat flow in these regions is largely explained by the presence of the flat slab. However, for the Pampean region in Chile, with an ~100-km-deep slab, our models predict that conductive cooling has not yet affected the surface heat flow. The low heat flow observed requires additional processes such as advective cooling from the infiltration of fluids released through dehydration of the flat slab.

The spatial association of accessory minerals with biotite in granitic rocks from the South Mountain Batholith, Nova Scotia, Canada

We use mineral liberation analysis (MLA) to quantify the spatial association of 15,118 grains of accessory apatite, monazite, xenotime, and zircon with essential biotite, and clustered with themselves, in a peraluminous biotite granodiorite from the South Mountain Batholith in Nova Scotia (Canada). A random distribution of accessory minerals demands that the proportion of accessory minerals in contact with biotite is identical to the proportion of biotite in the rock, and the binary touching factor (percentage of accessory mineral touching biotite divided by modal proportion of biotite) would be ~1.00. Instead, the mean binary touching factors for the four accessory minerals in relation to biotite are: apatite (5.06 for 11,168 grains), monazite (4.68 for 857 grains), xenotime (4.36 for 217 grains), and zircon (5.05 for 2876 grains). Shared perimeter factors give similar values. Accessory mineral grains that straddle biotite grain boundaries are larger than completely locked, or completely liberated, accessory grains. Only apatite-monazite clusters are significantly more abundant than expected for random distribution. The high, and statistically significant, binary touching factors and shared perimeter factors suggest a strong physical or chemical control on their spatial association. We evaluate random collisions in magma (synneusis), heterogeneous nucleation processes, induced nucleation in passively enriched boundary layers, and induced nucleation in actively enriched boundary layers to explain the significant touching factors. All processes operate during the crystallization history of the magma, but induced nucleation in passively and actively enriched boundary layers are most likely to explain the strong spatial association of phosphate accessories and zircon with biotite. In addition, at least some of the apatite and zircon may also enter the granitic magma as inclusions in grains of Ostwald-ripened xenocrystic biotite.

Multiple sediment incorporation events in a continental magmatic arc: Insight from the metasedimentary rocks of the northern North Cascades, Washington (USA)

The rheology and composition of arc crust and the overall evolution of continental magmatic arcs can be affected by sediment incorporation events. The exhumed Cretaceous–Eocene North Cascades arc exposes abundant metasedimentary rocks that were incorporated into the arc during multiple events. This study uses field relationships, detrital zircon geochronology, bulk rock geochemistry, geothermometry, and quartz­in­garnet geobarometry to distinguish approximate contacts and emplacement depths for different metasedimentary units to better understand their protolith incorporation history and impact on the arc. The Skagit Gneiss Complex is one of the main deep crustal units of the North Cascades arc. It includes metasedimentary rocks with distinct detrital zircon signatures: Proterozoic–Cretaceous (Group 1) or Triassic–Cretaceous (Group 2) zircon populations. Both metasedimentary groups achieved near­ peak metamorphic conditions of 640–800 °C and 5.5–7.9 kbar; several Group 2 samples reveal the higher pressures. A third group of metasedimentary rocks, which was previously interpreted as metamorphosed equivalents of backarc sediments (Group 3), exhibited unimodal Triassic or bimodal Late Jurassic–Early Cretaceous detrital zircon signatures and achieved near­peak conditions of 570–700 °C and 8.7–10.5 kbar. The combined field and analytical data indicate that protoliths of Group 1 and Group 2 metasedimentary rocks were successively deposited in a forearc basin and underthrusted into the arc as a relatively coherent body. Group 3 backarc sediments were incorporated into the arc along a transpressional step­over zone. The incorporation of both forearc and backarc sediments was likely facilitated by arc magmatism that weakened arc crust in combination with regional transpression.

Geochronology of the Wrangell Arc: Spatial-temporal evolution of slab-edge magmatism along a flat-slab, subduction-transform transition, Alaska-Yukon

The Wrangell Arc in Alaska (USA) and adjacent volcanic fields in the Yukon provide a long-term record of interrelations between flat-slab subduction of the Yakutat microplate, strike-slip translation along the Denali–Totschunda–Duke River fault system, and magmatism focused within and proximal to a Cretaceous suture zone. Detrital zircon (DZ) U-Pb (n = 2640) and volcanic lithic (DARL) 40Ar/39Ar dates (n = 2771) from 30 modern river sediment samples document the spatial-temporal evolution of Wrangell Arc magmatism, which includes construction of some of the largest Quaternary volcanoes on Earth. Mismatches in DZ and DARL date distributions highlight the impact of variables such as mineral fertility and downstream mixing/dilution on resulting provenance signatures. Geochronologic data document the initiation of Wrangell Arc magmatism at ca. 30–17 Ma along both sides of the Totschunda fault on the north flank of the Wrangell–St. Elias Mountains in Alaska, followed by southeastward progression of magmatism at ca. 17–10 Ma along the Duke River fault in the Yukon. This spatial-temporal evolution is attributable to dextral translation along intra-arc, strike-slip faults and a change in the geometry of the subducting slab (slab curling/steepening). Magmatism then progressed generally westward outboard of the Totschunda and Duke River faults at ca. 13–6 Ma along the southern flank of the Wrangell–St. Elias Mountains in Alaska and then northwestward from ca. 6 Ma to present in the western Wrangell Mountains. The 13 Ma to present spatial-temporal evolution is consistent with dextral translation along intra-arc, strike-slip faults and previously documented changes in plate boundary conditions, which include an increase in plate convergence rate and angle at ca. 6 Ma. Voluminous magmatism is attributed to shallow subduction-related flux melting and slab edge melting that is driven by asthenospheric upwelling along the lateral edge of the Yakutat flat slab. Magmatism was persistently focused within or adjacent to a remnant suture zone, which indicates that upper plate crustal heterogeneities influenced arc magmatism. Rivers sampled also yield subordinate Paleozoic–Mesozoic DZ and DARL age populations that reflect earlier episodes of magmatism within underlying accreted terranes and match magmatic flare-ups documented along the Cordilleran margin.

Ultrahigh-temperature granulite-facies metamorphism and exhumation of deep crust in a migmatite dome during late- to post- orogenic collapse and extension in the central Adirondack Highlands (New York, USA)

This study combines field observations, mineral and whole-rock geochemistry, phase equilibrium modeling, and U-Pb sensitive high-resolution ion microprobe (SHRIMP) zircon geochronology to investigate sillimanite-bearing felsic migmatites exposed on Ledge Mountain in the central Adirondack Highlands (New York, USA), part of an extensive belt of mid-crustal rocks comprising the hinterland of the Mesoproterozoic Grenville orogen. Phase equilibrium modeling suggests minimum peak metamorphic conditions of 960–1025 °C and 11–12.5 kbar during the Ottawan orogeny—significantly higher pressure-temperature conditions than previously determined—followed by a period of near-isothermal decompression, then isobaric cooling. Petrography reveals abundant melt-related microstructures, and pseudosection models show the presence of at least ~15%–30% melt during buoyancy-driven exhumation and decompression. New zircon data document late Ottawan (re)crystallization at ca. 1047 ± 5 to 1035 ± 2 Ma following ultrahigh-temperature (UHT) metamorphism and anatexis on the retrograde cooling path. Inherited zircon cores give a mean date of 1136 ± 5 Ma, which suggests derivation of these felsic granulites by partial melting of older igneous rocks. The ferroan, anhydrous character of the granulites is similar to that of the ca. 1050 Ma Lyon Mountain Granite and consistent with origin in a late- to post-Ottawan extensional environment. We present a model for development of a late Ottawan migmatitic gneiss dome in the central Adirondacks that exhumed deep crustal rocks including the Snowy Mountain and Oregon anorthosite massifs with UHT Ledge Mountain migmatites. Recognition of deep crustal meta-plutonic rocks recording UHT metamorphism in a migmatite gneiss dome has significant implications for crustal behavior in this formerly thickened orogen.

Late Cretaceous upper-crustal thermal structure of the Sevier hinterland: Implications for the geodynamics of the Nevadaplano

Crustal temperature conditions can strongly influence the evolution of deformation during orogenesis. The Sevier hinterland plateau in Nevada and western Utah (“Nevadaplano”) experienced a Late Cretaceous episode of shallow-crustal metamorphism and granitic magmatism. Here, we investigate the thermal history of the Nevadaplano by measuring peak thermal field gradients attained in the upper 10–20 km of the crust along an east-west transect through nine ranges in eastern Nevada and western Utah, by integrating Raman spectroscopy of carbonaceous material thermometry and published conodont alteration indices with reconstructed cross sections. Thermal field gradients of 29 ± 3 °C/km were obtained in the House and Confusion Ranges in westernmost Utah. The Deep Creek, Schell Creek, and Egan Ranges in easternmost Nevada yielded elevated gradients of 49 ± 7 °C/km, 36 ± 3 °C/km, and 32 ± 6 °C/km, respectively. Moving westward, the White Pine, Butte, Pancake, and Fish Creek Ranges exhibit gradients typically between ~20–30 °C/km. The elevated thermal gradients in easternmost Nevada are interpreted to have been attained during ca. 70–90 Ma granitic magmatism and metamorphism and imply possible partial melting at ~18 km depths. Our data are compatible with published interpretations of Late Cretaceous lithospheric mantle delamination under the Sevier hinterland, which triggered lower-crustal anatexis and the resulting rise of granitic melts. The lack of evidence for structures that could have accommodated deep burial of rocks in the nearby Northern Snake Range metamorphic core complex, combined with thermal gradients from adjacent ranges that are ~1.5–3 times higher than those implied by thermobarometry in the Northern Snake Range, further highlights the debate over possible tectonic overpressure in Cordilleran core complexes. Cross-section retro-deformation defines 73.4 ± 4.6 km (76 ± 8%) of extension across eastern Nevada and 15 km of shortening in the Eastern Nevada fold belt.

Early Pennsylvanian sediment routing to the Ouachita Basin (southeastern United States) and barriers to transcontinental sediment transport sourced from the Appalachian orogen based on detrital zircon U-Pb and Hf analysis

Carboniferous sediment dispersal from the Appalachian orogenic system (eastern United States) has become a topic of widespread interest. However, the actual pathways for continental-scale, east-to-west sediment transfer have not been documented. This study presents detrital zircon (DZ) U-Pb ages and Hf isotopic values from the Lower Pennsylvanian (Morrowan) Jackfork Group and Johns Valley Shale of the synorogenic Ouachita deepwater basin of Arkansas to document provenance and delineate the likely sediment-routing systems within the broader context of sediment dispersal across Laurentia. Twelve (12) DZ U-Pb age distributions are interpreted to indicate that sediments were derived from the Appalachians to the east and northeast, as well as the midcontinent region to the north. All samples display prominent ca. 500– 400 Ma, 1250–950 Ma, 1550–1300 Ma, and 1800–1600 Ma grains, consistent with ultimate derivation from the Appalachian, Grenville, Midcontinent, and Yavapai-Mazatzal provinces. DZ Hf values obtained from the Ouachita Basin are similar to published Hf values from Pennsylvanian samples in the Appalachian and Illinois Basins. Age distributions are generally consistent for seven samples collected from the Jackfork Group and Johns Valley Shale in the southern Ouachita Mountains through ~2400 m of stratigraphic section and are interpreted to indicate little change in provenance during the Morrowan in this part of the system. However, samples from the most northern and most source-proximal site in Little Rock, Arkansas, exhibit modest percentages of Appalachian ages and elevated contributions of Yavapai-Mazatzal ages when compared with samples collected farther to the south and west. We interpret differences between DZ signatures to indicate distinct sediment-routing pathways to the Ouachita Basin. We infer the strong Appalachian and Grenville signals to represent an axial system flowing through the Appalachian foredeep, whereas the more diverse signals represent a confluence of rivers from the northeast through the backbulge of southern Illinois and western Kentucky and from the north across the Arkoma shelf. Collectively, the Ouachita Basin represents a terminal sink for sediments derived from much of the eastern and central United States.

The rise and demise of deep accretionary wedges: A long-term field and numerical modeling perspective

Several decades of field, geophysical, analogue, and numerical modeling investigations have enabled documentation of the wide range of tectonic transport processes in accretionary wedges, which constitute some of the most dynamic plate boundary environments on Earth. Active convergent margins can exhibit basal accretion (via underplating) leading to the formation of variably thick duplex structures or tectonic erosion, the latter known to lead to the consumption of the previously accreted material and eventually the forearc continental crust. We herein review natural examples of actively underplating systems (with a focus on circum-Pacific settings) as well as field examples highlighting internal wedge dynamics recorded by fossil accretionary systems. Duplex formation in deep paleo–accretionary systems is known to leave in the rock record (1) diagnostic macro- and microscopic deformation patterns as well as (2) large-scale geochronological characteristics such as the downstepping of deformation and metamorphic ages. Zircon detrital ages have also proved to be a powerful approach to deciphering tectonic transport in ancient active margins. Yet, fundamental questions remain in order to understand the interplay of forces at the origin of mass transfer and crustal recycling in deep accretionary systems. We address these questions by presenting a suite of two-dimensional thermo-mechanical experiments that enable unravelling the mass-flow pathways and the long-term distribution of stresses along and above the subduction interface as well as investigating the importance of parameters such as fluids and slab roughness. These results suggest the dynamical instability of fluid-bearing accretionary systems causes either an episodic or a periodic character of subduction erosion and accretion processes as well as their topographic expression. The instability can be partly deciphered through metamorphic and strain records, thus explaining the relative scarcity of paleo–accretionary systems worldwide despite the tremendous amounts of material buried by the subduction process over time scales of tens or hundreds of millions of years. We finally stress that the understanding of the physical processes at the origin of underplating processes as well as the forearc topographic response paves the way for refining our vision of long-term plate-interface coupling as well as the rheological behavior of the seismogenic zone in active subduction settings.