极地冰盖粒雪密实化过程的研究进展

粒雪密实化过程在冰川物质平衡监测和冰芯古气候记录分析等冰川学研究中具有重要作用,近年来国内外学者在极地冰盖开展了大量工作并建立出一系列粒雪密实化模型。本文将典型的极地冰盖粒雪密实化模型分为经验模型（包括半经验模型）和物理模型两大类,综合分析和讨论了两类模型的建模机理、发展现状及不足之处。经验模型关注粒雪柱在密实化过程中的宏观演变,利用实测的粒雪深度-密度数据建立密实化速率方程。目前,该类模型主要研究了温度、积累率、融水渗透和再冻结过程对密实化速率的影响。物理模型聚焦于粒雪晶粒在密实化过程中的微观结构变化,利用晶界滑移、幂律蠕变和晶格扩散等物理机制解释晶粒的生长变形过程。目前,新雪-粒雪转化过程的研究主要基于晶界滑移理论,粒雪-冰川冰转化过程的研究主要基于球形粉末的压力烧结理论。总体而言,极地冰盖粒雪密实化过程的研究在过去几十年取得了长足的发展,但由于粒雪微观结构及其与宏观演变联系的认知匮乏,经验模型缺乏完善的湿粒雪密实化理论,且将稳态假设直接推广到瞬态情景存在一定问题;物理模型由于粒雪密实化的物理机制尚不明确,难以精确刻画粒雪演变的整个物理过程。     

国际南极冰盖与海平面变化研究述评

海平面上升是全球变暖的主要后果之一。尽管有少数冰川学家认为,气候变暖并不能确保雪积累量的显著增加,同时可能出现冰流的突然变化,因此南极冰盖在未来海平面变化中的作用存在很大的不确定性。但近几十年来南极半岛气温的急剧上升,已使大量的冰架崩解。冰架崩解并不对海平面产生真正的影响,但反映出南极洲气候与冰川存在急剧变化的可能。     

末次盛冰期气候模拟及青藏高原冰盖的可能影响

在国际古气候模拟比较计划设定的边界条件下(其中包括:海洋表面温度、陆地冰盖、CO2浓度、地球轨道参数等变化),本文汇报了利用中国科学院大气物理研究所的9层全球大气环流格点模式对末次盛冰期气候进行数值模拟的研究成果。结果表明:末次盛冰期全球年平均温度比现在低5.3℃,降雨量下降9%,其中大陆年平均降水量仅为目前的71%。通过与其他模式模拟结果和重建古气候资料的对比表明,该模式捕捉到了末次盛冰期干、冷的冰期气候特征。在此基础上,我们根据刘东生等的科研成果,初步探讨了青藏高原在末次盛冰期倘若存在一定范围的冰盖,那么该冰盖会对东亚区域古气候模拟结果产生怎样的冲击。     

南极冰盖GLIMMER模式移动边界试验研究

利用欧洲冰盖比较计划(EISMINT)的移动边界条件EISMINT-1测试集,通过三维有限差分南极冰盖GLIMMER模式,采用浅冰近似计算了一个定义在方形冰盖积分区域的流场特征量,研究了冰盖在大尺度长时间序列条件下对气候变化的反馈,考察了稳定态下冰盖演化的周期性行为.试验表明:GLIMMER模式在EISMINT-1移动边界条件下可以模拟出气候强迫下的冰盖演化趋势,在长时间尺度上可以刻画出气候变化的反馈机制.同时,对南极冰盖数值模拟面临的问题、GLIMMER模式以及EISMINT冰盖比较计划的内容和框架做了说明,对EISMINT比较计划目前的发展状况和遇到的问题做了分析,并且对GLIMMER模式下一步的研究方向做了一些探讨.     

冰期指数法模拟北半球冰盖演化的不确定性研究

IPCC第六次评估报告（AR6）显示,自20世纪起极地冰盖持续消融,全球海平面不断上升。目前对于地球冰盖未来的预测以及过去的演变历史尚不明确,而数值模拟能够提供一种有效的解决方案。在冰盖模拟研究中,冰期指数法可依据古气候代用指标将离散的气候强迫转化为连续的气候强迫,用于冰盖演变的瞬态模拟。基于该方法,利用2组（共6条）分别代表全球海平面和温度变化的代用指标,开展末次冰期旋回北半球冰盖的时空演变模拟研究,结果表明：(1)模拟的冰量演变特征受指数的变化趋势控制;(2)在指数变化特征（轨迹和变幅）相似时,千年尺度气候突变事件的存在会导致模拟的总冰量偏少;(3)在同一气候强迫下,不同指数模拟的最大冰盖范围受夏季0℃等温线的约束,同时指数的演变轨迹与变幅也会影响末次冰期盛冰期冰盖模拟的空间分布。因此,在利用冰期指数法开展冰盖瞬态模拟研究时,需根据关注的研究区域选取有代表性的指数并考虑古气候代用指标（即冰期指数）的不确定性对模拟结果的影响。     

冰盖下水流紊动特性试验研究

河冰是陆地冰冻圈的重要组成要素,在生消演变过程中影响水流运动规律,受冰凌影响河流流态可分为明流、完全冰封和部分冰封。为探究冰盖下水流的紊动特性,采用声学多普勒三维测速仪测量3种流态的恒定均匀流场,对比分析纵向时均流速、雷诺应力和紊动强度的分布规律。结果表明：完全冰封流纵向时均流速垂向上呈不对称“■”型分布,雷诺应力在垂向上线性分布,紊动强度在基于双层假定划分的冰盖层和床面层均符合指数分布律;部分冰封流纵向时均流速垂向分布兼具明流和完全冰封流的特性,存在横向动量交换,导致冰盖区、明流区和过渡区的雷诺应力和紊动强度在水深方向上也存在明显的差异。试验结果可为冰凌灾害防治和供水安全保障提供理论指导。     

冰盖下冰花浓度分布的数值模拟

封冻河道内的冰花分布规律对冰期水文及环境分析有着重要的意义. 视冰花为连续介质, 基于多相流理论和欧拉-欧拉模型方法, 建立了河渠垂向二维标准k-ε紊流数值模型. 在不同的水流流速、 冰花颗粒尺寸、 上游来冰量条件下, 对冰盖下冰花浓度场进行了模拟分析, 得到了冰盖下冰花分布特征. 结果表明: 水流流速越大, 冰花沿垂向分布范围越大, 冰花浓度分布曲线斜率越大, 冰花浓度沿垂向的均匀性越大, 同时浓度最大值越小; 冰花颗粒尺寸越大, 冰花沿垂向分布范围越大, 冰花浓度分布曲线斜率越大; 上游来冰量越大, 冰花沿垂向分布范围越大, 在同一水深其浓度值也越大, 冰花浓度沿垂向分布的均匀性越大.     

东南极冰盖/冰架变化监测与预测技术研究

南北极是研究全球变化的关键区域。＂十一五＂期间我国在南极地区系统开展了东南极冰盖/冰架变化监测与预测技术研究,对认识全球气候变化具有重要作用。通过项目实施,建立和发展了一批现场观测体系,发展了冰盖观测新技术并集成应用于东南极冰盖的综合观测,拓展了对冰穹A地区的新认识和新发现;在冰穹A边缘地区钻取的一支浅冰芯恢复了过去约780年的气候记录,记录了东南极地区存在小冰期的明显证据;发展了东南极冰盖积累和等时年层流动模型,研究在冰盖浅层、中层和深部的变化特征,反演了冰穹A地区的古积累率分布。本文概要介绍该项目基本情况。     

冰盖下水流速度分布的二维数值模拟分析

基于标准的κ-ε模型,利用有限体积法,对冰盖下水流速度分布进行了二维数值模拟和分析.结果表明:当冰盖和河床糙率比值变化时,水流速度剖面发生变化;冰盖下的最大流速值比速度对称分布时要大,且偏于光滑的一边;在相对光滑壁处近壁流速值大于相应对称分布时的流速值,该结果与实验室实测速度剖面情况吻合.数值试验表明,最大流速点位置的偏移量仅在有限的范围内,冰盖区、床面区及冰盖下水流平均速度相差较小,表明了工程应用中假定各区流速相等是可行的.     

1961—2017年南极冰盖近地面风时空变化研究

基于19个人工气象站1961—2017年风速风向实测数据对南极冰盖近地面风速时空变化和风向进行了分析。结果表明：近50年来,南极冰盖近地面各季节平均风速和年平均风速变化的空间模式基本一致。东南极0°~120° E沿海地区风速呈显著上升趋势;南极半岛气象站风速变化趋势各异,且变化速率相差甚大,但从区域平均结果来看,南极半岛年和季节平均风速均呈上升趋势。这与近几十来局地气温、气压变化及南半球环状模趋向于正位相发展有关。东南极受下降风和绕极东风影响,大部分地区盛行偏南风或偏东风,且频率较高风向稳定;而受天气活动影响,南极半岛风向复杂,主风向频率低,风向多变。     

丰满水库水温的原型观测及分析

为探究寒区冰封水库热状况的时空演变规律,于2010年6月~2011年8月对丰满水库的水温和冰情进行了原型观测。基于现场观测资料,阐述了丰满库区不同季节热结构的变化规律,重点分析了汛期洪水过程对分层结构的影响,并发现解冻初期水库表层浮力流动现象。同时,典型水温分层期1月和6月的坝前75~800m范围4条水温观测垂线基本一致,且同期电站下泄水温与坝前取水口对应高程平均水温差异较小,说明该水域范围内无三维水温效应。由于深层取水影响,丰满电站下泄水温过程较为稳定,但与坝址天然水温相比,5~7月、10~12月分别表现出月平均高达10.8℃、7.9℃的春夏低温水和秋冬高温水现象。此外,冰情监测结果表明,丰满水库冬季基本全库封冻,坝前最大冰厚约0.7m,且库区冰厚沿程出现不均匀分布的规律。     

草皮护坡水力糙率实验研究

生态护岸工程可以在保护河岸的同时对河道生态系统起到良好的维护作用,其中草本植物因其生长快、易存活成为生态护岸的重要工程材料。草皮护坡糙率实验对36组草皮试件的糙率进行了系统的实验研究。研究表明,草皮糙率的大致范围为0.020~0.090,主要与土壤类型、肥沃程度、水流条件等有关。通过实验数据分析发现:草皮糙率随生长期而增加;加筋草皮(基土中设置三维加筋材料的草皮)的糙率普遍大于天然草皮糙率。     

南水北调中线输水工程若干冰力学问题试验研究

通过在低温冰工程实验室进行的冰力学模型试验,研究了南水北调中线工程冰期输水中的若干冰力学问题,主要包括:冰的膨胀力、冰盖稳定性和浮冰对闸墩的撞击力.试验结果给出了冰的膨胀力变化典型时程曲线、不同厚度冰盖的膨胀力极值变化;不同厚度的冰盖随水位下降的最大塌落值;以及冰排对闸墩的撞击力的冰力时程曲线、冰力极值随冰速的变化曲线.这些试验成果可为南水北调中线工程冰期输水安全运行提供参考.     

环渤海海冰弯曲强度的试验测试及特性分析

在渤海冰区油气开发中,海冰物理力学性质对海洋结构设计、海冰动力学过程均有重要影响.目前,随着渤海冰区锥体海洋平台的增加,对海冰弯曲强度的研究具有重要意义.对环渤海沿岸9个测点的海冰弯曲强度进行了现场和室内测试,同时对影响海冰弯曲强度的海冰盐度、温度进行了测试;分析了莱州湾、辽东湾西岸和辽东湾东岸等不同测点海冰弯曲强度的...     

猫跳河流域土地覆被变化及其驱动力定量研究

应用RS、GIS以及GPS集成技术,以黔中喀斯特地区的猫跳河流域为对象,研究了这一地区从1973—2002年间土地利用／覆被的变化过程,并利用典型相关分析法研究了形成这种变化的驱动因素。研究表明,从1973—1990年,猫跳河流域耕地、建设用地、有林地、水域以及裸岩地呈增长趋势,而灌木林和灌草地趋于减少,流域生态环境质量下降。进入20世纪90年代后,有林地、灌草地、建设用地、水域持续增加,灌草地总量变化不大,耕地总量和裸岩地减少,流域生态环境质量趋于好转。定量分析表明,1990年以前,驱动建设用地变化的因素是总人口的变化,1990年以后则是农业产值;1990年以前驱动耕地变化的是农业人口的增长,1990年以后则和粮食产量的提高有关。      

渤海海冰断裂韧度试验

海冰的断裂韧度是近海结构冰荷载分析的重要参数,也是影响海冰断裂、重叠和堆积等动力行为的主要力学特性.2011-2012年冬季现场采集了渤海莱州湾的海冰试样,在室内低温环境下通过三点弯曲试验测得了不同温度、盐度和加载速率下的海冰断裂韧度KIC.结果表明,卤水体积(温度和盐度)和加载速率对KIC的影响比较显著.随温度和加载速率的降低,KIC值均会增大,且与卤水体积的平方根呈负指数关系.为极限破坏状态下海冰的断裂过程和冰荷载分析的参数选用提供了依据.     

渤海海冰断裂韧度试验

海冰的断裂韧度是近海结构冰荷载分析的重要参数,也是影响海冰断裂、重叠和堆积等动力行为的主要力学特性。2011—2012年冬季现场采集了渤海莱州湾的海冰试样,在室内低温环境下通过三点弯曲试验测得了不同温度、盐度和加载速率下的海冰断裂韧度K_IC。结果表明,卤水体积(温度和盐度)和加载速率对K_IC的影响比较显著。随温度和加载速率的降低,K_IC值均会增大,且与卤水体积的平方根呈负指数关系。为极限破坏状态下海冰的断裂过程和冰荷载分析的参数选用提供了依据。     

Inclusion of ice cover in a storm surge model for the Beaufort Sea

Storm surges in the Beaufort Sea present a severe problem for navigation as well as for offshore oil drilling activities. Influence of ice cover on storm surges in the Beaufort Sea is examined making use of a numerical model as well as a set of observations.The automated shallow-water model of Henry has been modified to incorporate ice cover and is adapted to the Beaufort Sea. The leading edge of the permanent ice is calculated from the loci of identifiable points. Generalized similarity theory is employed to compute wind stresses. Simulations are made using model-predicted ice concentrations and observed ice concentrations. Ice motion is relatively small in units of model grid distance (approximately 18 km) during surges. Spherical effects are important and should be included in future adaptations of the model. Comparison of the computed surges with observed surges for eight different events showed reasonable agreement.     

2008年夏季北冰洋海冰表面积雪特征初步分析

运用中国第三次北极考察期间观测的积雪资料, 分析了夏季北冰洋中心海域海冰表面积雪的空间分布特征. 结果表明: 积雪在垂直分布上呈现6种粒雪状态, 表层到底层依次为新降雪、风板、冰片、深霜、冻结状粗雪和渗浸冻结冰层, 积雪表面常被新降雪或厚度为2～3 cm的风板所覆盖. 考察区域积雪的平均密度为(304.01±29.00)kg·m^-3, 表层密度略低于次表层, 遵循雪的密实化原理. 积雪厚度, 雪水当量和新降雪皆具有由南向北递减的空间分布特征, 表明在海冰消融末期, 积雪的气候态分布主要是由降水量的多少决定的, 积雪的消融和蒸发并非海冰表面积雪评估需考虑的首要因素. 雪温随深度的增加而增加, 具有明显的垂直梯度, 观测区积雪表面的平均温度为(-2.01～±0.96)℃, 比海冰/积雪界面的温度高得多.     

饱和粉土冻结过程中的水分迁移试验研究

改变边界温度和土样高度,对饱和粉土进行了冻结试验,研究其水分迁移、水分重分布、冻胀和冰透镜体的发展规律。试验结果表明：当温度稳定时,水分持续迁移到冻结锋面附近,含水率急剧增大,易形成冰透镜体。饱和粉土冰透镜体几何形态较为规则,无枝状交叉结构,已冻土为整体状且无网状裂隙。冻结过程中存在起始冻胀时间,在起始冻胀时间之前,土中水分被排出,冻胀发生之后水分补给到冻土中,且补给水分产生的冻胀量和总冻胀量数值接近。土样高度影响水分迁移量和冻胀量,土样越高,冻胀量越小,冻土含水率增量越小,但水分增量分布区域越分散,且起始冻胀时间越长。