A viscoelastic-plastic constitutive model with Mohr-Coulomb yielding criterion for sea ice dynamics

A new viscoelastic-plastic （VEP） constitutive model for sea ice dynamics was developed based on continuum mechanics. This model consists of four components： Kelvin-Vogit viscoelastic model, Mohr-Coulomb yielding criterion, associated normality flow rule for plastic rehololgy, and hydrostatic pressure. The numerical simulations for ice motion in an idealized rectangular basin were made using smoothed particle hydrodynamics （SPH） method, and compared with the analytical solution as well as those based on the modified viscous plastic（VP） model and static ice jam theory. These simulations show that the new VEP model can simulate ice dynamics accurately. The new constitutive model was further applied to simulate ice dynamics of the Bohai Sea and compared with the traditional VP, and modified VP models. The results of the VEP model are compared better with the satellite remote images, and the simulated ice conditions in the JZ20-2 oil platform area were more reasonable.

A modified discrete element model for sea ice dynamics

Considering the discontinuous characteristics of sea ice on various scales,a modified discrete element model（DEM） for sea ice dynamics is developed based on the granular material rheology.In this modified DEM,a soft sea ice particle element is introduced as a self-adjustive particle size function.Each ice particle can be treated as an assembly of ice floes,with its concentration and thickness changing to variable sizes under the conservation of mass.In this model,the contact forces among ice particles are calculated using a viscous-elastic-plastic model,while the maximum shear forces are described with the Mohr-Coulomb friction law.With this modified DEM,the ice flow dynamics is simulated under the drags of wind and current in a channel of various widths.The thicknesses,concentrations and velocities of ice particles are obtained,and then reasonable dynamic process is analyzed.The sea ice dynamic process is also simulated in a vortex wind field.Taking the influence of thermodynamics into account,this modified DEM will be improved in the future work.

Mechanical and numerical models for sea ice dynamics on small-meso scale

On small-meso scale, the sea ice dynamic characteristics are quite different from that on large scale. To model the sea ice dynamics on small-meso scale, a new elastic-viscous-plastic （EVP） constitutive model and a hybrid Lagrangian- Eulerian （HLE） numerical method are developed based on continuum theory. While a modified discrete element model （DEM） is introduced to model the ice cover at discrete state. With the EVP constitutive model, the numerical simulation for ice ridging in an idealized rectangular basin is carried out and the results are comparable with the analytical solution of jam theory. Adopting the HLE numerical model, the sea ice dynamic process is simulated in a vortex wind field. The furthering application of DEM is discussed in details for modeling the discrete distribution of sea ice. With this study, the mechanical and numerical models for sea ice dynamics can be improved with high precision and computational efficiency.

Modified PIC Method for Sea Ice Dynamics

The sea ice cover displays various dynamical characteristics such as breakup, rafting, and ridging under external forces. To model the ice dynamic process accurately, the effective numerical modeling method should be established. In this paper, a modified particle-in-cell （PIC） method for sea ice dynamics is developed coupling the finite difference （FD） method and smoothed particle hydrodynamics （SPH）. In this method, the ice cover is first discretized into a series of lagrangian ice particles which have their own sizes, thicknesses, concentrations and velocities. The ice thickness and concentration at Eulerian grid positions are obtained by interpolation with the Gaussian function from their surrounding ice particles. The momentum of ice cover is solved with FD approach to obtain the Eulerian cell velocity, which is used to estimate the ice particle velocity with the Gaussian function also. The thickness and concentration of ice particles are adjnsted with particle mass density and smooth length, which are adjusted with the redistribution of ice particles. With the above modified PIC method, numerical simulations for ice motion in an idealized rectangular basin and the ice dynamics in the Bohai Sea are carried out. These simulations show that this modified PIC method is applicable to sea ice dynamics simulation.

DRUCKER-PRAGER YIELD CRITERIA IN VISCOELASTIC-PLASTIC CONSTITUTIVE MODEL FOR THE STUDY OF SEA ICE DYNAMICS

Based on the characteristics of sea ice drifting and ridging at meso-small scale, the Drucker-Prager （D-P） yield criteria was introduced into the Viscoelastic-Plastic （VEP） constitutive model for the study of sea ice dynamics. In this model, the Kelvin-Vogit viscoelastic model was adopted in the elastic stage, and the associated normal flow rule was used in the plastic stage. Using the VEP model, the sea ice ridging process was simulated in an idealized rectangular basin, and the simulation results show that the simulated ice ridge thickness is consistent with the analytical solution. Moreover, the VEP model with the D-P yield criteria was also applied for the sea ice simulation of Bohai Sea, and the ice thickness, concentration, velocity, and ice stress were obtained in 48 h. The simulated thickness distributions agree well with the satellite images. The singular problem in the Mohr^7oulomb （M-C） yield criteria was overcome by the D-P yield criteria, and the computational efficiency was also improved. In the numerical simulations described above, the smoothed particle hydrodynamics was applied.

DynIceData:a gridded ice-water classification dataset at short-time intervals based on observations from multiple satellites over the marginal ice zone

High-resolution observations of short-term changes in sea ice are critical to understanding ice dynamics and also provide important information used in advice to shipping,especially in the Arctic.Although individual satellite sensors provide periodic sea ice obser-vations with spatial resolutions of tens of meters,information regarding changes that occur over short time intervals of minutes or hours is limited.In this study,a gridded ice-water classification dataset with a high temporal resolution was developed based on observations acquired by multiple satellite sensors in the Marginal Ice Zone(MIZ).This dataset-DynIceData-which combines Sentinel-1 Synthetic Aperture Radar(SAR)data with Gaofen-3(GF-3)SAR and SDGSAT-1 thermal infrared imagery was used to obtain observations of the MIZ with a range of temporal resolutions ran-ging from minutes to tens of hours.The areas of the Arctic covered include the Kara Sea,Beaufort Sea,and Greenland Sea during the period from August 2021 to August 2022.Object-oriented segmen-tation and thresholding were used to obtain the ice-water classifi-cation map from Sentinel-1 and GF-3 SAR image pairs and Sentinel-1 SAR and SDGSAT-1 thermal image pairs.The time interval between the images in each pair ranged from 1 minute to 68 hours.Ten-kilometer grid sample granules with a spatial resolution of 25 m for the GF-3 SAR data and 30 m for the SDGSAT-1 thermal data were used.The classification was verified as having an overall accuracy of at least 95.58%.The DynIceData dataset consists of 7338 samples,which could be used as reference data for further research on rapid changes in sea ice patterns at different short time scales and provide support for research on thermodynamic and dynamic models of sea ice in combination with other environmen-tal data,thus potentially improving the accuracy of sea ice forecast-ing using Artificial Intelligence.The dataset can be accessed at https://doi.org/10.57760/sciencedb.j00001.00784.