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.

A package of changes will be described for the 32-km 45-layer Mesoscale Eta model (also known as the "Early" Eta at 00z and 12z and hereafter referred to as simply the Meso Eta). While it had been running there since November, the Meso Eta was officially implemented on the new NCEP IBM SP supercomputer in January 2000. The model is run 4 times per day at 00Z, 06Z, 12Z, and 18Z with identical resolution. Each run now goes out to 48 hours and each run is initialized from the fully cycled Eta Data Assimilation System (EDAS) of 12 hour duration. AWIPS output for the off-time (06 and 18Z) runs is only available out to 33 and 30 hours, respectively, due to the previous time format of those runs. (The Eta was run at 29 km/45 level resolution at 03Z and 18Z out to 33 and 30 hours, respectively, and those runs were uniquely referred to as the "Meso Eta."). An overview of the Meso Eta can be found here.

This study examines the performance of the NCEP Global Forecast System (GFS) surface layer parameterization scheme for strongly stable conditions over land in which turbulence is weak or even disappears because of high near-surface atmospheric stability. Cases of both deep snowpack and snow-free conditions are investigated. The results show that decoupling and excessive near-surface cooling may appear in the late afternoon and nighttime, manifesting as a severe cold bias of the 2-m surface air temperature that persists for several hours or more. Concurrently, because of negligible downward heat transport from the atmosphere to the land, a warm temperature bias develops at the first model level. The authors test changes to the stable surface layer scheme that include introduction of a stability parameter constraint that prevents the land–atmosphere system from fully decoupling and modification to the roughness-length formulation. GFS sensitivity runs with these two changes demonstrate the ability of the proposed surface layer changes to reduce the excessive near-surface cooling in forecasts of 2-m surface air temperature. The proposed changes prevent both the collapse of turbulence in the stable surface layer over land and the possibility of numerical instability resulting from thermal decoupling between the atmosphere and the surface. The authors also execute and evaluate daily GFS 7-day test forecasts with the proposed changes spanning a one-month period in winter. The assessment reveals that the systematic deficiencies and substantial errors in GFS near-surface 2-m air temperature forecasts are considerably reduced, along with a notable reduction of temperature errors throughout the lower atmosphere and improvement of forecast skill scores for light and medium precipitation amounts.

The purpose of this study is to evaluate the components of the land surface water budget in the four land surface models (Noah, SAC‐Sacramento Soil Moisture Accounting Model, (VIC) Variable Infiltration Capacity Model, and Mosaic) applied in the newly implemented National Centers for Environmental Prediction (NCEP) operational and research versions of the North American Land Data Assimilation System version 2 (NLDAS‐2). This work focuses on monthly and annual components of the water budget over 12 National Weather Service (NWS) River Forecast Centers (RFCs). Monthly gridded FLUX Network (FLUXNET) evapotranspiration (ET) from the Max‐Planck Institute (MPI) of Germany, U.S. Geological Survey (USGS) total runoff ( Q ), changes in total water storage (d S /d t , derived as a residual by utilizing MPI ET and USGS Q in the water balance equation), and Gravity Recovery and Climate Experiment (GRACE) observed total water storage anomaly (TWSA) and change (TWSC) are used as reference data sets. Compared to these ET and Q benchmarks, Mosaic and SAC (Noah and VIC) in the operational NLDAS‐2 overestimate (underestimate) mean annual reference ET and underestimate (overestimate) mean annual reference Q . The multimodel ensemble mean (MME) is closer to the mean annual reference ET and Q . An anomaly correlation (AC) analysis shows good AC values for simulated monthly mean Q and d S /d t but significantly smaller AC values for simulated ET. Upgraded versions of the models utilized in the research side of NLDAS‐2 yield largely improved performance in the simulation of these mean annual and monthly water component diagnostics. These results demonstrate that the three intertwined efforts of improving (1) the scientific understanding of parameterization of land surface processes, (2) the spatial and temporal extent of systematic validation of land surface processes, and (3) the engineering‐oriented aspects such as parameter calibration and optimization are key to substantially improving product quality in various land data assimilation systems.

The NCEP Climate Forecast System Reanalysis (CFSR) uses the NASA Land Information System (US) to create its land surface analysis: the NCEP Global Land Data Assimilation System (GLDAS). Comparing to the previous two generations of NCEP global reanalyses, this is the first time a coupled land-atmosphere data assimilation system is included in a global reanalysis. Global observed precipitation is used as direct forcing to drive the land surface analysis, rather than the typical reanalysis approach of using precipitation assimilating from a background atmospheric model simulation. Global observed snow cover and snow depth fields are used to constrain the simulated snow variables. This paper describes 1) the design and implementation of GLDAS/LIS in CFSR. 2) the forcing of the observed global precipitation and snow fields, and 3) preliminary results of global and regional soil moisture content and land surface energy and water budgets closure. With special attention made during the design of CFSR GLDAS/LIS, all the source and sink terms in the CFSR land surface energy and water budgets can be assessed and the total budgets are balanced. This is one of many aspects indicating improvements in CFSR from the previous NCEP reanalyses.

The Noah model is a land surface model of the National Centers for Environmental Prediction. It has been widely used in regional coupled weather and climate models (i.e. Weather Research and Forecasting Model, Eta Mesoscale Model) and global coupled weather and climate models (i.e. National Centers for Environmental Prediction Global Forecast System, Climate Forecast System). Therefore, its continued improvement and development are keys to enhancing our weather and climate forecast ability and water and energy flux simulation accuracy. North American Land Data Assimilation System phase 1 (NLDAS‐1) experiments indicated that the Noah model exhibited substantial bias in latent heat flux, total runoff and land skin temperature during the warm season, and such bias can significantly affect coupled weather and climate models. This paper presents a study to improve the Noah model by adding model parameterization processes such as including seasonal factor on leaf area index and root distribution and selecting optimal model parameters. We compared simulated latent heat flux, mean annual runoff and land skin temperature from the Noah control and test versions with measured latent heat flux, land surface skin temperature, mean annual runoff and satellite‐retrieved land surface skin temperature. The results show that the test version significantly reduces biases in latent heat, total runoff and land skin temperature simulation. The test version has been used for the NLDAS phase 2 (NLDAS‐2) to produce 30‐year water flux, energy flux and state variable products to support the US drought monitor of National Integrated Drought Information System. Copyright © 2012 John Wiley & Sons, Ltd.

Results are presented from the second phase of the multiinstitution North American Land Data Assimilation System (NLDAS‐2) research partnership. In NLDAS, the Noah, Variable Infiltration Capacity, Sacramento Soil Moisture Accounting, and Mosaic land surface models (LSMs) are executed over the conterminous U.S. (CONUS) in realtime and retrospective modes. These runs support the drought analysis, monitoring and forecasting activities of the National Integrated Drought Information System, as well as efforts to monitor large‐scale floods. NLDAS‐2 builds upon the framework of the first phase of NLDAS (NLDAS‐1) by increasing the accuracy and consistency of the surface forcing data, upgrading the land surface model code and parameters, and extending the study from a 3‐year (1997–1999) to a 30‐year (1979–2008) time window. As the first of two parts, this paper details the configuration of NLDAS‐2, describes the upgrades to the forcing, parameters, and code of the four LSMs, and explores overall model‐to‐model comparisons of land surface water and energy flux and state variables over the CONUS. Focusing on model output rather than on observations, this study seeks to highlight the similarities and differences between models, and to assess changes in output from that seen in NLDAS‐1. The second part of the two‐part article focuses on the validation of model‐simulated streamflow and evaporation against observations. The results depict a higher level of agreement among the four models over much of the CONUS than was found in the first phase of NLDAS. This is due, in part, to recent improvements in the parameters, code, and forcing of the NLDAS‐2 LSMs that were initiated following NLDAS‐1. However, large inter‐model differences still exist in the northeast, Lake Superior, and western mountainous regions of the CONUS, which are associated with cold season processes. In addition, variations in the representation of sub‐surface hydrology in the four LSMs lead to large differences in modeled evaporation and subsurface runoff. These issues are important targets for future research by the land surface modeling community. Finally, improvement from NLDAS‐1 to NLDAS‐2 is summarized by comparing the streamflow measured from U.S. Geological Survey stream gauges with that simulated by four NLDAS models over 961 small basins.

The augmented Noah land surface model described in the first part of the two-part series was evaluated here over global river basins. Across various climate zones, global-scale tests can reveal a model's weaknesses and strengths that a local-scale testing cannot. In addition, global-scale tests are more challenging than local-and catchment-scale tests. Given constant model parameters (e. g., runoff parameters) across global river basins, global-scale tests are more stringent. We assessed model performance against various satellite and ground-based observations over global river basins through six experiments that mimic a transition from the original Noah LSM to the fully augmented version. The model shows transitional improvements in modeling runoff, soil moisture, snow, and skin temperature, despite considerable increase in computational time by the fully augmented Noah-MP version compared to the original Noah LSM. The dynamic vegetation model favorably captures seasonal and spatial variability of leaf area index and green vegetation fraction. We also conducted 36 ensemble experiments with 36 combinations of optional schemes for runoff, leaf dynamics, stomatal resistance, and the beta factor. Runoff schemes play a dominant and different role in controlling soil moisture and its relationship with evapotranspiration compared to ecological processes such as the beta factor, vegetation dynamics, and stomatal resistance. The 36-member ensemble mean of runoff performs better than any single member over the world's 50 largest river basins, suggesting a great potential of land-based ensemble simulations for climate prediction.

This first paper of the two-part series describes the objectives of the community efforts in improving the Noah land surface model (LSM), documents, through mathematical formulations, the augmented conceptual realism in biophysical and hydrological processes, and introduces a framework for multiple options to parameterize selected processes (Noah-MP). The Noah-MP's performance is evaluated at various local sites using high temporal frequency data sets, and results show the advantages of using multiple optional schemes to interpret the differences in modeling simulations. The second paper focuses on ensemble evaluations with long-term regional (basin) and global scale data sets. The enhanced conceptual realism includes (1) the vegetation canopy energy balance, (2) the layered snowpack, (3) frozen soil and infiltration, (4) soil moisture-groundwater interaction and related runoff production, and (5) vegetation phenology. Sample local-scale validations are conducted over the First International Satellite Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE) site, the W3 catchment of Sleepers River, Vermont, and a French snow observation site. Noah-MP shows apparent improvements in reproducing surface fluxes, skin temperature over dry periods, snow water equivalent (SWE), snow depth, and runoff over Noah LSM version 3.0. Noah-MP improves the SWE simulations due to more accurate simulations of the diurnal variations of the snow skin temperature, which is critical for computing available energy for melting. Noah-MP also improves the simulation of runoff peaks and timing by introducing a more permeable frozen soil and more accurate simulation of snowmelt. We also demonstrate that Noah-MP is an effective research tool by which modeling results for a given process can be interpreted through multiple optional parameterization schemes in the same model framework.

Accurately representing complex land-surface processes balancing complexity and realism remains one challenge that the weather modelling community is facing nowadays. In this study, a photosynthesis-based Gas-exchange Evapotranspiration Model (GEM) is integrated into the Noah land-surface model replacing the traditional Jarvis scheme for estimating the canopy resistance and transpiration. Using 18-month simulations from the High Resolution Land Data Assimilation System (HRLDAS), the impact of the photosynthesis-based approach on the simulated canopy resistance, surface heat fluxes, soil moisture, and soil temperature over different vegetation types is evaluated using data from the Atmospheric Radiation Measurement (ARM) site, Oklahoma Mesonet, 2002 International H2O Project (IHOP_2002), and three Ameriflux sites. Incorporation of GEM into Noah improves the surface energy fluxes as well as the associated diurnal cycle of soil moisture and soil temperature during both wet and dry periods. An analysis of midday, average canopy resistance shows similar day-to-day trends in the model fields as seen in observed patterns. Bias and standard deviation analyses for soil temperature and surface fluxes show that GEM responds somewhat better than the Jarvis scheme, mainly because the Jarvis approach relies on a parametrised minimum canopy resistance and meteorological variables such as air temperature and incident radiation. The analyses suggest that adding a photosynthesis-based transpiration scheme such as GEM improves the ability of the land-data assimilation system to simulate evaporation and transpiration under a range of soil and vegetation conditions.

To examine the impact from land model upgrades and different land initializations on the National Centers for Environmental Prediction (NCEP)'s Climate Forecast System (CFS), extensive T126 CFS experiments are carried out for 25 summers with 10 ensemble members using the old Oregon State University (OSU) land surface model (LSM) and the new Noah LSM. The CFS using the Noah LSM, initialized in turn with land states from the NCEP-Department of Energy Global Reanalysis 2 (GR-2), Global Land Data System (GLDAS), and GLDAS climatology, is compared to the CFS control run using the OSU LSM initialized with the GR-2 land states. Using anomaly correlation as a primary measure, the summer-season prediction skill of the CFS using different land models and different initial land states is assessed for SST, precipitation, and 2-m air temperature over the contiguous United States (CONUS) on an ensemble basis.Results from these CFS experiments indicate that upgrading from the OSU LSM to the Noah LSM improves the overall CONUS June-August (JJA) precipitation prediction, especially during ENSO neutral years. Such an enhancement in CFS performance requires the execution of a GLDAS with the very same Noah LSM as utilized in the land component of the CFS, while improper initializations of the Noah LSM using the GR-2 land states lead to degraded CFS performance. In comparison with precipitation, the land upgrades have a relatively small impact on both of the SST and 2-m air temperature predictions.

To provide quality‐improved and consistent real‐time global green vegetation fraction (GVF) data products that are suitable for use in operational numerical weather, climate, and hydrological models, necessary processing steps are applied to the output data stream from the advanced very high resolution radiometer (AVHRR)‐based NOAA operational global vegetation index (GVI) system. This paper reviewed the NOAA GVI data and described the algorithm to derive weekly updated real‐time GVF from the normalized difference vegetation index (NDVI). The methodology description focuses on algorithm justification in an operational production context. The described algorithm was implemented in the global vegetation processing system (GVPS). The new global GVF data sets include the multiyear GVF weekly climatology and the real‐time weekly GVF. Compared to the old 5 year GVF monthly climatology currently used in the operational National Centers for Environmental Prediction (NCEP)/Environmental Modeling Center (EMC) weather and climate models, the new data sets provide an overall higher vegetation value, real‐time surface vegetation information, and numerous other improvements. The new GVF data set quality was partially assured by validation against Moderate Resolution Imaging Spectroradiometer (MODIS) NDVI at a few EOS land validation core sites and comparison with another independently processed NDVI data set. Impact of the new GVF data sets in numerical weather prediction (NWP) model was investigated using EMC mesoscale model simulations and concluded overall positive.

A negative snow water equivalent (SWE) bias in the snow model of the Noah land surface scheme used in the NCEP suite of numerical weather and climate prediction models has been noted by several investigators. This bias motivated a series of offline tests of model extensions and improvements intended to reduce or eliminate the bias. These improvements consist of changes to the model's albedo formulation that include a parameterization for snowpack aging, changes to how pack temperature is computed, and inclusion of a provision for refreeze of liquid water in the pack. Less extensive testing was done on the performance of model extensions with alternate areal depletion parameterizations. Model improvements were evaluated through comparisons of point simulations with National Resources Conservation Service (NRCS) Snowpack Telemetry (SNOTEL) SWE data for deep-mountain snowpacks at selected stations in the western United States, as well as simulations of snow areal extent over the conterminous United States (CONUS) domain, compared with observational data from the NOAA Interactive Multisensor Snow and Ice Mapping System (IMS). The combination of snow-albedo decay and liquid-water refreeze results in substantial improvements in the magnitude and timing of peak SWE, as well as increased snow-covered extent at large scales. Modifications to areal snow depletion thresholds yielded more realistic snow-covered albedos at large scales.

Abstract Hydrological processes are strongly coupled with atmospheric processes related, for example, to precipitation and temperature, and a coupled atmosphere–land surface system is required for a meaningful hydrological forecast. Since the atmosphere is a chaotic system with limited predictability, ensemble forecasts offer a practical tool to predict the future state of the coupled system in a probabilistic fashion, potentially leading to a more complete and informative hydrologic prediction. As ensemble forecasts with coupled meteorological–hydrological models are operationally running at major numerical weather prediction centers, it is currently possible to produce a gridded streamflow prognosis in the form of a probabilistic forecast based on ensembles. Evaluation and improvement of such products require a comprehensive assessment of both components of the coupled system. In this article, the atmospheric component of a coupled ensemble forecasting system is evaluated in terms of its ability to prov...