North American Midcontinent Rift Research Articles

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.

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

High geomagnetic field intensity recorded by anorthosite xenoliths requires a strongly powered late Mesoproterozoic geodynamo

Obtaining estimates of Earth's magnetic field strength in deep time is complicated by nonideal rock magnetic behavior in many igneous rocks. In this study, we target anorthosite xenoliths that cooled and acquired their magnetization within ca. 1,092 Ma shallowly emplaced diabase intrusions of the North American Midcontinent Rift. In contrast to the diabase which fails to provide reliable paleointensity estimates, the anorthosite xenoliths are unusually high-fidelity recorders yielding high-quality, single-slope paleointensity results that are consistent at specimen and site levels. An average value of ∼83 ZAm2 for the virtual dipole moment from the anorthosite xenoliths, with the highest site-level values up to ∼129 ZAm2, is higher than that of the dipole component of Earth's magnetic field today and rivals the highest values in the paleointensity database. Such high intensities recorded by the anorthosite xenoliths require the existence of a strongly powered geodynamo at the time. Together with previous paleointensity data from other Midcontinent Rift rocks, these results indicate that a dynamo with strong power sources persisted for more than 14 My ca. 1.1 Ga. These data are inconsistent with there being a progressive monotonic decay of Earth's dynamo strength through the Proterozoic Eon and could challenge the hypothesis of a young inner core. The multiple observed paleointensity transitions from weak to strong in the Paleozoic and the Proterozoic present challenges in identifying the onset of inner core nucleation based on paleointensity records alone.

Sources of Hydrothermal Fluids Inferred from Oxygen and Carbon Isotope Composition of Calcite, Keweenaw Peninsula Native Copper District, Michigan, USA

The Mesoproterozoic North American Midcontinent Rift hosts the world’s largest accumulation of native copper in Michigan’s Keweenaw Peninsula. During a regional metamorphogenic-hydrothermal event, native copper was deposited along with spatially zoned main-stage minerals in a thermal high. This was followed by deposition of late-stage minerals including minor copper sulfide. Inferences from the oxygen and carbon isotopic composition of main-stage hydrothermal fluids, as calculated from 296 new and compiled isotopic measurements on calcite, are consistent with existing models that low-sulfur saline native copper ore-forming fluids were dominantly derived by burial metamorphic processes from the very low sulfur basalt-dominated rift fill at depth below the native copper deposits. Co-variation of oxygen and carbon isotopic compositions are consistent with mixing of metamorphic-derived fluids with two additional isotopically different fluids. One of these is proposed to be evolved seawater that provided an outside source of salinity. This fluid mixed at depth and participated in the formation of a well-mixed hybrid metamorphic-dominated ore-forming fluid. Secondary Ion Mass Spectrometry in-situ isotopic analyses of calcite demonstrate a high degree of variability within samples that is attributed to variable degrees of shallow mixing of the hybrid ore-forming fluid with sulfur-poor, reduced evolved meteoric water in the zone of precipitation. The oxygen and carbon isotopic compositions of 100 new and compiled measurements on late-stage calcite are mostly isotopically different than the main-stage hydrothermal fluids. The late-stage hydrothermal fluids are interpreted as various proportions of mixing of evolved meteoric water, main-stage hybrid ore-forming fluid, and shallow, evolved seawater in the relatively shallow zone of precipitation.

Final inversion of the Midcontinent Rift during the Rigolet Phase of the Grenvillian Orogeny

Abstract Despite being a prominent continental-scale feature, the late Mesoproterozoic North American Midcontinent Rift did not result in the break-up of Laurentia, and subsequently underwent structural inversion. The timing of inversion is critical for constraining far-field effects of orogenesis and processes associated with the rift's failure. The Keweenaw fault in northern Michigan (USA) is a major thrust structure associated with rift inversion; it places ca. 1093 Ma rift volcanic rocks atop the post-rift Jacobsville Formation, which is folded in its footwall. Previous detrital zircon (DZ) U-Pb geochronology conducted by laser ablation–inductively coupled plasma–mass spectrometry (LA-ICP-MS) assigned a ca. 950 Ma maximum age to the Jacobsville Formation and led researchers to interpret its deposition and deformation as postdating the ca. 1090–980 Ma Grenvillian Orogeny. In this study, we reproduced similar DZ dates using LA-ICP-MS and then dated 19 of the youngest DZ grains using high-precision chemical abrasion–isotope dilution–thermal ionization mass spectrometry (CA-ID-TIMS). The youngest DZ dated by CA-ID-TIMS at 992.51 ± 0.64 Ma (2σ) redefines the maximum depositional age of the Jacobsville Formation and overlaps with a U-Pb LA-ICP-MS date of 985.5 ± 35.8 Ma (2σ) for late-kinematic calcite veins within the brecciated Keweenaw fault zone. Collectively, these data are interpreted to constrain deposition of the Jacobsville Formation and final rift inversion to have occurred during the 1010–980 Ma Rigolet Phase of the Grenvillian Orogeny, following an earlier phase of Ottawan inversion. Far-field deformation propagated >500 km into the continental interior during the Ottawan and Rigolet phases of the Grenvillian Orogeny.

Synchronous Emplacement of the Anorthosite Xenolith‐Bearing Beaver River Diabase and One of the Largest Lava Flows on Earth

AbstractNew geochronologic and paleomagnetic data from the North American Midcontinent Rift (MCR) reveal the synchronous emplacement of the Beaver River diabase, the anorthosite xenoliths within it, and the Greenstone Flow—one of the largest lava flows on Earth. A U‐Pb zircon date of 1091.83 0.21 Ma (2) from one of the anorthosite xenoliths is consistent with the anorthosite cumulate forming as part of the MCR and provides a maximum age constraint for the Beaver River diabase. Paired with the minimum age constraint of a cross‐cutting Silver Bay intrusion (1091.61 0.14 Ma; 2), these data tightly bracket the age of the Beaver River diabase to be 1091.7 0.2 Ma (95% CI), coeval with the eruption of the Greenstone Flow (1091.59 0.27 Ma; 2)—which is further supported by indistinguishable tilt‐corrected paleomagnetic pole positions. Geochronological, paleomagnetic, mineralogical and geochemical data are consistent with a hypothesis that the Beaver River diabase was the feeder system for the Greenstone Flow. The large areal extent of the intrusives and large estimated volume of the volcanics suggest that they represent a rapid and voluminous ca. 1,092 Ma magmatic pulse near the end of the main stage of MCR magmatism.

Mesoproterozoic porphyry copper mineralization at Mamainse Point, Ontario, Canada in the context of Midcontinent rift metallogeny

The Jogran copper prospect in southern Ontario, Canada, is located on the northeastern shoulder of the North American Midcontinent rift system where it is centered on two small, closely spaced quartz monzonite porphyry intrusions emplaced into Neoarchean mafic metavolcanic rocks. Notwithstanding an atypical alteration zoning pattern, the distinctive alteration styles and veinlet types confirm its porphyry copper affiliation. Two samples of molybdenite from chalcopyrite-bearing quartz veinlets in the two centers returned Re-Os ages of 1101 ± 5 and 1094 ± 5 Ma, showing that the mineralization took place during the main stage of extension and magmatic activity that gave rise to the Midcontinent rift and associated large igneous province (LIP). The nearby copper-bearing Tribag breccia pipes and Coppercorp veins have been shown previously to also be intrusion related, and likely formed penecontemporaneously with Jogran. The Jogran and nearby copper mineralization overlapped temporally with generation of low-grade magmatic copper-nickel-PGE deposits in the Duluth intrusive complex but predated the strata-bound native copper and chalcocite deposits in the Keweenaw Peninsula, which are believed to have formed from basinal fluids during rift inversion. Although such magmatic and strata-bound copper deposits typify intracontinental rift metallogeny, porphyry copper deposits are normally confined to convergent margin magmatic arcs, thereby making the rift-LIP setting of Jogran exceptional and possibly unique.

Rapid emplacement of massive Duluth Complex intrusions within the North American Midcontinent Rift

Abstract The Duluth Complex (Minnesota, USA) is one of the largest mafic intrusive complexes on Earth. It was emplaced as the Midcontinent Rift developed in Laurentia’s interior during an interval of magmatism and extension from ca. 1109 to 1084 Ma. This duration of magmatic activity is more protracted than is typical for large igneous provinces interpreted to have formed from decompression melting of upwelling mantle plumes. While the overall duration was protracted, there were intervals of more voluminous magmatism. New 206Pb/238U zircon dates for the anorthositic and layered series of the Duluth Complex constrain these units to have been emplaced ca. 1096 Ma in <1 m.y. (duration of 500 ± 260 k.y.). Comparison of paleomagnetic data from these units with Laurentia’s apparent polar wander path supports this interpretation. This rapid emplacement bears similarities to the geologically short duration of well-dated large igneous provinces. These data support hypotheses that call upon the co-location of lithospheric extension and anomalously hot upwelling mantle. This rapid magmatic pulse occurred >10 m.y. after initial magmatism following >20° of latitudinal plate motion. A likely scenario is one in which upwelling mantle encountered the base of Laurentian lithosphere and flowed via “upside-down drainage” to locally thinned lithosphere of the Midcontinent Rift.

Integrated geophysical analysis provides an alternate interpretation of the northern margin of the North American Midcontinent Rift System, Central Lake Superior

The Midcontinent Rift System (MRS) is a 1.1 Ga sequence of voluminous basaltic eruptions and multiple intrusions followed by widespread sedimentation that extends across the Midcontinent and northern Great Lakes region of North America. Previous workers have commonly used seismic-reflection data (Great Lakes International Multidisciplinary Program on Crustal Evolution [GLIMPCE] line A) to demonstrate that the northern rift margin in central Lake Superior developed as a normal growth fault that was structurally inverted to a reverse fault during a compressional event after rifting had ended. A prominent, curvilinear aeromagnetic anomaly that extends from Isle Royale, Michigan, to Superior Shoal in central Lake Superior, Ontario (the IR-SS anomaly), is commonly presented as a manifestation of this reverse fault. We have integrated multidisciplinary geophysical analyses (seismic-reflection, seismic-refraction, aeromagnetic, and gravity), physical-property information (density, magnetic susceptibility and remanence, and compressional-wave velocity), and geologic concepts to develop an alternate interpretation of the rift margin along GLIMPCE line A, where it intersects the IR-SS anomaly. Our new model indicates that a normal fault is the dominant structure at the northern rift margin along line A, contrary to the original rift-margin paradigm, which asserts that compressional structures are the dominant features preserved today. Integral to this alternate model is a newly interpreted, prerift sedimentary basin intruded by sills in northern Lake Superior. Our alternate model of the northern rift margin has implications for interpreting the style, scale, and timing of extension, rift-related intrusion, and compression during development of the MRS.

Failed rifting and fast drifting: Midcontinent Rift development, Laurentia’s rapid motion and the driver of Grenvillian orogenesis

Author(s): Swanson-Hysell, NL; Ramezani, J; Fairchild, LM; Rose, IR | Abstract: The late Mesoproterozoic was a time of large-scale tectonic activity both in the interior and on the margins of Laurentia-most notably the development of the Midcontinent Rift and the Grenvillian orogeny. Volcanism within the North American Midcontinent Rift between ca. 1109 and 1083 Ma, as well as other contemporaneous volcanism within Laurentia, has provided an opportunity to develop extensive paleomagnetic data sets spanning this time period. These data result in an apparent polar wander path (APWP) for Laurentia that goes from a high-latitude apex known as the Logan Loop into a swath known as the Keweenawan Track. A longstanding challenge of these data was the appearance of asymmetry between relatively steep reversed polarity directions from older rift rocks and relatively shallow normal polarity directions from younger rift rocks. This asymmetry was used to support an interpretation that there were large non-dipolar components to the geomagnetic field at the time. Recent data sets support the interpretation that this directional change was progressive and therefore a result of very rapid motion of Laurentia from high to low latitudes rather than a stepwise change across non-dipolar reversals. We present high-precision U-Pb dates from Midcontinent Rift volcanics that result in an improved chronostratigraphic framework for rift volcanics and unconformities that improves correlations as well as constraints on rift development. We use these dates in volcanostratigraphic context to temporally constrain a new compilation of Midcontinent Rift paleomagnetic poles. These paleomagnetic poles include new data from the North Shore Volcanic Group, Minnesota, USA and the Osler Volcanic Group, Ontario, Canada. The U-Pb dates constrain the rate of implied plate motion more precisely than has previously been possible. We apply a novel Bayesian approach to assess the rate of implied plate motion through inverting for paleomagnetic Euler poles. If the path is to be explained by a single Euler pole these inversions reveal that motion of the continent exceeded 27 cm/yr. The path is particularly well-explained by a model wherein there is continuous true polar wander in addition to rapid plate motion that changes direction and slows ca. 1096 Ma. Laurentia's movement from high to low latitudes resulted in collisional tectonics on its leading margin which could be associated with such a change in plate motion. We propose that upwelling of the Keweenawan mantle plume was associated with an avalanche of subducted slab material and associated downwelling that drove fast plate motion. This fast plate motion was followed by the Grenvillian orogeny from ca. 1090 to 980 Ma. Prolonged collisional orogenesis could have been sustained due to this strong convective cell that therefore played an integral role in the assembly of the supercontinent Rodinia.

Seismic Imaging of the North American Midcontinent Rift Using S-to-P Receiver Functions.

North America's ~1.1‐Ga failed Midcontinent Rift (MCR) is a striking feature of gravity and magnetic anomaly maps across the continent. However, how the rift affected the underlying lithosphere is not well understood. With data from the Superior Province Rifting Earthscope Experiment and the USArray Transportable Array, we constrain three‐dimensional seismic velocity discontinuity structure in the lithosphere beneath the southwestward arm of the MCR using S‐to‐P receiver functions. We image a velocity increase with depth associated with the Moho at depths of 33–40 ± 4 km, generally deepening toward the east. The Moho amplitude decreases beneath the rift axis in Minnesota and Wisconsin, where the velocity gradient is more gradual, possibly due to crustal underplating. We see hints of a deeper velocity increase at 61 ± 4‐km depth that may be the base of underplating. Beneath the rift axis further south in Iowa, we image two distinct positive phases at 34–39 ± 4 and 62–65 ± 4 km likely related to an altered Moho and an underplated layer. We image velocity decreases with depth at depths of 90–190 ± 7 km in some locations that do not geographically correlate with the rift. These include a discontinuity at depths of 90–120 ± 7 km with a northerly dip in the south that abruptly deepens to 150–190 ± 7 km across the Spirit Lake Tectonic Zone provincial suture. The negative phases may represent a patchy, frozen‐in midlithosphere discontinuity feature that likely predates the MCR and/or be related to lithospheric thickness.

Paleomagnetism and Geochemistry of ~1144‐Ma Lamprophyre Dikes, Northwestern Ontario: Implications for the North American Polar Wander and Plate Velocities

AbstractWe present new paleomagnetic and geochemical data from a suite of the ~1144‐Ma ultramafic lamprophyre dikes that outcrop in the Canadian Shield northeast of Lake Superior (Ontario, Canada). Nineteen of 22 sampled dikes yielded consistent characteristic remanent magnetization directions of normal (n = 5) and reversed (n = 14) polarity. The primary origin of characteristic remanent magnetization is bolstered by positive baked contact tests and a reversal test. The group mean direction (D = 306.4°, I = 72.1°, α95 = 5.5°, N = 19) obtained from the lamprophyre dikes is statistically indistinguishable from the group mean direction (D = 297.4°, I = 65.5°, α95 = 8.3°, N = 8) previously reported for the nearly coeval ~1142‐Ma Abitibi dikes. The geochemistry of the lamprophyre dikes suggests strong affinity with magmas derived from ocean island basalt‐type mantle sources, consistent with the mantle plume hypothesis for the formation of the ~1.1‐Ga North American Midcontinent Rift. The similarity in age, trend, paleomagnetism, and geochemistry indicates that the lamprophyre and Abitibi dike suites represent the earliest magmatic event associated with the commencement of rifting. The combined mean direction (D = 303.1°, I = 70.2°, α95 = 4.5°, N = 27) corresponds to a paleomagnetic pole at Plat = 55.8°N, Plong = 220.0°E (A95 = 7.3°). The new pole merits the highest classification on the Q‐scale of paleomagnetic reliability and represents a key pole defining the North American apparent polar wander path during the late Mesoproterozoic. Combined with high‐quality data from the ~1108‐Ma Coldwell Complex, our data indicate an equatorward motion of Laurentia at 3.8 ± 1.4 cm/year, comparable with the present‐day velocities of continental plates, before switching to extremely rapid motion between ~1108 and ~1099 Ma.

Emplacement of ultramafic-carbonatite intrusions along reactivated North American mid-continent rift structures

Understanding the various emplacement pathways of ultramafic‑carbonatite intrusive complexes and mapping of their global distribution are in an early stage. Discerning the pattern of their emplacement can reveal information about earth's dynamics and guide exploration for associated economic deposits. Here we apply multiple spatial analysis methods and synthetic aperture radar (SAR) data to assess the correlation between outcrops of ultramafic‑carbonatite intrusives and regional craton-scale fractures and assign targets for further exploration. Kernel density cluster, two-point azimuth, and the Fry method two-dimension spatial analyses are applied to 58 known outcrops of intrusives in the Avon Volcanic District (AVD), southeast Missouri USA. The results of the spatial analysis show linear trends among the outcrops with three primary orientations: northwest, northeast, and east-west and a strong correlation with regional fault zones. SAR data processing and analysis guided by the spatial results reveals rhombic lineament patterns with sides that parallel the northwest-northeast emplacement patterns of the AVD. These rhombic structures are related to possible extensional stepover duplexes, likely formed during Precambrian North American plate rifting along transfer faults. Reactivation of the Precambrian structures during high angle dip slip and wrench faulting that was penecontemporaneous with emplacement of the AVD in the Devonian likely provided pathways for the mantle-derived intrusives. These findings can guide exploration efforts and help us understand the emplacement of a unique intracontinental intrusive complex.

The end of Midcontinent Rift magmatism and the paleogeography of Laurentia

Paleomagnetism of the North American Midcontinent Rift provides a robust paleogeographic record of Laurentia (cratonic North America) from ca. 1110 to 1070 Ma, revealing rapid equatorward motion of the continent throughout rift magmatism. Existing age and paleomagnetic constraints on the youngest rift volcanic and sedimentary rocks have been interpreted to record a slowdown of this motion as rifting waned. We present new paleomagnetic and geochronologic data from the ca. 1090–1083 Ma “late-stage” rift volcanic rocks exposed as the Lake Shore Traps (Michigan), the Schroeder-Lutsen basalts (Minnesota), and the Michipicoten Island Formation (Ontario). The paleomagnetic data allow for the development of paleomagnetic poles for the Schroeder-Lutsen basalts (187.8°E, 27.1°N; A 95 = 3.0°, N = 50) and the Michipicoten Island Formation (174.7°E, 17.0°N; A 95 = 4.4°, N = 23). Temporal constraints on late-stage paleomagnetic poles are provided by high-precision, 206 Pb- 238 U zircon dates from a Lake Shore Traps andesite (1085.57 ± 0.25 Ma; 2s internal errors), a Michipicoten Island Formation tuff (1084.35 ± 0.20 Ma) and rhyolite (1083.52 ± 0.23 Ma), and a Silver Bay aplitic dike from the Beaver Bay Complex (1091.61 ± 0.14 Ma), which is overlain by the Schroeder-Lutsen basalt flows. These Michipicoten Island Formation dates are the youngest yet obtained from Midcontinent Rift volcanic rocks and indicate that rift magmatism was active for at least 25 m.y. The addition of these late-stage paleomagnetic poles to the Laurentian apparent polar wander path suggests that rapid motion of Laurentia continued throughout the entirety of rift volcanism.

Distinct crustal structure of the North American Midcontinent Rift from P wave receiver functions

AbstractEighty‐two broadband seismic stations of the Superior Province Rifting Earthscope Experiment (SPREE) collected 2.5 years of continuous seismic data in the area of the high gravity anomaly associated with the Midcontinent Rift (MCR). Over 100 high‐quality teleseismic earthquakes were used for crustal P wave receiver function analysis. Our analysis reveals that the base of the sedimentary layer is shallow outside the MCR, thickens near the flanks where gravity anomalies are low, and shallows again in the MCR's center where the gravity anomalies peak. This pattern is similar to that found from local geophysical studies and is consistent with reverse faulting having accompanied the cessation of rifting at 1.1 Ga. Intermittent intracrustal boundaries imaged by our analysis might represent the bottom of the MCR's mostly buried dense volcanic layers. Outside the MCR, the Moho is strong, sharp, and relatively flat, both beneath the Archean Superior Province and the Proterozoic terranes to its south. Inside the MCR, two weaker candidate Mohos are found at depths up to 25 km apart in the rift's center. The intermediate layer between these discontinuities tapers toward the edges of the MCR. The presence of this transitional layer is remarkably consistent along the strike of the MCR, including beneath its jog in southern Minnesota, near the Belle Plaine Fault. We interpret these results as evidence for extensive underplating as a defining characteristic of the rift, which remains continuous along the Minnesota jog, where due to its orientation, it is minimally affected by the reverse faulting that characterizes the NNE striking parts of the rift.

Nanoscale constraints on porosity generation and fluid flow during serpentinization

Research Article| February 01, 2016 Nanoscale constraints on porosity generation and fluid flow during serpentinization Benjamin M. Tutolo; Benjamin M. Tutolo * 1Department of Earth Sciences, University of Oxford, Oxford OX1 3AN, UK2Department of Earth Sciences, University of Minnesota, Minneapolis, Minnesota 55455, USA *E-mail: benjamin.tutolo@earth.ox.ac.uk Search for other works by this author on: GSW Google Scholar David F.R. Mildner; David F.R. Mildner 3National Institute of Standards and Technology Center for Neutron Research, Gaithersburg, Maryland 20899, USA Search for other works by this author on: GSW Google Scholar Cedric V.L. Gagnon; Cedric V.L. Gagnon 3National Institute of Standards and Technology Center for Neutron Research, Gaithersburg, Maryland 20899, USA4Department of Materials Science and Engineering, University of Maryland, College Park, Maryland 20742, USA Search for other works by this author on: GSW Google Scholar Martin O. Saar; Martin O. Saar 2Department of Earth Sciences, University of Minnesota, Minneapolis, Minnesota 55455, USA5Department of Earth Sciences, ETH-Zürich, 8006 Zurich, Switzerland Search for other works by this author on: GSW Google Scholar William E. Seyfried, Jr. William E. Seyfried, Jr. 2Department of Earth Sciences, University of Minnesota, Minneapolis, Minnesota 55455, USA Search for other works by this author on: GSW Google Scholar Author and Article Information Benjamin M. Tutolo * 1Department of Earth Sciences, University of Oxford, Oxford OX1 3AN, UK2Department of Earth Sciences, University of Minnesota, Minneapolis, Minnesota 55455, USA David F.R. Mildner 3National Institute of Standards and Technology Center for Neutron Research, Gaithersburg, Maryland 20899, USA Cedric V.L. Gagnon 3National Institute of Standards and Technology Center for Neutron Research, Gaithersburg, Maryland 20899, USA4Department of Materials Science and Engineering, University of Maryland, College Park, Maryland 20742, USA Martin O. Saar 2Department of Earth Sciences, University of Minnesota, Minneapolis, Minnesota 55455, USA5Department of Earth Sciences, ETH-Zürich, 8006 Zurich, Switzerland William E. Seyfried, Jr. 2Department of Earth Sciences, University of Minnesota, Minneapolis, Minnesota 55455, USA *E-mail: benjamin.tutolo@earth.ox.ac.uk Publisher: Geological Society of America Received: 10 Sep 2015 Revision Received: 25 Nov 2015 Accepted: 30 Nov 2015 First Online: 09 Mar 2017 Online Issn: 1943-2682 Print Issn: 0091-7613 © 2016 Geological Society of America Geology (2016) 44 (2): 103–106. https://doi.org/10.1130/G37349.1 Article history Received: 10 Sep 2015 Revision Received: 25 Nov 2015 Accepted: 30 Nov 2015 First Online: 09 Mar 2017 Cite View This Citation Add to Citation Manager Share Icon Share Facebook Twitter LinkedIn MailTo Tools Icon Tools Get Permissions Search Site Citation Benjamin M. Tutolo, David F.R. Mildner, Cedric V.L. Gagnon, Martin O. Saar, William E. Seyfried; Nanoscale constraints on porosity generation and fluid flow during serpentinization. Geology 2016;; 44 (2): 103–106. doi: https://doi.org/10.1130/G37349.1 Download citation file: Ris (Zotero) Refmanager EasyBib Bookends Mendeley Papers EndNote RefWorks BibTex toolbar search Search Dropdown Menu toolbar search search input Search input auto suggest filter your search All ContentBy SocietyGeology Search Advanced Search Abstract Field samples of olivine-rich rocks are nearly always serpentinized—commonly to completion—but, paradoxically, their intrinsic porosity and permeability are diminishingly low. Serpentinization reactions occur through a coupled process of fluid infiltration, volumetric expansion, and reaction-driven fracturing. Pores and reactive surface area generated during this process are the primary pathways for fluid infiltration into and reaction with serpentinizing rocks, but the size and distribution of these pores and surface area have not yet been described. Here, we utilize neutron scattering techniques to present the first measurements of the evolution of pore size and specific surface area distribution in partially serpentinized rocks. Samples were obtained from the ca. 2 Ma Atlantis Massif oceanic core complex located off-axis of the Mid-Atlantic Ridge and an olivine-rich outcrop of the ca. 1.1 Ga Duluth Complex of the North American Mid-Continent Rift. Our measurements and analyses demonstrate that serpentine and accessory phases form with their own, inherent porosity, which accommodates the bulk of diffusive fluid flow during serpentinization and thereby permits continued serpentinization after voluminous serpentine minerals fill reaction-generated porosity. You do not have access to this content, please speak to your institutional administrator if you feel you should have access.

Geochemical Stratigraphy of the Keweenawan Midcontinent Rift Volcanic Rocks with Regional Implications for the Genesis of Associated Ni, Cu, Co, and Platinum Group Element Sulfide Mineralization

The development of the North American Midcontinent rift can be understood in the context of a tectonic and magmatic event in which the far-field effects of continental drift influenced magmatism produced by a mantle plume. We report and interpret new data for samples collected from stratigraphically controlled sequences of basaltic rocks in the Osler Volcanic Group and the Mamainse Point Volcanic Group in Ontario, Canada. We confirm that the earliest phase of rift development involved primitive melts that produced basalts in the Lower Formation of the Osler Volcanic Group with Fo74–80 olivine (1,490–2,280 ppm Ni). The earliest basalts are primitive in both the Mamainse Point Volcanic Group and Osler Volcanic Group sections, with elevated MgO, Cr, and Ni; these rocks also have high TiO2 and high Gd/Yb, which are features that indicate low degrees of melting of a garnet-bearing lherzolite mantle. The Central Formation A and B basalts of the Osler Volcanic Group that overlie the Lower Formation lavas record an upward increase and then decrease in the degree of contamination of basaltic rocks with normal Ni and platinum group element (PGE) abundance levels that require no sulfide saturation event to accompany contamination. The Upper Formation of the Osler Volcanic Group is characterized by a limited number of PGE-depleted flows; although these basalts are highly fractionated and typically develop textures with stellate to laminated plagioclase feldspar, they are usually not significantly contaminated. The Lower and Upper formations of the Mamainse Point Volcanic Group show a diversity in geochemistry consistent with contributions from different mantle sources, but contamination is restricted to a small number of flows with normal Ni and PGE abundance levels. There is no indication of a regional metal depletion signature on the scale of the Nadezhdinsky Formation at Noril’sk, but there is evidence that local S saturation events were triggered in the unexposed conduits feeding some of the upper sequence flows in response to fractionation of the magma. Massive to heavy disseminated Ni-Cu-PGE sulfide mineralization in the Midcontinent rift is hosted by small primitive intrusions that were formed early in the evolution of the rift and are typically associated with transtensional structures in the cratonic blocks adjacent to the main rift. Examples include the Eagle deposit in Michigan, the Tamarack mineralization in Minnesota, and the Current Lake Complex in Ontario. The American examples are associated with splays of the Great Lakes tectonic zone and the Canadian example is associated with the Quetico fault zone. The intrusions share the geochemical fingerprints of the more primitive lowermost volcanic rocks of the Midcontinent rift and they are broadly contemporaneous with the development of the lower stratigraphy of the Osler Volcanic Group and Mamainse Point Volcanic Group. Although the low-grade but large-tonnage Cu-Ni-PGE deposits of the Duluth were formed significantly later than the small primitive intrusions, their host rocks share the same lithophile trace element signature as these early intrusions. The parental magmas that produced all of these intrusions were derived from deeper mantle depths, where garnet was a stable phase, than the magmas that produced the Central Formation basalts and some of the Upper Formation basalts investigated in this study.

Evolution of the Main Zone at the Marathon Cu-PGE Sulfide Deposit, Midcontinent Rift, Canada: Spatial Relationships in a Magma Conduit Setting

Various ultramafic Ni-Cu-platinum group element (PGE) deposits associated with the North American Midcontinent rift have been attributed to formation in a magma conduit setting, whereby PGE concentration is controlled by various fluid dynamic processes. The Marathon Cu-PGE sulfide deposit located within the Midcontinent rift-related Coldwell Alkaline Complex has been classified as a gabbro-associated contact-type deposit; however, both magmatic and hydrothermal processes have been proposed to account for the significant concentration of PGE. In light of the growing field of evidence for magma conduit-type settings, this study comprised a comprehensive geochemical investigation of the complicated crosscutting gabbroic to ultramafic units in the immediate vicinity of the Marathon deposit; and a thorough three-dimensional investigation of the distribution of Cu and Pd within the Main mineralized zone. The main objectives of this study were to test the applicability of the magma conduit deposit model to the Marathon deposit and to identify key exploration criteria for use elsewhere in the Coldwell Alkaline Complex. Mineralization is hosted by the Two Duck Lake gabbro, a 4-km-long and 250-m-thick unit of the Marathon Series. The Marathon Series is the latest of three magmatic series that make up the 1- to 2-km-thick Eastern Gabbro Suite, which wraps around the eastern and northern margin of the Coldwell Alkaline Complex. The three magmatic series are shown here to have distinct trace element signatures that enable reliable discrimination of potentially sulfide and PGE-bearing units of the Marathon Series from the barren rocks of either the Fine-Grained or Layered Series. At the Marathon deposit, sulfides consist of disseminated chalcopyrite, pyrrhotite, and minor bornite and occur within the Main, Footwall, and Hanging-wall zones and in the PGE-enriched W Horizon. This paper focused on sulfides located within the Main zone, including the keel-shaped feeder channel that continues downdip to over 550-m depth. The spatial distribution of Cu, Pd, and Cu/Pd were examined in relation to a three-dimensional surface model for the footwall contact; in a vertical profile through the Main zone; and in a longitudinal section that cuts the feeder channel. There are three important observations: (1) trends for elevated Cu and Pd are parallel to numerous troughs and ridges in the footwall, (2) Cu, Pd, and Cu/Pd varies up section in a saw-toothed pattern from high to low values, and (3) the proportion of high Cu/Pd sulfides is greatest within the thickest accumulations of sulfides within the feeder channel. Evaluation of interelement relationships between Cu and Pd and between Pd and Ir, Rh, Pt, and Au for mineralization within the Main zone indicate positive associative, but nonlinear behavior for all elements. Briefly, the data show nonlinear correlations between Cu and Pd in which Cu/Pd decreases with increasing Pd; and coherent but nonlinear behavior for Ir, Rh, Pt, and Pd in which Pd/Ir, Pd/Rh, Pd/Pt, and Pd/Au all increase with increasing Pd. The observed variation in Cu/Pd is consistent with a magmatic model calculated by others for deposits in the Duluth Complex, in which sulfides accumulated in a closed system from a melt with mantlelike Cu/Pd and an elevated silicate to sulfide ratio. The observed variations in Pd/Ir, Pd/Rh, and Pd/Pt are consistent with R factor fractionation related to differences in the relative partition coefficients between sulfide and silicate melts, and rule out the possibility that processes such as fractionation of sulfide melt by monosulfide solid solution (mss) or redistribution of metals during hydrothermal alteration played a significant role in the mineralizing event. The Two Duck Lake gabbro and associated sulfides of the Marathon deposit are proposed to have formed by multiple injections of plagioclase crystal mush that carried droplets of sulfide liquid along a conduit system that was controlled by radial and ring fault structures in the Coldwell Alkaline Complex. The accumulation of sulfides was controlled by flow dynamic processes within the magma channels, but Cu/Pd was controlled by local proportions of silicate melt to sulfide liquid. Key characteristics of the deposit that are critical to exploration elsewhere in the Coldwell Alkaline Complex include the following: (1) the recognition that gabbroic to ultramafic intrusions of the Marathon Series are the host for Cu and PGE mineralization, (2) the distribution of Cu/Pd data within sulfide occurrences are useful as vectors toward the feeder channel, (3) topographic lineaments are indicators of potential mineralized feeder zones, and (4) oxide- and apatite-rich, irregularly shaped gabbroic to ultramafic pods are potential indicators of an underlying feeder channel.