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.     

UNCERTAINTY IN THE SPECIFICATION OF SURFACE CHARACTERISTICS: A STUDY OF PREDICTION ERRORS IN THE BOUNDARY LAYER

The effects of uncertainty in the specification of surface characteristics on simulated atmospheric boundary layer (ABL) processes and structure were investigated using a one-dimensional soil-vegetation-boundary layer model. Observational data from the First International Satellite Land Surface Climatology Project Field Experiment were selected to quantify prediction errors in simulated boundary-layer parameters. Several numerical 12-hour simulations were performed to simulate the convective boundary-layer structure, starting at 0700 LT 6 June 1987.In the control simulation, measured surface parameters and atmospheric data were used to simulate observed boundary-layer processes. In the remaining simulations, five surface parameters – soil texture, initial soil moisture, minimum stomatal resistance, leaf area index, and vegetation cover – were varied systematically to study how uncertainty in the specification of these surface parameters affects simulated boundary-layer processes.The simulated uncertainty in the specification of these five surface parameters resulted in a wide range of errors in the prediction of turbulent fluxes, mean thermodynamic structure, and the depth of the ABL. Under certain conditions uncertainty in the specifications of soil texture and minimum stomatal resistance had the greatest influence on the boundary-layer structure. A lesser but still moderately strong effect on the simulated ABL resulted from (1) a small decrease (4%) in the observed initial soil moisture (although a large increase [40%] had only a marginal effect), and (2) a large reduction (66%) in the observed vegetation cover. High uncertainty in the specification of leaf area index had only a marginal impact on the simulated ABL. It was also found that the variations in these five surface parameters had a negligible effect on the simulated horizontal wind fields. On the other hand, these variations had a significant effect on the vertical distribution of turbulent heat fluxes, and on the predicted maximum boundary-layer depth, which varied from about 1400–2300 m across the 11 simulations. Thus, uncertainties in the specification of surface parameters can significantly affect the simulated boundary-layer structure in terms of meteorological and air quality model predictions.     

The Role of Planetary Boundary-Layer Parameterizations in the Air Quality of an Urban Area with Complex Topography

The effect of different planetary boundary-layer (PBL) parameterization schemes on the spatial distribution of atmospheric pollution over the complex topography of the greater Athens area is investigated. Four PBL schemes originally implemented in a numerical meteorological model and a fifth one simulating the urban effect are examined. Two different atmospheric conditions are analyzed; a typical summer and a typical winter pollution episode. The relative importance of chemical and physical processes of the pollution predictions is discussed using process analysis. It is revealed that, for primary pollutants, a local scheme seems more adequate to represent the maximum observed concentrations while, completely different in structure, a non-local scheme reproduces the mean observed values in the basin. Concerning secondary pollutants, peak concentration differences, due to the different PBL schemes, are smoothed out. Nevertheless, the PBL scheme selection shapes the horizontal and the vertical extension of maximum values. The non-local and semi non-local schemes are superior to the others, favouring strong vertical mixing and transport towards the surface. The stronger turbulence accommodated effectively by the semi non-local urban scheme enhances ozone production along the sea-breeze axis and preserves the high ozone concentrations during the nighttime hours in the urban core.     

A new design of the LAPS land surface scheme for use over and through heterogeneous and non-heterogeneous surfaces: Numerical simulations and tests

Numerical simulations and tests with the recently redesigned land–air parameterization scheme (LAPS) are presented. In all experiments, supported either by one-point micrometeorological, 1D or 3D simulations, the attention has been directed to: (1) comparison of simulation outputs, expressing the energy transfer over and through heterogeneous and non-heterogeneous surfaces, versus observations and (2) analysis of uncertainties occurring in the solution of the energy balance equation at the land–air interface. To check the proposed method for aggregation of albedo, “propagating hole” sensitivity tests with LAPS over a sandstone rock grid cell have been performed with the forcing meteorological data for July 17, 1999 in Baxter site, Philadelphia (USA). Micrometeorological and biophysical measurements from the surface experiments conducted over crops and apple orchard in Serbia, Poland, Austria and France were used to test the operation of LAPS in calculating surface fluxes and canopy environment temperatures within and above plant covers of different densities. In addition, sensitivity tests with single canopy covers over the Central Europe region and comparison against the observations taken from SYNOP data using 3D simulations were made. Validation of LAPS performances over a solid surface has been done by comparison of 2 m air temperature observations against 5-day simulations over the Sahara Desert rocky ground using 3D model. To examine how realistically the LAPS simulates surface processes over a heterogeneous surface, we compared the air temperature measured at 2 m and that predicted by the 1D model with the LAPS as the surface scheme. Finally, the scheme behaviour over urban surface was tested by runs over different parts of a hypothetical urban area. The corresponding 1D simulations were carried out with an imposed meteorological dataset collected during HAPEX-MOBILHY experiment at Caumont (France). The quantities predicted by the LAPS compare well with the observations and the various subcomponents of the scheme appear to operate realistically.     

Numerical Modelling of the Effects of Vegetation and Environmental Conditions on the Lake Breeze

The influence of vegetation and environmental conditions on the lake breeze and associated boundary-layer turbulence structure has been studied using a two-dimensional nonhydrostatic, compressible mesoscale model coupled with the SiB2 land-surface scheme. The results show that the impacts of vegetation on the lake effects are dependent on the environmental conditions, such as soil wetness and background wind, as well as vegetation characteristics. Both soil wetness and background wind play important roles in modifying lake effects on boundary-layer turbulence and the lake breeze, while the effects of vegetation type are secondary compared to the other factors. Without background wind, and under the same soil wetness, the maximum horizontal windspeed of the lake breeze is insensitive to the type of vegetation. Soil wetness can greatly affect both the maximum horizontal windspeed and the maximum vertical velocities of the lake breeze. With background wind, the lake-breeze circulations, upward motion regions, and boundary-layer turbulence structure all change markedly. A weaker background wind can strengthen the lake breeze, while stronger background wind suppresses the lake breeze circulations. The distribution of sensible and latent heat fluxes is also very sensitive to the soil wetness and background wind. However, for the same soil wetness (0.25 and 0.4 were chosen), there is only a small difference in the distribution of sensible and latent heat fluxes between the bare soil and vegetated soil or between the types of vegetated soils.     

A regional climate model simulation over West Africa: parameterization tests and analysis of land-surface fields

The West African Monsoon has been simulated with the regional climate model PROMES, coupled to the land-surface model ORCHIDEE and nested in ECMWF analysis, within AMMA-EU project. Three different runs are presented to address the influence of changes in two parameterizations (moist convection and radiation) on the simulated West African Monsoon. Another aim of the study is to get an insight into the relationship of simulated precipitation and 2-m temperature with land-surface fluxes. To this effect, data from the AMMA land-surface model intercomparison project (ALMIP) have been used. In ALMIP, offline simulations have been made using the same land-surface model than in the coupled simulation presented here, which makes ALMIP data particularly relevant for the present study, as it enables us to analyse the simulated soil and land-surface fields. The simulation of the monsoon depends clearly on the two analysed parameterizations. The inclusion of shallow convection parametrization affects the intensity of the simulated monsoon precipitation and modifies some dynamical aspects of the monsoon. The use of a fractional cloud-cover parameterization and a more complex radiation scheme is important for better reproducing the amplitude of the latitudinal displacement of the precipitation band. This is associated to an improved simulation of the surface temperature field and the easterly jets. However, the parameterization changes do not affect the timing of the main rainy and break periods of the monsoon. A better representation of downward solar radiation is associated with a smaller bias in the surface heat fluxes. The comparison with ALMIP land-surface and soil fields shows that precipitation and temperature biases in the regional climate model simulation are associated to certain biases in land-surface fluxes. The biases in soil moisture seem to be driven by atmospheric biases as they are strongly affected by the parameterization changes in atmospheric processes.     

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

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

A simple parameterization of sub-grid scale open water for climate models

The effects of small fractions ( < 30%) of open water covering a grid element are currently neglected even in atmospheric general circulation models (AGCMs) which incorporate complex land surface parameterization schemes. Here, a method for simulating sub-grid scale open water is proposed which permits any existing land surface model to be modified to account for open water. This new parameterization is tested as an addition to an advanced land surface scheme and, as expected, is shown to produce general increases in the surface latent heat flux at the expense of the surface sensible heat flux. Small changes in temperature are associated with this change in the partitioning of available energy which is driven by an increase in the wetness of the grid element. The sensitivity of the land surface to increasing amounts of open water is dependent upon the type of vegetation represented. Dense vegetation (with a high leaf area index) is shown to complicate the apparently simple model sensitivity and indicates that more advanced methods of incorporating open water into AGCMs need to be considered and compared against the parameterization suggested here. However, the sensitivity of one land surface model to incorporating open water is large enough to warrant consideration of its incorporation into climate models.     

陆面模式中土壤湿度影响蒸散参数化方案的评估

本文将四个常见陆面模式CLM3.5（Community Land Model Version 3.5）、Noah_LSM（The Noah Land Surface Model）、VIC（Variable Infiltration Capacity）以及SSiB（The Simplified Simple Biosphere Model）中土壤湿度影响蒸散的参数化方案进行简化,并利用实验观测资料对不同参数化方案进行评估,探究不同陆面模式对土壤湿度与蒸散关系的模拟差异,从而为提高模式的模拟能力提供依据。结果表明,（1）CLM与SSiB中计算土壤湿度影响裸土蒸发的参数化方案较Noah_LSM和VIC更接近真实的物理过程,同时CLM与SSiB模式中土壤湿度对蒸发的影响程度较Noah_LSM和VIC大;而对于下垫面有植被条件下的蒸散而言,CLM中包含了植被光合作用、呼吸作用等生物物理学过程,与实际情况更为接近,并且CLM与SSiB中土壤湿度对植被蒸散的影响程度大于VIC,Noah_LSM最低;（2）根据干旱区、半干旱区、半湿润区以及湿润区各站点的分析可知,CLM、SSiB与Noah_LSM中土壤湿度影响蒸散的参数化方案的拟合效果较VIC好,同时在部分站点CLM与SSiB的参数化方案稍优于Noah_LSM。区域之间比较说明,四个模式对干旱半干旱区的模拟效果明显较半湿润区和湿润区好。     

Role of land-surface temperature feedback on model performance for the stable boundary layer

At present a variety of boundary-layer schemes is in use in numerical models and often a large variation of model results is found. This is clear from model intercomparisons, such as organized within the GEWEX Atmospheric Boundary Layer Study (GABLS). In this paper we analyze how the specification of the land-surface temperature affects the results of a boundary-layer scheme, in particular for stable conditions. As such we use a well established column model of the boundary layer and we vary relevant parameters in the turbulence scheme for stable conditions. By doing so, we can reproduce the outcome for a variety of boundary-layer models. This is illustrated with the original set-up of the second GABLS intercomparison study using prescribed geostrophic winds and land-surface temperatures as inspired by (but not identical to) observations of CASES-99 for a period of more than two diurnal cycles. The model runs are repeated using a surface temperature that is calculated with a simple land-surface scheme. In the latter case, it is found that the range of model results in stable conditions is reduced for the sensible heat fluxes, and the profiles of potential temperature and wind speed. However, in the latter case the modelled surface temperatures are rather different than with the original set-up, which also impacts on near-surface air temperature and wind speed. As such it appears that the model results in stable conditions are strongly influenced by non-linear feedbacks in which the magnitude of the geostrophic wind speed and the related land-surface temperature play an important role.     

不同陆面模式对我国地表温度模拟的适用性评估

基于CLDAS大气驱动数据驱动CLM3.5陆面模式和3种不同参数化方案下的Noah-MP陆面模式模拟得到的地表温度,利用中国气象局2009-2013年2000多个国家级地面观测站地表温度进行质量评估。结果表明:从时间分布看,模拟地表温度与观测的偏差及均方根误差均呈季节性波动;从空间分布看,模拟地表温度与观测的偏差及均方根误差在中国东部地区相对于中国西部地区更小。选择Noah-MP陆面模式3种不同参数化方案模拟结果进行对比,结果表明:Noah-MP模式的非动态植被方案不变时,考虑植被覆盖度的二流近似辐射传输方案的Noah-MP陆面模式模拟的地表温度优于考虑太阳高度角和植被三维结构的二流近似辐射传输方案Noah-MP陆面模式模拟的地表温度;选择动态植被方案的Noah-MP陆面模式模拟的地表温度优于选择非动态植被方案的Noah-MP陆面模式;总体而言,考虑动态植被方案的Noah-MP陆面模式模拟的地表温度优于其他两种参数化方案的Noah-MP陆面模式以及CLM3.5陆面模式模拟的地表温度。     

Disaggregation of screen-level variables in a numerical weather prediction model with an explicit simulation of subgrid-scale land-surface heterogeneity

The earth’s surface is characterized by small-scale heterogeneity attributable to variability in land cover, soil characteristics and orography. In atmospheric models, this small-scale variability can be partially accounted for by the so-called mosaic approach, i.e., by computing the land-surface processes on a grid with an explicit higher horizontal resolution than the atmosphere. The mosaic approach does, however, not account for the subgrid-scale variability in the screen-level atmospheric parameters, part of which might be related to land-surface heterogeneity itself. In this study, simulations with the numerical weather prediction model COSMO are shown, employing the mosaic approach together with a spatial disaggregation of the atmospheric forcing by the screen-level variables to the subgrid-scale. The atmospheric model is run with a 2.8?km horizontal grid resolution while the land surface processes are computed on a 400-m horizontal grid. The disaggregation of the driving atmospheric variables at screen-level is achieved by a three-step statistical downscaling with rules learnt from high-resolution fully coupled COSMO simulations, where both, atmosphere and surface, were simulated on a 400-m grid. The steps encompass spline interpolation of the grid scale variables, conditional regression based on the high-resolution runs, and an optional stochastic noise generator which restores the variability of the downscaled variables. Simulations for a number of case studies have been carried out, with or without mosaic surface representation and with or without atmospheric disaggregation, and evaluated with respect to the surface state variables and the turbulent surface exchange fluxes of sensible and latent heat. The results are compared with the high-resolution fully coupled COSMO simulations. The results clearly demonstrate the high importance of accounting for subgrid-scale surface heterogeneity. It is shown that the atmospheric disaggregation leads to clear additional improvements in the structures of the two-dimensional surface state variable fields, but to only marginally impacts on the simulation of the turbulent surface exchange fluxes. A detailed analysis of these results identifies strongly correlated errors in atmospheric and surface variables in the mosaic approach as the main reason for the latter. The effects of these errors largely cancel out in the flux parameterization, and thus explain the comparably good results for the fluxes in the mosaic approach without atmospheric disaggregation despite inferior performance for the surface state variables themselves. Inserting noise in the disaggregation scheme leads to a deterioration of the results.     

A comparison of boundary layer diffusion schemes in unstable conditions over land

We compare the results of a local and a nonlocal scheme for vertical diffusion in the atmospheric boundary layer with observations at the 200 m tower at Cabauw. This is done for a 12 h period during daytime on 31 May 1978, which is characterised by strong insolation, clear skies, moderately strong winds and weak advection. The local diffusion scheme uses an eddy diffusivity determined independently at each point along the vertical based on local vertical gradients of wind and virtual potential temperature, similar to the usual approach in atmospheric models. The nonlocal scheme determines an eddy diffusivity profile based on a diagnosed boundary-layer height and a turbulent velocity scale. It also incorporates nonlocal (vertical) transport effects for heat and moisture. The boundary-layer diffusion schemes are forced with the locally observed fluxes for heat and moisture. The outputs of the scheme are compared with the observed mean structure along the Cabauw tower, and the radiosonde profile at a nearby location (De Bilt). Overall, the nonlocal scheme transports moisture away from the surface more rapidly than the local scheme, and deposits the moisture at higher levels. The local scheme tends to saturate the lowest model levels unrealistically in comparison with the observations. We also compare the outputs of the two diffusion schemes with the results of a transilient model simulation. Subsequently, we study the impact on the model behaviour by varying important parameters in both diffusion schemes and we investigate the sensitivity to uncertainty in the environmental conditions. Finally, we study the interaction of the diffusion schemes with a simple surface flux scheme.     

数值模式中的大气边界层参数化方案综述

大气边界层对下垫面和自由大气之间的热量、动量和物质交换以及全球辐射的收支平衡有着非常重要的作用,数值模式中对边界层的模拟需要通过参数化方案表现出来。综述了数值模式中常用的边界层参数化方案,从闭合框架以及边界层内的重要物理过程这两个方面对各类方案进行了总结,并讨论了边界层参数化方案的发展方向。     

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.     

Uncertainties linked to land-surface processes in climate change simulations

 The impact of climate change on the hydrology of continental surfaces is critical for human activities but the response of the surface to this perturbation may also affect the sensitivity of the climate. This complex feedback is simulated in general circulation models (GCMs) used for climate change predictions by their land-surface schemes. The present study attempts to quantify the uncertainty associated with these schemes and what impact it has on our confidence in the simulated climate anomalies. Four GCMs, each coupled to two different land-surface schemes, are used to explore the spectrum of uncertainties. It is shown that, in this sample, surface processes have a significant contribution to our ability to predict surface temperature changes and perturbations of the hydrological cycle in an environment with doubled greenhouse gas concentration. The results reveal that the uncertainty introduced by land-surface processes in the simulated climate is different from its impact on the sensitivity of GCMs to climate change, indeed an alteration of the surface parametrization with little impact on model climate can affect sensitivity significantly. This result leads us to believe that the validation of land-surface schemes should not be limited to the current climate but should also cover their sensitivity to variations in climatic forcing. Received: 24 June 1999 / Accepted: 20 April 2000     

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.     

一个陆面过程参数化模式与 MM5的耦合

在法国陆面过程模式的基础上,为了表示冠层叶片遮挡对水分蒸发阻抗的影响,在植被覆盖部分引入了遮盖因子,然后将这个修正的陆面过程参数化模式耦合到MM5模式中。耦合后的模式模拟出了因为降水造成土壤湿度变化和植被覆盖动态作用对地面通量的影响,而原MM5模式模拟结果则没有反映上述动态变化对地面通量的作用。原MM5模式和耦合模式模拟了1993年8月17日到20日以内蒙古半干旱草原为中心的中尺度区域的气象场,模拟结果和IMGRASS预观测资料进行了对比,对比说明新的陆面过程模式提高了MM5模式对地面通量和边界层各物理量（风、温、湿）的模拟精度。     

用于气候研究的雪盖模型参数化方案敏感性研究

为了得到一个适用于气候研究简化的季节性雪盖模式最佳方案,必须对雪盖内部的重要物理过程、其与上大气相互作用、相应模型的参数化方案和有关的参数选取以及模型的分层结构进行深入研究.利用作者的雪盖模型(SAST),对其中的一些关键性过程的有关参数化方案(如压实、相变、融化雪水流动及分层方案考虑等)及关键的参数(如雪面反照率、有效热传导系数及持水能力等)进行了分析和敏感性试验,得到若干有意义的结论,为雪盖模式改进提供有用的结论.