Ediacaran-Ordovician tectonic and geodynamic drivers of Great Unconformity exhumation on the southern Canadian Shield

Abstract

The Great Unconformity erosion surface between Paleozoic sedimentary rocks and Archean-Proterozoic basement has been related to erosion occurring during Snowball Earth glaciations, eustatic sea level fluctuations, and a range of tectonic and/or geodynamic mechanisms. Each class of mechanism predicts distinct timing and spatial patterns of exhumation. Snowball Earth glacial erosion is limited to 717-635 Ma and concentrated in narrow ice streams on continental margins. Sea-level related erosion is unconstrained in time but also spatially limited to continental margins. Tectonic and geodynamic mechanisms, in contrast, can result in exhumation distributed more broadly in time and space. We combine new zircon and apatite (U-Th)/He thermochronology data (ZHe, AHe) with independent paleodepth information from a continental interior location in southeastern Ontario, Canada to constrain the timing, magnitude, and regional pattern of exhumation associated with the Great Unconformity along a ∼650 km-long transect across the southern Canadian Shield. Here, the unconformity is defined by Middle Ordovician carbonates atop Archean-Proterozoic basement. ZHe analyses for seven basement samples display a range of dates from 960 ± 20 Ma to 37.5 ± 0.9 Ma that correlate negatively with radiation damage. AHe dates are ∼300-200 Ma and generally consistent across samples regardless of radiation damage. Independent evidence supports emplacement of both the 590 +2/-1 Ma Grenville dikes and the 577 ± 1 Ma Callander Complex in the study region at depths ≥6 km. The combined data require ≥6 km of exhumation between ca. 590-577 Ma and 470 Ma in the middle of our transect, with multi-km erosion up to ≥5 km elsewhere in the study area, well after the ca. 717-635 Ma Snowball Earth glaciations. These outcomes expand the spatial extent of an Ediacaran to early-Paleozoic exhumation signal previously inferred elsewhere across the Shield to ∼1.1 million km2. Thick Ediacaran successions on the Laurentian margins are complementary depositional signals of this erosion. The enormous spatial extent of Great Unconformity exhumation across the continental interior of the Canadian Shield is incompatible with glacial erosion and eustatic sea-level change as the primary causes. Instead, we attribute exhumation to tectonic and geodynamic mechanisms, which may include isostatic rebound, dynamic topography, plume activity associated with the Central Iapetus Magmatic Province, rifting during opening of the Iapetus Ocean, and development of the Transcontinental Arch.