下边界条件对多年冻土温度场变化数值模拟的影响.docxVIP

下边界条件对多年冻土温度场变化数值模拟的影响 摘要：在气候变暖背景下，北半球多年冻土呈现不同程度的退化趋势，冻土升温、活动层增厚、地下冰消融改变了区域工程地质条件、地形地貌，不仅对寒区环境和工程稳定性造成潜在的威胁，还影响着这些地区的气候、水文和生态过程。因此，准确评估和预估多年冻土热状况的变化具有重要科学和实践意义。现有用于模拟多年冻土热状况的各类模式重点考虑了近地表温度场变化对多年冻土的影响，主要集中于对气温和浅表层物理过程和参数化方案等改进和优化，而对于下边界条件设置对多年冻土热状况模拟的影响少有讨论。基于一维热传导冻土模型，以五道梁地区的多年冻土为研究对象，通过设置不同的下边界方案进行模拟实验，定量评估百年尺度气候变化下不同下边界条件对多年冻土温度场变化数值模拟的影响。结果表明：近地表层（3m）的温度场完全由年际气候变化决定，浅层（3～15m）及中层（15～30m）的多年冻土温度场受下边界条件的影响逐渐显著，深层（30m）的地温对百年尺度气候变化的响应不仅与气候变化的幅度有关，还与多年冻土相变热的多少有直接的关系。下边界条件不恰当的设置方式会对大尺度的气候变化下多年冻土消融程度的计算造成较大的影响，进而可能对深层地温乃至多年冻土区面积变化造成严重的误判。因此，开展百年尺度多年冻土温度场变化模拟时，应采用深层或多年冻土底板以下融土层的稳定地热流作为下边界条件。 0 引言多年冻土约占北半球陆地面积的25%瞬态数值冻土模式是研究多年冻土温度场变化的一个有效手段，因为其能灵活描述多年冻土剖面的异质性、几何形状和边界条件因此，本文拟利用藏北高原冰冻圈站在青藏高原五道梁地区架设的综合观测站的长期地温监测资料，基于一维热传导冻土模型，以重建过去50年的地表温度序列以及未来两种升温情景作为上边界驱动条件，在控制其他变量一致的基础上，通过设计不同的下边界方案进行模拟实验，以定量评估下边界条件对多年冻土地温场数值模拟的影响。分析结果可为冻土模拟研究中下边界条件的选取提供一定参考。1 研究区与数据青藏高原五道梁地区（35°13′N、93°05′E，海拔4 783m）位于昆仑山和唐古拉山脉之间的可可西里地区的东边（图1）。该地区位于高原北部连续多年冻土区，年均气温-5.1℃，年降水量为284.8mm，下垫面为盖度稀疏的高寒荒漠草原。其气候特征与多年冻土特征有较好的代表性藏北高原冰冻圈站在五道梁架设了综合观测站，对五道梁气候变化、地气物质能量交换以及多年冻土的动态变化进行了长期的监测及研究2 多年冻土数值模型2.1 考虑相变的一维热传导控制方程假设土壤中的热量传输完全由傅里叶定律所决定，忽视热对流及热辐射的影响。模型的热传导控制方程表示为式中：c通过显热容法考虑考虑土壤中的相变过程以及冻土融土热参数的差异。c式中：下标f与u分别表示土壤的冻结与融化状态；L为单位体积冰水相变潜热（334.54MJ·m从五道梁观测站现场采集的岩芯土壤样品进行实验分析，主要包括土壤质地、干容重和含水量等。土壤样品通过直径5cm、高5cm的采样器从探坑或钻孔岩芯中采集。土壤质地由德国UGT科学仪器公司生产的SEDIMAT4-12粒度分析仪进行分析，其测量原理为湿筛法结合移液管法。干容重和含水/冰量采用烘干法（60℃，24h）测定。模型中所有的土壤热参数基于以上测定的土壤性质资料，参考《冻土地区建筑地基基础设计规范》选取。模型中土壤参数取值见表1。2.2 上边界条件为充分考虑历史气候对多年冻土模拟的影响，本文提取位于五道梁的国家气象站（观测数据（http：///）中相应时间段的日气温数据，与综合观测场的日地表温度观测数据进行线性拟合，关系式为式中：GST为五道梁综合观测场的实测地表温度；X为五道梁国家气象台站的历史气温序列。然后，根据式（4）与五道梁的国家气象站自1965年的历史日气温数据序列重建过去50年（1965—2015年）的日地表温度序列（图2）。上边界的地表温度变化T用下面正弦方程进行拟合。式中：T此外，根据第5次耦合模式比较计划（CMIP5)11个气候模式的预测结果，2011—2100年中国气温升温率为0.006～0.063℃·a2.3 数值求解空间步长设置为多年冻土垂直剖面5m以上取0.05m,5m以下取0.5m；时间步长设置为1d。热传导偏微分方程通过MATLAB中ode23t函数（自由内插法的梯形法则）求解。3 实验设计3.1 9种不同下边界方案现有冻土数值模拟计算中，下边界条件的设置差异主要包括设置深度不同和设置方式不同。综合现有主要的冻土模式对于下边界条件的设置，本文设计3个下边界深度：10m（大气环流模式中的冻土层一般设置的范围3.2 实验步