Search

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 (dS/dt, 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 dS/dt 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.

This paper compares the annual and monthly components of the simulated energy budget from the North American Land Data Assimilation System phase 2 (NLDAS-2) with reference products over the domains of the 12 River Forecast Centers (RFCs) of the continental United States (CONUS). The simulations are calculated from both operational and research versions of NLDAS-2. The reference radiation components are obtained from the National Aeronautics and Space Administration Surface Radiation Budget product. The reference sensible and latent heat fluxes are obtained from a multitree ensemble method applied to gridded FLUXNET data from the Max Planck Institute, Germany. As these references are obtained from different data sources, they cannot fully close the energy budget, although the range of closure error is less than 15%formean annual results. The analysis here demonstrates the usefulness of basin-scale surface energy budget analysis for evaluating model skill and deficiencies. The operational (i.e., Noah, Mosaic, and VIC) and research (i.e., Noah-I and VIC4.0.5) NLDAS-2 land surface models exhibit similarities and differences in depicting basin-averaged energy components. For example, the energy components of the five models have similar seasonal cycles, but with different magnitudes. Generally, Noah and VIC overestimate (underestimate) sensible (latent) heat flux over several RFCs of the eastern CONUS. In contrast, Mosaic underestimates (overestimates) sensible (latent) heat flux over almost all 12 RFCs. The research Noah-I and VIC4.0.5 versions show moderate-to-large improvements (basin and model dependent) relative to their operational versions, which indicates likely pathways for future improvements in the operational NLDAS-2 system.

In the middle to high latitude and alpine regions, the seasonal snow pack can dominate the surface energy and water budgets due to its high albedo, low thermal conductivity, high emissivity, considerable spatial and temporal variability, and ability to store and then later release a winters cumulative snowfall (Cohen, 1994; Hall, 1998). With this in mind, the snow drought across the U.S. has raised questions about impacts on water supply, ski resorts and agriculture. Knowledge of various snow pack properties is crucial for short-term weather forecasts, climate change prediction, and hydrologic forecasting for producing reliable daily to seasonal forecasts. One potential source of this information is the multi-institution North American Land Data Assimilation System (NLDAS) project (Mitchell et al., 2004). Real-time NLDAS products are used for drought monitoring to support the National Integrated Drought Information System (NIDIS) and as initial conditions for a future NCEP drought forecast system. Additionally, efforts are currently underway to assimilate remotely-sensed estimates of land-surface states such as snowpack information into NLDAS. It is believed that this assimilation will not only produce improved snowpack states that better represent snow evolving conditions, but will directly improve the monitoring of drought.

Snow cover area affects snowmelt, soil moisture, evapotranspiration, and ultimately streamflow. For the Distributed Model Intercomparison Project - Phase 2 Western basins, we assimilate satellite-based fractional snow cover area (fSCA) from the Moderate Resolution Imaging Spectroradiometer, or MODIS, into the National Weather Service (NWS) SNOW-17 model. This model is coupled with the NWS Sacramento Heat Transfer (SAC-HT) model inside the National Aeronautics and Space Administration's (NASA) Land Information System. SNOW-17 computes fSCA from snow water equivalent (SWE) values using an areal depletion curve. Using a direct insertion, we assimilate fSCAs in two fully distributed ways: 1) we update the curve by attempting SWE preservation, and 2) we reconstruct SWEs using the curve. The preceding are refinements of an existing simple, conceptually-guided NWS algorithm. Satellite fSCA over dense forests inadequately accounts for below-canopy snow, degrading simulated streamflow upon assimilation during snowmelt. Accordingly, we implement a below-canopy allowance during assimilation. This simplistic allowance and direct insertion are found to be inadequate for improving calibrated results, still degrading them as mentioned above. However, for streamflow volume for the uncalibrated runs, we obtain: (1) substantial to major improvements (64-81 %) as a percentage of the control run residuals (or distance from observations), and (2) minor improvements (16-22 %) as a percentage of observed values. We highlight the need for detailed representations of canopy-snow optical radiative transfer processes in mountainous, dense forest regions if assimilation-based improvements are to be seen in calibrated runs over these areas.

Model simulated soil moisture fields are often biased due to errors in input parameters and deficiencies in model physics. Satellite derived soil moisture estimates, if retrieved appropriately, represent the spatial mean of soil moisture in a footprint area, and can be used to reduce model bias (at locations near the surface) through data assimilation techniques. While assimilating the retrievals can reduce model bias, it can also destroy the mass balance enforced by the model governing equation because water is removed from or added to the soil by the assimilation algorithm. In addition, studies have shown that assimilation of surface observations can adversely impact soil moisture estimates in the lower soil layers due to imperfect model physics, even though the bias near the surface is decreased. In this study, an ensemble Kalman filter (EnKF) with a mass conservation updating scheme was developed to assimilate the actual value of Advanced Microwave Scanning Radiometer (AMSR-E) soil moisture retrievals to improve the mean of simulated soil moisture fields by the Noah land surface model. Assimilation results using the conventional and the mass conservation updating scheme in the Little Washita watershed of Oklahoma showed that, while both updating schemes reduced the bias in the shallow root zone, the mass conservation scheme provided better estimates in the deeper profile. The mass conservation scheme also yielded physically consistent estimates of fluxes and maintained the water budget. Impacts of model physics on the assimilation results are discussed.

The accurate portrayal of the hydrological cycle is extremely important in land surface modeling. Central to this effort is the treatment of snow, as errors in the representation of this quantity can impact practically all other modeled quantities through alterations in the water and energy balances. Although land surface model (LSM) simulations can benefit from the assimilation of snow cover and snow depth observations, they can be negatively impacted if such observations contain errors or if a model bias exists in the simulation of surface or soil temperatures. Both cases may lead to excessive melting or growth of snow packs, and to large alterations in both the energy and water balances. Such problems in the snow assimilation process, made evident by the repeated melting and replenishing of snow pack over significant areas of the United States, exists in the Eta Data Assimilation System and is a product of the EDAS system's direct insertion assimilation of snow data. Occurring on a 24 hour cycle, the repeated melting infuses the soil column with a large quantity of water that upsets the hydrological cycle. In an effort to quantify the impacts of such errors in snow assimilation on water and energy budgets, a series of Mosaic LSM simulations were performed over the 12 month period covering October 1998 to October 1999.