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.

Notes:

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