[1] 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.

Scientists from several institutions and with different research backgrounds have worked together to develop a prototype modular land model for weather forecasting and climate studies. This model is now available for public use and further development.C limate and weather forecasting models require the energy, water, and momentum fluxes across the land-atmosphere interface to be specified. Various land surface parameterizations (LSPs), ranging from the simple bucket-type LSP in the 1960s to the current soil-vegetation-atmosphere interactive LSP, have been developed in the past four decades to calculate these fluxes. The Project for Intercomparison of Land Surface Parameterization Schemes (PILPS) has demonstrated that, even with the same atmospheric forcing data and similar land surface parameters, different LSPs still give significantly different surface fluxes and soil wetness, partly because of the differences in the formulations of individual processes and coding architectures in participant models . On the other hand, most LSPs share many common components, suggesting the need to develop a publicly available common land model with a modular structure that could facilitate the exploration of new issues, less repetition of past efforts, and sharing of improvements and refinements contributed by different groups.The Common Land Model (CLM) effort dates back to the mid-1990s and has evolved through various workshops and e-mail correspondence. The initial motivation was to provide a framework for a truly community-developed land component of the National Center for Atmospheric Research (NCAR) Community Climate System Model (CCSM). Interest in applying CLM came from the Goddard Space Flight Center (GSFC) Data Assimilation Office (DAO), which was implementing the Mosaic model (Koster and Suarez 1992), and the Center for Ocean-Land-Atmosphere Studies (COLA) scientists, who were revising their Simplified Simple Biosphere Model (SSiB) (Xue et al. 1991). We also established ties to groups performing carbon cycle and ecological modeling.In developing CLM, we attempted to combine the best features of three existing successful and relatively

[1] The Community Land Model version 3 (CLM3) is the land component of the Community Climate System Model (CCSM). CLM3 has energy and water biases resulting from deficiencies in some of its canopy and soil parameterizations related to hydrological processes. Recent research by the community that utilizes CLM3 and the family of CCSM models has indicated several promising approaches to alleviating these biases. This paper describes the implementation of a selected set of these parameterizations and their effects on the simulated hydrological cycle. The modifications consist of surface data sets based on Moderate Resolution Imaging Spectroradiometer products, new parameterizations for canopy integration, canopy interception, frozen soil, soil water availability, and soil evaporation, a TOPMODEL-based model for surface and subsurface runoff, a groundwater model for determining water table depth, and the introduction of a factor to simulate nitrogen limitation on plant productivity. The results from a set of offline simulations were compared with observed data for runoff, river discharge, soil moisture, and total water storage to assess the performance of the new model (referred to as CLM3.5). CLM3.5 exhibits significant improvements in its partitioning of global evapotranspiration (ET) which result in wetter soils, less plant water stress, increased transpiration and photosynthesis, and an improved annual cycle of total water storage. Phase and amplitude of the runoff annual cycle is generally improved. Dramatic improvements in vegetation biogeography result when CLM3.5 is coupled to a dynamic global vegetation model. Lower than observed soil moisture variability in the rooting zone is noted as a remaining deficiency.

[1] Groundwater interacts with soil moisture through the exchanges of water between the unsaturated soil and its underlying aquifer under gravity and capillary forces. Despite its importance, groundwater is not explicitly represented in climate models. This paper developed a simple groundwater model (SIMGM) by representing recharge and discharge processes of the water storage in an unconfined aquifer, which is added as a single integration element below the soil of a land surface model. We evaluated the model against the Gravity Recovery and Climate Experiment (GRACE) terrestrial water storage change (DS) data. The modeled total water storage (including unsaturated soil water and groundwater) change agrees fairly well with GRACE estimates. The anomaly of the modeled groundwater storage explains most of the GRACE DS anomaly in most river basins where the water storage is not affected by snow water or frozen soil. For this reason, the anomaly of the modeled water table depth agrees well with that converted from the GRACE DS in most of the river basins. We also investigated the impacts of groundwater dynamics on soil moisture and evapotranspiration through the comparison of SIMGM to an additional model run using gravitational free drainage (FD) as the model's lower boundary condition. SIMGM produced much wetter soil profiles globally and up to 16% more annual evapotranspiration than FD, most obviously in arid-to-wet transition regions.

[1] This paper develops a simple TOPMODEL-based runoff parameterization (hereinafter SIMTOP) for use in global climate models (GCMs) that improves the runoff production and the partitioning of runoff between surface and subsurface components. SIMTOP simplifies the TOPMODEL runoff formulations in two ways: (1) SIMTOP represents the discrete distribution of the topographic index as an exponential function, not as a three-parameter gamma distribution; this change improves the parameterization of the fractional saturated area, especially in mountainous regions. (2) SIMTOP treats subsurface runoff as a product of an exponential function of the water table depth and a single coefficient, not as a product of several parameters that depend on topography and soil properties; this change facilitates applying TOPMODEL-based runoff schemes on global scale. SIMTOP is incorporated into the National Center for Atmospheric Research (NCAR) Community Land Model version 2.0 (CLM 2.0). SIMTOP is validated at a watershed scale using data from the Sleepers River watershed in Vermont, USA. It is also validated on a global scale using the monthly runoff data from the University of New Hampshire Global Runoff Data Center (UNH-GRDC). SIMTOP performs favorably when compared to the baseline runoff formulation used in CLM2.0. Realistic simulations can be obtained using two distinct saturated hydraulic conductivity (K sat ) profiles. These profiles include (1) exponential decay of K sat with depth (as is typically done in TOPMODEL-based runoff schemes) and (2) the definition of K sat using the soil texture profile data (as is typically done in climate models) and the concordant reduction of the gravitational drainage from the bottom of the soil column.

The presence of ice in soil dramatically alters soil hydrologic and thermal properties. Despite this important role, many recent studies show that explicitly including the hydrologic effects of soil ice in land surface models degrades the simulation of runoff in cold regions. This paper addresses this dilemma by employing the Community Land Model version 2.0 (CLM2.0) developed at the National Center for Atmospheric Research (NCAR) and a simple TOPMODEL-based runoff scheme (SIMTOP). CLM2.0/ SIMTOP explicitly computes soil ice content and its modifications to soil hydrologic and thermal properties. However, the frozen soil scheme has a tendency to produce a completely frozen soil (100% ice content) whenever the soil temperature is below 0°C. The frozen ground prevents infiltration of snowmelt or rainfall, thereby resulting in earlier-and higher-than-observed springtime runoff. This paper presents modifications to the above-mentioned frozen soil scheme that produce more accurate magnitude and seasonality of runoff and soil water storage. These modifications include 1) allowing liquid water to coexist with ice in the soil over a wide range of temperatures below 0°C by using the freezing-point depression equation, 2) computing the vertical water fluxes by introducing the concept of a fractional permeable area, which partitions the model grid into an impermeable part (no vertical water flow) and a permeable part, and 3) using the total soil moisture (liquid water and ice) to calculate the soil matric potential and hydraulic conductivity. The performance of CLM2.0/SIMTOP with these changes has been tested using observed data in cold-region river basins of various spatial scales. Compared to the CLM2.0/SIMTOP frozen soil scheme, the modified scheme produces monthly runoff that compares more favorably with that estimated by the University of New Hampshire-Global Runoff Data Center and a terrestrial water storage change that is in closer agreement with that measured by the Gravity Recovery and Climate Experiment (GRACE) satellites.

Earth's terrestrial near-subsurface environment can be divided into relatively porous layers of soil, intact regolith, and sedimentary deposits above unweathered bedrock. Variations in the thicknesses of these layers control the hydrologic and biogeochemical responses of landscapes. Currently, Earth System Models approximate the thickness of these relatively permeable layers above bedrock as uniform globally, despite the fact that their thicknesses vary systematically with topography, climate, and geology. To meet the need for more realistic input data for models, we developed a high-resolution gridded global data set of the average thicknesses of soil, intact regolith, and sedimentary deposits within each 30 arcsec (1 km) pixel using the best available data for topography, climate, and geology as input. Our data set partitions the global land surface into upland hillslope, upland valley bottom, and lowland landscape components and uses models optimized for each landform type to estimate the thicknesses of each subsurface layer. On hillslopes, the data set is calibrated and validated using independent data sets of measured soil thicknesses from the U.S. and Europe and on lowlands using depth to bedrock observations from groundwater wells in the U.S. We anticipate that the data set will prove useful as an input to regional and global hydrological and ecosystems models.

The mapped rivers and streams of the contiguous United States are available in a geographic information system (GIS) dataset called National Hydrography Dataset Plus (NHDPlus). This hydrographic dataset has about 3 million river and water body reaches along with information on how they are connected into networks. The U.S. Geological Survey (USGS) National Water Information System (NWIS) provides streamflow observations at about 20 thousand gauges located on the NHDPlus river network. A river network model called Routing Application for Parallel Computation of Discharge (RAPID) is developed for the NHDPlus river network whose lateral inflow to the river network is calculated by a land surface model. A matrix-based version of the Muskingum method is developed herein, which RAPID uses to calculate flow and volume of water in all reaches of a river network with many thousands of reaches, including at ungauged locations. Gauges situated across river basins (not only at basin outlets) are used to automatically optimize the Muskingum parameters and to assess river flow computations, hence allowing the diagnosis of runoff computations provided by land surface models. RAPID is applied to the Guadalupe and San Antonio River basins in Texas, where flow wave celerities are estimated at multiple locations using 15-min data and can be reproduced reasonably with RAPID. This river model can be adapted for parallel computing and although the matrix method initially adds a large overhead, river flow results can be obtained faster than with the traditional Muskingum method when using a few processing cores, as demonstrated in a synthetic study using the upper Mississippi River basin.