As the sea‐ice modeling community is shifting to advanced numerical frameworks, developing new sea‐ice rheologies, and increasing model spatial resolution, ubiquitous deformation features in the Arctic sea ice are now being resolved by sea‐ice models. Initiated at the Forum for Arctic Modeling and Observational Synthesis, the Sea Ice Rheology Experiment (SIREx) aims at evaluating state‐of‐the‐art sea‐ice models using existing and new metrics to understand how the simulated deformation fields are affected by different representations of sea‐ice physics (rheology) and by model configuration. Part 1 of the SIREx analysis is concerned with evaluation of the statistical distribution and scaling properties of sea‐ice deformation fields from 35 different simulations against those from the RADARSAT Geophysical Processor System (RGPS). For the first time, the viscous‐plastic (and the elastic‐viscous‐plastic variant), elastic‐anisotropic‐plastic, and Maxwell‐elasto‐brittle rheologies are compared in a single study. We find that both plastic and brittle sea‐ice rheologies have the potential to reproduce the observed RGPS deformation statistics, including multi‐fractality. Model configuration (e.g., numerical convergence, atmospheric representation, spatial resolution) and physical parameterizations (e.g., ice strength parameters and ice thickness distribution) both have effects as important as the choice of sea‐ice rheology on the deformation statistics. It is therefore not straightforward to attribute model performance to a specific rheological framework using current deformation metrics. In light of these results, we further evaluate the statistical properties of simulated Linear Kinematic Features in a SIREx Part 2 companion paper.

The cumulative Greenland freshwater flux anomaly has exceeded 5,000 km³ since the 1990s. The volume of this surplus freshwater is expected to cause substantial freshening in the North Atlantic. Analysis of hydrographic observations in the subpolar seas reveals freshening signals in the 2010s. The sources of this freshening are yet to be determined. In this study, the relationship between the surplus Greenland freshwater flux and this freshening is tested by analyzing the propagation of the Greenland freshwater anomaly and its impact on salinity in the subpolar North Atlantic based on observational data and numerical experiments with and without the Greenland runoff. A passive tracer is continuously released during the simulations at freshwater sources along the coast of Greenland to track the Greenland freshwater anomaly. Tracer budget analysis shows that 44% of the volume of the Greenland freshwater anomaly is retained in the subpolar North Atlantic by the end of the simulation. This volume is sufficient to cause strong freshening in the subpolar seas if it stays in the upper 50–100 m. However, in the model the anomaly is mixed down to several hundred meters of the water column resulting in smaller magnitudes of freshening compared to the observations. Therefore, the simulations suggest that the accelerated Greenland melting would not be sufficient to cause the observed freshening in the subpolar seas and other sources of freshwater have contributed to the freshening. Impacts on salinity in the subpolar seas of the freshwater transport through Fram Strait and precipitation are discussed.

Simulating sea ice drift and deformation in the Arctic Ocean is still a challenge because of the multiscale interaction of sea ice floes that compose the Arctic Sea ice cover. The Sea Ice Rheology Experiment (SIREx) is a model intercomparison project of the Forum of Arctic Modeling and Observational Synthesis (FAMOS). In SIREx, skill metrics are designed to evaluate different recently suggested approaches for modeling linear kinematic features (LKFs) to provide guidance for modeling small‐scale deformation. These LKFs are narrow bands of localized deformation that can be observed in satellite images and also form in high resolution sea ice simulations. In this contribution, spatial and temporal properties of LKFs are assessed in 36 simulations of state‐of‐the‐art sea ice models and compared to deformation features derived from the RADARSAT Geophysical Processor System. All simulations produce LKFs, but only very few models realistically simulate at least some statistics of LKF properties such as densities, lengths, or growth rates. All SIREx models overestimate the angle of fracture between conjugate pairs of LKFs and LKF lifetimes pointing to inaccurate model physics. The temporal and spatial resolution of a simulation and the spatial resolution of atmospheric boundary condition affect simulated LKFs as much as the model's sea ice rheology and numerics. Only in very high resolution simulations (≤2 km) the concentration and thickness anomalies along LKFs are large enough to affect air‐ice‐ocean interaction processes.

Accelerating since the early 1990s, the Greenland Ice Sheet mass loss exerts a significant impact on thermohaline processes in the sub‐Arctic seas. Surplus freshwater discharge from Greenland since the 1990s, comparable in volume to the amount of freshwater present during the Great Salinity Anomaly events, could spread and accumulate in the sub‐Arctic seas, influencing convective processes there. However, hydrographic observations in the Labrador Sea and the Nordic Seas, where the Greenland freshening signal might be expected to propagate, do not show a persistent freshening in the upper ocean during last two decades. This raises the question of where the surplus Greenland freshwater has propagated. In order to investigate the fate, pathways, and propagation rate of Greenland meltwater in the sub‐Arctic seas, several numerical experiments using a passive tracer to track the spreading of Greenland freshwater have been conducted as a part of the Forum for Arctic Ocean Modeling and Observational Synthesis effort. The models show that Greenland freshwater propagates and accumulates in the sub‐Arctic seas, although the models disagree on the amount of tracer propagation into the convective regions. Results highlight the differences in simulated physical mechanisms at play in different models and underscore the continued importance of intercomparison studies. It is estimated that surplus Greenland freshwater flux should have caused a salinity decrease by 0.06–0.08 in the sub‐Arctic seas in contradiction with the recently observed salinification (by 0.15–0.2) in the region. It is surmised that the increasing salinity of Atlantic Water has obscured the freshening signal.

Five quantitative methodologies (metrics) that may be used to assess the skill of sea ice models against a control field are analyzed. The methodologies are Absolute Deviation, Root‐Mean‐Square Deviation, Mean Displacement, Hausdorff Distance, and Modified Hausdorff Distance. The methodologies are employed to quantify similarity between spatial distribution of the simulated and control scalar fields providing measures of model performance. To analyze their response to dissimilarities in two‐dimensional fields (contours), the metrics undergo sensitivity tests (scale, rotation, translation, and noise). Furthermore, in order to assess their ability to quantify resemblance of three‐dimensional fields, the metrics are subjected to sensitivity tests where tested fields have continuous random spatial patterns inside the contours. The Modified Hausdorff Distance approach demonstrates the best response to tested differences, with the other methods limited by weak responses to scale and translation. Both Hausdorff Distance and Modified Hausdorff Distance metrics are robust to noise, as opposed to the other methods. The metrics are then employed in realistic cases that validate sea ice concentration fields from numerical models and sea ice mean outlook against control data and observations. The Modified Hausdorff Distance method again exhibits high skill in quantifying similarity between both two‐dimensional (ice contour) and three‐dimensional (ice concentration) sea ice fields. The study demonstrates that the Modified Hausdorff Distance is a mathematically tractable and efficient method for model skill assessment and comparison providing effective and objective evaluation of both two‐dimensional and three‐dimensional sea ice characteristics across data sets.

Concurrent measurements by satellite tracked drifters of different hull and drogue configurations and coastal high-frequency radar reveal substantial differences in estimates of the near-surface velocity. These measurements are important for understanding and predicting material transport on the ocean surface as well as the vertical structure of the near-surface currents. These near-surface current observations were obtained during a field experiment in the northern Gulf of Mexico intended to test a new ultra-thin drifter design. During the experiment, thirty small cylindrical drifters with 5 cm height, twenty-eight similar drifters with 10 cm hull height, and fourteen drifters with 91 cm tall drogues centered at 100 cm depth were deployed within the footprint of coastal High-Frequency (HF) radar. Comparison of collocated velocity measurements reveals systematic differences in surface velocity estimates obtained from the different measurement techniques, as well as provides information on properties of the drifter behavior and near-surface shear. Results show that the HF radar velocity estimates had magnitudes significantly lower than the 5 cm and 10 cm drifter velocity of approximately 45% and 35%, respectively. The HF radar velocity magnitudes were similar to the drogued drifter velocity. Analysis of wave directional spectra measurements reveals that surface Stokes drift accounts for much of the velocity difference between the drogued drifters and the thin surface drifters except during times of wave breaking.