The crust surrounding SNO+

This web page provides the numerical 3D geophysical models of the crustal portion surrounding the SNO+ detector, a liquid scintillator detector located at SNOLAB underground physics laboratory (46.47° N, 81.20° W, Sudbury, Ontario, Canada).

The model of the LOcal Crust (LOC) belonging to the 6° x 4° (∼ 440 km x 460 km) region centered at SNOLAB is provided with a resolution of 1 km x 1 km x 0.1 km.

A high-resolution model (1 km x 1 km x 1 km) is available for the Close Upper Crust (CUC) corresponding to the 50 km x 50 km region centered at the detector.

Detailed information about these models are published in:

Huang, Y., V. Strati, F. Mantovani, S. B. Shirey and W. F. McDonough - Regional study of the Archean to Proterozoic crust at the Sudbury Neutrino Observatory (SNO+), Ontario: Predicting the geoneutrino flux. Geochemistry, Geophysics, Geosystems 15(10): 3925 - 3944 (2014). DOI: 10.1002/2014GC005397

Strati, V., S. A. Wipperfurth, M. Baldoncini, W. F. McDonough and F. Mantovani - Perceiving the Crust in 3-D: A Model Integrating Geological, Geochemical, and Geophysical Data. Geochemistry, Geophysics, Geosystems 18(12): 4326-4341 (2017). DOI: 10.1002/2017GC007067

The model of the LOC can be requested by email to radioactivity@fe.infn.it.

The model of the CUC can be downloaded at this link:

The Local Crust

The 3D geophysical model of the LOC was constructed from compiled geological and geophysical information. Crustal cross sections obtained from refraction and reflection seismic surveys were used to characterize the crust and assign uncertainties to its structure.

The Upper Crust was divided into seven dominant lithologic units on the basis of regional geological features (characteristic lithologies, metamorphism, tectonic events, U and Th abundances). Some interpreted crustal cross sections based upon seismic experiments were used to define the contacts between the units.


Lithology/GeologyThickness [km]Volume [106 km3]Volume [%]Density [g/cm3]Mass [1018 kg]
Tonalite/Tonalite gneiss (Wawa-Abitibi)16.62.5160.62.73 ± 0.086.9 ± 0.5
Central Gneiss Belt (Grenville Province)14.51.2530.22.73 ± 0.083.4 ± 0.2
(Meta)volcanic rocks (Abitibi subprovince)5.50.122.92.84 ± 0.140.34 ± 0.01
Paleozoic sedimentary rocks1.10.061.32.62 ± 0.190.16 ± 0.01
Granite or granodiorite (Wawa-Abitibi)5.20.092.22.67 ± 0.020.24 ± 0.00
Huronian Supergroup, Sudbury Basin4.40.112.72.69 ± 0.040.30 ± 0.01
Sudbury Igneous Complex6.10.0060.12.8 ± 0.10.02 ± 0.00

The Close Upper Crust

The 3D geophysical model of the CUC comprises 9 geological units defined considering the surface exposure described in the 1:250,000 scale Bedrock Geology of Ontario, conveniently simplified according to the spatial resolution of the available information about crustal structure.

The contacts between the units are defined combining multiple geological and geophysical inputs: (i) the contacts of the simplified geological map, (ii) the map of depth of the top of the middle crust of the LOC model, the 2.5-D geological modelsa along six profiles used for constructing the 3-D model reported in Olaniyan et al. (2015), and (v) five virtual cross sections derived from the LOC model. The surface topography for the CUC region uses the digital elevation modelb produced by the Shuttle Radar Topographic Mission.


LabelReservoirVolume [103 km3]Volume [%]Density [g/cm3]Mass [1015 kg]
GTGneissic Tonalite suite29.69 ± 1.4063.72.73 ± 0.0881.05 ± 5.01
HIHuronian Supergroup and minor Intrusions10.52 ± 0.4922.62.75 ± 0.0428.93 ± 1.79
NGNorite-gabbro, Sudbury Igneous Complex2.64 ± 0.125.72.83 ± 0.107.47 ± 0.46
GNGranophyre, Sudbury Igneous Complex1.43 ± 0.073.12.70 ± 0.103.86 ± 0.32
OPOnaping Fm., Whitewater Group0.94 ± 0.042.02.77 ± 0.042.60 ± 0.22
GFGrenville Front Tectonic Zone rocks0.83 ± 0.041.82.73 ± 0.082.27 ± 0.12
OWOnwatin Fm., Whitewater Group0.30 ± 0.010.62.68 ± 0.040.80 ± 0.05
CMChelmsford Fm., Whitewater Group0.23 ± 0.010.52.75 ± 0.040.62 ± 0.05
CTCartier Granite0.04 ± 0.0020.12.65 ± 0.020.11 ± 0.01

aOlaniyan, O., R. S. Smith and B. Lafrance (2015). "Regional 3D geophysical investigation of the Sudbury Structure." Interpretation 3(2): SL63-SL81.

bJarvis, A., H. Reuter, A. Nelson and E. Guevara (2008). "Hole-filled SRTM for the globe Version 4, available from the CGIAR-CSI SRTM 90m Database (http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp)."