Floes Research Articles

Seasonal evolution of the sea ice floe size distribution in the Beaufort Sea from 2 decades of MODIS data

Abstract. Arctic sea ice cover evolves seasonally from large plates separated by long, linear leads in the winter to a mosaic of smaller sea ice floes in the summer. The interplay between physical and thermodynamic mechanisms during this process ultimately creates the observed sea ice floe size distribution (FSD), which is an important metric for characterizing the sea ice cover and assessing model performance. Historically, the FSD has been studied at fixed locations over short periods, leaving a gap in our understanding of the spatial and temporal evolution of the FSD at large scales. Here, we present an automated framework for image segmentation, allowing the identification and labeling of individual ice floes in Moderate Resolution Imaging Spectroradiometer (MODIS) data. Using this algorithm, we automatically process and segment 4861 images, identifying more than 9.4 million floes over 23 years. The extracted characteristics of the floes – including area, perimeter, and orientation – evolve throughout the spring and summer in the Beaufort Sea. We find seasonal patterns of decreasing mean floe areas, increasing FSD power law slopes, and increasing variability in the floe orientation as the summer progresses.

CSEU-Net: ConvNeXt-SE-U-Net for river ice floe segmentation using unmanned aerial vehicle grayscale remote sensing images

CSEU-Net: ConvNeXt-SE-U-Net for river ice floe segmentation using unmanned aerial vehicle grayscale remote sensing images

Angular variation of SAR polarimetric parameters over multiyear ice

ABSTRACT Backscatter and derived parameters from Synthetic Aperture Radar (SAR) may vary across the swath as the incident angle increases from the near to the far range, especially in the ScanSAR mode. Previous studies characterize this angular variation in the form of linear regression for commonly used ice types, namely first-year (FYI) and multi-year (MYI). The present study took advantage of the appearance of three grounded multi-year ice floes, viewed at different radar incidence angles, in a series of 22 RADARSAT-2 fully-polarimetric images, acquired over the Resolute Passage in the Canadian Arctic, during the period from 21 October to 28 December 2017. The ice properties remained invariant during the observation period; hence, variation of radar parameters of MYI can be attributed to change of incidence angle only. Quantification of angular variation of the selected SAR parameters was performed. They include conventional orthogonal backscatter coefficients, derived parameters from the incoherent decomposition of Cloude-Pottier, Yamaguchi, and Touzi, as well as a set of parameters from the compact polarimetry mode. Results specify the parameters that reveal variation with incidence angle and present the linear regressions equations. The explanation for the behaviour of each parameter, in terms of showing or not showing an angular trend, is offered. Results will be useful when using the polarimetric parameters in the sea ice classification scheme.

Three-dimensional discrete element simulations on pressure ridge formation

Abstract. This study presents the first three-dimensional discrete element method simulations of pressure ridge formation. Pressure ridges are an important feature of the sea-ice cover, as they contribute to the mechanical thickening of ice and likely limit the strength of sea ice at large scales. We validate the simulations against laboratory-scale experiments, confirming their accuracy in predicting ridging forces and ridge geometries. Then we demonstrate that Cauchy–Froude scaling applies for translating laboratory-scale results on ridging to full-scale scenarios. We show that non-simultaneous failure, where an ice floe fails at distinct locations across the ridge length, is required for an accurate representation of the ridging process. This process cannot be described by two-dimensional simulations. We also find a linear relationship between the ridging forces and the ice thickness, contrasting with earlier results in the literature obtained by two-dimensional simulations.

Non-smooth discrete element method analysis of channel width's effect on ice resistance in broken ice field

A numerical model based on NDEM is used to investigate the effect of channel width on broken-ice resistance in the channel. The model incorporates a limited impulse method to simulate the ice crushing failure. The accuracy of the numerical model is benchmarked against experimental data from model tests conducted on the Araon icebreaker at KRISO. The simulations successfully capture the characteristic effects of channel width on broken-ice resistance observed in the experiment. Beyond benchmarking of the numerical model, this study analyses the force chain characteristics in the broken ice field, revealing the crucial role of force chains in the resistance increase caused by channel width constraints. Building upon the established numerical model, further sensitivity analyses are conducted to investigate the correlation between channel width effect and other relevant factors such as ice concentration, ship velocity and ice thickness. Among these factors, it is found that ice concentration emerges as the predominant determinant of the critical channel width. The increase in ice resistance attributable to the channel width is exacerbated under the condition of thicker ice floes and lower ship velocity, where the occurrence of ice floe rotation is less frequent and the corresponding force chains are stronger and more stable.

High-resolution regional sea-ice model based on the discrete element method with boundary conditions from a large-scale model for ice drift

Abstract Understanding sea-ice dynamics at the floe scale is crucial to comprehend polar climate systems. While continuum models are commonly used to simulate large-scale sea-ice dynamics, they have limitations in accurately representing sea-ice behaviour at small scales. DEMs, on the other hand, are well-suited for modelling the behaviour of individual ice floes but face limitations due to computational constraints. To address the limitations of both approaches while combining their strengths, we explored the feasibility of using a DEM within a continuum model, where the latter provides boundary conditions for a rectangular high-resolution DEM domain. This paper presents a feasibility study where a discrete model of a 100 × 100 km2 icefield was created using high-resolution optical satellite imagery. Sea-ice dynamics were simulated in the DEM considering environmental forces and integrating large-scale ice-drift velocities as boundary conditions. Model predictions were compared with satellite observations for ice drift and deformation parameters. This numerical approach showed potential for offering accurate, high-resolution predictions of sea ice, particularly in coastal areas and near islands, and may find applications in ice navigation and climate studies. However, further development of the DEM, along with upgrades to the coupled ocean models providing data for the ice component, may be necessary. Additionally, challenges remain to develop a two-way coupling between the DEM and a continuum model, which may be needed to improve the accuracy of large-scale simulations.

A theoretical model of excess attenuation of acoustic signals propagating under cracked sea-ice landscapes.

The Arctic sheet is transitioning from a continuous cover of thick multi-year ice to a fragmented landscape of thin young ice. If the type of acoustic transmission allows repetitive interaction of rays with the sea surface, in the fragmented scenario acoustic rays will undergo a random sequence of reflections from water or sea-ice interfaces. Calm sea conditions in the water channels between the ice floes (leads) and the smooth, flat surface of the young ice bottom reduce scattering due to interface roughness, resulting only in scattering due to inhomogeneity in surface reflectivity. Using an idealized framework, this study investigates the extent to which the mid- to high-frequency underwater acoustic propagation is altered due to repetitive interactions of acoustic signals with a sea surface consisting of a random distribution of ice sheets and leads. An expression for the coherent field (the acoustic field averaged over an ensemble of realizations of sea-ice distributions) was derived from theory. Any deviation from a homogeneous surface condition (either by randomly adding ice slabs in a free ice surface or by including leads in a fully ice-covered sea surface) leads to an excess attenuation of the coherent field. Results are validated by numerical simulations.

A Probabilistic Approach to Correct for Oblique Biases in Sea Ice Parametrization from Shipborne Cameras

Abstract This work investigates the perception bias caused by ice floes occluding open water regions behind them when analyzing sea ice from shipborne cameras. We propose a method to model this effect by predicting an occlusion probability for a given surface region based on a stochastic representation of the sea ice fields. We further show how the proposed framework can be used to design an unbiased estimator for image-based sea ice concentration by correcting for these occlusions. The occlusion model and proposed sea ice concentration estimator have been verified through simulated experiments. Results from the experiments showed that predictions from the occlusion model align with the simulated observations and that the proposed estimator remains unbiased for all environments tested. The contributions of this study can be used in algorithms for automated identification and characterization of ice environments from surface ships and are thus an essential step towards improved autonomy in Arctic navigation.

Notes on Towed Self-Propulsion Experiments with Simulated Managed Ice in Traditional Towing Tanks

Efficiency estimation of a propeller behind a vessel’s hull while sailing through ice floes, together with the ship’s resistance to motion, is a key factor in designing the power plant and determining the safety measures of a ship. This paper encloses the results from the experiments conducted at the CEHINAV towing tank, which consisted of analyzing the influence of the concentration at the free surface of artificial blocks, simulating ice, in propeller–block interactions. Thrust and torque were measured for a towed self-propelled ship model through simulated broken ice blocks made of paraffin wax. Three block concentrations of different block sizes and three model speeds were studied during the experimentation. Open-water self-propulsion tests and artificial broken ice towed self-propulsion tests are shown and compared in this work. The most relevant observations are outlined at the end of this paper, as well as some guidelines for conducting artificial ice-towed self-propulsion tests in traditional towing tanks.

Calibration of a hybrid sea ice model during an expedition to the Arctic

AbstractWe describe the acquisition and evaluation of data during the ArcWatch II expedition into the Arctic Ocean on board of the icebreaker Polarstern. The data is used to calibrate a hybrid sea ice model that focuses on modeling the marginal ice zone, which is the transition zone between pack ice and the open ocean. We describe the experimental setup used to track the motion of individual ice floes as well as the measurement of the wind and the ocean current. In the literature, the motion of an ice floe in the free‐drift regime is often estimated by Nansen's rule of thumb, which assumes the sea ice drifts at of the wind speed, deflected up to in clockwise direction. Analyzing the first datasets, we noticed that observed small ice floes in free drift often have a larger relative angle to the wind direction. Based on numerical simulation we quantify that these larger angles are related to the geometry of larger sea ice floes in the vicinity of the observed ice floe. These large icefloes influence the ocean current, which in turn, impacts the motion of the small ice floes in the surrounding.

Submicron Organic Aerosol Types in the Summertime Arctic: Mixing State, Geographic Distribution, and Drivers

AbstractDuring the 2017 summertime Arctic cruise observation campaigns, we measured over 290,000 individual submicron particles and clustered them into two inorganic classes (dominated by sea salt, accounting for 38.6% by number fraction) and five organic classes (dominated by natural and anthropogenic organics, 61.4%), presenting a distinct difference in geographic distribution. In the high Arctic and marginal ice zone (81.1–84.6°N) compared with the low Arctic (Chukchi Sea, Svalbard, and Iceland, <80°N), ocean‐derived organic aerosols were more prevalent (73.6% vs. 37.1%). Specifically, we found sharp contrasts in the geographic distributions of OC‐Ca (organics internally mixed with calcium, 29.0% vs. 9.4%) and OC‐S (organics internally mixed with sulfate, 3.2% vs. 21.4%). Utilizing an explainable machine learning technique, we inferred that OC‐Ca was driven by wind‐blown sea ice and/or sea ice floes and/or bubble bursting within sea ice leads under low wind speed conditions in the high Arctic, while OC‐S tended to associate with elemental carbon, sulfate, and higher temperatures, potentially originating from combustion emissions at low latitude regions.

Multiscale mushy layer model for Arctic marginal ice zone dynamics

Perhaps the most dynamic component of the Arctic sea ice cover is the marginal ice zone (MIZ), the transitional region between dense pack ice to the north and open ocean to the south. It widens by a factor of four while seasonally migrating more than 1600 km poleward in the Bering–Chukchi Sea sector, impacting climate dynamics, ecological processes, and human accessibility to the Arctic. Here we showcase a transformative mathematical modeling approach to understanding changes in MIZ location and width, focusing on their seasonal cycles as observed by satellites. We view the MIZ as a liquid-solid phase transition region, or mushy zone, on the scale of the Arctic Ocean. Invoking the physics of phase changes, the MIZ is modeled as a dynamic, multiscale composite material layer; this model captures 96% of the annual cycle of MIZ location and 78% of the annual cycle of MIZ width. Temperature in the upper ocean is described by a nonlinear heat equation with effective parameters obtained using homogenization theory for a random medium of ice floes in a sea water host. Observations and simulations together indicate that MIZ location closely tracks the below-ice 273 K isotherm while the width of the MIZ follows vertical heat flux convergence, but with a three-week lag.

The force action of ice cakes to stationary and floating objects

In freezing seas during winter periods between open water and cover by ice, transitional zones are formed from fragments of level ice fields (ice cakes). The quantity and sizes of these fragments decrease as they approach to open water. The ice cakes are drifting due to action of winds and sea currents. During storm periods, being on the stormy surface of the seas, ice cakes are capable of exerting the dynamic loads on stationary and floating objects. At the initial moment of contact of ice floes with various obstacles, significant local pressures are realized оn small contact areas, which represent a danger to the integrity of the structures of such obstacles. In the article considered an approach to solving the task of the dynamic impact by fragments of level ice fields on the supporting parts of offshore hydraulic structures that will be function in the conditions of open sea during severe storms in the winter and spring periods.

Research on the ice-resistance performance of conical wind turbine foundation in brash ice fields

Due to the significantly higher wind energy density in high-latitude sea areas compared to other regions, the construction of offshore wind farms is shifting towards these areas. However, in high-latitude sea areas, large amounts of sea ice float on the surface during winter, making ice loads a critical factor affecting the structural safety of wind turbines. For monopile flexible wind turbines, installing an ice-breaking cone near the waterline is an effective measure to withstand the threat of sea ice. To clarify the ice resistance performance of conical wind turbine foundations in brash ice zones, this study proposes a coupled Computational Fluid Dynamics-Discrete Element Method (CFD-DEM). This method is used to calculate the ice load results for a cylindrical riser, and the results were found to be in good agreement with the ice load experiment results from China University of Petroleum's wave flume (CUPWF). This validation confirms the rationality of the method. Taking a 5 MW monopile wind turbine installed with an ice-breaking cone as the research object, the time history of dynamic ice loads and the displacement response of ice-induced vibration during the interaction between the wind turbine and brash ice are simulated. Considering the water level variations caused by tides, the differences in ice resistance performance between the upright cone and the inverted cone of the ice-breaking cone are analyzed. Finally, the impact of changes in cone angle on the ice resistance performance of the wind turbine foundation is determined. The results indicate that under the impact of brash ice, the upright cone of the ice-breaking cone has a higher ice resistance capacity compared to the inverted cone. With changes in cone angle, the trends of ice loads in the horizontal and vertical directions vary oppositely. However, the overall ice resistance capacity of the wind turbine foundation decreases as the cone angle increases. This study can provide a reference for the rational design of conical wind turbines in cold regions.

Revealing the Freezing-Induced Alteration in Microplastic Behavior and Its Implication for the Microplastics Released from Seasonal Ice.

Ice can serve as a significant temporary repository and conveyance mechanism for microplastics (MPs). MPs present in the water column can become entrapped within developing ice formations, subsequently being sequestered and transported by ice floes. With changing temperatures, MPs stored in ice can be released back into the environment, while freezing conditions can alter the properties of MPs, ultimately affecting the fate of MPs in the environment. Freezing of MPs in freshwater ice results in the aggregation of MP particles due to physical compression, leading to an increase in particle size once the MPs are released from the ice. The freezing-induced aggregation enhances buoyancy effects, accelerating the settling/rising velocity of MPs in water. Additionally, freezing can lead to enhanced surface wetting alterations, thus improving the dispersion of hydrophobic MPs. The presence of salt in the water can mitigate the effect of freezing on MPs due to the formation of a brine network within the ice structure, which reduces the pressure on MPs entrapped by ice. In cold regions, numerous MPs undergo freezing and thawing, re-entering the water column.

An improved physical information network for forecasting the motion response of ice floes under waves

An improved physical information network for forecasting the motion response of ice floes under waves

Sea Ice Deformation Is Not Scale Invariant Over Length Scales Greater Than a Kilometer

AbstractIn March and April 2021, buoys were deployed in the Beaufort Sea, Arctic Ocean, to measure sea‐ice horizontal deformation over spatial scales that had not been previously achieved. Geodetic‐quality position measurements allowed measurements of strain‐rate over lengths from about 200 m to 2 km. Conventional ice‐drifters extended spatial coverage up to about 100 km. Past studies find there is multi‐fractal behavior for horizontal sea‐ice deformation from 10 to 1,000 km. Our results demonstrate that such behavior does not hold when including spatial scales below 10 km. We find that sea‐ice deformation is not scale invariant between the scale of individual sea‐ice floes and aggregates of floes. Therefore, we cannot expect the same physical laws or forcing to describe sea‐ice kinematics over these regimes, nor can we assume log‐log linear behavior for mean deformation. Using this scaling behavior as a metric to validate models that resolve sea ice floes and their interactions is hence not recommended.

Ice-Jam Investigations along the Oder River Based on Satellite and UAV Data

The Oder River, situated along the border between Poland and Germany, is regularly affected by ice-jam events and their associated hazards, such as a sudden rise in water level and the endangerment to flood-protection infrastructure. The existing databases on past ice-jam events lack substantial information considering ice formation, blockage origins or the spatiotemporal evolution of the ice cover needed for a comprehensive understanding of relevant ice processes. Within this study, the evaluation of satellite and Uncrewed Aerial Vehicle (UAV) data was carried out in order to analyze the capabilities of enhancing river ice information in the study area. Satellite imagery was proven to be a valuable source of investigating ice-jam phenomena on all scales, leading to the identification of initial ice-jam locations, surveying spatiotemporal ice cover evolution or monitoring the maximum ice-cover extent. A simplified approach for river ice classification of satellite radar data using the K-Means Cluster Analysis is introduced, enabling the differentiation between river ice formations. Based on UAV data taken in this study, workflows were presented, allowing for measurements of ice floe velocities and the localization of flooded and ice-covered flow control structures.

Numerical investigation of riser vibration induced by ice-water two-phase flow in icy regions by SPH method

In icy regions, the vibration induced by broken ice poses a hazard to risers, in addition to vortex-induced vibration (VIV). A numerical investigation has been conducted to explore the dynamic response of risers in seas with fragmented ice floes. Assuming the broken ice as a single rigid floating body bonded by particles, the flow with water and ice is considered a two-phase flow, which is established using the Smooth Particle Hydrodynamics (SPH) method based on the quasi-fluid model. Verification of the model was achieved by simulating the three-dimensional flow around a surface-piercing cylinder and VIV. Accurate predictions were made regarding hydrodynamic forces, uplift values of the free liquid surface at the front stagnation point of the cylinder, and riser motion. Furthermore, the SPH model was expanded to encompass the realm of ice-water structural interaction (IWSI) combined with the basic principle of bidirectional fluid-structure coupling, facilitating the simulation of riser vibration within ice-water flows. The results indicate that the ice load remained the principal contributor to the force on the vibrating riser, while the riser exhibited irregular motion, demonstrating a strong coupling between forces and riser movement.

An improved current retrieval method for ice-edge eddies from multi-sensor remote sensing image sequences

ABSTRACT Eddies in the Marginal Ice Zones (MIZ) contribute to an enhanced melting of sea ice by forcing contact between the sea ice and the warmer water off the ice edge. The horizontal current field can be retrieved by tracking the movement of ice floes between sequential satellite images. However, the traditional method of current retrieval imposes stringent requirements on the time intervals between sequential images, and the performance of the traditional method diminishes rapidly with increasing time intervals. In this paper, an improved method of retrieving currents from multi-sensor remote sensing images is proposed, which can significantly improve the accuracy of current retrieval. The method first extracts feature from remote sensing images based on texture features and a support vector machine model and then extracts prior information about the ice-edge eddy from the first two remote sensing images. Under the guidance of prior information, horizontal currents of the eddy are retrieved by tracking the movement of sea ice from sequential remote sensing images. The proposed method was used to retrieve currents from Sentinel-1 and Sentinel-2 sequential images over the MIZ of the Fram Strait. Comparisons with the traditional method show encouraging results that the proposed method can significantly reduce the number of mismatches and improve the reliability of current vectors.