Moho Depth and Bulk Crustal Properties in Northern Quebec and Labrador

Northern Quebec and Labrador lie at the heart of the Laurentian landmass and preserve over 3 billion years of continental evolution. In this region the Archean Superior and Nain cratons are surrounded by Paleoproterozoic orogens such as New-Quebec, Trans-Hudson and Torngat, as well as the younger Grenville orogen to the SE. Study of crustal structure in this region provides valuable information on the assembly of the North American continent. We use data from 8 seismic stations installed in summer 2011 as part of the QUiLLE (Quebec-Labrador Lithospheric Experiment) project to investigate crustal structure, using receiver function analysis. The data set covers 5 years (2011-2016) for most of the stations, comprising several hundred events of magnitude ≥5 and epicentral distance 30-90°. After initial data processing and quality control, several tens of events per station were used in an H-κ stacking analysis to estimate Moho depth and bulk crustal properties. Some stations show significant complexity in their receiver functions, leading to inconclusive H-κ results, but the majority show a consistent Moho signal from which crustal parameters are successfully extracted. Crustal thickness varies from 33 to 49 km, with the thickest crust associated with the Trans-Hudson orogen in the Ungava region of northernmost Quebec and the thinnest beneath the central Labrador coast. Vp/Vs ratios (κ) lie in the range 1.71-1.86, with the majority of values consistent with granite-gneiss-tonalite bulk crustal compositions. The receiver functions are combined with surface-wave group velocity data to model the crustal structures in more detail beneath each station, allowing us to investigate crustal layering, Moho complexity and lateral heterogeneity.