• 中文核心期刊
  • 中国科技核心期刊
  • ISSN 1007-6336
  • CN 21-1168/X

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于海冰热力学原理的营口海域冰厚演变规律及工程应用

马玉贤 谭红建 于福祥 陈雪 袁帅 许宁 史文奇

引用本文:
Citation:

基于海冰热力学原理的营口海域冰厚演变规律及工程应用

    作者简介: 马玉贤(1989-),男,山西大同人,博士研究生,主要从事海冰物理力学性质研究,E-mail:yxma@nmemc.org.cn;
    通讯作者: 许 宁(1983-),女,河北秦皇岛人,研究员,主要从事海冰物理力学性质研究,E-mail:nxu@nmemc.org.cn
  • 基金项目: 国家重点研发计划项目(2017YFA0604901);大连市高层次人才创新项目(2019RJ07);辽宁省民生科技计划(2021JH2/10300004)
  • 中图分类号: P753

The evolution law and engineering application of ice thickness in Yingkou sea area based on the thermodynamic of sea ice

  • 摘要: 本文依据营口气象站1951年以来的气象观测数据,结合热力学模式HIGHTSI,重现了1951-2018年冬季(67 a)海冰的生长与融化过程。对计算结果进行分析,多年平均冰厚在2月初较为对称,多年最大冰厚约为0.24 cm。自1951年以来最大海冰厚度与平均海冰厚度的减小速率分别为0.255 cm/a和0.076 cm/a,海冰厚度随时间的演变规律为:海冰厚度越大,发生概率最大值对应的年份越早,即大于40 cm海冰厚度的最大概率发生于20世纪50年代,小于10 cm海冰厚度的最大概率发生于20世纪90年代。1951-2018年,每10年的代际海冰厚度分布概率分析结果表明,20世纪70年代为持续的常冰年冰情,20世纪90年代为持续的轻冰年,各海冰厚度的概率呈连续减小趋势;其余代际年份的冰情呈现轻冰年与重冰年的波动变化,各海冰厚度的概率呈减小-增大-减小的规律。本文定量描述了气候变化对营口冰情的影响,为海洋工程设计的海冰参数选取奠定了数据基础。
  • 图 1  2018年雷达观测数据与计算结果(黑色竖线是沿岸观测的平整冰冰厚)

    Figure 1.  Seasonal evolution of modeled and observed ice thickness for 2018 ice season( The vertical bar is the standard deviation of observed ice thickness)

    图 2  营口气象站冬季平均气温趋势

    Figure 2.  The average temperature change trend of Yingkou

    图 3  1951-2018年海冰厚度计算结果;黑色线为平均海冰厚度[30]

    Figure 3.  The ice thicknesses for 1951-2018. The thick black line is the modeled seasonal average ice thicknesses

    图 4  海冰厚度演变过程图

    Figure 4.  Graph of the evolution of ice thickness

    图 5  代际海冰厚度发生概率

    Figure 5.  Probability of decade ice thickness

    表 1  HIGHTSI模型的主要组成部分和本文参数选取[30]

    Table 1.  Major components of HIGHTSI model and parameters used for this study

    方程冰表面热量和质量平衡冰-雪温度冰底热量和质量平衡
    参数 入射短波辐射
    表面反照率
    I(z) :穿透短波辐射
    入射与反射长波辐射
    感热通量与潜热通量
    表面透射率
    Fc:传递到表面的热通量
    Fm:表面融化的热通量
    Tsfc:表面温度
    T:温度
    t:时间
    z:垂直坐标
    ρ:密度
    c:比热容
    k:热传导系数
    q(z,t):吸收的太阳辐射
    si:分别代表雪与冰
    hi:海冰厚度
    Li:冻结潜热
    Fw:冰底热通量
    本文参数设置 ki= 1.5 m−1
    αs,I:参数化方案(Briegleb, et al., 2004)
    0.97
    k0 = 2.03 W m−1 K−1
    ρi = 910 kg m−3
    ci = 2093 J kg−1 K−1
    Fw= 2 W/m2
    Tf= −1.4 ℃
    Li= 0.33 × 106 J kg−1
    输入 风速;气温;相对湿度;云量系数;短波辐射(参数化)
    模型设置 步长:1 h;垂直分辨率:20层;初始海冰厚度:5 cm
    输出 海冰厚度;冰温
    下载: 导出CSV
  • [1] HUANG Y, SHI Q Z, SONG A. Model test study of the interaction between ice and a compliant vertical narrow structure[J]. Cold Regions Science and Technology, 2007, 49(2): 151-160. doi: 10.1016/j.coldregions.2007.01.004
    [2] YANG G J. Bohai sea ice conditions[J]. Journal of Cold Regions Engineering, 2000, 14(2): 54-67. doi: 10.1061/(ASCE)0887-381X(2000)14:2(54)
    [3] LI Z J, KANG J C, PU Y B. Fabric and crystal characteristics of Bohai and Arctic sea ice[J]. Chinese Journal of Polar Science, 2002, 13(1): 23-28.
    [4] 马玉贤, 关 湃, 许 宁, 等. 辽东湾海洋工程可靠性设计的海冰参数选取[J]. 海洋工程, 2019, 37(3): 136-142.
    [5] LI N, MAKI T, HAYAKAWA S, et al. Relationship between sea ice thickness and temperature in Bohai Sea of China[J]. Journal of the Faculty of Agriculture, Kyushu University, 2005, 50(1): 165-173. doi: 10.5109/4634
    [6] MAYKUT G A, UNTERSTEINER N. Some results from a time-dependent thermodynamic model of sea ice[J]. Journal of Geophysical Research, 1971, 76(6): 1550-1575. doi: 10.1029/JC076i006p01550
    [7] 程 斌. 一维海冰热力模式的守恒型差分格式和数值模拟[J]. 海洋通报, 1996, 15(4): 8-16.
    [8] CHENG B, LAUNIANEN J, VIHMA T. Modelling of superimposed ice formation and sub-surface melting in the Baltic Sea[J]. Geophysica, 2003, 39(1/2): 31-50.
    [9] CHENG B, VIHMA T, PIRAZZINI R, et al. Modelling of superimposed ice formation during the spring snowmelt period in the Baltic Sea[J]. Annals of Glaciology, 2006, 44: 139-146. doi: 10.3189/172756406781811277
    [10] CHENG B, ZHANG Z H, VIHMA T, et al. Model experiments on snow and ice thermodynamics in the Arctic Ocean with CHINARE 2003 data[J]. Journal of Geophysical Research, 2008, 113(C9): C09020.
    [11] CHENG B, VIHMA T, RONTU L, et al. Evolution of snow and ice temperature, thickness and energy balance in Lake Orajärvi, northern Finland[J]. Tellus A:Dynamic Meteorology and Oceanography, 2014, 66(1): 21564. doi: 10.3402/tellusa.v66.21564
    [12] CHENG B, VIHMA T, PALO T, et al. Observation and modelling of snow and sea ice mass balance and its sensitivity to atmospheric forcing during spring and summer 2007 in the Central Arctic[J]. Advances in Polar Science, 2021, 32(4): 312-326.
    [13] WANG C X, CHENG B, WANG K G, et al. Modelling snow ice and superimposed ice on landfast sea ice in Kongsfjorden, Svalbard[J]. Polar Research, 2015, 34: 20828. doi: 10.3402/polar.v34.20828
    [14] LAUNIAINEN J, CHENG B. Modelling of ice thermodynamics in natural water bodies[J]. Cold Regions Science and Technology, 1998, 27(3): 153-178. doi: 10.1016/S0165-232X(98)00009-3
    [15] YAN Y, UOTILA P, HUANG K Y, et al. Variability of sea ice area in the Bohai Sea from 1958 to 2015[J]. Science of the Total Environment, 2020, 709: 136164. doi: 10.1016/j.scitotenv.2019.136164
    [16] MÄKYNEN M, KARVONEN J, CHENG B, et al. Operational service for mapping the Baltic sea landfast ice properties[J]. Remote Sensing, 2020, 12(24): 4032. doi: 10.3390/rs12244032
    [17] KARVONEN J, CHENG B, VIHMA T, et al. A method for sea ice thickness and concentration analysis based on SAR data and a thermodynamic model[J]. The Cryosphere, 2012, 6(6): 1507-1526. doi: 10.5194/tc-6-1507-2012
    [18] 吴辉碇, 白 珊. 海冰动力学过程的数值模拟[J]. 海洋学报, 1998, 20(02): 1-13.
    [19] HUANG Y. Model test study of the nonsimultaneous failure of ice before wide conical structures[J]. Cold Regions Science and Technology, 2010, 63(3): 87-96. doi: 10.1016/j.coldregions.2010.06.004
    [20] SODHI D S. Crushing failure during ice——structure interaction[J]. Engineering Fracture Mechanics, 2001, 68(17/18): 1889-1921.
    [21] SCHULSON E M. Brittle failure of ice[J]. Engineering Fracture Mechanics, 2001, 68(17/18): 1839-1887.
    [22] 季顺迎, 岳前进. 渤海海冰特征厚度分析[J]. 海洋学报, 2000, 22(6): 117-123.
    [23] 季顺迎, 岳前进, 毕祥军. 辽东湾JZ20–2海域海冰参数的概率分布[J]. 海洋工程, 2002, 20(3): 39-43. doi: 10.3969/j.issn.1005-9865.2002.03.007
    [24] 吴辉碇. 渤海海冰设计作业条件[M]. 北京: 海洋出版社, 2001.
    [25] 段忠东, 杨冀东, 欧进萍, 等. 渤海及黄海北部年极值冰厚概率分布的统计分析[J]. 黄渤海海洋, 2000, 18(3): 34-40.
    [26] 刘永青, 李辉辉, 史文奇, 等. 单层平整冰厚设计条件推算方法的比较——以辽东湾JZ20-2海域为例[J]. 海洋环境科学, 2017, 36(6): 941-946.
    [27] GONG D Y, KIM S J, HO C H. Arctic Oscillation and ice severity in the Bohai Sea, East Asia[J]. International Journal of Climatology, 2007, 27(10): 1287-1302. doi: 10.1002/joc.1470
    [28] OUYANG L X, HUI F M, ZHU L X, et al. The spatiotemporal patterns of sea ice in the Bohai Sea during the winter seasons of 2000–2016[J]. International Journal of Digital Earth, 2019, 12(5): 893-909.
    [29] YAN Y, SHAO D D, GU W, et al. Multidecadal anomalies of Bohai Sea ice cover and potential climate driving factors during 1988–2015[J]. Environmental Research Letters, 2017, 12(9): 094014. doi: 10.1088/1748-9326/aa8116
    [30] MA Y X, CHENG B, XU N, et al. Long-term ice conditions in Yingkou, a coastal region northeast of the Bohai Sea, between 1951/1952 and 2017/2018: Modeling and Observations[J]. Remote Sensing, 2022, 14(1): 182. doi: 10.3390/rs14010182
  • [1] 韩蕾赵芮刘昭贾宝林姜祖林冯少姝梁玉波 . 大田软海绵酸和鳍藻毒素时间分辨荧光免疫层析试纸条的研制与应用. 海洋环境科学, 2022, 41(5): 783-790. doi: 10.12111/j.mes.2021-x-0194
    [2] 付东洋栾虹刘大召张莹丁又专王文芳李薛蒋城飞 . 珠江口冬春季悬浮泥沙浓度遥感反演模式分析. 海洋环境科学, 2016, 35(4): 600-604. doi: 10.13634/j.cnki.mes.2016.04.019
    [3] 贺琪孔令斌赵丹枫黄冬梅杜艳玲 . 基于限定motif关联规则的海洋热浪事件模式发现研究. 海洋环境科学, 2022, 41(6): 937-946. doi: 10.12111/j.mes.2022-x-0093
    [4] 周凤霞袁柳婷陆旋陈法锦朱庆梅孟亚飞 . 湛江湾表层沉积物不同形态磷的变化特征. 海洋环境科学, 2022, 41(4): 509-518. doi: 10.12111/j.mes.20210095
    [5] 孙鑫任景玲刘家琦房瑞雪蒋增杰 . 桑沟湾不同形态锰的分布、季节变化及影响因素. 海洋环境科学, 2023, 42(4): 511-522, 533. doi: 10.12111/j.mes.2022-x-0194
    [6] 王心哲谭琳琳范剑超 . 2010-2020年马来西亚马登红树林遥感监测与景观格局变化分析. 海洋环境科学, 2022, 41(4): 603-609. doi: 10.12111/j.mes.2021-x-0039
    [7] 徐冠球谭晓璇屠建波石海明何荣刘洋 . 天津大神堂牡蛎礁保护区海域海水水质变化趋势分析与评价. 海洋环境科学, 2022, 41(4): 554-562. doi: 10.12111/j.mes2021-x-0051
    [8] 于兵刘子洲翟方国顾艳镇吴文凡 . 2020年夏季威海瑜泰海洋牧场底层海水溶解氧的日变化研究. 海洋环境科学, 2022, 41(4): 563-571. doi: 10.12111/j.mes.2021-x-0070
    [9] 龚照辉池源扎边元冰瑜高建华 . 江苏省海岸带近30年人类干扰指数时空变化特征. 海洋环境科学, 2022, 41(6): 821-829. doi: 10.12111/j.mes.2021-x-0297
    [10] 罗新正孟宪伟 . 基于相对海平面变化的广西茅尾海红树林面积稳定性评估. 海洋环境科学, 2022, 41(6): 881-887. doi: 10.12111/j.mes.2022-x-0014
    [11] 姚慧平樊彦国禹定峰周燕杨雷高皜卞晓东赵春燕 . 2003-2020年黄海水体透明度时空变化特征及影响因素研究. 海洋环境科学, 2023, 42(2): 262-270. doi: 10.12111/j.mes.2022-x-0252
    [12] 方草宁志铭韦芬夏荣林杨斌 . 基于LOICZ收支模型的涠洲岛近海营养盐浓度及结构变化的关键影响因素研究. 海洋环境科学, 2022, 41(5): 689-696. doi: 10.12111/j.mes.20210040
    [13] 陈智杰曾志陈楚汉汤军健 . 滨海电厂温排水影响下的温升平面分布和垂向变化研究——以罗源湾可门电厂为例. 海洋环境科学, 2023, 42(1): 111-121. doi: 10.12111/j.mes.2022-x-0081
  • 加载中
图(5)表(1)
计量
  • 文章访问数:  1675
  • HTML全文浏览量:  877
  • PDF下载量:  92
出版历程
  • 收稿日期:  2022-01-03
  • 录用日期:  2022-04-24
  • 网络出版日期:  2022-06-06
  • 刊出日期:  2022-12-20

基于海冰热力学原理的营口海域冰厚演变规律及工程应用

    作者简介:马玉贤(1989-),男,山西大同人,博士研究生,主要从事海冰物理力学性质研究,E-mail:yxma@nmemc.org.cn
    通讯作者: 许 宁(1983-),女,河北秦皇岛人,研究员,主要从事海冰物理力学性质研究,E-mail:nxu@nmemc.org.cn
  • 1. 国家海洋环境监测中心, 辽宁 大连 116023
  • 2. 自然资源部第三海洋研究所, 福建 厦门 361005
  • 3. 大连测控技术研究所, 辽宁 大连 116013
基金项目: 国家重点研发计划项目(2017YFA0604901);大连市高层次人才创新项目(2019RJ07);辽宁省民生科技计划(2021JH2/10300004)

摘要: 本文依据营口气象站1951年以来的气象观测数据,结合热力学模式HIGHTSI,重现了1951-2018年冬季(67 a)海冰的生长与融化过程。对计算结果进行分析,多年平均冰厚在2月初较为对称,多年最大冰厚约为0.24 cm。自1951年以来最大海冰厚度与平均海冰厚度的减小速率分别为0.255 cm/a和0.076 cm/a,海冰厚度随时间的演变规律为:海冰厚度越大,发生概率最大值对应的年份越早,即大于40 cm海冰厚度的最大概率发生于20世纪50年代,小于10 cm海冰厚度的最大概率发生于20世纪90年代。1951-2018年,每10年的代际海冰厚度分布概率分析结果表明,20世纪70年代为持续的常冰年冰情,20世纪90年代为持续的轻冰年,各海冰厚度的概率呈连续减小趋势;其余代际年份的冰情呈现轻冰年与重冰年的波动变化,各海冰厚度的概率呈减小-增大-减小的规律。本文定量描述了气候变化对营口冰情的影响,为海洋工程设计的海冰参数选取奠定了数据基础。

English Abstract

  • 渤海海域每年冬季都有海冰生成,季节性海冰对涉海经济活动有重要影响。针对海冰的影响,科研人员开展了大量调查和研究工作[1]。随着温室气体的快速排放,气候变化引起的海冰冰情变化也逐渐被重视。渤海海域表征冰情严重程度的主要参数包括浮冰范围与海冰厚度[2-3]。受连续数据获取方式的限制,当前关于海冰浮冰范围的研究较多,以海冰厚度为对象的研究较少。浮冰范围仅在宏观尺度上对冰情严重程度进行评价,无法通过内部状态评估海冰的生消过程。为明确气候变化对冰情演变的影响,有必要剔除动力过程对冰情的影响,着重分析热力要素变化引起的冰情演变。

    当前,以鲅鱼圈雷达海冰观测(近20年)为代表的业务化海冰连续观测[4],由于数据序列长度不足等问题,不能有力支撑当前气候变化引起的海冰厚度演变研究工作的开展。海冰数值模型可以解决观测数据的不足问题。现有研究认为,海冰单层冰冰厚与天气条件有很强的联系[5],依据可靠的气象观测数据,可以得到有价值的长期海冰厚度数据。Maykut和Untersteiner[6](MU71模型)研发了首个海冰厚度方向的热力学模型,成功模拟了海冰随季节的生消过程,其模式框架一直沿用至今。程斌[7]利用积分插值法构造出守恒型差分格式,数值结果表明,积分插值法构造的守恒型差分格式计算稳定。Cheng和Launiainen等[8-14]在已有海冰热力学模式基础上发展了一维高分辨率雪冰热力学模式(HIGHTSI),给出较为完整的气-冰耦合参数化方案,通过积雪层、冰层的高分辨率进行计算,提高了模式计算的精度。20世纪80年代末90年代初,研究者建立了渤海海冰数值预报系统,渤海海冰物理和模拟论文陆续发表[15-18]。中国和芬兰在1990年冬季联合开展鲅鱼圈近岸海域气象、海冰冰温等参数的现场调查,获得了8 d逐小时连续数据。Cheng等[14]基于这些数据,进行了HIGHTSI模式冰温的检验。

    海冰载荷一直是有冰海域海洋工程结构关注的重点。一般而言,引起海洋工程结构构件失效的原因有两种:由冰力极值引起的强度失效与由交变冰力引起的疲劳失效 [19-21]。海冰厚度是极值冰力取值、交变冰力周期与幅值的控制性因素。因此,海冰厚度的合理取值是评估海洋工程结构、海冰灾害的重要内容。为确定冰情数据,研究者基于多年的实测气象资料对海冰厚度进行了大量的数值推算,季顺迎等[22-23]基于辽东湾 JZ20-2海域1996-2000年海冰定点观测数据与数值模拟结果,拟合出该海域平整冰冰厚的概率分布,发现海冰厚度服从对数正态分布。吴辉碇[24]认为不同冰情等级的发生概率服从正态分布,并给出不同重现期各冰区的设计冰厚。段忠东等[25]采用平稳二项随机过程模型和组合概率方法,统计分析了年最大海冰厚度的概率分布及参数,并推算了若干年重现期的极值海冰厚度。 然而,有关资料显示,与20世纪相比,近20年我国海域冰情偏轻[26-29],因此有必要对气候变化背景下海冰厚度的概率密度变化进行分析。

    本文选取辽东湾冰情最重的营口海域作为研究对象,使用1951年以来的气象观测数据,结合一维高分辨率热力学模式(HIGHTSI),重现1951-2018年营口海域平整冰的冰厚数据。以上述计算结果为基础,对气候变化背景下海冰厚度演变规律及其对工程的影响进行评估。

    • 渤海海冰主要分布于辽东湾、渤海湾、莱州湾,其中辽东湾冰情最严重。由于渤海湾与莱州湾冬季有冰天数较少,气候变化引起的冰情变化不明显。因此,本文选取辽东湾冰情最重的营口海域为研究对象。营口位于辽东湾顶部东侧,附近海域平均水深为5 m,海水盐度约为28,冰情为每年12月初至次年3月底。

    • (1)气象数据

      气象数据源于中国气象网《中国地面气候资料日值数据集》(http://data.cma.cn/)。根据各省上报的全国地面月报信息化文件和《全国地面气候资料(1961-1990)统计方法》《地面气象观测规范》,整编统计得到数据集。数据集为中国756个基本、基准地面气象观测站及自动站1951年到现在的日值数据集,要素包括:平均气温、日最高气温、日最低气温、平均相对湿度、最小相对湿度、平均风速、最大风速、最大风速的风向、极大风速、极大风速的风向、日照时数。本文选择距离鲅鱼圈雷达站最近的营口气象站数据进行分析。

      (2)鲅鱼圈海冰观测数据

      为研究渤海海冰的生消运移规律,国家海洋局在营口鲅鱼圈墩台山建立了海冰定点观测系统,通过雷达监测系统对海冰要素进行连续同步测量。海冰信息包括海冰厚度、类型和密集度等,其中,海冰厚度由海冰冰型反演而来。现场海冰观测主要针对雷达覆盖范围内的平整冰,海冰厚度代表能见范围内海冰的一般厚度(能见范围内分布最多的平整冰厚度)。海冰厚度观测结果虽然不是固定海冰的观测,但本文近似认为雷达覆盖范围内平整冰厚度与垂向热力学过程模拟结果具有一致性,观测数据可以为冰厚生长过程与平整冰冰厚最大值提供评价依据。本文以2018年冬季,每日8时、14时、20时海冰观测数据作为海冰厚度分析的基础资料。

    • 一维高分辨率雪/冰热力学数值模式(HIGHTSI)可模拟雪/冰生消过程,计算雪/冰的表面温度、内部温度和厚度变化。HIGHTSI模式已用于渤海、南极湖、北极海、波罗的海海冰和芬兰湖冰等的热力学数值模拟,现为芬兰气象局对波罗的海海冰的日常预报业务模式。模式采用的关键过程、模型参数、外强迫项及输出结果见表1。基于上述设置,本文对营口海域海冰进行热力学计算,时间为1951-2018年(67 a)每年12月1日至次年3月31日。

      方程冰表面热量和质量平衡冰-雪温度冰底热量和质量平衡
      参数 入射短波辐射
      表面反照率
      I(z) :穿透短波辐射
      入射与反射长波辐射
      感热通量与潜热通量
      表面透射率
      Fc:传递到表面的热通量
      Fm:表面融化的热通量
      Tsfc:表面温度
      T:温度
      t:时间
      z:垂直坐标
      ρ:密度
      c:比热容
      k:热传导系数
      q(z,t):吸收的太阳辐射
      si:分别代表雪与冰
      hi:海冰厚度
      Li:冻结潜热
      Fw:冰底热通量
      本文参数设置 ki= 1.5 m−1
      αs,I:参数化方案(Briegleb, et al., 2004)
      0.97
      k0 = 2.03 W m−1 K−1
      ρi = 910 kg m−3
      ci = 2093 J kg−1 K−1
      Fw= 2 W/m2
      Tf= −1.4 ℃
      Li= 0.33 × 106 J kg−1
      输入 风速;气温;相对湿度;云量系数;短波辐射(参数化)
      模型设置 步长:1 h;垂直分辨率:20层;初始海冰厚度:5 cm
      输出 海冰厚度;冰温

      表 1  HIGHTSI模型的主要组成部分和本文参数选取[30]

      Table 1.  Major components of HIGHTSI model and parameters used for this study

    • 本文以2018年冬季(常冰年)鲅鱼圈雷达观测站海冰厚度观测数据为基础,验证计算结果。鲅鱼圈附近海冰观测内容有:逐时雷达观测海冰厚度数据和每周1次现场沿岸观测。上述观测数据与本文计算结果比对如图1所示,可以发现,数值计算结果可以模拟初冰期与盛冰期的海冰情况,但对融冰期海冰的模拟存在偏差。这是因为,本文计算采用的是一维热力学模式,没有考虑海冰动力过程。现场观测显示,在初冰期与盛冰期,观测范围内分布大范围平整冰,且布满能见海面;在融冰期,离岸风导致平整冰快速向深水区输运(离开观测区域),此时观测到的海冰大多为重新成为浮冰的近岸固定冰,造成浮冰冰厚观测结果明显大于计算结果。可见,一维热力学模式可以模拟营口海域平整冰的生消过程。

      图  1  2018年雷达观测数据与计算结果(黑色竖线是沿岸观测的平整冰冰厚)

      Figure 1.  Seasonal evolution of modeled and observed ice thickness for 2018 ice season( The vertical bar is the standard deviation of observed ice thickness)

    • 作为海冰生消过程的重要条件,气温一直是海冰冰情发展的控制性热力条件。本文使用资料可查的营口气象站(鲅鱼圈附近)数据,对长序列的气温变化过程进行阐述。一般而言,辽东湾初冰日均在12月1日后,融冰日在3月31日前,所以营口气象站气象数据为12月1日至次年3月31日,计算结果如图2所示。从图2可以发现,气象站冬季气温呈明显增加趋势,年际升温趋势均为0.033 ℃/a(M-K趋势检验),升温趋势在1990年之后尤为明显。对冬季可能有海冰生成的月份进行月平均值趋势分析,12月与3月气温增加趋势相对较小,均为0.02 ℃/a;1月平均升温速率次之,约为0.034 ℃/a;2月升温趋势最明显,约为0.055 ℃/a。

      图  2  营口气象站冬季平均气温趋势

      Figure 2.  The average temperature change trend of Yingkou

    • 图3给出了67个冰季模型海冰厚度的季节演变及平均值,可以发现,季节海冰厚度演变表现出较大的年际变化,并且近几年海冰厚度与冰期长度呈明显减小趋势。所有年份海冰厚度的平均结果分布较为对称(图3黑色线),最大平均海冰厚度发生在2月初,平均24 cm。对每年计算结果的最大海冰厚度和平均海冰厚度进行Theil-Sen斜率分析,发现两者均呈减少趋势,减小速率分别为0.255 cm/a和0.076 cm/a。

      图  3  1951-2018年海冰厚度计算结果;黑色线为平均海冰厚度[30]

      Figure 3.  The ice thicknesses for 1951-2018. The thick black line is the modeled seasonal average ice thicknesses

    • 图3可以发现,营口海域海冰厚度年际变化波动较大,不利于海冰厚度演变规律的研究。为此,本文选择10 a的海冰厚度计算结果为样本,明确不同海冰厚度区间随时间的演变规律。以1955年度为例,使用1951-1960年冬季海冰厚度计算结果为样本库,分别计算[0~10 cm]、[10~20 cm]、[20~30 cm]、[30~40 cm]、[40~50 cm]、[50~60 cm]海冰厚度的发生概率。依据上述计算方法,获得1955-2013年冬季10 a滑动海冰厚度的发生概率。为直观展现海冰厚度的演变规律,本文对10~60 cm海冰厚度的分布概率进行插值,绘制1955-2014年期间10~60 cm海冰厚度的概率变化云图,结果见图4a;因为小于10 cm海冰厚度发生概率明显大于其他区间,所以本文单独绘制其发生概率的演变曲线,结果见图4b

      图  4  海冰厚度演变过程图

      Figure 4.  Graph of the evolution of ice thickness

      图4可以发现,20世纪70年代之前,大冰厚(大于35 cm视为极端海冰厚度情景)的发生概率逐渐减小,基本转移至20~30 cm;20世纪80年代,20 cm以上的海冰厚度发生概率逐渐减小,小于10 cm海冰的发生概率逐渐增大;20世纪90年代,气温连续维持高位,冰情较轻,小于10 cm海冰厚度的发生概率达到最大,其他海冰厚度发生概率维持最小;进入21世纪,气温波动变化较小,20~30 cm海冰厚度区间的发生概率略有增加,10 cm以下海冰发生的概率略有减小。若把20世纪90年代异常升温数据剔除掉,则20世纪60年代后海冰厚度呈连续减小趋势。

      对不同海冰厚度区间发生概率随时间的变化进行研究,可以发现,(1)当海冰厚度大于50 cm时,发生的概率逐渐减少,在20世纪90年代降为0,最大发生概率为20世纪50年代。(2)当海冰厚度为40~50 cm时,海冰厚度呈逐渐减小趋势(除2000年),最大发生概率为20世纪50年代。(3)当海冰厚度为30~40 cm时,发生概率呈先增大后减小的趋势,最大发生概率为20世纪60年代。(4)当海冰厚度处于20~30 cm时,发生概率呈先增大后减小再增大的趋势,最大发生概率为20世纪70年代,最小概率为21世纪前10年。(5)当海冰厚度处于10~20 cm时,发生概率呈先增大后减小的趋势,最大发生概率为20世纪70年代。(6)当海冰厚度小于10 cm时,除20世纪90年代明显增大外,小冰厚(小于10 cm)的发生概率随时间呈增大趋势,最大发生概率为1990年之后。

      自1950年以来,海冰厚度的演变规律为:冰厚越大,其发生概率的时间也越早。大于40 cm海冰的最大发生概率为20世纪50年代,30~40 cm海冰的最大发生概率为20世纪60年代,10~30 cm海冰的最大发生概率为20世纪70年代,小于10 cm海冰的最大发生概率为20世纪90年代。以上结果进一步证实了,随气温的升高,10年季大冰厚发生概率逐渐向小冰厚转移,并最终演化至小于10 cm。

    • 海洋工程结构使用累积损伤法计算结构寿命,离不开冰载荷概率密度函数的准确阐述。为明确气候变化背景下海冰概率密度函数的变化,本文分别对20世纪50年代、60年代、70年代、80年代、90年代,21世纪00年代、10年代,以10 cm为步长给出各自的分布概率,结果见图5。从图5可以发现,21世纪10年代的海冰厚度主要集中在30 cm以下。

      图  5  代际海冰厚度发生概率

      Figure 5.  Probability of decade ice thickness

      以海冰厚度为自变量对分布概率进行分析,发现20世纪50年代、60年代、80年代,21世纪00年代与10年代均表现为随海冰厚度增大,各厚度区间的概率呈减小-增大-减小的规律,也就是说,海冰厚度的分布概率表现出两个峰值(小于10 cm海冰厚度的发生概率与30 cm左右海冰厚度的分布概率);70年代、90年代期间冰情整体稳定,各海冰厚度区间的分布概率表现出连续减小的趋势。可见,一般情况下10年内冰情的严重程度呈现轻冰年与重冰年的波动变化,而70年代表现为持续的常冰年冰情、90年代表现为持续的轻冰年冰情。

      将本文计算结果与当前可查的两个海冰概率密度函数对比讨论。文献[23]依据1996-2000年(4 a)的观测与计算结果,给出锦州海域概率密度函数。这一概率密度函数的最大值在5 cm左右,且函数只有一个峰值(缺30 cm处的峰值);这一结果与本文计算结果的第一个峰值相当,与90年代海冰厚度分布概率结果类似;导致这一结果的原因是文献[23]所选取的海冰厚度样本库恰好为冰情较轻的时期,样本库缺乏冰情偏重年份的海冰厚度样本。文献[4]以鲅鱼圈雷达海冰观测站近20 a(2001-2018年)的实测海冰厚度数据进行统计分析,给出鲅鱼圈海域的海冰厚度概率密度函数;这一概率密度函数的最大值为30 cm左右,同样仅有一个峰值(缺5 cm处峰值),与本文计算结果的第二个峰值相当,与同时期模拟海冰厚度的概率分布结果类似;导致这一结果的原因是,鲅鱼圈海冰观测在每年小冰厚时期的观测数据存在缺失(观测基本从1月开始,此时海冰厚度已达10 cm以上),导致文献[4]同样仅有一个峰值。

    • (1)所有年份平均最大海冰厚度发生于2月初,平均最大海冰厚度约为24 cm。对每年计算结果的最大海冰厚度和平均海冰厚度进行Theil-Sen斜率分析,发现二者均呈减少趋势,减小速率分别为0.255 cm/a和0.076 cm/a。

      (2)自1950年以来,海冰厚度演变规律为:厚度越大,其发生概率的时间也越早。大于40 cm海冰厚度的最大发生概率为20世纪50年代,30~40 cm海冰厚度的最大发生概率为20世纪60年代,10~30 cm海冰厚度最大发生概率为20世纪70年代,小于10 cm海冰厚度最大发生概率为20世纪90年代。

      (3)给出1951年以来10年季海冰厚度的分布概率:一般情况下,10年内冰情的严重程度呈轻冰年与重冰年的波动变化,各海冰厚度区间的概率呈先增大后减小的规律;而20世纪70年代为持续的常冰年冰情,90年代为持续的轻冰年冰情,各海冰厚度呈连续减小趋势。

参考文献 (30)

目录

    /

    返回文章