David Gareth Babb

and 8 more

The loss of multiyear sea ice (MYI) in the Arctic Ocean is a significant change that affects all facets of the Arctic environment. Using a lagrangian ice age product we examine MYI loss and quantify the annual MYI area budget from 1980-2021 as the balance of export, melt and replenishment. Overall, MYI area declined at 72,500 km^2/yr, however a majority of the loss occurred during two stepwise reductions that interrupt an otherwise balanced budget and resulted in northward contractions of the MYI pack. First, in 1989, a change in atmospheric forcing led to a +56% anomaly in MYI export through Fram Strait. The second occurred from 2006-2008 with anomalously high melt (+25%) and export (+23%) coupled with low replenishment (-8%). In terms of trends, melt has increased since 1989, particularly in the Beaufort Sea, export has decreased since 2008 due to reduced MYI coverage north of Fram Strait, and replenishment has increased over the full time series due to a negative feedback that promotes seasonal ice survival at higher latitudes exposed by MYI loss. However, retention to older MYI has significantly declined, transitioning the MYI pack towards younger MYI that is less resilient than previously anticipated and could soon elicit another stepwise reduction. We speculate that future MYI loss will be driven by increased melt and reduced replenishment, both of which are enhanced with continued warming and will one day render the Arctic Ocean free of MYI, a change that will coincide with a seasonally ice-free Arctic Ocean.
In the summer of 2020, ESA changed the orbit of CryoSat-2 to align periodically with NASA’s ICESat-2 mission, a campaign known as CRYO2ICE, which allows for near-coincident CryoSat-2 and ICESat-2 observations in space and time over the Arctic until summer 2022, where the CRYO2ICE Antarctic campaign was initiated. This study investigates the Arctic CRYO2ICE radar and laser freeboards acquired by CryoSat-2 and ICESat-2, respectively, during the winter seasons of 2020–2022, and derives snow depths from their differences along the orbits. Along-track snow depth observations can provide high-resolution snow depth distributions which are vital for air-ice-ocean heat and momentum transfer, understanding light transmission, and snow-ice-interactions. Generally, ICESat-2 is backscattered at a surface above the elevation of the CryoSat-2 signal. CRYO2ICE snow depths are thinner than the daily model- or passive-microwave-based snow depth composites used for comparison, with differences being most pronounced in the Atlantic and Pacific Arctic. Satellite-derived and model-based snow estimates show similar seasonal accumulation over FYI, but CRYO2ICE has limited seasonal accumulation over MYI which is linked to a slow increase in ICESat-2, and to some extent CryoSat-2, freeboards. We present a first estimation of along-track snow depth estimates with average uncertainty of 9 +/- 3 cm for 7-km segments, with random and systematic contributions of 7 and 4 cm. These observations show the potential for along-track dual-frequency observations of snow depth from the future Copernicus mission CRISTAL; but they also highlight uncertainties in radar penetration and the correlation length scales of snow topography that still require further research. 
A growing number of studies are concluding that the resilience of the Arctic sea ice cover in a warming climate is essentially controlled by its thickness. Satellite radar and laser altimeters have allowed us to routinely monitor sea ice thickness across most of the Arctic Ocean for several decades. However, a key uncertainty remaining in the sea ice thickness retrieval is the error on the sea surface height (SSH) which is conventionally interpolated at ice floes from a limited number of lead observations along the altimeter’s orbital track. Here, we use an objective mapping approach to determine sea surface height from all proximal lead samples located on the orbital track and from adjacent tracks within a neighborhood of 10s of kilometers. The patterns of the SSH signal’s zonal, meridional, and temporal decorrelation length scales are obtained by analyzing the covariance of historic CryoSat-2 Arctic lead observations, which match the scales obtained from an equivalent analysis of high-resolution sea ice-ocean model fields. We use these length scales to determine an optimal SSH and error estimate for each sea ice floe location. By exploiting leads from adjacent tracks, we can increase the SSH precision estimated at orbital crossovers by a factor of three. In regions of high SSH uncertainty, biases in CryoSat-2 sea ice freeboard can be reduced by 25% with respect to coincident airborne validation data. The new method is not restricted to a particular sensor or mode, so it can be generalized to all present and historic polar altimetry missions.