Abstract We identify and discuss 57 magnetic anomaly pattern domains spanning the Circum-Arctic. The domains are based on analysis of a new Circum-Arctic data compilation. The magnetic anomaly patterns can be broadly related to general geodynamic classification of the crust into stable, deformed (magnetic and nonmagnetic), deep magnetic high, oceanic and large igneous province domains. We compare the magnetic domains with topography/bathymetry, regional geology, regional free air gravity anomalies and estimates of the relative magnetic ‘thickness’ of the crust. Most of the domains and their geodynamic classification assignments are consistent with their topographic/bathymetric and geological expression. A few of the domains are potentially controversial. For example, the extent of the Iceland Faroe large igneous province as identified by magnetic anomalies may disagree with other definitions for this feature. Also the lack of definitive magnetic expression of oceanic crust in Baffin Bay, the Norwegian–Greenland Sea and the Amerasian Basin is at odds with some previous interpretations. The magnetic domains and their boundaries provide clues for tectonic models and boundaries within this poorly understood portion of the globe.

A new isostatic residual gravity map of the conterminous United States presents continent-wide gravity data in a form that can be readily used, with geologic information and other geophysical data, in studies of the composition and structure of the continental crust. This map was produced from the gridded gravity data on which the recently released Gravity Anomaly Map of the United States is based. About 1 million onland and 0.8 million offshore gravity observations interpolated to a 4- by 4-km grid serve as the basis for both maps. The Airy-Heiskanen model of isostatic compensation of topography applied to topographic and bathymetric data averaged over 5- by 5-min compartments was used to remove, to first order, the large, long-wavelength Bouguer gravity anomalies caused by deep density distributions that support topographic loads. The parameters used in the Airy-Heiskanen model were topographic density, 2.67 g/cm 3 ; sea-level crustal thickness, 30 km; and density contrast across the base of the model crust, 0.35 g/cm 3 . Many of the conspicuous short-wavelength anomalies (widths less than several hundred kilometers) on the isostatic residual gravity map correlate with mapped or near-surface geologic features, and primarily reflect shallow-density distributions rather than any departures from isostatic equilibrium. In general, gravity highs occur over (1) mafic igneous bodies emplaced in rift or magmatic arc settings or as isolated intrusions controlled by structures; (2) accreted slices of mafic oceanic, island-arc, or transitional crust; and (3) uplifted crystalline basement. Gravity lows are found over (1) thick bodies of felsic intrusive or extrusive rocks; (2) sedimentary deposits in extensional, convergent, or transform settings; and (3) depressed crystalline basement. Anomalies with widths as much as 1,000 km or more also appear to reflect crustal properties in many cases—several broad gravity highs are associated with crust having a high average seismic wave velocity, and comparable broad gravity lows occur over areas of low average seismic velocity. Alternative ways of viewing the isostatic residual gravity data provide additional information about density distributions in the crust. The first-vertical derivative map accentuates gravity anomalies over shallow, abrupt density changes at the expense of those resulting from deeper or more gradual density transitions. The maximum horizontal gradient map contains information about the locations of pronounced density boundaries. Two-dimensional spectral analysis of the gravity data provides a quantitative means for identifying dominant fabrics in the gravity field and for distinguishing various terranes from each other. Neither Bouguer nor isostatic residual gravity anomalies are particularly well suited for practical modeling of deep structure in conjunction with deep seismic information. However, a scheme in which the entire Earth outside the area of interest is approximated by laterally homogeneous layers and isostatically compensated topography, and in which the area of interest is modeled using the seismic constraints applied in a two-and-one-half-dimensional geometry, holds promise for exploiting useful features of both the Bouguer and isostatic residual gravity anomalies.