Search

Many regions in Africa and the Middle East are vulnerable to drought and to water and food insecurity, motivating agency efforts such as the U.S. Agency for International Development’s (USAID) Famine Early Warning Systems Network (FEWS NET) to provide early warning of drought events in the region. Each year these warnings guide life-saving assistance that reaches millions of people. A new NASA multimodel, remote sensing–based hydrological forecasting and analysis system, NHyFAS, has been developed to support such efforts by improving the FEWS NET’s current early warning capabilities. NHyFAS derives its skill from two sources: (i) accurate initial conditions, as produced by an offline land modeling system through the application and/or assimilation of various satellite data (precipitation, soil moisture, and terrestrial water storage), and (ii) meteorological forcing data during the forecast period as produced by a state-of-the-art ocean–land–atmosphere forecast system. The land modeling framework used is the Land Information System (LIS), which employs a suite of land surface models, allowing multimodel ensembles and multiple data assimilation strategies to better estimate land surface conditions. An evaluation of NHyFAS shows that its 1–5-month hindcasts successfully capture known historic drought events, and it has improved skill over benchmark-type hindcasts. The system also benefits from strong collaboration with end-user partners in Africa and the Middle East, who provide insights on strategies to formulate and communicate early warning indicators to water and food security communities. The additional lead time provided by this system will increase the speed, accuracy, and efficacy of humanitarian disaster relief, helping to save lives and livelihoods.

We evaluate the impact of Gravity Recovery and Climate Experiment data assimilation (GRACE-DA) on seasonal hydrological forecast initialization over the United States, focusing on groundwater storage. GRACE-based terrestrial water storage (TWS) estimates are assimilated into a land surface model for the 2003–16 period. Three-month hindcast (i.e., forecast of past events) simulations are initialized using states from the reference (no data assimilation) and GRACE-DA runs. Differences between the two initial hydrological condition (IHC) sets are evaluated for two forecast techniques at 305 wells where depth to water table measurements are available. Results show that using GRACE-DA-based IHC improves seasonal groundwater forecast performance in terms of both RMSE and correlation. While most regions show improvement, degradation is common in the High Plains, where withdrawals for irrigation practices affect groundwater variability more strongly than the weather variability, which demonstrates the need for simulating such activities. These findings contribute to recent efforts toward an improved U.S. drought monitoring and forecast system.

The current study provides the irrigation water estimate based on incorporation of satellite-derived irrigation scheme and crop datasets into the NASA-Land Information System Framework (LISF) in India. NOAH 3.3 land surface model within NASA-LISF was run at 0.05-degree resolution for nine years from 2011 to 2019. The irrigation scheme accurately captures the seasonality and the two growing seasons that is December-March and August-November. The MODIS leaf area index product helps to regulate the seasonality and estimated irrigation amount and timing is based on 50% depletion of soil moisture at the field capacity in the rootzone. The results show that the evapotranspiration (ET) and latent heat flux (LE) have increased significantly in the cropped region with improvement in correlation with the MODIS ET and LE products. The study also shows an improvement in soil moisture simulation at the test sites (Varanasi and Gujarat). Besides, successfully demonstrating the irrigation timing and quantity, the present study can also be relevant to hydrological and energy fluxes studies of areas that still lack proper quantification of agricultural practices utilizing irrigation.

Representation of irrigation in Earth System Models has advanced over the past decade, yet large uncertainties persist in the effective simulation of irrigation practices, particularly over locations where the on-ground practices and climate impacts are less reliably known. Here we investigate the utility of assimilating remotely sensed vegetation data for improving irrigation water use and associated fluxes within a land surface model. We show that assimilating optical sensor-based leaf area index estimates significantly improves the simulation of irrigation water use when compared to the USGS ground reports. For heavily irrigated areas, assimilation improves the evaporative fluxes and gross primary production (GPP) simulations, with the median correlation increasing by 0.1-1.1 and 0.3-0.6, respectively, as compared to the reference datasets. Further, bias improvements in the range of 14-35 mm mo

Terrestrial water storage (TWS) provides important information on terrestrial hydroclimate and may have value for seasonal forecasting because of its strong persistence. We use the NASA Hydrological Forecast and Analysis System (NHyFAS) to investigate TWS forecast skill over Africa and assess its value for predicting vegetation activity from satellite estimates of leaf area index (LAI). Forecast skill is high over East and Southern Africa, extending up to 3–6 months in some cases, with more modest skill over West Africa. Highest skill generally occurs during the dry season or beginning of the wet season when TWS anomalies from the previous wet season are most likely to carry forward in time. In East Africa, this occurs prior to and during the transition into the spring “Long Rains” from January–March, while in Southern Africa this period of highest skill starts at the beginning of the dry season in April and extends through to the start of the wet season in October. TWS is highly and positively correlated with LAI, and a logistic regression model shows high cross-validation skill in predicting above or below normal LAI using TWS. Combining the LAI regression model with the NHyFAS forecasts, 1-month lead LAI predictions have high accuracy over East and Southern Africa, with reduced but significant skill at 3-month leads over smaller sub-regions. This highlights the potential value of TWS as an additional source of information for seasonal forecasts over Africa, with direct applications to some of the most vulnerable agricultural regions on the continent.

Recent annual trends of precipitation and terrestrial water storage (TWS) in West Africa have been increasing over the past decade. Despite a significant impact of soil moisture on the TWS in West Africa, there is little research on the recent spatial and temporal behaviors of surface soil moisture (SSM) along with the hydrological trends and variability in West Africa. In this study, we assimilate TWS estimates from the Gravity Recovery and Climate Experiment (GRACE) mission into the Catchment Land Surface Model (CLSM) and evaluate its impacts on SSM simulations for the years, 2002–2017. The evaluation is performed using reference datasets: the African Monsoon Multidisciplinary Analysis (AMMA) in situ soil moisture observations, three currently available microwave satellite SSM observations from the Advanced Scatterometer (ASCAT), the Soil Moisture Ocean Salinity (SMOS), and the Soil Moisture Active Passive (SMAP) satellites and also the triple collocation analysis (TCA). Overall, modeled SSM shows good agreement with the reference datasets in terms of the anomaly SSM correlations. However, both modeled and ASCAT SSM are limited in their representation of the drying rates, as observed by ground observations, SMOS and SMAP estimates. Further, GRACE data assimilation results in improved SSM simulations in the humid regions with large annual TWS variability. This study demonstrates the utility of land data assimilation to inform hydrological conditions in West Africa, where soil moisture monitoring is necessary for water resource and livestock management.

Terrestrial hydrologic trends over the conterminous United States are estimated for 1980–2015 using the National Climate Assessment Land Data Assimilation System (NCA-LDAS) reanalysis. NCA-LDAS employs the uncoupled Noah version 3.3 land surface model at 0.125° × 0.125° forced with NLDAS-2 meteorology, rescaled Climate Prediction Center precipitation, and assimilated satellite-based soil moisture, snow depth, and irrigation products. Mean annual trends are reported using the nonparametric Mann–Kendall test at p < 0.1 significance. Results illustrate the interrelationship between regional gradients in forcing trends and trends in other land energy and water stores and fluxes. Mean precipitation trends range from +3 to +9 mm yr−1 in the upper Great Plains and Northeast to −1 to −9 mm yr−1 in the West and South, net radiation flux trends range from +0.05 to +0.20 W m−2 yr−1 in the East to −0.05 to −0.20 W m−2 yr−1 in the West, and U.S.-wide temperature trends average about +0.03 K yr−1. Trends in soil moisture, snow cover, latent and sensible heat fluxes, and runoff are consistent with forcings, contributing to increasing evaporative fraction trends from west to east. Evaluation of NCA-LDAS trends compared to independent data indicates mixed results. The RMSE of U.S.-wide trends in number of snow cover days improved from 3.13 to 2.89 days yr−1 while trend detection increased 11%. Trends in latent heat flux were hardly affected, with RMSE decreasing only from 0.17 to 0.16 W m−2 yr−1, while trend detection increased 2%. NCA-LDAS runoff trends degraded significantly from 2.6 to 16.1 mm yr−1 while trend detection was unaffected. Analysis also indicated that NCA-LDAS exhibits relatively more skill in low precipitation station density areas, suggesting there are limits to the effectiveness of satellite data assimilation in densely gauged regions. Overall, NCA-LDAS demonstrates capability for quantifying physically consistent, U.S. hydrologic climate trends over the satellite era.

Research in meteorological prediction on sub-seasonal to seasonal (S2S) timescales has seen growth in recent years. Concurrent with this, demand for seasonal drought forecasting has risen. While there is obvious synergy between these fields, S2S meteorological forecasting has typically focused on low resolution global models, while the development of drought can be sensitive to the local expression of weather anomalies and their interaction with local surface properties and processes. This suggests that downscaling might play an important role in the application of meteorological S2S forecasts to skillful forecasting of drought. Here, we apply the Generalized Analog Regression Downscaling (GARD) algorithm to downscale meteorological hindcasts from the NASA Goddard Earth Observing System (GEOS) global S2S forecast system. Downscaled meteorological fields are then applied to drive offline simulations with the Catchment Land Surface Model (CLSM) to forecast United States Drought Monitor (USDM) style drought indicators derived from simulated surface hydrology variables. We compare the representation of drought in these downscaled hindcasts to hindcasts that are not downscaled, using the North American Land Data Assimilation System Phase 2 (NLDAS-2) dataset as an observational reference. We find that downscaling using GARD improves hindcasts of temperature and temperature anomalies, but the results for precipitation are mixed and generally small. Overall, GARD downscaling led to improved hindcast skill for total drought across the Contiguous United States (CONUS), and improvements were greatest for extreme (D3) and exceptional (D4) drought categories.

The majority of people in East Africa rely on the agro-pastoral system for their livelihood, which is highly vulnerable to droughts and flooding. Agro-pastoral droughts are endemic to the region and are considered the main natural hazard that contributes to food insecurity. Drought begins with rainfall deficit, gradually leading to soil moisture deficit, higher land surface temperature, and finally impacts to vegetation growth. Therefore, monitoring vegetation conditions is essential in understanding the progression of drought, potential effects on food security, and providing early warning information needed for drought mitigation decisions. Because vegetation processes couple the land and atmosphere, monitoring of vegetation conditions requires consideration of both water provision and demand. While there is consensus in using either the Normalized Difference Vegetation Index (NDVI) or evapotranspiration (ET) for vegetation monitoring, a comprehensive assessment optimizing the use of both has not yet been done. Moreover, the evaluation methods for understanding the relationships between NDVI and ET for vegetation monitoring are also limited. Taking these gaps into account we have developed a framework to optimize vegetation monitoring using both NDVI and ET by identifying where they perform the best by using triple collocation and cross-correlation methods. We estimated the random error structure in Moderate Resolution Imaging Spectroradiometer (MODIS) NDVI; ET from the Operational Simplified Surface Energy Balance (SSEBop) model; and ET from land surface models (LSMs). LSM ET and SSEBop ET have been found to be better indicators for vegetation monitoring during extreme drought events, while NDVI could provide better information on vegetation condition during wetter than normal conditions. The random error structures of these variables suggest that LSM ET is most likely to provide important information for vegetation monitoring over low and high ends of the vegetation fraction areas. Over moderate vegetative areas, any of these variables could provide important vegetation information for drought characterization and food security assessments. While this study provides a framework for optimizing vegetation monitoring for drought and food security assessments over East Africa, the framework can be adopted to optimize vegetation monitoring over any other drought and food insecure region of the world.

Though coarse in spatial resolution, the nearly all weather measurements from passive microwave sensors can help in improving the spatio‐temporal coverage of optical and thermal infrared sensors for monitoring vegetation changes on the land surface. This study demonstrates the use of vegetation optical depth (VOD) retrievals from the Soil Moisture Active Passive mission for capturing the vegetation alterations from the recent 2019 to 2020 Australian bushfires and drought. The impact of vegetation disturbances on terrestrial water budget is examined by assimilating the VOD retrievals into a dynamic phenology model. The results demonstrate that assimilating VOD observations lead to improved simulation of evapotranspiration, runoff, and soil moisture states. The study also demonstrates that the vegetation changes from the 2019 to 2020 Australian drought and fires led to significant modifications in the partitioning of evaporative and runoff fluxes, resulting in increased bare soil evaporation, reduced transpiration, and higher runoff.

This work provides an envisioned overview of scientific collaboration among multiple United States agencies including the National Aeronautics and Space Administration (NASA), U.S. Army Engineer Research and Development Center (ERDC), Oak Ridge National Laboratory (ORNL), and National Geospatial-Intelligence Agency (NGA) for the integration of existing data and model capabilities to support global scale water security applications. The primary objective is to develop a high-resolution, operational streamflow and flood forecasting system at the global scale, leveraging multiple process-based models, remote sensing data assimilation, and high-performance computing techniques. We present a preliminary case study that demonstrates the integration of the modeling framework using NASA’s Land Information System (LIS), ERDC’s Streamflow Prediction Tool (SPT), and ORNL’s GPU-accelerated 2D flood model (TRITON). Using the high-resolution terrain data from NGA, a historic flood event that occurred in March 2019 at Offutt Air Force Base in Nebraska, USA, was simulated on ORNL’s supercomputer, Summit. This benchmark test case is used to validate the modeling framework and to help establish a roadmap for the expanded modeling efforts at the global scale. In a broader sense, the proposed infrastructure will enable decision-makers to address issues such as transboundary water conflicts, flood and drought monitoring, and sustainable water resources management and to study their impacts on human, water-energy and natural systems in the short, medium and long term.

To support global drought early warning, the Climate Hazards Center (CHC) at the University of California, Santa Barbara developed the Climate Hazards center InfraRed Precipitation with Stations (CHIRPS) dataset, in collaboration with the US Geological Survey and NASA SERVIR. Specifically designed to support early warning applications, CHIRPS has high a spatial resolution (0.05°), a long period of record (1981 to the near present), and relatively low latencies. Here we will describe a brief formal analysis of distributional bias in CHIRPS2.0. This analysis reveals, as expected, that CHIRPS2.0 means are very similar to observed station data. However, a closer look suggests that low precipitation values are underestimated and high values are over-estimated in the CHIRPS2.0. We describe a potential correction for this below.

The Noah land surface model with multiple parameterization options (Noah-MP) includes a routine for the dynamic simulation of vegetation carbon assimilation and soil carbon decomposition processes. To use remote sensing observations of vegetation to constrain simulations from this model, it is necessary first to understand the sensitivity of the model to its parameters. This is required for efficient parameter estimation, which is both a valuable way to use observations and also a first or concurrent step in many state-updating data assimilation procedures. We use variance decomposition to assess the sensitivity of estimates of sensible heat, latent heat, soil moisture, and net ecosystem exchange made by certain standard Noah-MP configurations that include the dynamic simulation of vegetation and carbon to 43 primary user-specified parameters. This is done using 32 years’ worth of data from 10 international FluxNet sites. Findings indicate that there are five soil parameters and six (or more) vegetation parameters (depending on the model configuration) that act as primary controls on these states and fluxes.

Improved understanding of the water balance in the Blue Nile is of critical importance because of increasingly frequent hydroclimatic extremes under a changing climate. The intercomparison and evaluation of multiple land surface models (LSMs) associated with different meteorological forcing and precipitation datasets can offer a moderate range of water budget variable estimates. In this context, two LSMs, Noah version 3.3 (Noah3.3) and Catchment LSM version Fortuna 2.5 (CLSMF2.5) coupled with the Hydrological Modeling and Analysis Platform (HyMAP) river routing scheme are used to produce hydrological estimates over the region. The two LSMs were forced with different combinations of two reanalysis-based meteorological datasets from the Modern-Era Retrospective analysis for Research and Applications datasets (i.e., MERRA-Land and MERRA-2) and three observation-based precipitation datasets, generating a total of 16 experiments. Modeled evapotranspiration (ET), streamflow, and terrestrial water storage estimates were evaluated against the Atmosphere-Land Exchange Inverse (ALEXI) ET, in-situ streamflow observations, and NASA Gravity Recovery and Climate Experiment (GRACE) products, respectively. Results show that CLSMF2.5 provided better representation of the water budget variables than Noah3.3 in terms of Nash-Sutcliffe coefficient when considering all meteorological forcing datasets and precipitation datasets. The model experiments forced with observation-based products, the Climate Hazards group Infrared Precipitation with Stations (CHIRPS) and the Tropical Rainfall Measuring Mission (TRMM) Multi-Satellite Precipitation Analysis (TMPA), outperform those run with MERRA-Land and MERRA-2 precipitation. The results presented in this paper would suggest that the Famine Early Warning Systems Network (FEWS NET) Land Data Assimilation System incorporate CLSMF2.5 and HyMAP routing scheme to better represent the water balance in this region.

To assess growing season conditions where ground based observations are limited or unavailable, food security and agricultural drought monitoring analysts rely on publicly available remotely sensed rainfall and vegetation greenness. There are also remotely sensed soil moisture observations from missions like the European Space Agency (ESA) Soil Moisture and Ocean Salinity (SMOS) and NASA's Soil Moisture Active Passive (SMAP), however these time series are still too short to conduct studies that demonstrate the utility of these data for operational applications, or to provide historical context for extreme wet or dry events. To promote the use of remotely sensed soil moisture in agricultural drought and food security monitoring, we use East Africa as a case study to evaluate the quality of a 30+ year time series of merged active-passive microwave soil moisture from the ESA Climate Change Initiative (CCI-SM). Compared to the Normalized Difference Vegetation index (NDVI) and modeled soil moisture products, we found substantial spatial and temporal gaps in the early part of the CCI-SM record, with adequate data coverage beginning in 1992. From this point forward, growing season CCI-SM anomalies were well correlated (R0.5) with modeled, seasonal soil moisture, and in some regions, NDVI. We use correlation analysis and qualitative comparisons at seasonal time scales to show that remotely sensed soil moisture can add information to a convergence of evidence framework that traditionally relies on rainfall and NDVI in moderately vegetated regions.

The effective applications of land surface model (LSM) 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 pre-processor 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, 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 data sets and data analytics algorithms to aid land surface modelling and data assimilation.

Accurate determination of snow conditions is important for several water management applications, partly because of the significant influence of snowmelt on seasonal streamflow prediction. This article examines an approach using snow cover area (SCA) observations as snow detection constraints during the assimilation of snow depth retrievals from passive microwave sensors. Two different SCA products [the Interactive Multisensor Snow and Ice Mapping System (IMS) and the Moderate Resolution Imaging Spectroradiometer (MODIS)] are employed jointly with the snow depth retrievals from a variety of sensors for data assimilation in the Noah land surface model. The results indicate that the use of MODIS data is effective in obtaining added improvements (up to 6% improvement in aggregate RMSE) in snow depth fields compared to assimilating passive microwave data alone, whereas the impact of IMS data is small. The improvements in snow depth fields are also found to translate to small yet systematic improvements in streamflow estimates, especially over the western United States, the upper Missouri River, and parts of the Northeast and upper Mississippi River. This study thus demonstrates a simple approach for exploiting the information from SCA observations in data assimilation.

In snowmelt-driven river systems, it is critical to enable reliable predictions of the spatiotemporal variability in seasonal snowpack to support local and regional water management. Previous studies have shown that assimilating satellite-station blended snow depth data sets can lead to improved snow predictions, which however do not always translate into improved streamflow predictions, especially in complex mountain regions. In this study, we explore how an existing optimal interpolation-based blending strategy can be enhanced to reduce biases in satellite snow depth products for improving streamflow predictions. Two major new considerations are explored, including: (1) incorporating terrain aspect and (2) incorporating areal snow coverage information. The methodology is applied to the bias reduction of the Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E) snow depth estimates, which are then assimilated into the Noah land surface model via the ensemble Kalman Filtering (EnKF) for streamflow predictions in the Upper Colorado River Basin. Our results indicate that using only observations from low-elevation stations such as the Global Historical Climatology Network (GHCN) in the bias correction can lead to underestimation in streamflow, while using observations from high-elevation stations (e.g., the Snow Telemetry (SNOTEL) network) along with terrain aspect is critically important for achieving reliable streamflow predictions. Additionally incorporating areal snow coverage information from the Moderate Resolution Imaging Spectroradiometer (MODIS) can slightly improve the streamflow results further.

The Soil Moisture Active Passive (SMAP) mission will provide soil moisture data with unprecedented accuracy, resolution, and coverage, enabling models to better track agricultural drought and estimate yields. In turn, this information can be used to shape policy related to food and water from commodity markets to humanitarian relief efforts. New data alone, however, do not translate to improvements in drought and yield forecasts. New tools will be needed to transform SMAP data into agriculturally meaningful products. The objective of this study is to evaluate the possibility and efficiency of replacing the rainfall-derived soil moisture component of a crop water stress index with SMAP data. The approach is demonstrated with 0.1°-resolution, ~10-day microwave soil moisture from the European Space Agency and simulated soil moisture from the Famine Early Warning Systems Network Land Data Assimilation System. Over a West Africa domain, the approach is evaluated by comparing the different soil moisture estimates and their resulting Water Requirement Satisfaction Index values from 2000 to 2010. This study highlights how the ensemble of indices performs during wet versus dry years, over different land-cover types, and the correlation with national-level millet yields. The new approach is a feasible and useful way to quantitatively assess how satellite-derived rainfall and soil moisture track agricultural water deficits. Given the importance of soil moisture in many applications, ranging from agriculture to public health to fire, this study should inspire other modeling communities to reformulate existing tools to take advantage of SMAP data.

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 FEWS NET’s 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. Machine-accessible metadata file describing the reported data (ISA-Tab format)

The accurate knowledge of soil moisture and snow conditions is important for the skillful characterization of agricultural and hydrologic droughts, which are defined as deficits of soil moisture and streamflow, respectively. This article examines the influence of remotely sensed soil moisture and snow depth retrievals toward improving estimates of drought through data assimilation. Soil moisture and snow depth retrievals from a variety of sensors (primarily passive microwave based) are assimilated separately into the Noah land surface model for the period of 1979–2011 over the continental United States, in the North American Land Data Assimilation System (NLDAS) configuration. Overall, the assimilation of soil moisture and snow datasets was found to provide marginal improvements over the open-loop configuration. Though the improvements in soil moisture fields through soil moisture data assimilation were barely at the statistically significant levels, these small improvements were found to translate into subsequent small improvements in simulated streamflow. The assimilation of snow depth datasets were found to generally improve the snow fields, but these improvements did not always translate to corresponding improvements in streamflow, including some notable degradations observed in the western United States. A quantitative examination of the percentage drought area from root-zone soil moisture and streamflow percentiles was conducted against the U.S. Drought Monitor data. The results suggest that soil moisture assimilation provides improvements at short time scales, both in the magnitude and representation of the spatial patterns of drought estimates, whereas the impact of snow data assimilation was marginal and often disadvantageous.

Current satellite-based remote-sensing approaches are largely incapable of estimating precipitation over snow cover. This note reports a proof-of-concept study of a new satellite-based approach to the estimation of precipitation over snow-covered surfaces. The method is based on the principle that precipitation can be inferred from the changes in the snow water equivalent of the snowpack. Using satellite-based snow water equivalent measurements, we derived daily precipitation amounts for the northern hemisphere for three snow-accumulation seasons, and evaluated these against independent reference datasets. The new precipitation estimates captured realistic-looking storm events over largely un-instrumented regions. However, the data are noisy and, on a seasonal scale, the amount of precipitation is believed to be underestimated. Nevertheless, current uncertainty in snow measurements, albeit large (50–100%), is still lower than direct precipitation measurements over snow (100–140%) and therefore this approa...

An evapotranspiration (ET) ensemble composed of 36 land surface model (LSM) experiments and four diagnostic datasets (GLEAM, ALEXI, MOD16, and FLUXNET) is used to investigate uncertainties in ET estimate over five climate regions in West Africa. Diagnostic ET datasets show lower uncertainty estimates and smaller seasonal variations than the LSM-based ET values, particularly in the humid climate regions. Overall, the impact of the choice of LSMs and meteorological forcing datasets on the modeled ET rates increases from north to south. The LSM formulations and parameters have the largest impact on ET in humid regions, contributing to 90% of the ET uncertainty estimates. Precipitation contributes to the ET uncertainty primarily in arid regions. The LSM-based ET estimates are sensitive to the uncertainty of net radiation in arid region and precipitation in humid region. This study serves as support for better determining water availability for agriculture and livelihoods in Africa with earth observations and land surface models.

Eleven years of daily 500 m gridded Terra Moderate Resolution Imaging Spectroradiometer (MODIS) (MOD10A1) snow cover fraction (SCF) data are evaluated in terms of snow presence detection in Colorado and Washington states. The SCF detection validation study is performed using in-situ measurements and expressed in terms of snow and land detection and misclassification frequencies. A major aspect addressed in this study involves the shifting of pixel values in time due to sensor viewing angles and gridding artifacts of MODIS sensor products. To account for this error, 500 m gridded pixels are grouped and aggregated to different-sized areas to incorporate neighboring pixel information. With pixel aggregation, both the probability of detection (POD) and the false alarm ratios increase for almost all cases. Of the false negative (FN) and false positive values (referred to as the total error when combined), FN estimates dominate most of the total error and are greatly reduced with aggregation. The greatest POD increases and total error reductions occur with going from a single 500 m pixel to 3×3-pixel averaged areas. Since the MODIS SCF algorithm was developed under ideal conditions, SCF detection is also evaluated for varying conditions of vegetation, elevation, cloud cover and air temperature. Finally, using a direct insertion data assimilation approach, pixel averaged MODIS SCF observations are shown to improve modeled snowpack conditions over the single pixel observations due to the smoothing of more error-prone observations and more accurately snow-classified pixels. Copyright © 2012 John Wiley & Sons, Ltd.

Climate change is expected to accelerate the hydrologic cycle, increase the fraction of precipitation that is rain, and enhance snowpack melting. The enhanced hydrological cycle is also expected to increase snowfall amounts due to increased moisture availability. These processes are examined in this paper in the Colorado Headwaters region through the use of a coupled high-resolution climate–runoff model. Four high-resolution simulations of annual snowfall over Colorado are conducted. The simulations are verified using Snowpack Telemetry (SNOTEL) data. Results are then presented regarding the grid spacing needed for appropriate simulation of snowfall. Finally, climate sensitivity is explored using a pseudo–global warming approach. The results show that the proper spatial and temporal depiction of snowfall adequate for water resource and climate change purposes can be achieved with the appropriate choice of model grid spacing and parameterizations. The pseudo–global warming simulations indicate enhanced snowfall on the order of 10%–25% over the Colorado Headwaters region, with the enhancement being less in the core headwaters region due to the topographic reduction of precipitation upstream of the region (rain-shadow effect). The main climate change impacts are in the enhanced melting at the lower-elevation bound of the snowpack and the increased snowfall at higher elevations. The changes in peak snow mass are generally near zero due to these two compensating effects, and simulated wintertime total runoff is above current levels. The 1 April snow water equivalent (SWE) is reduced by 25% in the warmer climate, and the date of maximum SWE occurs 2–17 days prior to current climate results, consistent with previous studies.

Snow depletion curves (SDC) are functions that are used to show the relationship between snow covered area and snow depth or water equivalent. Previous snow cover data assimilation (DA) studies have used theoretical SDC models as observation operators to map snow depth to snow cover fraction (SCF). In this study, a new approach is introduced that uses snow water equivalent (SWE) observations and satellite-based SCF retrievals to derive SDC relationships for use in an Ensemble Kalman filter (EnKF) to assimilate snow cover estimates. A histogram analysis is used to bin the SWE observations, which the corresponding SCF observations are then averaged within, helping to constrain the amount of data dispersion across different temporal and regional conditions. Logarithmic functions are linearly regressed with the binned average values, for two U.S. mountainous states: Colorado and Washington. The SDC-based logarithmic functions are used as EnKF observation operators, and the satellite-based SCF estimates are assimilated into a land surface model. Assimilating satellite-based SCF estimates with the observation-based SDC shows a reduction in SWE-related RMSE values compared to the model-based SDC functions. In addition, observation-based SDC functions were derived for different intra-annual and physiographic conditions, and landcover and elevation bands. Lower SWE-based RMSE values are also found with many of these categorical observation-based SDC EnKF experiments. All assimilation experiments perform better than the open-loop runs, except for the Washington region&rsquo, s 2004&ndash, 2005 snow season, which was a major drought year that was difficult to capture with the ensembles and observations.

This paper presents results of four high resolution simulations of annual snowfall over Colorado, U.S.A. The results are verified using SNOTEL data. Sensitivity to model resolution is also explored. The results show that proper spatial and temporal depictions of snowfall adequate for water resource and climate change purposes can be achieved with the appropriate choice of model resolution and physical parameterizations.

Many studies have developed snow process understanding by exploring the impact of snow model complexity on simulation performance. This paper revisits this topic using several recently developed land surface models, including the Simplified Simple Biosphere Model (SSiB); Noah; Variable Infiltration Capacity (VIC); Community Land Model, version 3 (CLM3); Snow Thermal Model (SNTHERM); and new field measurements from the Cold Land Processes Field Experiment (CLPX). Offline snow cover simulations using these five snow models with different physical complexity are performed for the Rabbit Ears Buffalo Pass (RB), Fraser Experimental Forest headquarters (FHQ), and Fraser Alpine (FA) sites between 20 September 2002 and 1 October 2003. These models simulate the snow accumulation and snowpack ablation with varying skill when forced with the same meteorological observations, initial conditions, and similar soil and vegetation parameters. All five models capture the basic features of snow cover dynamics but show remarkable discrepancy in depicting snow accumulation and ablation, which could result from uncertain model physics and/or biased forcing. The simulated snow depth in SSiB during the snow accumulation period is consistent with the more complicated CLM3 and SNTHERM; however, early runoff is noted, owing to neglected water retention within the snowpack. Noah is consistent with SSiB in simulating snow accumulation and ablation at RB and FA, but at FHQ, Noah underestimates snow depth and snow water equivalent (SWE) as a result of a higher net shortwave radiation at the surface, resulting from the use of a small predefined maximum snow albedo. VIC and SNTHERM are in good agreement with each other, and they realistically reproduce snow density and net radiation. CLM3 is consistent with VIC and SNTHERM during snow accumulation, but it shows early snow disappearance at FHQ and FA. It is also noted that VIC, CLM3, and SNTHERM are unable to capture the observed runoff timing, even though the water storage and refreezing effects are included in their physics. A set of sensitivity experiments suggest that Noah’s snow simulation is improved with a higher maximum albedo and that VIC exhibits little improvement with a larger fresh snow albedo. There are remarkable differences in the vegetation impact on snow simulation for each snow model. In the presence of forest cover, SSiB shows a substantial increase in snow depth and SWE, Noah and VIC show a slight change though VIC experiences a later onset of snowmelt, and CLM3 has a reduction in its snow depth. Finally, we observe that a refined precipitation dataset significantly improves snow simulation, emphasizing the importance of accurate meteorological forcing for land surface modeling.

Despite recent advances in land surface modeling 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-of-the-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 1° spatial resolution were used to run 14 LSMs. Precipitation datasets that have been rescaled to match monthly Global Precipitation Climatology Project (GPCP) and Global 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 and Gravity Recovery and Climate Experiment (GRACE) TWS estimates 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−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.

[1] Eight years (2002–2010) of Advanced Microwave Scanning Radiometer–EOS (AMSR-E) snow water equivalent (SWE) retrievals and Moderate Resolution Imaging Spectroradiometer (MODIS) snow cover fraction (SCF) observations are assimilated separately or jointly into the Noah land surface model over a domain in Northern Colorado. A multiscale ensemble Kalman filter (EnKF) is used, supplemented with a rule-based update. The satellite data are either left unscaled or are scaled for anomaly assimilation. The results are validated against in situ observations at 14 high-elevation Snowpack Telemetry (SNOTEL) sites with typically deep snow and at 4 lower-elevation Cooperative Observer Program (COOP) sites. Assimilation of coarse-scale AMSR-E SWE and fine-scale MODIS SCF observations both result in realistic spatial SWE patterns. At COOP sites with shallow snowpacks, AMSR-E SWE and MODIS SCF data assimilation are beneficial separately, and joint SWE and SCF assimilation yields significantly improved root-mean-square error and correlation values for scaled and unscaled data assimilation. In areas of deep snow where the SNOTEL sites are located, however, AMSR-E retrievals are typically biased low and assimilation without prior scaling leads to degraded SWE estimates. Anomaly SWE assimilation could not improve the interannual SWE variations in the assimilation results because the AMSR-E retrievals lack realistic interannual variability in deep snowpacks. SCF assimilation has only a marginal impact at the SNOTEL locations because these sites experience extended periods of near-complete snow cover. Across all sites, SCF assimilation improves the timing of the onset of the snow season but without a net improvement of SWE amounts.

The Environmental Protection Agency (EPA) has estimated that over 20,000 water bodies within the United States do not meet water quality standards. Ninety percent of the impairments are typically caused by nonpoint sources. One of the regulations in the Clean Water Act of 1972 requires States to monitor the Total Maximum Daily Load (TMDL), or the amount of pollution that can be carried by a water body before it is determined to be "polluted", for any watershed in the U.S.. In response to this mandate, the EPA developed Better Assessment Science Integrating Nonpoint Sources (BASINS) as a Decision Support Tool (DST) for assessing pollution and to guide the decision making process for improving water quality. One of the models in BASINS, the Hydrological Simulation Program -- Fortran (HSPF), computes daily stream flow rates and pollutant concentration at each basin outlet. By design, precipitation and other meteorological data from weather stations serve as standard model input. In practice, these stations may be unable to capture the spatial heterogeneity of precipitation events especially if they are few and far between. An attempt was made to resolve this issue by substituting station data with NASA modified/NOAA precipitation data. Using these data within HSPF, stream flow was calculated for seven watersheds in the Chesapeake Bay Basin during low flow periods, convective storm periods, and annual flows. In almost every case, the modeling performance of HSPF increased when using the NASA-modified precipitation data, resulting in better stream flow statistics and, ultimately, in improved water quality assessment.

[1] This paper describes a spin-up experiment conducted over South America using the Simplified Simple Biosphere (SSiB) land surface model to study the process of model adjustment to atmospheric forcing data. The experiment was carried out as a precursor to the use of SSiB in a South American Land Data Assimilation System (SALDAS). The results from an 11 year long recursive simulation using three different initial conditions of soil wetness (control, wet and dry) are examined. The control run was initiated by interpolation of the NCEP/DOE Global Reanalysis-2 (NCEP/DOE R-2) soil moisture data set. In each case the time required for the model to reach equilibrium was calculated. The wet initialization leads to a faster adjustment of the soil moisture field, followed by the control and then the dry initialization. Overall, the final spin-up states using the SSiB-based SALDAS are generally wetter than both the NCEP/DOE R-2 and the Centro de Previsao do Tempo e Estudos Climaticos (CPTEC–Brazilian Center for Weather Forecast and Climate Studies) operational initial soil moisture states, consequently modeled latent heat is higher and sensible heat lower in the final year of simulation when compared with the first year. Selected regions, i.e., in semiarid northeastern Brazil, the transition zone to the south of the Amazon tropical forest, and the central Andes were studied in more detail because they took longer to spin up (up to 56 months) when compared with other areas (less than 24 months). It is shown that there is a rapid change in the soil moisture in all layers in the first 2 months of simulation followed by a subsequent slow and steady adjustment: This could imply there are increasing errors in medium range simulations. Spin-up is longest where frozen soil is present for long periods such as in the central Andes.

NASA/GSFC has developed with other groups a Land Data Assimilation System (LDAS) to output water and energy budgets for the primary purpose of improving weather and climate prediction. However, LDAS water and energy cycle outputs also may be coupled with other information to help with a wide range of water resources applications. For example, LDAS results may be used for water availability and quality, agricultural management and forecasting, assessment and prediction of snowmelt runoff, and flood and drought impact and prediction. Specifically, LDAS uses various satellites and ground based observations within a land surface modeling and data assimilation framework to produce optimal output fields of terrestrial energy, water and carbon fluxes. Current land surface outputs are gridded at 1/4° resolution globally and 1/8° for North America with work in progress to convert to a 1-km global grid. Integrated modeling, observations and data assimilation at various spatial and temporal scales helps LDAS to quantify terrestrial water, energy and biogeochemical processes. LDAS applications described in this paper are aimed at improving weather and seasonal forecasts. In addition, we also summarize the use of LDAS data to assist critical needs specified by the U.S. Bureau of Reclamation water resources management for selected basins in the western U.S.

To improve the efficiency of water management and irrigation scheduling in the Rio Grande River basin, the U.S. Bureau of Reclamation helped create the Agricultural Water Resources Decision Support (AWARDS) system. Through its Evapotranspiration (ET) Toolbox interface, the AWARDS system provides guidance to local farmers on when and where to deliver water to the crops. The ET Toolbox is based on water usage estimates (evapotranspiration and open water evaporation) on a grid cell basis (4 kmtimes4 km). Currently, crop water use estimates are determined using a modified-Penman ET approach. To improve upon this parameterization, we use the Community Land Model (CLM2.0) within the LDAS system downscaled to a 1 km grid cell resolution. Our work aims to improve evapotranspiration and soil moisture estimates in the Rio Grande River basin through the improvement of the CLM2.0 parameterization of surface processes. We specifically intend to assimilate up to four land surface temperature (LST) observations per day from the Terra/Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) instruments, at 1 km resolution, into the CLM2.0 land surface model. To verify the performance of the assimilation approach we deploy flux towers at two different sites. We focus our ground data collection on riparian and agricultural areas. Ultimately, our results should improve the daily forecasts of vegetation water requirements and, when integrated into operational decision support tools, aid water resource managers in making flood and drought assessments and predictions