Magnetics

Description: Gridded magnetic anomaly data from various sources.

Datasets:

ADMAP-1 data accessed from http://admap.kopri.re.kr/databases.html

Golynsky, Alexander V., Fausto Ferraccioli, Jong Kuk Hong, Dmitry A. Golynsky, Ralph R. B. von Frese, Duncan A. Young, Donald D. Blankenship, et al. “ADMAP2 Magnetic Anomaly Map of the Antarctic - Links to Grid (Grd) Files.” In Supplement to: Golynsky, AV et al. (2018): New Magnetic Anomaly Map of the Antarctic. Geophysical Research Letters, 45(13), 6437-6449, Https://Doi.Org/10.1029/2018GL078153. PANGAEA, 2018. https://doi.org/10.1594/PANGAEA.892723.

LCS-1: A high-resolution global model of the lithospheric magnetic field derived from CHAMP and Swarm satellite observations

Accessed from https://www.spacecenter.dk/files/magnetic-models/LCS-1/ version: scalar anomaly at Earth’s surface (WGS84 ellipsoid) for spherical harmonic degrees n=14-185

Associated papers:

Golynsky, Alexander, Massimo Chiappini, Detlef Damaske, Fausto Ferraccioli, Carol A. Finn, Takemi Ishihara, Hyung Rae Kim, et al. “ADMAP — A Digital Magnetic Anomaly Map of the Antarctic.” In Antarctica: Contributions to Global Earth Sciences, edited by Dieter Karl Fütterer, Detlef Damaske, Georg Kleinschmidt, Hubert Miller, and Franz Tessensohn, 109–16. Berlin, Heidelberg: Springer, 2006. https://doi.org/10.1007/3-540-32934-X_12.

Golynsky, A. V., F. Ferraccioli, J. K. Hong, D. A. Golynsky, R. R. B. von Frese, D. A. Young, D. D. Blankenship, et al. “New Magnetic Anomaly Map of the Antarctic: ADMAP2.” Geophysical Research Letters 45, no. 13 (July 16, 2018): 6437–49. https://doi.org/10.1029/2018GL078153.

Olsen, N., D. Ravat, C. C. Finlay, and L. K. Kother, LCS-1: A high-resolution global model of the lithospheric magnetic field derived from CHAMP and Swarm satellite observations, Geophys. J. Int., 211, 1461–1477, doi:10.1093/gji/ggx381, 2017

import polartoolkit as ptk

version_names = [
    "admap1",
    "admap2",
    # "admap2_gdb", # doesn't give grid so can't plot
    "LCS-1",
]

grids = []
for name in version_names:
    data = ptk.fetch.magnetics(
        version=name,
        # available options
        # region,
        # spacing,
        # registration ("g" for gridline or "p" for pixel),
        hemisphere="south",
    )
    grids.append(data)
    print(f"Info for {name}")
    _ = ptk.get_grid_info(data, print_info=True)
    print("##########")

Info for admap1
grid spacing: 5000.0 m
grid region: (-3330000.0, 3330000.0, -3330000.0, 3330000.0)
grid zmin: -936.875366211
grid zmax: 1766.1373291
grid registration: g
##########
Info for admap2
grid spacing: 1500.0 m
grid region: (-3423000.0, 3426000.0, -3424500.0, 3426000.0)
grid zmin: -2827.82373047
grid zmax: 4851.73828125
grid registration: g
##########
Info for LCS-1
grid spacing: 10000.0 m
grid region: (-3500000.0, 3500000.0, -3500000.0, 3500000.0)
grid zmin: -316.829925537
grid zmax: 601.757263184
grid registration: g
##########

fig = ptk.subplots(
    grids,
    region=ptk.regions.antarctica,
    fig_title="Magnetic Anomaly Data",
    titles=version_names,
    cbar_label="nT",
    margins=".5c/1.3c",
    coast=True,
    cmap="balance+h0",
    hemisphere="south",
    robust=True,
    hist=True,
)
fig.show(dpi=200)