Wave-ice dynamical interaction: a numerical model and its application

In this paper, an ice floe inner stress caused by the wave-induced bending moment is derived to estimate the stress failure of ice floe. The strain and stress failures are combined to establish a wave-induced ice yield scheme. We added ice stress and strain failure module in the Finite-Volume Community Ocean Model (FVCOM), which already includes module of ice-induced wave attenuation. Thus a fully coupled wave-ice dynamical interaction model is established based on the ice and wave modules of FVCOM. This model is applied to reproduce the ice and wave fields of the breakup events observed during the second Sea Ice Physics and Ecosystem Experiment (SIPEX-2) voyage. The simulation results show that by adopting the combined wave-induced ice yield scheme, the model can successfully predict the ice breakup events, which the strain failure model is unable to predict. By comparing the critical significant wave height deduced from strain and stress failure schemes, it is concluded that the ice breakup is caused by the strain failure when wave periods are shorter than a threshold value, while the stress failure is the main reason for the ice breakup when wave periods are longer than the threshold value. Neglecting either of these two ice-break inducement mechanisms could overestimate the ice floe size, and thus underestimate the velocity of the ice lateral melt and increase the error of simulation of polar ice extent.     

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

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

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…     

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.     

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…     

Discrete elementmodeling of ice loads on ship hulls in broken ice fields

Ice loads on a ship hull affect the safety of the hull structure and the ship maneuvering performance in ice-covered regions. A discrete element method （DEM） is used to simulate the interaction between drifting ice floes and a moving ship. The pancake ice floes are modelled with three-dimensional （3-D） dilated disk elements considering the buoyancy, drag force and additional mass induced by the current. The ship hull is modelled with 3D disks with overlaps. Ice loads on the ship hull are determined through the contact detection between ice floe element and ship hull element and the contact force calculation. The influences of different ice conditions （current velocities and directions, ice thicknesses, concentrations and ice floe sizes） and ship speeds are also examined on the dynamic ice force. The simulated results are compared qualitatively well with the existing field data and other numerical results. This work can be helpful in the shil3 structure design and the navigation securitv in ice-covered fields.     

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.     

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

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

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

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

浮冰界面融化速率参数化方案的实验室研究

融冰季节时天然浮冰表面、底面和侧向融化共存,三者融化速率是底面大于侧向,侧向大于表面。而且浮冰尺寸越小,侧向速率占比越高。为了解决将小尺度浮冰块尺度指标计入融化参数化方案,在低温环境实验室无辐射、无流速、控制气温和水温条件下对天然海冰和人工冻结淡水冰的圆盘试样,开展了不同初始水温和不同初始直径的圆盘试样融化过程实验,获得了圆盘试样直径、厚度和质量融化过程。依据这些实验数据,构建试样直径厚度比这一新指标,通过物理分析和数学统计手段,建立了海冰和淡水冰试样表面、底面融化速率同温度梯度,侧向融化速率和温差以及直径厚度比的关系式。这些关系式能够应用于天然直径100 m范围内浮冰的融化参数化方案,期望能解决北冰洋海冰和入海口近岸淡水冰夏季融化季节能量和质量平衡数值模拟的需求。     

北极薄冰-煎饼冰区域中不同重力波驱动形成海流的敏感性分析

A coupled ocean-ice-wave model is used to study ice-edge jet and eddy genesis during surface gravity wave dissipation in a frazil-pancake ice zone. With observational data from the Beaufort Sea, possible wave dissipation processes are evaluated using sensitivity experiments. As wave energy dissipated, energy was transferred into ice floe through radiation stress. Later, energy was in turn transferred into current through ocean-ice interfacial stress. Since most of the wave energy is dissipated at the ice edge, ice-edge jets, which contained strong horizontal shear, appeared both in the ice zone and the ocean. Meanwhile, the wave propagation direction determines the velocity partition in the along-ice-edge and cross-ice-edge directions, which in turn determines the strength of the along-ice-edge jet and cross-ice-edge velocity. The momentum applied in the along-ice-edge(cross-ice-edge)direction increased(decreased) with larger incident angle, which is favorable condition for producing stronger mesoscale eddies, vice versa. The dissipation rate increases(decreases) with larger(smaller) wavenumber, which enhances(reduces) the jet strength and the strength of the mesoscale eddy. The strong along-ice-edge jet may extend to a deep layer(> 200 m). If the water depth is too shallow(e.g., 80 m), the jet may be largely dampened by bottom drag, and no visible mesoscale eddies are found. The results suggest that the bathymetry and incident wavenumber(magnitude and propagation direction) are important for wave-driven current and mesoscale eddy genesis.     

2004～2005年冬季南堡-曹妃甸海域海冰监测报告

南堡-曹妃甸海域是渤海湾冰情严重海域之一。本文根据海上平台、陆地测点、航空监测、船舶调查及卫星遥感等有关海冰监测资料及相邻海域有关海冰监测资料,给出了该海域本年度冰期、冰型、冰厚、结冰范围、浮冰最大外缘线及海冰时空分布变化等状况,综合分析得出,本年度该海域的冰情状况接近于正常年份。     

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.     

海冰重叠动力条件分析

由于重叠冰厚度较大,给冰区海洋运输以及海洋生产作业等带来较大影响,为了得到科学客观的重叠冰厚度,需要对重叠冰的形成条件进行深入研究。利用理论分析和数值模拟方法,分析了波浪、风和流对渤海海冰重叠的影响。通过理论计算分析发现,波浪是导致海冰断裂的主要因素;采用离散元方法对海冰重叠过程进行了数值模拟,模拟结果表明海冰在外界动力作用下,海冰初始发生屈曲变形断裂,随后产生重叠。综合分析表明,海冰在同时具备波浪的垂向波动作用以及风或流水平作用力下,冰层比较容易发生重叠,波浪对海冰重叠的作用不可忽略,海冰的重叠长度取决于风和流作用的大小,其作用力越大,重叠长度越大。     

Evolution of the thickness distribution of ice in the sea of Japan

Simple mathematical models of the thermodynamics of the ice-cover thickness and the thermodynamics of the area of an individual ice floe are proposed. The equations of the models allow for an explicit consideration of the spatial boundedness of the seawater region containing ice covers. A kinetic model of the evolution of the distributions of ice area and ice thickness is formulated on the basis of the gas-dynamic theory. Integration of the equations of this model over the surface areas of individual floes gives the thickness distribution of ice areas. Several special cases are studied analytically. The adequacy of the models is assessed. The results of simulations are presented.     

Comparisons of passive microwave remote sensing sea ice concentrations with ship-based visual observations during the CHINARE Arctic summer cruises of 2010–2018

In order to apply satellite data to guiding navigation in the Arctic more effectively, the sea ice concentrations(SIC)derived from passive microwave(PM) products were compared with ship-based visual observations(OBS)collected during the Chinese National Arctic Research Expeditions(CHINARE). A total of 3 667 observations were collected in the Arctic summers of 2010, 2012, 2014, 2016, and 2018. PM SIC were derived from the NASA-Team(NT), Bootstrap(BT) and Climate Data Record(CDR) algorithms based on the SSMIS sensor, as well as the BT,enhanced NASA-Team(NT2) and ARTIST Sea Ice(ASI) algorithms based on AMSR-E/AMSR-2 sensors. The daily arithmetic average of PM SIC values and the daily weighted average of OBS SIC values were used for the comparisons. The correlation coefficients(CC), biases and root mean square deviations(RMSD) between PM SIC and OBS SIC were compared in terms of the overall trend, and under mild/normal/severe ice conditions. Using the OBS data, the influences of floe size and ice thickness on the SIC retrieval of different PM products were evaluated by calculating the daily weighted average of floe size code and ice thickness. Our results show that CC values range from 0.89(AMSR-E/AMSR-2 NT2) to 0.95(SSMIS NT), biases range from-3.96%(SSMIS NT) to 12.05%(AMSR-E/AMSR-2 NT2), and RMSD values range from 10.81%(SSMIS NT) to 20.15%(AMSR-E/AMSR-2 NT2). Floe size has a significant influence on the SIC retrievals of the PM products, and most of the PM products tend to underestimate SIC under smaller floe size conditions and overestimate SIC under larger floe size conditions. Ice thickness thicker than 30 cm does not have a significant influence on the SIC retrieval of PM products. Overall, the best(worst) agreement occurs between OBS SIC and SSMIS NT(AMSR-E/AMSR-2 NT2) SIC in the Arctic summer.     

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的一致性最差。     

基于油气平台的海冰雷达监测图像处理及冰速测量

对海冰的运动规律进行精确、连续和长周期的实时监测有助于海冰热力学和动力学的研究,也可保障冰区生产活动的安全进行。针对辽东湾海冰的运动特点和工程需求,在JZ20-2油气平台上建立了海冰雷达监测系统。采用数字图像处理技术对海冰雷达监测图像进行了分析和软件开发,可对海冰密集度、速度和冰块面积等海冰参数进行提取。采用该海冰雷达监测系统和数字图像处理软件,在2011-2012年冬季对该海域的海冰运动规律进行了全冰期的连续监测,在此基础上重点对海冰速度的雷达图像监测结果进行了分析,讨论了海冰速度场分布以及连续48 h的变化过程。以上结果为海冰的生消运移规律研究和油气作业区的海冰管理工作提供了可靠的现场监测数据。对海冰雷达现场监测及数字图像处理中的问题及改进方法进行了讨论。     

一种冰-海洋模式的热力耦合方案

冰与海洋的热力耦合对冰与海洋环流的模拟有极其重要的影响,是冰-海洋相互作用的一个重要方面.对其精确确定需要详细考虑冰-海洋界面附近的湍流过程,这在长时期的模拟特别是气候模拟中,常受到技术条件的限制.过去的研究常常假设冰下海洋混合层的温度为冰点,特别是在单纯冰模式的模拟中,但考虑海冰漂移和冰点变化的效应时,这一假设是不精确的.因此,弱化冰下海洋混合层温度为冰点的约束,不考虑详细的冰-海洋界面和海洋混合层的湍流过程,根据冰-海洋耦合系统的能量收支关系,设计了一个简化的冰-海洋热力耦合方案.对该方案引起的海洋混合层适应、热力结构和海冰发展的影响进行了分析,并将其用于全球冰海洋耦合模式的数值试验,结果表明,在大气热力强迫下该耦合方案既可使冰区混合层海水温度向冰点适应,又使冰边缘带海水温度与冰点保持明显差异,能够较好地反映冰-海洋热力相互作用.利用该耦合方案构造的全球冰-海洋耦合模式模拟的海冰范围及季节变化与实际观测非常接近.     

渤海和北极海冰组构及晶体结构特征分析

渤海一年冰的组构和晶体结构表现出的冰的非连续生长速率使得冰层底部具有气泡含量不同的许多分层;静水生长的灰冰具有典型的粒状冰和柱状冰;融化-重冻结冰和重叠冰具有特定的晶体结构.位于北冰洋72°24.037'N,153°33.994'W,长2.2m的冰心样为三年冰,并且发现其内有一种新的晶体类型,定义为碎屑凝聚冰.由位于北冰洋74°58.614'N,160°31.830'W,长4.86m的冰心样的晶体结构分析发现它过去为冰脊,后经融化-冻结改造形成冰丘的证据.