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

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

Numerical sea ice prediction in China

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

渤海海冰漂移数值研究

本文利用海冰－海洋动力耦合模式对渤海典型天气形势下海冰动力过程作了数值模拟，海冰模式建立在动量和质量守恒基础上，忽略了海冰变化的热力过程。海冰厚度被划分为三类：堆积冰、平整冰和开阔冰，冰的形变由一个厚度的重新分布约束条件确定。海洋模式是一个二维风暴潮模式．同时考虑了潮汐的作用。风场资料来自于沿岸气象观测站每日四次的风观测，计算网格是十分之一经纬度。主要的分析和模拟是在大气和海洋共同作用下海冰的漂流特点。观测结果比较表明，数值模拟结果基本上反映了该海域流冰的漂移性质，同时也可为短期冰情预报提供有益的参考。     

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

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

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

在对海冰漂移动力学分析基础上,利用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.     

基于Mohr-Coulomb准则的黏弹-塑性海冰动力学本构模型

针对中小尺度下海冰动力作用过程中的漂移、重叠和堆积特征,在连续介质力学基础上建立了一个黏弹塑性海冰动力学本构模型.该模型主要包括四部分,即海冰在屈服前的Kelvin Vogit黏弹模型、海冰塑性屈服的Mohr Coulomb准则、屈服后相关联的正交流动法则和影响海冰强度的静水压力.采用黏弹塑性本构模型对规则海域内的海冰堆积过程进行了数值试验,其结果与传统黏塑性模型、改进的黏塑性模型和经典冰坝理论的解析解相符,可合理地模拟海冰堆积的动力过程.为进一步检验该黏弹塑性本构模型的适用性,对渤海海冰的动力过程进行了48h的数值模拟,计算结果与海冰卫星遥感现场观测资料一致.通过与传统的黏塑性和改进的黏塑性模型结果的对比分析,进一步验证了黏弹塑性模型在中小尺度海冰动力学模拟中的可靠性.在以上海冰堆积的数值试验和渤海海冰动力学模拟中均采用了光滑质点流体动力学(SPH)方法.     

渤海冰-海洋耦合模式——I．模式和参数研究

在国内外冰-海洋耦合模式研究基础上，根据渤海水文、气象和冰情特点，以国家海洋环境预报中心的渤海海冰预报模式和POM海洋模式为基础，开发了一个冰-海洋耦合模式．在该耦合模式中，冰和海洋之间的动量和热量交换是双向的，冰厚和冰密集度的变化不仅由冰表面和冰底的热收支决定，还由开阔水的表面热收支决定．侧重阐述了耦合模式的动力和热力学过程的耦合，并对模式中一些热力参数进行了讨论．     

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

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

质点-网格海冰模式在渤海的数值预报试验

在海冰动力学和热力学的研究基础上,开发应用于渤海的质点-网格海冰模式.该模式采用质点-网格法,有效地避免了传统海冰模式的数值扩散问题.该模式采用冰厚分布函数取代传统渤海业务海冰预报模式中平整冰、堆积冰和开阔水的三类海冰分型.应用质点-网格海冰模式对2003年～2004年冬季渤海冰情进行了业务化逐日数值预报试验,并对预报结果进行统计检验和分析比较.     

A coupled ice-ocean model for the Bohai Sea I. Study on model and parameter

1 IntroductionAs an important compartment of global cli-mate system, the development and variation ofsea ice directly affect the global atmosphere andocean environment, especially the deep watersformation and circulation. Therefore, it be-comes one of the keys of the global climatestudy to explore the coupling physical processof sea ice to atmosphere and ocean in the cli-mate system and to set up the global coupledairiceocean system. The studies on the in-teraction of iceocean and its coupled…     

A coupled ice-ocean model for the Bohai Sea Ⅰ. Study on model and parameter

According to the earlier international studies on the coupled ice-ocean model and the hydrology, meteorology, and ice features in the Bohai Sea, a coupled ice-ocean model is developed based on the National Marine Environment Forecast Center's (NMEFC) numerical forecasting ice model of the Bohai Sea and the Princeton ocean model (POM).In the coupled model, the transfer of momentum and heat between ocean and ice is two-way, and the change of ice thickness and concentration depends on heat budget not only at the surface and bottom of ice, but also at the surface of open water between ices. The dynamic and thermodynamic coupling process is expatiated emphatically. Some thermodynamic parameters are discussed as well.     

A coupled ice-ocean model for the Bohai Sea Ⅰ. Study on model and parameter

According to the earlier international studies on the coupled ice-ocean model and the hydrology, meteorology, and ice features in the Bohai Sea, a coupled ice-ocean model is developed based on the National Marine Environment Forecast Center‘ s (NMEFC) numerical forecasting ice model of the Bohai Sea and the Princeton ocean model (POM). In the coupled model, the transfer of momentum and heat between ocean and ice is two-way, and the change of ice thickness and concentration depends on heat budget not only at the surface and bottom of ice, but also at the surface of open water between ices. The dynamic and thermodynamic coupling process is expatiated emphatically. Some thermodynamic parameters are discussed as well.     

渤海冰-海洋耦合模式Ⅰ.模式和参数研究

在国内外冰-海洋耦合模式研究基础上,根据渤海水文、气象和冰情特点,以国家海洋环境预报中心的渤海海冰预报模式和POM海洋模式为基础,开发了一个冰-海洋耦合模式.在该耦合模式中,冰和海洋之间的动量和热量交换是双向的,冰厚和冰密集度的变化不仅由冰表面和冰底的热收支决定,还由开阔水的表面热收支决定.侧重阐述了耦合模式的动力和热力学过程的耦合,并对模式中一些热力参数进行了讨论.     

Monitoring regional sea ice of China's Bohai Sea by SSM/Iscattering indexes

The Specific Sensor Microwave/Imager (SSM/I) data over China's Bohai Sea from 1993 to 1996 has been studied. Some empirical algorithms of data calibration/validation for sea ice, such as the method of the Atmospheric Environment Science/York University temperature brightness index (AES/YORK TBI) and ocean scattering index (OSI), are not well applicable to detection of regional sea ice such as the midlatitude Bohai Sea area. This paper presents a method of the scattering index (SI) and polarization index (PI) to identify regional Bohai Sea ice. Numerical results of vector radiative transfer for a model of a layer of sea ice over ocean simulate the relationship between SIs/PIs and sea ice depth and other variations. It is employed to categorize the graph of SIs/PIs for detection of sea ice     

2005～2006年冬季渤海海冰数值预报

本文概述2005-2006年冬季渤海冰情。介绍该冬季业务化海冰数值预报工作和应用PIC模式开展海冰数值后报试验,其业务化海冰数值预报的客观检验结果满足业务预报精度指标。     

渤海海冰现场监测的数字图像技术及其应用

在渤海冰区油气开发中,海冰给平台结构、油气运输和施工作业带来很大影响.油气作业区海冰参数精确、连续、实时的现场监测对分析油气开发的可靠性、检验海冰数值模式、校正海冰卫星遥感数据具有重要意义.针对渤海油气作业区的海冰运动和分布特性,通过数字图像技术对海冰的厚度、运动速度和密集度三个参数的提取进行了算法开发和软件研制.在2...     

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

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.     

星载SSM/Ⅰ微波遥感渤海海冰的辐射特征分析

根据我国渤海区域星载微波遥感ＤＭＳＰＳＳＭ／Ｉ１９９４～１９９６年的观察数据的研究,指出现有的星载微波遥感海冰的数据标定和验证ＴＢＩ和ＯＳＩ算法是不适用的。用２２和８５ＧＨｚ的散射指数ＳＩ和１９ＧＨｚ的极化指数ＰＩ,讨论了渤海区域海冰的主要辐射特征,用矢量辐射传输理论数值模拟了ＳＩ、ＰＩ与海冰覆盖及其厚度变化的定量关系,用ＳＩ－ＰＩ图规划了渤海区域性海冰的识别方法,给出了渤海海冰的空间分布及其辐射特征的时间序列变化。     

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

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