Magnetotelluric evidence for the formation of the layered Sask Craton by flat slab subduction

Long-period magnetotelluric (MT) data were collected at 56 locations over the Sask Craton in 2021 and 2022. The data were combined with existing broadband data and inverted to produce a 3-D resistivity model of the Sask Craton and Trans-Hudson Orogen (THO). The model reveals a number of northeast striking electrically conductive crustal structures that extend into the mantle lithosphere. In the mantle lithosphere, these conductors coalesce into a single large low resistivity anomaly in the depth range 70-85 km termed the Northern Sask Craton (NSC) conductor. The resistivity of the NSC conductor is attributed to sulfides deposited along an interface between a flat slab that was accreted to the base of the pre-THO Sask Craton lithosphere during closure of the Manikewan Ocean. Kimberlites have erupted along the margin of the NSC conductor. The boundary of the conductor likely represents deep-seated faults and mantle terrane boundaries formed during flat slab subduction that allowed the ascent of kimberlite melts. The resistivity of the northeast-trending conductors can be interpreted as due to graphite and sulfides precipitated by past fluid or melt flow during ocean closure and orogensis. A number of these conductors are located beneath known mineral districts and trends and may represent source pathways for regional base and precious metal deposits. Other conductors may represent possible, previously unknown, regions hosting mineralization. Many of these conductors are associated with major regional faults and shear zones, which may be deep-seated features that helped to guide both kimberlites and mineralizing fluids. Of the prominent northeast-trending conductors west of the Sask Craton, one corresponds to the previously reported North American Central Plains (NACP) conductor. The new model shows that this conductor abruptly terminates at 54 degrees N and is not observed farther south in the model. This shows that the NACP is not as spatially continuous as previously suggested, suggesting that the tectonic processes that formed the THO were not as uniform alongstrike as shown in existing tectonic models. The connection of one of the northeast-trending anomalies to the NSC conductor suggests that a previously unrecognized phase of east-dipping subduction may have occurred beneath the Sask Craton as the THO was formed.