基于中国北极考察布放浮标观测和数值模拟的海冰和积雪厚度变化分析

Sea ice and the snow pack on top of it were investigated using Chinese National Arctic Research Expedition(CHINARE) buoy data.Two polar hydrometeorological drifters,known as Zeno? ice stations,were deployed during CHINARE 2003.A new type of high-resolution Snow and Ice Mass Balance Arrays,known as SIMBA buoys,were deployed during CHINARE 2014.Data from those buoys were applied to investigate the thickness of sea ice and snow in the CHINARE domain.A simple approach was applied to estimate the average snow thickness on the basis of Zeno~ temperature data.Snow and ice thicknesses were also derived from vertical temperature profile data based on the SIMBA buoys.A one-dimensional snow and ice thermodynamic model(HIGHTSI) was applied to calculate the snow and ice thickness along the buoy drift trajectories.The model forcing was based on forecasts and analyses of the European Centre for Medium-Range Weather Forecasts(ECMWF).The Zeno~ buoys drifted in a confined area during 2003–2004.The snow thickness modelled applying HIGHTSI was consistent with results based on Zeno~ buoy data.The SIMBA buoys drifted from 81.1°N,157.4°W to 73.5°N,134.9°W in 15 months during2014–2015.The total ice thickness increased from an initial August 2014 value of 1.97 m to a maximum value of2.45 m before the onset of snow melt in May 2015;the last observation was approximately 1 m in late November2015.The ice thickness based on HIGHTSI agreed with SIMBA measurements,in particular when the seasonal variation of oceanic heat flux was taken into account,but the modelled snow thickness differed from the observed one.Sea ice thickness derived from SIMBA data was reasonably good in cold conditions,but challenges remain in both snow and ice thickness in summer.     

渤海模型冰运动特征的现场观测研究

Sea ice drift is mainly controlled by ocean currents, local wind, and internal ice stress. Information on sea ice motion, especially in situ synchronous observation of an ice velocity, a current velocity, and a wind speed, is of great significance to identify ice drift characteristics. A sea ice substitute, the so-called "modelled ice", which is made by polypropylene material with a density similar to Bohai Sea ice, is used to complete a free drift experiment in the open sea. The trajectories of isolated modelled ice, currents and wind in the Bohai Sea during non-frozen and frozen periods are obtained. The results show that the currents play a major role while the wind plays a minor role in the free drift of isolated modelled ice when the wind is mild in the Bohai Sea. The modelled ice drift is significantly affected by the ocean current and wind based on the ice–current–wind relationship established by a multiple linear regression. The modelled ice velocity calculated by the multiple linear regression is close to that of the in situ observation, the magnitude of the error between the calculated and observed ice velocities is less than12.05%, and the velocity direction error is less than 6.21°. Thus, the ice velocity can be estimated based on the observed current velocity and wind speed when the in situ observed ice velocity is missing. And the modelled ice of same thickness with a smaller density is more sensitive to the current velocity and the wind speed changes. In addition, the modelled ice drift characteristics are shown to be close to those of the real sea ice, which indicates that the modelled ice can be used as a good substitute of real ice for in situ observation of the free ice drift in the open sea, which helps solve time availability, safety and logistics problems related to in situ observation on real ice.     

HY-2卫星扫描微波辐射计数据反演北极海冰漂移速度

本文基于最大互相关法,利用海洋二号（HY-2）卫星扫描微波辐射计37 GHz通道多时相垂直极化亮温数据,获取了北极海冰漂移速度。采用2012年和2013年国际北极浮标计划海冰现场观测数据,对利用微波辐射计亮温资料反演的冬季北极海冰漂移速度进行了定量验证,结果表明:流速和流向均方根误差分别为1.12 cm/s和16.37°,从一定程度上说明了HY-2卫星扫描微波辐射计亮温数据反演海冰漂移速度的可行性。此外,使用美国国防气象卫星F-17搭载的专用微波成像仪91 GHz通道垂直极化亮温,采用高斯拉普拉斯滤波方法进行处理,结合最大互相关法反演的海冰漂移速度,优于法国海洋开发研究院海冰漂移速度产品。     

2016年南极中山站固定冰冰厚观测分析

极区海冰是全球气候系统的重要组成部分，南极的固定冰普遍存在于其沿海地区，中山站周边固定冰一般在11月中下旬达到最厚。海冰厚度是海冰的重要参数之一，2016年在南极中山站附近3个站点(S1、S2、S3站点)共布放了4套温度链浮标，包括1套SIMBA (Snow and Ice Mass Balance Array)温度链浮标和3套太原理工大学温度链浮标(TY温度链浮标)，SIMBA温度链浮标每天观测4次，TY温度链浮标每小时观测1次。利用浮标观测的温度剖面以及海冰和海水间不同介质温度差异计算得到海冰厚度。在S3站点，同时布放了SIMBA温度链浮标和TY温度链浮标。温度链浮标计算冰厚和人工钻孔观测冰厚比较结果显示，S1站点TY温度链浮标计算的海冰厚度平均误差和均方根误差分别为3.3 cm和14.7 cm，S2站点和S3站点分别为6.6 cm、6.9 cm以及4.0 cm、4.8 cm。S3站点的SIMBA温度链浮标计算冰厚和人工观测冰厚的平均误差和均方根误差为8.2 cm和9.7 cm。因而S3站点TY温度链浮标计算的海冰厚度更接近人工观测的结果。进一步对Stefan定律海冰生长模型进行对比，模型计算得到的海冰生长率为0.1～0.8 cm/d，生长率快于TY温度链浮标的结果，且受积雪影响明显。相比于卫星遥感反演冰厚的误差和观测时段的限制以及有限的人工观测，2种温度链浮标未来对于中山站附近海冰的长期监测均有重要的应用价值。     

1979-2012年北极海冰运动学特征初步分析

利用美国冰雪数据中心(NSIDC)发布的海冰速度和范围数据,本文分析了1979—2012年间北极海冰的运动学特征,以及北极海冰运动与分布范围演变之间的关系。结合欧洲中期天气预报中心(ECMWF)发布的2007和2012年高分辨率的气压场、风场数据,探讨了北极风场和气压场与海冰运动、辐散辐合和海冰面积的关系。结果表明,在1979-2012年间北极海冰平均运动速度呈显著增强的趋势,冬季海冰平均运动速度增加趋势明显强于夏季;北极、波弗特-楚科奇海域和弗拉姆海峡的冬、夏季海冰平均运动速度的增加率分别为2.1%/a和1.7%/a、2.0%/a和1.6%/a以及4.9%/a和2.2%/a。1979-2012年北极海冰平均运动速度和范围的相关性为-0.77,二者存在显著的负相关关系。北极冬季和夏季风场的长期变化趋势与海冰平均运动速度的变化趋势一致,冬季和夏季的相关系数分别为0.50和0.48。风场和气压场对海冰的运动、辐散及重新分布发挥着重要作用。2007年夏季,第234~273天波弗特海域一直被高压系统控制,波弗特涡旋加强,使得波弗特海域海冰聚集在北极中央区;顺时针的风场促使海冰向格陵兰岛和加拿大北极群岛以北聚合。2012年,白令海峡和楚科奇海域处于低压和高压系统的交界处,盛行偏北风,海冰从北极东部往西部输运,加拿大海盆的多年海冰因离岸运动而辐散,向楚科奇海域的海冰输运增加,受太平洋入流暖水影响,移入此区域的海冰加速融化,从而加剧海冰的减少。     

Observed and modelled sea‐ice drift response to wind forcing in the northern Baltic Sea

The wind dependence of sea‐ice motion was studied on the basis of ice velocity and wind observations, and weather model output. The study area was a transition zone between open water and the ice‐covered ocean in the northern Baltic Sea. In the centre of the basin the sea‐ice motion was highly wind‐dependent and the linear relationship between the wind and the drift velocities explained 80% of the drift's variance. On the contrary, the wind‐drift dependence was low near the coast. The wind‐drift coherence was significant over a broader frequency range in the central part of the basin than for the coastal drift. The ice motion was simulated by a numerical model forced with five types of wind stress and with two types of current data, and the outcome was compared with the observed buoy drift. The wind and the wind‐induced surface current were the main factors driving the ice in the basin's centre, while internal ice stresses were of importance in the shear zone near the fast ice edge. The best wind forcing was achieved by applying a method dependent on atmospheric stability and ice conditions. The average air–ice drag coefficient was 1.4×10⁻³ with the standard deviation of 0.2×10⁻³. The improvement brought about by using an accurate wind stress was comparable with that achieved by raising the model grid resolution from 18 km to 5 km.     

Mechanisms controlling the temporal variation of the sea ice edge in the Southern Ocean

This paper examines the mechanism controlling the short time-scale variation of sea ice cover over the Southern Ocean. Sea ice concentration and ice velocity datasets derived from images of the Defense Meteorological Satellite Program (DMSP) Special Sensor Microwave Imager (SSM/I) are employed to reveal this mechanism. The contribution of both dynamic and thermodynamic processes to the change in ice edge location is examined by comparing the meridional velocity of ice edge displacement and sea ice drift. In the winter expansion phase, the thermodynamic process of new ice production off the ice edge plays an important role in daily advances of ice cover, whereas daily retreats are mostly due to southward ice drift. On the other hand, both advance and retreat of ice edges in the spring contraction phase are mostly caused by the dynamic process of the ice drift. Based on the above mechanism and the linear relation between the degree of ice production at the ice edge and northward wind speed, the seasonal advance of ice cover can be roughly reproduced using the meridional velocity of ice drift at the ice edge.     

海水物质平衡浮标对北冰洋中心区海冰温度与物质平衡的观测

利用中国第3次北极科学考察所布放海冰物质平衡浮标(Ice Mass-Balance buoy,IMB)的观测数据,分析了北冰洋中心区多年冰2008年8月-2009年7月温度与物质平衡的变化特征.结果表明,冰温廓线呈现明显的季节变化,秋季降温过程从海冰表面开始向冰体内部传播.海冰底部的生长/消融率受海水温度控制,随水温的...     

Numerical simulation for dynamical processes of sea ice

ＮｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｆｏｒｄｙｎａｍｉｃａｌｐｒｏｃｅｓｓｅｓｏｆｓｅａｉｃｅＷｕＨｕｉｄｉｎｇ，ＢａｉＳｈａｎ，ＺｈａｎｇＺｈａｎｈａｉａｎｄＬｉＧｕｏｑｉｎｇ（ＲｅｃｅｉｖｅｄＭａｙ１６，１９９６；ａｃｃｅｐｔｅｄＪａｎｕａｒｙ...     

渤海冰漂移对海面风场、潮流场的响应

在对海冰漂移动力学分析基础上,利用MODIS资料,采用MCC方法获取渤海大范围冰覆盖的海域冰速场,并利用NCEP风速资料和潮流资料进行回归分析,得到渤海冰漂移速度与风速和流速的关系.利用MODIS和NOAA/AVHRR资料获取的渤海冰速资料的综合分析显示:渤海海冰运动,除受盛行风控制外,还受到复杂的海岸地形、流和冰内应力的共同作用,所得到的大范围海冰运动规律和多年历史观测资料分析结果基本一致,并清楚地显示了冰边缘带海冰运动的特征,弥补了局地、单站海冰观测的局限性.     

夏季加拿大海盆海冰边缘区声体积后向散射强度研究

声信号的体积后向散射强度是声传播过程中一个关键的参数。海冰边缘区的声体积后向散射强度研究对深入认识北极声场环境有着十分重要的意义。本文利用中国第六次北极科学考察获取的数据资料研究了海冰边缘区声体积后向散射强度特性。结果表明:加拿大海盆海冰边缘区是声体积后向散射强度的明显过渡区。无冰海面(海冰密集度小于15%)海洋深层水的声体积后向散射强度明显大于密集海冰区域的海水(海冰密集度大于50%)。讨论了声体积后向散射强度与海冰融化之间的关系,造成融冰区声体积后向散射强度增大的原因是水下悬浮泥沙、浮游生物等悬浮物质增加。根据海冰密集海域的海水后向散射强度弱的特点,对北极下放式声学多普勒测流仪(LADCP)观测的设置提出建议。     

大直径单桩风机基础冰荷载模型试验研究

针对渤海某海域以单桩结构为支撑的海上风电系统,对大直径直立桩风机基础进行了一系列静冰载荷模型试验研究。首先,针对目标海域海冰调研结果确定多个冰速工况,对3 MW及4 MW装机功率对应的两种不同直径的单桩基础开展静冰载荷模型试验;随后在试验现象及试验结果的分析中重点关注了冰排在大直径结构前的破坏模式及破坏特点;最终,通过对比模型试验冰力极值试验结果与多规范冰载荷计算结果,确定大直径直立桩静冰载荷计算规范的合理选择,并根据试验结果对直立桩静冰载荷计算方法进行了经验参数修正。得出的相关结论可为渤海海域大直径单桩式风机基础冰载荷的工程估算提供参考。     

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

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

北极各海域海冰覆盖范围的变化特征

Sea ice in the Arctic has been reducing rapidly in the past half century due to global warming.This study analyzes the variations of sea ice extent in the entire Arctic Ocean and its sub regions.The results indicate that sea ice extent reduction during 1979–2013 is most significant in summer,following by that in autumn,winter and spring.In years with rich sea ice,sea ice extent anomaly with seasonal cycle removed changes with a period of 4–6 years.The year of 2003–2006 is the ice-rich period with diverse regional difference in this century.In years with poor sea ice,sea ice margin retreats further north in the Arctic.Sea ice in the Fram Strait changes in an opposite way to that in the entire Arctic.Sea ice coverage index in melting-freezing period is an critical indicator for sea ice changes,which shows an coincident change in the Arctic and sub regions.Since 2002,Region C2 in north of the Pacific sector contributes most to sea ice changes in the central Aarctic,followed by C1 and C3.Sea ice changes in different regions show three relationships.The correlation coefficient between sea ice coverage index of the Chukchi Sea and that of the East Siberian Sea is high,suggesting good consistency of ice variation.In the Atlantic sector,sea ice changes are coincided with each other between the Kara Sea and the Barents Sea as a result of warm inflow into the Kara Sea from the Barents Sea.Sea ice changes in the central Arctic are affected by surrounding seas.     

渤海海冰特征厚度分析

通过海冰生消机理和数值试验,讨论了渤海海冰特征厚度的存在条件;对不同厚度的海冰表面温度、冰面热量收支、冰面下热传导和太阳辐射透射量进行了对比分析,分析了渤海海冰向特征冰厚的动态演化过程;在不同气温、风速、相对湿度和海洋热通量等气象和水文条件下,对渤海特征冰厚进行了计算;讨论了海冰生消的动态平衡过程,分析了1997/1998年冬季辽东湾JZ20-2海域实测冰厚与特征冰厚的相互关系。对渤海特征冰厚分析将有助于渤海海冰数值模拟工作的开展和对不同重现期设计冰厚的推算。     

海冰动力学过程的数值模拟

讨论了海冰动力学性质并阐述决定海冰漂移的动量平衡,冰脊和水道形成及确定冰应力与形变、强度之间关系的海冰流变学.提出了模拟海冰动力学过程的数值模式,模式中冰厚分布由开阔水、平整冰和堆积冰3种要素表示.在这3要素的预报方程中引入形变函数,采用一种参数化方法模拟冰脊和水道.为了表示冰内相互作用,将海冰作为一种非线性粘性可压缩物质,采用粘-塑性本构关系.本文还概述和讨论了模式中所采用的数值方法,应用此模式模拟了渤海、波罗的海的波的尼亚湾和拉布拉多海的冰漂移.渤海冰漂移模拟结果明显地显示出潮周期变化,还模拟了渤海的冰脊和水道,进行了海冰流变学参数的敏感性试验.并将此冰模式与大气模式和边界层模式联接,给出渤海海冰预报结果.     

夏季北极密集冰区范围确定及其时空变化研究

研究夏季北极密集冰区的范围变化是了解北极海冰融化过程的重要手段。密集冰区与海冰边缘区之间没有明确的分界线, 海冰密集度在两者之间平滑过渡, 确定密集冰区范围就需确定一个密集度阈值。文中依据分辨率为6.25 km的AMSR-E遥感数据, 发现不同密集度阈值所围范围在密集冰区边缘处的减小存在由快变慢的过程, 同时与周围格点的密集度差异变化在该处最为显著, 对这两个特征进行统计分析, 获得的阈值同为89%, 具有明确的物理意义和合理性。以此为基础, 运用腐蚀算法剔除海冰边缘区, 同时结合连通域法排除小范围密集冰的影响, 进而确定密集冰区的范围。结果表明, 2002-2006年密集冰区覆盖范围较大, 年际变化较小, 2007年以后明显减小, 2010年与2011年相继出现最小值, 其中2011年的范围最小值仅为2006年的64%。密集冰区范围的变化不同于海冰覆盖范围, 是具有独立特性的海冰变化参数, 反映出高密集度海冰区域的变化特征。海冰的融化与海冰边缘区的变化是导致密集冰区范围发生变化的两个主要因素, 受动力学因素的影响, 海冰边缘区发生伸展或收缩, 发生密集冰区与海冰边缘区互相转化。本文提出了一种研究北极海冰变化的新思路, 密集冰区覆盖范围的减小表明北极中央区域高密集度海冰正持续减少。     

Numerical sea ice prediction in China

ＮｕｍｅｒｉｃａｌｓｅａｉｃｅｐｒｅｄｉｃｔｉｏｎｉｎＣｈｉｎａＷｕＨｕｉｄｉｎｇ,ＢａｉＳｈａｎ,ＺｈａｎｇＺｈａｎｈａｉ１（ＲｅｃｅｉｖｅｄＳｅｐｔｅｍｂｅｒ１２,１９９６;ａｃｃｅｐｔｅｄＪｕｎｅ５,１９９７）Ａｂｓｔｒａｃｔ──Ａｄｙｎａｍｉ...     

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

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.     

基于质点-网格模式的海冰厚度变化过程数值模拟

根据渤海冰情,在海冰动力学和热力学研究基础上,应用一种质点-网格海冰模式于渤海海域.该模式采用了质点-网格法,有效地避免了传统模式的数值扩散问题.该模式采用了冰厚分布函数,用多种类型冰代替用于渤海业务预报的平整冰、堆积冰和开阔水3-level模式.进行理想场的数值试验,模拟冰厚变化动力过程.还使用该模式和业务预报模式对于实际渤海冰情进行了不同个例的预报试验,发现该模式在提高冰外缘线预报精度方面有一定的优势.