Parametrization of peatland hydraulic properties for the Canadian land surface scheme

Abstract

A hydraulic parametrization is developed for peatland environments in the Canadian Land Surface Scheme (CLASS). Three ‐wetland soil classes account for the typical variation in the hydraulic characteristics of the uppermost 0.5 m of organic soils. Review of the literature reveals that saturated hydraulic conductivity varies from a median of 1.0 × 10^?7m/s in deeply humified sapric peat to 2.8 × 10^?4 m/s in relatively undecomposed fibric peat. Average pore volume fraction ranges from 0.83 to 0.93. Parameters have been designed for the soil moisture characteristic curves for fibric, hemic and sapric peat using the Campbell (1974) equation employed in CLASS, and the van Genuchten (1980) formulation. There is little difference in modelled soil moisture between the two formulations within the range of conditions normally found in peatlands. Validation of modelled water table depth and peat temperature is performed for a fen in northern Québec and a bog in north‐central Minnesota. The new parametrization results in a more realistic simulation of these variables in peatlands than the previous version of CLASS, in which unrealistic mineral soil “equivalents “ were used for wetland soil climate modelling.     

陆面过程中冠层四流辐射传输模式的模拟性能检验

利用实际观测资料检验新发展的植被冠层四流辐射传输模式的模拟性能.将浙江大学田间观测数据用于模式试验,其结果表明:在可见光波段,因叶片散射率值较小,四流模式和二流模式模拟结果差别不大,模拟的冠层反照率都接近观测值.这表明在可见光波段,二流模式已经能够较好模拟芯层反照率,四流模式能够提高的精度范围有限;在近红外波段,因叶片散射率值较大,两个模式模拟结果差别较明显,四流模式模拟的冠层反照率相对二流模式的模拟结果更接近观测值.利用第2次雪模式比较计划SNOWMIP资料和改进的10层陆面模式BATS进行耦合试验,将四流模式以及二流模式均耦合到该陆面模式中,其目的是为考察四流模式对陆面模式模拟地表反射太刚短波辐射通量的影响.耦合后得到的结果与SNOWMIP中加拿大BERMS草地站和森林站的观测资料进行了对比.对比结果表明,采用网流模式、二流模式、BATS原辐射传输模式后,耦合四流模式的陆面模式模拟地表反射太阳短波辐射通量最接近观测值.说明采用四流模式能够改善陆面模式对地表反射太阳短波辐射通啦的模拟.     

Towards closing the vertical water balance in Canadian atmospheric models: Coupling of the land surface scheme class with the distributed hydrological model watflood

Abstract

Second generation land surface schemes are the subject of much development activity among atmospheric modellers. This work is aimed at, among other things, improving the representation of the soil water balance in order to simulate, more properly, exchanges with the atmosphere and to permit the use of model output to generate streamflow for model validation. The Canadian development program is centred on CLASS, the Canadian Land Surface Scheme, developed at Environment Canada. This paper focuses on the improvement of hydrology in CLASS. This was accomplished by designing a two‐way interface to WATFLOOD, a distributed hydrologic model developed at the University of Waterloo. The two models share many features, which facilitated the coupling procedure.

The interface retains the three‐layer vertical moisture budget representation in CLASS but adds three horizontal runoff possibilities. Runoff from the surface water follows Manning's equation for overland flow. Interflow is generated from the near‐surface soil layer using a parametrization of Richard's equation and base flow is produced by Darcian flow from the bottom of layer 3. An approximation of the internal topography of grid elements is used to supply horizontal gradients for the runoff components.

Tests are in progress in four Canadian study areas. Initial results are presented for the summer of 1993 for the Saugeen River in southwestern Ontario. The new scheme produces realistic hydrographs, whereas the old scheme did not. Bare ground evaporation is reduced by about 17% as a consequence of reduced water availability in layer 1. Evapotranspiration is not affected because the rooting depth extends into layer 3, in which soil moisture does not change appreciably with the new scheme. These results suggest that the new scheme improves the representation of streamflow in WATFLOOD/CLASS and of the soil moisture budget in CLASS. Work is in progress to validate this result over basins, such as the BOREAS study watersheds, where both runoff and evapotranspiration measurements are available.     

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     

A scheme for pixel-scale aerodynamic surface temperature over hilly land

Hilly-land satellite pixel-scale aerodynamic surface temperatures (AdST) are investigated using LAS (Large Aperture Scintillometer) and meteorological observations during 21-22 May 2001, indicating that the calculated temperatures are predominantly subject to estimated roughness lengths and, to a less extent, to estimated Bowen ratios, with errors to within 3.0 K between the AdST calculations and hilly radiometric surface temperatures retrieved from satellite data with the split window model. The errors depend heavily on the model used and the zenith angles and azimuth of the satellite and sun with respect to the observational site.     

Evaluation of an urban land surface scheme over a tropical suburban neighborhood

A two-step statistical downscaling method has been reviewed and adapted to simulate twenty-first-century climate projections for the Gulf of Fonseca (Central America, Pacific Coast) using Coupled Model Intercomparison Project (CMIP5) climate models. The downscaling methodology is adjusted after looking for good predictor fields for this area (where the geostrophic approximation fails and the real wind fields are the most applicable). The method’s performance for daily precipitation and maximum and minimum temperature is analysed and revealed suitable results for all variables. For instance, the method is able to simulate the characteristic cycle of the wet season for this area, which includes a mid-summer drought between two peaks. Future projections show a gradual temperature increase throughout the twenty-first century and a change in the features of the wet season (the first peak and mid-summer rainfall being reduced relative to the second peak, earlier onset of the wet season and a broader second peak).     

A study of the performance of a land surface process model (LSPM)

In this paper, a reliable Land-Surface Process Model (LSPM), which is a new version of the LPM of Ji and Hu (1989), is described. The LSPM has been validated with experimental data measured at two stations in the Po Valley (Northern Italy).     

The impact of implementing the bare essentials of surface transfer land surface scheme into the BMRC GCM

This study describes the first order impacts of incorporating a complex land-surface scheme, the bare essentials of surface transfer (BEST), into the Australian Bureau of Meteorology Research Centre (BMRC) global atmospheric general circulation model (GCM). Land seasonal climatologies averaged over the last six years of integrations after equilibrium from the GCM with BEST and without BEST (the control) are compared. The modeled results are evaluated with comprehensive sources of data, including the layer-cloud climatologies from the international satellite cloud climatology project (ISCCP) data from 1983 to 1991 and the surface-observed global data of Warren et al., a five-year climatology of surface albedo estimated from earth radiation budget experiment (ERBE) top-of-the-atmosphere (TOA) radiatioe fluxes, global grid point datasets of precipitation, and the climatological analyses of surface evaporation and albedo. Emphasis is placed on the surface evaluation of simulations of landsurface conditions such as surface roughness, surface albedo and the surface wetness factor, and on their effects on surface evaporation, precipitation, layer-cloud and surface temperature. The improvements due to the inclusion of BEST are: a realistic geographical distribution of surface roughness, a decrease in surface albedo over areas with seasonal snow cover, and an increase in surface albedo over snow-free land. The simulated reduction in surface evaporation due, in part, to the biophysical control of vegetation, is also consistent with the previous studies. Since the control climate has a dry bias, the overall simulations from the GCM with BEST are degraded, except for significant improvements for the northern winter hemisphere because of the realistic vegetation-masking effects. The implications of our results for synergistic developments of other aspects of model parameterization schemes such as boundary layer dynamics, clouds, convection and rainfall are discussed.     

Improving the vegetation dynamic simulation in a land surface model by using a statistical-dynamic canopy interception scheme

Canopy interception of incident precipitation, as a critical component of a forest＇s water budget, can affect the amount of water available to the soil, and ultimately vegetation distribution and function. In this paper, a statistical-dynamic approach based on leaf area index and statistical canopy interception is used to parameterize the canopy interception process. The statistical-dynamic canopy interception scheme is implemented into the Community Land Model with dynamic global vegetation model （CLM-DGVM） to improve its dynamic vegetation simulation. The simulation for continental China by the land surface model with the new canopy interception scheme shows that the new one reasonably represents the precipitation intercepted by the canopy. Moreover, the new scheme enhances the water availability in the root zone for vegetation growth, especially in the densely vegetated and semi-arid areas, and improves the model＇s performance of potential vegetation simulation.     

Local climates simulated by two generations of Canadian GCM land surface schemes

Abstract

The performance of two Canadian land surface schemes of widely differing complexity is compared and contrasted in a pair of year‐long simulations using the GCM developed at Atmospheric Environment Service, Canada. The old land surface model incorporates the force‐restore method for soil temperatures and the bucket approximation for soil moisture; the new model, CLASS (Canadian Land Surface Scheme) features three soil layers, an explicitly modelled snow layer, a thermally separate vegetation canopy, and physically‐based calculations of heat and moisture transfers between all of the land surface components and the atmosphere.

It was reported in previous papers that compared with observations, the old scheme tends to generate a climate which is characterized by anomalously high precipitation rates and cold temperatures over land. In this paper, by reference to field measurements and to the energy fluxes and temperatures generated by the two models at local scales, the hypotheses earlier postulated as to the underlying reasons for this are validated. The main factor contributing to the climate anomalies observed with the old scheme is found to be its generation of excessive evaporation rates; this is caused by the fact that the evaporation rate is never directly energy‐limited, the fact that the scaling of the evaporation rale with decreasing soil moisture content underestimates the effect of vegetation stomatal resistance, and the fact that the evaporation rate over bare soil depends not on the surface soil moisture, but on the moisture content of whole modelled soil column. The cold surface temperatures are additionally attributed to systematic errors incurred by the forward‐stepping temperature scheme, and to the fact that soils subjected to subzero temperature forcing in the winter are modelling as freezing completely. Finally, the inability of the old scheme to simulate partially frozen soils means that it proves unable to handle either shallow frost penetration at temperature latitudes, or the development of an active layer in permafrost.     

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

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

Design of land-air parameterization scheme (LAPS) for modelling boundary layer surface processes

Summary A land-air parametrization scheme (LAPS) describes mass, energy and momentum transfer between the land surface and the atmosphere. The scheme is designed as a software package which can be run as part of an atmospheric model or a stand-alone scheme. A single layer approach is chosen for the physical and biophysical scheme background. The scheme has six prognostic variables: two temperatures (one for the canopy vegetation and one for soil surface), one interception storage, and three soil moisture storage variables. The scheme's upper boundary conditions are: air temperature, water vapour pressure, wind speed, radiation and precipitation at some reference level within the atmospheric boundary layer. The sensible and latent heat are calculated using resistance representation. The evaporation from the bare soil is parametrized using the scheme. The soil part is designed as a three-layer model which is used to describe the vertical transfer of water in the soil.The performances of the LAPS scheme were tested using the results of meteorological measurements over a maize field at the experimental site De Sinderhoeve (The Netherlands). The predicted partitioning of the absorbed radiation into sensible and latent heat fluxes is in good agreement with observations. Also, the predicted leaf temperature agrees quite well with the observed values.With 9 Figures     

IAP94陆面过程模式在淮河流域的验证试验

利用ＩＡＰ９４陆面过程模式对淮河流域不同季节（春、夏、秋）、不同下垫面（旱地、森林、水田）进行了单点数值模拟试验，结果表明ＩＡＰ９４不但能较好地模拟出观测的地气间各种能量通量以及植被的冠层温度，同时还很好地模拟出整层土壤含水量的变化趋势，从而证实了ＩＡＰ９４陆面过程模式具有正确描述东亚半湿润季风区陆气相互作用的能力。另外比较结果还揭示了水田表面强烈蒸发作用的重要性，需要在陆过程模式中予以重视。     

陆面模拟中植被辐射传输参数化方案研究

在冠层二流辐射传输模式基础上新发展了一个描述太阳短波辐射在植被中传输的冠层四流辐射传输模式.冠层四流辐射传输模式是在大气辐射传输理论的基础上得到一组描述短波辐射在植被中传输过程的冠层辐射传输基本方程,引进大气中求解辐射传输方程的四流近似解法,并求得冠层四流辐射传输方程的解析解.方程中各项参量能够反映叶子或冠层特殊的几何和光学特征.冠层向上、向下辐射通量取决于冠层散射相函数、叶子在入射光方向投影面积、单个叶子反射率和透射率、叶面积指数以及直射光入射太阳高度角等.四流模式计算叶子水平倾角时对太阳短波辐射的反照率,与二流模式结果比较可以验证模式的理论推导和建模都是正确的:计算结果的比较,表明四流模式在水平叶角分布时计算的冠层反照率与二流模式结果一致,同时直射光从任何太阳高度角入射的冠层反照率结果也一致,从而证明发展的冠层四流辐射传输模式是成功的.模拟试验中将两种模型同时耦合到同一个陆面过程模式中进行比较试验,结果表明,冠层四流辐射传输模式能够得到更精确的植被反照率,从而使得陆面模式计算的地表吸收的净太阳辐射通量更接近于观测值.     

两次暴雨过程模拟对陆面参数化方案的敏感性研究

选取发生在江西和福建境内的两次暴雨个例,利用NCEP再分析资料在对暴雨发生前、后的环境场和物理量场进行诊断和对比分析的基础上,采用中尺度模式WRF V3.3,通过数值模拟探讨了陆面过程对两次暴雨过程的可能影响及其相关的物理过程。结果表明,2012年5月12日江西大暴雨主要受大尺度环流和中尺度天气系统影响,具有范围大、持续时间长等特点,属于大尺度降水为主的暴雨;而2011年8月23日福建暴雨发生在副热带高压控制下的午后,局地下垫面强烈的感热和潜热通量使低层大气不稳定性增强,触发了此次对流性降水为主的暴雨。通过资料诊断分析,可以判断陆面过程对福建暴雨个例的影响程度明显强于江西暴雨个例。通过关闭地表通量试验发现,陆面过程对暴雨模拟十分重要,尤其是对于该个例中对流性降水的发生起到关键性的作用。通过陆面参数化方案的敏感性试验发现,两次暴雨过程对陆面参数化方案均较为敏感。江西暴雨对陆面过程的敏感性主要体现在对流降水的模拟上,而福建暴雨则体现在大尺度降水的模拟方面,即福建暴雨对陆面参数化方案的敏感性强于江西暴雨。敏感性产生机制与降水类型关系紧密,大尺度降水对陆面过程的敏感性主要来源于不同参数化模拟的中高空对流系统的差异,而对流降水的敏感性则与不同参数化模拟的地表通量的差异有关。通过陆面参数的扰动试验进一步发现,相比于地表粗糙度和最小叶孔阻抗,土壤孔隙度和地表反照率则是影响对流降水对陆面过程敏感的关键因子,这在本质上与地表通量是否受到扰动有关。地表通量较风场而言,受扰动引起变化的空间范围广、时间响应快,变化具有明显规律性。所得结果可为深入理解陆面过程影响暴雨等天气过程和改进数值模式对暴雨的模拟能力提供一定的参考。     

Formulation and verification of a land surface parameterization for atmospheric models

The need for a well-defined lower boundary condition for atmospheric numerical models is well documented. This paper describes the formulation of a land surface parameterization, which will be used in atmospheric boundary-layer and mesoscale numerical models. The land surface model has three soil layers for the prediction of soil moisture and soil temperature. Model soil properties depend on soil texture and moisture content. A homogeneous distribution of vegetation is also included, so that transpiration may be included, as well as the interception of precipitation by vegetation elements. The simulated vegetation also affects the mean surface albedo and roughness characteristics.First ISLSCP Field Experiment (FIFE) data are used to verify the model. Three cases during the growing season were chosen, each case having different amounts of vegetation cover. Stand alone simulations, where observations of atmospheric and radiation variables are input to the land surface model, were performed. These simulations show that the model is able to reproduce observed surface energy budgets and surface temperatures reasonably well. The RMS differences between modeled and obsered turbulent fluxes of heat and moisture are quite comparable to those reported by more detailed land surface models.     

The greening of the McGill Paleoclimate Model. Part I: Improved land surface scheme with vegetation dynamics

The formulation of a new land surface scheme (LSS) with vegetation dynamics for coupling to the McGill Paleoclimate Model (MPM) is presented. This LSS has the following notable improvements over the old version: (1) parameterization of deciduous and evergreen trees by using the models climatology and the output of the dynamic global vegetation model, VECODE (Brovkin et al. in Ecological Modelling 101:251–261 (1997), Global Biogeochemical Cycles 16(4):1139, (2002)); (2) parameterization of tree leaf budburst and leaf drop by using the models climatology; (3) parameterization of the seasonal cycle of the grass leaf area index; (4) parameterization of the seasonal cycle of tree leaf area index by using the time-dependent growth of the leaves; (5) calculation of land surface albedo by using vegetation-related parameters, snow depth and the models climatology. The results show considerable improvement of the models simulation of the present-day climate as compared with that simulated in the original physically-based MPM. In particular, the strong seasonality of terrestrial vegetation and the associated land surface albedo variations are in good agreement with several satellite observations of these quantities. The application of this new version of the MPM (the green MPM) to Holocene millennial-scale climate changes is described in a companion paper, Part II.

---