融入冻融界面动态的陆面模式构建及与CAS-ESM耦合机制探究

融入冻融界面动态的陆面模式构建及与CAS-ESM耦合机制探究一、引言1.1研究背景陆面过程作为地球系统的重要组成部分，在全球气候变化研究中占据着关键地位。陆面模式则是描述土壤内部水热传输以及陆地-大气交界面物质和能量交换的物理、生化及生态过程的数值模型，其对陆面过程的准确模拟对于理解地球系统的运行机制至关重要。陆面过程涉及能量平衡、水分循环、动量交换以及生物地球化学循环等多个方面。例如，土壤温度和湿度的变化会影响地表的能量通量，进而影响大气的温度和运动；植被的蒸腾作用则是水分循环的重要环节，同时也会对大气的水汽含量和降水产生影响。此外，陆面过程还与生物地球化学循环密切相关，如碳、氮等元素在陆地生态系统中的循环过程会受到陆面过程的影响，进而影响大气中温室气体的浓度。冻融界面变化是陆面过程中的一个重要现象，它对陆面过程有着深远的影响。在冻土地区，土壤的冻融过程会导致土壤物理性质的改变，如土壤孔隙度、渗透率等的变化，从而影响土壤水分和热量的传输。例如，冻结的土壤会阻碍水分的下渗和蒸发，使得水分在地表积聚，增加地表径流的产生；而融化的土壤则会使水分更容易下渗，增加土壤的含水量。冻融过程还会影响土壤中微生物的活动和有机质的分解，进而影响土壤的肥力和碳循环。CAS-ESM（ChineseAcademyofSciencesEarthSystemModel）作为中国科学院研发的地球系统模式，在地球系统模拟中发挥着重要作用。它集成了大气环流、海洋环流、陆面过程、海冰、气溶胶和大气化学、植被动力学、陆地生物地球化学、海洋生物地球化学等多个分系统模式，能够全面地模拟地球系统的各个组成部分及其相互作用。通过与其他地球系统模式的比较和验证，CAS-ESM在模拟全球气候和环境变化方面具有较高的准确性和可靠性，为研究地球系统的演变规律和预测未来气候变化提供了有力的工具。1.2研究目的与意义本研究旨在构建一个能够准确考虑冻融界面变化的陆面模式，并将其与CAS-ESM进行有效耦合，以提高地球系统模拟的精度和可靠性，为深入理解气候变化机制提供更为坚实的理论基础和技术支持。冻融界面变化对陆面过程有着深远的影响，然而，目前的陆面模式在模拟这一过程时仍存在诸多不足。许多陆面模式对冻融界面的处理过于简单，未能充分考虑土壤冻融过程中物理性质的变化，导致对陆面水热传输的模拟存在较大误差。这不仅影响了对陆面过程的准确描述，也降低了地球系统模式对气候变化的模拟和预测能力。本研究通过构建考虑冻融界面变化的陆面模式，能够更准确地模拟土壤冻融过程对陆面水热传输的影响，从而提高陆面模式的模拟精度。通过将该陆面模式与CAS-ESM耦合，可以使地球系统模式更加全面地考虑陆面过程与其他圈层的相互作用，进而提升对气候变化的模拟和预测能力。构建考虑冻融界面变化的陆面模式并与CAS-ESM耦合，有助于我们更深入地理解气候变化机制。陆面过程与大气、海洋等圈层之间存在着复杂的相互作用，这些相互作用在气候变化中起着关键作用。冻融界面变化会影响土壤水分和热量的传输，进而影响大气的温度和湿度，对大气环流和气候产生影响。通过本研究，可以更准确地模拟这些相互作用，揭示气候变化的内在机制，为应对气候变化提供科学依据。这也有助于我们更好地理解陆面过程在地球系统中的作用，为地球系统科学的发展做出贡献。在全球气候变化的背景下，提高地球系统模拟的精度和可靠性具有重要的现实意义。准确的地球系统模拟可以为政府制定应对气候变化的政策提供科学依据，帮助决策者更好地了解气候变化的趋势和影响，从而制定出更加有效的应对策略。在农业领域，了解气候变化对土壤水分和温度的影响，可以帮助农民合理安排种植计划，提高农作物的产量和质量；在水资源管理方面，准确的气候预测可以帮助政府合理规划水资源，应对可能出现的干旱和洪涝等灾害。因此，本研究对于应对气候变化、保障人类社会的可持续发展具有重要的意义。1.3研究现状与问题分析陆面模式的发展历程是一个不断演进和完善的过程。自20世纪60年代末至70年代，第一代陆面模式诞生，其基于简单的“水桶”假定，对土壤含水量、地表蒸发和地表径流进行简单描述，虽能初步模拟陆面基本水分收支，但过于简化的假设无法准确反映真实陆面过程的复杂性。70年代末到90年代，第二代陆面模式引入植被生理过程，将植被冠层近似处理为薄叶片，建立了更细致的陆面参数化方案，能够较为真实地描述植被在陆面过程中的作用，以及植被蒸腾和冠层辐射传输对陆面水分和能量收支的影响，像BATS和SiB模式便是这一时期的典型代表，极大提升了对陆面过程中植被相关影响的模拟能力。90年代以后，第三代陆面模式引入植被光合作用的生物化学过程，不仅能模拟陆面水热交换，还能刻画地表碳通量和大气中CO₂浓度变化，为植物动态生长响应气候变化的生态模式研究奠定基础，NCARLSM和SiB2是该阶段的重要模式，开启了陆面模式对碳循环模拟的新篇章。90年代末以后发展的第四代陆面模式，如CoLM、CLM、JULES、CABLE和Noah-MP等，加入大量新子模块，对陆面物理、化学和生物等过程的描述更加精细化，涵盖的过程也更加完备，能从更多维度、更细致地模拟陆面复杂过程。在冻融界面变化研究方面，目前已取得了一定的进展。学者们通过野外观测、实验研究和数值模拟等多种手段，对土壤冻融过程中的物理机制有了更深入的认识。研究发现，土壤的冻融过程会导致土壤孔隙结构的改变，进而影响土壤的水力传导率和热导率。冻结过程中，土壤中的水分结冰膨胀，会使土壤孔隙增大；而融化过程中，土壤孔隙又会因冰的融化而发生变化。这种孔隙结构的变化对土壤水分和热量的传输有着重要影响。有研究表明，在冻土区，土壤冻结时，水分会向冻结锋面迁移，导致冻结锋面附近土壤含水量增加，而融化时，水分则会重新分布，影响土壤的干湿状况。冻融过程还会对土壤中微生物的活动和有机质的分解产生影响。在低温冻结条件下，微生物的活性受到抑制，有机质的分解速率减缓；而在土壤融化时，微生物活性增强，有机质分解加速，这会影响土壤中碳、氮等元素的循环，进而对陆地生态系统的碳收支产生影响。在一些高纬度地区，冻土中储存着大量的有机碳，随着气候变暖，冻土融化，这些有机碳可能会被分解并释放到大气中，成为重要的温室气体排放源，进一步加剧全球气候变化。尽管在冻融界面变化研究方面取得了上述成果，但当前陆面模式在模拟冻融界面变化及其对陆面过程的影响时仍存在一些问题。许多陆面模式对冻融界面的处理过于简单，未能充分考虑土壤冻融过程中物理性质的动态变化，导致对陆面水热传输的模拟存在较大误差。一些模式在模拟土壤冻融过程时，没有考虑土壤孔隙结构变化对水力传导率和热导率的影响，使得模拟结果与实际情况存在偏差，从而影响了对陆面水分和热量收支的准确模拟。在陆面模式与CAS-ESM耦合方面，也存在一些挑战。陆面模式与CAS-ESM中其他分系统模式之间的耦合机制尚不完善，导致各分系统之间的相互作用不能得到准确的模拟。在陆面与大气的耦合过程中，能量和物质的交换过程可能存在参数化不合理的问题，影响了整个地球系统模式对气候变化的模拟能力。陆面模式中一些参数的不确定性也会对耦合后的模拟结果产生影响，如何准确确定这些参数，提高陆面模式与CAS-ESM耦合的精度，是当前需要解决的重要问题。综上所述，当前陆面模式在模拟冻融界面变化及其与CAS-ESM耦合方面仍存在不足，亟待进一步改进和完善。本研究旨在针对这些问题，构建考虑冻融界面变化的陆面模式，并实现其与CAS-ESM的有效耦合，以提高地球系统模拟的精度和可靠性。二、陆面模式基础与冻融界面变化2.1陆面模式的基本原理与分类陆面模式作为描述土壤内部水热传输以及陆地-大气交界面物质和能量交换的物理、生化及生态过程的数值模型，在地球系统科学研究中扮演着重要角色。其基本原理基于对陆面过程中各种物理、生化和生态机制的数学描述，通过建立一系列的方程和参数化方案，来模拟陆面与大气之间的相互作用。陆面模式主要由物理过程、生化过程和生态过程三个部分构成。物理过程主要描述植被、土壤及雪盖等下垫面的水、热传输以及与大气进行水分、热量、动量和辐射的交换。在水热传输方面，通过能量守恒方程和质量守恒方程来描述土壤中热量和水分的运动。热量传输主要通过热传导和对流进行，而水分传输则涉及土壤水分的入渗、蒸发、蒸腾以及地下水的运动等。辐射交换过程则考虑了太阳辐射的吸收、反射和散射，以及地表长波辐射的发射和吸收，这些过程对地表能量平衡有着重要影响。生化过程主要描述碳、氮等化学元素在陆地生态系统中的循环，包括植被的光合作用以及土壤的呼吸等。光合作用是植被利用光能将二氧化碳和水转化为有机物和氧气的过程，这一过程不仅影响着植被的生长和发育，还对大气中二氧化碳的浓度有着重要的调节作用。土壤呼吸则是土壤中微生物和植物根系将有机物氧化分解，释放出二氧化碳的过程，它是陆地生态系统碳循环的重要组成部分。生态过程主要描述时间尺度较长的生态系统的结构、组成和成分的动态变化以及生态系统的演替等。随着时间的推移，生态系统中的植被类型、生物量以及物种组成等都会发生变化，这些变化会影响陆面过程中的能量和物质交换，进而影响整个地球系统的运行。陆面模式的发展经历了多个阶段，不同代际的陆面模式在结构和功能上存在显著差异。第一代陆面模式诞生于20世纪60年代末至70年代，基于简单的“水桶”假定，对土壤含水量、地表蒸发和地表径流进行简单描述。这类模式假设土壤为垂直方向均匀分布且容水量达15厘米的“水桶”，根据水分质量守恒的假设，用空气动力学总体输送公式和几个均匀的陆地表面参数（如反照率、粗糙度）来计算土壤含水量、地表蒸发和地表径流。这种简单的模型虽然能够初步模拟陆面的基本水分收支，但由于其过于简化的假设，无法准确反映真实陆面过程的复杂性，对土壤水热传输的模拟精度较低。70年代末到90年代发展起来的第二代陆面模式，在土壤水、热传输方案的基础上，引入了植被生理过程。将植被冠层近似地处理为一个薄的叶片，基于一系列可以直接观测到的陆面参数，建立起关于植被覆盖上空辐射、水分、热量和动量交换以及土壤中水热传输过程的参数化方案。这些模式能较为真实地描述植被在陆面过程中的作用，并能够细致地刻画植被蒸腾及植被冠层辐射传输对陆面水分和能量收支的影响。像罗伯特・狄金森研发的BATS（BiosphereAtmosphereTransferScheme）和皮尔斯・塞勒斯研制的SiB（SimpleBiosphereModel）便是这一时期的典型代表。BATS模式考虑了植被冠层对辐射的截留和散射，以及植被蒸腾对水分和能量的影响，能够更准确地模拟陆面的能量和水分收支。SiB模式则进一步完善了植被生理过程的描述，对植被光合作用和呼吸作用进行了参数化，提高了对陆面碳循环的模拟能力。90年代以后，第三代陆面模式逐渐发展起来，引入了植被吸收二氧化碳进行光合作用的生物化学过程。这类模式不仅能更加真实地刻画陆地-大气交界面的水、热交换过程，而且还可以模拟地表碳通量和大气中二氧化碳浓度的变化，从而为植物动态生长并响应气候变化的生态模式研究打下基础。美国国家大气研究中心研发的陆面模型（NCARLSM）和SiB的升级版本SiB2是该阶段的代表。NCARLSM模式考虑了植被的光合作用、呼吸作用以及土壤中碳的分解和积累等过程，能够模拟陆地生态系统中碳的循环和储存。SiB2模式则在SiB的基础上，进一步改进了对植被生理过程的描述，提高了对陆面碳通量和大气中二氧化碳浓度变化的模拟精度。90年代末以后发展的第四代陆面模式，如CoLM（CommonLandModel）、CLM（CommunityLandModel）、JULES（theJointUKLandEnvironmentSimulator）、CABLE（theCommunityAtmosphereBiosphereLandExchangeModel）和Noah-MP（Noah-MultiparameterizationLandSurfaceModel）等，加入了大量新的子模块，对陆面物理、化学和生物等过程的描述更加精细化，涵盖的过程也更加完备。CoLM模式考虑了土壤冻融过程对水热传输的影响，以及植被动态生长对陆面过程的反馈作用，能够更全面地模拟陆面过程。CLM模式则进一步完善了对陆面碳、氮循环的描述，考虑了土壤中微生物的活动和营养元素的循环，提高了对陆地生态系统功能的模拟能力。JULES模式注重对植被生态过程的模拟，考虑了植被的物候变化、生产力和碳分配等过程，能够更准确地模拟植被对气候变化的响应。CABLE模式则强调了陆面与大气之间的物质和能量交换，对陆面过程中的水文循环、能量平衡和生物地球化学循环进行了综合模拟。Noah-MP模式则在Noah模式的基础上，增加了多个参数化方案，提高了对不同下垫面条件下陆面过程的模拟能力。不同代际的陆面模式在模拟能力和应用范围上逐渐拓展和深化。随着对陆面过程认识的不断深入，陆面模式的结构和功能也在不断改进和完善，以更好地模拟复杂的陆面过程及其与大气的相互作用。2.2冻融界面变化的过程与影响因素冻融界面变化是一个复杂的物理过程，涉及到土壤中水分的相变、热量的传递以及土壤物理性质的改变。当土壤温度降至冰点以下时，土壤中的水分开始结冰，体积膨胀，导致土壤孔隙结构发生变化。随着温度的进一步降低，冻结锋面向下推进，冻融界面逐渐加深。在这个过程中，土壤的导热率、热容量和水力传导率等物理性质也会发生显著变化。土壤的导热率在冻结过程中会增大，因为冰的导热率远高于水和空气，这使得热量更容易在土壤中传导。而热容量则会减小，因为冰的热容量相对较低。土壤的水力传导率会大幅降低，因为孔隙被冰填充，阻碍了水分的流动。当土壤温度回升至冰点以上时，冰开始融化，水分重新分布，土壤孔隙结构再次发生改变。融化过程中，土壤的导热率和热容量会逐渐恢复到未冻结状态，而水力传导率则会随着孔隙的恢复而增加。冻融界面变化受到多种因素的影响，其中气温是最为关键的因素之一。气温的高低直接决定了土壤是否会发生冻结和融化，以及冻结和融化的程度。在寒冷的季节，当气温持续低于冰点时，土壤会逐渐冻结，冻融界面不断加深；而在温暖的季节，气温升高，土壤开始融化，冻融界面逐渐上升。气温的波动也会对冻融界面变化产生影响。频繁的气温升降会导致土壤反复冻结和融化，加速冻融界面的变化过程，对土壤结构和水分传输产生更为复杂的影响。在高纬度地区，冬季气温极低且持续时间长，土壤冻结深度较大，冻融界面变化相对缓慢；而在中低纬度地区，气温波动较大，土壤冻融过程更为频繁，冻融界面变化也更为活跃。降水对冻融界面变化也有着重要影响。降水为土壤提供了水分来源，而土壤水分含量是影响冻融过程的关键因素之一。充足的降水会增加土壤的含水量，使得土壤在冻结时能够形成更多的冰，从而影响冻融界面的变化。在降水较多的地区，土壤含水量较高，冻结时冰的形成量较大，冻融界面的变化可能更为剧烈。降水的时间分布也会影响冻融界面变化。如果降水主要集中在冬季，且在土壤冻结之前，那么土壤含水量会在冻结前增加，导致冻结过程中冰的形成量增加，冻融界面加深。相反，如果降水在土壤冻结后发生，可能会导致积雪覆盖增加，对土壤起到保温作用，减缓冻融界面的变化。积雪覆盖是影响冻融界面变化的另一个重要因素。积雪具有良好的隔热性能，能够有效地阻挡土壤与大气之间的热量交换。在冬季，积雪覆盖在土壤表面，就像给土壤盖上了一层棉被，使得土壤温度相对稳定，不易受到外界气温变化的影响。当积雪厚度较大时，土壤表面的温度会保持在相对较高的水平，从而减缓土壤的冻结速度，使冻融界面变浅。而在春季，随着气温升高，积雪开始融化，融化的雪水会渗入土壤，增加土壤的含水量，进一步影响冻融界面的变化。在一些高海拔地区，积雪覆盖时间长、厚度大，对冻融界面变化的影响尤为显著。积雪的存在可以延长土壤的冻结期，减少土壤的冻融循环次数，对土壤的物理性质和生态系统产生深远的影响。植被覆盖也会对冻融界面变化产生影响。植被通过蒸腾作用调节土壤水分含量，同时其根系和枯枝落叶层可以改变土壤的物理性质，进而影响冻融过程。植被的蒸腾作用会消耗土壤中的水分，降低土壤的含水量，从而减少土壤冻结时冰的形成量，影响冻融界面的变化。植被的根系可以增加土壤的孔隙度，改善土壤的通气性和透水性，这也会对土壤的冻融过程产生影响。植被的枯枝落叶层可以起到保温和保水的作用，类似于积雪覆盖，能够减缓土壤的冻融速度，使冻融界面变化相对稳定。在森林地区，茂密的植被覆盖使得土壤受到的外界影响较小，冻融界面变化相对较为缓和；而在草原地区，植被覆盖相对较稀疏，土壤的冻融过程可能会受到更多外界因素的影响，冻融界面变化更为复杂。土壤质地和结构也是影响冻融界面变化的重要因素。不同质地的土壤，如砂土、壤土和黏土，其孔隙大小、分布和连通性不同，这会影响土壤中水分的存储和传输，进而影响冻融过程。砂土的孔隙较大，通气性和透水性好，但保水性较差，土壤中的水分容易流失，在冻结时冰的形成量相对较少，冻融界面变化相对较快。黏土的孔隙较小，保水性强，但通气性和透水性较差，土壤中的水分不易排出，冻结时冰的形成量较大，冻融界面变化相对较慢。壤土的性质介于砂土和黏土之间，其冻融界面变化的特点也相对适中。土壤的结构，如团聚体的大小和稳定性，也会影响土壤的冻融过程。团聚体结构良好的土壤，孔隙分布均匀，有利于水分和热量的传输，冻融界面变化相对较为稳定；而结构松散的土壤，孔隙分布不均匀，水分和热量传输不畅，冻融界面变化可能更为复杂。地形因素，如海拔、坡度和坡向，也会对冻融界面变化产生影响。海拔高度决定了气温和降水的分布，随着海拔的升高，气温降低，降水增加，这会导致高海拔地区的土壤更容易冻结，冻融界面更深。坡度和坡向则影响太阳辐射的接收和地表径流的分布。阳坡接收的太阳辐射多，温度较高，土壤冻结程度相对较轻，冻融界面较浅；而阴坡接收的太阳辐射少，温度较低，土壤冻结程度相对较重，冻融界面较深。坡度较大的地区，地表径流速度快，土壤水分容易流失，冻结时冰的形成量相对较少，冻融界面变化相对较快；而坡度较小的地区，土壤水分相对稳定，冻融界面变化相对较慢。冻融界面变化是一个受多种因素共同影响的复杂过程。深入了解这些因素对冻融界面变化的影响机制，对于准确模拟陆面过程中的土壤水热传输具有重要意义。2.3冻融界面变化对陆面模式关键过程的影响2.3.1对土壤水热传输的影响冻融过程对土壤孔隙结构有着显著的影响。当土壤中的水分冻结时，水变成冰，体积膨胀约9%，这会对土壤孔隙产生挤压作用。在冻结初期，土壤孔隙中的水分开始结冰，冰晶体逐渐生长，占据孔隙空间，使得土壤孔隙被撑大。随着冻结程度的加深，冰的体积不断增大，孔隙结构进一步改变，原本连通的孔隙可能会被冰隔断，形成孤立的孔隙，导致土壤孔隙的连通性降低。这种孔隙结构的变化对土壤水分传输产生了重要影响。在水分传输方面，土壤孔隙结构的改变直接影响了土壤的水力传导率。水力传导率是衡量土壤中水分传输能力的重要参数，它与土壤孔隙的大小、形状和连通性密切相关。在冻结过程中，由于孔隙被冰填充，孔隙空间减小，连通性变差，土壤的水力传导率大幅降低。研究表明，当土壤冻结时，水力传导率可能会降低几个数量级。这使得水分在土壤中的运动变得困难，水分的下渗和蒸发过程受到阻碍。在冻结的土壤中，水分难以向下渗透，导致地表积水增加；同时，由于水分蒸发需要通过孔隙扩散到大气中，孔隙结构的改变使得水分蒸发速率减慢。在融化过程中，冰逐渐融化成水，土壤孔隙结构再次发生变化。冰的融化使得孔隙空间恢复，连通性增强，土壤的水力传导率逐渐增大。水分开始重新分布，原本被冻结在土壤中的水分得以释放，下渗和蒸发过程逐渐恢复。融化过程中，土壤孔隙的变化并非完全恢复到冻结前的状态，可能会留下一些残留的孔隙结构变化，这仍然会对水分传输产生一定的影响。冻融过程对土壤热量传输同样有着重要影响。土壤的热导率是影响热量传输的关键参数，它反映了土壤传导热量的能力。在冻结过程中，由于冰的热导率远高于水和空气，土壤的热导率会增大。这使得热量在土壤中的传导速度加快，更容易从高温区域向低温区域传递。当土壤表面温度较低时，热量能够更快地从土壤深层传导到表面，导致土壤表面温度下降更快，进一步促进土壤的冻结。融化过程中，随着冰的融化，土壤的热导率逐渐恢复到未冻结状态。由于融化过程中水分的相变需要吸收热量，这会消耗土壤中的热量，使得土壤温度升高缓慢。在融化初期，土壤中的冰开始融化，吸收大量的潜热，导致土壤温度在一段时间内保持相对稳定，即使外界气温升高，土壤温度的上升也较为缓慢。这种热量传输的变化对土壤的温度分布和冻融界面的移动有着重要影响。土壤冻融过程中，水分和热量的传输并非孤立进行，而是相互耦合的。水分的相变会伴随着热量的吸收和释放，从而影响土壤的热量传输；而热量的传输又会影响土壤的温度分布，进而影响水分的冻结和融化过程。在冻结过程中，水分结冰释放出潜热，使得土壤温度升高，减缓了土壤的冻结速度；而在融化过程中，冰融化吸收潜热，导致土壤温度降低，减缓了融化速度。这种水热耦合作用使得土壤冻融过程更加复杂，对陆面模式中土壤水热传输的模拟提出了更高的要求。2.3.2对地表能量平衡的影响冻融界面变化对地表反照率有着显著影响。地表反照率是指地表反射太阳辐射的能力，它是影响地表能量平衡的重要参数之一。在冻土地区，当土壤处于冻结状态时，地表通常被积雪或冰层覆盖，积雪和冰层具有较高的反照率，一般在0.5-0.9之间。这使得大量的太阳辐射被反射回大气中，减少了地表对太阳辐射的吸收，从而降低了地表获得的能量。在冬季，高反照率的积雪和冰层覆盖使得地表能量输入减少，地表温度降低，进一步维持了土壤的冻结状态。随着气温升高，土壤开始融化，积雪和冰层逐渐减少，地表反照率也随之降低。裸露的土壤和植被的反照率相对较低，一般在0.1-0.3之间。这使得地表对太阳辐射的吸收增加，获得的能量增多，地表温度升高，促进了土壤的进一步融化。在春季，随着积雪和冰层的融化，地表反照率的降低使得地表能量输入增加，加速了冻土的融化过程。冻融界面变化还会对地表热通量产生重要影响。地表热通量包括感热通量和潜热通量，它们是地表与大气之间能量交换的重要方式。在冻结过程中，土壤水分结冰，潜热释放到大气中，使得潜热通量增大。由于土壤冻结导致水力传导率降低，水分蒸发减少，潜热通量的增加幅度相对有限。而感热通量则由于地表温度降低，与大气的温度差减小，使得感热通量减小。这使得地表向大气的能量传输以潜热通量为主，导致大气的湿度增加，温度变化相对较小。在融化过程中，冰融化吸收潜热，使得潜热通量减小。随着土壤水分的增加和蒸发的增强，潜热通量又会逐渐增大。感热通量则由于地表温度升高，与大气的温度差增大，使得感热通量增大。这使得地表向大气的能量传输以感热通量和潜热通量共同作用为主，导致大气的温度和湿度都发生变化。在春季，随着冻土的融化，潜热通量和感热通量的变化使得大气的温度和湿度升高，对大气环流和气候产生影响。地表能量平衡是指地表获得的能量与损失的能量之间的平衡关系。冻融界面变化通过影响地表反照率和热通量，对地表能量平衡产生重要影响。在冻结过程中，地表反照率升高，太阳辐射吸收减少，潜热通量增大，感热通量减小，使得地表能量损失增加，能量平衡向负方向移动，地表温度降低。在融化过程中，地表反照率降低，太阳辐射吸收增加，潜热通量和感热通量的变化使得地表能量输入增加，能量平衡向正方向移动，地表温度升高。地表能量平衡的变化又会反过来影响冻融界面的变化。当地表能量平衡向负方向移动时，地表温度降低，促进土壤的冻结，使冻融界面加深；当地表能量平衡向正方向移动时，地表温度升高，促进土壤的融化，使冻融界面上升。这种相互作用形成了一个复杂的反馈机制，对冻土地区的气候和生态系统产生深远影响。如果气候变暖导致冻土融化，地表反照率降低，地表能量平衡向正方向移动，进一步加速冻土的融化，形成一个正反馈机制，可能导致冻土的快速退化和生态系统的改变。2.3.3对植被生长和生态系统的影响冻融过程对植被根系活动有着重要影响。在冻土地区，土壤的冻结和融化过程会改变土壤的物理性质，如土壤的硬度、孔隙度和水分含量等，这些变化直接影响着植被根系的生长和发育。在冻结期，土壤中的水分结冰，体积膨胀，导致土壤硬度增加，孔隙度减小。这使得植被根系在生长过程中受到更大的阻力，难以深入土壤中获取养分和水分。冻结的土壤还会阻碍根系的呼吸作用，影响根系的正常生理功能。一些研究表明，在冻结期，植被根系的生长速度明显减缓，根系的活力也会下降。在融化期，土壤中的冰逐渐融化，土壤硬度降低，孔隙度增大，水分含量增加。这为植被根系的生长提供了更有利的条件，根系能够更容易地在土壤中伸展和生长，吸收更多的养分和水分。融化期土壤中的微生物活动也会增强，有助于土壤中有机质的分解和养分的释放，进一步促进植被根系的生长。然而，融化期土壤水分的增加也可能导致土壤通气性变差，根系缺氧，对植被生长产生不利影响。冻融过程还会影响植被的生长周期。在冻土地区，气候条件较为恶劣，植被的生长季节相对较短。冻融过程的变化会直接影响植被的生长周期。如果春季气温升高较快，冻土融化提前，植被的生长季节可能会延长，有利于植被的生长和发育。相反，如果春季气温升高缓慢，冻土融化延迟，植被的生长季节可能会缩短，影响植被的生长和产量。冻融过程的异常变化，如频繁的冻融循环，也会对植被的生长周期产生负面影响。频繁的冻融循环会导致土壤结构的破坏，影响土壤的水分和养分供应，使得植被生长受到抑制，生长周期紊乱。植被是生态系统的重要组成部分，冻融过程对植被生长的影响必然会对生态系统的结构和功能产生重要作用。在生态系统结构方面，冻融过程的变化可能导致植被类型和群落结构的改变。随着气候变暖，冻土融化，一些原本适应寒冷环境的植被可能会逐渐减少，而一些适应温暖环境的植被可能会侵入，导致植被类型的更替和群落结构的变化。这种变化会影响生态系统中生物的多样性和物种间的相互关系，进而改变生态系统的结构。在生态系统功能方面，冻融过程的变化会影响生态系统的物质循环和能量流动。植被的生长和代谢过程是生态系统物质循环和能量流动的重要环节。冻融过程对植被生长的影响会直接影响植被的光合作用、呼吸作用和蒸腾作用等生理过程，进而影响生态系统的物质循环和能量流动。在冻结期，植被生长缓慢，光合作用减弱，生态系统中碳的固定和积累减少；而在融化期，植被生长旺盛，光合作用增强，生态系统中碳的固定和积累增加。冻融过程还会影响土壤中微生物的活动和有机质的分解，进一步影响生态系统的物质循环和能量流动。三、考虑冻融界面变化的陆面模式构建3.1现有陆面模式中冻融过程描述的不足在当前众多的陆面模式中，对冻融过程的描述存在着一系列的缺陷，这些不足在一定程度上限制了陆面模式对陆面过程的准确模拟，影响了对气候变化的理解和预测。在参数化方案方面，许多陆面模式对土壤冻融过程中的关键参数化方案不够完善。土壤导热率是描述土壤热量传输能力的重要参数，其准确计算对于模拟土壤温度变化至关重要。在实际的土壤冻融过程中，土壤的导热率会随着土壤含水量、孔隙度以及冰含量的变化而显著改变。在冻结过程中，土壤中的水分结冰，冰的导热率远高于水和空气，使得土壤的导热率增大；而在融化过程中，冰融化成水，土壤导热率又会发生相应的变化。当前大部分数值模式中土壤导热率采用Johansen方案及其派生方案进行计算，虽然这些方案在一定程度上能够描述土壤导热率的变化，但对于土壤冻融过程中复杂的物理变化考虑不足。例如，大部分方案没有充分考虑砾石和有机质对土壤导热率的影响。相对于普通的矿物质土粒而言，砾石具有高导热率和低热容，有机质具有低导热率和高热容，它们对热量在土壤中的传输及土壤温度垂直分布有不同的影响。已有研究表明，Balland–Arp方案考虑了砾石和有机质对土壤导热率的影响，能更好地刻画土壤冻融过程中土壤导热率变化的连续性，但该方案在实际应用中尚未得到广泛采用。土壤水力参数的计算也存在类似问题。土壤的水力传导率和持水率等参数对于模拟土壤水分的运动和分布起着关键作用。在冻融过程中，土壤孔隙结构的变化会导致水力参数发生显著改变。土壤冻结时，孔隙被冰填充，水力传导率大幅降低，水分难以在土壤中传输；而融化时，孔隙结构恢复，水力传导率增加。目前的陆面模式在计算这些参数时，往往没有充分考虑冻融过程中孔隙结构的动态变化，导致对土壤水分传输的模拟存在偏差。一些模式采用的导水率模型没有考虑土壤冻结对水分下渗的阻抗作用，使得在冻土区的水文过程模拟效果不佳。虽然一些改进的方案，如导水阻抗方案考虑了土壤冻结对土壤水分下渗的阻抗作用，改善了对冻土区水文过程的模拟效果，但这些方案在不同陆面模式中的应用还不够广泛，且在实际应用中仍存在一些问题需要进一步解决。在过程描述方面，现有陆面模式对冻融过程的物理机制描述不够全面和准确。土壤冻融过程伴随着水分的相变和能量的转化，这是一个复杂的水热耦合过程。目前的陆面模式在描述这一过程时，往往将水分和热量的传输视为独立的过程，没有充分考虑它们之间的相互作用。在冻结过程中，水分结冰会释放潜热，这部分潜热会影响土壤的温度分布，进而影响冻融界面的移动；而在融化过程中，冰融化吸收潜热，同样会对土壤的温度和水分分布产生影响。许多陆面模式没有准确地刻画这种水热耦合机制，导致对冻融过程的模拟与实际情况存在较大差异。现有陆面模式对等温框架的处理也存在问题。在等温框架的数值模式中，通常通过模拟每层土壤中间深度的冻融过程代表该模式分层的整体特征，这种处理方式导致对冻融深度的严重高估或低估，尤其是对厚度较大的模式深层土壤。实际的冻融锋面并不是水平均匀分布的，而是受到土壤质地、水分含量、温度梯度等多种因素的影响，呈现出复杂的形态。现有陆面模式没有充分考虑这些因素对冻融锋面的影响，使得对冻融深度的模拟存在较大误差。虽然冻融锋面计算方案的提出和应用在一定程度上减小了这种模拟偏差，但在实际应用中，该方案仍需要进一步优化和完善，以提高对冻融深度的模拟精度。现有陆面模式对土壤盐分导致土壤水冰点降低这一事实考虑不足。在一些盐碱地地区，土壤中含有大量的盐分，这些盐分的存在会降低土壤水的冰点，使得土壤在低于0℃时仍可能保持液态。然而，绝大多数数值模式没有考虑这一因素，导致在模拟这些地区的冻融过程时出现偏差。虽然大部分数值模式考虑了土壤有机质对土壤水热传输的影响，但是模式中对土壤有机质含量及垂直分布的考虑与植被根系的生长状态脱节。植被根系的生长会影响土壤有机质的分布和分解，进而影响土壤的水热传输过程。现有陆面模式没有充分考虑这种相互关系，使得对土壤水热传输的模拟不够准确。现有陆面模式模拟的土壤深度不足并且下边界通量为零的假定不符合实际情况。在实际的土壤中，深层土壤对浅层土壤的热力学状态有着重要影响，且下边界存在着水分和热量的交换。现有陆面模式往往忽略了这些因素，导致对土壤水热传输的模拟不够全面和准确。现有陆面模式在描述冻融过程时存在诸多不足，需要进一步改进和完善参数化方案，更全面准确地描述冻融过程的物理机制，以提高陆面模式对冻融界面变化及其对陆面过程影响的模拟能力。3.2改进的冻融过程参数化方案设计3.2.1基于物理过程的参数化方法为了更准确地描述冻融过程，本研究从土壤水热物理特性出发，提出了一套全新的基于物理过程的参数化公式和算法。在土壤导热率计算方面，充分考虑了土壤中砾石和有机质对导热率的影响。相较于普通的矿物质土粒，砾石具有高导热率和低热容，有机质具有低导热率和高热容，它们对热量在土壤中的传输及土壤温度垂直分布有着不同的作用。本研究采用了Balland–Arp方案，该方案能够更好地刻画土壤冻融过程中土壤导热率变化的连续性。具体而言，Balland–Arp方案通过引入砾石含量和有机质含量的参数，建立了与土壤导热率之间的数学关系。假设土壤导热率为\lambda，砾石含量为G，有机质含量为O，通过一系列实验和理论推导，得到如下关系：\lambda=f(G,O)，其中f是一个包含砾石和有机质影响的函数表达式，它考虑了不同成分在土壤中的比例以及它们各自的导热特性，从而能够更准确地反映土壤在冻融过程中导热率的变化。在土壤水力参数计算方面，本研究改进了导水率模型，充分考虑了冻融过程中土壤孔隙结构的动态变化对水力传导率和持水率的影响。土壤冻结时，孔隙被冰填充，水力传导率大幅降低，水分难以在土壤中传输；融化时，孔隙结构恢复，水力传导率增加。为了描述这一过程，引入了一个与孔隙冰含量相关的参数I，建立新的水力传导率K计算公式：K=K_0\timesg(I)，其中K_0是未冻结状态下的水力传导率，g(I)是一个关于孔隙冰含量I的函数，它反映了冰含量对水力传导率的影响程度。当I=0时，g(I)=1，此时水力传导率为未冻结状态下的值；随着I的增加，g(I)逐渐减小，水力传导率也相应降低，从而准确地描述了冻融过程中水力传导率的变化。针对土壤水热耦合过程，本研究建立了考虑水分相变和能量转化的水热耦合模型。土壤冻融过程伴随着水分的相变和能量的转化，是一个复杂的水热耦合过程。在冻结过程中，水分结冰会释放潜热，这部分潜热会影响土壤的温度分布，进而影响冻融界面的移动；融化过程中，冰融化吸收潜热，同样会对土壤的温度和水分分布产生影响。为了描述这一过程，建立了包含能量守恒方程和水分守恒方程的水热耦合模型。能量守恒方程考虑了土壤中热量的传导、对流以及水分相变所释放或吸收的潜热；水分守恒方程则考虑了土壤水分的入渗、蒸发、蒸腾以及在冻融过程中的相变。通过联立求解这两个方程，能够更准确地模拟土壤冻融过程中的水热耦合现象，揭示水热相互作用的内在机制。3.2.2考虑空间异质性的参数化策略不同地形、土壤类型等因素会导致冻融过程的空间异质性显著，为了更准确地模拟这种空间差异，本研究设计了一套能够反映空间异质性的冻融参数化方案。对于地形因素，考虑到海拔高度、坡度和坡向对冻融过程的影响。海拔高度决定了气温和降水的分布，随着海拔的升高，气温降低，降水增加，这会导致高海拔地区的土壤更容易冻结，冻融界面更深。坡度和坡向则影响太阳辐射的接收和地表径流的分布。阳坡接收的太阳辐射多，温度较高，土壤冻结程度相对较轻，冻融界面较浅；阴坡接收的太阳辐射少，温度较低，土壤冻结程度相对较重，冻融界面较深。坡度较大的地区，地表径流速度快，土壤水分容易流失，冻结时冰的形成量相对较少，冻融界面变化相对较快；坡度较小的地区，土壤水分相对稳定，冻融界面变化相对较慢。为了考虑这些因素，在参数化方案中引入了海拔高度H、坡度S和坡向A等参数。根据不同的地形条件，对土壤导热率、水力传导率等参数进行修正。对于高海拔地区，适当调整土壤导热率和热容量参数，以反映低温环境下土壤的热物理性质变化；对于不同坡度和坡向的区域，调整太阳辐射的吸收和地表径流的参数，从而更准确地模拟不同地形条件下的冻融过程。在土壤类型方面，不同质地的土壤，如砂土、壤土和黏土，其孔隙大小、分布和连通性不同，这会影响土壤中水分的存储和传输，进而影响冻融过程。砂土的孔隙较大，通气性和透水性好，但保水性较差，土壤中的水分容易流失，在冻结时冰的形成量相对较少，冻融界面变化相对较快；黏土的孔隙较小，保水性强，但通气性和透水性较差，土壤中的水分不易排出，冻结时冰的形成量较大，冻融界面变化相对较慢；壤土的性质介于砂土和黏土之间，其冻融界面变化的特点也相对适中。为了考虑土壤类型的影响，在参数化方案中对不同土壤类型的水力传导率、持水率等参数进行了差异化设置。根据土壤质地分类，建立不同土壤类型的参数库，在模拟过程中根据实际土壤类型选择相应的参数，以准确反映不同土壤类型下的冻融过程特征。对于砂土，设置较高的水力传导率和较低的持水率；对于黏土，设置较低的水力传导率和较高的持水率；对于壤土，参数取值介于两者之间。本研究还考虑了植被覆盖对冻融过程的影响。植被通过蒸腾作用调节土壤水分含量，同时其根系和枯枝落叶层可以改变土壤的物理性质，进而影响冻融过程。植被的蒸腾作用会消耗土壤中的水分，降低土壤的含水量，从而减少土壤冻结时冰的形成量，影响冻融界面的变化；植被的根系可以增加土壤的孔隙度，改善土壤的通气性和透水性，这也会对土壤的冻融过程产生影响；植被的枯枝落叶层可以起到保温和保水的作用，类似于积雪覆盖，能够减缓土壤的冻融速度，使冻融界面变化相对稳定。为了考虑植被覆盖的影响，在参数化方案中引入了植被覆盖度C、叶面积指数LAI等参数。根据植被覆盖度和叶面积指数的不同，调整土壤水分蒸发、热量传输等参数，以反映植被对冻融过程的调节作用。当植被覆盖度较高时，适当降低土壤水分蒸发速率，增加土壤的保温效果，从而减缓冻融界面的变化速度。通过以上考虑地形、土壤类型和植被覆盖等因素的空间异质性参数化策略，能够更准确地模拟不同区域的冻融过程，提高陆面模式对复杂下垫面条件下冻融界面变化的模拟能力。3.3模式验证与评估3.3.1验证数据的选取与来源为了全面、准确地验证考虑冻融界面变化的陆面模式的性能，本研究精心选取了多种类型的数据，包括实测数据和遥感数据，这些数据涵盖了不同的时间和空间范围，为模式验证提供了丰富的信息。实测数据主要来源于多个地面观测站，这些观测站分布在不同的气候区域和地形条件下，能够代表不同的陆面状况。在高纬度冻土区，选取了位于西伯利亚的多个观测站，这些站点长期监测土壤温度、湿度、气温、降水等气象要素以及土壤的冻融状态。在青藏高原地区，也设立了多个观测站，该地区地势高，气候寒冷，冻土分布广泛，是研究冻融界面变化的理想区域。这些观测站的数据时间跨度从数年到数十年不等，为研究冻融过程的长期变化提供了数据支持。例如，在某观测站，自2000年开始，就对土壤温度和湿度进行了详细的观测，每天记录多个深度的土壤温度和湿度数据，为验证模式在不同季节和年际变化中的模拟能力提供了有力的数据支撑。遥感数据则主要来源于卫星观测，具有覆盖范围广、时间分辨率高的特点，能够提供大面积的陆面信息。利用Landsat系列卫星的遥感影像，获取地表反照率、植被覆盖度等信息。Landsat卫星每隔16天对同一地区进行一次观测，能够提供长时间序列的地表信息，有助于分析地表参数的变化趋势。还利用了MODIS（Moderate-ResolutionImagingSpectroradiometer）卫星数据，获取地表温度、积雪覆盖等信息。MODIS卫星每天能够对全球大部分地区进行观测，其数据时间分辨率高，能够及时反映地表状况的变化。通过对这些遥感数据的处理和分析，可以得到不同区域的地表参数，用于与模式模拟结果进行对比验证。除了地面观测站和卫星遥感数据，还收集了一些其他来源的数据，如再分析数据。再分析数据是通过将观测数据与数值模式相结合，经过同化处理得到的，能够提供全球范围内的气象要素和陆面参数。欧洲中期天气预报中心（ECMWF）的ERA-Interim再分析数据，包含了全球的气温、降水、土壤湿度等信息，这些数据可以作为模式验证的参考，用于评估模式在大尺度空间上的模拟能力。在空间范围上，验证数据覆盖了全球多个冻土区和季节性冻土区，包括北半球的高纬度地区、青藏高原、阿尔卑斯山脉等。这些区域的冻融过程具有不同的特点，通过对这些区域的数据验证，可以全面评估模式在不同地理条件下的性能。在时间范围上，数据涵盖了多年的观测，能够反映冻融过程的年际变化和长期趋势。从2000年到2020年的观测数据，用于分析模式在不同年份和季节的模拟准确性，从而评估模式对冻融界面变化的长期模拟能力。通过多源数据的整合，构建了一个全面、系统的验证数据集，为准确评估考虑冻融界面变化的陆面模式提供了坚实的数据基础。这些数据的合理选取和有效利用，能够充分检验模式在不同条件下的模拟能力，为模式的改进和优化提供有力的支持。3.3.2验证指标与方法为了科学、客观地评估模式模拟结果的准确性，本研究选取了土壤温度、湿度等关键变量，并运用多种统计分析方法进行验证。在土壤温度方面，选用均方根误差（RMSE）、平均绝对误差（MAE）和相关系数（R）作为主要验证指标。均方根误差能够反映模拟值与观测值之间的偏差程度，其计算公式为：RMSE=\sqrt{\frac{1}{n}\sum_{i=1}^{n}(T_{sim,i}-T_{obs,i})^2}，其中n为样本数量，T_{sim,i}为第i个模拟值，T_{obs,i}为第i个观测值。均方根误差的值越小，说明模拟值与观测值越接近，模式的模拟精度越高。平均绝对误差则是模拟值与观测值偏差的绝对值的平均值，其计算公式为：MAE=\frac{1}{n}\sum_{i=1}^{n}|T_{sim,i}-T_{obs,i}|，它能够直观地反映模拟值与观测值之间的平均偏差大小。相关系数用于衡量模拟值与观测值之间的线性相关性，其取值范围在-1到1之间，R越接近1，表示两者之间的线性相关性越强，模式的模拟结果与观测数据的变化趋势越一致。在土壤湿度方面，同样采用均方根误差、平均绝对误差和相关系数作为验证指标。其计算公式与土壤温度的计算类似，只是将变量T替换为土壤湿度S。除了这些指标，还引入了偏差（Bias）来评估模式对土壤湿度的模拟偏差情况，偏差的计算公式为：Bias=\frac{1}{n}\sum_{i=1}^{n}(S_{sim,i}-S_{obs,i})，偏差为正值表示模拟值大于观测值，偏差为负值表示模拟值小于观测值。在评估过程中，采用了多种统计分析方法。将模式模拟结果与实测数据进行逐点对比，计算各项验证指标，以直观地了解模式在每个观测点的模拟准确性。还对不同区域和不同时间尺度的数据进行统计分析，例如，将研究区域划分为不同的子区域，分别计算每个子区域的验证指标，以分析模式在不同地理条件下的模拟性能差异；对不同季节和年份的数据进行统计分析，以评估模式在不同时间尺度上的模拟稳定性。为了进一步评估模式对冻融界面变化的模拟能力，还采用了空间分布对比和时间序列分析的方法。在空间分布对比中，将模式模拟的冻融界面位置与实际观测的冻融界面位置进行对比，分析模式对冻融界面空间分布的模拟准确性。通过绘制模拟和观测的冻融界面分布图，直观地展示两者之间的差异，并计算空间相关系数等指标来定量评估差异程度。在时间序列分析中，对比模式模拟的冻融过程时间序列与实际观测的时间序列，分析模式对冻融过程的时间变化特征的模拟能力。通过绘制时间序列图，观察模拟和观测的冻融开始时间、结束时间以及冻融持续时间等特征的一致性，并计算时间相关系数等指标来评估模拟的准确性。通过综合运用上述验证指标和方法，能够全面、深入地评估考虑冻融界面变化的陆面模式在模拟土壤温度、湿度以及冻融界面变化等方面的准确性和可靠性，为模式的改进和优化提供科学依据。3.3.3结果分析与讨论对比改进前后陆面模式的模拟结果，发现改进后的模式在模拟冻融界面变化及相关陆面过程方面展现出了显著的优势，同时也存在一些有待完善的地方。在土壤温度模拟方面，改进后的模式在精度上有了明显提升。以某高纬度冻土区观测站的数据为例，改进前的模式模拟的土壤温度均方根误差（RMSE）在10cm深度处为3.5℃，改进后降低至2.0℃；平均绝对误差（MAE）从2.8℃减小到1.5℃；相关系数（R）从0.75提高到0.85。这表明改进后的模式能够更准确地模拟土壤温度的变化，与实际观测值更为接近。从空间分布上看，改进后的模式对土壤温度的空间变化特征模拟得更加准确。在青藏高原地区，改进前的模式在模拟土壤温度时，存在明显的冷偏差，尤其是在冬季，模拟的土壤温度比实际观测值偏低较多；而改进后的模式能够较好地捕捉到土壤温度的空间差异，减少了冷偏差现象，使模拟结果更符合实际情况。在土壤湿度模拟方面，改进后的模式同样表现出色。在某季节性冻土区，改进前模式模拟的土壤湿度均方根误差在表层土壤为0.12m³/m³，改进后降低至0.08m³/m³；平均绝对误差从0.09m³/m³减小到0.06m³/m³；相关系数从0.7提高到0.8。偏差也得到了有效改善，改进前模式模拟的土壤湿度存在正偏差，即模拟值普遍大于观测值，而改进后的模式偏差明显减小，更接近实际观测值。改进后的模式对土壤湿度在冻融过程中的动态变化模拟得更加准确。在土壤冻结和融化过程中，改进前的模式对土壤湿度的变化趋势模拟不够准确，导致模拟的土壤湿度与实际观测值存在较大偏差；而改进后的模式能够较好地模拟土壤湿度在冻融过程中的变化，准确地反映出水分在土壤中的迁移和相变过程。在冻融界面变化模拟方面，改进后的模式取得了显著的进步。改进后的模式能够更准确地模拟冻融界面的位置和移动过程。通过与实际观测的冻融界面位置对比，发现改进前的模式对冻融界面深度的模拟存在较大误差，尤其是在深层土壤，误差可达20-30cm；而改进后的模式将误差控制在10cm以内，大大提高了模拟精度。改进后的模式对冻融界面变化的时间特征模拟也更加准确。在模拟冻融开始时间和结束时间时，改进前的模式与实际观测值的偏差可达1-2周，而改进后的模式偏差缩短至3-5天，能够更及时地反映冻融过程的时间变化。改进后的模式仍存在一些不足。在极端气候条件下，如强降雪或快速升温事件，模式的模拟能力有待提高。在一次强降雪事件中，改进后的模式虽然能够模拟出土壤温度和湿度的变化趋势，但在变化幅度上与实际观测值存在一定偏差，对冻融界面的响应速度也略显滞后。在复杂地形和下垫面条件下，模式的模拟精度也会受到影响。在山区，由于地形复杂，土壤质地和植被覆盖变化较大，模式在模拟土壤水热传输和冻融界面变化时，仍存在一定的误差，需要进一步改进参数化方案以适应复杂的地形和下垫面条件。改进后的考虑冻融界面变化的陆面模式在模拟土壤温度、湿度以及冻融界面变化等方面取得了显著的进展，能够更准确地反映陆面过程的实际情况，但在极端气候条件和复杂地形下仍需进一步优化和完善，以提高其模拟能力和可靠性。四、陆面模式与CAS-ESM耦合机制4.1CAS-ESM系统概述CAS-ESM作为中国科学院研发的地球系统模式，采用模块化框架，集成了多个分系统模式，能够全面地模拟地球系统的各个组成部分及其相互作用。其系统结构涵盖了大气、海洋、陆面、海冰、气溶胶和大气化学、植被动力学、陆地生物地球化学以及海洋生物地球化学等多个关键领域，这些分系统模式通过耦合器紧密连接，实现了地球各圈层间物质和能量的交换，共同构成了一个复杂而庞大的地球系统模拟体系。大气环流模式是CAS-ESM的重要组成部分，它主要用于模拟大气的运动和变化。该模式通过求解大气动力学和热力学方程组，来描述大气中的各种物理过程，如大气的水平和垂直运动、热量传输、水汽输送等。在模拟大气运动时，考虑了地球自转、地形地貌、太阳辐射等多种因素的影响。地球自转产生的科里奥利力会改变大气的运动方向，使得大气在北半球向右偏转，在南半球向左偏转；地形地貌的起伏会影响大气的流动，山脉会阻挡气流，导致气流上升或绕流，从而影响降水和气温的分布；太阳辐射是大气运动的主要能源，不同地区接收到的太阳辐射强度不同，导致大气温度和气压的差异，进而驱动大气的运动。通过对这些因素的综合考虑，大气环流模式能够准确地模拟大气环流的形成和演变，为研究气候变化提供重要的基础。海洋环流和海冰模式则专注于模拟海洋的环流和海冰的变化。海洋环流是海洋中大规模的水体运动，它对全球气候有着重要的调节作用。该模式通过考虑海水的温度、盐度、密度等因素，以及风应力、热通量、淡水通量等强迫作用，来模拟海洋环流的形成和变化。在模拟海洋环流时，考虑了海洋中不同深度的水温差异，水温较高的海水密度较小，会向上流动，而水温较低的海水密度较大，会向下流动，这种密度差异驱动了海洋环流的形成。海冰模式则主要描述海冰的生长、融化、漂移等过程，以及海冰与海洋和大气之间的相互作用。海冰的存在会影响海洋与大气之间的热量和水分交换，海冰的高反照率会反射大量的太阳辐射，减少海洋对太阳辐射的吸收，从而影响海洋的温度和环流。陆面及水文过程模式是描述陆面与大气之间物质和能量交换的关键模式，包括土壤水分、热量传输，植被生长和蒸散等过程。在土壤水分传输方面，考虑了降水、蒸发、下渗、地表径流等因素，通过建立水分平衡方程来模拟土壤水分的动态变化。在热量传输方面，考虑了太阳辐射、地表反照率、感热通量、潜热通量等因素，通过能量守恒方程来模拟土壤温度的变化。植被生长和蒸散过程则考虑了植被的生理特性、光照、温度、水分等因素，通过建立植被生长模型和蒸散模型来模拟植被的生长和水分消耗。该模式对研究陆面过程和气候变化具有重要意义，它能够为大气环流模式提供准确的下垫面条件，同时也能够反映气候变化对陆面生态系统的影响。植被动力学模式主要用于模拟植被的动态变化，包括植被的生长、死亡、演替等过程。该模式考虑了植被的生理特性、环境因素（如温度、水分、光照、土壤养分等）以及人类活动的影响，通过建立植被生长模型和生态系统模型来模拟植被的动态变化。在模拟植被生长时，考虑了植被的光合作用、呼吸作用、蒸腾作用等生理过程，以及这些过程与环境因素之间的相互关系。植被动力学模式能够反映植被对气候变化的响应，以及植被变化对生态系统和大气环境的影响，对于研究生态系统的稳定性和可持续性具有重要意义。气溶胶和大气化学模式用于模拟大气中的气溶胶和化学物质的分布和变化，以及它们与大气物理过程的相互作用。气溶胶是悬浮在大气中的固态或液态颗粒物，它们对大气的辐射平衡、云的形成和降水等过程都有着重要的影响。该模式考虑了气溶胶的来源、传输、转化和清除等过程，以及气溶胶与大气化学物质（如二氧化硫、氮氧化物、挥发性有机物等）之间的化学反应。通过模拟这些过程，能够了解气溶胶和大气化学物质对气候变化和空气质量的影响，为制定环境保护政策提供科学依据。陆地生物地球化学和海洋生物地球化学模式分别描述陆地和海洋生态系统中碳、氮、磷等元素的循环过程。在陆地生物地球化学模式中，考虑了植被的光合作用、呼吸作用、凋落物分解、土壤微生物活动等过程，以及这些过程对碳、氮、磷等元素循环的影响。在海洋生物地球化学模式中，考虑了海洋浮游生物的生长、死亡、分解，以及海洋中化学物质的溶解、沉淀、吸附等过程，以及这些过程对海洋生态系统中碳、氮、磷等元素循环的影响。这些模式对于研究全球碳循环和气候变化具有重要意义，它们能够揭示生态系统在调节气候变化中的作用，为预测未来气候变化提供重要的参考。这些分系统模式在CAS-ESM中相互关联、协同工作。大气环流模式为海洋环流和海冰模式提供大气强迫，如风速、气温、降水等；海洋环流和海冰模式则为大气环流模式提供海洋表面的温度、湿度和海冰覆盖等边界条件。陆面及水文过程模式与大气环流模式之间存在着密切的物质和能量交换，陆面的蒸发和蒸腾为大气提供水汽，而大气的降水和热量又影响着陆面的水分和热量平衡。植被动力学模式通过影响陆面的反照率、粗糙度和蒸散等参数，对大气环流和陆面过程产生影响；同时，大气环流和陆面过程的变化也会影响植被的生长和分布。气溶胶和大气化学模式与大气环流模式相互作用，气溶胶和大气化学物质的分布会影响大气的辐射平衡和物理过程，而大气环流的变化又会影响气溶胶和大气化学物质的传输和扩散。陆地生物地球化学和海洋生物地球化学模式与其他分系统模式之间也存在着复杂的相互关系，它们通过碳、氮、磷等元素的循环，影响着大气和海洋的化学成分，进而影响气候变化。CAS-ESM通过集成多个分系统模式，构建了一个全面、复杂的地球系统模拟体系，能够深入研究地球系统各圈层之间的相互作用和气候变化的机制，为全球气候变化研究提供了强大的工具。4.2耦合方案设计4.2.1耦合方式与接口选择陆面模式与CAS-ESM的耦合方式主要有单向耦合和双向耦合两种，它们各自具有不同的特点和适用场景，在耦合过程中，接口的选择至关重要，它直接影响着数据交换的效率和准确性。单向耦合是一种较为简单的耦合方式，在这种方式下，陆面模式仅接收来自CAS-ESM中大气模式的驱动数据，如气温、降水、辐射等，然后根据这些数据进行陆面过程的模拟，但模拟结果不会反馈给大气模式。这种耦合方式的优点是计算相对简单，计算效率较高，因为不需要考虑陆面模式对大气模式的反馈作用，减少了计算量。在一些对陆面过程反馈作用不太敏感的研究中，单向耦合可以快速地得到陆面过程的模拟结果，为研究提供初步的数据支持。单向耦合也存在明显的局限性，由于忽略了陆面过程对大气的反馈，可能会导致模拟结果与实际情况存在一定偏差。陆面的蒸发和蒸腾会影响大气的湿度和能量平衡，如果不考虑这些反馈，可能会使大气模式对降水和气温的模拟出现误差。双向耦合则充分考虑了陆面模式与大气模式之间的相互作用。在双向耦合中，陆面模式不仅接收大气模式的驱动数据，其模拟结果，如地表反照率、感热通量、潜热通量等，也会实时反馈给大气模式，作为大气模式下一时间步的边界条件。这种耦合方式能够更真实地反映地球系统中陆面与大气之间的复杂相互作用，提高模拟结果的准确性。在研究气候变化对陆面生态系统的影响时，双向耦合可以考虑到陆面生态系统变化对大气的反馈，从而更准确地预测气候变化的趋势。双向耦合的计算过程相对复杂，需要在陆面模式和大气模式之间频繁地进行数据交换和计算，对计算资源的要求较高。由于陆面过程和大气过程的时间尺度不同，如何合理地协调两者之间的数据交换频率和计算步长，是双向耦合中需要解决的关键问题。在接口选择方面，数据交换接口的设计是实现陆面模式与CAS-ESM有效耦合的关键环节。常用的接口包括基于文件的接口和基于内存共享的接口。基于文件的接口是一种较为传统的接口方式，它通过文件系统来实现陆面模式与CAS-ESM之间的数据交换。在这种接口方式下，陆面模式将模拟结果输出到文件中，大气模式从文件中读取这些数据作为边界条件。这种接口的优点是实现相对简单，兼容性好，几乎所有的陆面模式和大气模式都可以通过文件系统进行数据交换。它也存在一些缺点，如数据交换速度较慢，容易受到文件系统读写速度的限制；文件的存储和管理也会占用一定的存储空间和系统资源。基于内存共享的接口则是一种更为高效的接口方式，它利用计算机内存的共享机制，实现陆面模式与CAS-ESM之间的数据直接交换。在这种接口方式下，陆面模式和大气模式可以直接访问共享内存中的数据，无需通过文件系统进行读写操作，大大提高了数据交换的速度和效率。基于内存共享的接口还可以减少数据的冗余存储，提高系统的整体性能。这种接口的实现相对复杂，需要对陆面模式和大气模式的代码进行一定的修改，以支持内存共享机制；它对计算机硬件和操作系统的要求也较高，需要具备足够的内存和支持内存共享的操作系统。在实际应用中，需要根据具体的研究需求和计算资源情况，选择合适的耦合方式和接口。对于一些对计算效率要求较高、对陆面过程反馈作用不太敏感的研究，可以采用单向耦合和基于文件的接口；而对于一些对模拟准确性要求较高、需要考虑陆面与大气相互作用的研究，则应采用双向耦合和基于内存共享的接口。还可以结合不同的接口方式，发挥它们各自的优势，以实现陆面模式与CAS-ESM的高效、准确耦合。4.2.2数据传输与共享机制陆面模式与CAS-ESM之间的数据传输频率、格式和共享策略对于实现准确的地球系统模拟至关重要，它们直接影响着耦合系统的稳定性和模拟结果的准确性。数据传输频率是一个关键参数，它决定了陆面模式和CAS-ESM之间数据更新的及时性。不同的陆面过程和大气过程具有不同的时间尺度，因此需要根据具体情况合理确定数据传输频率。对于大气过程，其变化较为迅速，时间尺度较短，通常需要较高的数据传输频率。大气中的对流活动、降水过程等在短时间内会发生显著变化，因此大气模式可能需要每小时甚至更短时间接收一次陆面模式传来的数据，以准确反映陆面过程对大气的影响。而陆面过程相对较为缓慢，时间尺度较长，例如土壤水分的变化、植被的生长等过程在数天甚至数周内才会有明显变化，因此陆面模式向大气模式传输数据的频率可以相对较低，一般可以每天或数天传输一次。在实际应用中，为了协调不同时间尺度的过程，常采用时间步长嵌套的方法。大气模式采用较短的时间步长进行计算，陆面模式采用较长的时间步长进行计算。在陆面模式计算的每个时间步内，根据大气模式的时间步长，将陆面模式的输出数据进行插值，以提供给大气模式在不同时间步的边界条件。这样既保证了大气模式能够及时获取陆面过程的信息，又避免了因数据传输过于频繁而增加计算负担。数据格式的选择也十分重要，它直接影响数据的传输效率和准确性。常用的数据格式包括二进制格式和NetCDF（NetworkCommonDataForm）格式。二进制格式是一种高效的数据存储和传输格式，它以二进制编码的方式存储数据，占用空间小，读写速度快。由于二进制格式是一种机器特定的格式，不同的计算机系统可能对其解析方式不同，因此在跨平台数据交换时可能会出现兼容性问题。NetCDF格式则是一种自描述的、与平台无关的数据格式，它能够方便地存储和管理多维数组数据，并且提供了丰富的元数据信息，用于描述数据的属性、单位、坐标等。NetCDF格式具有良好的兼容性和可扩展性，被广泛应用于地球科学领域的数据存储和交换。在陆面模式与CAS-ESM的耦合中，采用NetCDF格式可以确保数据在不同模式之间的准确传输和共享，同时便于数据的管理和分析。数据共享策略是确保陆面模式与CAS-ESM之间数据有效传输和利用的重要保障。为了实现高效的数据共享，常采用分布式存储和并行计算的策略。分布式存储是将数据分散存储在多个存储节点上，通过网络进行数据访问和传输。这种方式可以提高数据的存储容量和访问速度，同时增强数据的可靠性和可扩展性。在陆面模式与CAS-ESM耦合系统中，将陆面模式和大气模式的输入输出数据分别存储在不同的存储节点上，通过高速网络实现数据的共享和交换。并行计算则是利用多处理器或多核处理器的计算能力，同时进行多个任务的计算。在耦合系统中，陆面模式和大气模式可以分别在不同的处理器上进行并行计算，提高计算效率。通过合理的任务分配和数据通信机制，实现陆面模式和大气模式之间的协同计算，确保数据的及时传输和共享。为了保证数据的质量和一致性，还需要建立数据质量控制和验证机制。在数据传输过程中，对数据进行质量检查，确保数据的准确性和完整性。对于异常数据或错误数据，及时进行处理和纠正。定期对耦合系统的模拟结果进行验证，与实际观测数据进行对比分析，评估模拟结果的准确性，根据验证结果对数据传输和共享策略进行调整和优化。陆面模式与CAS-ESM之间的数据传输与共享机制是一个复杂的系统工程，需要综合考虑数据传输频率、格式和共享策略等多个因素，以实现高效、准确的数据交换和共享，为地球系统模拟提供可靠的数据支持。4.3耦合过程中的关键技术问题及解决方法在陆面模式与CAS-ESM的耦合过程中，不可避免地会遇到一系列关键技术问题，这些问题若不妥善解决，将严重影响耦合系统的稳定性和模拟结果的准确性。能量不平衡是耦合过程中常见的问题之一。陆面模式与CAS-ESM中的大气模式在能量计算和传输过程中，由于各自的参数化方案、物理过程描述以及计算精度的差异，可能会导致能量在耦合界面处出现不平衡的情况。这种能量不平衡会随着时间的推移不断积累，进而影响整个耦合系统的稳定性和模拟结果的可靠性。在模拟地表能量平衡时，陆面模式计算的感热通量和潜热通量与大气模式接收的能量可能存在差异，这可能是由于陆面模式对植被蒸腾和土壤水分蒸发的参数化方案与大气模式对水汽输送和能量吸收的处理方式不一致所导致的。为了解决能量不平衡问题，采用能量再分配的方法进行调整。通过在耦合界面处建立能量平衡方程，对陆面模式和大气模式之间的能量通量进行实时监测和计算。当发现能量不平衡时，根据能量守恒定律，对能量通量进行合理的再分配。可以将多余的能量从能量通量较大的一方转移到能量通量较小的一方，以实现能量的平衡。为了确保能量再分配的合理性和准确性，还可以引入反馈机制，根据能量不平衡的程度和方向，动态调整再分配的比例和方式。如果能量不平衡主要是由于感热通量的差异导致的，则可以适当增加感热通量较小一方的感热通量，同时减少感热通量较大一方的感热通量，以达到能量平衡的目的。时间步长不一致也是耦合过程中需要解决的重要问题。陆面模式和CAS-ESM中的其他分系统模式，如大气模式、海洋模式等，由于各自所描述的物理过程的时间尺度不同，通常采用不同的时间步长进行计算。大气模式需要较短的时间步长来准确模拟快速变化的大气运动和天气过程，而陆面模式则可以采用相对较长的时间步长来模拟相对缓慢的陆面过程。这种时间步长的不一致会给耦合过程带来困难，可能导致数据传输和计算的不协调，影响模拟结果的准确性。为了解决时间步长不一致的问题，采用时间步长嵌套的方法。以大气模式的时间步长为基准，将陆面模式的计算结果在大气模式的每个时间步内进行插值，以提供与大气模式时间步长相匹配的数据。在陆面模式计算的每一个时间步内，根据大气模式的时间步长，将陆面模式输出的土壤温度、湿度、感热通量、潜热通量等数据进行线性插值或其他合适的插值方法，得到在大气模式每个时间步对应的数值。这样，大气模式在每个时间步都能获取到与陆面模式时间尺度相协调的数据，从而实现陆面模式与大气模式在不同时间步长下的有效耦合。还可以采用多重时间步长算法，即根据不同分系统模式的时间尺度特点，分别采用不同的时间步长进行计算，然后通过适当的算法将这些不同时间步长的计算结果进行整合，以实现整个耦合系统的稳定运行。变量尺度不匹配也是耦合过程中需要关注的问题。陆面模式和CAS-ESM中的其他分系统模式在空间分辨率和垂直分层等方面可能存在差异，这会导致变量在耦合过程中的尺度不匹配。陆面模式的空间分辨率可能较粗，而大气模式的空间分辨率较细，当陆面模式向大气模式传输地表参数时，可能会因为分辨率的差异而导致信息丢失或不准确。陆面模式和大气模式在垂直分层上的差异也会影响变量的传输和耦合效果。为了解决变量尺度不匹配问题，采用空间插值和垂直插值的方法。在空间插值方面，当陆面模式向大气模式传输数据时，根据大气模式的空间分辨率，采用合适的空间插值算法，如双线性插值、样条插值等，将陆面模式的粗分辨率数据插值到大气模式的细分辨率网格上，以保证数据在空间尺度上的一致性。在垂直插值方面，根据陆面模式和大气模式的垂直分层结构，采用线性插值或其他合适的方法，对变量在垂直方向上进行插值，以实现垂直尺度的匹配。还可以通过建立变量转换关系，对不同尺度下的变量进行转换和调整，以确保变量在耦合过程中的兼容性和准确性。通过针对能量不平衡、时间步长不一致和变量尺度不匹配等关键技术问题采取相应的解决方法，可以有效提高陆面模式与CAS-ESM耦合系统的稳定性和模拟结果的准确性，为地球系统模拟研究提供更可靠的工具。五、耦合模式的应用与案例分析5.1案例选取与模拟设置为了深入评估考虑冻融界面变化的陆面模式与CAS-ESM耦合模式的性能，本研究选取了青藏高原地区作为典型案例区域。青藏高原被誉为“世界屋脊”和“亚洲水塔”，是全球海拔最高、面积最大的高原，其独特的地理环境和气候条件使其成为研究冻融界面变化及其对陆面过程影响的理想区域。该地区冻土分布广泛，约占中国冻土总面积的70%，冻融过程复杂且对气候变化极为敏感。在时间跨度上，选择了2000-2010年这一时间段进行模拟研究。这一时期涵盖了不同的气候条件和季节变化，能够全面反映耦合模式在不同时间尺度上对冻融界面变化和陆面过程的模拟能力。在这十年间，青藏高原地区经历了气温的波动上升、降水模式的变化以及冻土的动态演变，为研究耦合模式对气候变化的响应提供了丰富的素材。模拟设置方面，在耦合模式中，陆面模式的水平分辨率设置为0.5°×0.5°，这种分辨率能够较好地捕捉青藏高原地区复杂的地形和下垫面特征，准确模拟冻融界面变化及其对陆面过程的影响。大气模式的水平分辨率设置为1.0°×1.0°，该分辨率在保证计算效率的同时，能够合理地模拟大气环流的变化，为陆面模式提供准确的大气驱动数据。海洋模式的水平分辨率设置为1/3°×1/3°，能够精确地模拟海洋环流和海温的变化，考虑海洋对气候的调节作用。垂直方向上，陆面模式分为10层，能够详细地描述土壤温度和湿度在不同深度的变化；大气模式分为30层，能够准确地模拟大气温度、湿度和气压在垂直方向上的分布；海洋模式分为40层，能够全面地反映海洋温度、盐度和密度在不同深度的变化。边界条件的设置对于模拟结果的准确性至关重要。在模拟过程中，大气边界条件采用NCEP/NCAR再分析资料，该资料提供了全球范围内的大气温度、湿度、气压、