To prepare for the next-generation North American Land Data Assimilation System (NLDAS), three advanced land surface models [LSMs; i.e., Community Land Model, version 4.0 (CLM4.0); Noah LSM with multiphysics options (Noah-MP); and Catchment LSM-Fortuna 2.5 (CLSM-F2.5)] were run for the 1979-2014 period within the NLDAS-based framework. Unlike the LSMs currently executing in the operational NLDAS, these three advanced LSMs each include a groundwater component. In this study, the model simulations of monthly terrestrial water storage anomaly (TWSA) and its individual water storage components are evaluated against satellite-based and in situ observations, as well as against reference reanalysis products, at basinwide and statewide scales. The quality of these TWSA simulations will contribute to determining the suitability of these models for the next phase of the NLDAS. Overall, it is found that all three models are able to reasonably capture the monthly and interannual variability and magnitudes of TWSA. However, the relative contributions of the individual water storage components to TWSA are very dependent on the model and basin. A major contributor to the TWSA is the anomaly of total column soil moisture content for CLM4.0 and Noah-MP, while the groundwater storage anomaly is the major contributor for CLSM-F2.5. Other water storage components such as the anomaly of snow water equivalent also play a role in all three models. For each individual water storage component, the models are able to capture broad features such as monthly and interannual variability. However, there are large intermodel differences and quantitative uncertainties, which are motivating follow-on investigations in the NLDAS Science Testbed developed by the NASA and NCEP NLDAS teams.

This paper evaluates multi-decadal simulations of the UCLA version of Climate Forecast System version 2, in which the default Noah land surface model has been replaced with the Simplified Simple Biosphere Model version-2. To examine the influence of the atmosphere–ocean (AO) interaction on the variability, two different simulations were conducted: one with interactive ocean component, and the other constrained by the prescribed sea surface temperature. We evaluate the mean seasonal climatology of precipitation and temperature, along with the model’s ability to reproduce atmospheric variability at different scales over the globe, including extratropical modes of atmospheric variability, and long-term trends of global and hemispheric temperature and regional precipitation. Here, we particularly selected two monsoon regions, East Asia and West Africa, where the simulation of multi-decadal variations which has heretofore been a challenging task, to examine decadal variation of monsoon precipitation. In general, temperature anomaly trends were better captured than those of precipitation in both simulations. Results suggest that the AO interaction, represented as latent heat flux, contributes to improve reproducibility of global-wide climatology, extratropical modes of atmospheric variability, and variability in the multi-decadal climate simulation, as well as for inter-decadal variability of the East Asian summer monsoon.

We tested four land surface parameterization schemes against long-term (5 months) area-averaged observations over the 15 km × 15 km First International Satellite Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE) area. This approach proved to be very beneficial to understanding the performance and limitations of different land surface models. These four surface models, embodying different complexities of the evaporation/hydrology treatment, included the traditional simple bucket model, the simple water balance (SWB) model, the Oregon State University (OSU) model, and the simplified Simple Biosphere (SSiB) model. The bucket model overestimated the evaporation during wet periods, and this resulted in unrealistically large negative sensible heat fluxes. The SWB model, despite its simple evaporation formulation, simulated well the evaporation during wet periods, but it tended to underestimate the evaporation during dry periods. Overall, the OSU model ably simulated the observed seasonal and diurnal variation in evaporation, soil moisture, sensible heat flux, and surface skin temperature. The more complex SSiB model performed similarly to the OSU model. A range of sensitivity experiments showed that some complexity in the canopy resistance scheme is important in reducing both the overestimation of evaporation during wet periods and underestimation during dry periods. Properly parameterizing not only the effect of soil moisture stress but also other canopy resistance factors, such as the vapor pressure deficit stress, is critical for canopy resistance evaluation. An overly simple-canopy resistance that includes only soil moisture stress is unable to simulate observed surface evaporation during dry periods. Given a modestly comprehensive time-dependent canopy resistance treatment, a rather simple surface model such as the OSU model can provide good area-averaged surface heat fluxes for mesoscale atmospheric models. Copyright 1996 by the American Geophysical Union.

Infiltration in soils is a key process that partitions precipitation at the land surface into surface runoff and water that enters the soil profile. We reviewed the basic principles of water infiltration in soils and we analyzed approaches commonly used in land surface models (LSMs) to quantify infiltration as well as its numerical implementation and sensitivity to model parameters. We reviewed methods to upscale infiltration from the point to the field, hillslope, and grid cell scales of LSMs. Despite the progress that has been made, upscaling of local-scale infiltration processes to the grid scale used in LSMs is still far from being treated rigorously. We still lack a consistent theoretical framework to predict effective fluxes and parameters that control infiltration in LSMs. Our analysis shows that there is a large variety of approaches used to estimate soil hydraulic properties. Novel, highly resolved soil information at higher resolutions than the grid scale of LSMs may help in better quantifying subgrid variability of key infiltration parameters. Currently, only a few LSMs consider the impact of soil structure on soil hydraulic properties. Finally, we identified several processes not yet considered in LSMs that are known to strongly influence infiltration. Especially, the impact of soil structure on infiltration requires further research. To tackle these challenges and integrate current knowledge on soil processes affecting infiltration processes into LSMs, we advocate a stronger exchange and scientific interaction between the soil and the land surface modeling communities.

The water-balance components of 16 Soil-Vegetation Atmospheric Transfer (SVAT) schemes were evaluated by comparing predicted and observed streamflow, predicted evapotranspiration and evapotranspiration inferred from an atmospheric moisture budget analysis, and soil moisture storage changes for a seven-year period (1980-1986) using data from the Red-Arkansas River basins of the Southern Great Plains of the USA. The evaluations support the following suggestions: (a) The mean annual runoff of all models follows, at least generally, the strong climatic East-West gradient of precipitation, although most models predict too much runoff in the dry part of the basin. (b) The mean monthly storage change tends to be underestimated, even though all models capture reasonably well the seasonality of the evapotranspiration. (c) The wide range of conceptualizations used for generation of surface and subsurface runoff strongly affect runoff generation on seasonal, and shorter, time scales. Model responses to summer precipitation ranged from almost no summer runoff (one model) to the (more common) situation of persistent overprediction of summer runoff, especially in the driest part of the basin. (d) All models tended to underpredict evapotranspiration in summer and overpredict in winter. (e) Model-derived mean seasonal cycles of changes in soil moisture storage are qualitatively similar to those inferred from observations, but most models do not predict the decrease in April soil moisture storage and the increase in October that is inferred from observations.

The energy components of sixteen Soil-Vegetation Atmospheric Transfer (SVAT) schemes were analyzed and intercompared using 10 years of surface meteorological and radiative forcing data from the Red-Arkansas River basin in the Southern Great Plains of the United States. Comparisons of simulated surface energy fluxes among models showed that the net radiation and surface temperature generally had the best agreement among the schemes. On an average (annual and monthly) basis, the estimated latent heat fluxes agreed (to within approximate estimation errors) with the latent heat fluxes derived from a radiosonde-based atmospheric budget method for slightly more than half of the schemes. The sensible heat fluxes had larger differences among the schemes than did the latent heat fluxes, and the model-simulated ground heat fluxes had large variations among the schemes. The spatial patterns of the model-computed net radiation and surface temperature were generally similar among the schemes, and appear reasonable and consistent with observations of related variables, such as surface air temperature. The spatial mean patterns of latent and sensible heat fluxes were less similar than for net radiation, and the spatial patterns of the ground heat flux vary greatly among the 16 schemes. Generally, there is less similarity among the models in the temporal (interannual) variability of surface fluxes and temperature than there is in the mean fields, even for schemes with similar mean fields.

Sixteen land-surface schemes participating in the project for the Intercomparison of Land-surface Schemes (PILPS) Phase 2(c) were run using 10 years (1979-1988) of forcing data for the Red-Arkansas River basins in the Southern Great Plains region of the United States. Forcing data (precipitation, incoming radiation and surface meteorology) and land-surface characteristics (soil and vegetation parameters) were provided to each of the participating schemes. Two groups of runs are presented. (1) Calibration-validation runs, using data from six small catchments distributed across the modeling domain. These runs were designed to test the ability of the schemes to transfer information about model parameters to other catchments and to the computational grid boxes. (2) Base-runs, using data for 1979-1988, designed to evaluate the ability of the schemes to reproduce measured energy and water fluxes over multiple seasonal cycles across a climatically diverse, continental-scale basin. All schemes completed the base-runs but five schemes chose not to calibrate. Observational data (from 1980-1986) including daily river flows and monthly basin total evaporation estimated through an atmospheric budget analysis, were used to evaluate model performance. In general, the results are consistent with earlier PILPS experiments in terms of differences among models in predicted water and energy fluxes. The mean annual net radiation varied between 80 and 105 W m-2 (excluding one model). The mean annual Bowen ratio varied from 0.52 to 1.73 (also excluding one model) as compared to the data-estimated value of 0.92. The run-off ratios varied from a low of 0.02 to a high of 0.41, as compared to an observed value of 0.15. In general, those schemes that did not calibrate performed worse, not only on the validation catchments, but also at the scale of the entire modeling domain. This suggests that further PILPS experiments on the value of calibration need to be carried out.

Abstract: The prediction skill for precipitation anomalies in late spring and summer months—a significant component of extreme climate events—has remained stubbornly low for years. This paper presents a new idea that utilizes information on boreal spring land surface temperature/subsurface temperature (LST/SUBT) anomalies over the Tibetan Plateau (TP) to improve prediction of subsequent summer droughts/floods over several regions over the world, East Asia and North America in particular. The work was performed in the framework of the GEWEX/LS4P Phase I (LS4P-I) experiment, which focused on whether the TP LST/SUBT provides an additional source for subseasonal-to-seasonal (S2S) predictability. The summer 2003, when there were severe drought/flood over the southern/northern part of the Yangtze River basin, respectively, has been selected as the focus case. With the newly developed LST/SUBT initialization method, the observed surface temperature anomaly over the TP has been partially produced by the LS4P-I model ensemble mean, and 8 hotspot regions in the world were identified where June precipitation is significantly associated with anomalies of May TP land temperature. Consideration of the TP LST/SUBT effect has produced about 25–50% of observed precipitation anomalies in most hotspot regions. The multiple models have shown more consistency in the hotspot regions along the Tibetan Plateau-Rocky Mountain Circumglobal (TRC) wave train. The mechanisms for the LST/SUBT effect on the 2003 drought over the southern part of the Yangtze River Basin are discussed. For comparison, the global SST effect has also been tested and 6 regions with significant SST effects were identified in the 2003 case, explaining about 25–50% of precipitation anomalies over most of these regions. This study suggests that the TP LST/SUBT effect is a first-order source of S2S precipitation predictability, and hence it is comparable to that of the SST effect. With the completion of the LS4P-I, the LS4P-II has been launched and the LS4P-II protocol is briefly presented.

Subseasonal-to-seasonal (S2S) precipitation prediction in boreal spring and summer months, which contains a significant number of high-signal events, is scientifically challenging and prediction skill has remained poor for years. Tibetan Plateau (TP) spring observed surface temperatures show a lag correlation with summer precipitation in several remote regions, but current global land-atmosphere coupled models are unable to represent this behavior due to significant errors in producing observed TP surface temperatures. To address these issues, the Global Energy and Water Exchanges (GEWEX) program launched the "Impact of Initialized Land Temperature and Snowpack on Subseasonal-to-Seasonal Prediction"(LS4P) initiative as a community effort to test the impact of land temperature in high-mountain regions on S2S prediction by climate models: more than 40 institutions worldwide are participating in this project. After using an innovative new land state initialization approach based on observed surface 2-m temperature over the TP in the LS4P experiment, results from a multimodel ensemble provide evidence for a causal relationship in the observed association between the Plateau spring land temperature and summer precipitation over several regions across the world through teleconnections. The influence is underscored by an out-of-phase oscillation between the TP and Rocky Mountain surface temperatures. This study reveals for the first time that high-mountain land temperature could be a substantial source of S2S precipitation predictability, and its effect is probably as large as ocean surface temperature over global "hotspot"regions identified here; the ensemble means in some "hotspots"produce more than 40% of the observed anomalies. This LS4P approach should stimulate more follow-on explorations.

Abstract. Sub-seasonal to seasonal (S2S) prediction, especially the prediction of extreme hydroclimate events such as droughts and floods, is not only scientifically challenging but has substantial societal impacts. Motivated by preliminary studies, the Global Energy and Water Exchanges (GEWEX)/Global Atmospheric System Study (GASS) has launched a new initiative called Impact of initialized Land Surface temperature and Snowpack on Sub-seasonal to Seasonal Prediction (LS4P), as the first international grass-root effort to introduce spring land surface temperature (LST)/subsurface temperature (SUBT) anomalies over high mountain areas as a crucial factor that can lead to significant improvement in precipitation prediction through the remote effects of land/atmosphere interactions. LS4P focuses on process understanding and predictability, hence it is different from, and complements, other international projects that focus on the operational S2S prediction. More than forty groups worldwide have participated in this effort, including 21 Earth System Models, 9 regional climate models, and 7 data groups. This paper overviews the history and objectives of LS4P, provides the first phase experimental protocol (LS4P-I) which focuses on the remote effect of the Tibetan Plateau, discusses the LST/SUBT initialization, and presents the preliminary results. Multi-model ensemble experiments and analyses of observational data have revealed that the hydroclimatic effect of the spring LST in the Tibetan Plateau is not limited to the Yangtze River basin but may have a significant large-scale impact on summer precipitation and its S2S prediction. LS4P models are unable to preserve the initialized LST anomalies in producing the observed anomalies largely for two main reasons: i) inadequacies in the land models arising from total soil depths which are too shallow and the use of simplified parameterizations which both tend to limit the soil memory; and ii) reanalysis data, that are used for initial conditions, have large discrepancies from the observed mean state and anomalies of LST over the Tibetan Plateau. Innovative approaches have been developed to largely overcome these problems.