Midcontinent Rift Research Articles

One Billion Years of Stability in the North American Midcontinent Following Two‐Stage Grenvillian Structural Inversion

AbstractThe North American craton interior preserves a >1 Ga history of near surface processes that inform ongoing debates regarding timing and drivers of continental‐scale deformation and erosion associated with far‐field orogenesis. We tested various models of structural inversion on a major segment of the Midcontinent Rift along the Douglas Fault (DF) in northern Wisconsin, which accommodated ≳10 km of total vertical displacement. U‐Pb detrital zircon and vein calcite Δ47/U‐Pb thermochronometry from the hanging wall constrain the majority of uplift (≳8.5 km) and deformation to 1052–1036 Ma during the Ottawan phase of the Grenvillian orogeny. Combined U‐Pb zircon dates, Δ47/U‐Pb calcite thermochronometry, and field data that document syn‐ to early post‐depositional deformation in the footwall constrain a second stage of uplift (1–1.5 km) ca. 995–980 Ma during the Rigolet phase of the Grenvillian orogeny. A minor phase of Appalachian far‐field orogenesis is associated with minimal thrust reactivation. Our combined analyses identified the 995–980 Ma Bayfield Group as a Grenvillian foreland basin with an original thickness 0.5–2 km greater than currently preserved. By quantifying flexural loading and other subsidence mechanisms along the Douglas Fault, we identify dynamic subsidence as a mechanism that could be consistent with the development of late‐Grenvillian transcontinental fluvial systems. Minimal post‐Grenvillian erosion (0.5–2 km) in this part of the craton interior has preserved the Bayfield Group and equivalent successions, limiting the magnitude of regional erosion that can be attributed to Neoproterozoic glaciation.

Mantle plume trail beneath the ca. 1.1 Ga North American Midcontinent Rift revealed by magnetotelluric data.

Whilst the 1.1 Ga North American Midcontinent Rift (MCR) system is formed in association with the Keweenaw mantle plume, the absence of a northern third rift arm or aulacogen (a general characteristic of mantle plumes) has previously not been well understood. To help clarify this unusual plume-rift relationship and to better establish the region affected by the Keweenaw mantle plume, we present the first electrical resistivity model of the MCR derived from 3D inversion of EarthScope USArray and Lithoprobe magnetotelluric (MT) data, extending northwards into the Archean Superior Province. Our model shows a prominent highly conductive anomaly trending NW-SE at the base of Western Superior's cratonic lithospheric mantle, cross-cutting and extending for over 300km on both sides of the western rift branch. We propose that this anomaly reflects the ancient signature of a plume trail, resulting from metasomatism and/or partial melting of the sulfide-rich basal lithospheric mantle during impingement of the Keweenaw mantle plume.

Crustal and uppermost mantle structures of the North American Midcontinent Rift revealed by joint full-waveform inversion of ambient-noise data and teleseismic P waves

The Midcontinent Rift (MCR), hosting several world-class ore deposits, is the fossil remnant of a massive Mesoproterozoic rifting event (1.1 Ga) that did not lead to the formation of an ocean basin. To better understand the lithospheric processes associated with the rifting stage and its subsequent failure, we developed a novel full-waveform joint inversion method using ambient noise data and teleseismic P waves for this seismically inactive region. We apply this approach to three years (2011-2013) of seismic recordings from the Superior Province Rifting EarthScope Experiment (SPREE) (~12 km average station spacing) and the USArray Transportable Array (~70 km average station spacing), and obtain a new 3D high-resolution Vs model down to 100 km depth, as well as Vp and density models down to 60 km depth. The model shows major velocity anomalies in agreement with previous seismic studies for the western arm of the MCR. In particular, we observe high density (2.8-3.0 g/cm3), Vp (6.3-6.5 km/s), and Vs (3.6-3.7 km/s) structures in the shallow upper crust within the rift, likely associated with volcanic rocks. Similar to a previously identified underplated layer, we also observe extensive normal-to-high Vs (3.8-4.2 km/s) along the whole rift axis and Vp (6.8-7.5 km/s) beneath the northern segment of the rift within the lower crust. However, the Vs and Vp values are lower than average for typical underplated materials. We suggest that this underplated layer may represent a combination of different intrusive rock types (e.g., gabbro, anorthosite) developed during magma differentiation processes, or contamination of the mafic magma by surrounding crustal material, or intrusions of sills.

Are the Fe-rich-clay veins in the igneous rock of the Kansas (USA) Precambrian crust of magmatic origin?

The North American Mid-Continent rift (MCR) is a 1.1 Ga aborted rift, which has recently become an area of intense focus for energy resource exploration following the report of H2 emissions. To document the nature of the producing rocks, we conducted a multi-scale study on preserved drill-core samples from the DR1-A well located in the same area as the H2-producing wells in Kansas. We showed that this well reaches an unmapped part of the MRS composed of fayalite-bearing monzo-diorites in which we identified atypical veins of iddingsite, a complex mixing of Fe-rich phyllosilicates. Combining scanning and transmission electron microscopy (SEM and TEM) with scanning transmission X-ray microscopy (STXM) allowed us to differentiate at least two types of sub-micro-meter Fe-rich veins. A central reduced vein cutting the fayalite, other mafic minerals, the plagioclase, and to a lesser extent the alkali feldspar, contains a complex mix of serpentine, chlorite, and mica. Furthermore, a border vein, with a higher degree of Fe-oxidation, is found to be contained only within the fayalite. This external vein mainly contains iron and silicon, together with a few percent of potassium and calcium, and can be divided into two sub-veins composed of Fe3+-rich interstratified chlorite-smectite and Fe3+-rich serpentine in direct contact with the fayalite. Textures and microstructures of these phyllosilicates suggest that they have crystallized from a late magmatic and differentiated fluid, which precipitation produced the central vein together with the exsolution of an H2O-enriched fluid phase. This exsolved fluid, chemically far from equilibrium with the fayalite, appears to have induced a deuteric alteration of the fayalite, leading to the crystallization of the external veins enriched in ferric iron. These observations bring new perspectives on the history of formation of iron-rich clay minerals, which may, somehow, be related to H2 production.

Genesis of Fe–Ti Oxide-Bearing Ultramafic Intrusions in the Duluth Complex, Minnesota, USA

Abstract The Duluth Complex is a large mafic intrusive system located in northeastern Minnesota emplaced as part of the 1.1-Ga Midcontinent Rift. Several Fe–Ti oxide-bearing ultramafic intrusions are hosted along the Western Margin of the Duluth Complex, and are discordant bodies present in a variety of geometries, hosted in multiple rock types, and dominated by peridotite, pyroxenite, and semi-massive to massive Fe–Ti oxide rock types. Their origin has been debated, and here we present geochemical evidence and modeling that supports a purely magmatic origin for the Titac and Longnose Fe–Ti oxide-bearing ultramafic intrusions. Ilmenite and titanomagnetite textures indicate a protracted cooling process, and δ34S values of sulfides reveal little assimilation of the footwall Virginia Formation, a fine-grained pelitic unit that contains sulfide-rich bands. We model the crystallization of a hypothetical parental magma composition to the host intrusion of Longnose using Rhyolite-MELTS and demonstrate that the accumulation of Fe–Ti oxides in the discordant intrusions cannot be explained by density-driven segregation of crystallized Fe–Ti oxides. Instead, we show that the development of silicate liquid immiscibility, occurring by the unmixing of the silicate melt into conjugate Si- and Fe-rich melts, can result in the effective segregation and transportation of the Fe-rich melt. The Fe-rich melt is ~2 orders of magnitude less viscous than the Si-rich melt, allowing the Fe-rich melt to be more effectively segregated and transported in the mush regime (crystallinities >50%). This suggests that viscosity, in addition to density, plays a significant role in forming the discordant Fe–Ti oxide-bearing ultramafic intrusions. We propose a genetic model that could also be responsible for the Fe–Ti oxide-rich layers or bands that are hosted within the igneous stratigraphy of mafic intrusions of the Duluth Complex.

Geochemical assessment of mineral sequestration of carbon dioxide in the midcontinent rift

AbstractThis study examines the potential of Midcontinent Rift rocks to facilitate long‐term CO2 sequestration by providing the necessary Ca and Mg for carbonate mineralization. Surface samples were collected from the Oronto and Bayfield‐Jacobsville Groups around Lake Superior and used for petrography and X‐ray diffraction to determine their mineral composition. Also, X‐ray fluorescence was also used to assess their bulk chemical composition. The samples were then exposed to CO2 and deionized water in Teflon‐lined vessels at 90°C, and the resulting leachate fluids were analyzed for the cation released during the testing. SEM microscopy was used to examine the samples for potential mineralization of carbonate minerals. The Oronto Group sediments consist primarily of feldspathic to feldspathic lithic arenites with a chlorite‐dominated matrix, and the primary porosity is blocked by calcite and hematite cement. The Bayfield–Jacobsville sequences are porous quartz arenites to feldspathic quartz arenites that do not contain significant accumulation of Ca‐, Mg‐, and Fe‐bearing minerals. The leachate fluids obtained from Oronto Group samples exhibit a maximum Ca release rate (5.2 × 10−4 mole/cm2.day), indicating rapid calcite cement dissolution and increased porosity and permeability. SEM/EDS microanalysis revealed areas where pore‐filling calcite was preferentially dissolved. Longer‐term rock‐water reactions resulted in induced carbonate mineralization, as evidenced by calcite crystals observed in a sample reacted for 102 days. © 2024 Society of Chemical Industry and John Wiley & Sons, Ltd.

Facies analysis and depositional model of the Midcontinent Rift System in Kansas, USA

AbstractThe Midcontinent Rift System of North America is one of the oldest continental rifts but rifting ceased before continental breakup. The southern segment of the Midcontinent Rift System lies in Kansas, USA, where the stratigraphic succession and rift evolution are largely unknown. This study analysed the rift basin infill in this part of the Midcontinent Rift System to propose a depositional model. The Precambrian rift succession was described in discontinuous cores drilled in the Texaco Noel Poersch#1 well in Washington County. Sixteen lithofacies were identified and grouped into four different facies associations (fluvial, aeolian, lacustrine and alluvial fan). Overall, the studied succession comprises continental deposits accumulated dominantly in alluvial and aeolian settings, with the intermittent development of lacustrine systems. The proposed depositional models for the available core intervals indicate cyclic patterns of overfilled and underfilled phases within the rift basin. These changes in the accommodation‐to‐supply ratio were controlled by tectonism and probably modulated by climate during evolution in the syn‐rift phase. This study advances our understanding of variations across the Midcontinent Rift System.

High-precision U-Pb geochronology links magmatism in the Southwestern Laurentia large igneous province and Midcontinent Rift

Abstract The Southwestern Laurentia large igneous province (SWLLIP) comprises voluminous, widespread ca 1.1 Ga magmatism in the southwestern United States and northern Mexico. The timing and tempo of SWLLIP magmatism and its relationship to other late Mesoproterozoic igneous provinces have been unclear due to difficulties in dating mafic rocks at high precision. New precise U-Pb zircon dates for comagmatic felsic segregations within mafic rocks reveal distinct magmatic episodes at ca. 1098 Ma (represented by massive sills in Death Valley, California, the Grand Canyon, and central Arizona) and ca. 1083 Ma (represented by the Cardenas Basalts in the Grand Canyon and a sill in the Dead Mountains, California). The ca. 1098 Ma magmatic pulse was short-lived, lasting 0.25−0.24+0.67 m.y., and voluminous and widespread, evidenced by the ≥100 m sills in Death Valley, the Grand Canyon, and central Arizona, consistent with decompression melting of an upwelling mantle plume. The ca. 1083 Ma magmatism may have been generated by a secondary plume pulse or post-plume lithosphere extension. The ca. 1098 Ma pulse of magmatism in southwestern Laurentia occurred ~2 m.y. prior to an anomalous renewal of voluminous melt generation in the Midcontinent Rift of central Laurentia that is recorded by the ca. 1096 Ma Duluth Complex layered mafic intrusions. Rates of lateral plume spread predicted by mantle plume lubrication theory support a model where a plume derived from the deep mantle impinged near southwestern Laurentia, then spread to thinned Midcontinent Rift lithosphere over ~2 m.y. to elevate mantle temperatures and generate melt. This geodynamic hypothesis reconciles the close temporal relationships between voluminous magmatism across Laurentia and provides an explanation for that anomalous renewal of high magmatic flux within the protracted magmatic history of the Midcontinent Rift.

End-Mesoproterozoic (ca. 1.08 Ga) epeiric seaway of the Nonesuch Formation, Wisconsin and Michigan, USA

Abstract The Nonesuch Formation and related sedimentary units of the Oronto Group, southern Lake Superior region, midwestern United States, are commonly held to have been deposited in a lacustrine rift basin within interior continental Laurentia. Here, we present new sedimentologic and stratigraphic evidence that shows a marine influence on deposition. Tidally influenced shallow-marine sandstone and evaporitic, sandy and muddy tidal flat facies pass upward into fine-grained estuarine and sandy turbidite deposits, which are sharply overlain by mixed sandy and muddy tidal flat and floodplain deposits. These observations are evidence that the lower Oronto Group was deposited in an epeiric seaway, one of several such seaways that developed during the final amalgamation of Rodinia at a time of globally high sea level. Retrogradational-aggradational-progradational-degradational stratal architecture records changes in the relative balance between generation of accommodation space and sedimentation rates, which we interpret to reflect the combined influence of Grenvillian Ottawan phase tectonic subsidence and thermal subsidence from earlier Midcontinent Rift magmatism. We use this revised stratigraphic framework to show that the geochemical proxies of the Nonesuch Formation are tied closely to sedimentary facies and reflect intrabasinal redox heterogeneity rather than global anoxia at the end-Mesoproterozoic. Further, our sedimentology shows that the microfossils recovered from the Nonesuch rocks are primarily associated with tidal flat facies. The combined influence of marine and local nonmarine conditions must be considered when invoking the Nonesuch Formation, or similar marine-influenced interior basin deposits, as global analogues.

Joint Inversion of SPREE Receiver Functions and Surface Wave Dispersion Curves for 3‐D Crustal and Upper Mantle Structure Beneath the U.S. Midcontinent Rift

AbstractBroadband seismograms from the EarthScope Transportable Array and Superior Province Rifting EarthScope Experiment (SPREE) deployments are used to map the crust and uppermost mantle structures beneath the failed Midcontinent Rift (MCR) of Minnesota/Wisconsin, USA. The results suggest the existence of a variable zone of mafic underplating that is up to 20 km thick (40–60 deep). We jointly invert receiver functions and Rayleigh wave dispersion curves to quantify the region's crustal and mantle shear‐wave velocity structure. Basin sediment thicknesses are mildly asymmetric about the rift axis, with thickest regions immediately beneath the rift. 3‐D modeling shows anomalous lower crust and crust‐mantle transitions beneath the MCR. Sub‐MCR crustal thicknesses are generally >50 km with lower crust Vs of 4.0–4.2 km/s. Away from the MCR, the crust is typically ∼40 km thick. Strong variations in apparent crustal thickness are found along the MCR, increasing significantly in places. An additional layer of shear velocities intermediate between typical lower crust and upper mantle velocities (4.1–4.6 km/s) exists beneath most of the MCR which is thickest beneath the rift axis and pinches out away from the rift. This structure corroborates previous proposals of the presence of an underplated layer near the Moho. Results cannot distinguish between different mechanisms of emplacement (e.g., mafic interfingering within a subsequently down‐dropped lower crust vs. development of a high‐density pyroxenitic residuum at the top of the mantle). Also observed are anomalously high (>4.7 km/s) sub‐rift shear‐wave velocities at ∼70–90‐km depths, suggesting the presence of cold, depleted upper mantle material.

Post-Archean Nb-REE-U enrichment in the Superior craton recorded in metasomatised mantle rocks erupted in the 1.1 Ga Midcontinental Rift event

Post-Archean Nb-REE-U enrichment in the Superior craton recorded in metasomatised mantle rocks erupted in the 1.1 Ga Midcontinental Rift event

Reconstructing the paleoenvironment of an oxygenated Mesoproterozoic shoreline and its record of life

The Nonesuch Formation microbiota provide a window into ca. 1075 Ma life within the interior of ancient North America. The Nonesuch water body formed following the cessation of widespread volcanism within the Midcontinent Rift as the basin continued to subside. In northern Michigan and Wisconsin, USA, the Copper Harbor Conglomerate records terrestrial alluvial fan and fluvial plain environments that transitioned into subaqueous lacustrine deposition of the Nonesuch Formation. These units thin toward a paleotopographic high associated with the Brownstone Falls angular unconformity. Due to these “Brownstone Highlands,” we were able to explore the paleoenvironment laterally at different depths in contemporaneous deposits. Rock magnetic data constrain that when the lake was shallow, it was oxygenated as evidenced by an oxidized mineral assemblage. Oxygen levels were lower at greater depth—in the deepest portions of the water body, anoxic conditions are recorded. An intermediate facies in depth and redox between these endmembers preserves detrital magnetite and hematite, which can be present in high abundance due to the proximal volcanic highlands. This magnetic facies enabled the development of a paleomagnetic pole based on both detrital magnetite and hematite that constrains the paleolatitude of the lake to 7.1 ± 2.8°N. Sediments of the intermediate facies preserve exquisite organic-walled microfossils, with microfossils being less diverse to absent in the anoxic facies where amorphous organic matter is more likely to be preserved. The assemblage of cyanobacteria and eukaryotes (both photoautotrophs and heterotrophs) lived within the oxygenated waters of this tropical Mesoproterozoic water body.

Fluid Flow, Alteration, and Timing of Cu-Ag Mineralization at the White Pine Sediment-Hosted Copper Deposit, Michigan, USA

Abstract White Pine, located in Michigan’s Upper Peninsula, is an archetypal sediment-hosted stratiform copper deposit. The Midcontinent rift system is one of only seven basins globally that host a giant sediment-hosted stratiform copper deposit. Despite many similarities with other deposits of this type, White Pine displays some important differences, including the late Mesoproterozoic age, a thick basalt sequence, an apparent lack of evaporites, and a lacustrine depositional setting. This study analyzes paleofluid flow related to the formation of White Pine and places a particular emphasis on structural and diagenetic fluid pathways. Most of the ore is located in a 30-m-wide zone spanning the Copper Harbor Formation red beds and the overlying Nonesuch Formation shales. Sedimentation of these units was accompanied by subtle synsedimentary faulting. Premineralization phases include calcite concretions and nodules, illite and hematite grain coatings, isopachous chlorite rims, emplacement of liquid petroleum (now pyrobitumen), and bleaching. Mineralization introduced native copper into the footwall sandstones and a zoned suite of native copper and sulfur-poor copper sulfide minerals across a migrating redox front in the overlying shales where copper minerals nucleated on authigenic and detrital chlorite grains. Postmineralization phases include quartz cement, calcite cement, and calcite veins that partially overlapped inversion of synsedimentary faults. Contrary to previous studies, we identified evidence for only one phase of mineralization. An Re-Os chalcocite age of 1067 ± 11 Ma places mineralization 11 to 17 m.y. after host-rock deposition. Sulfide δ34S values of –14.0 to 29.9‰ suggest an important contribution from sour gas and thermochemical sulfate reduction of seawater. Carbon (δ13C) and oxygen (δ18O) isotope compositions of five calcite generations range from –15.1 to –1.3‰ and 10.4 to 41.3‰, respectively, and record early meteoric pore water displaced by later seawater. White Pine is both a sediment-hosted stratiform copper deposit and a paleo-oil field. Synsedimentary faults controlled the sedimentological character of the upper Copper Harbor Formation, and together these imparted a strong control on fluid flow and later diagenetic processes. Early oxidized meteoric fluids were displaced by liquid petroleum and sour gas, which were in turn succeeded by metal-rich but sulfate-poor oxidized seawater. Burial compaction during deposition of the overlying Freda Formation drove fluids through White Pine due to its situation on a paleotopographic high near the basin margin. Mineralization occurred at ~125°C at depths of ~2.0 km and spanned incipient basin inversion related to the distal effects of Grenvillian orogenesis. The hightenor copper mineral assemblage is the product of an abundant supply of metal from basaltic volcanic detritus in the Copper Harbor Formation and low seawater sulfate concentrations in late Mesoproterozoic oceans. This demonstrates that viable sediment-hosted stratiform copper systems can form when a readily leachable metal source rock is present, even if hypersaline and sulfate-rich brines are not.

Characterizing the supra- and subsolidus processes that generated the Current PGE-Cu-Ni deposit, Thunder Bay North Intrusive Complex, Canada: insights from trace elements and multiple S isotopes of sulfides.

The online version contains supplementary material available at 10.1007/s00126-023-01193-9.

Controls on the Stratiform Copper Mineralization in the Western Syncline, Upper Peninsula, Michigan

The Western Syncline hosts reduced-facies, or Kupferschiefer-type, sedimentary rock-hosted stratiform Cu deposits (SSC) in the lowermost meters of the Nonesuch Formation, which is part of a thick section of clastic sedimentary rocks that comprise the upper fill of the Mesoproterozoic Midcontinent Rift of North America. Located in the Porcupine Mountains Cu district in Upper Peninsula, Michigan, these blind deposits were discovered in 1956, but are not yet developed, although recent renewed interest may result in near-term production. The deposits are distinguished by their relatively undeformed nature and lack of superposed hydrothermal events. Prior to lithification, chalcocite mineralization replaced diagenetic pyrite within two discrete tabular, albeit discontinuous, potential orebodies referred to as the lower Cu-bearing sequence (LCBS) and the upper Cu-bearing sequence (UCBS). The Top Cu Zone transgresses lithologic boundaries, suggesting that a limited volume of Cu-bearing fluids moved vertically upwards through the unlithified stratigraphy, since reductant pyritic rocks above this zone are essentially barren of Cu. The total Cu inventory that has a reasonable expectation of economic extraction is 3678 M lbs. of Cu with 15.3 M oz. of byproduct Ag. When a cutoff grade of 0.9% Cu over a minimum thickness of 2 m is applied to justify an underground room-and-pillar mine, the LCBS and UCBS are not continuous over the Western Syncline. Sedimentology is the first-order control of potential ore and its continuity; dark-gray shales and siltstones deposited under low-energy, anoxic conditions are preferred host rocks, whose thickness must be >2 m to be potential ore since host-rock thickness determines economic viability of extraction. Furthermore, stratigraphy influences the time constraints on mineralization as the lithification process impedes vertical permeability and thus the flow of Cu-bearing fluids upward through the unlithified section. Syn-sedimentary tectonic movements, likely along pre-existing buried faults, are a third-order control as the thickness of host rocks is enhanced under such conditions. Therefore, an understanding of the depositional and tectonic history throughout the Western Syncline is fundamental to understanding the limits of possible economic exploitation and to optimizing ore extraction.

The ~1.1 Ga St. Ignace Island complex, Northern Ontario, Canada: Evidence for magma mixing and crustal melting in the generation of Midcontinent Rift-related bimodal magmas and implications for regional metallogeny

Abstract The St. Ignace Island complex in Northern Ontario is a package of dominantly felsic rocks emplaced within the upper portions of the Osler Volcanic rocks of the ~1.1 Ga Midcontinent Rift System. The Osler volcanic rocks are predominantly tholeiitic basalts intercalated with rare interflow sediments and rhyolites. The St. Ignace Island complex is a ~26 km2 stock with a felsic core of quartz-feldspar-phyric rhyolites and dacites and an outer ring of anorthosite and gabbro. Textures at a variety of scales within the rocks of the complex show clear evidence of the mingling and mixing of partially crystallized mafic and felsic liquids. Two multigrain (zircon/baddeleyite) fractions from a sample of the gabbro define a Discordia line with an upper intercept date of 1107±8.9 Ma. The core of the complex consists of dacites and rhyolites with similar REE abundances with negative Nb anomalies whereas the surrounding mafic rocks are gabbros to monzogabbros that are less LREE-enriched than the felsic rocks but with similar HREE. Felsic units have a narrow range of 87Sr/86Sri (0.7032-0.7045) and 143Nd/144Ndi (0.51051-0.51057) whereas the mafic end members have similar 87Sr/86Sri (0.7040-0.7061) but more radiogenic 143Nd/144Ndi (0.51067-0.51085). Very well-preserved silicate melt inclusions (MI), many completely glassy, were observed in quartz, clinopyroxene and some plagioclase phenocrysts from the complex. These represent some of the oldest unrecrystallized silicate melt inclusions described to date. Melt inclusions within quartz from the felsic volcanics are broadly rhyolitic in composition whereas MI from plagioclase in the mafic volcanics range from basalt to basaltic andesite; these felsic and mafic melt compositions are interpreted to represent the end-member liquids in the system and bulk rock analyses affirm mixtures of the two. Concentrations of Cu and Ag (in both mafic and felsic MI), and Mo (in felsic MI), are up to an order of magnitude higher in the mafic and felsic MI than in continental crust. Bulk rock metal concentrations are also significantly lower than in the MI, suggesting that the melt inclusions may preserve pre-eruptive metal tenors that were subsequently modified by sulfide saturation, degassing, or post-solidus hydrothermal alteration. The whole rock and MI geochemistry of the St. Ignace complex are broadly similar to the Central Osler Group and, given the broad similar ages, suggests they may have been derived from a similar mantle source, but distinct from the source of rhyolites in the Black Bay Peninsula. The negative Nb anomalies and negative εNd values for the St. Ignace complex are consistent with mixing with older continental crust during ascent and emplacement. The rocks of the St. Ignace Island complex likely formed as the result of emplacement of a large mafic magma chamber at the base of the Osler volcanic pile that triggered partial melting to generate the rhyolite end members. The felsic melts ascended to shallower levels in the crust where they mixed with mafic magmas derived directly from the deeper chamber. Generally, melt inclusions in the complex have very high Cu and Ag contents, similar to those observed in arc-related and extremely oxidized early rift-related rocks and may account for the world-class volcano-sediment-hosted Cu-(Ag) deposits within the rift and the presence of small porphyry-style deposits.

A tectonostratigraphic framework for the late Mesoproterozoic Bylot basins of Laurentia

The Bylot basins of northeastern Canada and northwestern Greenland comprise the Borden, Aston-Hunting, Fury and Hecla, and Thule basins. This system of late Mesoproterozoic ( c. 1.27–1.0 Ga) sedimentary basins preserves an important record of present day northeastern Laurentia coincident with the emplacement of the Mackenzie large igneous province, the Shawinigan and Ottawan phases of the Grenville Orogeny, and the development of the Midcontinent Rift. However, establishing correlations between the sedimentary successions of the Bylot basins has been hindered by the absence of robust chronostratigraphic constraints. As a result, the degree to which these basins were interconnected, whether they share a common tectonostratigraphic history, and how their sedimentary patterns relate to regional tectonic events remain open questions. Recent Re–Os geochronology from organic-rich strata has yielded depositional ages from the Borden (1048 and 1046 Ma) and Fury and Hecla (1087 Ma) basins, which we integrate with existing models for the depositional history of these basins to derive three tectonostratigraphic assemblages from the Bylot basins. We project our refined tectonostratigraphic framework for the Borden and Fury and Hecla successions to Greenland to establish a testable hypothesis for how the Thule Supergroup fits into this tectonostratigraphic picture.

New insights on Moho depth and regional lithospheric structure of central Oklahoma based on receiver function analysis from dense seismic networks

Knowledge of the midcontinent crustal structure of North America is crucial for understanding the evolutionary history of the ancient North American craton as Laurentia grew through accretion ∼1.5 Ga to 1 Ga. Although Oklahoma has been recognized as a tectonically stable region since the Phanerozoic, its crustal structure records the earlier formation of the Mazatzal and southern Granite-Rhyolite provinces 1.6 Ga to 1.4 Ga. We present results from teleseismic receiver function analysis applied to 221 events recorded on 169 broadband stations in central Oklahoma. Our findings include a Moho depth map of central Oklahoma based on stacked and depth-converted teleseismic P receiver functions. The results are interpreted together with gravimetric and magnetic datasets and a recently established seismic velocity model of the crust.The Moho map shows a generally flat crust-mantle boundary in central Oklahoma with an average depth of 43.5 km while we observe a sudden thickening on the crust of the northwestern part of Oklahoma where the Moho deepens to over 50 km depth. We also find a Mid-lithosphere discontinuity at the upper-most mantle in north-central Oklahoma, presented as a negative phase deepening southeastward from 60 km to 80 km. We further observe an intracrustal discontinuity at the Nemaha uplift and Anadarko shelf regions in a depth range of 17–30 km. The hypothesis of the Mid-continent Rift (MCR) extending into Oklahoma is examined in terms of the crustal structure and Moho depth variation revealed by receiver functions. We do not find evidence of Moho structure or lower-crustal underplay characteristics similar to what has been discovered in the northern part of MCR, but the intracrustal discontinuity that deepens towards the hypothesized MCR region suggests upper-crustal volcanics potentially caused by the extended expansion regime of the failed rift near the south-most termination.

Osmium isotopes record a complex magmatic history during the early stages of formation of the North American Midcontinent Rift — Implications for rift initiation

Osmium isotopes record a complex magmatic history during the early stages of formation of the North American Midcontinent Rift — Implications for rift initiation

Temporal evolution of mantle temperature and lithospheric thickness beneath the ∼1.1 Ga Midcontinent Rift, North America: Implications for rapid motion of Laurentia

The production of the magma volumes necessary for the formation of continental flood basalts is thought to result from extension-related decompression melting of thermo-chemically anomalous mantle. The relative importance of elevated mantle temperature or decompression associated with lithospheric thinning remains difficult to constrain because the temporal relationship between the initiation of extension and continental flood basalt is often ambiguous. The ∼1.1 Ga Midcontinent Rift (MCR) in North America provides an opportunity to probe these magma generation factors because the mantle thermo-chemical anomaly (Keweenaw plume) thought responsible for the continental flood basalts of the Keweenaw Large Igneous Province (LIP) intersected an existing rift. The role of elevated mantle temperatures in Keweenaw LIP magmatism remains ambiguous because a temporal decrease in MgO content in primitive lavas may reflect decreasing melting depths with time due to progressive lithospheric thinning, a temporal decrease in mantle temperature, or both. Furthermore, paleomagnetic studies show that Laurentia was moving rapidly during Keweenaw LIP magmatism, raising the question of how the plume remained connected to the MCR during 25 Myr of volcanism. Using samples from the most complete volcanic section preserved at Mamainse Point, Ontario, Canada, we constrain mantle potential temperature and lithospheric thickness by combining three major element thermobarometry models with forward modeling of rare earth element compositions. We find that the earliest stage of Keweenaw LIP magmatism (Early phase; 1109-1104 Ma) was governed by initially deeper melting (>85-100 km thick lithosphere) and high mantle potential temperatures (>1480-1630°C) related to the plume. Potential temperatures decreased to near-ambient conditions (∼1400-1445°C) during the latter portion of the Early phase, with these lower temperatures persisting throughout the Latent (1104-1098 Ma) and Main (1098-1090 Ma) phases. Major element-derived pressures of melt equilibration place broad constraints on lithospheric thickness, which decreased from ∼60-110 km at the onset of magmatism to 45-65 km by the Main phase. The decrease in primitive MgO contents thus reflects thinning lithosphere and decreasing mantle potential temperatures. We argue that the temporal decrease in mantle potential temperature reflects the severing of the connection between the MCR and Keweenaw plume after the Early phase of magmatism by the rapid plate motion. Our constraints on mantle potential temperature resolve the relationship between the MCR, the plume, Keweenaw LIP magmatism, and rapid plate motion, and indicate that an alternative mechanism, such as enhanced small-scale convection due to rapid plate motion and plume-lithosphere interaction, was required to generate the voluminous Main magmatic phase.