3-D Modeling of the Duluth Complex from geophysical data

By: , and 



The Mesoproterozoic Duluth Complex in northeastern Minnesota is one of the major plutonic components of the Midcontinent Rift System and hosts a variety of copper-nickel sulfide and platinum-group element deposits. The Duluth Complex is composed of a series of individual mafic and felsic intrusions emplaced 1110-1098 Ma within Paleoproterozoic sedimentary rocks of the Animikie basin and volcanic flows of the Midcontinent Rift. Prior work has included 2-D modeling and qualitative geologic interpretations of gravity and magnetic data (e.g., Chandler, 1990; Chandler and Ferderer, 1989), much of which is still preliminary (V. Chandler, written commun., 2020). Three-dimensional modeling has been limited, with only one 3-D model created using Bouguer gravity data constrained by seismic-reflection interpretations as part of a PhD thesis (Allen, 1994). Given the complex geology of the area, 3-D modeling is useful for providing a complete picture of the variable densities, susceptibilities, and electrical resistivities throughout the Duluth Complex and associated volcanic rocks as well as their depth extent beneath sedimentary cover. Models of these geophysical properties at depth enable more accurate geologic mapping in the subsurface which can lead to an improved understanding of the formation history of the Duluth Complex. In this study, we use aeromagnetic data acquired between 1979-1991 (Chandler, 2007), Bouguer gravity data collected since 1950 (Chandler and Lively, 2019), and magnetotelluric data collected in 2019 to create new 2-D and 3-D geophysical models of the Duluth Complex constrained by seismic reflection, geologic, and rock property data. An inversion of the Bouguer gravity data for thickness of the Duluth Complex using constant densities of 3110 kg/m3 and 2670 kg/m3 for the Duluth Complex and surrounding crustal rocks, respectively, results in thicknesses ranging from ~3-28 km for the Duluth Complex and related intrusions and volcanic rocks (Figure 1A). A 3-D model of the magnetotelluric data reveals low resistivity anomalies at ~5-10 km depth below the northern margin of the Duluth Complex and below the Greenwood Lake intrusion (Figure 1B). We expect to encounter low resistivities at depth associated with the Paleoproterozoic Animikie basin, which makes up the floor of the Duluth Complex, and therefore interpret these anomalies as either the base of the complex or as fragments of Animikie sediments interfingered with igneous intrusive rocks. Finally, 3-D voxel models of density and susceptibility illuminate the subsurface distribution of rock properties below the Duluth Complex which, in combination with resistivity and thickness models, can be used to create a 3-D geologic map of this area.

Study Area

Publication type Conference Paper
Publication Subtype Conference Paper
Title 3-D Modeling of the Duluth Complex from geophysical data
Year Published 2021
Language English
Publisher Institute on Lake Superior Geology
Contributing office(s) Geology, Geophysics, and Geochemistry Science Center
Description 2 p.
Larger Work Type Book
Larger Work Subtype Conference publication
Larger Work Title 67th Institute on Lake Superior Geology Proceedings
First page 52
Last page 53
Country United States
State Minnesota
Other Geospatial Duluth Complex
Google Analytic Metrics Metrics page
Additional publication details