海冰重叠动力条件分析

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

Aerial survey of a salt marsh: Ice rafting to the lower intertidal zone

An aerial photographic survey was conducted using an ultralight aircraft to assess both the extent and the directionality of ice rafted salt marsh pear in a northern New England salt marsh. To address the extent of ice rafting, a map was constructed from aerial photographs to determine the number, area and percentage of the tidal flat covered by accumulated ice rafted peat islands. Seven percent of the tidal flat or 117·4 m² of the salt marsh embayment was covered by 202 transported peat islands. To address the direction of ice rafting, an ice tagging experiment tested the hypothesis that salt marsh peat is transported from the high to the low intertidal. Eighty-eight percent of the recovered ice rafted peat islands were transported to lower intertidal elevations. Thus, the ice rafting process is an important physical disturbance in this north temperate salt marsh and has the potential to impact the salt marsh peat community.     

海冰动力过程的改进离散元模型及在渤海的应用

海冰的断裂、重叠和堆积等离散分布特性广泛地存在于极区和副极区的不同海域,并对海冰的生消、运移过程有着重要影响。针对海冰在不同尺度下的离散分布特点,发展海冰动力过程的离散元方法有助于完善海冰数值模式,提高海冰数值模拟的计算精度。为此,本文针对海冰生消运移过程中的非连续分布和形变特性,发展了适用于海冰动力过程的改进离散元模型(MDEM)。不同于传统离散元方法,该模型将海冰离散为具有一定厚度、尺寸和密集度的圆盘单元。海冰单元设为诸多浮冰块的集合体,其在运移和相互接触碰撞过程中,依照质量守恒发生单元尺寸、密集度和厚度的相应变化。基于海冰离散性和流变性的特点,该模型采用黏弹性接触本构模型计算单元间的作用力,并依据Mohr-Coulomb准则计算海冰法向作用下的塑性变形及切向摩擦力。为验证该模型的可靠性,本文对海冰在规则水域内的运移和堆积过程进行了分析,离散元计算结果与解析值相一致;此外,对旋转风场下海冰漂移规律的模拟进一步验证了本文方法的精确性。在此基础上,对渤海辽东湾的海冰动力过程进行了48h数值分析,计算结果与卫星遥感资料和油气作业区的海冰现场监测数据吻合良好。在下一步工作中将考虑海冰离散元模拟中的热力因素影响,发展具有冻结、断裂效应的海冰离散元模型,更精确地模拟海冰动力-热力耦合作用下的生消和运移过程。     

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

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.     

渤海海域冰水共存期的波能流密度推算

为给寒区海域的波浪能估算提供科学依据,提出一种合理推算冰水共存海域波浪条件及波能流密度的方法,该方法将海冰模型与水动力学模型耦合模拟得到的冰浓度以线性修正函数的方式纳入波浪模型的海面摩阻风速方程中,并基于MCT (model coupling toolkit)耦合器将海冰模型、水动力学模型与波浪模型进行实时耦合。基于该方法模拟了渤海冬季寒潮大风期间的海冰以及波能流密度的演化。模拟结果表明,在2012年2月5~8日寒潮大风期间,结冰区域占到渤海总面积的1/3,约有76%的渤海海域的平均波能流密度受海冰影响减小,其中辽东湾近岸的波能流密度平均受冰影响最多减小了100%,而渤海湾和莱州湾近岸受冰影响最多分别减小了60%和50%。即使是无冰覆盖的老铁山水道,其波能流密度的最大值也受冰影响减少了14%。耦合模拟可以更为准确地对渤海冬季的波能流密度分布进行评估,为波浪能发电厂选址提供依据。     

冰和直立结构相互作用的静冰力可靠性分析

我国渤海每年约有3个月的结冰期,冰与海洋结构的相互作用远较波浪为甚,因此有必要进行冰与海洋结构相互作用的可靠性分析.首先,采用热力学的方法根据每年实测的环境参数对海冰平整冰厚进行数值模拟计算,然后根据计算结果对不同重现期的极值冰厚采用极值统计法给出了极值分布,并与其他文献作比较.其次,以应变速率为海冰强度的主要影响参数,给出了海冰强度的概率分布.最后,以Schwarz公式为静冰力计算公式,结合得到的冰厚和冰强度的概率分布,采用映射变换法进行不同重现期静冰力的可靠度分析计算,并用蒙特卡罗重点抽样法进行模拟.     

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

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

西北航道夏季气旋登陆期间海冰对波浪的影响

西北航道是指从北大西洋经加拿大北极群岛进入北冰洋,再进入太平洋的航道,是连接大西洋和太平洋的捷径。为了探讨西北航道通航期极端天气条件下强风及海冰对波浪场的影响机制,建立并验证了考虑海冰影响下的西北航道风浪演化模型,并以2012年8月北极气旋登陆期间为例探讨西北航道通航期波浪特性及波能流密度的时空演化及其对风和海冰的响应。研究结果表明,北极夏季海冰大多分布于西北航道以北海域,而风向大部分集中在SSW(南偏西22.5°)至SW(南偏西45°),西北航道海冰的存在并不会引起有效风区的明显减少,也不会引起无冰海域波能流的明显减小(不超过5%)。但是,当风向变为北向风时,无冰海域波能流减小幅度最多高达62%。最后,综合海冰和波浪要素的时空分布,提出了极端天气条件下西北航道通航期的最佳适航路线,为西北航道的夏季安全通航提供了科学依据。     

Observations and modeling of the ice-ocean conditions in the coastal Chukchi and Beaufort Seas

The Chukchi and Beaufort Seas include several important hydrological features： inflow of the Pacific water, Alaska coast current （ ACC ）, the seasonal to perennial sea ice cover, and landfast ice ＇along the Alaskan coast. The dynamics of this coupled ice-ocean system is important for both regional scale oceanography and large-scale global climate change research. A mumber of moorings were deployed in the area by JAMSTEC since 1992, and the data revealed highly variable characteristics of the hydrological environment. A regional high-resolution coupled ice-ocean model of the Chukchi and Beaufort Seas was established to simulate the ice-ocean environment and unique seasonal landfast ice in the coastal Beaufort Sea. The model results reproduced the Beaufort gyre and the ACC. The depthaveraged annual mean ocean currents along the Beaufort Sea coast and shelf hreak compared well with data from four moored ADCPs, but the simulated velocity had smaller standard deviations, which indicate small-scale eddies were frequent in the region. The model resuits captured the sea,real variations of sea ice area as compared with remote sensing data, and the simulated sea ice velocity showed an ahnost stationary area along the Beaufort Sea coast that was similar to the observed landfast ice extent. It is the combined effects of the weak oceanic current near the coast, a prevailing wind with an onshore component, the opposite direction of the ocean current, and the blocking hy the coastline that make the Beaufort Sea coastal areas prone to the formation of landfast ice.     

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.     

渤海潮汐对冰作用的数值模拟

为研究渤海潮对海冰运动的作用,将渤海海冰模式与Ｂｌｕｍｂｅｒｇ的ＥＣＯＭ－ｓｉ模式通过动力作用联结,构成一种海冰一海洋动力耦合模式。应用上述耦合模式对渤海冬季潮流对冰的动力作用进行模拟研究,得出了渤海海冰运动的一些显著特征。     

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

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

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年,白令海峡和楚科奇海域处于低压和高压系统的交界处,盛行偏北风,海冰从北极东部往西部输运,加拿大海盆的多年海冰因离岸运动而辐散,向楚科奇海域的海冰输运增加,受太平洋入流暖水影响,移入此区域的海冰加速融化,从而加剧海冰的减少。     

北极冰海耦合模式对两种不同大气再分析资料响应的分析

本文中我们比较了Climate Forecast System Reanalysis(CFSR)高分辨率的再分析数据集和低分辨率的Japanese 25-year Reanalysis Project(JRA25)再分析数据集在向下短波辐射、向下长波辐射、10m风场、近地面气温、降水、湿度上的不同,发现二者差异最大的为降水数据,其次为向下短波辐射数据、向下长波辐射数据。用这两个数据集驱动同一冰海耦合模式,CFSR强迫的海冰、北冰洋中层水和加拿大海盆温盐结构与实测相比有很大差距,等密度面上的地转流速在加拿大海盆和欧亚海盆比JRA25强迫的结果高20%,同时等密度面的深度偏深、位温偏高,在弗拉姆海峡的流通量也比海洋再分析数据Simple Ocean Data Assimilation(SODA)偏多。CFSR的向下辐射数据更加接近实测,采用此数据的敏感性实验模拟结果与实测符合的更好。对于海冰的模拟,云量起着至关重要的作用,降水带来的淡水通量通过影响大西洋入流水携带的热量进而影响到冰区。此外,CFSR过量的降水也是二者对于北冰洋温盐结构、弗拉姆海峡流通量以及地转流强度模拟产生偏差的主要原因。尽管风场的分辨率不同,在海盆尺度上对于海冰和海水温盐结构的影响并不大。     

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.     

锥角对锥体结构抗冰性能影响的离散元分析

在寒区海洋工程中,锥体海洋结构的尺寸参数对其破冰性能具有重要影响。采用具有黏结-破碎功能的离散元方法模拟海冰与锥体结构的相互作用过程中海冰的破坏模式及冰载荷分析。该离散元方法的计算参数通过与渤海油气平台的实测数据对比进行了可靠性验证。在此基础上,当设计海域的潮差固定时,对不同锥角下风电单桩锥体结构的冰载荷和海冰破坏模式进行了离散元分析。计算结果表明:锥体冰载荷随锥角的增大而增大,海冰的平均断裂长度则随锥角的增大而减小;当锥角小于70°时海冰的破碎模式主要为弯曲破坏,而当锥角大于70°时海冰破坏模式则主要为挤压破碎。以上研究结果可为冰区海上风电单桩结构的抗冰锥设计提供参考。     

Experimental Investigation on Breakup of Ice Floe on Waves

Experimental investigation of for floe breakup on waves by use of non-refrigerated breakablematerials has been carried out.The incident wave heights for breakup of ice floe with different lengths,andthe influence of rigidity on reflection and transmission coeffcients are discussed.The experimental resultsshow that the ratio of the ice floe length,L_i,to the wavelength,L,is one of the significant factor affectingthe minimum wave height to cause fracture of ice floe,and another key factor is the ratio of the ice floethickness,h_i and L_i.     

BCC_CSM对北极海冰的模拟:CMIP5和CMIP6历史试验比较

本文利用北京气候中心气候系统模式(BCC_CSM)在最近两个耦合模式比较计划(CMIP5和CMIP6)的历史试验模拟结果,对北极海冰范围和冰厚的模拟性能进行了比较,结果表明:(1) CMIP6改善了CMIP5模拟海冰范围季节变化过大的问题,总体上更接近观测结果;(2)两个CMIP试验阶段中BCC_CSM模拟的海冰厚度都偏小,但CMIP6试验对夏季海冰厚度过薄问题有所改进。通过对影响海冰生消过程的冰面和冰底热收支的分析,我们探讨了上述模拟偏差以及CMIP6模拟结果改善的成因。分析表明,8?9月海洋热通量、向下短波辐射和反照率对模拟结果的误差影响较大,CMIP6试验在这些方面有较大改善;而12月至翌年2月,CMIP5模拟的北极海冰范围偏大主要是海洋热通量偏低所导致,CMIP6模拟的海洋热通量较CMIP5大,但北大西洋表层海流的改善才是巴芬湾附近海冰外缘线位置改善的主要原因。CMIP试验模拟的夏季海冰厚度偏薄主要是因为6?8月海洋热通量和冰面热收支都偏大,而CMIP6试验模拟的夏季海冰厚度有所改善主要是由于海洋热通量和净短波辐射的改善。海冰模拟结果的改善与CMIP6海冰模块和大气模块参数化的改进有直接和间接的关系,通过改变短波辐射、冰面反照率和海洋热通量,使BCC_CSM模式对北极海冰的模拟性能也得到有效提高。     

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

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

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.