陆面气象要素非均匀分布对模式计算的影响及其参数化初步探讨

为改进陆面过程的参数化，初步探讨了模式中陆面气象要素次网格尺度不均匀性分布对模拟计算结果的影响。应用平均化方法初步探讨陆面过程的参数化，并分析了气象要素次网格尺度不均匀性分布对计算结果的影响特征。结果表明，次网格尺度不均匀性分布对模式的计算结果会产生一定的偏差，可用平均化方法对陆面过程次网格尺度不均匀性分布进行初步参数化。     

尺度自适应大气边界层参数化改进及其对一次海雾的数值模拟研究

大气边界层湍流运动是地球大气运动最重要的能量输送过程之一。当数值模式分辨率接近活跃含能湍涡长度尺度时，湍流运动被部分解析，被称为“灰色区域”，传统的边界层方案不适合此时模式湍流问题的描述。为了提高模式边界层方案在包括“灰色区域”的不同网格尺度上的描述能力，适应不同分辨率模式的需要，在雷诺平均湍流理论基础上，修正Mellor-Yamada-Nakanishi-Niino (MYNN)方案湍流长度尺度参数和非局地湍流的参数表达，改进湿度和温度在“灰色区域”的湍流输送参数化及对网格尺度的自适应能力。利用改进的MYNN尺度自适应方案，分别采用3 km和1 km、1.5 km和0.5 km分辨率单向嵌套网格WRF中尺度模式，对2014年2月26日的一次黄海海雾过程进行模拟试验，检验不同分辨率下改进后的MYNN大气边界层参数化方案的合理性和对海雾的模拟效果。尺度自适应MYNN大气边界层参数化方案在千米级网格尺度上获得稳定、合理的湍流垂直输送计算结果。参照雾区卫星云图，不同分辨率模式低层云水混合比模拟结果具有稳定表现，模拟的雾区分布和温度、湿度等物理量结构都较好地再现了再分析“观测事实”，初步表明该参数化方案有较高的网格尺度自适应能力。      

边界层参数化方案在“灰色区域”尺度下的适用性评估

随着数值预报模式分辨率的提高,当模式网格距与含能湍涡的长度尺度相当时,模式动力过程可解析一部分湍流运动,而剩余的湍流运动仍需参数化,此时便产生了湍流参数化的“灰色区域”问题。对传统的PBL（Planetary Boundary Layer）方案在“灰色区域”下的适用性评估,是改进PBL方案以使其能够适应分辨率变化的前提和基础。本研究基于干对流边界层的大涡模拟试验,比较了WRF（Weather Research and Forecast Model）模式中四种常用的边界层参数化方案[YSU（Yonsei University）、MYJ（Mellor-Yamada-Janjic）、MYNN2.5（Mellor-Yamada-Nakanishi-Niino Level 2.5）、MYNN3）]在“灰色区域”尺度下的表现。研究表明,混合层内总热通量对所使用的参数化方案和水平分辨率均不敏感。不同参数化方案中次网格与网格通量的比例表现出对水平网格距不同的依赖性。局地PBL方案（MYJ、MYNN2.5）在混合层内的平均位温随网格距减小而增大,次网格通量随网格距减小而减小,较参考湍流场对次网格通量有所低估。YSU方案的非局地项几乎不随水平格距改变而变化,对次网格通量的表征并未表现出较强的分辨率依赖性,且过强的非局地次网格输送使混合层内温度层结呈弱稳定,抑制了可分辨湍流输送,不易于激发次级环流。MYNN3方案的非局地次网格通量（负梯度输送项）随网格距减小而减小,使其对次网格通量的表征具有较好的分辨率依赖性。PBL方案在“灰色区域”尺度下的适用性与具体分辨率有关。以分辨率500 m为例,四种PBL方案中不存在一种最佳方案,能对边界层的热力结构和湍流统计特征均有准确的描述。     

海冰在大气环流模式中的重要作用

文章简要综述了次网格尺度海冰非均匀性对大气环流模式性能的影响；南极冰在全球环流和短期气候变化中的作用；以及模式中不同的海冰反照率参数化对地表温度和辐射的影响等研究结果．说明海冰对极地海洋和大气的能量收支及短期气候变化有重要作用，不同的海冰参数化方案对气候模拟结果有重要影响．     

次网格地形坡度坡向参数化及其对区域气候模拟的影响

随着数值模式水平分辨率的提高，模式下垫面特征的描述更加仔细，针对东亚地区复杂的地形和植被特征，发展适合非均匀下垫面的地表通量参数化方案，对改进数值模式和刻画高原附近地区复杂地形动力和热力。效应的能力非常必要。本文通过计算次网格地形坡度、坡向参数及其对到达地面短波辐射通量的影响，提出次网格地形热力效应的参数化方案，改进数值模式中复杂地形区域地面热量平衡的计算。利用p—σ区域气候模式的数值试验结果表明，次网格地形热力效应参数化方案的引入，对东亚地区夏季气候尤其是降水的模拟有明显的改进。     

陆面特征非均一作用参数化及其对区域气候影响的数值模拟试验

本文首先改进了陆地下垫面特征非均一性的次网格尺度参数化方法，然后利用三维地气耦合的区域气候模式，设计一系列值试验，研究了下垫面特征改变对区域气候环境变化的影响，主要分析了陆地表面特征变化对我国苏南附近地区夏季温度化的影响，结果表明采用地下垫面特征非均一作用的次网格尺度参数化方法对于改进数值模拟结果的质量有一定效果。     

降水气候强迫下非均匀地表区域平均径流的一种参数化方案

将任一中尺度区域的平均瞬间径流率考虑为区域平均降水量和地表土壤层水分渗透垦的余项.根据降水量在地理空间上分布的实测资料拟合其空间概率密度函数(PDF),并结合土壤入渗物理过程的数学描述及其经验公式,精确估计出地表土壤渗透率及其空间分布,由此建立区域地表径流率的统计-动力学估计方案.换言之,区域内地表产流率可视为区域平均降水量与区域平均的土壤下渗量之差值,而区域内土壤的平均下渗量又町分为非饱和区和饱和区两部分的下渗量来分别计算.就陆面水分循环的物理过程而言,地表入渗现象是在一定的下垫面特性基础上,由一定的水分供应源而形成的.根据大气降水向地表层输送水分的物理过程,在满足植被表层覆盖需水(截流水)和地表层土壤人渗水基础上,多余的降水量才会形成地表径流.凶此,推求地表产流率的主要关键在于地表土壤层需水量.为此奉文根据土壤水分通量方程推导出水分入渗公式.又从描述土壤水分和降水的空间PDF出发,推导出非均匀土壤含水量及降水气候强迫所形成的次网格尺度区域平均径流率计算公式.利用长江三角洲地区1996年降水量和土壤特性等实测资料建立区域平均地表径流率的估计公式,并对其影响凶素进行敏感性试验.结果表明,该方法与用Mosaic方法计算的区域径流率(或产流率)结果十分接近.由此可见,该文提出的降水气候强迫下非均匀地表区域平均径流的这种统计-动力参数化方案,具有相当的可靠性与可行性.     

利用尺度自适应大气边界层湍流参数化方案对一次陆地浓雾的数值模拟

为提高大气数值模式的模拟能力,改进大气边界层水汽、热量湍流输送计算和大雾天气的模拟效果,选用WRF三维非静力模式,采用具有局地和非局地垂直湍流尺度自适应计算能力的MYNN_SA参数化方案,对2017年12月28—29日我国华北—江淮地区的大范围浓雾过程进行了数值模拟研究,探讨从中尺度到灰区尺度分辨率范围,模式的尺度自适应大气边界层湍流参数化方案对稳定大气边界层发展、湍流输送和大雾天气模拟的影响。利用中国地面气象站观测资料和ERA5再分析数据,在接近灰区尺度的网格分辨率上,利用尺度自适应大气边界层湍流MYNN_SA参数化方案较之中尺度参数化MYNN方案,可明显改善次网格湍流输送计算,以及陆地浓雾的强度、空间分布和时间演变特征,可更精确地模拟云水混合比、逆温层和雾区的垂直结构。     

海冰在大气环流模型中的重要作用

文章简要综述了次网格尺度海非均匀性大对大气环流模式性能的影响，南极冰在全球环流和短期气候变化中的作用，以及模式中不同的海冰反照率参数化对地表温度和辐射的影响等研究结果，说明海冰对极地海洋和大气的能量收支及短期气候变化有重要作用，不同的海洋参数化方案对气候模拟结果有重要影响。     

基于降水气候强迫的一种地表径流估计方法

文中从地表水分平衡的物理机制出发 ,引进降水概率统计分布理论 ,推导出一种由降水气候强迫形成的次网格尺度非均匀径流率计算方法。应用于mosaic方案 ,可进一步计算区域平均径流及产流和汇流。试验证明 ,对于不同物理性质的地表而言 ,由于其土壤入渗能力的差异 ,相同降水气候强迫所能生成的径流量及其相对比值是不同的 ,例如干旱区与湿润区就有很大差异。同时 ,不同空间分布概率的降水量 ,不同站点的地表性质、土壤水力学条件等物理因素的千变万化 ,可使径流特征的空间分布具有很大变异性。因上述各种因素的综合影响 ,由降水气候强迫所形成的地表径流具有特定的非均匀分布是必然的。文中用实测资料验证了其可靠性与可行性     

Implementation of a Surface Runoff Model with Horton and Dunne Mechanisms into the Regional Climate Model RegCM_NCC

A surface runoff parameterization scheme that dynamically represents both Horton and Dunne runoff generation mechanisms within a model grid cell together with a consideration of the subgrid-scaie soil heterogeneity, is implemented into the National Climate Center regional climate model （RegCM_NCC）. The effects of the modified surface runoff scheme on RegCMANCC performance are tested with an abnormal heavy rainfall process which occurred in summer 1998. Simulated results show that the model with the original surface runoff scheme （noted as CTL） basically captures the spatial pattern of precipitation, circulation and land surface variables, but generally overestimates rainfall compared to observations. The model with the new surface runoff scheme （noted as NRM） reasonably reproduces the distribution pattern of various variables and effectively diminishes the excessive precipitation in the CTL. The processes involved in the improvement of NRM-simulated rainfall may be as follows： with the new surface runoff scheme, simulated surface runoff is larger, soil moisture and evaporation （latent heat flux） are decreased, the available water into the atmosphere is decreased; correspondingly, the atmosphere is drier and rainfall is decreased through various processes. Therefore, the implementation of the new runoff scheme into the RegCMANCC has a significant effect on results at not only the land surface, but also the overlying atmosphere.     

Sensitivity of a GCM simulation to subgrid infiltration and surface runoff

A subgrid parameterization of infiltration and surface runoff was evaluated using a land surface model coupled to an atmospheric general circulation model. Averaged over 5 year simulations, the subgrid parameterization resulted in significantly less infiltration of water into the soil compared to a simulation without subgrid hydrologic processes. As a result, the soils were drier, latent heat flux decreased, and surface air temperature increased. These results are consistent with other studies of subgrid hydrologic parameterizations, which also resulted in drier soils, decreased latent heat flux, and warmer surface temperatures. Several river basins were studied in detail. In the Amazon and Lena basins, the subgrid parameterization resulted in better annual runoff compared to observed annual river flow; surface air temperature was unchanged in the Amazon and better, compared to observations, in the Lena. In the Ob, Yenisey, and Amur basins, the subgrid parameterization resulted in too much annual runoff; July surface air temperature was unchanged or worse (Amur). Annual runoff for the Mississippi basin was better with the subgrid parameterization, but July surface air temperature was worse. These results suggest the utility of subgrid hydrologic parameterizations vary among river basins depending on the relative importance of Horton and Dunne runoff and the geologic factors affecting runoff generation.     

具有Horton及Dunne机制的径流模型在VIC模型中的应用(英)

地表径流主要由蓄满(Dunne)和超渗产流(Horton)机制产生;土壤性质的空间变异性、前期土壤水、地形及降水的空间变异性导致不同的径流机制。在研究区域或模型网格内,蓄满产流及超渗产流可能同时出现,缺乏考虑任何一种机制以及土壤性质的次网格空间变率可能导致地表径流的过高或过低估计,从而影响土壤水的计算。利用Philip入渗公式用于时间压缩逼近(TCA)给出了一种径流参数化方法,该方法可以动态实现模型网格中的Horton及Dunne产流机理,它考虑了土壤空间变异性对Horton和Dunne径流的影响。该径流模型应用到基于水文原理的陆面过程模型VIC,在淮河流域及美国宾西法尼亚州的一个流域进行了测试,结果表明:新的参数化方法对地表径流和土壤水分含量的分配起着重要作用,对于改进径流和土壤水的模拟有重要意义。     

基于WRF模式的清江流域降雨-径流模拟研究

基于WRF模式,采用4层嵌套方案,选取3种积云参数化方案和7种微物理方案组成21种组合,对清江流域2016—2018年6—10月6次典型降雨事件进行数值预报,结合CMORPH卫星-地面自动站-雷达三源融合降水产品,采用TS评分和FSS评分,分析不同分辨率和云微物理方案的降雨预报效果;基于较优组合方案的WRF模式与WRF-Hydro水文模式耦合进行径流模拟,分析WRF模式在水文模拟中的应用效果。结果表明:3 km和1 km分辨率对降雨中心位置及强度预报的差别不大,对降雨落区都有较好的预报能力;在积云参数化方案中,KF方案和BMJ方案的降雨预报效果优于GF方案;在微物理方案中,WSM3、WSM5、WSM6、Thompson方案的预报结果与融合数据有较好的一致性;基于较优组合方案BMJ_WSM3,将WRF模式与WRF-Hydro模式耦合,耦合模式能较好地模拟洪水过程,径流模拟相关系数都在0.67以上,且NSE最高可达0.79。      

Improving a subgrid runoff parameterization scheme for climate models by the use of high resolution data derived from satellite observations

In this study it is shown that the availability of a very high resolution dataset of land surface characteristics leads to the improvement of a surface runoff parameterization scheme. The improved parameterization scheme was developed for application in global and regional climate models and is a further development of the Arno scheme that is widely used in climate models. Here, surface runoff is computed as infiltration excess from a "bucket" type reservoir which takes the subgrid variability of soil saturation within a model gridbox into account. Instead of prescribing a distribution of subgrid scale soil water capacities as in the original Arno scheme, the array of high resolution soil water capacities taken from a global 1 km dataset of land surface parameters is used to obtain individual fractional saturation curves for each model gridbox. From each saturation curve, the three parameters (a shape parameter describing the shape of the subgrid distribution of soil water capacities, subgrid minimum and maximum soil water capacity) required in the modified formulation of the scheme are derived via optimization. As in the original Arno scheme applied in the ECHAM general circulation model and the REMO regional climate model, topography variations will influence the distribution of saturated subgrid areas within a model gridbox. At most gridboxes the net effect of these changes is such that more runoff is produced for high soil water contents and less runoff for low soil water contents. A validation of simulated discharge computed with a simplified land surface scheme applied to reanalysis data of the European Centre for Medium-Range Weather Forecasts and a hydrological discharge model has shown that these changes lead to a more realistic simulation of the annual cycle of discharge for several catchments. In particular this could be shown for the Yangtze Kiang and Amur catchments where adequate input data are available.     

Calculations of river-runoff in the GISS GGM: impact of a new land-surface parameterization and runoff routing model on the hydrology of the Amazon River

This study examines the impact of a new land-surface parameterization and a river routing scheme on the hydrology of the Amazon basin, as depicted by the NASA/Goddard Institute of Space Studies (GISS) global climate model (GCM). The more physically realistic land surface scheme introduces a vegetation canopy resistance and a six-layer soil system. The new routing scheme allows runoff to travel from a river's headwater to its mouth according to topography and other channel characteristics and improves the timing of the peak flow. River runoff is examined near the mouth of the Amazon and for all of its sub-basins. With the new land-surface parameterization, river run-off increases significantly and is consistent with that observed in most basins and at the mouth. The representation of the river hydrology in small basins is not as satisfactory as in larger basins. One positive impact of the new land-surface parameterization is that it produces more realistic evaporation over the Amazon basin, which was too high in the previous version of the GCM. The realistic depiction of evaporation also affects the thermal regime in the lower atmosphere in the Amazon. In fact, the lower evaporation in some portions of the basin reduces the cloudiness, increases the solar radiation reaching the ground, increases the net radiation at the surface, and warms the surface as compared to observations. Further GCM improvement is needed to obtain a better representation of rainfall processes.     

Testing a GCM land surface scheme against catchment-scale runoff data

A GCM land surface scheme was used, in off-line mode, to simulate the runoff, latent and sensible heat fluxes for two distinct Australian catchments using observed atmospheric forcing. The tropical Jardine River catchment is 2500 km² and has an annual rainfall of 1700 mm y^–1 while the Canning River catchment is 540 km², has a Mediterranean climate (annual rainfall of 800 mm y^–1) and is ephemeral for half the year. It was found that the standard version of a land surface scheme developed for a GCM, and initialised as for incorporation into a GCM, simulated similar latent and sensible heat fluxes compared to a basin-scale hydrological model (MODHYDROLOG) which was calibrated for each catchment. However, the standard version of the land surface scheme grossly overestimated the observed peak runoff in the wet Jardine River catchment at the expense of runoff later in the season. Increasing the soil water storage permitted the land surface scheme to simulate observed runoff quite well, but led to a different simulation of latent and sensible heat compared to MODHYDROLOG. It is concluded that this 2-layer land surface scheme was unable to simulate both catchments realistically. The land surface scheme was then extended to a three-layer model. In terms of runoff, the resulting control simulations with soil depths chosen as for the GCM were better than the best simulations obtained with the two-layer model. The three-layer model simulated similar latent and sensible heat for both catchments compared to MODHYDROLOG. Unfortunately, for the ephemeral Canning River catchment, the land surface scheme was unable to time the observed runoff peak correctly. A tentative conclusion would be that this GCM land surface scheme may be able to simulate the present day state of some larger and wetter catchments but not catchments with peaky hydrographs and zero flows for part of the year. This conclusion requires examination with a range of GCM land surface schemes against a range of catchments. Received: 9 June 1995 / Accepted: 4 April 1996     

陆面过程参数化对宁波地区雷暴过程模拟的影响

利用耦合复杂程度不同的陆面过程参数化方案(Noah、RUC)的新一代中尺度数值模式WRF,对2006年6月24日发生在宁波地区的一次典型的雷暴过程进行了数值模拟试验。结果表明:对于雷暴发生前期近地面热力、动力场的特征,Noah方案的模拟较为逼真,RUC方案没有反映出下垫面覆盖的多样性以及城市下垫面的影响,城乡之间差异不明显;Noah方案模拟的雷暴启动、发展过程与观测较为一致,RUC方案较好地描述出了演变过程中的关键阶段(3次合并过程);由于参数化所考虑的要素和物理过程存在一定差异,Noah方案在对降水的强度、降水中心位置的模拟方面具有一定优势;雷暴的持续时间对陆面过程参数化方案的选择比较敏感,两个方案所模拟的雷暴过程持续时间不同程度地长于实际雷暴持续时间;无论是哪种下垫面覆盖类型,白天Noah方案模拟的感热通量均大于RUC方案,而Noah方案模拟的潜热通量均小于RUC方案。