Ice Velocity

Description: Gridded estimates of ice velocity from the MEaSUREs Phase-Based Greenland Ice Velocity Map.

Requires an EarthData login, see link for how to configure this.

Dataset:

National Snow and Ice Data Center. “MEaSUREs Multi-Year Greenland Ice Sheet Velocity Mosaic, Version 1,” October 21, 2020. https://nsidc.org/data/nsidc-0670/versions/1.

Associated paper:

N.A.

%%capture
%load_ext autoreload
%autoreload 2

from polartoolkit import fetch, maps, utils

data = fetch.ice_vel(
    hemisphere="north",
    # available options
    # region,
    # spacing,
    # registration ("g" for gridline or "p" for pixel),
)
data

<xarray.DataArray (y: 10920, x: 6020)> Size: 263MB
array([[nan, nan, nan, ..., nan, nan, nan],
       [nan, nan, nan, ..., nan, nan, nan],
       [nan, nan, nan, ..., nan, nan, nan],
       ...,
       [nan, nan, nan, ..., nan, nan, nan],
       [nan, nan, nan, ..., nan, nan, nan],
       [nan, nan, nan, ..., nan, nan, nan]], dtype=float32)
Coordinates:
  * x        (x) float64 48kB -6.45e+05 -6.448e+05 ... 8.595e+05 8.598e+05
  * y        (y) float64 87kB -6.402e+05 -6.405e+05 ... -3.37e+06 -3.37e+06

xarray.DataArray

• y: 10920
• x: 6020

• nan nan nan nan nan nan nan nan ... nan nan nan nan nan nan nan nan

  array([[nan, nan, nan, ..., nan, nan, nan],
         [nan, nan, nan, ..., nan, nan, nan],
         [nan, nan, nan, ..., nan, nan, nan],
         ...,
         [nan, nan, nan, ..., nan, nan, nan],
         [nan, nan, nan, ..., nan, nan, nan],
         [nan, nan, nan, ..., nan, nan, nan]], dtype=float32)

• Coordinates: (2)
  • x
    (x)
    float64
    -6.45e+05 -6.448e+05 ... 8.598e+05

    array([-645000., -644750., -644500., ...,  859250.,  859500.,  859750.])

  • y
    (y)
    float64
    -6.402e+05 -6.405e+05 ... -3.37e+06

    array([ -640250.,  -640500.,  -640750., ..., -3369500., -3369750., -3370000.])

• Indexes: (2)
  • x
    PandasIndex

    PandasIndex(Index([-645000.0, -644750.0, -644500.0, -644250.0, -644000.0, -643750.0,
           -643500.0, -643250.0, -643000.0, -642750.0,
           ...
            857500.0,  857750.0,  858000.0,  858250.0,  858500.0,  858750.0,
            859000.0,  859250.0,  859500.0,  859750.0],
          dtype='float64', name='x', length=6020))

  • y
    PandasIndex

    PandasIndex(Index([ -640250.0,  -640500.0,  -640750.0,  -641000.0,  -641250.0,  -641500.0,
            -641750.0,  -642000.0,  -642250.0,  -642500.0,
           ...
           -3367750.0, -3368000.0, -3368250.0, -3368500.0, -3368750.0, -3369000.0,
           -3369250.0, -3369500.0, -3369750.0, -3370000.0],
          dtype='float64', name='y', length=10920))

• Attributes: (0)

_ = utils.get_grid_info(data, print_info=True)

grid spacing: 250.0 m
grid region: (-645000.0, 859750.0, -3370000.0, -640250.0)
grid zmin: 0.00136277475394
grid zmax: 12906.4941406
grid registration: g

fig = maps.plot_grd(
    data,
    title="Ice Velocity",
    cmap="plasma",
    coast=True,
    cbar_label="meters/year",
    hist=True,
    robust=True,
    grd2cpt=True,
    hemisphere="north",
)

fig.show(dpi=200)

/home/mattd/polartoolkit/src/polartoolkit/maps.py:1626: UserWarning: getting max/min values from grid/points, if cpt_lims were used to create the colorscale, histogram will not properly align with colorbar!
  add_colorbar(