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.

import polartoolkit as ptk

data = ptk.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]],
      shape=(10920, 6020), 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]],
        shape=(10920, 6020), dtype=float32)

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

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

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

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

• 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)

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

grid spacing: 250.0 m
grid region: (-645125.0, 859875.0, -3369875.0, -640375.0)
grid zmin: 0.00136277475394
grid zmax: 12906.4941406
grid registration: p

fig = ptk.plot_grid(
    data,
    title="Ice Velocity",
    cmap="magma",
    cpt_lims=(0.001, 3),
    cpt_log=True,
    cbar_log=True,
    coast=True,
    cbar_label="meters/year",
    hemisphere="north",
)

fig.show(dpi=200)