Comparisons with CryoSat-2 and PIOMAS sea ice thickness estimates
Comparisons with CryoSat-2 and PIOMAS sea ice thickness estimates#
Summary: In this notebook, we compare our ICESat-2-derived sea ice thickness estimates with various estimates derived from ESA’s CryoSat-2 radar altimeter, a combined ICESat-2/CryoSat-2 product and PIOMAS. See the CryoSat-2 wrangling notebook for more insight into the different CryoSat-2 products.
Doesn’t include spatial comparison stats, but would be pretty trivial to add that in - working on it!
Also lots more to be done in terms of assessing consistency between CS-2 and IS-2 and where/why that breaks down.
Version history: Version 1 (01/01/2022)
Import notebook dependencies#
import xarray as xr import pandas as pd import numpy as np import itertools import pyproj from netCDF4 import Dataset import scipy.interpolate # Helper function for reading the data from the bucket from utils.read_data_utils import read_book_data from utils.plotting_utils import compute_gridcell_winter_means, interactiveArcticMaps, interactive_winter_mean_maps, interactive_winter_comparison_lineplot # Plotting # Plotting dependencies import cartopy.crs as ccrs from textwrap import wrap import hvplot.xarray import holoviews as hv import matplotlib.pyplot as plt from matplotlib.axes import Axes from cartopy.mpl.geoaxes import GeoAxes GeoAxes._pcolormesh_patched = Axes.pcolormesh # Helps avoid some weird issues with the polar projection %config InlineBackend.figure_format = 'retina' import matplotlib as mpl mpl.rcParams['figure.dpi'] = 150 # Sets figure size in the notebook # Remove warnings to improve display import warnings warnings.filterwarnings('ignore')