基于陆面模式Noah-MP的不同参数化方案在半干旱区的适用性

针对陆面模式Noah-MP对兰州大学半干旱气候与环境观测站（SACOL）2009年8月地表热通量模拟值偏差大的问题,通过分析相关物理过程和模拟试验来探究偏差的来源,并确定合适的参数化方案:采用Chen97方案计算感热输送系数可以改善感热通量的模拟;采用Jarvis气孔阻抗方案能增大植被蒸腾,改进模式对潜热通量的模拟效果,同时也使热通量在感热和潜热间的分配比例合理;采用LP92方案可减小土壤蒸发阻抗并有利于土壤蒸发,使得模式对潜热通量的模拟效果变好。不同参数化方案的组合试验表明:同时采用2组或3组新的参数化方案组合可以进一步减小模拟的地表感热和潜热通量的均方根误差,但是土壤湿度和温度的模拟效果并没有同步改善。     

青藏高原夏季土壤有机质及砾石影响水热 传输特性的数值模拟

本文针对青藏高原部分地区土壤有机质和砾石含量较高的特点,在前人工作的基础上,发展了一个新的参数化方案以描述土壤有机质和砾石对土壤导热率、导水率的影响。通过对通用陆面模式CoLM中的土壤水、热参数化方案以及地表蒸发阻抗三方面的逐步改进,对青藏高原藏东南站和纳木错站两种不同下垫面进行单点数值模拟分析。对比原方案与最终优化方案的模拟结果表明:采用新方案的CoLM模式对藏东南站土壤湿度模拟性能明显提高,平均偏差减小到0.04,而对纳木错站浅层20 cm以上土壤湿度的模拟偏差略微增大。新方案在藏东南站对土壤内部温度的模拟改善较为显著,平均偏差减小了0.2℃;而在纳木错站40 cm以上有所改进。新参数化方案较好地模拟了两个观测站表面能量通量的时间变化,纳木错站7、8月份的潜热通量改进尤为明显,比原方案减少大约20 W m-2,与观测结果较为接近。     

根系吸水过程参数化方案对青藏高原陆面过程模拟的影响研究

根系吸水过程对地表能量平衡和水循环起着重要作用，目前不同的根系吸水过程参数化方案对青藏高原陆面过程模拟的影响尚不明确，探讨相关参数化方案的影响，可以为今后建立陆面过程模式根系参数化方案提供参考。本文利用2010年6月1日至9月30日青藏高原玛曲站的观测资料作为大气强迫资料，驱动BCC_AVIM模式（北京气候中心陆面模式）引入不同的根系吸水过程参数化方案，对玛曲站2010年6月1日至9月30日时段感热通量、潜热通量、土壤温度、土壤含水量等要素进行数值模拟，分析根系吸水过程参数化方案对青藏高原地区陆面过程的影响。模式中有关根系吸水过程的参数化方案主要分为根分布模型和土壤水分对根系有效性函数两类，根分布模型用Jackson方案、Schenk方案替换，土壤水分对根系有效性函数用Li方案、LSM1.0方案、CLM4.5方案替换。对比结果表明:不同的根系吸水过程参数化方案对土壤温度、土壤含水量的模拟影响较小，对感热通量、潜热通量模拟影响较大，尤其对冠层蒸腾量模拟差异显著，相关参数化方案的变动直接影响冠层蒸腾量。两类方案模拟的差异受降水的影响，在多雨期，根分布对比方案与原模式方案模拟的感热、潜热通量间存在较大差异；在少雨期，土壤水分对根系有效性函数对比方案与原模式方案模拟的感热、潜热通量间存在较大差异。     

Improving CLM4.5 simulations of land–atmosphere exchange during freeze–thaw processes on the Tibetan Plateau

Soil is heterogeneous and has different thermal and hydraulic properties, causing varied behavior in heat and moisture transport. Therefore, soil has an important effect on land–atmosphere interactions. In this study, an improved soil parameterization scheme that considers gravel and organic matter in the soil was introduced into CLM4.5 (Community Land Model). By using data from the Zoige and Madoi sites on the Tibetan Plateau, the ability of the model to simultaneously simulate the duration of freeze–thaw periods, soil temperature, soil moisture, and surface energy during freeze–thaw processes, was validated. The results indicated that: (1) the new parameterization performed better in simulating the duration of the frozen, thawing, unfrozen, and freezing periods; (2) with the new scheme, the soil thermal conductivity values were decreased; (3) the new parameterization improved soil temperature simulation and effectively decreased cold biases; (4) the new parameterization scheme effectively decreased the dry biases of soil liquid water content during the freezing, completely frozen, and thawing periods, but increased the wet biases during the completely thawed period; and (5) the net radiation, latent heat flux, and soil surface heat flux of the Zoige and Madoi sites were much improved by the new organic matter and thermal conductivity parameterization.     

Simulation of the Bare Soil Surface Energy Balance at the Tongyu Reference Site in Semiarid Area of North China

The performance of a 1-D soil model in a semiarid area of North China was investigated using observational data from a cropland station at the Tongyu reference site of the Coordinated Enhanced Observing Period (CEOP) during the non-growing period, when the ground surface was covered with bare soil. Comparisons between simulated and observed soil surface energy balance components as well as soil temperatures and water contents were conducted to validate the soil model. Results show that the soil model could produce good simulations of soil surface temperature, net radiation flux and sensible heat flux against observed values with the RMSE of 1.54oC, 7.71 W m-2 and 27.79 W m-2, respectively. The simulated volumetric soil water content is close to the observed values at various depths with the maximal difference between them being 0.03. Simulated latent heat and ground heat fluxes have relatively lower errors in relative to net radiation and sensible heat flux. In conclusion, the soil model has good capacity to simulate the bare soil surface energy balance at the Tongyu cropland station and needs to be further tested in longer period and at more sites in semiarid areas.     

基于CLM4.5的高寒草地辐射收支和水热交换的数值模拟研究

利用那曲高寒气候环境观测研究站本部BJ观测点2014年6-8月的近地层观测资料,结合CLM4.5陆面模型,探究空气动力学粗糙长度、叶面积指数、植被覆盖度和热力学粗糙长度参数化方案的改变对陆面能水平衡的模拟结果产生的影响,并且探讨了粗糙度及植被状态指数影响陆面能水平衡模拟性能的机制.结果表明:(1)CLM4.5默认的热力...     

Evaluation of the profile and the resistance method for estimation of surface fluxes of momentum,sensible and latent heat

Summary Hourly lysimetric and micrometeorological data taken over a grass surface at the Meteorological Research Unit, Cardington U.K. have been analysed. A temperature difference and measurements of wind speed at only one height, combined with an independently estimated effective roughness length allowed sensible heat and momentum fluxes determination by the profile method on an hourly basis. The estimates are compared with direct measurements of sensible heat and friction velocity obtained by the eddy correlation method. The sensible and latent heat fluxes are also modelled by the resistance method. Equations based on the Monin—Obukhov similarity theory are used to account for stability effects through various forms of parameterization Aerodynamic and surface resistances, necessary for the Penman—Monteith equation are calculated from routinely measured meteorological data. The profile method for estimation of sensible heat flux and friction velocity is found to work excellently on the discussed daytime experimental data which correspond mainly to near neutral or slightly unstable conditions.Surface latent and sensible heat fluxes can also be described very well by the resistance method. A slightly better estimate of the sensible heat flux is achieved when stability corrections are taken into account. On the contrary Penman-Monteith equation for estimating latent heat flux is insensitive to adjustments for atmospheric stability.The comparison of the various methods leads to the establishment of empirical relationships which correlate various quantities such as soil heat flux, resistances, evapotranspiration etc. to routinely measured meteorological data.With 8 Figures     

Bulk Formulation of the Surface Heat Flux

An interpretive literature survey examines different approachesfor applying the bulk aerodynamic formulato predict the surface heat flux. The surface heat flux is often predicted in terms of the surface radiation temperature, which is also used to predict the upward longwave radiation and the heat flux into the soil. In models, the thermal roughness length based on the surface radiation temperature (radiometric roughness length) is often specified to be smaller than the roughness length for momentum for a number of distinct reasons. The definition of the radiometric roughness length depends on the way that the surface temperature is measured, the choice of stability functions and displacement height and inclusion of any additional resistances.Using airborne eddy correlation data collected over eight different sites including bare soil, crops and grassland and several types of forests, the radiometric roughness length is found to vary by orders of magnitude in a manner that is difficult to formulate. Alternatively, we evaluate the approach where the thermal roughness length is equated with the better behaved roughness length for momentum and the corresponding aerodynamic surface temperature is modelled in terms of the surface radiation temperature, solar radiation, and vegetation index. The influence of wind speed and soil moisture on the difference between the aerodynamic and surface radiation temperatures is also examined.     

Sensitivity study and validation of a land surface parameterization using the HAPEX-MOBILHY data set

A simple parameterization of land surface processes, amenable to the structure of a two-layer soil model, including a representation of the vegetation, has been designed for use in meteorological models. Prior to implementation in a mesoscale model, it is necessary to check the components and to verify the good working order of the parameterization as a whole. The aims of this paper then are: (i) evaluation and a sensitivity study of the various components of the model, specifying the needed accuracy for the parameters; (ii) micrometeorological validation of the model against the HAPEX-MOBILHY data set.First, we present the basic scheme. The focus is on the parameterization of surface resistance, and especially on its relationship with soil moisture.A sensitivity study is then performed through a set of one-dimensional simulations which allow a full interaction between the ground and the atmosphere. Above bare ground, it is shown that both soil texture and initial moisture greatly influence the outcome of the simulation. Latent heat flux ranges from that associated with potential evaporation through a switch-like behavior to that of dry soil. Next, the effects of transpiring vegetation canopies on the physical processes involved and the surface energy balance are examined. The sensitivity of the latent heat flux to changes in the soil and canopy parameters is emphazised; the major influence of the initial mean soil moisture and of the vegetation cover is pointed out. Finally, the evolution of the boundary layer in response to various surface conditions is studied.A validation of the land surface scheme is conducted through daily cycles during cloudless days. Simulated turbulent fluxes are successfully compared to micrometeorological measurements over a maize field at different growth stages. Over a pine forest, the correct simulation of the turbulent fluxes is obtained with an adequate parameterization of the surface resistance accounting for the atmospheric moisture deficit.     

Evaluation of the CAM and PX Surface Layer Parameterization Schemes for Momentum and Sensible Heat Fluxes Using Observations

In this study, the performances of the Community Atmosphere Model (CAM) and Pleim–Xiu (PX) surface layer parameterization schemes are investigated by using field observations. The parameterization schemes are evaluated against continuous momentum and sensible heat flux observations measured at two flat and homogeneous grassland sites in the suburb of Nanjing, eastern China. The observations were conducted from 30 December 2014 to 18 April 2017 at Jiangxinzhou and from 9 February 2015 to 26 March 2018 at Jiangning. It is found that the momentum flux is overall in good agreement with the observation, and the sensible heat flux is overestimated. The parameterizations of the momentum and sensible heat fluxes well capture the diurnal and seasonal patterns seen in the observations at the two sites. At Jiangxinzhou, the PX parameterization underestimates the momentum flux throughout the day and the CAM parameterization slightly overestimates it around the noon, while they underestimate the momentum flux throughout the year. The two parameterizations overestimate the sensible heat flux in the daytime as well as over the entire year. At Jiangning, the two parameterizations overestimate the momentum flux throughout the day and the sensible heat flux in the daytime, and overestimate both of them over the entire year. The two parameterizations are not significantly different from each other in reproducing the turbulent fluxes at the same site, while they perform differently at the two sites in terms of statistics. In addition, the parameterized fluxes increase with increased roughness length.     

Bare-ground surface heat and water exchanges under dry conditions: Observations and parameterization

A simplified land-surface parameterization is tested against bare-soil data collected during the EFEDA experiment conducted in Spain in June 1991. A complete data set, made up of soil properties as well as hydrological and atmospheric measurements, is described and discussed. The 11-day data set is characterized by very dry conditions and high surface temperatures during the day. Large values of sensible and soil heat fluxes and small values of surface evaporation (≈1 mm/day) were observed. This data set was modelled, leading to the following conclusions:

1. In the model, the parameterization provides values of the soil thermal properties and subsequently of the predicted soil heat fluxes which are overestimated when compared with the observations.
2. Following the literature, a value of the ratio between the roughness lengths for momentumZ _oand heatZ _ohof close to 10 for fairly homogeneous areas of bare soil and vegetation is used. This value leads to a fair prediction of the surface temperature. If the roughness lengths were taken to be equal, as is often assumed in atmospheric modelling, a poorer prediction results.
3. Finally, the vapor phase transfer mode is found dominant close to the surface and a modified parameterization including this effect is proposed. It allows a fair prediction of both surface evaporation and near-surface water content.

     

Further Development of the Vegetated Urban Canopy Model Including a Grass-Covered Surface Parametrization and Photosynthesis Effects

The vegetated urban canopy model (VUCM), which includes parametrizations of urban physical processes for artificial surfaces and vegetated areas in an integrated system, has been further developed by including physical processes associated with grass-covered surfaces in urban pervious surfaces and the photosynthesis effects of urban vegetation. Using measurements made from three urban/suburban sites during the BUBBLE field campaign in 2002, the model’s performance in modelling surface fluxes (momentum flux, net radiation, sensible and latent heat fluxes and storage heat flux) and canopy air conditions (canopy air temperature and specific humidity) was critically evaluated for the non-precipitation and the precipitation days. The observed surface fluxes at the urban/suburban sites were significantly altered by precipitation as well as urban vegetation. Especially, the storage heat at urban surfaces and underlying substrates varied drastically depending on weather conditions while having an important role in the formation of a nocturnal urban surface layer. Unlike the nighttime canopy air temperature that was largely affected by the storage-heat release, the daytime canopy air conditions were highly influenced by the vertical turbulent exchange with the overlying atmosphere. The VUCM well reproduced these observed features in surface fluxes and canopy air conditions at all sites while performing well for both the non-precipitation and the precipitation days. The newly implemented parametrizations clearly improved the model’s performance in the simulation of sensible and latent heat fluxes at the sites, more noticeably at the suburban site where the vegetated area fraction is the largest among the sites. Sensitivity analyses for model input parameters in VUCM elucidated the relative importance of the morphological, aerodynamic, hydrological and radiative/thermal properties in modelling urban surface fluxes and canopy air conditions for daytime and nighttime periods. These results suggest that the VUCM has great potential for urban atmospheric numerical modelling for a range of cities and weather conditions in addition to having a better physical basis in the representation of urban vegetated areas and associated physical processes.     

Impact of Atmospheric Surface-layer Parameterizations in the new Land-surface Scheme of the NCEP Mesoscale Eta Model

We tested three atmospheric surface-layer parameterization schemes (Mellor-Yamadalevel 2, Paulson, and modified Louis), both ina 1-D mode in the new NCEP land-surface scheme against long-term FIFE and HAPEX observations, and in a coupled 3-D mode withthe NCEP mesoscale Eta model. The differences inthese three schemes and the resulting surface exchange coefficients do not, in general, lead to significant differences in model simulated surface fluxes, skin temperature, andprecipitation, provided the same treatment of roughness length for heat is employed.Rather, the model is more sensitive to the choice of the roughness length for heat. To assess the latter, we also tested two approaches to specifythe roughness length for heat: 1) assuming the roughness length for heat is a fixed ratio of the roughness length for momentum, and 2) relating this ratio to the roughness Reynolds number as proposed by Zilitinkevich.Our 1-D column model sensitivity tests suggested that the Zilitinkevich approach can improve the surface heat fluxand skin temperature simulations. A long-term test with the NCEP mesoscaleEta model indicated that this approach can also reduce forecast precipitation bias. Based on these simulations, in January 1996 we operationally implemented the Paulsonscheme with the new land-surface scheme of the NCEP Eta model, along with the Zilitinkevich formulation to specify the roughness length for heat.     

Evaluation of the Urban Tile in MOSES using Surface Energy Balance Observations

The UK Met Office has introduced a new scheme for its urban tile in MOSES 2.2 (Met Office Surface Exchange Scheme version 2.2), which is currently implemented within the operational Met Office weather forecasting model. Here, the performance of the urban tile is evaluated in two urban areas: the historic core of downtown Mexico City and a light industrial site in Vancouver, Canada. The sites differ in terms of building structures and mean building heights. In both cases vegetation cover is less than 5%. The evaluation is based on surface energy balance flux measurements conducted at approximately the blending height, which is the location where the surface scheme passes flux data into the atmospheric model. At both sites, MOSES 2.2 correctly simulates the net radiation, but there are discrepancies in the partitioning of turbulent and storage heat fluxes between predicted and observed values. Of the turbulent fluxes, latent heat fluxes were underpredicted by about one order of magnitude. Multiple model runs revealed MOSES 2.2 to be sensitive to changes in the canopy heat storage and in the ratio between the aerodynamic roughness length and that for heat transfer (temperature). Model performance was optimum with heat capacity values smaller than those generally considered for these sites. The results suggest that the current scheme is probably too simple, and that improvements may be obtained by increasing the complexity of the model.     

Parameterization of the thermal impacts of sub-grid orography on numerical modeling of the surface energy budget over East Asia

Summary A parameterization scheme for the thermal effects of subgrid scale orography is incorporated into a regional climate model (developed at Nanjing University) and its impact on modeling of the surface energy budget over East Asia is evaluated. This scheme includes the effect of terrain slope and orientation on the computation of solar and infrared radiation fluxes at the surface, as well as the surface sensible and latent heat fluxes. Calculations show that subgrid terrain parameters alter the diurnal cycle and horizontal distributions of surface energy budget components. This effect becomes more significant with increased terrain slope, especially in winter. Due to the inclusion of the subgrid topography, the surface area of a model grid box changes over complex terrain areas. Numerical experiments, with and without the subgrid scale topography scheme, show that the parameterization scheme of subgrid scale topography modifies the distribution of the surface energy budget and surface temperature around the Tibetan Plateau. Comparisons with observations indicate that the subgrid topography scheme, implemented in the climate model, reproduces the observed detailed spatial temperature structures at the eastern edge of the Tibetan Plateau and reduces the tendency to overestimate precipitation along the southern coastal areas of China in summer.     

On using mixed-layer transport parameterizations with radiometric surface temperature for computing regional scale sensible heat flux

Recent mixed-layer formulations for computing large-scale surface energy fluxes under daytime convective conditions do not require the estimation of surface-layer parameters, such as the roughness lengths for momentum and heat. This greatly simplifies approaches using operational satellite measurements of surface temperature for computing the surface energy balance at regional scales because the surface roughness parameters are not well known for many landscapes. The utility of such mixed-layer formulations is tested using data from several recent multidisciplinary field experiments (HAPEX-MOBILHY, FIFE and Monsoon 90). The results indicate that specific mixed-layer formulations adequately simulate surface sensible heat fluxes in the grassland and shrubland sites. However, use of the original values of proposed empirical coefficients for the forested site yield poor results. This is probably due to the fact that the forested site has significantly different surface geometry and associated distribution of temperature among the surface components (especially the relative importance of soil background temperatures) compared to the other sites. Therefore, the relationship between aerodynamic and radiometric surface temperature may have greatly differed between the forested site and the other locations. However, differences in aerodynamic roughness between the experimental sites were not correlated with changes required in the values of the coefficients. Instead, a two-source model which makes the distinction between aerodynamic and radiative temperature is proposed, as a means to determine which surface properties significantly affect the magnitude of the mixed-layer coefficients.     

MM5和ETA相似理论近地层方案对农田下垫面通量模拟比较研究

近地层湍流通量计算对于中尺度数值模式有重要意义, 湍流通量的参数化是当前大气边界层研究的重要课题之一。选择青藏高原东缘大理观象台边界层通量观测系统, 离线测试了WRF区域模式中的两种常用的近地层参数化方案(MM5相似理论非迭代方案A和ETA 相似理论迭代方案B), 并将参数化方案计算结果与边界层铁塔涡动相关法的观测值进行对比分析。在大理观象台观测场不同植被随季节交替的状况下, 根据边界层铁塔4层高度风速拟合, 发现近地层空气动力学粗糙度随季节变化特征明显。将拟合的空气动力学粗糙度输入模式参数化方案进行通量计算。结果表明:稳定度是影响近地层参数化方案精度的重要因素, 在不稳定条件下方案B低估了动量通量, 方案A优于方案B, 而在稳定条件下方案A低估了动量通量, 方案B优于方案A, 两种方案总体来看误差不大。对于大理边界层通量观测场地农田植被交替的环境条件, 不同季节下垫面植被类型的差异, 以及植被的稀疏对近地层参数化方案湍流通量计算结果的精度有显著影响。方案B考虑了空气动力学粗糙度z₀和热量粗糙度z_0h的差异, 不稳定条件下感热通量计算结果在裸土或稀少植被条件下明显优于方案A。针对方案B不稳定条件下感热通量计算结果在裸土下垫面仍出现高估的现象, 在使用了(Zeng, et al, 1998)提出的对于使用辐射地表温度在裸土下垫面时的订正方法后, 计算结果也有明显改善。     

黄土高原塬区地表辐射和热量平衡观测与分析

利用2005年夏季黄土高原塬区陆面过程野外试验(LOPEX05)的观测资料,初步分析了甘肃平凉黄土高原塬区地表辐射收支和热量平衡特征。结果显示,黄土高原塬区地面长波辐射大于大气长波辐射,典型晴天、阴天和雨天情况下两者平均差值分别为65,25和8 W.m-2;对于地气能量交换各个分量而言,黄土高原塬上和塬下在相同下垫面下的差别不大,但裸地和有植被的下垫面差别很明显;在白天,潜热在净辐射中所占的比重较大,其次是感热,最后是土壤热通量。对能量平衡中的储存项如热通量板上层土壤的热存储和植被冠层存储进行了估算,结果表明,土壤的热储存项在-30~70W.m-2之间,而植被的热能储存项在-10~25 W.m-2之间。在考虑估算的存储项之后,能量平衡散布图斜率由0.68提高到0.79,相关系数R由0.90提高到0.93,两者分别提高了11.0%和3.0%,并对能量不平衡有明显的改进,说明能量储存项在地表能量闭合中必须考虑。     

MM5和ETA相似理论近地层方案对农田下垫面通量模拟比较研究

近地层湍流通量计算对于中尺度数值模式有重要意义,湍流通量的参数化是当前大气边界层研究的重要课题之一.选择青藏高原东缘大理观象台边界层通量观测系统,离线测试了WRF区域模式中的两种常用的近地层参数化方案(MM5相似理论非迭代方案A和ETA相似理论迭代方案B),并将参数化方案计算结果与边界层铁塔涡动相关法的观测值进行对比分析.在大理观象台观测场不同植被随季节交替的状况下,根据边界层铁塔4层高度风速拟合,发现近地层空气动力学粗糙度随季节变化特征明显.将拟合的空气动力学粗糙度输入模式参数化方案进行通量计算.结果表明:稳定度是影响近地层参数化方案精度的重要因素,在不稳定条件下方案B低估了动量通量,方案A优于方案B,而在稳定条件下方案A低估了动量通量,方案B优于方案A,两种方案总体来看误差不大.对于大理边界层通量观测场农田植被交替的环境条件,不同季节下垫面植被类型的差异,以及植被的稀疏对近地层参数化方案湍流通量计算结果的精度有显着影响.方案B考虑了空气动力学粗糙度z₀和热量粗糙度z_0h的差异,不稳定条件下感热通量计算结果在裸土或稀少植被条件下明显优于方案A.针对方案B不稳定条件下感热通量计算结果在裸土下垫面仍出现高估的现象,使用了Zeng等1998年提出的用辐射地表温度订正裸土下垫面感热能量方法后,计算结果也有明显改善.