Pedro Viterbo

The land surface-atmosphere interaction : A review based on observational and global modeling perspectives

This review discusses the land-surface-atmosphere interaction using observations from two North American field experiments (First International Satellite Land Surface Climatology Project Field Experiment (FIFE) and Boreal Ecosystem Atmosphere Study (BOREAS)) and the application of research data to the improvement of land surface and boundary layer parameterizations in the European Centre for Medium-Range Weather Forecast (ECMWF) global forecast model. Using field data, we discuss some of the diurnal and seasonal feedback loops controlling the net surface radiation and its partition into the surface sensible and latent heat fluxes and the ground heat flux. We consider the impact on the boundary layer evolution and show the changes in the diurnal cycle with soil moisture in midsummer. We contrast the surface energy budget over the tropical oceans with that over both dry and wet land surfaces in summer. Results from a new ECMWF model with four predicted soil layers illustrate the interaction between the soil moisture reservoir, evaporation and precipitation on different timescales and space scales. An analysis of an ensemble of 30-day integrations for July 1993 (the month of the Mississippi flood) showed a large sensitivity of the monthly precipitation pattern (and amount) to different initial soil moisture conditions. Short-range forecasts with old and new land surface and boundary layer schemes showed that the new scheme produced much better precipitation forecasts for the central United States because of a more realistic thermodynamic structure, which in turn resulted from improved evaporation in an area that is about 1-day upstream. The results suggest that some predictability exists in the extended range as a result of the memory of the soil moisture reservoir. We also discuss briefly the problem of soil moisture initialization in a global forecast model and summarize recent experience with nudging of soil moisture at ECMWF and improvements in the surface energy budget coming from the better prediction of clouds.

An improved land surface parameterization scheme in the ECMWF model and its validation

Abstract A new version of the ECMWF land surface parameterization scheme is described. It has four prognostic layers in the soil for temperature and soil moisture, with a free drainage and a zero heat flux condition at the bottom as a boundary condition. The scheme has been extensively tested in stand-alone mode with the help of long observational time series from three different experiments with different climatological regimes: the First ISLSCP (International Satellite Land Surface Climatology Project) Field Experiment in the United States, Cabauw in the Netherlands, and the Amazonian Rainforest Meteorological Experiment in Brazil. The emphasis is on seasonal timescales because it was felt that the main deficiencies in the old ECMWF land surface scheme were related to its capability of storing precipitation in spring and making it available for evaporation later in the year. It is argued that the stand-alone testing is particularly important, because it allows one to isolate problems in the land surface...

The Anomalous Rainfall over the United States during July 1993: Sensitivity to Land Surface Parameterization and Soil Moisture Anomalies

Abstract This paper discusses the sensitivity of short- and medium-range precipitation forecasts for the central United States to land surface parametrization and soil moisture anomalies. Two forecast systems with different land surface and boundary layer schemes were running in parallel during the extreme rainfall events of July 1993. One forecast system produces much better precipitation forecasts due to a more realistic thermodynamic structure resulting from improved evaporation in an area that is about 1 day upstream from the area of heaviest rain. The paper also discusses two ensembles of 30-day integrations for July 1993. In the first ensemble, soil moisture is initialized at field capacity (100% availability); in the second ensemble it is at 25% of soil moisture availability. It is shown that the moist integrations produce a much more realistic precipitation pattern than the dry integrations. These results suggest that there may be some predictive skill in the monthly range related to the time-scal...

A Revised Hydrology for the ECMWF Model: Verification from Field Site to Terrestrial Water Storage and Impact in the Integrated Forecast System

Abstract The Tiled ECMWF Scheme for Surface Exchanges over Land (TESSEL) is used operationally in the Integrated Forecast System (IFS) for describing the evolution of soil, vegetation, and snow over the continents at diverse spatial resolutions. A revised land surface hydrology (H-TESSEL) is introduced in the ECMWF operational model to address shortcomings of the land surface scheme, specifically the lack of surface runoff and the choice of a global uniform soil texture. New infiltration and runoff schemes are introduced with a dependency on the soil texture and standard deviation of orography. A set of experiments in stand-alone mode is used to assess the improved prediction of soil moisture at the local scale against field site observations. Comparison with basin-scale water balance (BSWB) and Global Runoff Data Centre (GRDC) datasets indicates a consistently larger dynamical range of land water mass over large continental areas and an improved prediction of river runoff, while the effect on atmospheric...

Basin scale estimates of evapotranspiration using GRACE and other observations

Evapotranspiration is integral to studies of the Earth system, yet it is difficult to measure on regional scales. One estimation technique is a terrestrial water budget, i.e., total precipitation minus the sum of evapotranspiration and net runoff equals the change in water storage. Gravity Recovery and Climate Experiment (GRACE) satellite gravity observations are now enabling closure of this equation by providing the terrestrial water storage change. Equations are presented here for estimating evapotranspiration using observation based information, taking into account the unique nature of GRACE observations. GRACE water storage changes are first substantiated by comparing with results from a land surface model and a combined atmospheric-terrestrial water budget approach. Evapotranspiration is then estimated for 14 time periods over the Mississippi River basin and compared with output from three modeling systems. The GRACE estimates generally lay in the middle of the models and may provide skill in evaluating modeled evapotranspiration.

Multimodel estimate of the global terrestrial water balance: setup and first results

Six land surface models and five global hydrological models participate in a model intercomparison project [Water Model Intercomparison Project (WaterMIP)], which for the first time compares simulation results of these different classes of models in a consistent way. In this paper, the simulation setup is described and aspects of the multimodel global terrestrial water balance are presented. All models were run at 0.58 spatial resolution for the global land areas for a 15-yr period (1985–99) using a newly developed global meteorological dataset. Simulated global terrestrial evapotranspiration, excluding Greenland and Antarctica, ranges from 415 to 586 mm yr 21 (from 60 000 to 85 000 km 3 yr 21 ), and simulated runoff ranges from 290 to 457 mm yr 21 (from 42 000 to 66 000 km 3 yr 21 ). Both the mean and median runoff fractions for the land surface models are lower than those of the global hydrological models, although the range is wider. Significant simulation differences between land surface and global hydrological models are found to be caused by the snow scheme employed. The physically based energy balance approach used by land surface models generally results in lower snow water equivalent values than the conceptual degreeday approach used by global hydrological models. Some differences in simulated runoff and evapotranspiration are explained by model parameterizations, although the processes included and parameterizations used are not distinct to either land surface models or global hydrological models. The results show that differences between models are a major source of uncertainty. Climate change impact studies thus need to use not only multiple climate models but also some other measure of uncertainty (e.g., multiple impact

The WFDEI meteorological forcing data set: WATCH Forcing Data methodology applied to ERA-Interim reanalysis data

The WFDEI meteorological forcing data set has been generated using the same methodology as the widely used WATCH Forcing Data (WFD) by making use of the ERA-Interim reanalysis data. We discuss the specifics of how changes in the reanalysis and processing have led to improvement over the WFD. We attribute improvements in precipitation and wind speed to the latest reanalysis basis data and improved downward shortwave fluxes to the changes in the aerosol corrections. Covering 1979–2012, the WFDEI will allow more thorough comparisons of hydrological and Earth System model outputs with hydrologically and phenologically relevant satellite products than using the WFD.

Evaluation of forest snow processes models (SnowMIP2)

Thirty-three snowpack models of varying complexity and purpose were evaluated across a wide range of hydrometeorological and forest canopy conditions at five Northern Hemisphere locations, for up t ...

The Rhône-Aggregation Land Surface Scheme Intercomparison Project: An Overview

The Rhone-Aggregation (Rhone-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 Rhone 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.

Simulation of high-latitude hydrological processes in the Torne-Kalix basin: PILPS Phase 2(e) 1: Experiment description and summary intercomparisons

Abstract Twenty-one land-surface schemes (LSSs) participated in the Project for Intercomparison of Land-surface Parameterizations (PILPS) Phase 2(e) experiment, which used data from the Torne–Kalix Rivers in northern Scandinavia. Atmospheric forcing data (precipitation, air temperature, specific humidity, wind speed, downward shortwave and longwave radiation) for a 20-year period (1979–1998) were provided to the 21 participating modeling groups for 218 1/4° grid cells that represented the study domain. The first decade (1979–1988) of the period was used for model spin-up. The quality of meteorologic forcing variables is of particular concern in high-latitude experiments and the quality of the gridded dataset was assessed to the extent possible. The lack of sub-daily precipitation, underestimation of true precipitation and the necessity to estimate incoming solar radiation were the primary data concerns for this study. The results from two of the three types of runs are analyzed in this, the first of a three-part paper: (1) calibration–validation runs—calibration of model parameters using observed streamflow was allowed for two small catchments (570 and 1300 km2), and parameters were then transferred to two other catchments of roughly similar size (2600 and 1500 km2) to assess the ability of models to represent ungauged areas elsewhere; and 2) reruns—using revised forcing data (to resolve problems with apparent underestimation of solar radiation of approximately 36%, and certain other problems with surface wind in the original forcing data). Model results for the period 1989–1998 are used to evaluate the performance of the participating land-surface schemes in a context that allows exploration of their ability to capture key processes spatially. In general, the experiment demonstrated that many of the LSSs are able to capture the limitations imposed on annual latent heat by the small net radiation available in this high-latitude environment. Simulated annual average net radiation varied between 16 and 40 W/m2 for the 21 models, and latent heat varied between 18 and 36 W/m2. Among-model differences in winter latent heat due to the treatment of aerodynamic resistance appear to be at least as important as those attributable to the treatment of canopy interception. In many models, the small annual net radiation forced negative sensible heat on average, which varied among the models between −11 and 9 W/m2. Even though the largest evaporation rates occur in the summer (June, July and August), model-predicted snow sublimation in winter has proportionately more influence on differences in annual runoff volume among the models. A calibration experiment for four small sub-catchments of the Torne–Kalix basin showed that model parameters that are typically adjusted during calibration, those that control storage of moisture in the soil column or on the land surface via ponding, influence the seasonal distribution of runoff, but have relatively little impact on annual runoff ratios. Similarly, there was no relationship between annual runoff ratios and the proportion of surface and subsurface discharge for the basin as a whole.