[1] A global model is developed for both induced and remanent magnetizations in the terrestrial lithosphere. The model is compared with, and well-described by, Ørsted satellite observations. Interpretation of the observations over North America suggests that the large total field anomalies, associated with spherical harmonic degrees 15 -26 and centered over Kentucky and the south-central United States, are the manifestations of the magnetic edges of the southern boundaries of cratonic North America. The techniques and models developed here may be of use in defining other cratonic boundaries.INDEX TERMS: 1545 Geomagnetism and Paleomagnetism: Spatial variations (all harmonics and anomalies); 1219 Geodesy and Gravity: Local gravity anomalies and crustal structure; 7218 Seismology: Lithosphere and upper mantle; 9350 Information Related to Geographic Region: North America
Background[2] The launch of the Ørsted high-precision geomagnetic field satellite [Neubert et al., 2001] has invigorated efforts to understand the magnetic field of the earth's lithosphere. Early attempts [Langel and Hinze, 1998] to model the lithospheric field relied on forward and inverse approaches over local regions. After the realization that much of the lithospheric magnetic signal might be obscured by overlap with the long-wavelength magnetic field from the core [Meyer et al., 1985], recent work has explored the potential of global forward [Cohen and Achache, 1994; Dyment and ArkaniHamed, 1998] and inverse models [Purucker et al., 1998]. Our work elaborates on but differs from previous work which 1) included remanent magnetizations associated only with the Cretaceous Quiet Zones [Cohen and Achache, 1994;Purucker et al., 1998], 2) did not consider induced magnetizations at all [Dyment and Arkani-Hamed, 1998], 3) used a model of induced magnetization that had fewer geologic and geophysical inputs [Purucker et al., 1998].
A New Global Magnetization Model[3] The global model of induced magnetization is based on an estimate of the volume of the magnetic crust and its magnetic susceptibility. We assume that induced magnetizations are restricted to the crust [Wasilewski and Mayhew, 1992] and utilize a global seismic tomography model [Nataf and Ricard, 1996] for estimating crustal thickness. The model also contains a tectonicbased subdivision of the crust into three categories, each of which has an associated geotherm. These geotherms, when coupled with an assumption about the magnetic mineral(s) responsible for the bulk of the magnetization, allow for the calculation of a depth to the Curie isotherm. We assume here that the magnetic mineral is magnetite or low-Ti magnetite. The magnetic layer thickness is calculated as the thickness of the igneous crust above the magnetite Curie isotherm. We utilize a sediment thickness model [Laske and Masters, 1997] to account for the presence of effectively non-magnetic sediment or sedimentary rock which serves to decrease the effective magnetic layer thickness. Although we calculate our induced magnetization mod...