Sea Ice Thickness Distribution Research Articles

Combining a New Parameterization Scheme of Oceanic Heat Flux with Thickness Assimilation to Improve Sea Ice Forecast Accuracy in Liaodong Bay

Liaodong Bay is one of the lowest latitude areas with seasonal sea ice cover in the Northern Hemisphere. Sea ice forecasting faces challenges in accuracy due to its low thickness. Therefore, a novel parameterization scheme for oceanic heat flux was developed to forecast sea ice thickness accurately. Application of the parameterization scheme for oceanic heat flux in the HIGHTSI model significantly improved the forecasting accuracy of sea ice thickness before the severe ice period. During this period, the RMSE of sea ice thickness measured on the JZ9–3 and the JZ20–2 oil platforms decreased by 0.53 cm and 2.90 cm compared to previous schemes, respectively. In order to improve the forecasting accuracy during the severe and melting ice periods, the observed and retrieved sea ice thickness was nudged into the model. The simulated results demonstrated a good agreement with monitored sea ice thickness distribution. During the entire season with sea ice cover, the R-squared values between simulated and retrieved sea ice thickness in the core area of Liaodong Bay reached 0.76. Furthermore, this study revealed a relatively strong correlation between air temperature and ice temperature on the following day. The proposed scheme provides a valuable approach to improve the forecasting accuracy for the areas with low thickness in the sea ice numerical models.

Assessment of Arctic sea ice and surface climate conditions in nine CMIP6 climate models

ABSTRACT The observed retreat and anticipated further decline in Arctic sea ice holds strong climate, environmental, and societal implications. In predicting climate evolution, ensembles of coupled climate models have demonstrated appreciable accuracy in simulating sea-ice area trends throughout the historical period, yet individual climate models still show significant differences in accurately representing the sea-ice thickness distribution. To better understand individual model performance in sea-ice simulation, nine climate models were evaluated in comparison with Arctic satellite and reanalysis-derived sea-ice thickness data, sea-ice area records, and atmospheric reanalysis data of surface wind and air temperature. This assessment found that the simulated spatial distribution of historical sea-ice thickness varies greatly between models and that several key limitations persist among models. Primarily, most models do not capture the thickest regimes of multiyear ice present in the Wandel and Lincoln seas; those that do often possess erroneous positive bias in other regions such as the Laptev Sea or along the Eurasian Arctic Shelf. This analysis provides enhanced understanding of individual model historical simulation performance, which is critical in informing the selection of coupled climate model projections for dependent future modeling efforts.

Improvement of sea ice thermodynamics with variable sea ice salinity and melt pond parameterizations in an OGCM

Enhanced representation of sea ice processes in ocean modeling studies is required for advancing our understanding and prediction of the climate variability. In the present study, we improved the sea ice thermodynamics in an OGCM by introducing processes for variable sea ice salinity and melt ponds on sea ice. The former affects the latent heat and conductivity of sea ice as well as the salt budget in the ocean–sea ice system. We formulated salinification and desalination processes based on previous observational studies. Melt ponds enhance summertime surface melting by reducing the surface albedo. In our formulation of the melt pond thermodynamics, the shortwave radiation absorbed in the pond was partly released to the air and partly conducted to the underlying sea ice. By adopting these parameterizations, seasonal cycles of the sea ice salinity distributions in both hemispheres were generally reproduced. Simulated seasonal and interannual variations in the Arctic melt pond area and associated surface albedo reductions were basically consistent with observational estimates. The melt pond effect reduced the Arctic sea ice volume by about 1.5 × 103 km3 in the summer; half of this impact remained throughout the year. Our results suggest that the new parameterizations have the potential to improve the seasonal evolution of the sea ice thickness and concentration distributions, which would contribute to prediction studies.

The effects of assimilating a sub-grid-scale sea ice thickness distribution in a new Arctic sea ice data assimilation system

Abstract. In the past decade groundbreaking new satellite observations of the Arctic sea ice cover have been made, allowing researchers to understand the state of the Arctic sea ice system in greater detail than before. The derived estimates of sea ice thickness are useful but limited in time and space. In this study the first results of a new sea ice data assimilation system are presented. Observations assimilated (in various combinations) are monthly mean sea ice thickness and monthly mean sea ice thickness distribution from CryoSat-2 and NASA daily Bootstrap sea ice concentration. This system couples the Centre for Polar Observation and Modelling's (CPOM) version of the Los Alamos Sea Ice Model (CICE) to the localised ensemble transform Kalman filter (LETKF) from the Parallel Data Assimilation Framework (PDAF) library. The impact of assimilating a sub-grid-scale sea ice thickness distribution is of particular novelty. The sub-grid-scale sea ice thickness distribution is a fundamental component of sea ice models, playing a vital role in the dynamical and thermodynamical processes, yet very little is known of its true state in the Arctic. This study finds that assimilating CryoSat-2 products for the mean thickness and the sub-grid-scale thickness distribution can have significant consequences for the modelled distribution of the ice thickness across the Arctic and particularly in regions of thick multi-year ice. The assimilation of sea ice concentration, mean sea ice thickness and sub-grid-scale sea ice thickness distribution together performed best when compared to a subset of CryoSat-2 observations held back for validation. Regional model biases are reduced: the thickness of the thickest ice in the Canadian Arctic Archipelago (CAA) is decreased, but the thickness of the ice in the central Arctic is increased. When comparing the assimilation of mean thickness with the assimilation of sub-grid-scale thickness distribution, it is found that the latter leads to a significant change in the volume of ice in each category. Estimates of the thickest ice improve significantly with the assimilation of sub-grid-scale thickness distribution alongside mean thickness.

Seasonal Evolution of the Arctic Sea Ice Thickness Distribution

AbstractThe Thorndike et al. (1975, https://doi.org/10.1029/jc080i033p04501) theory of the ice thickness distribution, g(h), treats the dynamic and thermodynamic aggregate properties of the ice pack in a novel and physically self‐consistent manner. Therefore, it has provided the conceptual basis of the treatment of sea‐ice thickness categories in climate models. The approach, however, is not mathematically closed due to the treatment of mechanical deformation using the redistribution function ψ, the authors noting “The present theory suffers from a burdensome and arbitrary redistribution function ψ.” Toppaladoddi and Wettlaufer (2015, https://doi.org/10.1103/physrevlett.115.148501) showed how ψ can be written in terms of g(h), thereby solving the mathematical closure problem and writing the theory in terms of a Fokker‐Planck equation, which they solved analytically to quantitatively reproduce the observed winter g(h). Here, we extend this approach to include open water by formulating a new boundary condition for their Fokker‐Planck equation, which is then coupled to the observationally consistent sea‐ice growth model of Semtner (1976, https://doi.org/10.1175/1520-0485(1976)006<0379:amfttg>2.0.co;2) to study the seasonal evolution of g(h). We find that as the ice thins, g(h) transitions from a single‐ to a double‐peaked distribution, which is in agreement with observations. To understand the cause of this transition, we construct a simpler description of the system using the equivalent Langevin equation formulation and solve the resulting stochastic ordinary differential equation numerically. Finally, we solve the Fokker‐Planck equation for g(h) under different climatological conditions to study the evolution of the open‐water fraction.

Simple Hybrid Sea Ice Nudging Method for Improving Control Over Partitioning of Sea Ice Concentration and Thickness

AbstractTo assess the effect of ocean‐atmosphere coupling in the climate response to forced sea ice loss, the Polar Amplification Model Intercomparison Project protocol includes centennial coupled atmosphere‐ocean general circulation model simulations with imposed sea ice loss. The protocol, which specifies sea ice concentration and thickness distribution targets, does not prescribe a method for achieving them. Although different methods for imposing sea ice loss (or growth) in models have been documented, testing of the method‐dependence of the resulting climate responses has been limited. Achieving the targeted sea ice state has proven to be challenging using the so‐called ghost‐flux nudging method, which induces ice melt from below, as this method does not constrain the partitioning between thickness and concentration. We propose, describe and test a simple method that combines the advantages of direct sea ice nudging and ghost‐flux nudging. The hybrid nudging method better captures the partitioning between thickness and concentration while conserving total water content. We document how this novel sea ice constraining method reaches specific targets, enhances surface turbulent heat flux responses to sea ice loss, and induces tropospheric warming for both polar regions.

Effects of Increasing the Category Resolution of the Sea Ice Thickness Distribution in a Coupled Climate Model on Arctic and Antarctic Sea Ice Mean State

AbstractMany modern sea ice models used in global climate models represent the subgrid‐scale heterogeneity in sea ice thickness with an ice thickness distribution (ITD), which improves model realism by representing the significant impact of the high spatial heterogeneity of sea ice thickness on thermodynamic and dynamic processes. Most models default to five thickness categories. However, little has been done to explore the effects of the resolution of this distribution (number of categories) on sea‐ice feedbacks in a coupled model framework and resulting representation of the sea ice mean state. Here, we explore this using sensitivity experiments in CESM2 with the standard 5 ice thickness categories and 15 ice thickness categories. Increasing the resolution of the ITD in a run with preindustrial climate forcing results in substantially thicker Arctic sea ice year‐round. Analyses show that this is a result of the ITD influence on ice strength. With 15 ITD categories, weaker ice occurs for the same average thickness, resulting in a higher fraction of ridged sea ice. In contrast, the higher resolution of thin ice categories results in enhanced heat conduction and bottom growth and leads to only somewhat increased winter Antarctic sea ice volume. The spatial resolution of the ICESat‐2 satellite mission provides a new opportunity to compare model outputs with observations of seasonal evolution of the ITD in the Arctic (ICESat‐2; 2018–2021). Comparisons highlight significant differences from the ITD modeled with both runs over this period, likely pointing to underlying issues contributing to the representation of average thickness.

Thermodynamic and dynamic contributions to seasonal Arctic sea ice thickness distributions from airborne observations

Sea ice thickness is a key parameter in the polar climate and ecosystem. Thermodynamic and dynamic processes alter the sea ice thickness. The Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition provided a unique opportunity to study seasonal sea ice thickness changes of the same sea ice. We analyzed 11 large-scale (∼50 km) airborne electromagnetic sea thickness and surface roughness surveys from October 2019 to September 2020. Data from ice mass balance and position buoys provided additional information. We found that thermodynamic growth and decay dominated the seasonal cycle with a total mean sea ice thickness increase of 1.4 m (October 2019 to June 2020) and decay of 1.2 m (June 2020 to September 2020). Ice dynamics and deformation-related processes, such as thin ice formation in leads and subsequent ridging, broadened the ice thickness distribution and contributed 30% to the increase in mean thickness. These processes caused a 1-month delay between maximum thermodynamic sea ice thickness and maximum mean ice thickness. The airborne EM measurements bridged the scales from local floe-scale measurements to Arctic-wide satellite observations and model grid cells. The spatial differences in mean sea ice thickness between the Central Observatory (<10 km) of MOSAiC and the Distributed Network (<50 km) were negligible in fall and only 0.2 m in late winter, but the relative abundance of thin and thick ice varied. One unexpected outcome was the large dynamic thickening in a regime where divergence prevailed on average in the western Nansen Basin in spring. We suggest that the large dynamic thickening was due to the mobile, unconsolidated sea ice pack and periodic, sub-daily motion. We demonstrate that this Lagrangian sea ice thickness data set is well suited for validating the existing redistribution theory in sea ice models. Our comprehensive description of seasonal changes of the sea ice thickness distribution is valuable for interpreting MOSAiC time series across disciplines and can be used as a reference to advance sea ice thickness modeling.

An evaluation of the E3SMv1 Arctic ocean and sea-ice regionally refined model

Abstract. The Energy Exascale Earth System Model (E3SM) is a state-of-the-science Earth system model (ESM) with the ability to focus horizontal resolution of its multiple components in specific areas. Regionally refined global ESMs are motivated by the need to explicitly resolve, rather than parameterize, relevant physics within the regions of refined resolution, while offering significant computational cost savings relative to the respective cost of configurations with high-resolution (HR) everywhere on the globe. In this paper, we document results from the first Arctic regionally refined E3SM configuration for the ocean and sea-ice components (E3SM-Arctic-OSI), while employing data-based atmosphere, land, and hydrology components. Our aim is an improved representation of the Arctic coupled ocean and sea-ice state, its variability and trends, and the exchanges of mass and property fluxes between the Arctic and the sub-Arctic. We find that E3SM-Arctic-OSI increases the realism of simulated Arctic ocean and sea-ice conditions compared to a similar low-resolution E3SM simulation without the Arctic regional refinement in ocean and sea-ice components (E3SM-LR-OSI). In particular, exchanges through the main Arctic gateways are greatly improved with respect to E3SM-LR-OSI. Other aspects, such as the Arctic freshwater content variability and sea-ice trends, are also satisfactorily simulated. Yet, other features, such as the upper-ocean stratification and the sea-ice thickness distribution, need further improvements, involving either more advanced parameterizations, model tuning, or additional grid refinements. Overall, E3SM-Arctic-OSI offers an improved representation of the Arctic system relative to E3SM-LR-OSI, at a fraction (15 %) of the computational cost of comparable global high-resolution configurations, while permitting exchanges with the lower-latitude oceans that cannot be directly accounted for in Arctic regional models.

Arctic sea ice anomalies during the MOSAiC winter 2019/20

Abstract. During the winter of 2019/2020, as the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) project started its work, the Arctic Oscillation (AO) experienced some of its largest shifts, ranging from a highly negative index in November 2019 to an extremely positive index during January–February–March (JFM) 2020. The permanent positive AO phase for the 3 months of JFM 2020 was accompanied by a prevailing positive phase of the Arctic Dipole (AD) pattern. Here we analyze the sea ice thickness (SIT) distribution based on CryoSat-2/SMOS satellite-derived data augmented with results from the hindcast simulation by the fully coupled Regional Arctic System Model (RASM) from November 2019 through March 2020. A notable result of the positive AO phase during JFM 2020 was large SIT anomalies of up to 1.3 m that emerged in the Barents Sea (BS), along the northeastern Canadian coast and in parts of the central Arctic Ocean. These anomalies appear to be driven by nonlinear interactions between thermodynamic and dynamic processes. In particular, in the Barents and Kara seas (BKS), they are a result of enhanced ice growth connected with low-temperature anomalies and the consequence of intensified atmospherically driven sea ice transport and deformations (i.e., ice divergence and shear) in this area. The Davies Strait, the east coast of Greenland and the BS regions are characterized by convergence and divergence changes connected with thinner sea ice at the ice borders along with an enhanced impact of atmospheric wind forcing. Low-pressure anomalies that developed over the eastern Arctic during JFM 2020 increased northerly winds from the cold Arctic Ocean to the BS and accelerated the southward drift of the MOSAiC ice floe. The satellite-derived and simulated sea ice velocity anomalies, which compared well during JFM 2020, indicate a strong acceleration of the Transpolar Drift relative to the mean for the past decade, with intensified speeds of up to 6 km d−1. As a consequence, sea ice transport and deformations driven by atmospheric surface wind forcing accounted for the bulk of the SIT anomalies, especially in January 2020 and February 2020. RASM intra-annual ensemble forecast simulations with 30 ensemble members forced with different atmospheric boundary conditions from 1 November 2019 through 30 April 2020 show a pronounced internal variability in the sea ice volume, driven by thermodynamic ice-growth and ice-melt processes and the impact of dynamic surface winds on sea ice formation and deformation. A comparison of the respective SIT distributions and turbulent heat fluxes during the positive AO phase in JFM 2020 and the negative AO phase in JFM 2010 corroborates the conclusion that winter sea ice conditions in the Arctic Ocean can be significantly altered by AO variability.

A new state-dependent parameterization for the free drift of sea ice

Abstract. Free-drift estimates of sea ice motion are necessary to produce a seamless observational record combining buoy and satellite-derived sea ice motion vectors. We develop a new parameterization for the free drift of sea ice based on wind forcing, wind turning angle, sea ice state variables (thickness and concentration), and estimates of the ocean currents. Given the fact that the spatial distribution of the wind–ice–ocean transfer coefficient has a similar structure to that of the spatial distribution of sea ice thickness, we take the standard free-drift equation and introduce a wind–ice–ocean transfer coefficient that scales linearly with ice thickness. Results show a mean bias error of −0.5 cm s−1 (low-speed bias) and a root-mean-square error of 5.1 cm s−1, considering daily buoy drift data as truth. This represents a 35 % reduction of the error on drift speed compared to the free-drift estimates used in the Polar Pathfinder dataset (Tschudi et al., 2019b). The thickness-dependent transfer coefficient provides an improved seasonality and long-term trend of the sea ice drift speed, with a minimum (maximum) drift speed in May (October), compared to July (January) for the constant transfer coefficient parameterizations which simply follow the peak in mean surface wind stresses. Over the 1979–2019 period, the trend in sea ice drift in this new model is +0.45 cm s−1 per decade compared with +0.39 cm s−1 per decade from the buoy observations, whereas there is essentially no trend in a free-drift parameterization with a constant transfer coefficient (−0.09 cm s−1 per decade) or the Polar Pathfinder free-drift input data (−0.01 cm s−1 per decade). The optimal wind turning angle obtained from a least-squares fitting is 25∘, resulting in a mean error and a root-mean-square error of +3 and 42∘ on the direction of the drift, respectively. The ocean current estimates obtained from the minimization procedure resolve key large-scale features such as the Beaufort Gyre and Transpolar Drift Stream and are in good agreement with ocean state estimates from the ECCO, GLORYS, and PIOMAS ice–ocean reanalyses, as well as geostrophic currents from dynamical ocean topography, with a root-mean-square difference of 2.4, 2.9, 2.6, and 3.8 cm s−1, respectively. Finally, a repeat of the analysis on two sub-sections of the time series (pre- and post-2000) clearly shows the acceleration of the Beaufort Gyre (particularly along the Alaskan coastline) and an expansion of the gyre in the post-2000s, concurrent with a thinning of the sea ice cover and the observed acceleration of the ice drift speed and ocean currents. This new dataset is publicly available for complementing merged observation-based sea ice drift datasets that include satellite and buoy drift records.

Spatio-temporal distribution of sea-ice thickness using a machine learning approach with Google Earth Engine and Sentinel-1 GRD data

Knowledge of ice thickness and its distribution is of great interest to plan ship and offshore operations in ice areas. It is also a major direct factor in climate change. However, it is currently the most inaccurate sea-ice parameter according to a large majority of the scientific community. Satellite-borne synthetic aperture radar (SAR) data are highly valuable for monitoring daily ice-covered oceans. In this study, we develop a new technique based on a machine learning Random Forest (RF) regression approach using the combination of the in-situ thickness measurements and the backscattering information from Sentinel-1 to retrieve the level first-year ice (FYI) thickness. By applying the technique over the Sentinel-1 ground range detected (GRD) data set available in Google Earth Engine (GEE) over the Beaufort Sea spanning a time period from Apr 2015 to Sep 2018, a thickness of up to 1.5 m with a root mean square error (RMSE) of 22 cm is retrieved.

Towards a swath-to-swath sea-ice drift product for the Copernicus Imaging Microwave Radiometer mission

Abstract. Across spatial and temporal scales, sea-ice motion has implications for ship navigation, the sea-ice thickness distribution, sea-ice export to lower latitudes and re-circulation in the polar seas, among others. Satellite remote sensing is an effective way to monitor sea-ice drift globally and daily, especially using the wide swaths of passive microwave missions. Since the late 1990s, many algorithms and products have been developed for this task. Here, we investigate how processing sea-ice drift vectors from the intersection of individual swaths of the Advanced Microwave Scanning Radiometer 2 (AMSR2) mission compares to today's status quo (processing from daily averaged maps of brightness temperature). We document that the “swath-to-swath” (S2S) approach results in many more (2 orders of magnitude) sea-ice drift vectors than the “daily map” (DM) approach. These S2S vectors also validate better when compared to trajectories of on-ice drifters. For example, the RMSE of the 24 h winter Arctic sea-ice drift is 0.9 km for S2S vectors and 1.3 km for DM vectors from the 36.5 GHz imagery of AMSR2. Through a series of experiments with actual AMSR2 data and simulated Copernicus Imaging Microwave Radiometer (CIMR) data, we study the impact that geolocation uncertainty and imaging resolution have on the accuracy of the sea-ice drift vectors. We conclude by recommending that a swath-to-swath approach is adopted for the future operational Level-2 sea-ice drift product of the CIMR mission. We outline some potential next steps towards further improving the algorithms and making the user community ready to fully take advantage of such a product.

Tenuous Correlation between Snow Depth or Sea Ice Thickness and C- or X-Band Backscattering in Nunavik Fjords of the Hudson Strait

Radar penetration in brine-wetted snow-covered sea ice is almost nil, yet reports exist of a correlation between snow depth or ice thickness and SAR parameters. This article presents a description of snow depth and first-year sea ice thickness distributions in three fjords of the Hudson Strait and of their tenuous correlation with SAR backscattering in the C- and X-band. Snow depth and ice thickness were directly measured in three fjords of the Hudson Strait from 2015 to 2018 in April or May. Bayesian linear regression analysis was used to investigate their relationship with RADARSAT-2 (C-band) or TerraSAR-X (X-band). Polarimetric ratios and the Cloude–Pottier decomposition parameters were explored along with the HH, HV and VV bands. Linear correlations were generally no higher than 0.3 except for a special case in May 2018. The co-polarization ratio did not perform better than the backscattering coefficients.

Airborne mapping of the sub-ice platelet layer under fast ice in McMurdo Sound, Antarctica

Abstract. Basal melting of ice shelves can result in the outflow of supercooled ice shelf water, which can lead to the formation of a sub-ice platelet layer (SIPL) below adjacent sea ice. McMurdo Sound, located in the southern Ross Sea, Antarctica, is well known for the occurrence of a SIPL linked to ice shelf water outflow from under the McMurdo Ice Shelf. Airborne, single-frequency, frequency-domain electromagnetic induction (AEM) surveys were performed in November of 2009, 2011, 2013, 2016, and 2017 to map the thickness and spatial distribution of the landfast sea ice and underlying porous SIPL. We developed a simple method to retrieve the thickness of the consolidated ice and SIPL from the EM in-phase and quadrature components, supported by EM forward modelling and calibrated and validated by drill-hole measurements. Linear regression of EM in-phase measurements of apparent SIPL thickness and drill-hole measurements of “true” SIPL thickness yields a scaling factor of 0.3 to 0.4 and rms error of 0.47 m. EM forward modelling suggests that this corresponds to SIPL conductivities between 900 and 1800 mS m−1, with associated SIPL solid fractions between 0.09 and 0.47. The AEM surveys showed the spatial distribution and thickness of the SIPL well, with SIPL thicknesses of up to 8 m near the ice shelf front. They indicate interannual SIPL thickness variability of up to 2 m. In addition, they reveal high-resolution spatial information about the small-scale SIPL thickness variability and indicate the presence of persistent peaks in SIPL thickness that may be linked to the geometry of the outflow from under the ice shelf.

Sea Ice Thickness in the Western Ross Sea

AbstractUsing airborne measurements, we provide a first direct glimpse of the sea ice thickness distribution in the western Ross Sea, Antarctica, where the distinguishing sea ice process is the regular occurrence of the Ross Sea, McMurdo Sound, and Terra Nova Bay polynyas. Two flights in November 2017 over a length of 800 km reveal a heavily deformed ice regime with a mean thickness of 2.0 ± 1.6 m. Supported by satellite image analysis, we identify regional variability in ice thickness based on formation history. Sea ice thickness gradients are highest within 100 and 200 km of the Terra Nova Bay and McMurdo Sound polynyas, respectively, where the mean thickness of the thickest 10% of ice is 7.6 m. Overall, about 80% of the ice is heavily deformed, concentrated in ridges with thicknesses of 3.0–11.8 m. This is evidence that sea ice is much thicker than in the central Ross Sea.

Impact of the ice thickness distribution discretization on the sea ice concentration variability in the NEMO3.6–LIM3 global ocean–sea ice model

Abstract. This study assesses the impact of different sea ice thickness distribution (ITD) discretizations on the sea ice concentration (SIC) variability in ocean stand-alone NEMO3.6–LIM3 simulations. Three ITD discretizations with different numbers of sea ice thickness categories and boundaries are evaluated against three different satellite products (hereafter referred to as “data”). Typical model and data interannual SIC variability is characterized by K-means clustering both in the Arctic and Antarctica between 1979 and 2014. We focus on two seasons, winter (January–March) and summer (August–October), in which correlation coefficients across clusters in individual months are largest. In the Arctic, clusters are computed before and after detrending the series with a second-degree polynomial to separate interannual from longer-term variability. The analysis shows that, before detrending, winter clusters reflect the SIC response to large-scale atmospheric variability at both poles, while summer clusters capture the negative and positive trends in Arctic and Antarctic SIC, respectively. After detrending, Arctic clusters reflect the SIC response to interannual atmospheric variability predominantly. The cluster analysis is complemented with a model–data comparison of the sea ice extent and SIC anomaly patterns. The single-category discretization shows the worst model–data agreement in the Arctic summer before detrending, related to a misrepresentation of the long-term melting trend. Similarly, increasing the number of thin categories reduces model–data agreement in the Arctic, due to a poor representation of the summer melting trend and an overly large winter sea ice volume associated with a net increase in basal ice growth. In contrast, more thin categories improve model realism in Antarctica, and more thick ones improve it in central Arctic regions with very thick ice. In all the analyses we nonetheless identify no optimal discretization. Our results thus suggest that no clear benefit in the representation of SIC variability is obtained from increasing the number of sea ice thickness categories beyond the current standard with five categories in NEMO3.6–LIM3.

Impact of Arctic sea ice floe-scale anisotropy on airborne electromagnetic surveys

AbstractAirborne electromagnetic induction sensors have demonstrated their extensive capacities to measure sea-ice thickness distributions. However, biases can emerge when comparing these 1-D measurements to a broader 2-D regional scale due to the spatial anisotropy inherent to sea-ice cover. Automated processing of available sea-ice maps could significantly ease the decision on how to set up an optimised flight pattern, which would result in representative ice thickness numbers for the region. In this study, first we investigate the extent to which the sea-ice anisotropy can influence the representativeness of an airborne survey compared to the regional situation. Second, we propose a method to process sea-ice maps prior to flights to help preparing the most representative flight plan possible for the local area. The method is based on automated segmentation of radar satellite images and extensive simulation of flight transects over the image. The spatial analysis of these transects enables for the identification of the most representative survey trajectories for the area. The method was applied for seven different synthetic aperture radar satellite images over Arctic sea ice north of Svalbard. The results indicate that the proposed method improved the representativeness of the airborne survey by identifying the most suitable transect over the ice pack.

A Textural Approach to Improving Snow Depth Estimates in the Weddell Sea

The snow depth on Antarctic sea ice is critical to estimating the sea ice thickness distribution from laser altimetry data, such as from Operation IceBridge or ICESat-2. Snow redistributed by wind collects around areas of deformed ice and forms a wide variety of features on sea ice; the morphology of these features may provide some indication of the mean snow depth. Here, we apply a textural segmentation algorithm to classify and group similar textures to infer the distribution of snow using snow surface freeboard measurements from Operation IceBridge campaigns over the Weddell Sea. We find that texturally-similar regions have similar snow/ice ratios, even when they have different absolute snow depth measurements. This allows for the extrapolation of nadir-looking snow radar data using two-dimensional surface altimetry scans, providing a two-dimensional estimate of the snow depth with ∼22% error. We show that at the floe scale (∼180 m), snow depth can be directly estimated from the snow surface with ∼20% error using deep learning techniques, and that the learned filters are comparable to standard textural analysis techniques. This error drops to ∼14% when averaged over 1.5 km scales. These results suggest that surface morphological information can improve remotely-sensed estimates of snow depth, and hence sea ice thickness, as compared to current methods. Such methods may be useful for reducing uncertainty in Antarctic sea ice thickness estimates from ICESat-2.

Winter Arctic Sea Ice Thickness From ICESat‐2 Freeboards

AbstractNational Aeronautics and Space Administration's (NASA's) Ice, Cloud, and land Elevation Satellite‐2 (ICESat‐2) mission was launched in September 2018 with the primary goal of monitoring our rapidly changing polar regions. The sole instrument onboard, the Advanced Topographic Laser Altimeter System, is now providing routine, very high‐resolution, surface elevation data across the globe, including the Arctic and Southern oceans. In this study, we demonstrate our new processing chain for converting the along‐track ICESat‐2 sea ice freeboard product (ATL10) into sea ice thickness, focusing our initial efforts on the Arctic Ocean. For this conversion, we primarily make use of snow depth and density data from the NASA Eulerian Snow on Sea Ice Model. The coarse resolution (~100 km) snow data are redistributed onto the high‐resolution (approximately 30–100 m) ATL10 freeboards using relationships obtained from snow depth and freeboard data collected by NASA's Operation IceBridge mission. We present regional sea ice thickness distributions and highlight their seasonal evolution through our first winter season of data collection. We include ice thickness uncertainty estimates, while also acknowledging the limitations of these estimates. We generate a gridded monthly thickness product and compare this with various monthly sea ice thickness estimates obtained from European Space Agency's CryoSat‐2 satellite mission, with ICESat‐2 showing consistently lower thicknesses. Finally, we compare our February/March 2019 thickness estimates to ICESat February/March (19 February to 21 March) 2008 ice thickness estimates using the same input assumptions, which show an ~0.37 m or ~20% thinning across an inner Arctic Ocean domain in this 11‐year time period.