Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Ice Temperature Profiles

ESA
import xarray as xr
import cartopy.crs as ccrs
import hvplot.xarray

url = 'https://s3.waw4-1.cloudferro.com/EarthCODE/OSCAssets/polar_cube_datasets/ice_sheet_temperature/SM_TEST_MIR_ITUDP4_20130101T000000_20141231T000000_200_001_0.zarr/'
ds = xr.open_zarr(url)
ds
Loading...
Loading...
Loading...
Loading...
Loading...
lon = ((ds["longitude"] + 180) % 360) - 180

tice = ds["Tice"].assign_coords(longitude=lon, latitude=ds["latitude"])


tice = tice.rename("Tice (Kelvin)")

tice.hvplot.quadmesh(
    x="longitude",
    y="latitude",
    groupby="depth",
    geo=True,
    crs=ccrs.PlateCarree(),
    projection=ccrs.SouthPolarStereo(central_longitude=0, true_scale_latitude=-71),
    coastline=True,
    cmap="coolwarm",
    colorbar=True,
    clabel=tice.name,
    width=850,
    height=700,
    title="Ice temperature by depth",
)
Loading...
References
  1. Davison, B. J., Hogg, A. E., Gourmelen, N., Jakob, L., Wuite, J., Nagler, T., Greene, C. A., Andreasen, J., & Engdahl, M. E. (2023). Annual mass budget of Antarctic ice shelves from 1997 to 2021. Science Advances, 9(41). 10.1126/sciadv.adi0186