Description: 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.

Description: Knowledge of sea-ice thickness and volume depends on freeboard observations from satellite altimeters and in turn on information of snow mass and sea-ice density required for the freeboard-to-thickness conversion. These parameters, especially sea-ice density, are usually based on climatologies constructed from in situ observations made in the 1980s and earlier while contemporary and representative measurements are lacking. Our aim with this paper is to derive updated sea-ice bulk density estimates suitable for the present Arctic sea-ice cover and a range of ice types to reduce uncertainties in sea-ice thickness remote sensing. Our sea-ice density measurements are based on over 3000 km of high-resolution collocated airborne sea-ice and snow thickness and freeboard measurements in the western Arctic Ocean in 2017 and 2019. Sea-ice bulk density is derived assuming isostatic equilibrium for different ice types. Our results show higher average bulk densities for both first-year ice (FYI) and especially multi-year ice (MYI) compared to previous studies. In addition, we find a small difference between deformed and possibly unconsolidated FYI and younger MYI. We find a negative-exponential relationship between sea-ice bulk density and sea-ice freeboard and apply this parameterisation to one winter of monthly gridded CryoSat-2 sea-ice freeboard data. We discuss the suitability and the impact of the derived FYI and MYI bulk densities for sea-ice thickness retrievals and the uncertainty related to the indirect method of measuring sea-ice bulk density. The results suggest that retrieval algorithms be adapted to changes in sea-ice density and highlight the need of future studies to evaluate the impact of density parameterisation on the full sea-ice thickness data record.

Description: An integrated sensor platform including an inertial navigation system (INS) and a commercial airborne laser scanner (ALS) among other sensor was mounted in the cargo compartment in one of the Polarstern helicopters during MOSAiC. ALS data was acquired from more than 60 flights between October 2019 and September 2020 with a range of survey types intended to map changes of the sea ice surface during the full annual cycle at high spatial resolution and coverage. Here, we provide an overview of the collected data, the challenge of achieving centimeter elevation accuracy with a helicopter platform at high polar latitudes as well as the content and specifications of ALS data products. The high spatial resolution and repeated coverage of the larger area around Polarstern allow studying various surface features (e.g. pressure ridges, floes, melt ponds, snow drifts, etc.), their seasonal evolution, and their impact on atmosphere and ocean. Finally, we outline methods for planned applications, such as identifying individual floes and surface types using both measured freeboard and surface reflectance. Collocated helicopter-based optical and infrared imagery allow analyzing sea ice properties in further applications and to upscale comparable in-situ observations.

Description: Arctic sea ice is a mosaic of ice floes whose distribution greatly impacts the interaction of sea ice with the atmosphere and the ocean. However, we are still lacking physics to describe the complex interplay of ice floes that are a key characteristic of sea ice. Kilometer-scale satellite and buoy observations have been used to describe the deformation of sea-ice as a continuum of floes, whereas ice fracture of individual floes has been studied in laboratory and field experiments. In our contribution, we aim to bridge this gap and outline a framework to characterize sea-ice deformation at the floe-scale from observational data by studying the mechanical interaction of multiple identifiable floes. We use Sentinel SAR imagery and airborne laser scanner surveys acquired during the MOSAiC expedition to map ice floes in the larger area around Polarstern. This combination of data products allows us to describe the floe-size distribution from hundreds of kilometers down to the meter-scale. With the repeated coverage of both data products, ice motion is tracked and deformation estimates are derived. By combining both floe-size estimates and deformation rates we provide insights into how the floe composition changes in regions that were exposed to deformation. This relationship between floe sizes and mechanical redistribution could constrain parameterizations to improve the representation of sea-ice deformation and floe sizes in large-scale continuum models. Finally, we highlight the potential to use the presented floe and deformation information to initialise and evaluate floe-resolving sea-ice simulations using discrete element models.

Description: Sea-ice thickness is a key factor and indicator in understanding the impact of the global climate change. Deriving basin-wide sea-ice thickness estimates from satellite laser and radar altimetry relies on freeboard measurements. The freeboard-to-thickness conversion in turn requires information of snow mass and the density of the sea-ice layer that have unknown spatio-temporal variabilities and trends directly translating into the uncertainty of decadal sea-ice thickness data records. In addition, inter-mission biases arise from, e.g., different sensor types and frequencies as well as varying footprint sizes affected by surface roughness across regions and seasons. Therefore, carrying out validation and inter-calibration studies is crucial for reliable and continuous observation of the Earth’s cryosphere. To achieve this, it is beneficial to have simultaneous measurements of freeboard, snow depth, and sea-ice thickness, which provide reference data for both direct satellite observations and geophysical target parameters. Here, we present Alfred Wegener Institute’s (AWI) IceBird program, which is a series of fixed-wing aircraft campaigns to measure Arctic sea ice and to monitor its change. During two late-winter campaigns in the western Arctic Ocean in 2017 and 2019, we have carried out surveys with the unique scientific instrument configuration including an airborne laser scanner (ALS) for surface topography and freeboard measurements, a tethered electromagnetic induction sounding instrument (EM-Bird) for total (snow+ice) thickness measurements, and an ultrawideband frequency-modulated continuous-wave microwave radar to measure snow thickness. Therefore, we are able to observe all three bounding interfaces in the sea-ice–snow system in high resolution along survey tracks on regional scales. During the ship-based drift expedition Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) between October 2019 and September 2020, helicopter surveys were carried out in high spatio-temporal resolution throughout the year, including the polar night, to measure freeboard and roughness with the ALS both in local grid pattern and in larger scale. Coincident EM-Bird ice thickness data and information from snow measurements on the ground will help linking these parameters and monitor them and their effect on satellite retrievals for a full seasonal cycle. The individual parameters are important for describing and monitoring the state of the Arctic sea ice and validating retrievals from satellite data, but combined they offer further possibilities to characterise sea ice. By assuming isostatic equilibrium, we are able to estimate up-to-date bulk density values for different sea-ice types from the IceBird data and to derive a parametrisation of sea-ice bulk density based on sea-ice freeboard. These data allow us to explore spatio-temporal variations in sea-ice parameters observable from space and to evaluate the validity of the freeboard-to-thickness conversion in satellite altimetry through comparison against dedicated satellite overpasses and orbit collections.

Description: Discrete Element Models presented in literature have been mostly run in idealized configurations focussing on certain aspects of sea ice fracture and deformation. Realistic simulations with DEMs initialised, forced, and validated with observational data are limited due to the lack of sea ice observations at floe scale. In our presentation, we outline a possible framework to use MOSAiC observations to design realistic DEM configurations to simulate pack-ice. We present airborne laser scanner (ALS) data of the sea ice surface elevation retrieved during the MOSAiC campaign, from which floe size, shape, and thickness information in an 10x5 km domain can be derived. With a spatial resolution of 50cm the ALS offers the unique opportunity to map meter-scale sea ice floes in 25 records around the year can be derived. These floe observations could be used to initialise the DEM simulations. The simulations could be forced with drift observations derived from satellite or buoy observations along the boundaries. Finally, the simulations could be evaluated against the floe observations derived from the next ALS record.

Description: Arctic sea ice is a mosaic of ice floes whose distribution greatly impacts the interaction of sea ice with the atmosphere and the ocean. However, we are still lacking physics to describe the complex interplay of ice floes that are a key characteristic of sea ice. Kilometer-scale satellite and buoy observations have been used to describe the deformation of sea-ice as a continuum of floes, whereas ice fracture of individual floes has been studied in laboratory and field experiments. In our contribution, we aim to bridge this gap and outline a framework to characterize sea-ice deformation at the floe-scale from observational data by studying the mechanical interaction of multiple identifiable floes. We use Sentinel SAR imagery and airborne laser scanner surveys acquired during the MOSAiC expedition to map ice floes in the larger area around Polarstern. This combination of data products allows us to describe the floe-size distribution from hundreds of kilometers down to the meter-scale. With the repeated coverage of both data products, ice motion is tracked and deformation estimates are derived. By combining both floe-size estimates and deformation rates we provide insights into how the floe composition changes in regions that were exposed to deformation. This relationship between floe sizes and mechanical redistribution could constrain parameterizations to improve the representation of sea-ice deformation and floe sizes in large-scale continuum models. Finally, we highlight the potential to use the presented floe and deformation information to initialise and evaluate floe-resolving sea-ice simulations using discrete element models.

Description: In the polar oceans in winter, fractures and leads are the hotspots of exchange between the ocean and atmosphere which are otherwise well separated by sea ice. By altering the heat, gas, and momentum fluxes they play a crucial role in atmospheric, ecological, and oceanic processes. At the same time, leads represent a part of the present state of strain of the ice cover, opening up the possibility to study ice rheology. The transient nature of leads and their narrow appearance has set limits to the detection of leads from satellites. Different approaches using active and passive sensors from the microwave and infrared spectrum are employed so far to observe leads by means of satellite data. They make use of the strong contrast between leads and the surrounding ice pack in (i) surface temperature, (ii) microwave backscatter, (iii) emission or (iv) a change in ice drift speed. With the increasing availability of high-resolution SAR data for the Arctic, we explored the potential to use SAR derived sea ice deformation to estimate lead fractions. We calculated sea ice drift and divergence with a spatial resolution of 1.4 km from daily Sentinel-1 scenes. We obtained the divergence-based lead fraction of a region by summing up all positive divergence pixels multiplied by the respective time step length. We derived a second lead fraction product from the deformation fields that calculates the position of linear kinematic features (LKFs) first. The advantage is a skilled noise reduction, and a tracking algorithm of the deformation zones. We compared divergence- and LKF-based lead fractions to several other established lead fraction products in the Transpolar Drift along the drift track of the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) between October 2019 to April 2020. We used lead fractions from helicopter-borne infrared surveys at a grid resolution of 5 m, classified Sentinel-1 (SAR) scenes at 80 m, MODIS (thermal infrared) at 1 km, AMSR2 (passive microwaves) at 3.25 km, and CryoSat-2 (altimeter in Ku-band) at 12.5 km. Since the methods rely on different physical properties of the water and ice in leads and are affected by different constraints, derived mean lead fractions vary by 1-2 magnitudes between the products. For example, infrared, SAR and microwave radiometer-based algorithms do not only detect open-water leads but also leads with thin ice up to a certain thickness, which differs between the products. Common lead events were identified across products. The time series mostly indicated a phase of increased lead activity during freeze-up in autumn 2019 and spring 2020. We used the different lead fraction time series to estimate new ice formation in the leads and compared the results to ice thickness and oceanographic measurements obtained during the MOSAiC campaign. Results yield lower and upper bounds for ice formation and brine expulsion in and from leads. Due to the wide range of lead fractions obtained from different methods, we conclude that the specific lead fraction product must be chosen depending on research question. Divergence- and LKF-based lead fractions provide valuable information in addition to established lead fraction products at high spatial resolution and independent of cloud coverage.