基于水热变化的青藏高原土壤冻融过程研究进展

青藏高原近地层土壤冻融过程是高原地表最显著的陆面特征之一,也是判断冻土发育、存在以及反映气候变化的重要指标。近地层土壤昼夜、季节性的冻结、融化会导致青藏高原陆—气间能水平衡的变化甚至异常,从而显著影响高原地表水文过程、生态环境、碳氮循环以及高原及其周边区域的天气和气候系统。论文从观测、模拟以及对气候的影响3个角度来探讨1990年以来青藏高原土壤冻融过程的最新研究进展。结果表明：① 在一个完整的年冻融循环过程中,近地表各层土壤大体都经历了夏季融化期、春秋季融化—冻结期、冬季冻结期4个阶段。受局地因素的影响,不同站点的冻结或消融起止时间、速率、类型均有差异。② 多年冻土区和季节冻土区的日冻融循环过程差异较大,主要体现在日冻融循环持续时间上。③ 不同陆面模式都可以很好地抓住冻融过程中物理量的时空变化,但都需要针对高原陆面过程的特点进行参数化改进。④ 规避不稳定的迭代计算并根据热力学平衡方程确定冻融临界温度可以改进不合理的冻融参数化方案。基于已有研究回顾,发现增加高质量的观测站,利用卫星遥感等多种手段来反演高原土壤冻融过程以及加强陆面模式与区域气候模式和全球气候模式的耦合,并立足于高原冻融过程的特点发展相适应的参数化方案以及模拟结构的调整,能够有助于高原冻融过程的模拟。     

An improvement of soil temperature simulations on the Tibetan Plateau

The simulation of soil temperature on the Tibetan Plateau (TP) plays a dominant role in the performance of both global climate and numerical weather forecast models. To improve the simulation of soil temperature on the TP, the Johansen soil thermal conductivity parameterization scheme was introduced into Community Land Model 3.5 (CLM3.5) and Regional Climatic Model 4 (RegCM4). The improved CLM3.5 and RegCM4-CLM were utilized to conduct offline and regional simulation experiments on the TP. Comparison of the new and old schemes revealed that CLM3.5 provides high thermal conductivity parameters of mineral soil solid on the TP. The Johansen scheme is more practical for the TP than the soil thermal conductivity parameterization in CLM3.5. The simulation of soil temperature and liquid water content was improved in offline experiment. The improved parameterization scheme can also reduce the simulation error of soil temperature in winter throughout the entire TP.     

Numerical Simulation on Climate Effects of Freezing-Thawing Processes Using CCM3

A parameterization of soil freezing-thawing physics for use in the land-surface model of the National Center for Atmospheric Research(NCAR) Community Climate Model(CCM3) is developed and evaluated.The new parameterization scheme has improved the representation of physical processes in the existing land surface model.Numerical simulations using CCM3 with improved land-surface processes and with the original land-surface processes are compared against the NCEP reanalysis.It is found that the CCM3 version using the improved land surface model shows significant improvements in simulating precipitation in China during the summer season,the general circulation over East Asia,and wind fields over the Tibet Plateau.For the summer season,the improved model was able to better simulate the Indian summer monsoon components,including the mean northerly wind in the upper troposphere and mean southerly wind in the lower troposphere.     

沙漠下垫面地面温度及能量收支的BATS模式参数化改进

利用陆面过程模式BATS,引入地表发射率及两种大气发射率参数化方案,同时引入不同的地表粗糙度参数化方案,对比各种参数化方案对沙漠下垫面地面温度及能量收支的模拟状况。结果表明：采用Van Bavel等发展的地表发射率及Chung等发展的大气发射率方案可以明显改进地面温度及向上长波辐射的模拟,Chung等方案在夜间与正午的模拟效果更好,减小了1 ℃左右的地面温度模拟偏差,减小了10 W·m^-2左右的向上长波辐射模拟偏差。晴天地面温度及向上长波辐射的模拟结果优于阴天。利用Zhang等发展的裸土粗糙度参数化方案也会提高模式对地表感热通量模拟的准确性。     

Wind tunnel simulation of the effects of freeze–thaw cycles on soil erosion in the Qinghai–Tibet Plateau

Intense freezing and thawing actions occur in the Qinghai–Tibet Plateau because of its high elevation and cold temperature. The plateau's unique environment makes it easy to generate wind erosion under dry, windy weather conditions, resulting in the emergence of desertification. As a major form of freeze–thaw erosion, freeze–thaw and wind erosion is displayed prominently on the Qinghai–Tibet Plateau. Therefore, in this study, soil samples were collected from the surface of the plateau to undergo freeze–thaw and wind erosion simulation experiments. Results show that wind erosion strength increases with an increasing number of freeze–thaw cycles, water content in the freezing–thawing process, and the difference in freeze–thaw temperatures. Therefore, in the conditions of water participation, the main reason for the freeze–thaw and wind erosion in the Qinghai–Tibet Plateau is the damage to the soil structure by repeated, fierce freeze–thaw actions, and the sand-bearing wind is the main driving force for this process. The research results have theoretical significance for exploring the formation mechanism of freeze–thaw and wind erosion in the Qinghai–Tibet Plateau, and provide a scientific basis for freeze–thaw desertification control in the plateau.     

Characteristics of GLDAS soil-moisture data on the Tibet Plateau

In this paper, the applicability of soil-moisture (SM) datasets of GLDAS (Global Land Data Assimilation System) in an alpine region (Tibet Plateau, TP) is investigated. The relations and discrepancies between the GLDAS-NOAH SM (0~10 cm) and the observations are compared; the possible reasons for errors over the TP are explored. The results show that GLDAS SM biases mainly show up in errors of values in the nonfrozen period (April to October) and changes of SM along with the temperature, especially during the freezing-thawing process in the frozen period (November to March). The biases of GLDAS SM in the nonfrozen period are mainly caused by the GLDAS precipitation-forcing data. The errors of GLDAS SM in the frozen period are speculated to be induced by the freeze-thaw parameterization scheme in the land-surface model.     

陆面过程参数化模式的研究

在概述了近年来国内外陆面过程参数化研究的基础上,发展了一个改进的陆面过程参数化模式。该模式能被应用于地球表面不同土壤和植被类型的陆面过程的相互作用研究。在该模式中通过土壤结构参数和陆面过程湍流通量的参数化,可广泛应用于不同类型的区域。在模式中,大气-陆面系统包括3个子系统:近地层大气、冠层和土壤活跃层。土壤被分为4层,上两层分别相当于温度的日和季变化层,植被的根系活动也在这两层内,在最深层的温度和含水量是气候的平均值。发展的陆面过程参数化模式可应用于不同气候区域能量平衡和气候的模拟研究。对我国在该领域的研究,特别是在开发我国西北地区过程中的环境生态研究中具有重要的意义。     

模型驱动数据空间分辨率对模拟生态水文过程的影响

设定两种独立的调参方案,检验将高/低分辨率下的调参值应用到低/高分辨率中驱动模型模拟的可行性;同时分析驱动数据分辨率的差异对模型模拟生态水文过程的影响;此外,设置9种集水面积阈值,分析驱动数据分辨率与集水面积阈值对模拟结果的影响。结果表明两种调参方案下模拟的水文与生态数据相关性强、差异小,说明调参方案受驱动数据分辨率的影响不明显;分辨率对土地覆被和土壤类型面积影响小,而对高程和坡度影响大;驱动数据分辨率降低,模拟径流的均值减小、峰值增大、生态变量值增大;驱动数据分辨与集水面积阈值之间不存在交互作用;集水面积阈值增大,模拟径流的均值、峰值均减小,生态变量值也减小,但阈值变化对峰现时间无影响。     

DLM-FvCB模型光合能力关键参数季节性参数化方案及模拟效果评估

选取包含9种植被功能型的15个通量塔观测数据（60站点年）,采用控制变量法,分别在月尺度和年尺度上对25℃时最大羧化速率（V_cmax25）进行优化,并利用未参与参数反演的站点年数据集定量评估这2种参数化方案的模拟效果。研究发现：① 在参数的季节变异方面：9种植被功能型的V_cmax25均呈现明显的季节性波动。V_cmax25的波动幅度为：冬季>秋季>春季>夏季,不同植被类型季节变化幅度相差不大,而寒带植被V_cmax25季节变化幅度接近温带植被的2倍;② 对生态系统初级生产力（GPP）估算的影响方面：V_cmax25季节性参数化方案显著提高了GPP的模拟能力和模拟精度,其中森林和灌木冬季提升最显著（R²提高了35.7%,RMSE降低了23.24%）,春秋季次之,夏季最小。然而,对于C₃草地,DLM无论采取V_cmax25季节性参数化方案还是年尺度参数化方案,均为系统低估GPP。     

Characteristics of abrupt changes of snow cover and seasonal freeze-thaw layer in the Tibetan Plateau and their impacts on summer precipitation in China

In this paper, a variation series of snow cover and seasonal freeze-thaw layer from 1965 to 2004 on the Tibetan Plateau has been established by using the observation data from meteorological stations. The sliding T-test, M-K test and B-G algorithm are used to verify abrupt changes of snow cover and seasonal freeze-thaw layer in the Tibetan plateau. The results show that the snow cover has not undergone an abrupt change, but the seasonal freeze-thaw layer obviously witnessed a rapid degradation in 1987, with the frozen soil depth being reduced by about 15 cm. It is also found that when there is less snow in the plateau region, precipitation in South China and Southwest China increases. But when the frozen soil is deep, precipitation in most of China apparently decreases. Both snow cover and seasonal freeze-thaw layer on the plateau can be used to predict the summer precipitation in China. However, if the impacts of snow cover and seasonal freeze-thaw layer are used at the same time, the predictability of summer precipitation can be significantly improved. The significant correlation zone of snow is located in middle reaches of the Yangtze River covering the Hexi Corridor and northeastern Inner Mongolia, and the seasonal freeze-thaw layer exists in Mt. Nanling, northern Shannxi and northwestern part of North China. The significant correlation zone of simultaneous impacts of snow cover and seasonal freeze-thaw layer is larger than that of either snow cover or seasonal freeze-thaw layer. There are three significant correlation zones extending from north to south: the north zone spreads from Mt. Daxinganling to the Hexi Corridor, crossing northern Mt. Taihang and northern Shannxi; the central zone covers middle and lower reaches of the Yangtze River; and the south zone extends from Mt. Wuyi to Yunnan and Guizhou Plateau through Mt. Nanling.     

天山暴雨数值模拟的试验研究

利用新疆天山地区加密自动气象站资料、常规探空资料和NOAA卫星AMSU辐射亮温资料,应用中尺度模式(WRF-ARW)和(3DVAR)同化系统,对2007年7月16-17日发生在新疆东天山的暴雨天气过程进行了数值模拟试验。结果表明:①通过WRF模式对天山地区“2004-07-18”和“2010-06-22”两次暴雨过程的数值模拟试验, WRF模式在天山地区以Lin微物理过程方案,Grell积云对流参数化方案,MYJ边界层参数化方案,耦合Noah陆面模式的方案组合更优。②同化天山地区加密自动气象站资料、常规探空资料和AMSU-A/B辐射亮温资料,使初始场更接近于当时的实际大气状况,且同化后的模拟预报效果,无论是对降水落区还是对降水强度均较控制试验有明显的改善。③适当提高模式水平分辨率,有助于改进暴雨落区的预报强度。     

SIS海冰模式中两种盐度参数化方案的差异

详细介绍了Sea Ice Simulator(SIS)海冰模式中引进两种盐度参数化方案即等盐度方案和盐度廓线方案对海冰模拟所存在的差异。利用盐度廓线方案导出的表征盐度与海冰温度间关系的方程比等盐度方案多出一项,将该项定义为盐度差异项。盐度差异项对海冰厚度的热力作用表现为:在海冰厚度增长季节(11月到次年5月),盐度差异项通过升高海冰内部温度,抑制海冰增长;在消融的第一阶段(6—8月),盐度差异项通过升高海冰内部温度加快海冰消融;在消融的第二阶段(9—10月),盐度差异项通过降低海冰内部的温度抑制海冰消融。但尺度分析表明,盐度差异项要比方程中对海冰温度作用最大项小1—2个量级,如果采用一级近似,可以略去盐度差异项,因此盐度差异项对与海冰增长和消融影响很小。同时利用GFDL中心研制的冰-洋耦合模式Modular Ocean Model(MOM4),分别采用两种盐度参数化方案模拟北极海冰厚度和海冰密集度的季节性变化,模拟结果也表明两种方案模拟得到的海冰厚度和海冰密集度的季节性变化相差甚小。     

Simulated effect of soil freeze-thaw process on surface hydrologic and thermal fluxes in frozen ground region of the Northern Hemisphere

Soil freeze-thaw process is closely related to surface energy budget,hydrological activity,and terrestrial ecosystems.In this study,two numerical experiments(including and excluding soil freeze-thaw process)were designed to examine the effect of soil freeze-thaw process on surface hydrologic and thermal fluxes in frozen ground region in the Northern Hemisphere based on the state-of-the-art Community Earth System Model version 1.0.5.Results show that in response to soil freeze-thaw process,the area averaged soil temperature in the shallow layer(0.0175?0.0451 m)decreases by 0.35℃in the TP(Tibetan Plateau),0.69℃in CES(Central and Eastern Siberia),and 0.6℃in NA(North America)during summer,and increases by 1.93℃in the TP,2.28℃in CES and 1.61℃in NA during winter,respectively.Meanwhile,in response to soil freeze-thaw process,the area averaged soil liquid water content increases in summer and decrease in winter.For surface heat flux components,the ground heat flux is most significantly affected by the freeze-thaw process in both summer and winter,followed by sensible heat flux and latent heat flux in summer.In the TP area,the ground heat flux increases by 2.82 W/m2(28.5%)in summer and decreases by 3.63 W/m2(40%)in winter.Meanwhile,in CES,the ground heat flux increases by 1.89 W/m2(11.3%)in summer and decreases by 1.41 W/m2(18.6%)in winter.The heat fluxes in the Tibetan Plateau are more susceptible to the freeze-thaw process compared with the high-latitude frozen soil regions.Soil freeze-thaw process can induce significant warming in the Tibetan Plateau in winter.Also,this process induces significant cooling in high-latitude regions in summer.The frozen ground can prevent soil liquid water from infiltrating to deep soil layers at the beginning of thawing;however,as the frozen ground thaws continuously,the infiltration of the liquid water increases and the deep soil can store water like a sponge,accompanied by decreasing surface runoff.The influence of the soil freeze-thaw process on surface hydrologic and thermal fluxes varies seasonally and spatially.     

Evaluating effects of Dielectric Models on the surface soil moisture retrieval in the Qinghai-Tibet Plateau

Based on the measurement of L-band ground-based microwave radiometer(ELBARA-III type)in the Qinghai-Tibet Plateau and theτ-ωradiative transfer model,this research evaluated the effects of four soil dielectric models,i.e.,Wang-Schmugge,Mironov,Dobson,and Four-phase,on the L-band microwave brightness temperature simulation and soil moisture retrieval.The results show that with the same vegetation and roughness parameterization scheme,the four soil dielectric models display obvious differences in microwave brightness temperature simulation.When the soil moisture is less than 0.23 m3/m3,the simulated microwave brightness temperature in Wang-Schmugge model is significantly different from that of the other three models,with maximum differences of horizontal polarization and vertical polarization reaching 8.0 K and 4.4 K,respectively;when the soil moisture is greater than 0.23 m3/m3,the simulated microwave brightness temperature of Four-phase significantly exceeds that of the other three models;when the soil moisture is saturated,maximum differences in simulated microwave brightness temperature with horizontal polarization and vertical polarization are 6.1 K and 4.8 K respectively,and the four soil dielectric models are more variable in the microwave brightness temperature simulation with horizontal polarization than that with vertical polarization.As for the soil moisture retrieval based on the four dielectric models,the comparison study shows that,under the condition of horizontal polarization,Wang-Schmugge model can reduce the degree of retrieved soil moisture underestimating the observed soil moisture more effectively than other parameterization schemes,while under the condition of vertical polarization,the Mironov model can reduce the degree of retrieved soil moisture overestimating the observed soil moisture.Finally,based on the Wang-Schmugge model and FengYun-3C observation data,the spatial distribution of soil moisture in the study area is retrieved.     

2017年夏季秦岭降水的数值模拟及其空间分布

秦岭是中国气候的分界线,其高海拔地区降水数据的缺乏,限制了对秦岭高山气候及水资源变化的研究。为认知秦岭地区降水的空间分布,利用WRF模式,采用GF、KF和BMJ 3种对流参数化方案,对秦岭及周边地区2017年夏季降水进行了模拟。结果显示,模拟和卫星数据都揭示降水与地形呈剖面一致性,从南到北呈条带状分布的空间格局,在秦岭出现一条降水高值带,秦岭对水汽向北输送有明显的阻滞作用,造成南坡降水量显著大于北坡。模拟值要比卫星降水数据偏大,其中KF方案模拟降水偏多的主要原因为强烈的对流不稳定性导致对流降水的过高估计,GF方案则由于大气偏湿激发网格尺度降水,造成大尺度降水的模拟偏多,而BMJ方案模拟值与观测值最为接近。把模式分辨率提高到2 km,可显著改善模式对秦岭山区的模拟水平,但2 km高分辨率仍不足以完全显式解析积云对流过程,需要恰当的积云参数化的协同作用。     

RegCM4.6两种积云参数化方案在东亚模拟结果的评估

新一代区域气候模式RegCM4.6引进了Mix积云参数化方案,可以将之前版本中的Emanuel和Grell方案结合在一起,以弥补单个参数化方案的不足。利用2016年MODIS(Moderate-resolution imaging spectroradiometer)数据对RegCM4.6中Mix和Emanuel积云参数化方案模拟的东亚云量(Cloud fraction,CF)、冰水柱含量(Ice water path,IWP)和液水柱含量(Liquid water path,LWP)进行初步评估,计算了相关系数(r)、平均绝对误差(Mean absolute error,MAE)、平均偏差(Mean bias error,MBE)和均方根误差(Root mean square error,RMSE),以便为相关研究选取积云参数化方案提供参考依据。结果表明：(1)模拟的CF的MBE大致以胡焕庸线为界,西北部为轻微高估,东南部通常为低估。2种方案在夏季的模拟效果最好,冬季最差。Mix方案的MAE、MBE和RMSE的绝对值在四季普遍小于Emanuel方案。(2)模式明显低估了东亚的IWP,...     

WRF模式中3种闪电参数化方案的对比分析

利用WRF 模式V3.6.1 中3 种基于经验公式的闪电参数化方案对2013-07-05 江苏省内的1 次飑线过程进行数值模拟,对比3 种方案在模拟闪电位置和频数的特点,分析模拟结果的成因。结果表明：PR92w 方案能够很好地模拟闪电的分布位置;PR92z 方案的模拟结果和PR92w 方案正好相反,闪电频数与实测结果有很好的一致性但分布位置出现偏移;LPI 方案虽然能指示闪电分布区域,但对高闪电密度区的预报有一定误差。3 种方案的模拟结果表明：垂直风速和雷达回波是影响模式预报闪电频数和分布的重要因素。因此,增大网格模拟精度和同化雷达资料是改善PR92 方案模拟效果的重要途径。     

陆面水文—气候耦合模拟研究进展

陆面水文过程是全球/区域气候模式十分重要而又十分薄弱的环节。本文通过科学文献计量法,分析了陆面水文—气候耦合研究的发展状况及研究热点和趋势,并进一步对研究中存在的问题和挑战进行综述和探讨。现有多数气候模式中的陆面模式主要基于一维垂向结构设计,缺乏对流域尺度水文过程的精细描述,尤其缺乏下垫面人类活动影响的描述。因此,为了科学认识水文过程与气候变化的相互反馈作用机制,大量研究主要通过耦合流域水文模型与气候模式,研究不同时空尺度下水文过程变化的气候反馈效应。陆面过程模式中水文过程的改进和大尺度水文模型发展为陆面水文—气候耦合模拟奠定基础,在此基础上,陆面水文—气候耦合研究正从传统的单向耦合研究逐步发展为考虑气候—水文反馈的双向耦合研究。然而,双向耦合研究远未成熟,其问题集中表现为陆面水文—气候模型耦合过程如何匹配并提高系统稳定性、研制有效的尺度转换方案、完善参数化方案并评估参数不确定性、研制有效参数移植方法并提高模型适用性以及高分辨率甚至超分辨率模拟等方面,逐步解决上述问题并提高模拟精度是未来水文—气候耦合模拟研究的重要发展方向。     

藏东地区斜坡土壤冻融侵蚀力学机制及敏感性分析

冻融侵蚀是青藏高原草甸覆盖区的主要侵蚀方式,以气候条件一致的藏东地区斜坡表层土壤侵蚀为对象,基于区域地质条件和土体赋存特征,分析了土壤剥蚀输移的力学过程,探索了缓变的隐性因子和灾变的显性因子对冻融侵蚀的作用机制。结果表明：① 地表冻融侵蚀是自基岩风化和土壤演化起始,经历冻融拉裂破坏与沙土输移,到重新裸露基岩的一个循环演化过程;② 冻融侵蚀输移驱动力为受全因素控制的重力和径流动水压力,抗侵蚀力则由土壤二元结构特征和植被根系锚固力控制,多因素耦合作用下,冻融侵蚀演化具有“多因素,同表观”的特征;③ 基于土体极限平衡状态的破坏机制,冻融侵蚀分为牵引式的逐级崩解型和牵引—推移复合式的多级崩解型,两者在侵蚀强度和触发因素上均存在显著区别;④ 引入冻融侵蚀强度特征参数并建立其与构造、地形、植被、岩性等因子之间的关系式,提出基于小流域土壤侵蚀观测的冻融侵蚀强度估算公式。基于力学过程和土体破裂机制建立的侵蚀强度多因素估算方法,可有效减小区域适宜性产生的误差,并且为制定精细化的不同工况下的冻融侵蚀防治对策提供理论基础。     

基于GIS的城市用地布局多情景模拟与方案评价——以江苏省太仓市为例

用地布局一直是城市总体规划中的关键难题,以往的规划方案多凭借规划师的经验判断、简单的图层叠加和有限的公众参与生成。然而新的城乡规划法要求从多角度综合考虑城乡用地布局问题,更加注重规划过程的科学性和准确性,也就特别需要有一种适应多情景分析下的城市用地布局模拟与方案评价方法。文章以江苏省太仓市为例,借助GIS技术的强大空间分析功能,探讨了这种方法的技术流程：① 参考已有的大尺度城市未来模型,结合我国特别是研究区的区域特点,构建一个实用的城市未来模型(Urban Future Model,UFM);② 通过用地评价、战略归纳、情景模拟等步骤,生成若干个可选的用地布局模拟方案;③ 应用多目标达成矩阵法从粮食、生态、灾害等多个角度对这些方案进行综合评价,并根据评价结果选择一个最佳方案作为本轮总规用地布局的规划参考方案。