Geologic information

Description: Download fault locations, geologic unit boundaries, and source information for Antarctic from the GeoMAP compilation.

Dataset:

Cox, Simon Christopher, Belinda Smith Lyttle, Samuel Elkind, Christine S. Smith Siddoway, Paul Morin, Giovanni Capponi, Tamer Abu-Alam, et al. “The GeoMAP (v.2022-08) Continent-Wide Detailed Geological Dataset of Antarctica.” PANGAEA, 2023. https://doi.org/10.1594/PANGAEA.951482.

Associated paper:

Cox, Simon C., Belinda Smith Lyttle, Samuel Elkind, Christine Smith Siddoway, Paul Morin, Giovanni Capponi, Tamer Abu-Alam, et al. “A Continent-Wide Detailed Geological Map Dataset of Antarctica.” Scientific Data 10, no. 1 (May 18, 2023): 250. https://doi.org/10.1038/s41597-023-02152-9.

%%capture
%load_ext autoreload
%autoreload 2

from polartoolkit import fetch, maps, regions

version_names = [
    "faults",
    "units",
    # "sources", # not interesting to plot
    # "quality", # not interesting to plot
]


data = []
for name in version_names:
    d = fetch.geomap(
        version=name,
        # available options
        # region,
        region=regions.antarctica,
    )
    # only keep subset of largest units to speed up the plotting below
    if name == "units":
        d = d.sort_values(by="Shape_Area", ascending=False).copy()
        d = d.iloc[0 : int(len(d) / 4)]  # largest 1/4 of shapes
    data.append(d)

fig = maps.basemap(
    region=regions.minna_bluff,
    hemisphere="south",
    coast=True,
    imagery_basemap=True,
)
fig.add_geologic_units()
fig.add_faults()
fig.show()

Plotting faults

You can also use the function maps.add_faults for an easy way to specify and plot faults.

fig = maps.basemap(
    region=regions.antarctica,
    hemisphere="south",
    coast=True,
)

# with defaults
fig.add_faults()

fig.show(dpi=200)

fig = maps.basemap(
    region=regions.antarctica,
    hemisphere="south",
    coast=True,
)

# specifying fault types
fig.add_faults(
    faults_activity="inactive",
    faults_motion="normal",
    faults_exposure="inferred",
    pen="2p,red",
)

fig.show(dpi=200)