From 30 September to 2 October 1999 a workshop was held in Gif-sur-Yvette, France, with the central objective to develop a research strategy for the next 3–5 years, aiming at a systematic description of root functioning, rooting depth, and root distribution for modeling root water uptake from local and regional to global scales. The goal was to link more closely the weather prediction and climate and hydrological models with ecological and plant physiological information in order to improve the understanding of the impact that root functioning has on the hydrological cycle at various scales. The major outcome of the workshop was a number of recommendations, detailed at the end of this paper, on root water uptake parameterization and modeling and on collection of root and soil hydraulic data.

Abstract

The Project for the Intercomparison of Land-Surface Parameterization Schemes (PILPS) aims to improve understanding and modeling of land surface processes. PILPS phase 2(d) uses a set of meteorological and hydrological data spanning 18 yr (1966–83) from a grassland catchment at the Valdai water-balance research site in Russia. A suite of stand-alone simulations is performed by 21 land surface schemes (LSSs) to explore the LSSs’ sensitivity to downward longwave radiative forcing, timescales of simulated hydrologic variability, and biases resulting from single-year simulations that use recursive spinup. These simulations are the first in PILPS to investigate the performance of LSSs at a site with a well-defined seasonal snow cover and frozen soil. Considerable model scatter for the control simulations exists. However, nearly all the LSS scatter in simulated root-zone soil moisture is contained within the spatial variability observed inside the catchment. In addition, all models show a considerable sensitivity to longwave forcing for the simulation of the snowpack, which during the spring melt affects runoff, meltwater infiltration, and subsequent evapotranspiration. A greater sensitivity of the ablation, compared to the accumulation, of the winter snowpack to the choice of snow parameterization is found. Sensitivity simulations starting at prescribed conditions with no spinup demonstrate that the treatment of frozen soil (moisture) processes can affect the long-term variability of the models. The single-year recursive runs show large biases, compared to the corresponding year of the control run, that can persist through the entire year and underscore the importance of performing multiyear simulations.

Abstract

In the PILPS Phase 2a experiment, 23 land-surface schemes were compared in an off-line control experiment using observed meteorological data from Cabauw, the Netherlands. Two simple sensitivity experiments were also undertaken in which the observed surface air temperature was artificially increased or decreased by 2 K while all other factors remained as observed. On the annual timescale, all schemes show similar responses to these perturbations in latent, sensible heat flux, and other key variables. For the 2-K increase in temperature, surface temperatures and latent heat fluxes all increase while net radiation, sensible heat fluxes, and soil moistures all decrease. The results are reversed for a 2-K temperature decrease. The changes in sensible heat fluxes and, especially, the changes in the latent heat fluxes are not linearly related to the change of temperature. Theoretically, the nonlinear relationship between air temperature and the latent heat flux is evident and due to the convex relationship between air temperature and saturation vapor pressure. A simple test shows that, the effect of the change of air temperature on the atmospheric stratification aside, this nonlinear relationship is shown in the form that the increase of the latent heat flux for a 2-K temperature increase is larger than its decrease for a 2-K temperature decrease. However, the results from the Cabauw sensitivity experiments show that the increase of the latent heat flux in the +2-K experiment is smaller than the decrease of the latent heat flux in the −2-K experiment (we refer to this as the asymmetry). The analysis in this paper shows that this inconsistency between the theoretical relationship and the Cabauw sensitivity experiments results (or the asymmetry) is due to (i) the involvement of the β _g formulation, which is a function of a series stress factors that limited the evaporation and whose values change in the ±2-K experiments, leading to strong modifications of the latent heat flux; (ii) the change of the drag coefficient induced by the changes in stratification due to the imposed air temperature changes (±2 K) in parameterizations of latent heat flux common in current land-surface schemes. Among all stress factors involved in the β _g formulation, the soil moisture stress in the +2-K experiment induced by the increased evaporation is the main factor that contributes to the asymmetry.

Abstract

Twenty-one land surface schemes (LSSs) performed simulations forced by 18 yr of observed meteorological data from a grassland catchment at Valdai, Russia, as part of the Project for the Intercomparison of Land-Surface Parameterization Schemes (PILPS) Phase 2(d). In this paper the authors examine the simulation of snow. In comparison with observations, the models are able to capture the broad features of the snow regime on both an intra- and interannual basis. However, weaknesses in the simulations exist, and early season ablation events are a significant source of model scatter. Over the 18-yr simulation, systematic differences between the models’ snow simulations are evident and reveal specific aspects of snow model parameterization and design as being responsible. Vapor exchange at the snow surface varies widely among the models, ranging from a large net loss to a small net source for the snow season. Snow albedo, fractional snow cover, and their interplay have a large effect on energy available for ablation, with differences among models most evident at low snow depths. The incorporation of the snowpack within an LSS structure affects the method by which snow accesses, as well as utilizes, available energy for ablation. The sensitivity of some models to longwave radiation, the dominant winter radiative flux, is partly due to a stability-induced feedback and the differing abilities of models to exchange turbulent energy with the atmosphere. Results presented in this paper suggest where weaknesses in macroscale snow modeling lie and where both theoretical and observational work should be focused to address these weaknesses.

Abstract

The Protocol for the Analysis of Land Surface Models (PALS) Land Surface Model Benchmarking Evaluation Project (PLUMBER) illustrated the value of prescribing a priori performance targets in model intercomparisons. It showed that the performance of turbulent energy flux predictions from different land surface models, at a broad range of flux tower sites using common evaluation metrics, was on average worse than relatively simple empirical models. For sensible heat fluxes, all land surface models were outperformed by a linear regression against downward shortwave radiation. For latent heat flux, all land surface models were outperformed by a regression against downward shortwave radiation, surface air temperature, and relative humidity. These results are explored here in greater detail and possible causes are investigated. It is examined whether particular metrics or sites unduly influence the collated results, whether results change according to time-scale aggregation, and whether a lack of energy conservation in flux tower data gives the empirical models an unfair advantage in the intercomparison. It is demonstrated that energy conservation in the observational data is not responsible for these results. It is also shown that the partitioning between sensible and latent heat fluxes in LSMs, rather than the calculation of available energy, is the cause of the original findings. Finally, evidence is presented that suggests that the nature of this partitioning problem is likely shared among all contributing LSMs. While a single candidate explanation for why land surface models perform poorly relative to empirical benchmarks in PLUMBER could not be found, multiple possible explanations are excluded and guidance is provided on where future research should focus.

Abstract

In the Project for Intercomparison of Land-Surface Parameterization Schemes phase 2a experiment, meteorological data for the year 1987 from Cabauw, the Netherlands, were used as inputs to 23 land-surface flux schemes designed for use in climate and weather models. Schemes were evaluated by comparing their outputs with long-term measurements of surface sensible heat fluxes into the atmosphere and the ground, and of upward longwave radiation and total net radiative fluxes, and also comparing them with latent heat fluxes derived from a surface energy balance. Tuning of schemes by use of the observed flux data was not permitted. On an annual basis, the predicted surface radiative temperature exhibits a range of 2 K across schemes, consistent with the range of about 10 W m⁻² in predicted surface net radiation. Most modeled values of monthly net radiation differ from the observations by less than the estimated maximum monthly observational error (±10 W m⁻²). However, modeled radiative surface temperature appears to have a systematic positive bias in most schemes; this might be explained by an error in assumed emissivity and by models’ neglect of canopy thermal heterogeneity. Annual means of sensible and latent heat fluxes, into which net radiation is partitioned, have ranges across schemes of30 W m⁻² and 25 W m⁻², respectively. Annual totals of evapotranspiration and runoff, into which the precipitation is partitioned, both have ranges of 315 mm. These ranges in annual heat and water fluxes were approximately halved upon exclusion of the three schemes that have no stomatal resistance under non-water-stressed conditions. Many schemes tend to underestimate latent heat flux and overestimate sensible heat flux in summer, with a reverse tendency in winter. For six schemes, root-mean-square deviations of predictions from monthly observations are less than the estimated upper bounds on observation errors (5 W m⁻² for sensible heat flux and 10 W m⁻² for latent heat flux). Actual runoff at the site is believed to be dominated by vertical drainage to groundwater, but several schemes produced significant amounts of runoff as overland flow or interflow. There is a range across schemes of 184 mm (40% of total pore volume) in the simulated annual mean root-zone soil moisture. Unfortunately, no measurements of soil moisture were available for model evaluation. A theoretical analysis suggested that differences in boundary conditions used in various schemes are not sufficient to explain the large variance in soil moisture. However, many of the extreme values of soil moisture could be explained in terms of the particulars of experimental setup or excessive evapotranspiration.

Abstract

The Rhône-Aggregation (Rhône-AGG) Land Surface Scheme (LSS) intercomparison project is an initiative within the Global Energy and Water Cycle Experiment (GEWEX)/Global Land–Atmosphere System Study (GLASS) panel of the World Climate Research Programme (WCRP). It is a intermediate step leading up to the next phase of the Global Soil Wetness Project (GSWP) (Phase 2), for which there will be a broader investigation of the aggregation between global scales (GSWP-1) and the river scale. This project makes use of the Rhône modeling system, which was developed in recent years by the French research community in order to study the continental water cycle on a regional scale.

The main goals of this study are to investigate how 15 LSSs simulate the water balance for several annual cycles compared to data from a dense observation network consisting of daily discharge from over 145 gauges and daily snow depth from 24 sites, and to examine the impact of changing the spatial scale on the simulations. The overall evapotranspiration, runoff, and monthly change in water storage are similarly simulated by the LSSs, however, the differing partitioning among the fluxes results in very different river discharges and soil moisture equilibrium states. Subgrid runoff is especially important for discharge at the daily timescale and for smaller-scale basins. Also, models using an explicit treatment of the snowpack compared better with the observations than simpler composite schemes.

Results from a series of scaling experiments are examined for which the spatial resolution of the computational grid is decreased to be consistent with large-scale atmospheric models. The impact of upscaling on the domain-averaged hydrological components is similar among most LSSs, with increased evaporation of water intercepted by the canopy and a decrease in surface runoff representing the most consistent inter-LSS responses. A significant finding is that the snow water equivalent is greatly reduced by upscaling in all LSSs but one that explicitly accounts for subgrid-scale orography effects on the atmospheric forcing.

Abstract

The project Land-Use and Climate, Identification of Robust Impacts (LUCID) was conceived to address the robustness of biogeophysical impacts of historical land use–land cover change (LULCC). LUCID used seven atmosphere–land models with a common experimental design to explore those impacts of LULCC that are robust and consistent across the climate models. The biogeophysical impacts of LULCC were also compared to the impact of elevated greenhouse gases and resulting changes in sea surface temperatures and sea ice extent (CO2SST). Focusing the analysis on Eurasia and North America, this study shows that for a number of variables LULCC has an impact of similar magnitude but of an opposite sign, to increased greenhouse gases and warmer oceans. However, the variability among the individual models’ response to LULCC is larger than that found from the increase in CO2SST. The results of the study show that although the dispersion among the models’ response to LULCC is large, there are a number of robust common features shared by all models: the amount of available energy used for turbulent fluxes is consistent between the models and the changes in response to LULCC depend almost linearly on the amount of trees removed. However, less encouraging is the conclusion that there is no consistency among the various models regarding how LULCC affects the partitioning of available energy between latent and sensible heat fluxes at a specific time. The results therefore highlight the urgent need to evaluate land surface models more thoroughly, particularly how they respond to a perturbation in addition to how they simulate an observed average state.

Abstract

The 12 weather and climate models participating in the Global Land–Atmosphere Coupling Experiment (GLACE) show both a wide variation in the strength of land–atmosphere coupling and some intriguing commonalities. In this paper, the causes of variations in coupling strength—both the geographic variations within a given model and the model-to-model differences—are addressed. The ability of soil moisture to affect precipitation is examined in two stages, namely, the ability of the soil moisture to affect evaporation, and the ability of evaporation to affect precipitation. Most of the differences between the models and within a given model are found to be associated with the first stage—an evaporation rate that varies strongly and consistently with soil moisture tends to lead to a higher coupling strength. The first-stage differences reflect identifiable differences in model parameterization and model climate. Intermodel differences in the evaporation–precipitation connection, however, also play a key role.