Overview of the MOSAiC expedition: Snow and sea ice

Sea ice and its snow cover make the Arctic Ocean distinct from most other oceans; they control energy transfer and important interactions between the atmosphere and the ocean. Therefore, they require a realistic representation in Arctic climate models, which are needed to improve our understanding of causes and consequences of the currently observed dramatic changes of the Arctic climate system (Thoman et al., 2020). Although model projections are continuously improved and becoming more consistent in forecasting the further decline of the Arctic sea ice cover, they still differ considerably on the rate of the decline (AMAP, 2017).

After the pioneering drift of Fridtjof Nansen with his vessel Fram (1893–1896; Nansen, 1897), numerous expeditions were performed to study the properties and interactions of Arctic sea ice and its snow cover with the atmosphere and the ocean. Most of these studies were internationally coordinated activities using icebreakers. Furthermore, 41 Russian North Pole drifting ice camps were established between 1937 and 2015 (Frolov et al., 2005). Only a few studies were conducted observing the evolution of sea ice and snow cover during different seasons or over a full annual cycle. During the International Geophysical Year (1957–1958), sea ice heat and mass balance measurements were made at the Ice Station Alpha in the Beaufort Sea (Untersteiner, 1961). The Arctic Ice Dynamics Experiment focused on understanding ice mechanics and dynamics during the yearlong drifting ice camp in the Beaufort Sea from 1975 to 1976 (Untersteiner et al., 2007). The international Surface Heat Budget of the Arctic Ocean (SHEBA) field campaign was a yearlong drift of the Canadian Coast Guard research icebreaker Des Groseilliers in the Beaufort and Chukchi Seas from 1997 to 1998. SHEBA focused on understanding the ice albedo and cloud radiation feedback mechanisms to improve climate models (Perovich et al., 1999). During the International Polar Year 2007/2008, snow and sea ice processes were studied along the drift of the schooner Tara (Gascard et al., 2008; Nicolaus et al., 2010; Haas et al., 2011). Between 2002 and 2008, Canadian Arctic programs involved the overwintering of the Canadian Coast Guard research icebreaker Amundsen to study the role of sea ice (Fortier and Cochran, 2008; Barber et al., 2010). The multidisciplinary Norwegian young ICE (N-ICE2015) expedition using the Norwegian research vessel Lance was carried out from January to June 2015. N-ICE2015 focused on the transition processes from multiyear to younger and thinner sea ice in the Arctic Ocean. It spent nearly 6 months studying atmosphere-ice-ocean-ecosystem interactions in an ice pack dominated by relatively thin (<1.5 m) first year ice (FYI) and second year ice (SYI; Granskog et al., 2018).

Beyond these crewed programs, autonomous drifting buoys, such as those coordinated by the International Arctic Buoy Programme (Rigor et al., 2008), observe annually the atmospheric, snow, sea ice, and ocean properties of the Arctic sea ice. Sea ice observations and data from airborne campaigns build another important source of Arctic snow and sea ice data (e.g., Haas et al., 2010; Kwok et al., 2017; Lawrence et al., 2018; MacGregor et al., 2021), supplemented by sporadic measurements collected by adventurers (Gerland and Haas, 2011). Numerous satellite observations provide an Arctic-wide perspective of sea ice properties including extent, concentration, thickness, type, drift, snow cover, and others. Passive microwave satellites have provided a sea ice climatology over more than 40 years, but they lack details for many ice-related processes (Kern et al., 2019). To establish new satellite products, coordinated ground truth observations are required (Gerland et al., 2019; König et al., 2019).

As a result of past expeditions, the understanding of the role of snow and sea ice in the Arctic and its influence on global climate has improved in recent decades. However, there are still significant knowledge gaps (Webster et al., 2018; Gerland et al., 2019) in quantifying the interactions and feedback mechanisms between the physical, ecological, and biogeochemical processes. Furthermore, an accurate representation of sea ice effects on biology and biogeochemistry needs to be better understood. To contribute to solving these problems, the work of the MOSAiC snow and sea ice team (in short ICE team) aims to achieve 5 overarching goals:

• Characterize the properties of snow and ice cover and understand the processes that govern these properties,

• Determine the snow and sea ice mass and freshwater balances,

• Quantify the partitioning of solar radiation between the snow, the sea ice, and the ocean,

• Describe the spatial variability and temporal evolution of the snow and ice cover, and

• Integrate snow and sea ice measurements with the atmosphere, ocean, and biosphere of the coupled Arctic system.

To achieve these main goals, a coordinated, integrated, and interdisciplinary approach was implemented as part of the MOSAiC expedition. Section 2 introduces the MOSAiC expedition and explains the concept of the snow and sea ice research program. The work program and a synopsis of the resulting data sets are presented in Section 3. In Section 4, the methodological and conceptual advances, and exemplary results from 2 case studies during 1 week of winter and 1 week of summer observations are shown. They result in a general description of the seasonal changes of the ice pack and the associated challenges of consistent and integrated observations of the coupled system components. Finally, linkages to the different Arctic subsystems, numerical modeling, and satellite observations are discussed in Section 5. Section 6 provides the major conclusions of this paper.

The MOSAiC field measurements started with the departure of the German research ice breaker Polarstern (Alfred-Wegener-Institut Helmholtz-Zentrum für Polar- und Meeresforschung, 2017) from Tromsø, Norway, on September 20, 2019, and ended with the arrival of the vessel on October 12, 2020, in Bremerhaven, Germany. Key dates are summarized in Table 1 and a map is shown in Figure 1. The observational year was divided into 5 legs: Leg 1 included the setup of the first Central Observatory (CO1, ice camp plus installations on Polarstern) and was supported by the Russian research ice breaker Akademik Fedorov. Akademik Fedorov performed most of the deployment of the Distributed Network (DN). The DN consisted of a hierarchy of autonomous systems deployed on ice floes surrounding the CO at 3–40 km distance to measure atmospheric and oceanic lateral gradients and increase sampling the highly heterogeneous ice pack properties. The network consisted of more than 100 simple position nodes, 8 nodes of medium instrumentation (M-sites), and 3 nodes of large sensor suites (L-sites). A more detailed description of the DN is under development, led by B. Rabe. The winter Leg 2 and spring Leg 3 continued the work on CO1, before Polarstern had to end the (manned) Drift1 and leave the floe, for logistical reasons, on May 16, 2020. The autonomous stations on CO1 and in DN1 continued the (unmanned) drift and recorded data while the vessel performed logistical operations in Svalbard. Afterward, Polarstern returned to the original ice floe but at a different location some hundred meters away. Leg 4 continued the drift with the new CO2 over summer until the decay of the floe in Fram Strait (between Greenland and Svalbard; Figure 1) on July 31, 2020. This ended the Drift2 and only a few autonomous devices continued the drift of CO1, CO2, and DN1. After the final rotation of personnel, Polarstern travelled back into the ice and started Drift3 with the setup of CO3 and DN2 on August 21, 2020, near the North Pole. Drift3 ended on September 20, 2020, when the vessel started the return voyage. A publication with details on the expedition in general and its logistics is under development led by M. Rex.

To achieve the main goals of MOSAiC, a coordinated, integrated, and interdisciplinary approach was implemented, as a collaboration of 5 scientific teams: the atmosphere (ATMOS; Shupe et al., 2022) team, the snow and sea ice (ICE, this manuscript) team, the oceanography (OCEAN; Rabe et al., 2022) team, the ecosystem (ECO) team, and the biogeochemistry (BGC) team. The work on physical snow and sea ice processes and properties was organized in 13 tasks, which are described in the subsections of this chapter (see also Text S1 with Figure S1). All tasks comprise elements of observations and numerical modeling, while the balance between these elements differs between the tasks. Additional details on some methods are given in Text S2, including a list of short names and acronyms in Table S2:

1. General snow and ice observations

2. Snow measurements and sampling

3. Ice coring

4. Mass balance observations

5. Terrestrial laser scanning (TLS)

6. Optical measurements

7. Melt pond observations

8. Sea ice dynamics

9. Ridge observations

10. Transect measurements

11. Remotely operated vehicle (ROV) observations

12. Helicopter observations

13. On-ice remote sensing (RS).

The snow and sea ice measurements were distributed over the entire CO (Figure 2). A detailed description of all 3 COs with additional maps and details on installations is under preparation by the project team led by M. Nicolaus. While the setup of the CO changed between legs due to ice dynamics and ship relocation, its general layout is exemplified in the map in Figure 2, which shows the CO1 ice camp shortly before the temporary departure of Polarstern at the end of Leg 3.

First results with a focus on snow and sea ice conditions have been published, for example, on the history of the ice drift prior to the MOSAiC expedition (Krumpen et al., 2020), on the ice conditions from satellite observation along the drift track in previous years (Krumpen et al., 2021), on the representativeness of the selected ice floe (Belter et al., 2021), on the linkages to the general atmospheric circulation (Dethloff et al., 2021), on sea ice RS methods (Munoz-Martin et al., 2020; Stroeve et al., 2020; Semmling et al., 2021), and on the platelet ice accumulation in winter (Katlein et al., 2020).

The general snow and ice observations and specific activities were performed during all legs to document the general snow and ice conditions and to support data interpretation afterward. Mapping activities of discrete measurements, installations, and operations on the ice were also organized in this task. Its activities included:

1. Acquiring time-lapse photos of the surroundings of the vessel with a 360-degree panorama (Panomax) camera mounted above the crow’s nest. Each photo has a resolution of 15,680 × 2,048 pixels and photos were taken at 20-minute intervals (72 photos per day).

2. Recording standard sea ice observations daily (on transits hourly) from the bridge following the protocol of the Arctic Shipbourne Sea Ice Standardization Tool (Hutchings and Faber, 2018). The aim of these observations was to describe the surrounding conditions over a 10-minute period within a radius of 1.5 nautical miles around the vessel, including the polar night.

3. Documenting the on-ice work by establishing a local floe-based x, y coordinate system to provide the relative locations of the various sites (FloeNavi system), because measurements on sea ice have the particular challenge that the ice floes are in constant movement. In addition, geographic positions were recorded to allow reference to any other measurements, in particular airborne and satellites, that are not referenced to the floe.

4. Acquiring visible and thermal photography from helicopters and drones to provide information about sea ice surface conditions, to support mapping of the main floe, and to generate digital elevation models of the surface. Visual floe maps and digital elevation models were created by stitching photographs from single flights. The final images were geo-referenced to the local x-y coordinate system.

5. Upgrading, before the start of the drift expedition, the Electronic Charting System on board Polarstern, for receiving and displaying various satellite and weather data in near real time. The effort was focused on the charting of weather forecast products provided by the German Weather Service. In addition, model-based ice drift forecasts were included, as well as a combination of various new satellite products.

6. Documenting the on-ice work almost daily on board and synchronizing with land as the “MOSAiC LogBook” to inform the project consortium on details of the ongoing work.

The overarching goal of the snow measurements was to characterize the spatial and temporal variability of snow over different sea ice types (e.g., FYI, SYI, and ridges) during an entire annual cycle. During the melt season, techniques used to investigate the snow were applied to the uppermost layer of the melting ice cover, the so-called “surface scattering layer” (SSL). The snow measurements during MOSAiC were conducted with the aim to improve snow models, and subsequently to improve climate models and RS retrieval algorithms, and in support of all MOSAiC teams. Snow is one of the most insulative naturally occurring materials on Earth (Webster et al., 2018), and one of the most efficient reflectors of incident solar radiation (Warren, 1982). Snow on sea ice therefore acts as an effective barrier between the atmosphere and the ocean, preventing both heating of the upper ocean in summer and cooling of the ocean in winter. While the bulk physical properties of snow on sea ice (e.g., snow depth, snow density, snow thermal conductivity, snow reflectivity, albedo) and the resulting stratigraphy have been studied previously, data on the spatial heterogeneity of the optical and thermal properties of snow on sea ice and the relationship between these properties and their governing processes are limited. Snow measurements during MOSAiC ranged from microscale (micrometer) to local scale (approximately 1000 m) to address these properties across scales. The following novel instruments were used during MOSAiC to measure physical properties of snow:

• –  The Snow Micro Pen (SMP), a penetrometer (Proksch et al., 2015) for fast retrieval of the vertical profile of snow stratigraphy, snow density, and specific surface area, with a high vertical resolution (better than 1 mm). The SMP measures vertical profiles of snow layers at one location. By repeating samples over a large area, micro-, macro- and regional-scale properties can be linked. SMP was used on snow profiles and transects. Overall, 6,959 SMP profiles were measured.

• –  An X-ray microtomograph (micro-CT, Scanco Medical micro-CT 9; Calonne et al., 2014), which was installed on board Polarstern and allowed measurements of snow microstructure on quasi in situ samples. We used the micro-CT to measure the full snow profiles throughout the year, as well as on the SSL, lead and pond ice, and sea ice cores. Overall, 614 samples were analyzed with the micro-CT.

• –  Near Infrared (NIR) photography (Matzl and Schneebeli, 2006) to characterize the specific surface area of snow profiles as well as the snow and SSL surfaces. During MOSAiC, photographs at wavelength of 850 nm and 940 nm were taken for each profile/surface scene. Overall, some 600 profile/surface scenes were photographed using NIR.

• –  Structure from motion photograph sets to derive elevation data. Based on these data, high spatial resolution digital elevation models were created at each snow pit location, and sometimes along transects, of the surface and of the snow–ice interface. These “mini digital elevation models” will be used to determine small scale (millimeter) surface roughness that is relevant for optical modeling and satellite retrievals, and link data to coarser-scale digital elevation models obtained by terrestrial and airborne laser scanning. Overall, some 300 structure from motion sets were taken during MOSAiC.

Snow observations were collected mainly in specific areas assigned to snow measurements. Those areas included the prevalent types of sea ice, namely FYI (e.g., Snow1), SYI (e.g., Snow2), ridges, and new ice in leads, as highlighted by yellow patches in Figure 2. Additional measurements were performed across the COs to connect physical snow properties to relevant processes that were studied by other teams. The team members conducted about 600 snow pit measurements of different types:

• –  “Type C” snow pits characterized the vertical profile of the snowpack and consisted of SMP measurements only.

• –  “Type B” snow pits included density, snow water equivalent (SWE), temperature, and salinity measurements in addition to observations of Type C.

• –  “Type A” snow pits included the full suite of snow measurements, adding snow sampling (see below), micro-CT scans, SMP transect (to either side of pit), NIR photography, structure from motion pictures, and snow sampling for a wide range of constituents, in addition to observations of Type B and C.

For practical reasons, this perfect scheme of snow pit types was not always followed in every pit. We performed a comprehensive snow sampling program, which was centrally organized and coordinated across all scientific teams. Snow samples colocated with the physical measurements were taken for analysis of salinity, water and sulfur stable isotopes, major ions, halogens and halocarbons, beryllium, black carbon, marine sugars, organic acids, microplastics, and ice-nucleating particles. Snow sampling was combined with snow pit measurements whereby the frequency and number depended on the purpose of the snow sample. This approach ensured that all samples from each sampling event were taken at the same site(s) and in a consistent and reproducible manner facilitating the colocation with the physical properties of the snow. In addition, dielectric permittivity measurements were made at 50 MHz to invert snow wetness during melting conditions.

Figure 3 shows the temporal distribution of snow pits. Distributed observations of physical snow properties on transects and of snow depth (Section 2.10), and observations of snow surface position with terrestrial and airborne laser scanner (TLS Section 2.5 and ALS Section 2.12) will be used to upscale these point observations. Collectively, these measurements provide a comprehensive characterization of the snow cover, enabling synthesis studies on snow evolution, spatial variability, energy budget and heat fluxes, RS algorithm development, and toward improved parameterizations in climate models.

The primary aim of the ice coring work was to capture the seasonal evolution of sea ice physical properties over a full year, for both level FYI and level SYI at 5-cm vertical resolution, and to document the differences between FYI and older ice. This work was performed at 2 main coring sites (Figure S2 with Tables S2 and S3), one on FYI and the other on SYI, which was the oldest ice at the start of MOSAiC. We used the approach recommended by Oggier et al. (2020) to section the ice cores, relative to both the ice surface and ice bottom (Section 4.2), to better capture the processes at the surface and bottom interfaces. The physical properties affect energy and matter fluxes and ecosystem processes in sea ice and are also important to the RS of sea ice (Section 2.13). The ice coring work captures the seasonal evolution of level ice thickness, temperature, salinity, and density. At each coring event, temperature profiles were measured using a thermistor probe into holes drilled at 5 or 10 cm intervals within 2–3 min of core extraction. Two cores collected within 50 cm of the temperature core were cut in sections and stored in sealed containers for transport to the ship, after noting total length and fractures. Onboard, bulk salinity was measured on melted ice samples, and density was measured using the hydrostatic density method (Pustogvar and Kulyakhtin, 2016). The brine volume fraction can be derived from these measurements. Two additional cores were collected for textural and microstructure analysis: one for thin sections and a second for micro-CT scans (Section 2.2). Bottom and surface sections (in most cases the lower and uppermost 10 cm of ice cores) were imaged on board using the micro-CT on Polarstern. The remaining ice cores were archived for more detailed studies of texture and microstructure at home laboratories after the expedition. In the sunlit period, ice cores were also collected to measure the optical properties at the home lab (Section 2.6). Coring was coordinated with the snow task, and often the coring was associated with a simple snow pit (which typically included temperature, density, salinity, and sampling for oxygen isotopes; Section 2.2) supporting the interpretation of the ice core data. Coring was conducted in close collaboration with the ECO and BGC teams.

In addition to the time series data from the main coring sites (Figure S2), the BGC team observed the temporal evolution of both FYI and SYI at a number of additional sites (primarily salinity and temperature) coincidently with their own work. In addition, sediment laden ice cores were collected for analysis jointly with the BGC team, and an occasional ice core was collected for analysis of microplastics. These ice cores complement the observations made by the ICE team and provide more information on spatial variability. Coring was also done repeatedly in ridges (Section 2.9) and for ice mechanical properties (Section 2.8). Opportunistic sampling was done especially for thin young ice in leads and linked to RS work during some events and occasionally in the DN. The collected data set captures the seasonal evolution of FYI, from early growth and new saline ice to mature FYI and the rapid desalination in summer (Section 4). While the residual FYI that had survived the previous summer melt became SYI at the end of the year, the refreezing and new ice growth at the bottom of SYI could be observed until its summer decay. An overview of the ice coring activities by the ICE team is shown in Figure 3, and a list of existing ice cores is referenced in the data availability section.

The goal of the mass balance observations was to document the amount and spatial variability of ice growth during winter and surface and bottom ice melt during summer. Mass balance observations provide insight on the relative contributions of the atmosphere and ocean to changes in the sea ice thickness and snow depth. Multiple sites sampling new ice, FYI, SYI, deformed ice, undeformed ice, and ponded ice were included to comprehensively address the spatial variability of sea ice mass balance. The in situ mass balance measurements were made both autonomously and manually with 37 digital temperature chains, 4 sea ice mass balance buoys (IMB) of different types (e.g., Polashenski et al., 2011; Jackson et al., 2013), 20 Snow Buoys (Nicolaus et al., 2021), and 120 manual ice thickness gauge/ablation stake pairs. The temperature chains and IMBs recorded vertical profiles of temperature of the air through the snow and ice into the upper ocean, as well as (for some IMBs) the positions of surface and bottom interfaces using acoustic rangefinders. Analysis of temperature chain data manually and by algorithms are also used to determine the position of the air–snow or –ice and the ice–ocean interfaces (Hoppmann et al., 2015; Cheng et al., 2020). Manual measurements at the ablation stake sites also provided a time series of the sea ice surface and bottom position changes. Thicknesses measured from ice cores collection and laser scanning measurements contributed to the snow and sea ice mass balance measurements (Sections 2.3, 2.5, and 2.12). Together, these multiple methods of observing ice surface and bottom position are used to calculate time series of snow depth and ice thickness, ice growth, surface melt, and bottom melt. In combination with transects (Section 2.10) this approach provides a comprehensive view of the sea ice mass balance (in space and time) and the best local (thermodynamic) and best spatial (dynamic) information. To complement the mass balance observations, systems for atmosphere and ocean measurements were colocated on the same sites.

The goal of the TLS was to map the snow and ice surface topography and how it changed in time and space throughout the seasonal cycle. Using the eye-safe Riegl VZ1000 TLS during winter and the non-eye-safe Riegl VZ6000 TLS during summer allowed large aerial scans. Analysis of these data is being used to link these changes with other in situ measurements of snow and ice properties and spatially extend point measurements of snow and ice properties. TLS measurements were used to map an area of approximately 0.5 km² every 2 weeks with a centimeter-scale vertical resolution and decimeter-scale or better horizontal resolution. The same region was mapped with 42 scans during the drift with CO1 from October 18, 2019, to May 9, 2020, (Figure 3) and included most of the snow pit measurement sites, most of the mass balance observation sites as well as the southern loop of the transect, all of the albedo transects in spring, portions of the strain gauge installations, the area studied from beneath the ice at ROV sites (ROV 2.0 and 3.0 in Figure 2), the on-ice RS instruments, and the on-ice meteorological installations at Met City (Sections 2.2, 2.4, 2.6, 2.11, and 2.13). This colocation of the various methods will allow linked data processing, while the impact of surface destruction on the TLS results through other activities is minor given the very small fraction (<1%) of impacted surface area. Summertime TLS observations were more limited due to safety and logistics but included areas of snow pits, albedo, ponds, and other surface property measurements. The TLS observations are used to generate a time series of ice dynamics, changes in atmospheric form drag, and snow accumulation, redistribution, and melt. TLS measurements can be used to infer material properties throughout the mapped region by upscaling the vertical profile measurements of snow and ice physical properties made by other instruments.

The objective of the “optics” task was to quantify the partitioning of solar radiation within and beneath the ice cover during the sunlit season. The quantities measured included albedo, transmittance, and light extinction in snow, sea ice, and the underlying water. A variety of sensors (measuring either broadband, spectral, or photosynthetically active radiation) were used for measuring repeat-visit survey lines and grids in addition to fixed-point study sites. The albedo measurements included surface survey lines, multiple fixed site installations, unmanned aerial systems (e.g., the drones HELiX and Mavic and the Spectra Quadrocopter), helicopter-borne instrument platforms, and the tethered balloon. Transmittance measurements included repeat under-ice surveys using an ROV (Section 2.11) and site-specific measurements using an L-arm photosynthetically active radiation sensor. In-ice extinction and inherent optical property measurements were carried out at specific locations using profiling sensor arrays frozen into the ice (light chains and light harp) and a spatially resolved reflectance probe lowered into auger holes. Light extinction within the water column beneath the ice cover was estimated using sensors mounted on the ocean CTD-sensor (conductivity, temperature, depth), and the ROV was also equipped with a number of optical sensors.

Sensors were also deployed on radiation stations to autonomously collect albedo, transmittance, and water column radiation at fixed locations throughout the season. In addition to these time series, specific experiments were carried out to assess the response of spectral albedo to specific events or processes: (1) snow during early summer melt, (2) diurnal variability during both melt and freeze-up, (3) the removal and formation of a SSL, and (4) rain on new snow. Additionally, ice was sampled for investigation of the linkages of optical and structural properties to be carried out in freezer laboratories at home institutions (Section 2.3).

Ongoing investigations of the spatial heterogeneity of the surface albedo at various horizontal scales will be aided by the drone-based spectral and broadband albedo measured along grids and horizontal transects at various heights above the ice surface. The vertical profiles of albedo recorded from drones will be crucial for determining the contributions of various surface features (melt ponds, ridges, flat-ice, leads) to the areal-averaged albedo. These measurements enable scaling from ground-based observations (meter-to-tens-of-meters footprint), to satellite observations (tens-to-hundreds-of-meters footprint), and large-grid model outputs. Moreover, high spatial resolution hyperspectral albedo measurements collected from drone-mounted cameras will help relate the scattering properties of each surface feature to the snow and ice microstructure properties observed at the surface.

Melt ponds play a key role in energy exchange between atmosphere and ocean in the Arctic summer because they substantially reduce surface albedo of sea ice. Compared to bare ice, melt ponds reflect less solar energy to the atmosphere and transmit more solar energy into the ice underneath and the ocean, thus fostering further sea ice melt (Light et al., 2008; Nicolaus et al., 2012). At present, a synoptic view of pond evolution, distribution, and depth using airborne and satellite data is still limited. The main objective of the Ponds task was to provide field and airborne data (Section 2.12) of the physical pond characteristics (e.g., size and shape distributions) to improve the understanding of the spatiotemporal development during all stages of pond evolution. These characteristics include pond area fraction, size, shape, and number. To improve preprocessing of airborne and satellite data (i.e., atmospheric correction, correction of adjacency effects; König et al., 2019), the spectral signature of the different surface types on the floe were measured using hand-held water color spectroradiometers, a hyperspectral camera, and a goniometer. Up- and downwelling broadband irradiance data as well as multispectral imagery are available from drone surveys, helicopter flights, and satellite observations (Sections 2.6 and 2.12). To assist a proper reflectance calibration of water surfaces, we used a remotely controlled boat equipped with water color spectroradiometers. Above and below water spectral measurements were used to assess the optical properties of pond water (attenuation coefficient, inherent optical properties; König and Oppelt, 2020). In addition to the radiometer measurements, the field measurements included pond depth as well as chlorophyll content, salinity, and temperature of pond water. To assess the absorption coefficient of the dissolved and particulate matter, water samples were taken and analyzed in the Polarstern lab using a point-source integrating-cavity absorption meter. Micro-CT samples were processed from frozen melt ponds to derive optical properties.

One of the key questions of MOSAiC is how sea ice moves and deforms. The dynamics and mechanics of the ice around the Polarstern and its interaction with the ship, ocean, and atmosphere were monitored throughout the drift. We aimed to examine the spatial and temporal distribution of ice motion and ice stress, and to determine how they vary locally, regionally, and seasonally in relation to the physical properties of the ice. The CO observations were coordinated with satellite and airborne RS (Section 2.12) monitoring ice drift, ice morphology, and floe size. Drift and strain-rate (divergence, shear, and vorticity) of the ice pack were measured over multiple nested scales, from hundreds of kilometers in a wide array of buoys coordinated with the International Arctic Buoy Program, to ocean synoptic and mesoscale scale in the DN measured with more than 100 ice-deployed GPS position sensors, to ship-mounted ice radar (Oikkonen et al., 2017) providing dense target tracking over 7 km around the ship every few minutes, to high frequency and centimeter-level position accuracy from GNSS helicopters inertial navigation system stations and robotic theodolite strain observations. Ice motion was also monitored at 2 locations with seismic stations, providing information on vertical motion, and ice–wave interaction (Marsan et al., 2011). Three seismic stations, located at 3 points of the ice floe, recorded ice oscillations in 3 directions (horizontal and vertical) and transmitted the records to the central station installed on the ship. The signal recording frequency was 100 Hz. The graphical format of the waveform records allows to see ice events such as new fractures, compression or ridging. Internal ice stress was recorded with vibrating wire strain gauges (Cox and Johnson, 1983) at 17 sites in the CO. Twelve of these gauges had 1-minute sampling, the others 10-minute sampling, allowing for specific breaking events to be aligned with the stress time series. To estimate the ice stress on the Polarstern, stress panels were placed in the ice about 20 m from the ship parallel to “Void 92,” a specific compartment in the hull of Polarsternice-, which was instrumented with strain gauges. Throughout the drift, the mechanical properties of the ice (local strength, horizontal, and vertical uniaxial strength tests) were measured at a level SYI site. For these measurements, an ice-hydro-complex with a probe-indenter and a press was used. The probe-indenter was used to measure the local strength of the ice cover at different levels, and the press was used to measure the ice strength for uniaxial compression (horizontal and vertical). Salinity, temperature, density, and texture analysis of colocated ice cores allow characterization of the mechanical properties of ice (Section 2.3). A log of specific dynamic events (cracks, leads opening, ridging) was kept; Panomax camera and Ice Radar were used to identify the start time of each event to within 20-min precision. The full data set allows for specific case studies that can correlate ice fractures and the response of ice to ocean and wind forcing on scales beyond tens of meters.

Ridges are one of the least studied features in the Arctic ice pack because of logistical challenges in making measurements in these sometimes-massive chaotic formations of ice that pile up when ice floes collide. The main aim of the ridge work was to study the development of the geometry and of the thermomechanical properties of ridges during the MOSAiC drift. The 3 main phases of ridge evolution, initial formation, winter consolidation, and summer decay, were documented for the first time for selected ridges throughout a full seasonal cycle and compared to the evolution of level FYI and SYI. The main geometrical parameters characterizing the ridges were the thickness of the consolidated layer, the keel depth, and the macroporosity. These 3 parameters were determined by repeated manual drilling (Figure 3) and complemented by continuous temperature profiles obtained by IMBs and thermistor strings that were deployed vertically through specific ridges. Additionally, sea ice cores were taken from ridges to characterize the temporal and spatial evolution of sea ice temperature, salinity, and density. In addition to this baseline work, additional specialized measurements were performed to determine bottom topography, ocean properties, currents and keel-associated turbulence, as well as surface characteristics by snow and ice thickness transects, snow surveys, and laser scans of surface topography (Sections 2.2, 2.3, 2.5, 2.10, and 2.11). The characterization of the physical properties of ridges was of particular importance to coordinated ecological studies of sea ice ridges performed by the ECO team.

Making systematic and representative observations of physical sea ice parameters, their spatial variability, and temporal evolution in the CO throughout the full MOSAiC drift was the objective of the transect task. The sampling strategy consisted of repeat surveys of sea ice and snow thicknesses throughout all MOSAiC legs, similar to observations carried out from May to December during the Tara expedition in 2007 (Haas et al., 2011). These baseline observations were complemented with seasonal activities, specifically measurements of surface albedo and melt pond properties during summer and autumn. In addition, cross-task activities with joint measurements by mobile on-ice remote-sensing sensors (Section 2.13; L-band microwave radiometer and Ku/Ka-band radar) along the transect lines or transect observations at specific ridges (Section 2.9) were carried out to provide input data for algorithm development and temporal and spatial evolution at ridges.

The transect measurements were carried out by repeating dedicated loops typically 1 km or longer, except for Leg 5 which was shorter. The baseline observations included a broadband electromagnetic induction sensor (Geophex GEM-2, Hunkeler et al., 2016) for the measurement of the total thickness of the sea ice and snow/SSL layers and an automated snow/SSL depth probe (Snow Hydro MagnaProbe; Sturm and Holmgren, 2018) which in combination yield point measurements of snow/SSL and sea ice thicknesses. To accommodate measurements of snow/SSL and melt pond depths when both were present, the snow disc of the MagnaProbe was outfitted with floatation on its upper side to also measure pond depth. During winter, only limited FYI and mixed type transects were performed, as no observations on FYI were possible in early winter. These transects were not revisited on later legs as they became inaccessible due to deformation. During summer and autumn, additional transects were conducted parallel to albedo lines for relating the properties of snow, SSL, sea ice, and melt ponds to the spectral and broadband albedo measurements (Sections 2.6 and 2.7).

The 2 goals for the transect surveys were to achieve representative coverage of sea ice conditions in the CO and continuity of observations throughout the drift. However, due to the relocation of Polarstern after Legs 3 and 4 (Table 1), new transect loops were established in addition to the original transect loops. In the beginning of the MOSAiC expedition, the 2 original loops were laid out, termed “northern” and “southern” loop. The rationale behind the 2 loops was to survey the 2 dominant ice types in the CO: the northern loop on the older, thicker, and more deformed part of the “MOSAiC floe”; and the southern loop on the younger and thinner part that had widespread refrozen melt ponds. While these 2 loops were traversed continuously from October 2019 to early May 2020, the southern loop was often inaccessible after sea ice deformation events. The summer transect (CO2) was established in June 2020 on similar ice conditions as the prior transect loops and with small (approximately 300 m) overlap with the northern transect loop. The transect loop of Leg 5 (CO3) was fully independent of earlier transects due to the change in regional focus. A full overview of the dates, location, and acquisitions of all transect observations is referenced in the data availability section. In addition, the dates of all total ice plus snow thickness observations including data acquisitions outside the transect task are illustrated in Figure 3.

The main goal of the work with the ROV was to obtain a better understanding of the seasonal cycle with respect to (1) ice draft and bottom topography, (2) the light field beneath the ice, (3) the biophysical properties of the ice and uppermost ocean, and (4) the organisms living in and under the sea ice. All measurements focused on the spatial variability under sea ice and across different surface features and ice types. To accomplish these objectives, the ROV was operated 1–3 times per week amounting to a total of 83 dives. The ROV was equipped with a comprehensive suite of biophysical sensors, several cameras, and a mechanism to haul sampling nets under the ice (Katlein et al., 2017; Wollenburg et al., 2020). During summer, additional dives were performed with an underwater hyperspectral imager (UHI) to survey the 2D and 3D radiance under sea ice. In addition to the direct data sets from the ROV, the work was complemented by data sets from other ICE and ECO team tasks to allow the upscaling of point observations, observing links to ecosystem parameters, and estimating energy budgets combined with other observations.

The ROV was operated within an approximate 300-m radius around the ROV sites (Figures 2 and S3). Over the year, the ROV site was moved several times due to the dynamic icescape. Specific dive missions were defined and executed as part of the weekly schedule. Over most of the year, weekly multibeam sonar surveys and net hauls were carried out. Multibeam sonar surveys mapped the under-ice topography in a grid covering an area as large as possible at a depth of about 20 m (Section 4.1). Net hauls were conducted along linear or triangular transects at different depths for dive times of approximately 15 minutes. Once sufficient sunlight returned (mid-March), comprehensive optical dives were carried out on the ROV days (including downwelling radiation), in particular for grids and transects under FYI, SYI, leads, and ridges. This effort also included stable measurements at marker positions, vertical profiles, and measurements of upwelling radiation. Additionally, specific dive missions were added to locate lost devices, to document under-ice installations and to deploy and recover sediment traps.

The main objectives of the helicopter-borne sea ice surveys were, first, to document the temporal evolution of snow and sea ice cover on the CO floe and within the DN, from sea ice formation in autumn to sea ice melt in summer and, second, to document the spatial distribution and spatial variability of snow and sea ice properties within the DN and beyond. Key parameters were sea ice thickness distribution, surface topography, freeboard and floe size distribution, surface temperature and surface albedo, areal fraction of different sea ice types, and melt ponds as well as additional melt pond characteristics like pond-size distribution and bathymetry (Section 2.7).

To reach these goals, 3 different sensor packages were operated independently on different flights: (1) an airborne laser scanner (ALS, Riegl VQ-580), 2 RGB cameras (Canon EOS 1D Mark III) with wide-angle and fisheye lens, in winter an infrared (IR) camera (VarioCAM HD head 680, InfraTec) or in summer a VIS/NIR hyperspectral camera (Specim AisaEAGLE, 400 nm to 970 nm), a radiation thermometer (Heitronics KT19) and an upward and downward looking pyranometer (Kipp&Zonen CMP22), (2) EM-Bird for measuring the combined thickness distribution of sea ice and snow using airborne electromagnetic induction sounding (Haas et al., 2009), and (3) the towed system HELiPOD, which was operated by the BGC team for measurements of trace gases, energy and momentum fluxes, aerosol and sea ice surface properties. The laser scanner and camera configuration was operated during all legs on a total of 68 flights, the EM-Bird was operated during all legs, except in the midwinter darkness, on 23 flights, and the HELiPOD was operated in spring and summer on 5 flights.

There were 5 main mission types to investigate the temporal evolution and spatial variability of snow and sea ice properties: (1) local floe grid survey covering the CO floe and surrounding floes, (2) regional survey consisting of lines between the 3 L-site buoys and between these buoys and Polarstern, (3) butterfly survey consisting of triangles between Polarstern and 4 selected buoys of the DN to observe the relation between ice deformation and thickness change, and (4) mixed pattern of horizontal lines and vertical profiles mostly across the CO. Missions (1) and (4) were flown with the ALS and camera sensor package as well as with the EM-Bird and HELiPOD, mission (2) with the ALS and camera sensor package only, and mission (3) with the EM-Bird only. The maximum airtime of individual flights was about 2 hours.

Almost all Arctic-wide sea ice climate time series are based on microwave satellite observations (IPCC, 2013), which now span more than 40 years. However, the retrieved geophysical parameters like ice concentration, ice type, ice thickness, melt onset, or snow depth are not measured directly by the microwave radiometers and radars but inferred from the interaction of microwave radiation with sea ice and snow as well as seawater. Depending on frequency and physical properties (temperature, salinity, snow grain size and geometry, porosity, etc.) of the snow and sea ice (including melt ponds), microwave radiation can penetrate into the ice, get scattered, absorbed and emitted. Information of the microwave radiative transfer in snow and sea ice is needed to derive the observational quantities from satellite measurements.

The main objective of the ice RS task is to develop new methods to retrieve improved sea ice and snow parameters from satellites. This goal relates both to reducing and quantifying the uncertainties of existing satellite retrievals (e.g., ice concentration, sea ice thickness, and ice type) and to developing new methods for current and future satellite missions (e.g., snow depth for CRISTAL and CIMR, surface properties and ice types from Radarsat Constellation Mission and TerraSAR-X, melt pond depth using hyperspectral airborne data and Sentinel-2). Most measurements described here were in the microwave domain from 0.5 to 89 GHz using passive microwave radiometers and active radars. Here, the most pressing issue is to better understand how the microwave emissivity and backscatter depend on physical snow and ice properties and environmental conditions. Thus, RS instruments were codeployed with snow and ice measurements (Figures 2 and S4). In turn, radiative transfer modeling of snow and ice will be improved with the availability of the RS measurements. Further observations were carried out in the visual and IR domain to support the microwave measurements and to measure ice, snow and melt pond optical properties as well as surface temperature distribution and evolution. MOSAiC offered the unique opportunity to observe a full seasonal cycle and colocate measurements at different frequencies and polarizations. RS measurements were taken at a central RS site on the ice floe during all 5 legs, from the ship, and along transects (Section 2.10). An overview of the RS site is given in Figure S4. A detailed overview of satellite related work is under development led by G. Spreen. Sea ice conditions along the drift trajectory based on satellite observations and comparison to previous years are described in Krumpen et al. (2021).

Table 2 provides a short overview of the 14 RS instruments that were operated during MOSAiC, with their basic instrument parameters. In addition, examples of corresponding satellites or sensors are given. The 3 main instrument categories were microwave radiometers, radars, and other sensors like GNSS-R, IR, and cameras. These RS measurements on the ice correspond to a wide range of satellite sensors, which are also listed in the table. Due to ice dynamics and instrument problems not all instruments were operational all the time (see Figure 3B for observation periods). The primary measured parameters are emitted infrared and microwave radiation (brightness temperature from radiometers), backscattered (radars) and reflected (GNSS-R, nadir Ku/Ka radar) microwave radiation, and spectral reflectances (hyperspectral camera).

The work on snow and sea ice processes and properties required an integrated approach to carry out the observational program along the tasks outlined in Section 2. Performing consistent measurements throughout a full year was a major challenge and one of the main differences between MOSAiC and many previous snow and sea ice field studies. To guarantee consistent high-quality data collection between all legs, a particular preparation phase was performed (Text S1). This section describes the key processes, which were studied by the ICE team across the different tasks (Section 2). These processes are illustrated in Figure 4. While most processes were observed continuously over the full annual cycle, seasonal processes were only studied during parts of the drift (e.g., solar short-wave radiation interactions). The work covered multiple ice types and the highly variable snow distribution: new ice, FYI, and SYI. We worked on level, deformed, and ridged sea ice as well as in the marginal ice zone, leads, and melt ponds. In addition, we describe the measured parameters and the resulting data sets, which will be used to improve and develop parameterizations for a better understanding of the coupled Arctic sea ice system. The key parameters are summarized in Table 3 and their temporal coverage is shown in Figure 3. Methods and tasks are restricted to the work in the COs.

New ice formation in open water (mainly in leads) was observed through the entire autumn, winter, and spring, including observations of surface, bottom, lateral, and internal freezing. In October 2019, most measurements were performed on SYI that had survived the previous summer. However, parts of the ice camp, for example, the coring activities, were also established on younger sea ice that involved FYI. Thus, the evolution of both FYI and SYI during the initial freezing was observed and followed throughout the campaign. Furthermore, MOSAiC was able to observe the initial accumulation of snow in autumn. Snow accumulation, distribution, and redistribution were observed through the entire year with a particular focus on comparisons of snow processes between level FYI, level SYI, and deformed (ridged) ice in the CO and how these processes affect the growth of different ice types.

The novel observation of platelet ice under winter sea ice in the Arctic showed that platelet ice contributed to bottom ice growth during winter (Katlein et al., 2020). Over time, the ice pack became strongly deformed. After melt onset, the melt and decay of the snow and sea ice were observed, and contribution from surface, bottom, lateral, and internal melt processes were quantified. The decay of the original MOSAiC floe (CO1 and CO2) on July 31, 2020, showed strong linkages between dynamic and thermodynamic processes controlling the decay of the ice pack in the marginal ice zone.

In addition to the general time series of each ice type, the comparison of sea ice and snow mass balance for level and deformed ice was an important element along the transects. The thickness, consolidation, and decay of pressure ridges were examined by a combination of drilling, coring, and analysis of IMB and thermistor string data. The macroporosity of ridges was estimated from drillings, while the keel depth was investigated by drilling, ROV multibeam mapping, and IMB/thermistor data. Almost all tasks contributed in some way to the sea ice and snow mass balance work (Table 3).

The physical properties of sea ice were studied based on ice cores and samples on the microscale and also through surface and under-ice transects as well as larger scale airborne and dynamical studies (Table 3). The ice coring work of the ICE team focused on the time series of the main coring sites to capture the evolution of FYI, SYI, and deformed ice. Additional information on the spatial variability is available through the work related to the mechanical properties, the RS work and in particular through the BGC team, which surveyed additional sites across the COs and DN1. Optical properties were mostly studied in relation to energy budget studies (see below).

Snow on sea ice is the thermal and structural interface between sea ice and the atmosphere (Sturm and Massom, 2017). Thus, even though major parts of the snow work were performed by the ICE team, many linkages exist to the ATMOS team, and comprehensive parts of the snow work were carried out in strong collaboration with the ATMOS team (Shupe et al., 2022). Snow precipitation, rain fall, and snow redistribution were studied with common setups by both teams. Snow formation, snow metamorphism, and changes in the snowpack were studied in detail from micrometer scales to floe scales. Snow cover is the most (temporally and spatially) variable layer dominating the thermal and optical properties of the sea ice system. Due to the highly metamorphic nature of snow, its microstructure changed at a fast pace. Snow processes during MOSAiC were highly sensitive to quickly changing air temperatures, which led to strong temperature gradients in the snow, particularly during Legs 3 (spring) and 5 (autumn), and resulted in a high degree of snow metamorphism. The strong temperature gradients affected the thermodynamics, growth rate, and microstructure of the underlying sea ice and consequently strongly affected the temperature, salinity, and permeability of sea ice. Ultimately, the physical snow properties and processes build the link to light and nutrient availability in the sea ice and in the upper ocean, and significantly affect the polar marine ecosystem.

Incident solar shortwave radiation is partitioned into reflected, absorbed, and transmitted components during the sunlit season. This split determines the surface shortwave radiative budget, with effects on mass balance and physical properties of snow and sea ice, and contributes to ocean heat (Figure 4). The optical properties of the ice cover were investigated by a variety of methods. This variety was needed to constrain radiation budgets, provide some spatial and temporal overlap, achieve redundancy to help alleviate data gaps where instruments failed, and to provide insight into physics that are driven by or drive processes not well constrained by a single set of measurements. An example would be the utility of surveying both albedo and transmittance simultaneously to understand the full partitioning. While the work of the ICE team focused on studies of this partitioning within the snow and ice, the ATMOS team observed radiative transfer in the atmosphere as well as the full surface energy budget including radiative, turbulent, and conductive heat fluxes. The OCEAN team observed radiative heating and turbulent mixing in the ocean, using a combination of near-ice direct heat, salt, and momentum ocean flux measurements and profiling CTD and microstructure dissipation profiles (Rabe et al., 2022).

Conductive heat flux between atmosphere and ocean through the snow and sea ice pack was studied mostly based on thermistor chains (Section 2.4) and also as residuals from surface and basal mass balance calculations. These heat fluxes are crucial for quantifying melt and freeze processes at both interfaces and thus are key elements to close the atmospheric and oceanographic energy budgets. Heat fluxes through snow and sea ice strongly determine sea ice mass balance.

Fixed stations measured the temporal variability at single locations. Surveys measured both spatial and temporal variability. Fixed installations included automated radiation stations, the UHI, light chains hanging through the ice into the ocean below the ice, and a light harp, which is a novel in-ice light profile sensor that was colocated with a salt harp system. Optical property survey work included measurements of surface and drone-based albedo and through-ice light transmittance recorded on the ROV. Other measurements that did not use a repeat-survey sampling strategy and also were not tied to a fixed field of view on the ice, included helicopter and HELiPOD albedo, in-ice inherent optical property probe, and ice cores taken for dedicated laboratory assessment of inherent optical property profiles (Table 3).

The surfaces of snow and sea ice, which comprise the dominant properties for RS and large-scale linkages, changed seasonally. While our work during winter was strongly related to changes in snow accumulation, ridging, and surface roughness, surface melt processes were examined in detail over summer. Starting with melt onset, the melt dynamics of the snowpack and the formation and evolution of melt ponds were studied. In addition, the physical properties of melt ponds were investigated, mostly in connection to ecological and biogeochemical properties. Surface properties and melt ponds were linked to process studies on radiative and heat transfer, due to their large impact on the surface energy budget.

During summer and autumn, the flux of snow and ice meltwater into the ocean played a key role with immediate links to the ocean properties, the ecosystem, and biogeochemical fluxes. The formation and evolution of the meltwater layer in leads and under the sea ice was studied in detail in a joint effort across the teams (Section 4.2).

Studies of sea ice dynamics and mechanics were carried out during all legs in the CO. They give insights into the connections of sea ice material properties at the millimeter scale to the visco-plastic behavior of sea ice at the floe scale. Sea ice deformation at the floe scale was measured down to millimeters in the CO using the laser strain instrument. Ridge and lead formation were observed through the many dynamical events in the CO by, for example, airborne mapping by ALS, the EM-Bird, and aerial imagery. Over 120 position buoys in the DN contributed to a regional-scale data set of sea ice drift and displacement. The ice radar system onboard of Polarstern also mapped the dynamics within 5 nm around the vessel (Jäkel et al., 2021).

Ridges and leads had strong logistical implications for the work and fate in/of the ice camp. After ridge formation, the temporal development of the consolidated layer, the macroporosity and the keel depth were studied in a number of ridges using the described mass balance, ice coring, TLS, transect, and ROV methods (Sections 2.3, 2.4, 2.5, 2.10, and 2.11). In situ stresses were measured on the ice floe and will be correlated to atmospheric and oceanic forcing as well as the kinematics of ice drift at different scales.

The mechanical properties of the ice were studied through in situ borehole jack testing with 2 different jacks and through sampling ice cores and testing uniaxial strength in the field. The full annual cycle of data gives a unique data set that can be correlated with small-scale physical properties (temperature, density, salinity) and ice texture allowing for a quantification of small-scale (0.01 m) ice mechanical properties and their seasonal variability.

Measurements in the microwave domain (radiometers, radars, reflected GNSS signals) at different frequencies and polarizations were conducted simultaneously with the extensive physical ice and snow measurement program during the whole MOSAiC drift. The combination of all measured frequencies and polarizations allows a better understanding of the processes, as well as the volume fraction of emitting and scattering constituents and their geometry, which define the seasonal development of microwave signals. Key processes influencing the microwave radiation are, for example, snow metamorphism and ice lenses/layers in the snow, snow liquid water content, wicking of brine into snow and desalination of ice during summer, changes of ice thickness and the snow/ice temperature profiles (Ulaby and Long, 2014). These processes can cause major changes and fluctuations in the microwave signals and can add significant uncertainties in satellite sea ice retrieval algorithms. To develop new satellite retrievals and quantify uncertainties of existing retrievals, first, the processes that influence the microwave retrievals have to be understood better. Second, based on the improved process understanding, better microwave emission and scattering models (Tonboe, 2010; Picard et al., 2018) can be developed. They require a number of poorly constrained input parameters (e.g., correlation length scales or salinity profiles) for which the MOSAiC snow/sea ice physics measurements in combination with the microwave measurements are essential. The same microwave radiative transfer models (sometimes called forward operators) are also needed to assimilate microwave satellite observations in climate and weather prediction models.

The full seasonal cycle was captured by the multifrequency RS observations. Largest uncertainties in satellite retrievals occur when environmental conditions change, like during melt-refreeze cycles, warm air intrusions, or rain on snow events, which all were observed during MOSAiC and are of special importance.

This section presents first results from the MOSAiC drift. These results consist of methodological advances, which make use of technological developments (Section 4.1) as well as scientific results that contrast winter and summer snow and sea ice properties and conditions (Section 4.2). The last part (Section 4.3) discusses remaining challenges of the snow and sea ice observational program.

Figure 5 is called the “sea ice clock” and illustrates the sea ice and surface conditions in the 3 COs based on photographs taken by the panorama camera (Section 2.1). Each image represents the conditions of 1 month during the annual cycle: starting from thin and new ice in September/October, further ice growth with snow accumulation through the dark season until February, into the dense and cold spring ice pack until April, and then snow melt and melt pond formation, and the decay of the ice pack in summer. At the beginning of the drift, CO1 consisted of SYI with a highly compressed core of the floe and refrozen melt ponds over large parts of the floe. Over time new ice formed and the icescape became strongly deformed (see below).

Methodological and technological advances were implemented for the MOSAiC drift in all tasks, allowing a sharp improvement in observing sea ice and its snow cover for an entire year. The measurements benefited from development of sensor and instrument technology over the last decades, advances in data handling and processing, and improved coordination across tasks and disciplines. Technological advances resulted in data sets with higher spatial and temporal resolution as well as with reduced uncertainties. It was possible to obtain long time series with instrumentation that was earlier only used for case studies or dedicated measurements. Covering most of the year, these advanced technologies often turned into a backbone of the MOSAiC data set. The technological (including data processing) advances enabled the immediate use of measurements to plan additional measurements and aid decision-making. Selected examples are as follows:

1. Combining an eye-safe TLS during winter and a non-eye-safe TLS during summer allowed the same regions to be scanned in the CO approximately every 2 weeks (Section 2.5). A custom heated enclosure enabled scans throughout the winter in temperatures below –30°C. During summer, the Riegl VZ6000 allowed for challenging measurements of wet surfaces. A first version of data was processed on board for immediate analysis in support of other measurements (e.g., precisely locating the ROV transponders). Figure 6 shows the results of 2 TLS scans of the same piece of ice, containing level ice and a small SYI ridge, on January 19 and 25, 2020. In this way, the TLS data are used to quantify surface changes related to snow fall and redistribution, as shown in this example. Figure 6D illustrates that the snow deposition is concentrated in a snowdrift off the ridge (left of the peak). Similar results on snow accumulation in ridged areas are expected from the transect data across the COs, while the TLS data provide 2-dimensional accumulation details on the centimeter scale.

2. The micro-CT system (Section 2.2) operation in a freezer-lab container allowed for almost in situ imaging of samples, without the delays and likely changes in microstructure. Immediate scanning was of particular value for snow samples and the more delicate samples of the bottom of sea ice cores, because it ensured that the samples were as close to their original state as possible. Sea ice samples from ice cores were transported in insulated customized boxes that kept, through eutectic phase change cooling elements, temperatures close to the in situ seawater freezing temperature of approximately –2°C. On the ship the samples were centrifuged at this temperature (Weissenberger et al., 1992). This procedure ensured the conservation of in situ microstructure as well as sufficient contrast to retrieve pore network details by a subsequent micro-CT imaging (Section 4.3). The immediate onboard visualization of the 3D structure of snow and ice gave insights into relevant processes and allowed for targeted additional sampling, specifically with respect to the less studied snow–ice interface processes.

3. The snow pit program was carried out by instrument-based measurements, which resulted in a consistent time series with minimal human bias. Main advances resulted from standardizing the snow characterization using micro-CT, SMP, NIR, and structure from motion photography (Section 2.2). These novel quantitative measurement techniques were combined with the concept of defining different types of snow pits (e.g., relatively quick snow pits with fewer measurements versus more elaborate snow pits using the complete suite of measurements; Section 2.2), and thereby ensuring that a standard set of measurements was performed for each snow pit type during all legs.

4. Studies of dynamical and mechanical properties benefitted from the advance in autonomous platforms. In the CO, the laser strain array was operated to investigate sea ice deformation and mechanics on the millimeter scale. These data will give new insights into mechanical properties, in particular in combination with the measurements of stress and strain rates and high precision position monitoring (Section 2.8).

5. Optimizing and winterizing the ROV (Figure S3) and surface unit allowed us to operate the vehicle in all conditions and without interruptions related to sensors or technology during the full annual cycle. Advances in sensor technology allowed for more precise under-ice navigation and synchronous recording of all data streams. Specialized additional sensors, for example, the UHI or suction sampling systems, were integrated into the vehicle in the field. In addition, the reliable operation of the ROV allowed site selection and instrument maintenance under the sea ice. Particularly the sensitive manipulation skills of the ROV system allowed the deployment and recovery of sediment traps, inspection and maintenance of under-ice sensors, as well as sampling without the need for human divers. Figure 7 shows the light transmittance distribution under sea ice in July 2020. Compiling 6 dive days with focus on optical measurements makes it possible to quantify the increasing amount of radiation penetrating the sea ice as melt progresses, allowing for detailed analysis of the increasing spatial variability and the evolution of individual patches with high or low transmittance. The overall evolution shows that the initially unimodal distribution with a modal transmittance of 0.03 develops into a bimodal distribution with modes at 0.15 and 0.26. This shift in modes represents the evolution from wide-spread wet and melting snow and SSL to a white ice and melt pond pattern. For illustration, see the surface images of June to August in Figure 5.

6. The MOSAiC ICE field program strongly benefitted from advances in drone technology. Operation of (small) drones for systematic aerial photography and optical measurements (see below) advanced the study of surface processes. Advances in and miniaturization of both drone and sensors systems allowed for successful operation of these platforms at high latitudes, and accurate measurements of broadband and spectral irradiance. The latter was achieved in part through the stabilization of sensors using miniaturized gimbal systems that could be readily carried by the drones. Additionally, continued advancement of battery systems allowed for extended flight times for these platforms, increasing capabilities with respect to sampling area and achievable altitude. The drones supported the alignment of various measurements at the surface, for example, the optical measurements along the albedo lines (see Figure 8), as well as mapping of the surface at high spatial resolution to provide enhanced information over limited spatial coverage within the CO (Section 4.3).

7. Data from the ALS and the real-time navigation solution of the inertial navigation system of selected flights covering the CO were processed with a short delay of 1 or 2 days after the surveys. The processing steps included the estimation of ellipsoidal elevations of the ALS range measurements, the transformation of the geographical coordinates into a local Cartesian coordinate system using position, and true heading data from Polarstern as well as the merging of all ALS swaths into a single digital elevations model with a spatial resolution of 0.5 m. Visualizations of these models served as valuable science planning information in the absence of large-scale aerial photography during the winter season (Figures 2 and 9). More than 20 ALS surveys provided floe maps throughout MOSAiC, with the largest survey consisting of 3 coordinated flights in late December 2019 and covering approximately 70 km².

8. A newly developed sensor suite (Section 2.12) allowed to observe different important surface properties at the same time. The topography/roughness of the sea ice cover can be directly related to areal properties like ice thickness, melt pond coverage, and surface albedo. With a lateral resolution of 0.5 m and an elevation uncertainty between 2.5 cm (center) and 10 cm (edges), the laser scanner data allow digital elevation models to be derived with a precision never reached before on airborne operations from Polarstern. The TIR camera has a precision of 0.2 K and an accuracy of 1 K and allows detailed mapping of thin ice thickness distributions and lead coverage. The digital single lens reflex and hyperspectral cameras operated during daylight time allow retrieval of albedo and detailed snow and melt pond optical properties. Sea ice thickness surveys (EM-bird) in close temporal vicinity to the camera/ALS flights resulted in synergistic measurements that allow a novel aerial view of sea ice properties. The drones were also used for floe mapping to document and plan the work under daylight conditions, thus supplementing the ALS maps.

9. The deployment of novel and advanced buoy systems in the CO and DN strengthened the acquired time series particularly during the time when Polarstern had to leave the MOSAiC floe. In addition to more classical IMBs and Snow Buoys, emerging technologies such as the continuous measurements of the vertical profile of sea ice salinity (so called salt harp), in-ice and water column light and characterization (optics probe) of under-ice biological activity by the means of fluorometry and autonomous acoustic profiling were key contributions of the ICE team to the interdisciplinary efforts in the DN.

10. Novel on-ice RS instruments were developed for MOSAiC (Ku/Ka-radar, GNSS-R, HUTRAD and L-band microwave radiometers) to evaluate measurements of existing and upcoming satellite missions (Section 2.13). For the first time, a large set of 14 RS instruments, covering both active and passive measurements as well as a wide range of frequencies and polarizations, was operated on the same sea ice floe. Several high-resolution optical and synthetic aperture radar (SAR) satellite sensors performed specific acquisitions for the MOSAiC region. In particular for SAR, a unique data set from newly launched sensors (e.g., Radarsat Constellation Mission, SAOCOM) and wide range of frequencies (ALOS-2, TerraSAR-X, PAZ, COSMO-SkyMed, KOMPSAT-5) were acquired, which will allow novel multifrequency and full polarimetric sea ice analysis at unprecedented temporal resolution. These acquisitions required placing an order of the satellite scenes about 2 days in advance on a daily basis. This process was facilitated by the availability of a dedicated near-real-time drift forecast product for the MOSAiC floe, to have Polarstern as central as possible in the satellite scenes. The product was provided by the Year of Polar Prediction (YOPP; Jung et al., 2016) Sea Ice Drift Forecast Experiment (SIDFEx). For MOSAiC, a consensus ensemble forecast product, based on the different forecast systems, was used. The consensus forecasts were provided onboard Polarstern to support decision-making. Figure S5 exemplifies the consensus forecasts product, showing the forecast issued on February 24, 2020, and revealing that, while the drift over the subsequent 4 months toward Fram Strait was on average considerably faster than anticipated, the trajectory remained within the ensemble uncertainty margins.

Coordinating measurements and sampling strategies within the ICE team (across tasks) and with other teams allowed the required systematic and representative observations of multiple components of the coupled system to be obtained. Here, we present examples of such activities.

Probably the most striking and obvious advantage of MOSAiC was the interdisciplinary coordination across the scientific teams. Although most field measurements benefitted from this advantage, the most prominent example is the coordination of sea ice coring across disciplines. Traditionally physical scientists are interested in the ice–snow or ice–atmosphere interface and reference their work from the ice surface (i.e., the zero-reference level is set at the ice surface). However, most (sea ice) ecological studies are primarily interested in the highly porous underside (i.e., bottom) of the ice, where the majority of the biomass is found (Manes and Gradinger, 2009). Thus, they tend to reference their work from the ice bottom upward. This difference often makes direct comparisons of data collected with the 2 approaches difficult. During MOSAiC, the sectioning of ice cores for physical and biological sampling was optimized such that sectioning was initiated both from the top and bottom, and the middlemost part of an ice core was left with the “odd” length section. Colocation of coring with the ECO and BGC teams captures the seasonal evolution of both physical and biological properties for FYI and SYI. Approximately 1,500 ice cores were collected during approximately 30 visits at FYI sites and 25 visits at SYI sites. This coordination effort will allow better direct comparison of ice core data across teams.

Snow sampling was coordinated with other teams, usually collected by the team making snow pit measurements, independently of the purpose and later analysis of the snow samples. This strategy ensured that samples from each sampling event were taken at the same site(s) and same time and in a consistent and reproducible manner throughout the campaign. Special emphasis was given to sampling colocated with other measurements, for example, at the RS site. For interpretation of the RS measurements, the snow and ice physics measurements are critically needed.

The observations of low-salinity meltwater layers under the sea ice and in leads during summer resulted in a dedicated interdisciplinary program across all teams. Drill lines spanning in total 900 m were completed covering different ice types with the specific aim to measure the ice thickness and ocean temperature and salinity profiles directly beneath the ice and to map the presence of “false bottoms” formed at the meltwater–seawater interface. These 130 ice thickness observations, of which 30 have associated ocean temperature and salinity profiles, overlapped with existing long-term mass balance observations (stakes and thickness surveys). Additional melt pond, ice core and under-ice, and lead water samples were collected by the ICE, ECO, and BGC teams to better understand the meltwater sources, and the role of these meltwater layers in physical, chemical, and biological processes in summer.

MOSAiC offered the opportunity to colocate measurements of sea ice surface and bottom topography, the internal structure of pressure ridges and the biology associated with ridges. This approach allows a complete 3-dimensional view of ridges. Figure 9 shows the surface and bottom topography as well as a visual photograph of the pressure ridge called “Jaridge.” This ridge was located between the SYI closer to Polarstern and the FYI where the ROV and the FYI coring site were located (Figure S3). Jaridge was formed during deformation events primarily during winter (Leg 2) and included parts of the floe adjacent to the former dark site coring site. The Jaridge ridge consisted mostly of 20–40 cm thick ice blocks. This site was used for a large suite of coordinated measurements and installations: an IMB was deployed, several transects with drilling for ice thickness and consolidated layer thickness were undertaken, sea ice cores were collected (ICE and ECO), sediment traps were deployed with the ROV, ROV multibeam and the UHI was used to map ridged and level ice, biological samples were taken in the ridge, and airborne mapping was conducted with ALS and cameras. This major effort will provide a comprehensive and unique view on the evolution of sea ice ridges in the Arctic summer.

The evolution of the sea ice surface during the summer is a culmination of changes in the radiative energy and mass budget due to variation in snow cover, melt ponds, and bare ice. Capturing the heterogeneity in surface conditions, in both space and time, is key to understanding the seasonal evolution of surface albedo and thus surface energy balance. Surface measurements of snow and melt pond properties, and albedo (both spectral and broadband), were coordinated with aerial measurements of both surface topography and albedo. These observations also overlapped with TLS scans (Section 2.5). Micro-CT observations of the ice and snow provide novel insights into the evolution of the surface microstructure in relation to the observed radiation budget. In addition to the surface-based observations of ice and snow, surface albedo measurements and imagery of the surface were collected using drones. These activities included flights over the surface albedo lines to compare observations from surface pyranometers and spectral radiometers with those on the drones. To align measured albedos along the measurement survey line, drone GPS coordinates were corrected to account for ice drift, and the albedo was averaged every 2 m from 4 individual flight passes along the line. Mean drone-based albedo values and standard deviation measured along the “lemon drop” albedo line on July 06 are shown in Figure 8 (red line). The drone operating along this albedo line carried a modified Kipp and Zonen PR1-V broadband pyranometer, measuring wavelengths between 310 and 2700 nm. Surface and drone-based measurements show good agreement, with differences primarily resulting from the difference in sampling altitude (1 m for surface-based measurements and 15 m for drone-based observations), with the higher sensor aggregating the influence of surface type variability. This influence of sensor altitude on observed surface albedo is being investigated in a separate study.

Figure 10 shows winter and summer time-series measurements of several colocated RS instruments (Table 2 and Figure S4). During a storm event in November 2019, measurements by both radars and radiometers show changes in backscatter and brightness temperature, respectively. Due to the fact that air and surface temperatures reach a maximum of –5°C, melting plays no role in these changes, but snow metamorphism already can happen. Also, strong winds of up to 20 m s^–1 caused snow compaction and redistribution. The radars and radiometers show stronger changes for higher frequencies and partly opposite behavior depending on frequency. The 35 GHz radar backscatter decreases by close to 2 dB while at 15 GHz it increases by about 1 dB during the event. At the lowest frequency of 1.3 GHz, no changed radar response is observed during the event, which is in correspondence with the minimal changes in GNSS-R reflectivity operating at a similar frequency. During the storm event, the increase in magnitude of microwave brightness temperatures (TB) is increasing with frequency: for example, at 89 GHz TB is increasing by about 50 K, which corresponds to a strong emissivity change, while at 1.4 GHz almost no TB increase is observed. The peak of the storm on November 16, 2019, could not be captured by most of the RS instruments because a crack opening, due to strong ice dynamics, at the site caused a power outage.

In summer, on September 13, 2020, even stronger fluctuations in both backscatter and brightness temperatures occurred when the air temperature became positive and melting increased. Rain on snow additionally changed surface properties. GNSS-R reflectivity, both at the RS site and on Polarstern, increased during the event. On September 15, after the temperature dropped below zero again, measurements also stabilized. Similar to the winter event, the changes in magnitude and variability increased for higher frequencies both for the radars as well as the radiometers, with the exception of the 1.4 GHz radiometer. However, changes were much stronger in summer compared to winter (note the different y-scales in Figure 10).

Ice thickness influences the signal at very low microwave frequencies. Measurements by the 1.4 GHz L-band show lower brightness temperatures at a partly refrozen lead (light blue line) than for the thicker ice at the RS site (dark blue line; Figure 10). This effect is used to retrieve thin ice thickness from L-band satellite measurements from SMOS and SMAP (Kaleschke et al., 2016; Paţilea et al., 2019) and in future CIMR. The higher measurement variability at the lead likely is due to the changing ice conditions caused by ice dynamics.

Variability in brightness temperature and backscatter induced by environmental changes like those shown in these 2 case studies are not accounted for in current satellite sea ice retrievals for ice concentration, thickness, or snow depth and are causing uncertainties in the satellite data sets. In conclusion, during both events, lower frequencies were less affected by the environmental changes and thus should result in more stable satellite retrievals. By combining the on-ice RS measurements with the detailed snow and ice physics measurements, we can better understand and model the interactions of microwaves with snow and ice. Future work will improve satellite RS methods to obtain better sea ice satellite climate records.

To demonstrate the comprehensiveness of the snow and sea ice measurements during MOSAiC, 2 weeks of measurements are selected as case studies. The week starting January 20, 2020, (Leg 2) represents the work during winter darkness in the polar night and the week starting on July 06, 2020, (Leg 4) represents the work during summer. Figure S1 shows how the work from the different tasks was distributed over the respective week. Additional optical measurements were carried out under daylight conditions.

Sea ice, snow, and weather conditions differed significantly for both case study weeks. As expected for high Arctic conditions, the mean air temperature was –29.0°C during the week in January and 0.3°C during the week in July. Wind speed at Met City (Figure 2) was 4.6 m s^–1 and 6.4 m s^–1 for the weeks, respectively. In general, the winter period was characterized by a transition from weak winds with periodic thin, ice-dominated clouds to a couple of passing snow storms later in the week with increased winds. The summer week was similar in some regards, but with much higher temperatures, characterized by persistent low-level, liquid-dominated cloud and fog, with 2 short storms that produced light drizzle near the surface. A detailed description of the meteorological conditions at the CO is in preparation led by M. Shupe.

To contrast the different snow and sea ice conditions, Figure 11 illustrates various aspects of the surface and ice conditions. The figure shows the CO and its surrounding on composites from hundreds of aerial photographs taken from the Polarstern helicopters. During polar night, TIR photos were used to map the surface temperature, which indicates sea ice types based on differences in surface temperature. The map shows a dense ice pack with few deformation zones and elongated structures of thinner ice (higher temperatures) resulting from leads with new thin sea ice. During summer, visible photos were taken and give a view of the surface conditions. Toward the end of June and in July, the sea ice surface consisted of many distinct floes, widely covered with melt ponds. The observation over the year of the same ice pack and area will allow consideration of how much of the wintertime surface and structural features determine summer features, such as the melt pond distribution. In general, the observational concept of TIR and ALS flights during the dark season proved most useful (Figure 2). Both methods were used to map surface conditions and guide the on-ice measurement program as well as to support logistics through results on snow and sea ice surface properties, immediately processed on board Polarstern.

The drift and deformation of the ice pack was one of the most remarkable features of the MOSAiC field experiment, and more extensive than expected. In total, the ice camps covered distances of 2,354 km (Drift1 with CO1), 622 km (Drift2 with CO2), and 345 km (Drift3 with CO3) based on hourly position data during the manned observations. These tracks are shown in Figure 1, which also shows the continued drift of the COs after Polarstern had departed (Table 1). Satellite data indicate that the MOSAiC drift was about 25% faster than the climatological mean drift, and was caused by large-scale low-pressure anomalies that prevailed around the Barents-Kara-Laptev Sea region between January and March (Krumpen et al., 2021). Ice drift and deformation was tracked in the DN as well as in the CO. Figure 12 compares the ice drift speed with the wind speed for January and July and discriminates whether the floe was in free drift or moved with the closed ice pack. It reveals much higher drift speeds with a mean of 0.19 m s^–1 in July, when the floe was in free drift, while drift speeds with a mean of 0.10 m s^–1 were only about half that value in January, when the floe was embedded in the interior ice pack with other floes. In addition, the variability of drift speed was 1.5 times higher during free drift in July than in January. Additional analyses of this data set for the entire year will reveal details of the seasonality of different drift modes as well as dominating processes, as the start of inertial motions became more important around mid-June (Dethloff et al., 2021). Seasonal contrasts in drift and deformation patterns can also be seen in the ice radar animations over the entire CO1 (see reference to video animation in data section). Future work will analyze the interaction between ice fractures and wind forcing also on larger scales by combining the wind stress information across the DN.

Figure 13 shows the rate of deformation (⁠$divergence2+shear2$⁠) derived from 2 consecutive Sentinel-1 SAR images in a 200 × 200 km distance centered around the MOSAiC CO for winter (December 2019) and summer (June 2020) conditions (for method see Albedyll et al., 2021). Intersecting lines of strong deformation (linear kinematic features) are present both in winter and summer ice pack. Differences in winter and summer deformation, for example, in the degree of localization of the deformation visualized in the width of the linear kinematic features (Figure 13), are caused by changes in the ice pack strength. The seasonally varying deformation provides additional insights to the different sea ice conditions previously presented in Figure 11, and adds the larger scale to the airborne observational data. Additional analyses of spaceborne RS are described in Krumpen et al. (2021). The SAR-based ice deformation can, in future, be combined with the buoy observations from the DN and the ship-radar-based deformation to obtain data products of high spatial and temporal resolution. Thus, these observations may be used toward an improved understanding of ice dynamics on different spatial scales.

On the floe scale, Figure 14 describes the distribution of sea ice and snow thickness for January and July, and melt pond depth (for July only). The transect data (Section 2.10) shows a winter ice mode of 1.3 m and a thicker summer mode of 1.8 m (mean 2.7 m). The winter and summer modes of snow and sea ice thickness distributions demonstrate the effects of thermodynamic ice growth/decay and deformation throughout the seasonal cycle. The winter distribution features 2 distinct total (snow plus ice) thickness modes associated with a thinner young ice and a thicker mix of level/deformed FYI and SYI. In summer, in turn, the distribution exhibits 2 modes of level FYI and SYI accompanied by a broad tail of thicker ice formed by deformation processes during the preceding winter. The transect lines between winter and summer had only partial overlap. As one would expect, snow thickness was much reduced in summer compared to winter/spring, which is obvious from the reduction in modal snow thickness from 0.06 to 0.03 m. In July, the modal melt pond depth was 0.1 m, with maximum depth of 0.5 m. The large coverage of melt ponds in summer is shown in Figure 11B. Results from a helicopter electromagnetic ice thickness survey on July 01, 2020, (data not shown) revealed a thicker summer modal total thickness of 2.1 m, but a thinner mean total sea ice thickness of 2.5 m over the CO and the surrounding ice within a radius of 50 km.

The characteristics and evolution of FYI is exemplified in Figure 15 and data are given in Tables S2 and S3. In winter (December), the upper 15 cm from the ice surface consisted of granular ice (Figure 15A), while most of the core (19–65 cm) was composed of columnar ice. A 6-cm thick skeletal layer was apparent at the core bottom (65–71 cm) of the growing FYI. While the ice in December (Figure 15B) had the typical C-shape salinity profile and was cold due to the low atmospheric temperatures, the rapid warming and desalination of FYI in summer is shown with the data from early July (Figure 15C).

The surface properties of the topmost layer of sea ice showed large temporal and spatial variations throughout the MOSAiC observation period. The typical winter snowpack consisted of an often wind-packed new snow layer at the surface, with low density, large specific surface area and relatively small optical equivalent diameter (0.00–0.15 mm; Figure 16A). Due to the large temperature difference between the air and the ocean, the snowpack experienced large temperature gradients during most of the winter, which caused strong recrystallization and faceting into columnar snow structures with a strong geometric anisotropy (depth hoar). We observed the formation of euhedral crystals in the low-density regions (0.40–0.65 mm; Figure 16A), in contrast to subhedral crystals in the layers with higher densities (0.25–0.40 mm; Figure 16A). Our data show that the strong temperature gradients in the Arctic snowpack during winter cause a highly anisotropic microstructure that has significant effects on the thermal conductivity of snow and sea ice growth. In addition, they affect microwave properties for RS (Figure 10). The sea ice surface in summer (Figure 16B) is generally free from snow and consists of the SSL that is formed by melting and draining processes in the sea ice. In the field visually distinguishing between snow and the SSL is hard; however, data from the micro-CT highlight the microstructural differences between the SSL and snow (Figure 16B). The density and the optical equivalent diameter increased with depth, reflecting the transition from the drained and melted sea ice at the surface to values typical of frozen sea ice (approximately 700 kg m^–3) toward the bottom of our sample. The specific surface area was relatively consistent, showing an increase at the surface due to preferential melting of columnar ice. The strong anisotropy of the SSL reflects the preexisting structure of the original columnar sea ice crystals and brine channels and influences optical properties of the sea ice surface as well as RS retrievals.

The observational program of MOSAiC had to compromise in various aspects. The largest impact on the observations of the full annual cycle was that Polarstern had to leave the main floe (CO1) temporarily to exchange personnel and to supply the vessel. This absence resulted in a gap in manned observations between May 16 and June 19, 2020 (Table 1), a key phase of the annual cycle: the spring–summer transition and the early melt onset. The formation and decay of sea ice during the shoulder seasons is of increasing importance in the increasingly seasonal ice cover of the Arctic. Studying these processes will need dedicated expeditions to observe sea ice and snow processes of thin and very fragile ice covers in very close connection to oceanographic and atmospheric conditions. During MOSAiC, only autonomous measurements in the DN were able to bridge this gap to some extent and to obtain some key parameters continuously across this gap. Nevertheless, further studies of similar interdisciplinary complexity are needed to add more coordinated in situ observations of this phase, where small changes in timing have large impacts on the total energy budget and the seasonality of ecological processes (Nicolaus et al., 2012).

All observations are limited to one specific drift trajectory along the Transpolar Drift from 2019 to 2020. This leaves the challenge of upscaling and generalization, which will mostly rely on merging the in situ observations with numerical models (Section 5.2) and RS (Section 5.3). Although many relevant sea ice processes were very well covered, individual parts are likely missing or lacking details. The work on the snow and sea ice properties in the DN was, in the end, quite limited both in space and time. The need for a permanent ice camp and fully manned research icebreaker for most aspects of the field program excluded most observations of marginal ice zone processes.

MOSAiC demonstrated the benefit of improved coordination across disciplines, which allowed us to study linkages and to realize ambitious plans in integrating methods and disciplines. Section 5.1 describes linkages to the individual subsystems (including the work of the other 4 teams) and points to the comprehensive data sets, which will build the legacy of MOSAiC. The sampling strategy was guided by the needs of sea ice and climate models as well as for satellite RS. These linkages are discussed in Sections 5.2 and 5.3, also in connection to the DN, because they allow upscaling of observations relevant for climate models and RS.

Combined analyses of key parameters across all subsystems will allow unprecedented studies of the coupled Arctic system (Rabe et al., 2022; Shupe et al., 2022). The physical properties of snow and sea ice are crucial for interpreting biogeochemical and biological observations as well as for understanding atmospheric composition and aerosols (Figure 4). For example, the evolution of snow and sea ice microstructure (Figures 14 and 15), porosity and permeability, is known to affect the transfer of trace gases (e.g., Nomura et al., 2018) and aerosol particles across the atmosphere–ice–ocean boundary (Frey et al., 2020). The snow on sea ice and the sea ice both act as a chemical reservoir and reactor releasing a range of chemical trace gases, aerosol particles, and their precursors with significant impacts on tropospheric composition and oxidizing capacity, and potentially on climate via influencing clouds (Grannas et al., 2007). Observations of gas and particle fluxes across the atmosphere–snow–ice–ocean boundary were carried out at different scales by the BGC team and ATMOS team. Here, both snow and sea ice properties from different ice types as well as ice dynamics (especially formation of leads, e.g., Figure 13) are likely crucial factors for the observed fluxes. The close coupling of the snow and sea ice cover with atmospheric forcing (especially wind forcing) became obvious (Graham et al., 2017). The dynamics of the ice pack challenged the field program, but dynamics are critical for a better understanding of energy and mass exchanges between atmosphere, sea ice, and ocean. Here, we gave examples of dynamical processes with respect to the fast drift (Figure 12), the formation of leads and ridges (Figures 9 and 13), and their effects on the sea ice thickness distribution (Figure 14). While the sampling activity (coring, snow sampling) captures the seasonality, more rapid or transient changes (e.g., due to warm air advection) are more reliably covered by continuous measurements of autonomous systems; their combination provides new insights into the coupling between snow and sea ice physical properties and atmosphere–ocean exchanges.

The role of melt ponds and open water fraction was investigated intensively with respect to formation processes (Figure 5), energy budgets (Figures 7 and 8), and the impacts of their large-scale distribution (Figure 11) on the sea ice mass balance (Figure 14). Findings from these investigations will also link to concurrent ecological and biogeochemical studies. The small-scale surface heterogeneity and, in particular, the contrasts of surface temperature, trigger exchange processes between the Arctic system components. During low ice surface temperatures, the warmer open water areas or partly refrozen leads act as a source of latent and sensible heat, which modifies the atmospheric boundary layer and triggers atmospheric convection (e.g., Lüpkes et al., 2008; Schmale et al., 2021). The strong release of water vapor into the atmosphere in summer results in enhanced formation of low-level clouds and fog (Tjernström et al., 2012). While the leads represent a conduit for strong heat loss (Figure 11) and new ice formation during winter (Itkin et al., 2018), they are windows in the sunlit summer (Figure 7) period that allow more solar heating of the upper ocean (Taskjelle et al., 2017) and intensified pelagic primary production before melt ponds appear (Assmy et al., 2017).

The observations of sea ice mass balance provide an integrated measure of the energy balance, which is a result of the interactions and energy exchange with the atmosphere and the ocean. Here, both thermodynamic and dynamic forcing are relevant. The first observations of platelet ice in the central Arctic during winter (Katlein et al., 2020) suggest that this mechanism can also contribute to the mass balance of Arctic sea ice. The detailed observations of ocean heat and freshwater content and fluxes by the OCEAN team provide data sets to examine the interaction between upper ocean properties and sea ice under different ice regimes.

The detailed work on ridges and leads (Figure 9), in close coordination and collaboration with oceanographic and ecological aspects, represents a leap forward from simply investigating the level and modal snow and sea ice properties toward a much more comprehensive understanding of the whole ice pack. These measurements allow specific analyses of deformed and ridged ice (Gradinger et al., 2010; Fernández-Méndez et al., 2018; Shestov et al., 2018), for example, through comparisons of sea ice habitats, energy and mass fluxes, and contrasting processes of level ice with deformed ice, thin ice, and leads.

A synthesis of sea ice field measurements with results from numerical models involved in MOSAiC is expected to demonstrate how process models with different levels of complexity and coupling can help to advance our understanding and prediction of Arctic climate change and contribute to improved numerical weather and sea ice prediction. Some of the expected progress is related to:

• the assessment of unresolved, missing or overly simplified or underrepresented processes in the sea ice component of the fully coupled Arctic system;

• the exploration of inadequacies and possible improvements along the model atmosphere-sea ice and ocean-sea ice coupling channels;

• the analysis of model sensitivities to the scale-aware sea ice parameter space; and

• the placement of the MOSAiC sea ice measurements within a broader spatial (i.e., pan-Arctic) and temporal (i.e., interannual to multidecadal) context.

Major challenges remain in how to optimally utilize observational data in different models, including those related to:

• using Lagrangian point-observations for evaluation of gridded, space- and time-discrete model output;

• postprocessing heterogeneous, erroneous and inconsistent data sets into data products that may be useful for model initialization, forcing and evaluation;

• bridging different scales of in situ observations (e.g., from the CO) and remotely sensed measurements in order to up- or down-scale such data to a model grid cell or mesh; and

• homogenizing the heterogenous data sets of individual parameters into a common data product, for example, by interpolating to common time stamps or by classifying data according to the snow and sea ice states.

On the model side, future work will be needed to:

• evaluate sea ice models and their existing parameterizations in a process-oriented manner that benefits from an improved process understanding;

• use process models, large eddy simulations or direct numerical simulations to derive, test, and make available advanced parameterizations of processes impacting snow and sea ice states and variability;

• assimilate sea ice data to produce a better state estimation and forecast initial conditions for near-term forecasts; and

• derive improved sea ice forecasts and projections across timescales, including the consideration of internal and forced variability.

During the field phase of MOSAiC, different observations were realized to enable a process-oriented model evaluation and thus to assess and improve systematic model biases through improved process understanding and related parameterizations (e.g., Figure S5).

From the sea ice perspective, the most immediate connections, and thus model developments are currently expected from:

1. Improvements to the representation of snow on sea ice (Figures 6, 14, and 16). The evolution of the snow cover (fraction, depth, distribution, thermal properties) is often overly simplistic in various types of models, for example, in global coupled climate models (Chen et al., 2021; Webster et al., 2021) and models to reproduce reanalyses (Sato and Inoue, 2018; Batrak and Müller, 2019; Wang et al., 2019). Model sensitivity studies also indicate that the simulated climate is very sensitive to the snow representation (Urrego-Blanco et al., 2016). The new field observations cover all relevant processes (Figure 4) from freeze-up to melt, including RS observations on large scales. Particular improvements are expected for estimates of snow accumulation and redistribution over deformed ice (Liston et al., 2018), surface topography (Figures 6 and 9), thermal properties, and heat transfer between atmosphere and ocean, as well as internal snow processes (Wever et al., 2020).

2. Improvements in surface albedo and optical properties (Figures 7 and 8). The evolution of optical properties, along with other energy transfer terms, has huge implications for the total energy budget and is still one of the most critical aspects of sea ice modeling and tuning (e.g., Holland et al., 2012; Jäkel et al., 2019). Model studies indicate that the treatment of the albedo, including factors like melt ponds, modify the transient climate evolution (e.g., Holland et al., 2012). They also show the need for inclusion of specific surface albedo for different surface types and the factors that drive the evolution of those surface types over time. Based on the MOSAiC observational data sets, better implementation of albedo schemes and melt pond thermodynamics is expected, leading to improved heat fluxes over the heterogeneous pack ice and more realistic surface energy budgets. These improvements will also have important implications for atmospheric processes and how they are represented in models, due to the strong coupling between surface fluxes and boundary layer processes.

3. Improvements in sea ice dynamics (Figures 12 and 13). An improved inclusion of dynamical processes, including deformation and sea ice surface/bottom roughness (Figure 9) is critical for an accurate description of momentum fluxes across the atmosphere–ice–ocean interfaces and internal stress (Hutchings et al., 2011). Here, the formation and decay of ice pressure ridges and the sea ice floe size distribution are of particular importance to describe the atmosphere-ice and ice-ocean drag and internal ice stress. Both thermodynamical and dynamical processes are important for sea ice change (Dethloff et al., 2021). The main benefit of the MOSAiC observations is expected from the dedicated program on pressure ridges and leads in the CO (local scale) and the unique DN, which allows deriving dynamic parameters on a broad range of scales (e.g., ice radar, helicopter surveys, satellite RS).

To enable model evaluation and development of improved parameterizations, a main goal is to merge the existing observations into an explicit “MOSAiC standard forcing and benchmarking data set” for thermodynamic/dynamic sea ice models. A second goal is to advance coupling with ecological and biogeochemical processes, where the colocated measurements will be most important. They bridge across geophysical and biological aspects of snow and sea ice and allow for benchmarking simulations of causal relationships between environmental conditions and responses of different components of the Arctic sea ice ecosystem.

Over the last decades, assimilation of sea ice and snow data became an important tool to study sea ice and its role in the coupled climate system and to initialize forecasts. The field data will be used to improve the assimilation of sea ice concentration, sea ice thickness and snow depth in coupled ocean–sea ice models (Fritzner et al., 2019), in particular by calibrating and validating satellite-derived products and thus reducing uncertainties in the large-scale assimilation data sets (Figure 10). The impact of assimilating subgrid scale sea ice thickness distribution derived from Cryosat-2 and ICESat-2 observations will be assessed and support simulations with a stand-alone sea ice model.

The field phase of MOSAiC also demonstrated how model applications can support and guide the field experiment. Sea ice simulations and in particular ice drift forecasts were provided to Polarstern in near real time through the SIDFEx initiative (Figure S5). Supplementing the atmospheric forecasts, these provisions supported the on-ice measurement program, for example, in planning intensified observation periods. These data sets were mainly used for near-real-time verification of the model performance and for real time educational applications.

As part of the MOSAiC expedition, several airborne surveys were conducted with the AWI research aircraft Polar-5 and Polar-6 from Longyearbyen (Spitsbergen) between August 17 and September 17, 2020. The aircrafts are of type Basler BT-67 and were operated by the Canadian Company Kenn Borek Air Ltd. Calgary. The main goal of these survey flights was to extend the spatially limited observations of the atmosphere and sea ice in the vicinity of CO to a larger area. While the airborne RS from the Polarstern helicopter cover scales up to 100 km with a high, weekly temporal sampling, the aircrafts cover several hundreds of kilometers and can carry more instrumentation. When the airborne surveys began in August 2020, Polarstern was already on its way to the central Arctic to relocate to CO3 (Figure 1), and thus out of aircraft range. However, various DN buoys, positioned between 79°N and 81°N in early September, were still active and within range. These buoys allowed the sea ice conditions over the DN area to be recorded after Polarstern had departed. The surveys also complement earlier sea ice surveys made in Fram Strait between 2001 and 2018 (Belter et al., 2021). For this purpose, the aircraft Polar 6 was equipped with a number of sensors, including the EM-Bird, a laser scanner, and optical instruments. Together with the laser scanner, the optical camera recorded melt pond distribution, surface elevation, floe size distribution, and other surface properties. An insight into the different sensor systems and their specifications is given in Herber et al. (2021). The track of a flight over the DN made on September 02, 2020, is shown in Figure 1. Average ice thickness in the vicinity of the DN was 1.44 m, while the modal ice thickness was 0.93 m. Note that in parallel to the MOSAiC airborne campaign, sea ice and oceanographic surveys were carried out by the Norwegian research vessel Kronprins Haakon in Fram Strait, and a sea ice floe next to the vessel was overflown while in situ observations were taken. The second aircraft Polar 5 was equipped with in situ, RS, and basic meteorological instrumentation. The observations focused on characterizing the Arctic atmospheric boundary layer in conjunction with ocean-atmosphere interactions, clouds, solar and terrestrial radiation, and aerosols. More details about the airborne atmospheric survey program can be found in Herber et al. (2021).

Satellite observations can extend the local MOSAiC observations to regional and hemispheric scales. Long-term satellite time series and climate data records can put the MOSAiC observations into a temporal context and answer questions about how representative or unusual was the MOSAiC year (Dethloff et al., 2021; Krumpen et al., 2021). Routinely taken satellite acquisitions already provide a large collection of daily sea ice related properties from space, which will help and extend the MOSAiC data analysis: sea ice area, thickness, ice type, drift, albedo, snow depth, melt pond coverage and more. However, in addition, a large collection of satellite data was acquired specifically for MOSAiC, especially high-resolution SAR and optical data. During sunlight, observations of various optical sensors with different footprints are available. WorldView-2 (meter footprint), Sentinel-2 (tens of meter footprint), Sentinel-3, and MODIS (hundreds of meter footprints) allow the scaling of in situ to satellite observations at these scales. Several space agencies helped to acquire a unique data set of radar backscatter images at different spatial resolution, frequencies, and polarizations from several different SAR sensors (ALOS-2, COSMO-SkyMed, KOMPSAT-5, PAZ, Radarsat Constellation Mission, SAOCOM, Sentinel-1, TerraSAR-X). This help allowed monitoring the CO and DN and an unprecedented temporal resolution of several acquisitions per day. However, MOSAiC was drifting through the “pole hole” of many satellite sensors (Figure 1) for a significant amount of time and thus not all observations are available year-around. Retrieved quantities from these higher resolution satellites acquisitions will be, for example, floe size distributions, lead and ridge locations, ice types, ice deformation, melt pond, and albedo distributions. On the other hand, the MOSAiC in situ observations will help to improve and develop new satellite retrievals (descriptions in Sections 2 and 3 and examples in Section 4).

The snow and sea ice field program was designed to monitor all key parameters of the snow and ice system over one full annual cycle, while the intensity of observations over the year were adapted with respect to (1) deciphering sea ice evolution from autumn freeze-up to summer melt, (2) capturing key events and seasonal changes, and (3) accommodating targeted research on emerging science questions. The resulting data set will allow to better quantify the causes and consequences of the evolving and diminishing Arctic sea ice cover, as one of the central elements of the Arctic climate system. The first results show the active dynamics of the entire Arctic ice pack. MOSAiC experienced and observed rapid ice transformations and motions, and strong deformations along the drift in all seasons. These observations demonstrate the importance of dynamic processes in comparison to thermodynamic processes, as well as the role of ridges (deformed sea ice) and leads in comparison to level sea ice. Leads, in particular, link the ice pack to atmosphere-ocean exchange. The program also captured the evolution of the snow cover, across different ice types, including leads and deformed (ridged) ice. Strong temperature gradients during winter cause a highly anisotropic microstructure with significant effects on the thermal conductivity of snow. The large spatial and temporal variabilities of the snow pack impact sea ice growth and have to be considered in more detail, both in observations and in simulations. These details will allow better understanding of feedback processes between the ice and the atmosphere, including gas and particle exchange. The snow and ice measurements combined with the RS observations will lead to a better understanding of sea ice microwave emission and scattering, including their temporal variability which is needed to improve satellite data sets. In particular, the role of events like the passage of low-pressure systems (storms), warm air intrusions, rain and snow fall, lead and ridge formations can be studied using integrated data sets from all key parameters. Ice–ocean interface processes, such as the formation of platelet ice, the summer freshwater layer or the roughness of the ice–ocean interface, were observed in great detail and will likely guide upcoming research with respect to the changing Arctic sea ice cover.

Even though the MOSAiC snow and sea ice work was extremely comprehensive, some aspects could only be studied in a very limited way. For example, MOSAiC was only marginally able to study the transitions of sea ice from/into open water. Dedicated studies in the marginal ice zone may be able to shed more light on this topic in the future. Furthermore, the temporary departure of Polarstern from the ice due to logistical reasons resulted in a gap in many of the in situ observations during melt onset and transition into summer. Despite these limitations, we expect that the overall work will lead to a better process understanding of snow and sea ice and their linkages to atmosphere and ocean as well as improved forecast capabilities. Our expectation is that this work will strengthen research on the global coupled climate system, in particular with respect to the seasonally frozen Arctic Ocean.

All data in this manuscript are publicly available from online repositories. Note that most data sets contain raw or preliminary data, while advanced versions will become available in future. The data may be found under the following references: drift track data (Figure 1, Nicolaus et al., doi:10.1594/PANGAEA.937204), observational dates (Figure 4, Nicolaus et al., doi:10.5281/zenodo.5898517), panorama photographs (Figure 5, Nicolaus et al., doi:10.1594/PANGAEA.938534), TLS data (Figure 6, Clemens-Sewall et al., doi:10.18739/A27S7HT3B), ROV radiation data (Figure 7, Nicolaus et al., doi:10.1594/PANGAEA.935688), surface albedo data on ground (Figure 8, Smith et al., broadband data under doi:10.18739/A2KK94D36 and spectral data under doi:10.18739/A2FT8DK8Z) and from the HELiX drone (Figure 8, Calmer et al., doi:10.18739/A2GH9BB0Q), on-ice RS data (Figure 10, Spreen et al., doi:10.5281/zenodo.5725870), surface images from thermal infrared and true color (Figure 11, Thielke et al, doi:10.1594/PANGAEA.934666), drift speed data from Polarstern (Figure 12, Nicolaus et al., doi:10.1594/PANGAEA.937204), deformation data from SAR (Figure 13, von Albedyll et al, doi:10.5281/zenodo.5195366), sea ice thickness and snow depth distribution (Figure 14, Hendricks et al., doi:10.5281/zenodo.5155244), sea ice physical properties (Figure 15, in Tables S2 and S3) with a sea ice core overview (Granskog et al., doi:10.5281/zenodo.4719905), snow pack properties (Figure 16, Macfarlane et al., doi: 10.1594/PANGAEA.935934), and ship radar video sequence (Jäkel et al., doi:10.5446/52953).

The supplemental files for this article can be found as follows:

Text S1 Sea ice and snow work and concept

Text S2 Methods and field setup details

    Sea ice coring

    Remotely Operated Vehicle (ROV)

    On-ice remote sensing

    Sea ice drift forecast

    Technological challenges

    References

This work was carried out, and data used in this manuscript were produced as part of the international Multidisciplinary drifting Observatory for the Study of the Arctic Climate (MOSAiC) with the tag MOSAiC20192020. We thank all persons involved in the expedition of the research vessel Polarstern during MOSAiC in 2019–2020 (AWI_PS122_00) as listed in Nixdorf et al. (2021). We highly acknowledge the contribution by Tânia Casal, Claude Duguay, Hajo Eicken, Joel Johnson, Ron Kwok, Michael Lehning, Robert Rember, Mike Schwank, and John Yackel as principal investigators on different projects that contributed to financing and enabling this work. We are grateful to Anja Nicolaus, Annekathrin Jäkel, Dirk Notz, Frederike Fohlmeister, Jeremy Wilkinson, Lovro Valic, and Pedro de la Torre for supporting the field work by the provision of instrumentation, data handling, and expertise in setting up the snow and sea ice observations. We thank Marika Holland and Wieslaw Maslowski for providing comments and fruitful discussions in preparing this manuscript. The comments of Stephen Ackley and Jody Deming as editors as well as Pat Langhorne and an anonymous reviewer are highly appreciated, as they helped improve this long and complex manuscript significantly.

This work was funded by the following:

• – the German Federal Ministry of Education and Research (BMBF) through financing the Alfred-Wegener-Institut Helmholtz-Zentrum für Polar- und Meeresforschung (AWI) and the Polarstern expedition PS122 under the grant N-2014-H-060_Dethloff,

• – the AWI through its projects: AWI_ROV, AWI_ICE, AWI_SNOW, AWI_ECO. The AWI buoy program and ROV work were funded by the Helmholtz strategic investment Frontiers in Arctic Marine Monitoring (FRAM),

• – the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) through the Transregional Collaborative Research Centre TRR-172 “ArctiC Amplification: Climate Relevant Atmospheric and SurfaCe Processes, and Feedback Mechanisms (AC)³” (grant 268020496), the International Research Training Group 1904 ArcTrain (grant 221211316), the MOSAiCmicrowaveRS project (grant 420499875), the HELiPOD grant (LA 2907/11-1), and the SCASI (NI 1096/5-1 and KA 2694/7-1) and SnowCast (AR1236/1) projects,

• – the BMBF through the projects Diatom-ARCTIC (03F0810A), IceSense (BMBF 03F0866A and 03F0866B), MOSAiC3-IceScan (BMBF 03F0916A), NiceLABpro (BMBF 03F0867A), SSIP (01LN1701A), and SIDFExplore (03F0868A),

• – the German Federal Ministry for Economic Affairs and Energy through the project ArcticSense (BMWi 50EE1917A),

• – the US National Science Foundation (NSF) through the project PROMIS (OPP-1724467, OPP-1724540, and OPP-1724748), the buoy work (OPP-1723400), the MiSNOW (OPP-1820927), the snow transect work (OPP-1820927), the sea ice coring work (OPP-1735862), the HELiX drone operations (OPP-1805569), surface energy fluxes (OPP-1724551), Climate Active Trace Gases (OPP-1807496), and Reactive Gas Chemistry (OPP-1914781). The last 4 of these were also supported by the NOAA Physical Sciences Laboratory,

• – the European Union’s Horizon 2020 research and innovation program projects ARICE (grant 730965) for berth fees associated with the participation of the DEARice team and INTAROS (grant 727890) supporting the drone and albedo measurements,

• – the US Department of Energy Atmospheric Radiation Measurement (ARM) and Atmospheric System Research (ASR) programs (DE-SC0019251, DE-SC0021341),

• – the National Aeronautics and Space Administration (NASA) project 80NSSC20K0658,

• – the European Space Agency (ESA) MOSAiC microwave radiometer (EMIRAD2, ELBARA, HUTRAD), (EMIRAD2, ELBARA, HUTRAD), CIMRex (contract 4000125503/18/NL/FF/gp) and GNSS-R (P.O. 5001025474, C.N. 4000128320/19/NL/FF/ab) GNSS-R (contracts P.O. 5001025474 and C.N. 4000128320/19/NL/FF/ab) projects,

• – the Canadian Space Agency FAST project (grant no. 19FACALB08),

• – EUMETSAT support for microwave scatterometer measurements,

• – the Research Council of Norway through the projects HAVOC (grant no. 280292), SIDRiFT (grant no. 287871), and CAATEX (grant no. 280531),

• – the Fram Centre (Tromsø, Norway), from its flagship program on Arctic Ocean through the PHOTA project,

• – the UKRI Natural Environment Research Council (NERC) and BMBF, who jointly funded the Changing Arctic Ocean program (project Diatom Arctic, NE/R012849/1 and 03F0810A),

• – the UK Natural Environment Research Council (project SSAASI-CLIM grant NE/S00257X/1),

• – the Agencia Estatal de Investigación AEI of Spain (grant no. PCI2019-111844-2, RTI2018-099008-B-C22),

• – the Swedish Research Council (VR, grant no. 2018-03859),

• – the Swedish Polar Research Secretariat for berth fees for MOSAiC,

• – the Swiss Polar Institute project SnowMOSAiC,

• – the Werner-Petersen-Foundation for the development of a remotely operated floating platform (grant no. FKZ 2019/610).

All authors declare that they have no competing interests.

Contributed to conception and design: KD, AAF, AH, BR, BoL, BrL, CG, CH, CP, DKP, EC, ED, GB, GdB, GS, JKH, JS, JuR, KVH, LK, MAG, MDS, MN, MR, MSchn, MS, NO, RG, RKS, SG, SH, SSK, TK.

Contributed to acquisition of data: AAF, AJ, AM, AN, AR, ARM, AT, AnT, BAL, BC, BH, BR, CH, CK, DC, DVD, DK, DKP, DNW, DW, ED, ES, FL, FP, GB, GS, HJB, H-RH, IAR, IO M, JB, JKH, JS, JaH, JaR, JoH, JuR, KVH, LK, LT, LvA, MA, MAG, MDS, MF, MH, MJ, MM, MMS, MN, MO, MR, MSchi, MSchn, MT, MeW, NF, NK, OD, PA, PI, RC, RD, RL, RN, RP, RR, RTT, SA, SF, SH, TK, TP, TPS, VB, VN, VaS, VlS.

Contributed to analysis and interpretation of data: AJ, ARM, BoL, CG, CH, CK, DC, DVD, DKP, EC, FL, GB, GdB, GEL, GS, HFG, IO M, JH, JaH, JuR, KR, LT, LvA, MAG, MMS, MN, MO, MS, MSchn, MeW, NF, PA, PI, RC, RD, RKS, RP, RN, SH, SM, TK.

Drafted and/or revised the article: All authors.

Approved the submitted version for publication: All authors.