Search

Permanently frozen and seasonally frozen soils occur over a large portion of the Earth's land surface. Changes in the freeze/thaw state of the land surface reflects major changes in thermal and hydraulic properties as well as acting as a "switch" for many ecological processes. In short, soil freeze/thaw state is a fundamental land surface variable in the water and energy cycles, and it connects to the carbon cycle. Surface freeze/thaw state is observable by passive and active microwave sensors. For example, NASA's Soil Moisture Active Passive (SMAP) mission includes a freeze/thaw data product. Such satellite sensing offers routine all-season and all-weather global observations of soil freeze/thaw state with the application of suitable algorithms. We describe early finding from the SLAPex Freeze/Thaw campaign, believed to be the first airborne campaign of its type, focusing on soil freeze/thaw.

In flat homogenous terrain such as in Cambodia and Vietnam, the monsoon season brings significant and consistent flooding between May and November. To monitor flooding in the Lower Mekong region, the near real-time NASA Flood Extent Product (NASA-FEP) was developed using seasonal normalized difference vegetation index (NDVI) differences from the 250 m resolution Moderate Resolution Imaging Spectroradiometer (MODIS) sensor compared to daily observations. The use of a percentage change interval classification relating to various stages of flooding reduces might be confusing to viewers or potential users, and therefore reducing the product usage. To increase the product usability through simplification, the classification intervals were compared with other commonly used change detection schemes to identify the change classification scheme that best delineates flooded areas. The percentage change method used in the NASA-FEP proved to be helpful in delineating flood boundaries compared to other change detection methods. The results of the accuracy assessments indicate that the −75% NDVI change interval can be reclassified to a descriptive 'flood' classification. A binary system was used to simplify the interpretation of the NASA-FEP by removing extraneous information from lower interval change classes.

Soil freezing and thawing is an important process in the terrestrial water, energy, and carbon cycles, marking the change between two very different hydraulic, thermal, and biological regimes. NASA's Soil Moisture Active/Passive (SMAP) mission includes a binary freeze/thaw data product. While there have been ground-based remote sensing field measurements observing soil freeze/thaw at the point scale, and airborne campaigns that observed some frozen soil areas (e.g., BOREAS), the recently-completed SLAPex Freeze/Thaw (F/T) campaign is the first airborne campaign dedicated solely to observing frozen/thawed soil with both passive and active microwave sensors and dedicated ground truth, in order to enable detailed process-level exploration of the remote sensing signatures and in situ soil conditions. SLAPex F/T utilized the Scanning L-band Active/Passive (SLAP) instrument, an airborne simulator of SMAP developed at NASA's Goddard Space Flight Center, and was conducted near Winnipeg, Manitoba, Canada, in October/November, 2015. Future soil moisture missions are also expected to include soil freeze/thaw products, and the loss of the radar on SMAP means that airborne radar-radiometer observations like those that SLAP provides are unique assets for freeze/thaw algorithm development. This paper will present an overview of SLAPex F/T, including descriptions of the site, airborne and ground-based remote sensing, ground truth, as well as preliminary results.

Precipitation in the central U.S. decreases by about 25% during the seasonal transition from June to July, and this precipitation decrease has been observed to have intensified since 1979. Such an intensification could enhance future spring drought occurrences such as was the case in the 2012 "flash drought" in the Midwestern U.S., where conditions evolved quickly from being abnormally dry to exceptionally dry within a mere month from June to July. In this study, various atmospheric and land reanalysis datasets were analyzed to examine the trend calculated from 1979 to 2012 in the June-to-July seasonal transition. It was found that the change in precipitation deficit was accompanied by increased downward shortwave radiation flux and tropospheric subsidence, enhanced evaporative fraction, as well as an elevated planetary boundary layer height. The change in the tropospheric circulation encompassed an anomalous ridge over the western U.S. and a trough on either side; this wave-form circulation pattern is known to induce dry conditions in the central U.S. Possibly, the trends in the June-to-July seasonal shifts in precipitation, drought severity and tropospheric circulation intensified the 2012 "flash drought" in timing and extent. The knowledge of the trends allows one to anticipate the evolution of spring onset of drought into the summer.

This paper summarizes research related to the 2012 record drought in the central United States conducted by members of the NASA Energy and Water cycle Study (NEWS) Working Group. Past drought patterns were analyzed for signal coherency with latest drought and the contribution of long-term trends in the Great Plains low-level jet, an important regional circulation feature of the spring rainy season in the Great Palins. Long-term changes in the seasonal transition from rainy spring into dry summer were also examined. Potential external forcing from radiative processes, soil-air interactions, and ocean teleconnections were assessed as contributors to the intensity of the drought. The atmospheric Rossby wave activity was found to be a potential source of predictability for the onset of drought. A probabilistic model was introduced and evaluated for its performance in predicting drought recovery in the Great Plains.

Monitoring Earth's terrestrial water conditions is critically important to many hydrological applications such as global food production; assessing water resources sustainability; and flood, drought, and climate change prediction. These needs have motivated the development of pilot monitoring and prediction systems for terrestrial hydrologic and vegetative states, but to date only at the rather coarse spatial resolutions (approx.10-100 km) over continental to global domains. Adequately addressing critical water cycle science questions and applications requires systems that are implemented globally at much higher resolutions, on the order of 1 km, resolutions referred to as hyperresolution in the context of global land surface models. This opinion paper sets forth the needs and benefits for a system that would monitor and predict the Earth's terrestrial water, energy, and biogeochemical cycles. We discuss six major challenges in developing a system: improved representation of surface-subsurface interactions due to fine-scale topography and vegetation; improved representation of land-atmospheric interactions and resulting spatial information on soil moisture and evapotranspiration; inclusion of water quality as part of the biogeochemical cycle; representation of human impacts from water management; utilizing massively parallel computer systems and recent computational advances in solving hyperresolution models that will have up to 10(exp 9) unknowns; and developing the required in situ and remote sensing global data sets. We deem the development of a global hyperresolution model for monitoring the terrestrial water, energy, and biogeochemical cycles a grand challenge to the community, and we call upon the international hydrologic community and the hydrological science support infrastructure to endorse the effort.

Hydrological years 2006 (HY06, 10/2005-09/2006) and 2007 (HY07, 10/2006-09/2007) provide a unique opportunity to examine hydrological extremes in the central US because there are no other examples of two such highly contrasting precipitation extremes occurring in consecutive years at the Southern Great Plains (SGP) in recorded history. The HY06 annual precipitation in the state of Oklahoma, as observed by the Oklahoma Mesonet, is around 61% of the normal (92.84 cm, based on the 1921-2008 climatology), which results in HY06 the second-driest year in the record. In particular, the total precipitation during the winter of 2005-06 is only 27% of the normal, and this winter ranks as the driest season. On the other hand, the HY07 annual precipitation amount is 121% of the normal and HY07 ranks as the seventh-wettest year for the entire state and the wettest year for the central region of the state. Summer 2007 is the second-wettest season for the state. Large-scale dynamics play a key role in these extreme events. During the extreme dry period (10/2005-02/2006), a dipole pattern in the 500-hPa GH anomaly existed where an anomalous high was over the southwestern U.S. region and an anomalous low was over the Great Lakes. This pattern is associated with inhibited moisture transport from the Gulf of Mexico and strong sinking motion over the SGP, both contributing to the extreme dryness. The precipitation deficit over the SGP during the extreme dry period is clearly linked to significantly suppressed cyclonic activity over the southwestern U.S., which shows robust relationship with the Western Pacific (WP) teleconnection pattern. The precipitation events during the extreme wet period (May-July 2007) were initially generated by active synoptic weather patterns, linked with moisture transport from the Gulf of Mexico by the northward low level jet, and enhanced by the mesoscale convective systems. Although the drought and pluvial conditions are dominated by large-scale dynamic patterns, we have demonstrated that the two positive feedback processes during the extreme dry and wet periods found in this study play a key role to maintain and reinforce the length and severity of existing drought and flood events. For example, during the extreme dry period, with less clouds, LWP, PWV, precipitation, and thinner Cu cloud thickness, more net radiation was absorbed and used to evaporate water from the ground. The evaporated moisture, however, was removed by low-level divergence. Thus, with less precipitation and removed atmospheric moisture, more absorbed incoming solar radiation was used to increase surface temperature and to make the ground drier.

An ensemble Kalman filter (EnKF) is used in a suite of synthetic experiments to assimilate coarse-scale (25 km) snow water equivalent (SWE) observations (typical of satellite retrievals) into fine-scale (1 km) model simulations. Coarse-scale observations are assimilated directly using an observation operator for mapping between the coarse and fine scales or, alternatively, after disaggregation (re-gridding) to the fine-scale model resolution prior to data assimilation. In either case observations are assimilated either simultaneously or independently for each location. Results indicate that assimilating disaggregated fine-scale observations independently (method 1D-F1) is less efficient than assimilating a collection of neighboring disaggregated observations (method 3D-Fm). Direct assimilation of coarse-scale observations is superior to a priori disaggregation. Independent assimilation of individual coarse-scale observations (method 3D-C1) can bring the overall mean analyzed field close to the truth, but does not necessarily improve estimates of the fine-scale structure. There is a clear benefit to simultaneously assimilating multiple coarse-scale observations (method 3D-Cm) even as the entire domain is observed, indicating that underlying spatial error correlations can be exploited to improve SWE estimates. Method 3D-Cm avoids artificial transitions at the coarse observation pixel boundaries and can reduce the RMSE by 60% when compared to the open loop in this study.

The definition and derivation of a 5-year, 0.125deg, 3-hourly atmospheric forcing dataset for the South America continent is described which is appropriate for use in a Land Data Assimilation System and which, because of the limited surface observational networks available in this region, uses remotely sensed data merged with surface observations as the basis for the precipitation and downward shortwave radiation fields. The quality of this data set is evaluated against available surface observations. There are regional difference in the biases for all variables in the dataset, with biases in precipitation of the order 0-1 mm/day and RMSE of 5-15 mm/day, biases in surface solar radiation of the order 10 W/sq m and RMSE of 20 W/sq m, positive biases in temperature typically between 0 and 4 K, depending on region, and positive biases in specific humidity around 2-3 g/Kg in tropical regions and negative biases around 1-2 g/Kg further south.

Land surface models (LSMs) are computer programs, similar to weather and climate prediction models, which simulate the stocks and fluxes of water (including soil moisture, snow, evaporation, and runoff) and energy (including the temperature of and sensible heat released from the soil) after they arrive on the land surface as precipitation and sunlight. It is not currently possible to measure all of the variables of interest everywhere on Earth with sufficient accuracy and space-time resolution. Hence LSMs have been developed to integrate the available observations with our understanding of the physical processes involved, using powerful computers, in order to map these stocks and fluxes as they change in time. The maps are used to improve weather forecasts, support water resources and agricultural applications, and study the Earth's water cycle and climate variability. NASA's Global Land Data Assimilation System (GLDAS) project facilitates testing of several different LSMs with a variety of input datasets (e.g., precipitation, plant type). Precipitation is arguably the most important input to LSMs. Many precipitation datasets have been produced using satellite and rain gauge observations and weather forecast models. In this study, seven different global precipitation datasets were evaluated over the United States, where dense rain gauge networks contribute to reliable precipitation maps. We then used the seven datasets as inputs to GLDAS simulations, so that we could diagnose their impacts on output stocks and fluxes of water. In terms of totals, the Climate Prediction Center (CPC) Merged Analysis of Precipitation (CMAP) had the closest agreement with the US rain gauge dataset for all seasons except winter. The CMAP precipitation was also the most closely correlated in time with the rain gauge data during spring, fall, and winter, while the satellitebased estimates performed best in summer. The GLDAS simulations revealed that modeled soil moisture is highly sensitive to precipitation, with differences in spring and summer as large as 45% depending on the choice of precipitation input.

The U.S. Bureau of Reclamation identified the remote sensing of evapotranspiration (ET) as an important water flux for study and designated a test site in the Lower Colorado River basin. A consortium of groups will work together with the goal to develop more accurate and cost effective techniques using the enhanced spatial and temporal coverage afforded by remote sensing. ET is a critical water loss flux where improved estimation should lead to better management of Reclamation responsibilities. There are several areas where NASA satellite and modeling data may be useful to meet Reclamation's objectives for improved ET estimation. In this paper we outline one possible contribution to use NASA's data integration capability of the Land Information System (LIS) to provide a merger of observational (in situ and satellite) with physical process models to provide estimates of ET and other water availability outputs (e.g., runoff, soil moisture) retrospectively, in near real-time, and also providing short-term predictions.

The HYDROS mission is under development by NASA as part of its Earth System Science Pathfinder program. HYDROS is designed to provide global maps of the Earth's soil moisture and freeze/thaw state every 2-3 days, for weather and climate prediction, water and carbon cycle studies, natural hazards monitoring, and national security applications.

Water Management Applications is one of twelve elements in the Earth Science Enterprise National Applications Program. NASA Goddard Space Flight Center is supporting the Applications Program through partnering with other organizations to use NASA project results, such as from satellite instruments and Earth system models to enhance the organizations critical needs. The focus thus far has been: 1) estimating water storage including snowpack and soil moisture, 2) modeling and predicting water fluxes such as evapotranspiration (ET), precipitation and river runoff, and 3) remote sensing of water quality, including both point source (e.g., turbidity and productivity) and non-point source (e.g., land cover conversion such as forest to agriculture yielding higher nutrient runoff). The objectives of the partnering cover three steps of: 1) Evaluation, 2) Verification and Validation, and 3) Benchmark Report. We are working with the U.S. federal agencies including the Environmental Protection Agency (EPA), the Bureau of Reclamation (USBR) and the Department of Agriculture (USDA). We are using several of their Decision Support Systems (DSS) tools. This includes the DSS support tools BASINS used by EPA, Riverware and AWARDS ET ToolBox by USBR and SWAT by USDA and EPA. Regional application sites using NASA data across the US. are currently being eliminated for the DSS tools. The current NASA data emphasized thus far are from the Land Data Assimilation Systems WAS) and MODIS satellite products. We are currently in the first two steps of evaluation and verification validation. Water Management Applications is one of twelve elements in the Earth Science Enterprise s National Applications Program. NASA Goddard Space Flight Center is supporting the Applications Program through partnering with other organizations to use NASA project results, such as from satellite instruments and Earth system models to enhance the organizations critical needs. The focus thus far has been: 1) estimating water storage including snowpack and soil moisture, 2) modeling and predicting water fluxes such as evapotranspiration (ET), precipitation and river runoff, and 3) remote sensing of water quality, including both point source (e.g., turbidity and productivity) and non-point source (e.g., land cover conversion such as forest to agriculture yielding higher nutrient runoff). The objectives of the partnering cover three steps of 1) Evaluation, 2) Verification and Validation, and 3) Benchmark Report. We are working with the U.S. federal agencies the Environmental Protection Agency (EPA), the Bureau of Reclamation (USBR) and the Department of Agriculture (USDA). We are using several of their Decision Support Systems (DSS) tools. T us includes the DSS support tools BASINS used by EPA, Riverware and AWARDS ET ToolBox by USBR and SWAT by USDA and EPA. Regional application sites using NASA data across the US. are currently being evaluated for the DSS tools. The current NASA data emphasized thus far are from the Land Data Assimilation Systems (LDAS) and MODIS satellite products. We are currently in the first two steps of evaluation and verification and validation.

The HYDROS mission is under development by NASA as part of its Earth System Science Pathfinder (ESSP) program. HYDROS is designed to provide global maps of the Earth's soil moisture and freezel/thaw state every 2-3 days, for weather and climate prediction, water and carbon cycle studies, natural hazards monitoring, and national security applications. HYDROS uses a unique active and passive L-band microwave system that optimizes measurement accuracy, spatial resolution, and coverage. It provides measurements in nearly all weather conditions, regardless of solar illumination. The designs of the radar and radiometer electronics, antenna feedhorn and reflector, and science data system, are driven by specific mission and science objectives. These objectives impose requirements on the frequencies, polarizations, sampling, spatial resolution, and accuracy of the system. In this paper we describe the HYDROS mission requirements, baseline design, and measurement capabilities.

This is the first in a two-part series of manuscripts describing numerical experiments on the influence of 2-30 km striplike heterogeneity on wet and dry boundary layers coupled to the land surface. The strip-like heterogeneity is shown to dramatically alter the structure of the free-convective boundary layer by inducing significant organized circulations that modify turbulent statistics. The coupling with the land-surface modifies the circulations compared to previous studies using fixed surface forcing. Total boundary layer turbulence kinetic energy increases significantly for surface heterogeneity at scales between Lambda/z(sub i) = 4 and 9, however entrainment rates for all cases are largely unaffected by the strip-like heterogeneity.

We examine the influence of surface heterogeneity on boundary layers using a large-eddy simulation coupled to a land-surface model. Heterogeneity, imposed in strips varying from 2-30 km (1 less than lambda/z(sub i) less than 18), is found to dramatically alter the structure of the free convective boundary layer by inducing significant organized circulations. A conditional sampling technique, based on the scale of the surface heterogeneity (phase averaging), is used to identify and quantify the organized surface fluxes and motions in the atmospheric boundary layer. The impact of the organized motions on turbulent transport depends critically on the scale of the heterogeneity lambda, the boundary layer height zi and the initial moisture state of the boundary layer. Dynamical and scalar fields respond differently as the scale of the heterogeneity varies. Surface heterogeneity of scale 4 less than lamba/z(sub i) less than 9 induces the strongest organized flow fields (up, wp) while heterogeneity with smaller or larger lambda/z(sub i) induces little organized motion. However, the organized components of the scalar fields (virtual potential temperature and mixing ratio) grow continuously in magnitude and horizontal scale, as lambda/z(sub i) increases. For some cases, the organized motions can contribute nearly 100% of the total vertical moisture flux. Patch-induced fluxes are shown to dramatically impact point measurements that assume the time-average vertical velocity to be zero. The magnitude and sign of this impact depends on the location of the measurement within the region of heterogeneity.

This paper presents an approach to measuring all major components of the water cycle from space. The goal of the paper is to explore the concept of using a sensor-web of satellites to observe the global water cycle. The details of the required measurements and observation systems are therefore only an initial approach and will undergo future refinement, as their details will be highly important. Key elements include observation and evaluation of all components of the water cycle in terms of the storage of water-in the ocean, air, cloud and precipitation, in soil, ground water, snow and ice, and in lakes and rivers-and in terms of the global fluxes of water between these reservoirs. For each component of the water cycle that must be observed, the appropriate temporal and spatial scales of measurement are estimated, along with the some of the frequencies that have been used for active and passive microwave observations of the quantities. The suggested types of microwave observations are based on the heritage for such measurements, and some aspects of the recent heritage of these measurement algorithms are listed. The observational requirements are based on present observational systems, as modified by expectations for future needs. Approaches to the development of space systems for measuring the global water cycle can be based on these observational requirements.

The potential role of soil moisture initialization in seasonal forecasting is illustrated through ensembles of simulations with the NASA Seasonal-to-Interannual Prediction Project (NSIPP) model. For each boreal summer during 1997-2001, we generated two 16-member ensembles of 3-month simulations. The first, "AMIP-style" ensemble establishes the degree to which a perfect prediction of SSTs would contribute to the seasonal prediction of precipitation and temperature over continents. The second ensemble is identical to the first, except that the land surface is also initialized with "realistic" soil moisture contents through the continuous prior application (within GCM simulations leading up to the start of the forecast period) of a daily observational precipitation data set and the associated avoidance of model drift through the scaling of all surface prognostic variables. A comparison of the two ensembles shows that soil moisture initialization has a statistically significant impact on summertime precipitation and temperature over only a handful of continental regions. These regions agree, to first order, with regions that satisfy three conditions: (1) a tendency toward large initial soil moisture anomalies, (2) a strong sensitivity of evaporation to soil moisture, and (3) a strong sensitivity of precipitation to evaporation. The degree to which the initialization improves forecasts relative to observations is mixed, reflecting a critical need for the continued development of model parameterizations and data analysis strategies.

The Project for Intercomparison of Land Surface Parameterization Schemes (PILPS) Phase 2(e) showed that in cold regions the annual runoff production in Land Surface Schemes (LSSs) is closely related to the maximum snow accumulation, which in turn is controlled in large part by winter sublimation. To help further explain the relationship between snow cover, turbulent exchanges and runoff production, a simple equivalent model-(SEM) was devised to reproduce the seasonal and annual fluxes simulated by 13 LSSs that participated in PILPS Phase 2(e). The design of the SEM relates the annual partitioning of precipitation and energy in the LSSs to three primary parameters: snow albedo, effective aerodynamic resistance and evaporation efficiency. Isolation of each of the parameters showed that the annual runoff production was most sensitive to the aerodynamic resistance. The SEM was somewhat successful in reproducing the observed LSS response to a decrease in shortwave radiation and changes in wind speed forcings. SEM parameters derived from the reduced shortwave forcings suggested that increased winter stability suppressed turbulent heat fluxes over snow. Because winter sensible heat fluxes were largely negative, reductions in winter shortwave radiation imply an increase in annual average sensible heat.

The Global Land Data Assimilation System (GLDAS) is a global land parameterization that uses prescribed meteorology as forcing in order to determine regular gridded land surface states (temperature and moisture) and other properties (e.g. water and heat fluxes). In the present experiment, the assimilation of surface skin temperature is incorporated into the land parameterizations. The meteorological forcing was derived from the Goddard Earth Observing System (GEOS-3) Data Assimilation System (DAS) for the full year of 1998 GLDAS can use several land parameterizations, but here we use the Mosaic land surface model and the Common Land Model (CLM). TOVS surface temperature observations are assimilated into GLDAS. The TOVS observations are less frequent that observations used in previous experiments (ISCCP). The purpose of this presentation is to evaluate the impact of the TOVS assimilation on both Mosaic and CLM. We will especially consider the impact of coarse temporal observations on the assimilation and bias correction.

Accurate, spatially distributed surface temperatures are required for modeling evapotranspiration (ET) over agricultural fields under wide ranging conditions, including stressed and unstressed vegetation. Modeling approaches that use surface temperature observations, however, have the burden of estimating surface emissivities. Emissivity estimation, the subject of much recent research, is facilitated by observations in multiple thermal infrared bands. But it is nevertheless a difficult task. Using observations from a multiband thermal sensor, the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), estimated surface emissivities and temperatures are retrieved in two different ways: the temperature emissivity separation approach (TES) and the normalized emissivity approach (NEM). Both rely upon empirical relationships, but the assumed relationships are different. TES relies upon a relationship between the minimum spectral emissivity and the range of observed emissivities. NEM relies upon an assumption that at least one thermal band has a pre-determined emissivity (close to 1.0). The benefits and consequences of each approach will be demonstrated for two different landscapes: one in central Oklahoma, USA and another in southern New Mexico.

The Moderate Resolution Imaging Spectroradiometer (MODIS) on the Earth Observing System (EOS) Terra Mission began to produce data in February 2000. The Terra MODIS is in a sun-synchronous orbit going north to south in the daylight portion of the orbit crossing the equator at about 1030 hours local time. The spacecraft, instrument, and data systems are performing well and are producing a wide variety of data products useful for scientific and applications studies in relatively consistent fashion extending from November 2000 to the present. Within the approximately 40 MODIS data products, several are new and represent powerful and exciting capabilities such the ability to provide observations over the globe of fire occurrences, microphysical properties of clouds and sun-stimulated fluorescence from phytoplankton in the surface waters of the ocean. The remainder of the MODIS products exceed or, at a minimum, match the capabilities of products from heritage sensors such as, for example, the Advanced Very High Resolution Radiometer (AVHRR). Efforts are underway to provide data sets for the greater Earth science community and to improve access to these products at the various Distributed Active Archive Centers (DAAC's) or through Direct Broadcast (DB) stations. The MODIS instrument on the EOS Aqua mission should also be expected to be in orbit and functioning in the Spring of 2002. The Aqua spacecraft will operate in a sun-synchronous orbit going south to north in the daylight portion of the orbit crossing the equator at approximately 1330 hours local time. Subsequently the Aqua MODIS observations will substantially add to the capabilities of the Terra MODIS for environmental applications and global change studies.

On December 18, 1999, the Terra satellite was launched with a complement of five instruments including the Moderate Resolution Imaging Spectroradiometer (MODIS). Many geophysical products are derived from MODIS data including global snow-cover products. MODIS snow and ice products have been available through the National Snow and Ice Data Center (NSIDC) Distributed Active Archive Center (DAAC) since September 13, 2000. MODIS snow-cover products represent potential improvement to or enhancement of the currently-available operational products mainly because the MODIS products are global and 500-m resolution, and have the capability to separate most snow and clouds. Also the snow-mapping algorithms are automated which means that a consistent data set may be generated for long-term climate studies that require snow-cover information. Extensive quality assurance (QA) information is stored with the products. The MODIS snow product suite begins with a 500-m resolution, 2330-km swath snow-cover map which is then gridded to an integerized sinusoidal grid to produce daily and 8-day composite tile products. The sequence proceeds to a climate-modeling grid (CMG) product at about 5.6-km spatial resolution, with both daily and 8-day composite products. Each pixel of the CMG contains fraction of snow cover from 40 - 100%. Measured errors of commission in the CMG are low, for example, on the continent of Australia in the spring, they vary from 0.02 - 0.10%. Near-term enhancements include daily snow albedo and fractional snow cover. A case study from March 6, 2000, involving MODIS data and field and aircraft measurements, is presented to show some early validation work.

The Global Land Data Assimilation System (GLDAS) drives multiple uncoupled land surface models in order to produce optimal output fields of surface states in near-real time, globally, at 1/4 degree spatial resolution. These fields are then made available for coupled atmospheric model initialization and further research. One of the unique aspects of GLDAS is its ability to ingest both modeled and observation-derived forcing for running global scale land surface models. This paper compares results of runs forced by modeled and observed precipitation and shortwave radiation fields. Differences are examined and the impact of the observations on model skill is assessed.

Topographic effects on runoff generation have been documented observationally (e.g., Dunne and Black, 1970) and are the subject of the physically based rainfall-runoff model TOPMODEL (Beven and Kirkby, 1979; Beven, 1986a;b) and its extensions, which incorporate variable soil transmissivity effects (Sivapalan et al, 1987, Wood et al., 1988; 1990). These effects have been shown to exert significant control over the spatial distribution of runoff, soil moisture and evapotranspiration, and by extension, the latent and sensible heat fluxes

Some studies suggest that the proper initialization of soil moisture in a forecasting model may contribute significantly to the accurate prediction of seasonal precipitation, especially over mid-latitude continents. In order for the initialization to have any impact at all, however, two conditions must be satisfied: (1) the initial soil moisture anomaly must be "remembered" into the forecasted season, and (2) the atmosphere must respond in a predictable way to the soil moisture anomaly. In our previous studies, we identified the key land surface and atmospheric properties needed to satisfy each condition. Here, we tie these studies together with an analysis of an ensemble of seasonal forecasts. Initial soil moisture conditions for the forecasts are established by forcing the land surface model with realistic precipitation prior to the start of the forecast period. As expected, the impacts on forecasted precipitation (relative to an ensemble of runs that do not utilize soil moisture information) tend to be localized over the small fraction of the earth with all of the required land and atmosphere properties.

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.

Realistic representation of the land surface is crucial in global climate modeling (GCM). Recently, the Mosaic land-surface Model (LSM) has been driven off-line using GEOS DAS (Goddard Earth Observing System Data Assimilation System) atmospheric forcing, forming the Off-line Land-surface Global Assimilation (OLGA) system. This system provides a computationally efficient test bed for land surface data assimilation. Here, we validate the OLGA simulation of surface processes and the assimilation of ISCCP surface temperatures. Another component of this study as the incorporation of the Physical-space Statistical Analysis System (PSAS) into OLGA, in order to assimilate surface temperature observations from the International Satellite Cloud Climatology Project (ISCCP). To counteract the subsequent forcing of the analyzed skin temperature back to the initial state following the analysis. incremental bias correction (IBC) was included in the assimilation. The IBC scheme effectively removed the time mean bias, but did not remove him in the mean diurnal cycle. Therefore, a diurnal him correction (DBC) scheme was developed, where the time-dependent bias was modeled with a sine wave parameterization. In addition, quality control of the ISCCP data and anisotropic temperature correction were implemented in PSAS. Preliminary results showed a substantial impact from the inclusion of PSAS and DBC that was visible in the surface meteorology fields and energy budget. Also, the monthly mean diurnal cycle from the experiment closely matched the diurnal cycle from the observations.

Seasonal snow cover in South America was examined in this study using passive microwave satellite data from the Special Sensor Microwave Imagers (SSM/I) on board Defense Meteorological Satellite Program (DMSP) satellites. For the period from 1992-1998, both snow cover extent and snow depth (snow mass) were investigated during the winter months (May-August) in the Patagonia region of Argentina. Since above normal temperatures in this region are typically above freezing, the coldest winter month was found to be not only the month having the most extensive snow cover but also the month having the deepest snows. For the seven-year period of this study, the average snow cover extent (May-August) was about 0.46 million sq km and the average monthly snow mass was about 1.18 x 10(exp 13) kg. July 1992 was the month having the greatest snow extent (nearly 0.8 million sq km) and snow mass (approximately 2.6 x 10(exp 13) kg).

The rate of carbon accumulation by terrestrial plant communities in a process-level, mechanistic modeling is the difference of the rate of gross photosynthesis by a canopy (A(sub g)) and autotrophic respiration (R) of the stand. Observations for different biomes often show that R to be a large and variable fraction of A(sub g), ca. 35% to 75%, although other studies suggest the ratio of R and A(sub g) to be less variable. Here, R has been calculated according to the two compartment model as being the sum of maintenance and growth components. The maintenance respiration of foliage and living fine roots for different biomes has been determined objectively from observed nitrogen content of these organs. The sapwood maintenance respiration is based on pipe theory, and checked against an independently derived equation considering sapwood biomass and its maintenance coefficient. The growth respiration has been calculated from the difference of A(sub g) and maintenance respiration. The A(sub g) is obtained as the product of biome-specific radiation use efficiency for gross photosynthesis under unstressed conditions and intercepted photosynthetically active radiation, and adjusted for stress. Calculations have been done using satellite and ground observations for 36 consecutive months (1987-1989) over large contiguous areas (ca. 10(exp 5) sq km) of boreal forests, crop land, temperate deciduous forest, temperate grassland, tropical deciduous forest, tropical evergreen forest, tropical savanna, and tundra. The ratio of annual respiration and gross photosynthesis, (R/A(sub g)), is found to be 0.5-0.6 for temperate and cold adopted biome areas, but somewhat higher for tropical biome areas (0.6-0.7). Interannual variation of the fluxes is found to be generally less than 15%. Calculated fluxes are compared with observations and several previous estimates. Results of sensitivity analysis are presented for uncertainties in parameterization and input data. It is found that uncertainty in determining maintenance respiration for tropical biomes is such that R/A(sub g) for these biomes could be similar to that for temperate biomes.

This study presents two new modeling strategies. First, a methodology for representing the physical process of stormflow within a TOPMODEL framework is developed. In using this approach, discharge at quickflow time scales is simulated and a fuller depiction of hydrologic activity is brought about. Discharge of water from the vadose zone is permitted in a physically realistic manner without a priori assumption of the level within the soil column at which stormflow saturation can take place. Determination of the stormflow contribution to discharge is made using the equation for groundwater flow. No new parameters are needed. Instead, regions of near saturation that develop during storm events, producing vertical recharge, are allowed to contribute to soil column discharge. These stormflow contributions to river runoff, as for groundwater flow contributions, are a function of catchment topography and local hydraulic conductivity at the depth of these regions of near saturation. The second approach improves groundwater flow response through a reduction of porosity and field capacity with depth in the soil column. Large storm events are better captured and a more dynamic water table develops with application of this modified soil column profile (MSCP). The MSCP predominantly reflects soil depth differences in upland and lowland regions of a watershed. Combined, these two approaches - stormflow and the MSCP - provide a more accurate representation of the time scales at which soil column discharge responds and a more complete depiction of hydrologic activity. Storm events large and small are better simulated, and some of the biases previously evident in TOPMODEL simulations are reduced.

Moderate Resolution Imaging Spectroradiometer (MODIS) snow-cover maps have been available since September 13, 2000. These products, at 500-m spatial resolution, are available through the National Snow and Ice Data Center Distributed Active Archive Center in Boulder, Colorado. By the 2001-02 winter, 5-km climate-modeling grid (CMG) products will be available for presentation of global views of snow cover and for use in climate models. All MODIS snow-cover products are produced from automated algorithms that map snow in an objective manner. In this paper, we describe the MODIS snow products, and show snow maps from the fall of 2000 in North America.

On December 18, 1999, the Terra satellite was launched with a complement of five instruments including the Moderate Resolution Imaging Spectroradiometer (MODIS). Many geophysical products are derived from MODIS data including global snow-cover products. These products have been available through the National Snow and Ice Data Center (NSIDC) Distributed Active Archive Center (DAAC) since September 13, 2000. MODIS snow-cover products represent potential improvement to the currently available operation products mainly because the MODIS products are global and 500-m resolution, and have the capability to separate most snow and clouds. Also the snow-mapping algorithms are automated which means that a consistent data set is generated for long-term climates studies that require snow-cover information. Extensive quality assurance (QA) information is stored with the product. The snow product suite starts with a 500-m resolution swath snow-cover map which is gridded to the Integerized Sinusoidal Grid to produce daily and eight-day composite tile products. The sequence then proceeds to a climate-modeling grid product at 5-km spatial resolution, with both daily and eight-day composite products. A case study from March 6, 2000, involving MODIS data and field and aircraft measurements, is presented. Near-term enhancements include daily snow albedo and fractional snow cover.

The net primary productivity (C) or annual rate of carbon accumulation per unit ground area by terrestrial plant communities is the difference of the rate of gross photosynthesis (A(sub g)) and autotrophic respiration (R) per unit ground area. Although available observations show that R is a large and variable fraction of A(sub g), viz., 0.3 to 0.7, it is generally recognized that much uncertainties exist in this fraction due to difficulties associated with the needed measurements. Additional uncertainties arise when these measurements are extrapolated to regional or global land surface using empirical equations, for example, using regression equations relating C to mean annual precipitation and air temperature. Here, a process-based approach has been taken to calculate A(sub g) and R using satellite and ancillary data. A(sub g) has been expressed as a product of radiation use efficiency, magnitude of intercepted photosynthetically active radiation (PAR), and normalized by stresses due to soil water shortage and air temperature away from the optimum range. A biophysical model has been used to determine the radiation use efficiency from the maximum rate of carbon assimilation by a leaf, foliage temperature, and the fraction of diffuse PAR incident on a canopy. All meteorological data (PAR, air temperature, precipitation, etc.) needed for the calculation are derived from satellite observations, while a land use, land cover data (based on satellite and ground measurements) have been used to assess the maximum rate of carbon assimilation by a leaf of varied cover type based on field measurements. R has been calculated as the sum of maintenance and growth components. The maintenance respiration of foliage and live fine roots at a standard temperature of different land cover has been determined from their nitrogen content using field and satellite measurements, while that of living fraction of woody stem (viz., sapwood) from the seasonal maximum leaf area index as determined from satellite observations. These maintenance respiration values were then adjusted to that corresponding to air temperature according to a prescribed non-linear variation of respiration with temperature. The growth respiration has been calculated from the difference of Ag and maintenance respiration, according to the two-compartment model. The results of calculations will be reported for 36 consecutive months (1987-1989) over large contiguous areas (ca. 10(exp 5) sq km) Of agricultural land and tropical humid evergreen forests, and compared with available field data.

Land surface models (LSMs) play a critical role in the simulation of climate, for they determine the character of a large fraction of the atmosphere's lower boundary. The LSM partitions the net radiative energy at the land surface into sensible heat, latent heat, and energy storage, and it partitions incident precipitation water into evaporation, runoff, and water storage. Numerous modeling experiments and the existing (though very scant) observational evidence suggest that variations in these partitionings can feed back on the atmospheric processes that induce them. This land-atmosphere feedback can in turn have a significant impact on the generation of continental precipitation. For this and other reasons (including the role of the land surface in converting various atmospheric quantities, such as precipitation, into quantities of perhaps higher societal relevance, such as runoff), many modeling groups are placing a high emphasis on improving the treatment of land surface processes in their models. LSMs have evolved substantially from the original bucket model of Manabe et al. This evolution, which is still ongoing, has been documented considerably. The present paper also takes a look at the evolution of LSMs. The perspective here, though, is different - the evolution is considered strictly in terms of the 'balance' between the formulations of evaporation and runoff processes. The paper will argue that a proper balance is currently missing, largely due to difficulties in treating subgrid variability in soil moisture and its impact on the generation of runoff.

There are several hemispheric-scale satellite-derived snow-cover maps available, but none has been fully validated. For the period October 23 - December 25, 2000, we compare snow maps of North America derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the National Oceanic and Atmospheric Administration (NOAA) National Operational Hydrologic Remote Sensing Center (NOHRSC), which both rely on satellite data from the visible and near-infrared parts of the spectrum; we also compare MODIS and Defense Meteorological Satellite Program (DMSP) Special Sensor Microwave/Imager (SSM/I) passive-microwave snow maps. The maps derived from visible and near-infrared data are more accurate for mapping snow cover than are the passive-microwave-derived maps, however discrepancies exist as to the location and extent of the snow cover among those maps. The large (approx. 30 km) footprint of the SSM/I data and the difficulty in distinguishing wet and shallow snow from wet or snow-free ground, reveal differences up to 5.32 million sq km in the amount of snow mapped using MODIS versus SSM/I data. Algorithms that utilize both visible and passive-microwave data, which would take advantage of the all-weather mapping ability of the passive-microwave data, will be refined following the launch of the Advanced Microwave Scanning Radiometer (AMSR) in the fall of 2001.

Moderate Resolution Imaging Spectroradiometer (MODIS) snow-cover maps have been available since September 13, 2000. These products, at 500 m spatial resolution, are available through the National Snow and Ice Data Center Distributed Active Archive Center in Boulder, Colorado. By the 2001-02 winter, 5 km climate-modeling grid (CMG) products will be available for presentation of global views of snow cover and for use in climate models. All MODIS snow-cover products are produced from automated algorithms that map snow in an objective manner. In this paper, we describe the MODIS snow products, and show snow maps from the fall of 2000 in North America.

Land-atmosphere feedback, by which (for example) precipitation-induced moisture anomalies at the land surface affect the overlying atmosphere and thereby the subsequent generation of precipitation, has been examined and quantified with many atmospheric general circulation models (AGCMs). Generally missing from such studies, however, is an indication of the extent to which the simulated feedback strength is model dependent. Four modeling groups have recently performed a highly controlled numerical experiment that allows an objective inter-model comparison of land-atmosphere feedback strength. The experiment essentially consists of an ensemble of simulations in which each member simulation artificially maintains the same time series of surface prognostic variables. Differences in atmospheric behavior between the ensemble members then indicates the degree to which the state of the land surface controls atmospheric processes in that model. A comparison of the four sets of experimental results shows that feedback strength does indeed vary significantly between the AGCMs.

Assimilation of remotely sensed land surface observations into regional to global scale numerical models have the potential to significantly advance our ability, to assess, understand, and predict surface water, energy, and carbon cycles. This session seeks to assess the state-of-the-art in data assimilation methods for integrating land surface remote sensing and modeling, with a focus on practical applications and techniques. Assimilated land surface variables of interest include (but are not limited to, soil moisture, surface temperature, snowpack, streamflow, vegetation dynamics, and carbon storage. Contributions describing the development of practical land surface data assimilation methods, multivariate land surface data assimilation strategies, evaluation of the required accuracy and resolution of remote sensing observations, the effects of scale, process complexity, and uncertainty on data assimilation, and the optimal treatment of model and observation errors are encouraged.

The Mosaic Land-surface Model (LSM) has been included into the current GEOS Data Assimilation System (DAS). The LSM uses a more advanced representation of physical processes than previous versions of the GEOS DAS, including the representation of sub-grid heterogeneity of the land-surface through the Mosaic approach. As a first approximation, Mosaic assumes that all similar surface types within a grid-cell can be lumped together as a single'tile'. Within one GCM grid-cell, there might be 1 - 5 different tiles or surface types. All tiles are subjected to the grid-scale forcing (radiation, air temperature and specific humidity, and precipitation), and the sub-grid variability is a function of the tile characteristics. In this paper, we validate the LSM sub-grid scale variability (tiles) using a variety of surface observing stations from the Southern Great Plains (SGP) site of the Atmospheric Radiation Measurement (ARM) Program. One of the primary goals of SGP ARM is to study the variability of atmospheric radiation within a G,CM grid-cell. Enough surface data has been collected by ARM to extend this goal to sub-grid variability of the land-surface energy and water budgets. The time period of this study is the Summer of 1998 (June I - September 1). The ARM site data consists of surface meteorology, energy flux (eddy correlation and bowen ratio), soil water observations spread over an area similar to the size of a G-CM grid-cell. Various ARM stations are described as wheat and alfalfa crops, pasture and range land. The LSM tiles considered at the grid-space (2 x 2.5) nearest the ARM site include, grassland, deciduous forests, bare soil and dwarf trees. Surface energy and water balances for each tile type are compared with observations. Furthermore, we will discuss the land-surface sub-grid variability of both the ARM observations and the DAS.

Land-surface processes in a data assimilation system influence the lower troposphere and must be properly represented. With the recent incorporation of the Mosaic Land-surface Model (LSM) into the GEOS Data Assimilation System (DAS), the detailed land-surface processes require strict validation. While global data sources can identify large-scale systematic biases at the monthly timescale, the diurnal cycle is difficult to validate. Moreover, global data sets rarely include variables such as evaporation, sensible heat and soil water. Intensive field experiments, on the other hand, can provide high temporal resolution energy budget and vertical profile data for sufficiently long periods, without global coverage. Here, we evaluate the GEOS DAS against several intensive field experiments. The field experiments are First ISLSCP Field Experiment (FIFE, Kansas, summer 1987), Cabauw (as used in PILPS, Netherlands, summer 1987), Atmospheric Radiation Measurement (ARM, Southern Great Plains, winter and summer 1998) and the Surface Heat Budget of the Arctic Ocean (SHEBA, Arctic ice sheet, winter and summer 1998). The sites provide complete surface energy budget data for periods of at least one year, and some periods of vertical profiles. This comparison provides a detailed validation of the Mosaic LSM within the GEOS DAS for a variety of climatologic and geographic conditions.

Accurate soil moisture information has proved to be important to climate simulations and climate and weather forecasts. However, many difficulties exist that limit our understanding of soil moisture distribution and variability. One of them is the lack of accurate precipitation with appropriate spatial and temporal resolution. Precipitation as an input forcing to the land surface greatly influences soil moisture characteristics and variability. To improve precipitation data quality, an algorithm has been developed to generate a spatially and temporally continuous 3-hourly global precipitation data for the period of 1987 to present. This precipitation product is a combination of the precipitation from Special Sensor Microwave Imager (SSMI) with the Goddard Earth Observing System-1 Data Assimilation System (GEOS-1 DAS) employing a Physical-space Statistical Analysis System (PSAS). In this study we investigate the impact of this merged/analyzed precipitation data on the global soil moisture variability using an Off-line Land-surface GEOS Assimilation (OLGA) system. Two OLGA integrations starting from 1987 to 1993 are performed forced with the analyzed and GEOS-1 DAS precipitation respectively. We examine the spatial and temporal characteristics of soil moisture variability in response to the analyzed precipitation. The influence of this merged precipitation on the soil moisture variability and regional hydrological budget is estimated throughout the comparison with the results forced with the GEOS-1 DAS precipitation only. In the OLGA the sut@-grid scale horizontal heterogeneity is explicitly represented on the tile space. This provides a means to assess the role of the surface moisture heterogeneity in the interaction with the surface atmosphere and surface hydrological budget, and to validate OLGA results at tile space with in situ observation. ABRACOS (Anglo-Brazilian Amazonian Climate Observation Study), FIFE (First ISLSCP Field Experiment) I and HAPEX data will be used in the validation.

Recently, NASA Goddard Earth Observing System (GEOS-1) reanalysis data has been used to provide forcing for the Koster and Suarez Mosaic Land-surface Model (LSM). The LSM was integrated off-line at all global land points for the period of 1983 - 1995 by the Off-line Land- surface GEOS Assimilation system (OLGA). Here, we compare the interannual variability of OLGA, GEOS-1 and surface observing stations temperature and moisture. Particular attention is given to the United States because of the extreme seasons of 1988 and 1993. Furthermore, the comparison of OLGA is extended to include the analysis of data on the'tiles' (different surface types) in the Mosaic LSM. Results indicate that the GEOS-1 near-surface temperature and moisture reasonably represents the interannual variability in more normal years. However, OLGA also simulates the extreme drought and floods years well. The analysis of the tile information shows that the "Bare soil" surface type is most sensitive to the climate extremes. Off-line testing has provided valuable information on the performance of the Mosaic LSM prior to its incorporation into the new version of the GEOS Data Assimilation System and the integration of a new long reanalysis.

Several recent studies have shown that much can be learned by comparing grid-point data from a data assimilation system with in-situ observations from field experiments. While the surface heterogeneity is acknowledged in these studies, they lack quantitative representations of the influence of heterogeneity on the near-surface meteorology and surface hydrologic and energy balance. Here, we use the Betts and Ball FIFE site-averaged data. Standard deviations of the site-average will provide an estimate of the FIFE site heterogeneity. Recently, the Mosaic Land-Surface Model (LSM) has been incorporated into the Goddard Earth Observing System (GEOS) Data Assimilation System (DAS). The Mosaic LSM computes the surface energy and hydrologic balance for nine distinct surface types at each grid-point. Each surface type is proportionally weighted to determine the mean grid point properties. Hence, we can compare modeled and observed grid-point variability in addition to the mean properties. Also, assimilated data sets created with and without the LSM are compared. The results indicate the importance of including quantitative estimates of heterogeneity in the analysis of the land surface hydrology and energy balances in assimilation systems.

Land surface hydrology for the Off-line Land-surface GEOS Analysis (OLGA) system and Goddard Earth Observing System (GEOS-1) Data Assimilation System (DAS) has been examined using a river routing model. The GEOS-1 DAS land-surface parameterization is very simple, using an energy balance prediction of surface temperature and prescribed soil water. OLGA uses near-surface atmospheric data from the GEOS-1 DAS to drive a more comprehensive parameterization of the land-surface physics. The two global systems are evaluated using a global river routing model. The river routing model uses climatologic surface runoff from each system to simulate the river discharge from global river basins, which can be compared to climatologic river discharge. Due to the soil hydrology, the OLGA system shows a general improvement in the simulation of river discharge compared to the GEOS-1 DAS. Snowmelt processes included in OLGA also have a positive effect on the annual cycle of river discharge and source runoff. Preliminary tests of a coupled land-atmosphere model indicate improvements to the hydrologic cycle compared to the uncoupled system. The river routing model has provided a useful tool in the evaluation of the GCM hydrologic cycle, and has helped quantify the influence of the more advanced land surface model.

The land surface component of the hydrological cycle is fundamental to the overall functioning of the atmospheric and climate processes. Spatially and temporally variable rainfall and available energy, combined with land surface heterogeneity cause complex variations in all processes related to surface hydrology. The characterization of the spatial and temporal variability of water and energy cycles are critical to improve our understanding of land surface-atmosphere interaction and the impact of land surface processes on climate extremes. Because the accurate knowledge of these processes and their variability is important for climate predictions, most Numerical Weather Prediction (NWP) centers have incorporated land surface schemes in their models. However, errors in the NWP forcing accumulate in the surface and energy stores, leading to incorrect surface water and energy partitioning and related processes. This has motivated the NWP to impose ad hoc corrections to the land surface states to prevent this drift. A proposed methodology is to develop Land Data Assimilation schemes (LDAS), which are uncoupled models forced with observations, and not affected by NWP forcing biases. The proposed research is being implemented as a real time operation using an existing Surface Vegetation Atmosphere Transfer Scheme (SVATS) model at a 40 km degree resolution across the United States to evaluate these critical science questions. The model will be forced with real time output from numerical prediction models, satellite data, and radar precipitation measurements. Model parameters will be derived from the existing GIS vegetation and soil coverages. The model results will be aggregated to various scales to assess water and energy balances and these will be validated with various in-situ observations.

The Earth's land surface, including its biomass, is an integral part of the Earth's weather and climate system. Land surface heterogeneity, such as the type and amount of vegetative covering., has a profound effect on local weather variability and therefore on regional variations of the global climate. Surface conditions affect local weather and climate through a number of mechanisms. First, they determine the re-distribution of the net radiative energy received at the surface, through the atmosphere, from the sun. A certain fraction of this energy increases the surface ground temperature, another warms the near-surface atmosphere, and the rest evaporates surface water, which in turn creates clouds and causes precipitation. Second, they determine how much rainfall and snowmelt can be stored in the soil and how much instead runs off into waterways. Finally, surface conditions influence the near-surface concentration and distribution of greenhouse gases such as carbon dioxide. The processes through which these mechanisms interact with the atmosphere can be modeled mathematically, to within some degree of uncertainty, on the basis of underlying physical principles. Such a land surface model provides predictive capability for surface variables including ground temperature, surface humidity, and soil moisture and temperature. This information is important for agriculture and industry, as well as for addressing fundamental scientific questions concerning global and local climate change. In this study we apply a methodology known as tangent linear modeling to help us understand more deeply, the behavior of the Mosaic land surface model, a model that has been developed over the past several years at NASA/GSFC. This methodology allows us to examine, directly and quantitatively, the dependence of prediction errors in land surface variables upon different vegetation conditions. The work also highlights the importance of accurate soil moisture information. Although surface variables are predicted imperfectly due to inherent uncertainties in the modeling process, our study suggests how satellite observations can be combined with the model, through land surface data assimilation, to improve their prediction.

Meteorological measurements in the Walnut Gulch catchment in Arizona were used to synthesize a distributed, hourly-average time series of data across a 26.9 by 12.5 km area with a grid resolution of 480 m for a continuous 18-month period which included two seasons of monsoonal rainfall. Coupled surface-atmosphere model runs established the acceptability (for modelling purposes) of assuming uniformity in all meteorological variables other than rainfall. Rainfall was interpolated onto the grid from an array of 82 recording rain gauges. These meteorological data were used as forcing variables for an equivalent array of stand-alone Biosphere-Atmosphere Transfer Scheme (BATS) models to describe the evolution of soil moisture and surface energy fluxes in response to the prevalent, heterogeneous pattern of convective precipitation. The calculated area-average behaviour was compared with that given by a single aggregate BATS simulation forced with area-average meteorological data. Heterogeneous rainfall gives rise to significant but partly compensating differences in the transpiration and the intercepted rainfall components of total evaporation during rain storms. However, the calculated area-average surface energy fluxes given by the two simulations in rain-free conditions with strong heterogeneity in soil moisture were always close to identical, a result which is independent of whether default or site-specific vegetation and soil parameters were used. Because the spatial variability in soil moisture throughout the catchment has the same order of magnitude as the amount of rain failing in a typical convective storm (commonly 10% of the vegetation's root zone saturation) in a semi-arid environment, non-linearitv in the relationship between transpiration and the soil moisture available to the vegetation has limited influence on area-average surface fluxes.