Search

The focus of this chapter is progress in hydrology for the last 100 years. During this period, we have seen a marked transition from practical engineering hydrology to fundamental developments in hydrologic science, including contributions to Earth system science. The first three sections in this chapter review advances in theory, observations, and hydrologic prediction. Building on this foundation, the growth of global hydrology, land-atmosphere interactions and coupling, ecohydrology, and water management are discussed, as well as a brief summary of emerging challenges and future directions. Although the review attempts to be comprehensive, the chapter offers greater coverage on surface hydrology and hydrometeorology for readers of this American Meteorological Society (AMS) Monograph.

The advent of commercially available quantum computers has marked the beginning of quantum computing as a reality. Both quantum gate and annealing computers have been released by major computer hardware companies. In this work, the D-Wave 2XTM quantum annealing computer housed at the NASA Advanced Systems computational facility is investigated to accelerate Machine Learning (ML) for image registration. NASA collects large amounts of images over the globe remotely using space-based monitoring. Images of a fixed areas of the land surface are taken over time. Due to the orbit of the sensors, the viewing angles deviate slightly, and it is necessary to align or register the images precisely to create image time series over the land surface. Unaligned images can lead to substantial analysis errors. These time-series are then used in modeling Earth Systems models such as hydrological, weather, and carbon monitoring models. In this work, we consider the Moderate Resolution Image Spectrometer (MODIS) data collected by the NASA's terra satellite. Artificial Neural Networks (ANNs) is a natural fit for ML modelling of images. Several successes have been reported using machine learning related to image processing. We investigate the use of ML to register MODIS images. ANNs are investigated in combination with a Restricted Boltzmann Machines (RBM) as an auto-encoder. We will present results showing the accuracy and efficiency of this approach.The D-Wave 2XTM quantum annealer samples the ground-state wave-function of a spin-Ising systems with quadratic interactions between qubits and a Chimera connectivity. The system sits in a ~15 mK thermal bath. One can think of the system as being placed in the ground state initially and subject to thermal excitations governed by Boltzmann statistics. If this is assumed true, one can use the statistics from the D-Wave 2XTM to train RBMs. Generating statistics for training Boltzmann machines is an NP-hard problem and constitutes the largest compute cost. We investigate the use of the D-Wave 2XTM to accelerate the training of the RBMs in our ANNs and report on the results.

The effective applications of land surface models (LSMs) and hydrologic models pose a varied set of data input and processing needs, ranging from ensuring consistency checks to more derived data processing and analytics. This article describes the development of the Land surface Data Toolkit (LDT), which is an integrated framework designed specifically for processing input data to execute LSMs and hydrological models. LDT not only serves as a preprocessor to the NASA Land Information System (LIS), which is an integrated framework designed for multi-model LSM simulations and data assimilation (DA) integrations, but also as a land-surface-based observation and DA input processor. It offers a variety of user options and inputs to processing datasets for use within LIS and stand-alone models. The LDT design facilitates the use of common data formats and conventions. LDT is also capable of processing LSM initial conditions and meteorological boundary conditions and ensuring data quality for inputs to LSMs and DA routines. The machine learning layer in LDT facilitates the use of modern data science algorithms for developing data-driven predictive models. Through the use of an object-oriented framework design, LDT provides extensible features for the continued development of support for different types of observational datasets and data analytics algorithms to aid land surface modeling and data assimilation.

A Saharan air layer (SAL) event associated with a nondeveloping African easterly wave (AEW) over the main development region of the eastern Atlantic was sampled by the NASA Global Hawk aircraft on 24-25 August 2013 during the NASA Hurricane and Severe Storm Sentinel (HS3) campaign and was simulated with the NASA Unified Weather Research and Forecasting (NU-WRF) Model. Airborne, ground-based, and spaceborne measurements were used to evaluate the model performance. The microphysical and radiative effects of dust and other aerosols on the SAL structure and environment were investigated with the factor-separation method. The results indicate that relative to a simulation without dust-radiative and microphysical impacts, Saharan dust and other aerosols heated the SAL air mainly through shortwave heating by the direct aerosol-radiation (AR) effect, resulting in a warmer (up to 0.6 K) and drier (up to 5% RH reduction) SAL and maintaining the strong temperature inversion at the base of the SAL in the presence of predominant longwave cooling. Radiative heating of the dust accentuated a vertical circulation within the dust layer, in which air rose (sank) in the northern (southern) portions of the dust layer. Furthermore, above and to the south of the dust layer, both the microphysical and radiative impacts of dust tended to counter the vertical motions associated with the Hadley circulation, causing a small weakening and southward shift of convection in the intertropical convergence zone (ITCZ) and reduced anvil cloud to the north. Changes in moisture and cloud/precipitation hydrometeors were largely driven by the dust induced changes in vertical motion. Dust strengthened the African easterly jet by up to ~1ms(exp -1) at the southern edge of the jet, primarily through the AR effect, and produced modest increases in vertical wind shear within and in the vicinity of the dust layer. These modulations of the SAL and AEW environment clearly contributed to the nondevelopment of this AEW.

Accurate quantification of the terrestrial evapotranspiration (ET) components of plant transpiration (T), soil evaporation (E) and evaporation of the intercepted water (I) is necessary for improving our understanding of the links between the carbon and water cycles. Recent studies have noted that, among the modeled estimates, large disagreements exist in the relative contributions of T, E and I to the total ET. As these models are often used in data assimilation environments for incorporating and extending ET relevant remote sensing measurements, understanding the sources of inter-model differences in ET components is also necessary for improving the utilization of such remote sensing measurements. This study quantifies the contributions of two key factors explaining inter-model disagreements to the uncertainty in total ET: (1) contribution of the local partitioning and (2) regional distribution of ET. The analysis is conducted by using outputs from a suite of land surface models in the North American Land Data Assimilation System (NLDAS) configuration. For most of these models, transpiration is the dominant component of the ET partition. The results indicate that the uncertainty in local partitioning dominates the inter-model spread in modeled soil evaporation E. The inter-model differences in T are dominated by the uncertainty in the distribution of ET over the Eastern U.S. and the local partitioning uncertainty in the Western U.S. The results also indicate that uncertainty in the T estimates is the primary driver of total ET errors. Over the majority of the U.S., the contribution of the two factors of uncertainty to the overall uncertainty is non-trivial.

Multimodel ensembles are often used to produce ensemble mean estimates that tend to have increased simulation skill over any individual model output. If multimodel outputs are too similar, an individual LSM would add little additional information to the multimodel ensemble, whereas if the models are too dissimilar, it may be indicative of systematic errors in their formulations or configurations. The article presents a formal similarity assessment of the North American Land Data Assimilation System (NLDAS) multimodel ensemble outputs to assess their utility to the ensemble, using a confirmatory factor analysis. Outputs from four NLDAS Phase 2 models currently running in operations at NOAA/NCEP and four new/ upgraded models that are under consideration for the next phase of NLDAS are employed in this study. The results show that the runoff estimates from the LSMs were most dissimilar whereas the models showed greater similarity for root zone soil moisture, snow water equivalent, and terrestrial water storage. Generally, the NLDAS operational models showed weaker association with the common factor of the ensemble and the newer versions of the LSMs showed stronger association with the common factor, with the model similarity increasing at longer time scales. Trade-offs between the similarity metrics and accuracy measures indicated that the NLDAS operational models demonstrate a larger span in the similarity-accuracy space compared to the new LSMs. The results of the article indicate that simultaneous consideration of model similarity and accuracy at the relevant time scales is necessary in the development of multimodel ensemble.

Rainfall variability in the Tigris Euphrates headwaters is a result of interaction between topography and meteorological features at a range of spatial scales. Here, the Weather Research and Forecasting (WRF) Model, driven by the NCEP-DOE AMIP-II reanalysis (R-2), has been implemented to better understand these interactions. Simulations were performed over a domain covering most of the Middle East. The extended simulation period (1983 - 2013) enables us to study seasonality, interannual variability, spatial variability, and extreme events of rainfall. Results showed that the annual cycle of precipitation produced by WRF agrees much more closely with observations than does R-2. This was particularly evident during the transition months of April and October, which were further examined to study the underlying physical mechanisms. In both months, WRF improves representation of interannual variability relative to R-2, with a substantially larger benefit in April. This improvement results primarily from WRFs ability to resolve two low-level, terrain-induced flows in the region that are either absent or weak in R-2: one parallel to the western edge of the Zagros Mountains, and one along the east Turkish highlands. The first shows a complete reversal in its direction during wet and dry days, when flowing southeasterly it transports moisture from the Persian Gulf to the region, and when flowing northwesterly it blocks moisture and transports it away from the region. The second is more directly related to synoptic-scale systems and carries moist, warm air from the Mediterranean and Red Seas toward the region. The combined contribution of these flows explains about 50 of interannual variability in both WRF and observations for April and October precipitation.

Accurate specification of the model error covariances in data assimilation systems is a challenging issue. Ensemble land data assimilation methods rely on stochastic perturbations of input forcing and model prognostic fields for developing representations of input model error covariances. This article examines the limitations of using a single forcing dataset for specifying forcing uncertainty inputs for assimilating snow depth retrievals. Using an idealized data assimilation experiment, the article demonstrates that the use of hybrid forcing input strategies (either through the use of an ensemble of forcing products or through the added use of the forcing climatology) provide a better characterization of the background model error, which leads to improved data assimilation results, especially during the snow accumulation and melt-time periods. The use of hybrid forcing ensembles is then employed for assimilating snow depth retrievals from the AMSR2 (Advanced Microwave Scanning Radiometer 2) instrument over two domains in the continental USA with different snow evolution characteristics. Over a region near the Great Lakes, where the snow evolution tends to be ephemeral, the use of hybrid forcing ensembles provides significant improvements relative to the use of a single forcing dataset. Over the Colorado headwaters characterized by large snow accumulation, the impact of using the forcing ensemble is less prominent and is largely limited to the snow transition time periods. The results of the article demonstrate that improving the background model error through the use of a forcing ensemble enables the assimilation system to better incorporate the observational information.

Seasonal agricultural drought monitoring systems, which rely on satellite remote sensing and land surface models (LSMs), are important for disaster risk reduction and famine early warning. These systems require the best available weather inputs, as well as a long-term historical record to contextualize current observations. This article introduces the Famine Early Warning Systems Network (FEWS NET) Land Data Assimilation System (FLDAS), a custom instance of the NASA Land Information System (LIS) framework. The FLDAS is routinely used to produce multi-model and multi-forcing estimates of hydro-climate states and fluxes over semi-arid, food insecure regions of Africa. These modeled data and derived products, like soil moisture percentiles and water availability, were designed and are currently used to complement FEWSNETs operational remotely sensed rainfall, evapotranspiration, and vegetation observations. The 30+ years of monthly outputs from the FLDAS simulations are publicly available from the NASA Goddard Earth Science Data and Information Services Center (GES DISC) and recommended for use in hydroclimate studies, early warning applications, and by agro-meteorological scientists in Eastern, Southern, and Western Africa.

The objective of the North American Land Data Assimilation System (NLDAS) is to provide best available estimates of near-surface meteorological conditions and soil hydrological status for the continental United States. To support the ongoing efforts to develop data assimilation (DA) capabilities for NLDAS, the results of Gravity Recovery and Climate Experiment (GRACE) DA implemented in a manner consistent with NLDAS development are presented. Following previous work, GRACE terrestrial water storage (TWS) anomaly estimates are assimilated into the NASA Catchment land surface model using an ensemble smoother. In contrast to many earlier GRACE DA studies, a gridded GRACE TWS product is assimilated, spatially distributed GRACE error estimates are accounted for, and the impact that GRACE scaling factors have on assimilation is evaluated. Comparisons with quality-controlled in situ observations indicate that GRACE DA has a positive impact on the simulation of unconfined groundwater variability across the majority of the eastern United States and on the simulation of surface and root zone soil moisture across the country. Smaller improvements are seen in the simulation of snow depth, and the impact of GRACE DA on simulated river discharge and evapotranspiration is regionally variable. The use of GRACE scaling factors during assimilation improved DA results in the western United States but led to small degradations in the eastern United States. The study also found comparable performance between the use of gridded and basin averaged GRACE observations in assimilation. Finally, the evaluations presented in the paper indicate that GRACE DA can be helpful in improving the representation of droughts.

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.

A common set of statistical metrics has been used to summarize the performance of models or measurements-­ the most widely used ones being bias, mean square error, and linear correlation coefficient. They assume linear, additive, Gaussian errors, and they are interdependent, incomplete, and incapable of directly quantifying un­certainty. The authors demonstrate that these metrics can be directly derived from the parameters of the simple linear error model. Since a correct error model captures the full error information, it is argued that the specification of a parametric error model should be an alternative to the metrics-based approach. The error-modeling meth­odology is applicable to both linear and nonlinear errors, while the metrics are only meaningful for linear errors. In addition, the error model expresses the error structure more naturally, and directly quantifies uncertainty. This argument is further explained by highlighting the intrinsic connections between the performance metrics, the error model, and the joint distribution between the data and the reference.

The multiscale modeling framework (MMF), which replaces traditional cloud parameterizations with cloud-resolving models (CRMs) within a host atmospheric general circulation model (GCM), has become a new approach for climate modeling. The embedded CRMs make it possible to apply CRM-based cloud microphysics directly within a GCM. However, most such schemes have never been tested in a global environment for long-term climate simulation. The benefits of using an MMF to evaluate rigorously and improve microphysics schemes are here demonstrated. Four one-moment microphysical schemes are implemented into the Goddard MMF and their results validated against three CloudSat/CALIPSO cloud ice products and other satellite data. The new four-class (cloud ice, snow, graupel, and frozen drops/hail) ice scheme produces a better overall spatial distribution of cloud ice amount, total cloud fractions, net radiation, and total cloud radiative forcing than earlier three-class ice schemes, with biases within the observational uncertainties. Sensitivity experiments are conducted to examine the impact of recently upgraded microphysical processes on global hydrometeor distributions. Five processes dominate the global distributions of cloud ice and snow amount in long-term simulations: (1) allowing for ice supersaturation in the saturation adjustment, (2) three additional correction terms in the depositional growth of cloud ice to snow, (3) accounting for cloud ice fall speeds, (4) limiting cloud ice particle size, and (5) new size-mapping schemes for snow and graupel. Despite the cloud microphysics improvements, systematic errors associated with subgrid processes, cyclic lateral boundaries in the embedded CRMs, and momentum transport remain and will require future improvement.

Advances in satellite monitoring of the terrestrial water cycle have led to a concerted effort to assimilate soil moisture observations from various platforms into offline land surface models (LSMs). One principal but still open question is that of the ability of land data assimilation (LDA) to improve LSM initial conditions for coupled short-term weather prediction. In this study, the impact of assimilating Advanced Microwave Scanning Radiometer for EOS (AMSR-E) soil moisture retrievals on coupled WRF Model forecasts is examined during the summers of dry (2006) and wet (2007) surface conditions in the southern Great Plains. LDA is carried out using NASAs Land Information System (LIS) and the Noah LSM through an ensemble Kalman filter (EnKF) approach. The impacts of LDA on the 1) soil moisture and soil temperature initial conditions for WRF, 2) land-atmosphere coupling characteristics, and 3) ambient weather of the coupled LIS-WRF simulations are then assessed. Results show that impacts of soil moisture LDA during the spin-up can significantly modify LSM states and fluxes, depending on regime and season. Results also indicate that the use of seasonal cumulative distribution functions (CDFs) is more advantageous compared to the traditional annual CDF bias correction strategies. LDA performs consistently regardless of atmospheric forcing applied, with greater improvements seen when using coarser, global forcing products. Downstream impacts on coupled simulations vary according to the strength of the LDA impact at the initialization, where significant modifications to the soil moisture flux- PBL-ambient weather process chain are observed. Overall, this study demonstrates potential for future, higher-resolution soil moisture assimilation applications in weather and climate research.

Powered by the sun, water moves continuously between and through Earths oceanic, atmospheric, and terrestrial reservoirs. It enables life, shapes Earths surface, and responds to and influences climate change. Scientists measure various features of the water cycle using a combination of ground, airborne, and space-based observations, and seek to characterize it at multiple scales with the aid of numerical models. Over time our understanding of the water cycle and ability to quantify it have improved, owing to advances in observational capabilities, the extension of the data record, and increases in computing power and storage. Here we present some of the most recent estimates of global and continental ocean basin scale water cycle stocks and fluxes and provide examples of modern numerical modeling systems and reanalyses.Further, we discuss prospects for predicting water cycle variability at seasonal and longer scales, which is complicated by a changing climate and direct human impacts related to water management and agriculture. Changes to the water cycle will be among the most obvious and important facets of climate change, thus it is crucial that we continue to invest in our ability to monitor it.

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.

Seasonal streamflow forecasts, contingent on climate information, can be utilized to ensure water supply for multiple uses including municipal demands, hydroelectric power generation, and for planning agricultural operations. However, uncertainties in the streamflow forecasts pose significant challenges in their utilization in real-time operations. In this study, we systematically decompose various sources of errors in developing seasonal streamflow forecasts from two Land Surface Models (LSMs) (Noah3.2 and CLM2), which are forced with downscaled and disaggregated climate forecasts. In particular, the study quantifies the relative contributions of the sources of errors from LSMs, climate forecasts, and downscaling/disaggregation techniques in developing seasonal streamflow forecast. For this purpose, three month ahead seasonal precipitation forecasts from the ECHAM4.5 general circulation model (GCM) were statistically downscaled from 2.8deg to 1/8deg spatial resolution using principal component regression (PCR) and then temporally disaggregated from monthly to daily time step using kernel-nearest neighbor (K-NN) approach. For other climatic forcings, excluding precipitation, we considered the North American Land Data Assimilation System version 2 (NLDAS-2) hourly climatology over the years 1979 to 2010. Then the selected LSMs were forced with precipitation forecasts and NLDAS-2 hourly climatology to develop retrospective seasonal streamflow forecasts over a period of 20 years (1991-2010). Finally, the performance of LSMs in forecasting streamflow under different schemes was analyzed to quantify the relative contribution of various sources of errors in developing seasonal streamflow forecast. Our results indicate that the most dominant source of errors during winter and fall seasons is the errors due to ECHAM4.5 precipitation forecasts, while temporal disaggregation scheme contributes to maximum errors during summer season.

Despite recent advances in land surfacemodeling and remote sensing, estimates of the global water budget are still fairly uncertain. This study aims to evaluate the water budget of the Amazon basin based on several state-ofthe- art land surface model (LSM) outputs. Water budget variables (terrestrial water storage TWS, evapotranspiration ET, surface runoff R, and base flow B) are evaluated at the basin scale using both remote sensing and in situ data. Meteorological forcings at a 3-hourly time step and 18 spatial resolution were used to run 14 LSMs. Precipitation datasets that have been rescaled to matchmonthly Global Precipitation Climatology Project (GPCP) andGlobal Precipitation Climatology Centre (GPCC) datasets and the daily Hydrologie du Bassin de l'Amazone (HYBAM) dataset were used to perform three experiments. The Hydrological Modeling and Analysis Platform (HyMAP) river routing scheme was forced with R and B and simulated discharges are compared against observations at 165 gauges. Simulated ET and TWS are compared against FLUXNET and MOD16A2 evapotranspiration datasets andGravity Recovery and ClimateExperiment (GRACE)TWSestimates in two subcatchments of main tributaries (Madeira and Negro Rivers).At the basin scale, simulated ET ranges from 2.39 to 3.26 mm day(exp -1) and a low spatial correlation between ET and precipitation indicates that evapotranspiration does not depend on water availability over most of the basin. Results also show that other simulated water budget components vary significantly as a function of both the LSM and precipitation dataset, but simulated TWS generally agrees with GRACE estimates at the basin scale. The best water budget simulations resulted from experiments using HYBAM, mostly explained by a denser rainfall gauge network and the rescaling at a finer temporal scale.

Land surface microwave emissivity affects remote sensing of both the atmosphere and the land surface. The dynamical behavior of microwave emissivity over a very diverse sample of land surface types is studied. With seven years of satellite measurements from AMSR-E, we identified various dynamical regimes of the land surface emission. In addition, we used two radiative transfer models (RTMs), the Community Radiative Transfer Model (CRTM) and the Community Microwave Emission Modeling Platform (CMEM), to simulate land surface emissivity dynamics. With both CRTM and CMEM coupled to NASA's Land Information System, global-scale land surface microwave emissivities were simulated for five years, and evaluated against AMSR-E observations. It is found that both models have successes and failures over various types of land surfaces. Among them, the desert shows the most consistent underestimates (by approx. 70-80%), due to limitations of the physical models used, and requires a revision in both systems. Other snow-free surface types exhibit various degrees of success and it is expected that parameter tuning can improve their performances.

Uncertainties in the retrievals of microwaveland-surface emissivities are quantified over two types of land surfaces: desert and tropical rainforest. Retrievals from satellite-based microwave imagers, including the Special Sensor Microwave Imager, the Tropical Rainfall Measuring Mission Microwave Imager, and the Advanced Microwave Scanning Radiometer for Earth Observing System, are studied. Our results show that there are considerable differences between the retrievals from different sensors and from different groups over these two land-surface types. In addition, the mean emissivity values show different spectral behavior across the frequencies. With the true emissivity assumed largely constant over both of the two sites throughout the study period, the differences are largely attributed to the systematic and random errors inthe retrievals. Generally, these retrievals tend to agree better at lower frequencies than at higher ones, with systematic differences ranging 1%-4% (3-12 K) over desert and 1%-7% (3-20 K) over rainforest. The random errors within each retrieval dataset are in the range of 0.5%-2% (2-6 K). In particular, at 85.5/89.0 GHz, there are very large differences between the different retrieval datasets, and within each retrieval dataset itself. Further investigation reveals that these differences are most likely caused by rain/cloud contamination, which can lead to random errors up to 10-17 K under the most severe conditions.

The present generation of general circulation models (GCM) use parameterized cumulus schemes and run at hydrostatic grid resolutions. To improve the representation of cloud-scale moist processes and landeatmosphere interactions, a global, Multi-scale Modeling Framework (MMF) coupled to the Land Information System (LIS) has been developed at NASA-Goddard Space Flight Center. The MMFeLIS has three components, a finite-volume (fv) GCM (Goddard Earth Observing System Ver. 4, GEOS-4), a 2D cloud-resolving model (Goddard Cumulus Ensemble, GCE), and the LIS, representing the large-scale atmospheric circulation, cloud processes, and land surface processes, respectively. The non-hydrostatic GCE model replaces the single-column cumulus parameterization of fvGCM. The model grid is composed of an array of fvGCM gridcells each with a series of embedded GCE models. A horizontal coupling strategy, GCE4fvGCM4Coupler4LIS, offered significant computational efficiency, with the scalability and I/O capabilities of LIS permitting landeatmosphere interactions at cloud-scale. Global simulations of 2007e2008 and comparisons to observations and reanalysis products were conducted. Using two different versions of the same land surface model but the same initial conditions, divergence in regional, synoptic-scale surface pressure patterns emerged within two weeks. The sensitivity of largescale circulations to land surface model physics revealed significant functional value to using a scalable, multi-model land surface modeling system in global weather and climate prediction.

Land-atmosphere (L-A) interactions play a critical role in determining the diurnal evolution of both planetary boundary layer (PBL) and land surface temperature and moisture budgets, as well as controlling feedbacks with clouds and precipitation that lead to the persistence of dry and wet regimes. Recent efforts to quantify the strength of L-A coupling in prediction models have produced diagnostics that integrate across both the land and PBL components of the system. In this study, we examine the impact of improved specification of land surface states, anomalies, and fluxes on coupled WRF forecasts during the summers of extreme dry (2006) and wet (2007) land surface conditions in the U.S. Southern Great Plains. The improved land initialization and surface flux parameterizations are obtained through the use of a new optimization and uncertainty estimation module in NASA's Land Information System (LIS-OPT/UE), whereby parameter sets are calibrated in the Noah land surface model and classified according to a land cover and soil type mapping of the observation sites to the full model domain. The impact of calibrated parameters on the a) spinup of the land surface used as initial conditions, and b) heat and moisture states and fluxes of the coupled WRF simulations are then assessed in terms of ambient weather and land-atmosphere coupling along with measures of uncertainty propagation into the forecasts. In addition, the sensitivity of this approach to the period of calibration (dry, wet, average) is investigated. Finally, tradeoffs of computational tractability and scientific validity, and the potential for combining this approach with satellite remote sensing data are also discussed.

Land-atmosphere (L-A) Interactions playa critical role in determining the diurnal evolution of both planetary boundary layer (PBL) and land surface heat and moisture budgets, as well as controlling feedbacks with clouds and precipitation that lead to the persistence of dry and wet regimes. Recent efforts to quantify the strength of L-A coupling in prediction models have produced diagnostics that integrate across both the land and PBL components of the system. In this study, we examine the impact of improved specification of land surface states, anomalies, and fluxes on coupled WRF forecasts during the summers of extreme dry (2006) and wet (2007) land surface conditions in the U.S. Southern Great Plains. The improved land initialization and surface flux parameterizations are obtained through the use of a new optimization and uncertainty estimation module in NASA's Land Information System (US-OPT/UE), whereby parameter sets are calibrated in the Noah land surface model and classified according to a land cover and soil type mapping of the observation sites to the full model domain. The impact of calibrated parameters on the a) spinup of the land surface used as initial conditions, and b) heat and moisture states and fluxes of the coupled WRF Simulations are then assessed in terms of ambient weather and land-atmosphere coupling along with measures of uncertainty propagation into the forecasts. In addition, the sensitivity of this approach to the period of calibration (dry, wet, average) is investigated. Finally, tradeoffs of computational tractability and scientific validity, and the potential for combining this approach with satellite remote sensing data are also discussed.

Since June 2010, the NASA Short-term Prediction Research and Transition (SPoRT; Goodman et al. 2004; Darden et al. 2010; Stano et al. 2012; Fuell et al. 2012) Center has been generating a real-time Normalized Difference Vegetation Index (NDVI) and corresponding Green Vegetation Fraction (GVF) composite based on reflectances from NASA s Moderate Resolution Imaging Spectroradiometer (MODIS) instrument. This dataset is generated at 0.01 resolution across the Continental United States (CONUS), and updated daily. The goal of producing such a vegetation dataset is to improve over the default climatological GVF dataset in land surface and numerical weather prediction models, in order to have better simulations of heat and moisture exchange between the land surface and the planetary boundary layer. Details on the SPoRT/MODIS vegetation composite algorithm are presented in Case et al. (2011). Vegetation indices such as GVF and Leaf Area Index (LAI) are used by land surface models (LSMs) to represent the horizontal and vertical density of plant vegetation (Gutman and Ignatov 1998), in order to calculate transpiration, interception and radiative shading. Both of these indices are related to the NDVI; however, there is an inherent ambiguity in determining GVF and LAI simultaneously from NDVI, as described in Gutman and Ignatov (1998). One practice is to specify the LAI while allowing the GVF to vary both spatially and temporally, as is done in the Noah LSM (Chen and Dudhia 2001; Ek et al. 2003). Operational versions of Noah within several of the National Centers for Environmental Prediction (NCEP) global and regional modeling systems hold the LAI fixed, while the GVF varies according to a global monthly climatology. This GVF climatology was derived from NDVI data on the NOAA Advanced Very High Resolution Radiometer (AVHRR) polar orbiting satellite, using information from 1985 to 1991 (Gutman and Ignatov 1998; Jiang et al. 2010). Representing data at the mid-point of every month, the climatological dataset is on a grid with 0.144 (~16 km) spatial resolution and is distributed with the community WRF model (Ek et al. 2003; Jiang et al. 2010; Skamarock et al. 2008).

Land-atmosphere (L-A) interactions play a critical role in determining the diurnal evolution of both planetary boundary layer (PBL) and land surface heat and moisture budgets, as well as controlling feedbacks with clouds and precipitation that lead to the persistence of dry and wet regimes. Recent efforts to quantify the strength of L-A coupling in prediction models have produced diagnostics that integrate across both the land and PBL components of the system. In this study, we examine the impact of improved specification of land surface states, anomalies, and fluxes on coupled WRF forecasts during the summers of extreme dry and wet land surface conditions in the U.S. Southern Great Plains. The improved land initialization and surface flux parameterizations are obtained through calibration of the Noah land surface model using the new optimization and uncertainty estimation subsystem in NASA's Land Information System (LIS-OPT/UE). The impact of the calibration on the a) spinup of the land surface used as initial conditions, and b) the simulated heat and moisture states and fluxes of the coupled WRF simulations is then assessed. Changes in ambient weather and land-atmosphere coupling are evaluated along with measures of uncertainty propagation into the forecasts. In addition, the sensitivity of this approach to the period of calibration (dry, wet, average) is investigated. Results indicate that the offline calibration leads to systematic improvements in land-PBL fluxes and near-surface temperature and humidity, and in the process provide guidance on the questions of what, how, and when to calibrate land surface models for coupled model prediction.

Land-atmosphere (L-A) interactions play a critical role in determining the diurnal evolution of land surface and planetary boundary layer (PBL) temperature and moisture states and fluxes. In turn, these interactions regulate the strength of the connection between surface moisture and precipitation in a coupled system. To address model deficiencies, recent studies have focused on development of diagnostics to quantify the strength and accuracy of the land-PBL coupling at the process level. In this paper, a diagnosis of the nature and impacts of local land-atmosphere coupling (LoCo) during dry and wet extreme conditions is presented using a combination of models and observations during the summers of 2006 and 2007 in the U.S. southern Great Plains. A range of diagnostics exploring the links and feedbacks between soil moisture and precipitation is applied to the dry/wet regimes exhibited in this region, and in the process, a thorough evaluation of nine different land-PBL scheme couplings is conducted under the umbrella of a high-resolution regional modeling test bed. Results show that the sign and magnitude of errors in land surface energy balance components are sensitive to the choice of land surface model, regime type, and running mode. In addition, LoCo diagnostics show that the sensitivity of L-A coupling is stronger toward the land during dry conditions, while the PBL scheme coupling becomes more important during the wet regime. Results also demonstrate how LoCo diagnostics can be applied to any modeling system (e.g., reanalysis products) in the context of their integrated impacts on the process chain connecting the land surface to the PBL and in support of hydrological anomalies.

An assessment ofETestimates for current LDAS systems is provided along with current research that demonstrates improvement in LSM ET estimates due to assimilating satellite-based soil moisture products. Using the Ensemble Kalman Filter in the Land Information System, we assimilate both NASA and Land Parameter Retrieval Model (LPRM) soil moisture products into the Noah LSM Version 3.2 with the North American LDAS phase 2 CNLDAS-2) forcing to mimic the NLDAS-2 configuration. Through comparisons with two global reference ET products, one based on interpolated flux tower data and one from a new satellite ET algorithm, over the NLDAS2 domain, we demonstrate improvement in ET estimates only when assimilating the LPRM soil moisture product.

Land-atmosphere (L-A) interactions play a critical role in determining the diurnal evolution of land surface and planetary boundary layer (PBL) temperature and moisture states and fluxes. In turn, these interactions regulate the strength of the connection between surface moisture and precipitation in a coupled system. To address model deficiencies, recent studies have focused on development of diagnostics to quantify the strength and accuracy of the land- PBL coupling at the process-level. In this paper, a diagnosis of the nature and impacts of local land-atmosphere coupling (LoCo) during dry and wet extreme conditions is presented using a combination of models and observations during the summers of 2006 and 2007 in the U.S. Southern Great Plains. A range of diagnostics exploring the links and feedbacks between soil moisture and precipitation are applied to the dry/wet regimes exhibited in this region, and in the process a thorough evaluation of nine different land-PBL scheme couplings is conducted under the umbrella of a high-resolution regional modeling testbed. Results show that the sign and magnitude of errors in land surface energy balance components are sensitive to the choice of land surface model, regime type, and running mode. In addition, LoCo diagnostics show that the sensitivity of L-A coupling is stronger towards the land during dry conditions, while the PBL scheme coupling becomes more important during the wet regime. Results also demonstrate how LoCo diagnostics can be applied to any modeling system (e.g. reanalysis products) in the context of their integrated impacts on the process-chain connecting the land surface to the PBL and in support of hydrological anomalies.

The NASA developed Land Information System (LIS) is the Air Force Weather Agency's (AFWA) operational Land Data Assimilation System (LDAS) combining real time precipitation observations and analyses, global forecast model data, vegetation, terrain, and soil parameters with the community Noah land surface model, along with other hydrology module options, to generate profile analyses of global soil moisture, soil temperature, and other important land surface characteristics. (1) A range of satellite data products and surface observations used to generate the land analysis products (2) Global, 1/4 deg spatial resolution (3) Model analysis generated at 3 hours. AFWA recognizes the importance of operational benchmarking and uncertainty characterization for land surface modeling and is developing standard methods, software, and metrics to verify and/or validate LIS output products. To facilitate this and other needs for land analysis activities at AFWA, the Model Evaluation Toolkit (MET) -- a joint product of the National Center for Atmospheric Research Developmental Testbed Center (NCAR DTC), AFWA, and the user community -- and the Land surface Verification Toolkit (LVT), developed at the Goddard Space Flight Center (GSFC), have been adapted to operational benchmarking needs of AFWA's land characterization activities.

The NASA Short-term Prediction Research and Transition (SPoRT) Center has developed a Greenness Vegetation Fraction (GVF) dataset, which is updated daily using swaths of Normalized Difference Vegetation Index data from the Moderate Resolution Imaging Spectroradiometer (MODIS) data aboard the NASA-EOS Aqua and Terra satellites. NASA SPoRT started generating daily real-time GVF composites at 1-km resolution over the Continental United States beginning 1 June 2010. A companion poster presentation (Bell et al.) primarily focuses on impact results in an offline configuration of the Noah land surface model (LSM) for the 2010 warm season, comparing the SPoRT/MODIS GVF dataset to the current operational monthly climatology GVF available within the National Centers for Environmental Prediction (NCEP) and Weather Research and Forecasting (WRF) models. This paper/presentation primarily focuses on individual case studies of severe weather events to determine the impacts and possible improvements by using the real-time, high-resolution SPoRT-MODIS GVFs in place of the coarser-resolution NCEP climatological GVFs in model simulations. The NASA-Unified WRF (NU-WRF) modeling system is employed to conduct the sensitivity simulations of individual events. The NU-WRF is an integrated modeling system based on the Advanced Research WRF dynamical core that is designed to represents aerosol, cloud, precipitation, and land processes at satellite-resolved scales in a coupled simulation environment. For this experiment, the coupling between the NASA Land Information System (LIS) and the WRF model is utilized to measure the impacts of the daily SPoRT/MODIS versus the monthly NCEP climatology GVFs. First, a spin-up run of the LIS is integrated for two years using the Noah LSM to ensure that the land surface fields reach an equilibrium state on the 4-km grid mesh used. Next, the spin-up LIS is run in two separate modes beginning on 1 June 2010, one continuing with the climatology GVFs while the other uses the daily SPoRT/MODIS GVFs. Finally, snapshots of the LIS land surface fields are used to initialize two different simulations of the NU-WRF, one running with climatology LIS and GVFs, and the other running with experimental LIS and NASA/SPoRT GVFs. In this paper/presentation, case study results will be highlighted in regions with significant differences in GVF between the NCEP climatology and SPoRT product during severe weather episodes.

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.

Both NLDAS Phase 1 (1996-2007) and Phase 2 (1979-present) datasets have been evaluated against in situ observational datasets, and NLDAS forcings and outputs are used by a wide variety of users. Drought indices and drought monitoring from NLDAS were recently examined by Mo et al. (2010) and Sheffield et al. (2010). In this poster, we will present results analyzing NLDAS Phase 2 forcings and outputs for 3 North American Case studies being analyzed as part of the NOAA MAPP Drought Task Force: (1) Western US drought (1998- 2004); (2) plains/southeast US drought (2006-2007); and (3) Current Texas-Mexico drought (2011-). We will examine percentiles of soil moisture consistent with the NLDAS drought monitor.

Uncertainties in the retrievals of microwave land surface emissivities were quantified over two types of land surfaces: desert and tropical rainforest. Retrievals from satellite-based microwave imagers, including SSM/I, TMI and AMSR-E, were studied. Our results show that there are considerable differences between the retrievals from different sensors and from different groups over these two land surface types. In addition, the mean emissivity values show different spectral behavior across the frequencies. With the true emissivity assumed largely constant over both of the two sites throughout the study period, the differences are largely attributed to the systematic and random errors in the retrievals. Generally these retrievals tend to agree better at lower frequencies than at higher ones, with systematic differences ranging 1~4% (3~12 K) over desert and 1~7% (3~20 K) over rainforest. The random errors within each retrieval dataset are in the range of 0.5~2% (2~6 K). In particular, at 85.0/89.0 GHz, there are very large differences between the different retrieval datasets, and within each retrieval dataset itself. Further investigation reveals that these differences are mostly likely caused by rain/cloud contamination, which can lead to random errors up to 10~17 K under the most severe conditions.

The NASA developed Land Information System (LIS) is the Air Force Weather Agency's (AFWA) operational Land Data Assimilation System (LDAS) combining real time precipitation observations and analyses, global forecast model data, vegetation, terrain, and soil parameters with the community Noah land surface model, along with other hydrology module options, to generate profile analyses of global soil moisture, soil temperature, and other important land surface characteristics. (1) A range of satellite data products and surface observations used to generate the land analysis products (2) Global, 1/4 deg spatial resolution (3) Model analysis generated at 3 hours

Land-atmosphere (L-A) interactions play a critical role in determining the diurnal evolution of land surface and planetary boundary layer (PBL) temperature and moisture states and fluxes. In turn, these interactions regulate the strength of the connection between surface moisture and precipitation in a coupled system. To address deficiencies in numerical weather prediction and climate models due to improper treatment of L-A interactions, recent studies have focused on development of diagnostics to quantify the strength and accuracy of the land-PBL coupling at the process-level. In this study, a diagnosis of the nature and impacts of local land-atmosphere coupling (LoCo) during dry and wet extreme conditions is presented using a combination of models and observations during the summers of2006-7 in the U.S. Southern Great Plains. Specifically, the Weather Research and Forecasting (WRF) model has been coupled to NASA's Land Information System (LIS), which provides a flexible and high resolution representation and initialization of land surface physics and states. A range of diagnostics exploring the links and feedbacks between soil moisture and precipitation are examined for the dry/wet regimes of this region, along with the behavior and accuracy of different land-PBL scheme couplings under these conditions. Results demonstrate how LoCo diagnostics can be applied to coupled model components in the context of their integrated impacts on the process-chain connecting the land surface to the PBL and support of hydrological anomalies.

Model evaluation and verification are key in improving the usage and applicability of simulation models for real-world applications. In this article, the development and capabilities of a formal system for land surface model evaluation called the Land surface Verification Toolkit (LVT) is described. LVT is designed to provide an integrated environment for systematic land model evaluation and facilitates a range of verification approaches and analysis capabilities. LVT operates across multiple temporal and spatial scales and employs a large suite of in-situ, remotely sensed and other model and reanalysis datasets in their native formats. In addition to the traditional accuracy-based measures, LVT also includes uncertainty and ensemble diagnostics, information theory measures, spatial similarity metrics and scale decomposition techniques that provide novel ways for performing diagnostic model evaluations. Though LVT was originally designed to support the land surface modeling and data assimilation framework known as the Land Information System (LIS), it also supports hydrological data products from other, non-LIS environments. In addition, the analysis of diagnostics from various computational subsystems of LIS including data assimilation, optimization and uncertainty estimation are supported within LVT. Together, LIS and LVT provide a robust end-to-end environment for enabling the concepts of model data fusion for hydrological applications. The evolving capabilities of LVT framework are expected to facilitate rapid model evaluation efforts and aid the definition and refinement of formal evaluation procedures for the land surface modeling community.

The inherent coupled nature of earth s energy and water cycles places significant importance on the proper representation and diagnosis of land atmosphere (LA) interactions in hydrometeorological prediction models. However, the precise nature of the soil moisture precipitation relationship at the local scale is largely determined by a series of nonlinear processes and feedbacks that are difficult to quantify. To quantify the strength of the local LA coupling (LoCo), this process chain must be considered both in full and as individual components through their relationships and sensitivities. To address this, recent modeling and diagnostic studies have been extended to 1) quantify the processes governing LoCo utilizing the thermodynamic properties of mixing diagrams, and 2) diagnose the sensitivity of coupled systems, including clouds and moist processes, to perturbations in soil moisture. This work employs NASA s Land Information System (LIS) coupled to the Weather Research and Forecasting (WRF) mesoscale model and simulations performed over the U.S. Southern Great Plains. The behavior of different planetary boundary layers (PBL) and land surface scheme couplings in LIS WRF are examined in the context of the evolution of thermodynamic quantities that link the surface soil moisture condition to the PBL regime, clouds, and precipitation. Specifically, the tendency toward saturation in the PBL is quantified by the lifting condensation level (LCL) deficit and addressed as a function of time and space. The sensitivity of the LCL deficit to the soil moisture condition is indicative of the strength of LoCo, where both positive and negative feedbacks can be identified. Overall, this methodology can be applied to any model or observations and is a crucial step toward improved evaluation and quantification of LoCo within models, particularly given the advent of next-generation satellite measurements of PBL and land surface properties along with advances in data assimilation schemes.

Data assimilation is being increasingly used to merge remotely sensed land surface variables such as soil moisture, snow and skin temperature with estimates from land models. Its success, however, depends on unbiased model predictions and unbiased observations. Here, a suite of continental-scale, synthetic soil moisture assimilation experiments is used to compare two approaches that address typical biases in soil moisture prior to data assimilation: (i) parameter estimation to calibrate the land model to the climatology of the soil moisture observations, and (ii) scaling of the observations to the model s soil moisture climatology. To enable this research, an optimization infrastructure was added to the NASA Land Information System (LIS) that includes gradient-based optimization methods and global, heuristic search algorithms. The land model calibration eliminates the bias but does not necessarily result in more realistic model parameters. Nevertheless, the experiments confirm that model calibration yields assimilation estimates of surface and root zone soil moisture that are as skillful as those obtained through scaling of the observations to the model s climatology. Analysis of innovation diagnostics underlines the importance of addressing bias in soil moisture assimilation and confirms that both approaches adequately address the issue.

Land-Atmosphere coupling is typically designed and implemented independently for physical (e.g. water and energy) and chemical (e.g. biogenic emissions and surface depositions)-based models and applications. Differences in scale, data requirements, and physics thus limit the ability of Earth System models to be fully coupled in a consistent manner. In order for the physical-chemical-biological coupling to be complete, treatment of the land in terms of surface classification, condition, fluxes, and emissions must be considered simultaneously and coherently across all components. In this study, we investigate a coupling strategy for the NASA-Unified Weather Research and Forecasting (NU-WRF) model that incorporates the traditionally disparate fluxes of water and energy through NASA's LIS (Land Information System) and biogenic emissions through BEIS (Biogenic Emissions Inventory System) and MEGAN (Model of Emissions of Gases and Aerosols from Nature) into the atmosphere. In doing so, inconsistencies across model inputs and parameter data are resolved such that the emissions from a particular plant species are consistent with the heat and moisture fluxes calculated for that land cover type. In turn, the response of the atmospheric turbulence and mixing in the planetary boundary layer (PBL) acts on the identical surface type, fluxes, and emissions for each. In addition, the coupling of dust emission within the NU-WRF system is performed in order to ensure consistency and to maximize the benefit of high-resolution land representation in LIS. The impacts of those self-consistent components on' the simulation of atmospheric aerosols are then evaluated through the WRF-Chem-GOCART (Goddard Chemistry Aerosol Radiation and Transport) model. Overall, this ambitious project highlights the current difficulties and future potential of fully coupled. components. in Earth System models, and underscores the importance of the iLEAPS community in supporting improved knowledge of processes and innovative approaches for models and observations.

Land-atmosphere interactions play a critical role in determining the. diurnal evolution of both planetary boundary layer (PBL) and land surface temperature and moisture states. The degree of coupling between the land surface and PBL in numerical weather prediction and climate models remains largely unexplored and undiagnosed due to the complex interactions and feedbacks present across a range of scales. Further, uncoupled systems or experiments (e.g., the Project for Intercomparison of Land Parameterization Schemes, PILPS) may lead to inaccurate water and energy cycle process understanding by neglecting feedback processes such as PBL-top entrainment. In this study, a framework for diagnosing local land-atmosphere coupling (LoCo) is presented using a coupled mesoscale model with a suite of PBL and land surface model (LSM) options along with observations during the summers of 200617 in the U.S. Southern Great Plains. Specifically, the Weather Research and Forecasting (WRF) model has been coupled to NASA's Land Information System (LIS), which provides a flexible and high-resolution representation and initialization of land surface physics and states. A range of diagnostics exploring the links and feedbacks between soil moisture and precipitation are examined for the dry/wet extremes of this region, along with the sensitivity of PBL-LSM coupling to perturbations in soil moisture. As such, this methodology provides a potential pathway to study factors controlling local land-atmosphere coupling (LoCo) using the LIS-WRF system, which is serving as a testbed for LoCo experiments to evaluate coupling diagnostics within the community.

The degree of coupling between the land surface and PBL in NWP models remains largely undiagnosed due to the complex interactions and feedbacks present across a range of scales. In this study, a framework for diagnosing local land-atmosphere coupling (LoCo) is presented using a coupled mesoscale model with observations during the summers of 2006/7 in the U.S. Southern Great Plains. Specifically, the Weather Research and Forecasting (WRF) model has been coupled to NASA's Land Information System (LIS), which enables a suite of PBL and land surface model (LSM) options along provides a flexible and high-resolution representation and initialization of land surface physics and states. This coupling is one component of a larger project to develop a NASA-Unified WRF (NU-WRF) system. A range of diagnostics exploring the feedbacks between soil moisture and precipitation are examined for the dry/wet extremes, along with the sensitivity of PBL-LSM coupling to perturbations in soil moisture.

Land-atmosphere interactions play a critical role in determining the diurnal evolution of both planetary boundary layer (PBL) and land surface temperature and moisture states. The degree of coupling between the land surface and PBL in numerical weather prediction and climate models remains largely unexplored and undiagnosed due to the complex interactions and feedbacks present across a range of scales. Further, uncoupled systems or experiments (e.g., the Project for Intercomparison of Land Parameterization Schemes, PILPS) may lead to inaccurate water and energy cycle process understanding by neglecting feedback processes such as PBL-top entrainment. In this study, a framework for diagnosing local land-atmosphere coupling (LoCo) is presented using a coupled mesoscale model with a suite of PBL and land surface model (LSM) options along with observations during the summers of200617 in the U.S. Southern Great Plains. Specifically, the Weather Research and Forecasting (WRF) model has been coupled to NASA's Land Information System (LIS), which provides a flexible and high-resolution representation and initialization of land surface physics and states. A range of diagnostics exploring the links and feedbacks between soil moisture and precipitation are examined for the dry/wet extremes of this region, along with the sensitivity of PBL-LSM coupling to perturbations in soil moisture. As such, this methodology provides a potential pathway to study factors controlling local land-atmosphere coupling (LoCo) using the LIS-WRF system, which is serving as a testbed for LoCo experiments to evaluate coupling diagnostics within the community.

The Land Information System (LIS; http://lis.gsfc.nasa.gov) is a flexible land surface modeling framework that has been developed with the goal of integrating satellite-and ground-based observational data products and advanced land surface modeling techniques to produce optimal fields of land surface states and fluxes. As such, LIS represents a step towards the next generation land component of an integrated Earth system model. In recognition of LIS object-oriented software design, use and impact in the land surface and hydrometeorological modeling community, the LIS software was selected as a co-winner of NASA?s 2005 Software of the Year award.LIS facilitates the integration of observations from Earth-observing systems and predictions and forecasts from Earth System and Earth science models into the decision-making processes of partnering agency and national organizations. Due to its flexible software design, LIS can serve both as a Problem Solving Environment (PSE) for hydrologic research to enable accurate global water and energy cycle predictions, and as a Decision Support System (DSS) to generate useful information for application areas including disaster management, water resources management, agricultural management, numerical weather prediction, air quality and military mobility assessment. LIS has e volved from two earlier efforts -- North American Land Data Assimilation System (NLDAS) and Global Land Data Assimilation System (GLDAS) that focused primarily on improving numerical weather prediction skills by improving the characterization of the land surface conditions. Both of GLDAS and NLDAS now use specific configurations of the LIS software in their current implementations.In addition, LIS was recently transitioned into operations at the US Air Force Weather Agency (AFWA) to ultimately replace their Agricultural Meteorology (AGRMET) system, and is also used routinely by NOAA's National Centers for Environmental Prediction (NCEP)/Environmental Modeling Center (EMC) for their land data assimilation systems to support weather and climate modeling. LIS not only consolidates the capabilities of these two systems, but also enables a much larger variety of configurations with respect to horizontal spatial resolution, input datasets and choice of land surface model through "plugins". LIS has been coupled to the Weather Research and Forecasting (WRF) model to support studies of land-atmosphere coupling be enabling ensembles of land surface states to be tested against multiple representations of the atmospheric boundary layer. LIS has also been demonstrated for parameter estimation, who showed that the use of sequential remotely sensed soil moisture products can be used to derive soil hydraulic and texture properties given a sufficient dynamic range in the soil moisture retrievals and accurate precipitation inputs.LIS has also recently been demonstrated for multi-model data assimilation using an Ensemble Kalman Filter for sequential assimilation of soil moisture, snow, and temperature.Ongoing work has demonstrated the value of bias correction as part of the filter, and also that of joint calibration and assimilation.Examples and case studies demonstrating the capabilities and impacts of LIS for hydrometeorological modeling, assimilation and parameter estimation will be presented as advancements towards the next generation of integrated observation and modeling systems

Determining uncertainties in satellite-based multi-sensor quantitative precipitation estimates over land of fundamental importance to both data producers and hydro climatological applications. ,Evaluating TRMM-era products also lays the groundwork and sets the direction for algorithm and applications development for future missions including GPM. QPE uncertainties result mostly from the interplay of systematic errors and random errors. In this work, we will synthesize our recent results quantifying the error characteristics of satellite-based precipitation estimates. Both systematic errors and total uncertainties have been analyzed for six different TRMM-era precipitation products (3B42, 3B42RT, CMORPH, PERSIANN, NRL and GSMap). For systematic errors, we devised an error decomposition scheme to separate errors in precipitation estimates into three independent components, hit biases, missed precipitation and false precipitation. This decomposition scheme reveals hydroclimatologically-relevant error features and provides a better link to the error sources than conventional analysis, because in the latter these error components tend to cancel one another when aggregated or averaged in space or time. For the random errors, we calculated the measurement spread from the ensemble of these six quasi-independent products, and thus produced a global map of measurement uncertainties. The map yields a global view of the error characteristics and their regional and seasonal variations, reveals many undocumented error features over areas with no validation data available, and provides better guidance to global assimilation of satellite-based precipitation data. Insights gained from these results and how they could help with GPM will be highlighted.

Objective: The objective of this proposal is to provide a routine land surface modeling and data assimilation capability for GPM in order to provide global land surface states that are necessary to support physical precipitation retrieval algorithms over land. It is well-known that surface emission, particularly over the range of frequencies to be included in GPM, is sensitive to land surface states, including soil properties, vegetation type and greenness, soil moisture, surface temperature, and snow cover, density, and grain size. Therefore, providing a robust capability to routinely provide these critical land states is essential to support GPM-era physical retrieval algorithms over land.

The Land Information System (LIS; http://lis.gsfc.nasa.gov; Kumar et al., 2006; Peters- Lidard et al.,2007) is a flexible land surface modeling framework that has been developed with the goal of integrating satellite- and ground-based observational data products and advanced land surface modeling techniques to produce optimal fields of land surface states and fluxes. As such, LIS represents a step towards the next generation land component of an integrated Earth system model. In recognition of LIS object-oriented software design, use and impact in the land surface and hydrometeorological modeling community, the LIS software was selected ase co-winner of NASA's 2005 Software of the Year award. LIS facilitates the integration of observations from Earth-observing systems and predictions and forecasts from Earth System and Earth science models into the decision-making processes of partnering agency and national organizations. Due to its flexible software design, LIS can serve both as a Problem Solving Environment (PSE) for hydrologic research to enable accurate global water and energy cycle predictions, and as a Decision Support System (DSS) to generate useful information for application areas including disaster management, water resources management, agricultural management, numerical weather prediction, air quality and military mobility assessment. LIS has evolved from two earlier efforts North American Land Data Assimilation System (NLDAS; Mitchell et al. 2004) and Global Land Data Assimilation System (GLDAS; Rodell al. 2004) that focused primarily on improving numerical weather prediction skills by improving the characterization of the land surface conditions. Both of GLDAS and NLDAS now use specific configurations of the LIS software in their current implementations. In addition, LIS was recently transitioned into operations at the US Air Force Weather Agency (AFWA) to ultimately replace their Agricultural Meteorology (AGRMET) system, and is also used routinely by NOAA's National Centers for Environmental Prediction (NCEP)/Environmental Modeling Center (EMC) for their land data assimilation systems to support weather and climate modeling. LIS not only consolidates the capabilities of these two systems, but also enables a much larger variety of configurations with respect to horizontal spatial resolution, input datasets and choice of land surface model through "plugins,". As described in Kumar et al., 2007, and demonstrated in Case et al., 2008, and Santanello et al., 2009, LIS has been coupled to the Weather Research and Forecasting (WRF) model to support studies of land-atmosphere coupling the enabling ensembles of land surface states to be tested against multiple representations of the atmospheric boundary layer. LIS has also been demonstrated for parameter estimation as described in Peters-Lidard et al. (2008) and Santanello et al. (2007), who showed that the use of sequential remotely sensed soil moisture products can be used to derive soil hydraulic and texture properties given a sufficient dynamic range in the soil moisture retrievals and accurate precipitation inputs. LIS has also recently been demonstrated for multi-model data assimilation (Kumar et al., 2008) using an Ensemble Kalman Filter for sequential assimilation of soil moisture, snow, and temperature. Ongoing work has demonstrated the value of bias correction as part of the filter, and also that of joint calibration and assimilation. Examples and case studies demonstrating the capabilities and impacts of LIS for hydrometeoroogical modeling, assimilation and parameter estimation will be presented as advancements towards the next generation of integrated observation and modeling systems.

Land-atmosphere interactions play a critical role in determining the diurnal evolution of both planetary boundary layer (PBL) and land surface temperature and moisture states. The degree of coupling between the land surface and PBL in numerical weather prediction and climate models remains largely unexplored and undiagnosed due to the complex interactions and feedbacks present across a range of scales. Further, uncoupled systems or experiments (e.g., the Project for Intercomparison of Land Parameterization Schemes, PILPS) may lead to inaccurate water and energy cycle process understanding by neglecting feedback processes such as PBL-top entrainment. In this study, a framework for diagnosing local land-atmosphere coupling is presented using a coupled mesoscale model with a suite of PBL and land surface model (LSM) options along with observations during field experiments in the U. S. Southern Great Plains. Specifically, the Weather Research and Forecasting (WRF) model has been coupled to the Land Information System (LIS), which provides a flexible and high-resolution representation and initialization of land surface physics and states. Within this framework, the coupling established by each pairing of the available PBL schemes in WRF with the LSMs in LIS is evaluated in terms of the diurnal temperature and humidity evolution in the mixed layer. The co-evolution of these variables and the convective PBL is sensitive to and, in fact, integrative of the dominant processes that govern the PBL budget, which are synthesized through the use of mixing diagrams. Results show how the sensitivity of land-atmosphere interactions to the specific choice of PBL scheme and LSM varies across surface moisture regimes and can be quantified and evaluated against observations. As such, this methodology provides a potential pathway to study factors controlling local land-atmosphere coupling (LoCo) using the LIS-WRF system, which will serve as a testbed for future experiments to evaluate coupling diagnostics within the community.