冰脊形成过程的离散元模拟及影响因素分析

冰脊对极地船舶及海洋工程结构的冰载荷设计及冰区安全运行具有重要的影响,其几何形态是重要的考虑因素。为研究冰脊的几何形态,本文基于Voronoi切割算法构造扩展多面体海冰单元,并使其相互冻结形成平整冰;采用离散单元法模拟海冰在相对运动时挤压形成冰脊的动力过程,统计分析脊帆高度、龙骨宽度、龙骨深度、龙骨水平倾角等主要几何参数,并确定这些几何参数间的对应关系。将离散元模拟结果同冰脊现场观测资料对比分析以验证扩展多面体离散元方法模拟冰脊形成过程的可行性。在此基础上讨论冰厚、冰速、海冰强度等因素对冰脊形态的影响。本文采用离散元方法对冰脊形成过程及几何特性的数值分析可为深入开展冰脊的形成机理及其对船舶、海洋工程结构物的冰载荷分析提供参考依据。     

Some parameters in arctic sea ice dynamics from an expedition in the summer of 2003

1Introduction Seaiceoccupiesthemainpartofthesurfaceof theArcticOcean.ThefocusoftheSecondChineseNa- tionalArcticResearchExpedition(CHINAE-2003) wastounderstandthevariationsofarcticmarineenvi- ronmentsandtheseaiceeffectsontheclimatechanges ofglobalextent,inmiddleandlowerlatitudesareas, especiallyinChina.Therefore,thejointsea-ice-airob- servationforseaicestudieswasoneofthekeypro- jectsinCHINARE-2003.Theinvestigatedareacov- ered3000kmfromsouthtonorthand900kmfrom westtoeast.Seventemporali…     

GPU-Based DEM Simulations of Global Ice Resistance on Ship Hull During Navigation in Level Ice

The ice resistance on a ship hull affects the safety of the hull structure and the ship maneuvering performance in icecovered regions.In this paper,the discrete element method(DEM)is adopted to simulate the interaction between level ice and ship hull.The level ice is modeled with 3D bonded spherical elements considering the buoyancy and drag force of the water.The parallel bonding approach and the de-bonding criterion are adopted to model the freezing and breakage of level ice.The ship hull is constructed with rigid triangle elements.To improve computational efficiency,the GPU-based parallel computational algorithm was developed for the DEM simulations.During the interaction between the ship hull and level ice,the ice cover is broken into small blocks when the interparticle stress approaches the bonding strength.The global ice resistance on the hull is calculated through the contacts between ice elements and hull elements during the navigation process.The influences of the ice thickness and navigation speed on the dynamic ice force are analyzed considering the breakage mechanism of ice cover.The Lindqvist and Riska formulas for the determination of ice resistance on ship hull are employed to validate the DEM simulation.The comparison of results of DEM,Lindqvist,and Riska formula show that the DEM result is between those the Lindqvist formula and Riska formula.Therefore the proposed DEM is an effective approach to determine the ice resistance on the ship hull.This work can be aided in the hull structure design and the navigation operation in ice-covered fields.     

北极巴伦支海和拉普捷夫海海冰中细菌和微藻的垂直分布

通过计算机图像测定系统测定、计算获得海冰细菌的大小、面积和碳含量,研究了北极巴伦支海和拉普捷夫海的浮冰中细菌大小、丰度和生物量的垂直分布及其与海冰微藻生物量的关系.细菌在冰心中是非均匀分布的,至少有两种不同分布类型:第一类有单一的丰度峰值,在冰心底部或在冰心中部;第二类有两个丰度峰值.海冰中细菌的丰度为0.4×105～36.7×105个/cm3.不同类型海冰中的细菌大小变化极大,在多年冰中,最大的细菌与叶绿素a最大值出现在同一层,而在一年冰中,细菌大小几乎没有垂直变化.整柱冰心的细菌生物量变化为19.2～79.2mg/cm2,细菌与海冰微藻生物量之比为0.43～10.00.对固定冰比较和研究的结果发现,海冰冰心中细菌大小、丰度和生物量的垂直分布差别极大.据此,分析了对目前海冰研究采样方法的局限性,并提出了规范采样方法的设想.     

Arctic sea ice distribution in summer based on aerial photos

1Introduction TheArcticOceanisoneoftheimportantcold sourcesontheearth,whichaffectsglobalclimateand oceancirculationseriously.Itsinteractionwithglobal climatesystemisrepresentedbyseaice,whichisthe mainfeatureonthesurfaceoftheArcticOcean(Aa- gaard,etal.,1989).Firstly,seaiceplaysapivotalrole intheheatandmassbalanceonthesurfaceoftheArc- ticOcean.Seaicenotonlyobstructstheheatexchange betweenatmosphereandocean,butalsoreflectsthe mostofthelocalsolarradiationbacktotheatmo- spherebecauseofitshighalb…     

Improvement of short-term sea ice forecast in the southern Okhotsk Sea

In this study, a numerical model of 7-day forecast of sea ice produced by the Japan Meteorological Agency was improved by the following approaches. First, a new ice dynamic model was introduced: the distributed mass/discrete floe model. The model takes account of discrete characteristics of ice floes and well simulates the ice edge location at low computational cost. Secondly, the grid size was reduced to 5 × 5 km for the future high resolution forecasts. Next, the sea surface current data was examined because it significantly influences sea ice movement. We applied two new datasets of HINO and Okhotsk Ocean General Circulation Model (Okhotsk OGCM), which are estimated by numerical simulations, for the 7-day forecast of sea ice. Ice southward speed in January and the whorl formations in February and March were well reproduced with Okhotsk OGCM datasets. Finally, the ocean heat flux at the ice-ocean interface was refined. As a result, we achieved an ice edge error reduction from 30.8 km to 23.5 km.     

Effects of hull–ice friction coefficient on predictions of pack ice forces for moored offshore vessels

Designers and offshore operators frequently predict pack ice loading on offshore vessels by conducting scale model tests. One factor that can affect pack ice loading is the hull–ice friction coefficient. This research investigates the effect of hull–ice friction coefficients for a moored offshore vessel model and includes ice floe size and ice concentration as additional variables. A method of non-dimensional analysis is modified in order to deal with the multivariate nature of the new data. The resulting non-dimensional equation provides insight on relationships between the predicted pack ice force and the variables under investigation. The relationship between pack ice force and hull–ice friction coefficient is shown to be approximately a fourth root function, while the relationship for ice floe size, non-dimensionalised by the vessel beam, is approximately linear. The relationship between predicted pack ice force and ice concentration exists in a band bounded by cubic and sixth power curves. Applying the modified equation to the previous data sets shows the current analysis slightly improves the normalization of pack ice forces.     

渤海海冰漂移过程的数值模拟和试验

建立了一个包含潮流作用的准定常海冰动力学模式，利用实测风资料和计算的潮流场对辽东湾中部的冰块漂移过程进行数值模拟，模拟的冰块漂移过程和实况基本一致。表明模式具有反映冰漂移过程动力特征的能力。通过对各动力因子的数值试验，说明引入潮流作用的必要性，并分析了各动力因子在冰漂移过程中的作用。     

浮冰在规则波和双色波下的运动特性研究

为研究浮冰在波浪下的运动特性,以单个聚氯乙烯（PVC）塑料圆板模拟浮冰,进行了一系列波浪试验研究。在规则波试验基础上,同时研究了双色波条件下的浮冰运动特性。结果显示：规则波条件下浮冰运动幅频响应算子首先随无量纲波长的增大而增大,当无量纲波长大于4.0后,运动幅频响应算子随无量纲波长的变化不明显,且趋于一定值;在双色波条件下,对应频率组成成分的浮冰运动幅频响应算子与规则波条件下随波长的变化规律一致。根据部分规则波试验结果提出预测浮冰慢漂速度的经验公式,并用已有的试验和余下的规则波与双色波试验结果进行验证。结果表明,经验公式对规则波和双色波条件下的浮冰慢漂速度的误差在20%以内,预测结果吻合较好。     

Amodified 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.     

A linear Boltzmann equation to model wave scattering in the marginal ice zone

We present a linear Boltzmann equation to model wave scattering in the Marginal Ice Zone (the region of ocean which consists of broken ice floes). The equation is derived by two methods, the first based on Meylan et al. [Meylan, M.H., Squire, V.A., Fox, C., 1997. Towards realism in modeling ocean wave behavior in marginal ice zones. J. Geophys. Res. 102 (C10), 22981–22991] and second based on Masson and LeBlond [Masson, D., LeBlond, P., 1989. Spectral evolution of wind-generated surface gravity waves in a dispersed ice field. J. Fluid Mech. 202, 111–136]. This linear Boltzmann equation, we believe, is more suitable than the equation presented in Masson and LeBlond [Masson, D., LeBlond, P., 1989. Spectral evolution of wind-generated surface gravity waves in a dispersed ice field. J. Fluid Mech. 202, 111–136] because of its simpler form, because it is a differential rather than difference equation and because it does not depend on any assumptions about the ice floe geometry. However, the linear Boltzmann equation presented here is equivalent to the equation in Masson and LeBlond [Masson, D., LeBlond, P., 1989. Spectral evolution of wind-generated surface gravity waves in a dispersed ice field. J. Fluid Mech. 202, 111–136] since it is derived from their equation. Furthermore, the linear Boltzmann equation is also derived independently using the argument in Meylan et al. [Meylan, M.H., Squire, V.A., Fox, C., 1997. Towards realism in modeling ocean wave behavior in marginal ice zones. J. Geophys. Res. 102 (C10), 22981–22991]. We also present details of how the scattering kernel in the linear Boltzmann equation is found from the scattering by an individual ice floe and show how the linear Boltzmann equation can be solved straightforwardly in certain cases.     

2010-2018年北极夏季中国北极科学考察航行期间被动微波遥感海冰密集度与船基目视观测资料的比较

为了更有效地将卫星数据应用于北极航行导航，被动微波（PM）产品的海冰密集度（SIC）与从中国北极科学考察中收集到的船基目视观测（OBS）资料进行了比较。在2010、2012、2014、2016和2018年的北极夏季总共收集了3667组目测数据。PM SIC取自基于SSMIS传感器的NASA-Team（NT）、Bootstrap（BT）以及Climate Data Record（CDR）算法和基于AMSR-E/AMSR-2传感器的BT、enhanced NT（NT2）以及ARTIST Sea Ice（ASI）算法。使用PM SIC的日算术平均值和OBS SIC的日加权平均值进行比较。比较了PM SIC和OBS SIC之间的相关系数，偏差和均方根偏差，包括总体趋势以及在轻度/普通/严重冰况下的情况。使用OBS数据，浮冰尺寸和冰厚对不同PM产品SIC反演的影响可以通过计算浮冰尺寸编码和冰厚的日加权平均值来评估。我们的结果显示相关系数的范围为0.89（AMSR-E/AMSR-2 NT2）到0.95（SSMIS NT），偏差的范围为-3.96%（SSMIS NT）到12.05%（AMSR-E/AMSR-2），均方根偏差的范围为10.81%（SSMIS NT）到20.15%（AMSR-E/AMSR-2 NT2）。浮冰尺寸对PM产品的SIC反演有显著的影响，大多数PM产品倾向于在小浮冰尺寸情况下低估SIC，而在大浮冰尺寸情况下高估SIC。超过30 cm的冰厚对于PM产品的SIC反演没有明显影响。总体来看，在北极夏季，SSMIS NT SIC与OBS SIC之间有着最好的一致性，而AMSR-E/AMSR-2 NT2 SIC与OBS SIC的一致性最差。     

基于黏聚单元法计算分析海洋平台锥形立柱冰载荷

基于非线性有限元分析方法和LS-DYNA计算软件,使用黏聚单元法建立层冰的数值模型,并进行了层冰与锥形立柱相互作用的数值模拟,研究抗冰锥角、冰厚等参数对立柱所受冰载荷的影响。模拟结果与实测和试验结果吻合良好。研究结果表明:使用黏聚单元法建立的层冰模型可以较好地模拟出层冰的破碎和堆积过程,冰载荷计算结果与ISO冰力标准结果以及离散单元法计算结果吻合较好;锥形立柱的抗冰锥角大小对冰载荷影响明显,随着锥角增大,水平方向的最大冰力增大而竖直方向上最大冰力减小;锥形结构改变了碰撞过程中层冰的破坏形式;结构物所受冰力大小与层冰厚度呈正相关,且增加速度接近线性。     

<Emphasis Type="Italic">In-situ</Emphasis> measured primary productivity of ice algae in Arctic sea ice floes using a new incubation method

Recent changes in climate and environmental conditions have had great negative effects such as decreasing sea ice thickness and the extent of Arctic sea ice floes that support ice-related organisms. However, limited field observations hinder the understanding of the impacts of the current changes in the previously ice-covered regions on sea ice algae and other ice-related ecosystems. Our main objective in this study was to measure recent primary production of ice algae and their relative contribution to total primary production (ice plus pelagic primary production). In-situ primary productivity experiments with a new incubation system for ice algae were conducted in 3 sea ice cores at 2 different ice camps in the northern Chukchi Sea, 2014, using a ¹³C and ¹⁵N isotope tracer technique. A new incubation system was tested for conducting primary productivity experiments on ice algae that has several advantages over previous incubation methods, enabling stable carbon and nitrogen uptake experiments on ice algae under more natural environmental conditions. The vertical C-shaped distributions of the ice algal chl-a, with elevated concentrations at the top and bottom of the sea ice were observed in all cores, which is unusual for Arctic sea ice. The mean chl-a concentration (0.05 ± 0.03 mg chl-a m^?3) and the daily carbon uptake rates (ranging from 0.55 to 2.23 mg C m^?2 d^?1) for the ice algae were much lower in this study than in previous studies in the Arctic Ocean. This is likely because of the late sampling periods and thus the substantial melting occurring. Ice algae contributed 1.5–5.7% of the total particulate organic carbon (POC) contents of the combined euphotic water columns and sea ice floes. In comparison, ice algae contributed 4.8–8.6% to the total primary production which is greater than previously reported in the Arctic Ocean. If all of the ice-associated productions were included, the contributions of the sea ice floes to the total primary production would be greater in the Arctic Ocean and their importance would be greater in the arctic marine ecosystems.     

冰激振动及冰荷载动力特性的试验研究

回顾了解决海冰引起结构振动分析方法的要求，通过室内冰力模型试验得出了冰激振动的重要参数－冰力频率，幅值的特性。研究表明，在低冰速时，冰力频率与冰速对冰厚之比呈线性关系；高冰速时，冰力步与冰速对冰厚之比呈现恒定性，文章还根据试验数据以及统计分析，推荐了冰力幅值计算公式和取值范围。     

Relation of ice conditions to climate change in the Bohai Sea of China

1NTRODUcrIONThe bohai As is a seasonally ice-covered sea and is located in the lowest latitudes (37' -4l'N), where sea ice occurs. The bohai ffea is nearly enclosed by land in the south, the northand the west, and only connects to the Huanghai ffea through the bohai Strait in the east.The width of the strait is abeut l06 km. The boai ffea is very shallow basin with the meandepth of l8 m and the maximum depth of 78 m. The topography of the sea bottom and thecoastal regions has an importan…     

Some dynamics of spilled oil in a fractured ice field in Buzzards Bay, Massachusetts

A spill of #2 fuel oil occurred 28 January 1977 in a dynamic ice field in BBuzzards Bay, Massachusetts. An oil budget study revealed that within four days after the spill occurred, approximately 60% of the oil was incorporated into the fractured ice zones and onto the ice surface. The movement of the oil (initially spilled under the ice) into the fracture zones appears to be the result of water current below the ice and buoyancy spreading mechanisms. Areas of fracture zones and leads provided a surface outlet for the less dense oil as it moved under the ice. When the leads closed, the oil was trapped into tools on top of the ice at the floe edges. The further spread of the oil onto the surface of the ice floes appears to have been caused by the wind.     

引入拖曳系数参数化的海冰自由漂流模拟研究

Many interesting characteristics of sea ice drift depend on the atmospheric drag coefficient(C_a) and oceanic drag coefficient(C_w).Parameterizations of drag coefficients rather than constant values provide us a way to look insight into the dependence of these characteristics on sea ice conditions.In the present study,the parameterized ice drag coefficients are included into a free-drift sea ice dynamic model,and the wind factor α and the deflection angle θ between sea ice drift and wind velocity as well as the ratio of C_a to C_w are studied to investigate their dependence on the impact factors such as local drag coefficients,floe and ridge geometry.The results reveal that in an idealized steady ocean,C_a/C_w increases obviously with the increasing ice concentration for small ice floes in the marginal ice zone,while it remains at a steady level(0.2-0.25) for large floes in the central ice zone.The wind factor α increases rapidly at first and approaches a steady level of 0.018 when A is greater than 20%.And the deflection angle θ drops rapidly from an initial value of approximate 80° and decreases slowly as A is greater than20%without a steady level like α.The values of these parameters agree well with the previously reported observations in Arctic.The ridging intensity is an important parameter to determine the dominant contribution of the ratio of skin friction drag coefficient(C_s' /C_s) and the ratio of ridge form drag coefficient(C_r'/C_r) to the value of C_a/C_w,α,and θ,because of the dominance of ridge form drag for large ridging intensity and skin friction for small ridging intensity among the total drag forces.Parameterization of sea ice drag coefficients has the potential to be embedded into ice dynamic models to better account for the variability of sea ice in the transient Arctic Ocean.     

Sea ice volume and age: Sensitivity to physical parameterizations and thickness resolution in the CICE sea ice model

New dynamics parameterizations in Version 5 of the Los Alamos Sea Ice Model, CICE, feature an anisotropic rheology and variable drag coefficients. This study investigates their effect on Arctic sea ice volume and age simulations, along with the effects of several pre-existing model options: a parameter that represents the mean cumulative area of ice participating in ridging, the resolution of the ice thickness distribution, and the resolution of the vertical temperature and salinity profiles.By increasing shear stress between floes, the anisotropic rheology slows the ice motion, producing a thicker, older ice pack. The inclusion of variable drag coefficients, which depend on modeled roughness elements such as deformed ice and melt pond edges, leads to thinner ice and a more realistic simulation of sea ice age. Several feedback processes act to enhance differences among the runs. Notably, if less open water is produced mechanically through ice deformational processes, the simulated ice thins relative to runs with more mechanically produced open water. Thermodynamic processes can have opposing effects on ice age and volume; for instance, growth of new ice increases the volume while decreasing the age of the pack. Therefore, age data provides additional information useful for differentiating among process parameterization effects and sensitivities to other model parameters.Resolution of thicker ice types is crucial for proper modeling of sea ice volume, because the volume of ice in the thicker ice categories determines the total ice volume. Model thickness categories tend to focus resolution for thinner ice; this paper demonstrates that 5 ice thickness categories are not enough to accurately resolve the ice thickness distribution for simulations of ice volume.     

Influence of dynamic factors on ice rafting

Theoretical analysis and numerical simulation have been applied to researching the effect of wave, current and wind on the ice rafting. The analytical results show that wave is the key factor for ice fracture. The process of ice rafting is simulated using the discrete element method. The simulated results show that the distortion or flexure breakage happens firstly, then followed by ice rafting. The comprehensively analysed results show that the ice is easily rafted on the simultaneity action of wave, current and wind, and the action of wave can't be ignored. Moreover, the rafted ice length increases with the enhancive force of ocean current and wind.