Search

On June 23, 2023, the Advanced Information Systems Technology (AIST) program conducted grouped annual reviews of all its Earth System Digital Twins (ESDT) projects. This report regroups all these technical presentations.

Rain-on-snow (ROS) over snow-dominated regions such as High Mountain Asia (HMA) modulates snowmelt and runoff and is key contributor in influencing water availability and hazards (e.g., floods and landslides). We studied the trends in ROS in HMA over the past two decades from 2001 to 2018 using the land surface model Noah-MP driven by an ensemble precipitation data set. Our results show that changes in precipitation phase and rainfall are altering ROS. Because of the strong physical heterogeneity and atmospheric dynamics of HMA, ROS characteristics and trends are region-dependent and ROS occurs predominantly over the Indus, Ganges-Brahmaputra, and northwestern basins. In the Indus, ROS representing ∼5% of the annual precipitation and ∼20% of the annual snowmelt, has an increasing trend. This is contrary to the Ganges-Brahmaputra characterized by decreasing ROS trends, where it represents ∼11% of the annual precipitation and ∼60% of the annual snowmelt. In the northwestern basins, ROS has bidirectional trends due to elevation patterns and trends in rainfall, and it constitutes ∼5 to ∼10% of the annual precipitation. Increasing trends in ROS over Indus contribute to reducing the snowpack in late summer, with concerns of reduced water availability and increased groundwater exploitation. Similarly, because of its high amount and contribution to snowmelt, the decreasing ROS trends in the Ganges-Brahmaputra will have consequences of decreased recharge from the headwaters and exacerbated use of groundwater unless increasing trends in rainfall compensate for the decreasing snowmelt. These results provide new insights on ROS-driven changes in the hydrological cycle over HMA.

In a changing climate, the likelihood of hydrologic extremes has been increasing as climate change can impact both means and extremes4 of hydrologic cycle processes, potentially resulting in an increased frequency of floods in some regions and decreases in others. In a warming world, the physical processes that affect hydrologic response, such as rain-on snow runoff events, are also changing, such that the seasonality of streamflow has been shifting. The geography of rain-on-snow runoff events is predicted to move from low to high elevations. In addition to floods, there is also potential for an increase in dry extremes in a warming world with increased drought frequency and occurrences in many parts of the world. The increased frequency of drought and heatwave events is expected to have consequences such as escalating crop failures in future projection scenarios1. Thus, the consensus of literature shows that climate change is increasing the magnitude and frequency of extreme hydrologic events, and the human influence in many of these events is substantial.

In situ hydrologic monitoring over regions most susceptible to food insecurity can be a challenge in current times due to various socio-economic and political issues in combination with environmental factors such as ongoing famine or drought. Hydrologic monitoring and initializing forecasts based on remotely sensed and analyzed data can contribute significantly to early warning in such regions. Routine hydrologic forecasts, as provided by NASA’s Hydrologic Forecasting and Analysis System (NHyFAS), are a recent addition to early warning systems. A custom instance of NHyFAS, termed FLDAS-Forecast, is used by FEWS NET’s Land Data Assimilation System (FLDAS). The FLDAS-Forecast’s dynamic forecasting component was originally set up with Goddard Earth Observing System (GEOS) forecast inputs and has been recently expanded with precipitation forecast forcing from the North American Multi-Model Ensemble (NMME). This paper describes the improvements in seasonal hydrologic forecasts produced with this updated system. Evaluations in this study focus on soil moisture across southern Africa’s growing season. Soil moisture forecasts are benchmarked and evaluated relative to climatology-based forecasts and historic runs, which are driven by observation-based meteorological forcing fields, and they are verified with remotely sensed observations of soil moisture and vegetation. Through multiple deterministic and probabilistic skill assessments, we show that using the larger ensemble of NMME precipitation inputs in the forecast system results in higher quality hydrologic forecasts than are allowed by climatology- or GEOS-only-based forecasts. Further, the near-real-time NMME-based rootzone soil moisture forecasts were able to correctly predict developing drought conditions over southern Africa through late 2019 and into early 2020.

Flash droughts evolve and intensify rapidly under the influence of anomalous atmospheric conditions. In this study, we investigate the role of assimilating remotely sensed soil moisture (SM) and vegetation properties in capturing the evolution and impacts of two flash droughts in the Northern Great Plains. We find that during 2016 drought triggered by anomalously high temperatures and excessive evaporative demands, multivariate data assimilation (DA) of MODIS-derived leaf area index (LAI) and Soil Moisture Active Passive SM within Noah-Multiparameterization model helps capture elevated transpiration at onset. Assimilation of LAI particularly helped model the resulting rapid decline in SM during onset with as high as 10.0% steeper rate of decline compared to the simulation without any assimilation. Modeled-SM anomalies exhibit a 7.5% and 11.7% increase in similarity with Evaporative Stress Index (ESI) data and U.S. Drought Monitor (USDM) maps, respectively. In contrast, during 2017 flash drought driven by record-low precipitation during summers, SM assimilation resulted in largest rates of decline in rootzone SM, as large as 48.4% compared to results from no assimilation. Multivariate DA of SM and LAI results in 6.7% and 14.3% higher spatial similarity with ESI and USDM, respectively, and is necessary to model rapid intensification caused by anomalous precipitation deficits. This study elucidates the need to incorporate multiple observational constraints from remote sensing to effectively capture rapid onset rates, intensification, and severity of flash drought following different propagation mechanisms. This is fundamental for drought early detection to provide a wider window of response and implement efficient mitigation strategies.

Precipitation estimates are highly uncertain in complex regions such as High-Mountain Asia (HMA), where ground measurements are very difficult to obtain, and atmospheric dynamics poorly understood. Though gridded products derived from satellite-based observations and/or reanalysis can provide temporally and spatially distributed estimates of precipitation, there are significant inconsistencies in these products. As such, to date, there is little agreement in the community on the best and most accurate gridded precipitation product in HMA, which is likely area dependent because of HMA’s strong heterogeneities and complex orography. Targeting these gaps, this article presents the development of a consensus ensemble precipitation product using three gridded precipitation datasets (the Integrated Multi-satellitE Retrieals for Global Precipitation Measurement IMERG, the Climate Hazards group Infrared Precipitation with Stations CHIRPS, and the ECMWF Reanalysis ERA5) with a localized probability matched mean (LPM) approach. We evaluate the performance of the LPM estimate along with a simple ensemble mean (EM) estimate to overcome the differences and disparities of the three selected constituent products on long-term averages and trends in HMA. Our analysis demonstrates that LPM reduces the high biases embedded in the ensemble members and provides more realistic spatial patterns compared to EM. LPM is also a good alternative for merging data products with different spatio-temporal resolutions. By filtering disparities among the individual ensemble members, LPM overcomes the problem of a certain product performing well only in a particular area and provides a consensus estimate with plausible temporal trends.

This article describes a new precipitation analysis algorithm developed by NASA for time-sensitive operations at the United States Air Force. Implemented as part of the Land Information System—a land modeling and data assimilation software framework—this NASA–Air Force Precipitation Analysis (NAFPA) combines numerical weather prediction model outputs with rain gauge measurements and satellite estimates to produce global, gridded 3-h accumulated precipitation fields at approximately 10-km resolution. Input observations are subjected to quality control checks before being used by the Bratseth analysis algorithm that converges to optimal interpolation. NAFPA assimilates up to 3.5 million observations without artificial data thinning or selection. To evaluate this new approach, a multiyear reanalysis is generated and intercompared with eight alternative precipitation products across the contiguous United States, Africa, and the monsoon region of eastern Asia. NAFPA yields superior accuracy and correlation over low-latency (up to 14 h) alternatives (numerical weather prediction and satellite retrievals), and often outperforms high-latency (up to 3.5 months) products, although the details for the latter vary by region and product. The development of NAFPA offers a high-quality, near-real-time product for use in meteorological, land surface, and hydrological research and applications.

Ozone (O3) dry deposition is a major O3 sink. As a follow-up study of Huang et al. (2021), we quantify the impact of satellite soil moisture (SM) on model representations of this process when different dry-deposition parameterizations are implemented, based on which the implications for interpreting O3 air pollution levels and assessing the O3 impacts on human and ecosystem health are provided. The SM data from NASA's Soil Moisture Active Passive mission are assimilated into the Noah-Multiparameterization (Noah-MP) land surface model within the NASA Land Information System framework, semicoupled with Weather Research and Forecasting model with online Chemistry (WRF-Chem) regional-scale simulations covering the southeastern US. Major changes in the modeling system used include enabling the dynamic vegetation option, adding the irrigation process, and updating the scheme for the surface exchange coefficient. Two dry-deposition schemes are implemented, i.e., the Wesely scheme and a “dynamic” scheme, in the latter of which dry-deposition parameterization is coupled with photosynthesis and vegetation dynamics. It is demonstrated that, when the dynamic scheme is applied, the simulated O3 dry-deposition velocities vd and their stomatal and cuticular portions, as well as the total O3 fluxes Ft, are larger overall; vd and Ft are 2–3 times more sensitive to the SM changes due to the data assimilation (DA). Further, through case studies at two forested sites with different soil types and hydrological regimes, we highlight that, applying the Community Land Model type of SM factor controlling stomatal resistance (i.e., β factor) scheme in replacement of the Noah-type β factor scheme reduced the v(d) sensitivity to SM changes by ∼75 % at one site, while it doubled this sensitivity at the other site. Referring to multiple evaluation datasets, which may be associated with variable extents of uncertainty, the model performance of vegetation, surface fluxes, weather, and surface O3 concentrations shows mixed responses to the DA, some of which display land cover dependency. Finally, using model-derived concentration- and flux-based policy-relevant O3 metrics as well as their matching exposure–response functions, the relative biomass/crop yield losses for several types of vegetation/crops are estimated to be within a wide range of 1 %–17 %. Their sensitivities to the model's dry-deposition scheme and the implementation of SM DA are discussed.

The Middle East and North Africa (MENA) region has experienced more frequent and severe drought events in recent decades, leading to increasingly pressing concerns over already strained food and water security. An effective drought monitoring and early warning system is thus critical to support risk mitigation and management by countries in the region. Here we investigate the potential for assimilation of leaf area index (LAI) and soil moisture observations to improve the representation of the overall hydrological and carbon cycles and drought by an advanced land surface model. The results reveal that assimilating soil moisture does not meaningfully improve model representation of the hydrological and biospheric processes for this region, but instead it degrades the simulation of the interannual variation in evapotranspiration (ET) and carbon fluxes, mainly due to model weaknesses in representing prognostic phenology. However, assimilating LAI leads to greater improvement, especially for transpiration and carbon fluxes, by constraining the timing of simulated vegetation growth response to evolving climate conditions. LAI assimilation also helps to correct for the erroneous interaction between the prognostic phenology and irrigation during summertime, effectively reducing a large positive bias in ET and carbon fluxes. Independently assimilating LAI or soil moisture alters the categorization of drought, with the differences being greater for more severe drought categories. We highlight the vegetation representation in response to changing land use and hydroclimate as one of the key processes to be captured for building a successful drought early warning system for the MENA region.

A soil moisture retrieval assimilation framework is implemented across South Asia in an attempt to improve regional soil moisture estimation as well as to provide a consistent regional soil moisture dataset. This study aims to improve the spatiotemporal variability of soil moisture estimates by assimilating Soil Moisture Active Passive (SMAP) near-surface soil moisture retrievals into a land surface model. The Noah-MP (v4.0.1) land surface model is run within the NASA Land Information System software framework to model regional land surface processes. NASA Modern-Era Retrospective Analysis for Research and Applications (MERRA2) and Global Precipitation Measurement (GPM) Integrated Multi-satellitE Retrievals (IMERG) provide the meteorological boundary conditions to the land surface model. Assimilation is carried out using both cumulative distribution function (CDF)-corrected (DA-CDF) and uncorrected SMAP retrievals (DA-NoCDF). CDF matching is applied to correct the statistical moments of the SMAP soil moisture retrieval relative to the land surface model. Comparison of assimilated and model-only soil moisture estimates with publicly available in situ measurements highlights the relative improvement in soil moisture estimates by assimilating SMAP retrievals. Across the Tibetan Plateau, DA-NoCDF reduced the mean bias and RMSE by 8.4 % and 9.4 %, even though assimilation only occurred during less than 10 % of the study period due to frozen (or partially frozen) soil conditions. The best goodness-of-fit statistics were achieved for the IMERG DA-NoCDF soil moisture experiment. The general lack of publicly available in situ measurements across irrigated areas limited a domain-wide direct model validation. However, comparison with regional irrigation patterns suggested correction of biases associated with an unmodeled hydrologic phenomenon (i.e., anthropogenic influence via irrigation) as a result of SMAP soil moisture retrieval assimilation. The greatest sensitivity to assimilation was observed in cropland areas. Improvements in soil moisture potentially translate into improved spatiotemporal patterns of modeled evapotranspiration, although limited influence from soil moisture assimilation was observed on modeled processes within the carbon cycle such as gross primary production. Improvement in fine-scale modeled estimates by assimilating coarse-scale retrievals highlights the potential of this approach for soil moisture estimation over data-scarce regions.

This second paper of the two-part series focuses on demonstrating the impact of assimilating satellite-based snow cover and freeze/thaw observations into the hyper-resolution, offline terrestrial modeling system used for the High Mountain Asia (HMA) region from 2003 to 2016. To this end, this study systematically evaluates a total of six sets of 0.01° (∼1 km) model simulations forced by different precipitation forcings, with and without the dual assimilation scheme enabled, at point-scale, basin-scale, and domain-scale. The key variables of interest include surface net shortwave radiation, surface net longwave radiation, skin temperature, near-surface soil temperature, snow depth, snow water equivalent (SWE), and total runoff. First, the point-scale assessment is mainly conducted via evaluating against ground-based measurements. In general, the assimilation enabled estimates are better than no-assimilation counterparts. Second, the basin-scale runoff assessment demonstrates that across three snow-dominated basins, the assimilation enabled experiment yields systematic improvements in all goodness-of-fit statistics through mitigating the negative effects brought by the fixed long-term precipitation correction factors. For example, when forced by the bias-corrected precipitation, the assimilation-enabled experiment improves the bias by 69%, the root-mean-squared error by 30%, and the unbiased root-mean-squared error by 18% (relative to the no-assimilation counterpart). Finally, the domain-scale assessment is conducted via evaluating against satellite-based SWE and skin temperature products. Both sets of domain-scale analysis further corroborate the findings in the point-scale evaluations. Overall, this study suggests the benefits of the proposed multi-variate assimilation system in improving the cryospheric-hydrological process within a land surface model for use in HMA.

he NASA LIS/WRF-Hydro system is a coupled modeling framework that combines the modeling and data assimilation (DA) capabilities of the NASA Land Information System (LIS) with the multi-scale surface hydrological modeling capabilities of the WRF-Hydro model, both of which are widely used in both operations and research. This coupled modeling framework builds on the linkage between land surface models (LSMs), which simulate surface boundary conditions in atmospheric models, and distributed hydrologic models, which simulate horizontal surface and sub-surface flow, adding new land DA capabilities. In the present study, we employ this modeling framework in the Tuolumne River basin in central California. We demonstrate the added value of the assimilation of NASA Airborne Snow Observatory (ASO) snow water equivalent (SWE) estimates in the Tuolumne basin. This analysis is performed in both LIS as an LSM column model and LIS/WRF-Hydro, with hydrologic routing. Results demonstrate that ASO DA in the basin reduced snow bias by as much as 30% from an open-loop (OL) simulation compared to three independent datasets. It also reduces downstream streamflow runoff biases by as much as 40%, and improves streamflow skill scores in both wet and dry years. Analysis of soil moisture and evapotranspiration (ET) also reveals the impacts of hydrologic routing from WRF-Hydro in the simulations, which would otherwise not be resolved in an LSM column model. By demonstrating the beneficial impact of SWE DA on the improving streamflow forecasts, the article outlines the importance of such observational inputs for reservoir operations and related water management applications.

NASA's Advanced Information Systems Technology (AIST) Program is one of several Technology programs managed by the Earth Science Technology Office (ESTO) in the Earth Science Division (ESD). The AIST Program focuses on advanced information systems and novel computer science technologies that will be needed by NASA Earth Science in the next 5 to 10 years. New Observing Strategies is one the three main thrusts of the AIST Program. Each year, all the PI's present the technical advancements of their projects during a Grouped Technical Session.

This study assimilated Sentinel-1 C-band backscatter observations over snow-covered terrain into the Noah-Multiparameterization land surface model using support vector machine (SVM) regression and an ensemble Kalman filter to improve the modeled terrestrial snow mass estimates. The data assimilation (DA) experiment was conducted across Western Colorado from September 2016 to August 2017. As part of the DA experiments, the impact of a rule-based update was evaluated by comparing snow water equivalent (SWE) estimates via DA (with [ DAv1 ] and without [ DAv2 ] the rule-based update) against SNOTEL SWE measurements. Results confirmed that rule-based update helped minimize SVM controllability issues, and in turn, improved the accuracy of SWE estimates relative to both open loop (OL) and DAv2 . Comparison of SWE estimates from Sentinel-1 DAv1 against SNOTEL SWE revealed that 75% of stations showed improvements in bias and correlation coefficient relative to the OL. Assimilated SWE estimates also showed statistical improvements during both the snow accumulation and snow ablation periods. However, unbiased root mean square error showed a slight increase during the snow ablation period due to the large variability in the electromagnetic response of C-band backscatter over deep and/or wet snow. Improvement of the SWE estimates also resulted in improving river discharge estimates compared to in situ measurements. River discharge using Sentinel-1 DAv1 improved the Nash–Sutcliffe efficiency at all available stations. These results suggest that physically constrained SVM can serve as an efficient observation operator for snow mass DA through explicit consideration of the first-order C-band scattering mechanisms over different terrestrial snow conditions.

Identification of flood water extent from satellite images has historically relied on either synthetic aperture radar (SAR) or multi-spectral (MS) imagery. MS sensors are limited to cloud free conditions, whereas SAR imagery is plagued by noise-like speckle. Prior studies that use combinations of MS and SAR data to overcome individual limitations of these sensors have not fully examined sensitivity of flood mapping performance to different combinations of SAR and MS derived spectral indices or band transformations in color space. This study explores the use of diverse bands of Sentinel 2 (S2) through well-established water indices and Sentinel 1 (S1) derived SAR imagery along with their combinations to assess their capability for generating accurate flood inundation maps. The robustness in performance of S-1 and S-2 band combinations was evaluated using 446 hand labeled flood inundation images spanning across 11 flood events from Sen1Floods11 dataset which are highly diverse in terms of land cover as well as location. A modified K-fold cross validation approach is used to evaluate the performance of 32 combinations of S1 and S2 bands using a fully connected deep convolutional neural network known as U-Net. Our results indicated that usage of elevation information has improved the capability of S1 imagery to produce more accurate flood inundation maps. Compared to a median F1 score of 0.62 when using only S1 bands, the combined use of S1 and elevation information led to an improved median F1 score of 0.73. Water extraction indices based on S2 bands have a statistically significant superior performance in comparison to S1. Among all the band combinations, HSV (Hue, Saturation, Value) transformation of S2 bands provides a median F1 score of 0.9, outperforming the commonly used water spectral indices owing to HSV’s transformation’s superior contrast distinguishing abilities. Additionally, U-Net algorithm was able to learn the relationship between raw S2 based water extraction indices and their corresponding raw S2 bands, but not of HSV owing to relatively complex computation involved in the latter. Results of the paper establishes important benchmarks for the extension of S1 and S2 data-based flood inundation mapping efforts over large spatial extents.

Accurately representing land-atmosphere (L-A) interactions and coupling in NWP systems remains a challenge. New observations, incorporated into models via assimilation or calibration, hold the promise of improved forecast skill, but erroneous model coupling can hinder the benefits of such activities. To better understand model representation of coupled interactions and feedbacks, this study demonstrates a novel framework for coupled calibration of the Single Column Model (SCM) capability of the NASA Unified Weather Research and Forecasting (NU-WRF) system coupled to NASA’s Land Information System (LIS). The local land-atmosphere coupling (LoCo) process chain paradigm is used to assess the processes and connections revealed by calibration experiments. Two summer case studies in the U. S. Southern Great Plains are simulated in which LSM parameters are calibrated to diurnal observations of LoCo process chain components including 2-meter temperature, 2-meter humidity, surface fluxes (Bowen ratio), and PBL height. Results show a wide range of soil moisture and hydraulic parameter solutions depending on which L-A variable (i.e. observation) is used for calibration, highlighting that improvement in either SHP or ISM when not in tandem with the other can provide undesirable results. Overall, this work demonstrates that a process chain calibration approach can be used to assess L-A connections, feedbacks, strengths, and deficiencies in coupled models, as well as quantify the potential impact of new sources of observations of land-PBL variables on coupled prediction.

This study quantifies the level of observational accuracy required from a spaceborne light detection and ranging (LiDAR) snow depth retrieval mission for enabling beneficial impacts for snow estimation. The study is conducted over a region in Western Colorado using a suite of observing system simulation experiments(OSSEs).The Joint UK Land Environment Simulator, version 5.0 (JULES v5.0) is employed to simulate a suite of idealized LiDAR observations, considering a range of LiDAR snow depth retrieval errors, different hypothetical sensor swath widths, and the impact of cloud cover on observability. These simulated observations are then assimilated into the Noah land surface modelwith multi-parameterization options, version 3.6 (Noah-MP v3.6) model. This data assimilation setup is used to systematically evaluate the potential utility of LiDAR observations for improving modeled snow water equivalent (SWE)estimates and water budget variables such as runoff. Results from the OSSE runsshow that, in general, assimilation of synthetic LiDAR observations provide beneficial impacts when theLiDAR snow depth retrieval error standard deviation (σerror) is below 60 cm.Based on comparisons between the realistic (i.e., swath-limited and cloud-attenuated) case and the idealized (i.e., infinite swath width in the absence of cloud cover) case,this study concludes that observations with a conservative error standard deviation threshold of 40 cm (i.e., upper limit of the snow depth retrieval error that adds value to the SWE estimates viaassimilation) are needed for improving modeled snow estimates. More than a 33% reduction in SWEroot mean square errors and more than a 15% increase in correlation coefficientsare achieved when σ error ≤ 40 cmusing a 170-km sensor swath width in the presence of cloud attenuation effects. Further, the integrated hydrologic response, as represented by total (surface and subsurface) runoff estimates during the snow ablation season, are also enhanced when assimilating synthetic LiDAR snow depth retrievals with errors below this level.

In this study, we tested the impact of a revised set of soil, vegetation and land cover parameters on the performance of three different state-of-the-art land surface models (LSMs) within the NASA Land Information System (LIS). The impact of this revision was tested over the South American Dry Chaco, an ecoregion characterized by deforestation and forest degradation since the 1980s. Most large-scale LSMs may lack the ability to correctly represent the ongoing deforestation processes in this region, because most LSMs use climatological vegetation indices and static land cover information. The default LIS parameters were revised with improved soil parameters, satellite-based interannually varying vegetation indices (leaf area index and green vegetation fraction) instead of climatological vegetation indices, and yearly land cover information instead of static land cover. A relative comparison in terms of water budget components and “efficiency space” for various baseline and revised experiments showed that large regional and long-term differences in the simulated water budget partitioning relate to different LSM structures, whereas smaller local differences resulted from updated soil, vegetation and land cover parameters. Furthermore, the different LSM structures redistributed water differently in response to these parameter updates. A time-series comparison of the simulations to independent satellite-based estimates of evapotranspiration and brightness temperature (Tb) showed that no LSM setup significantly outperformed another for the entire region and that not all LSM simulations improved with updated parameter values. However, the revised soil parameters generally reduced the bias between simulated surface soil moisture and pixel-scale in situ observations and the bias between simulated Tb and regional Soil Moisture Ocean Salinity (SMOS) observations. Our results suggest that the different hydrological responses of various LSMs to vegetation changes may need further attention to gain benefits from vegetation data assimilation.

Though coarse in spatial resolution, the nearly all weather measurements from passive microwave sensors can help in improving the spatiotemporal coverage of optical and thermal infrared sensors for monitoring vegetation changes on the land surface. This study demonstrates the use of vegetation optical depth retrievals from the Soil Moisture Active Passive mission for capturing the vegetation alterations from the recent 2019-2020 Australian bushfires and drought. The impact of vegetation disturbance s on terrestrial water budget is examined by assimilating the vegetation optical depth retrievals into a dynamic phenology model. The results demonstrate that assimilating vegetation optical depth observations lead to improved simulation of evapotranspiration, runoff, and soil moisture states. The study also demonstrates that the vegetation changes from the 2019-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 study presents an evaluation of the impact of vegetation conditions on a land-surface model (LSM) simulation of agricultural drought. The Noah-MP LSM is used to simulate water and energy fluxes and states, which are transformed into drought categories using percentiles over the continental U.S. from 1979 to 2017. Leaf Area Index (LAI) observations are assimilated into the dynamic vegetation scheme of Noah-MP. A weekly operational drought monitor (the U.S. Drought Monitor) is used for the evaluation. The results show that LAI assimilation into Noah-MP’s dynamic vegetation scheme improves the model’s ability to represent drought, particularly over cropland areas. LAI assimilation improves the simulation of the drought category, detection of drought conditions, and reduces the instances of drought false alarms. The assimilation of LAI in these locations not only corrects model errors in the simulation of vegetation, but also can help to represent unmodeled physical processes such as irrigation towards improved simulation of agricultural drought.

This first paper of the two-part series focuses on demonstrating the accuracy of a hyper-resolution, offline terrestrial modeling system used for the High Mountain Asia (HMA) region. To this end, this study systematically evaluates four sets of model simulations at point scale, basin scale, and domain scale obtained from different spatial resolutions including 0.01° (∼1-km) and 0.25° (∼25-km). The assessment is conducted via comparisons against ground-based observations and satellite-derived reference products. The key variables of interest include surface net shortwave radiation, surface net longwave radiation, skin temperature, near-surface soil temperature, snow depth, snow water equivalent, and total runoff. In the evaluation against ground-based measurements, the superiority of the 0.01° estimates are mostly demonstrated across relatively complex terrain. Specifically, hyper-resolution modeling improves the skill in meteorological forcing estimates (except precipitation) by 9% relative to coarse-resolution estimates. The model forced by downscaled forcings in its entirety yields the highest skill in model output states as well as precipitation, which improves the skill obtained by coarse-resolution estimates by 7%. These findings, on one hand, corroborate the importance of employing the hyper-resolution versus coarse-resolution modeling in areas characterized by complex terrain. On the other hand, by evaluating four sets of model simulations forced with different precipitation products, this study emphasizes the importance of accurate hyper-resolution precipitation products to drive model simulations.

The Western Land Data Assimilation System (WLDAS) is a custom instance of the NASA Land Information System that combines land surface parameters, meteorological forcing data, and satellite products within a land surface model to produce daily estimates of the water and energy budget variables for the western United States. WLDAS was configured through discussions with partners, with the goal of groundwater sustainability planning for the state of California in mind. The publicly available output dataset has a 1-km grid resolution and spans 1979–present, which makes it suitable for water resources assessments. The data are also able to contextualize the recent drought events in California. Assimilation of Leaf Area Index, which is demonstrated herein to improve simulation over agricultural areas in California, specifically in terms of evapotranspiration in irrigation regions, will be included along with other data assimilation in subsequent releases of WLDAS.

South and Southeast Asia is subject to significant hydrometeorological extremes, including drought. Under rising temperatures, growing populations, and an apparent weakening of the South Asian monsoon in recent decades, concerns regarding drought and its potential impacts on water and food security are on the rise. Reliable sub-seasonal to seasonal (S2S) hydrological forecasts could, in principle, help governments and international organizations to better assess risk and act in the face of an oncoming drought. Here, we leverage recent improvements in S2S meteorological forecasts and the growing power of Earth Observations to provide more accurate monitoring of hydrological states for forecast initialization. Information from both sources is merged in a South and Southeast Asia sub-seasonal to seasonal hydrological forecasting system (SAHFS-S2S), developed collaboratively with the NASA SERVIR program and end-users across the region. This system applies the Noah-MultiParameterization (NoahMP) Land Surface Model (LSM) in the NASA Land Information System (LIS), driven by downscaled meteorological fields from the Global Data Assimilation System (GDAS) and Climate Hazards InfraRed Precipitation products (CHIRP and CHIRPS) to optimize initial conditions. The NASA Goddard Earth Observing System Model - sub-seasonal to seasonal (GEOS-S2S) forecasts, downscaled using the National Center for Atmospheric Research (NCAR) General Analog Regression Downscaling (GARD) tool and quantile mapping, are then applied to drive 5-km resolution hydrological forecasts to a 9-month forecast time horizon. Results show that the skillful predictions of root zone soil moisture can be made one to two months in advance for forecasts initialized in rainy seasons and up to 8 months when initialized in dry seasons. The memory of accurate initial conditions can positively contribute to forecast skills throughout the entire 9-month prediction period in areas with limited precipitation. This SAHFS-S2S has been operationalized at the International Centre for Integrated Mountain Development (ICIMOD) to support drought monitoring and warning needs in the region.

The model configuration employs the modified International Geosphere Biosphere Programme (IGBP) MODIS 20 category landcover data, soil parameters derived from the International Soil Reference and Information Centre, and the Shuttle Radar Topography Mission based elevation, slope, and aspect data. Statistical downscaling approaches are used to transform the coarse resolution MERRA2 meteorological inputs to 1km. The input meteorological fields of air temperature, humidity, surface pressure, wind, downward shortwave radiation, and downward longwave radiation are downscaled to 1km by adjusting for terrain differences in elevation, slope, and aspect. The high resolution monthly precipitation climatology from WorldClim is used to spatially disaggregate input MERRA2 precipitation to 1km. The initial conditions for the model simulations are generated from a long spinup of NoahMP starting from year 2000. All model integrations and evaluations are conducted using the NASA Land Information System and the Land surface Verification Toolkit.

Climate variability is an important driver of irrigation water use in many regions. Efforts to anticipate climate change impacts on future water availability can benefit from understanding how irrigation water demand has responded to these drivers to date. Here we apply satellite derived data, meteorological reanalysis, an advanced land surface model, and available state level reports to quantify irrigation demand sensitivities to temperature and precipitation across the Contiguous United States, for the period of 2002-2017. As expected, strong negative correlations are found between precipitation and irrigation withdrawals, both simulated and reported. Temperature sensitivities, however, vary by region and season, as do the interactive effects of temperature and precipitation on irrigation. Climate-induced irrigation variability is largest in transitional climate zones. These transitional zones are generally separate from the regions where rates of irrigation withdrawals are greatest, such that climate-induced variability in irrigation demand represents a water resource consideration that is distinct from chronic over pumping.

The utility of hydrologic land surface models (LSMs) can be enhanced by using information from observational platforms, but mismatches between the two are common.This study assesses the degree to which model agreement with observations (observability) is affected by two mechanisms in particular: 1) physical incongruities between the support volumes being characterized and 2) inadequateor inconsistent parameterizations of physical processes. The Noah and Noah-MP LSMs by default characterize surface soil moisture (SSM) in the top 10 cm of the soil column. This depth may be reasonable when comparing against soil moisture from in situ probes centered at 5 cm,but it is notably different from the 5 cm (or less) sensing depth of NASA’s Soil Moisture Active Passive (SMAP) satellite mission. These depth inconsistencies are examined by using thinner model layers in the Noah and Noah-MP LSMs and comparing resultant simulations to in situ and SMAP soil moisture. In addition, a forward radiative transfer model to simulate microwave brightness temperatures (Tbs) is used to facilitate direct comparisons of LSM-based and SMAP-based L-band Tb retrievals. Observability is quantified using Kolmogorov-Smirnov distance values, calculated from empirical cumulative distribution functions of SSM and Tb time series. Experiment results depend on the particular subspace being analyzed (SSM or Tb). This study concludes that therole of increasedsoil layer discretizationon LSM observability is secondary to the influence of component parameterizations, the effects of which dominate systematic differences with observations

Traditional watershed modeling often overlooks the role of vegetation dynamics. There is also little quantitative evidence to suggest that increased physical realism of vegetation dynamics in process-based models improves hydrology and water quality predictions simultaneously. In this study, we applied a modified Soil and Water Assessment Tool (SWAT) to quantify the extent of improvements that the assimilation of remotely sensed Leaf Area Index (LAI) would convey to streamflow, soil moisture, and nitrate load simulations across a 16,860 km2 agricultural watershedin the midwestern United States. We modified the SWAT source code to automatically override the model’s built-in semiempirical LAI with spatially distributed and temporally continuous estimates from Moderate Resolution Imaging Spectroradiometer (MODIS). Compared to a “basic” traditional model with limited spatial information, our LAI assimilation model (i) significantly improved daily streamflow simulations during medium-to-low flow conditions, (ii) provided realistic spatial distributions of growing season soil moisture, and (iii) substantially reproduced the long-term observed variability of daily nitrate loads. Further analysis revealed that the overestimation or underestimation of LAI imparted a proportional cascading effect on how the model partitions hydrologic fluxes and nutrient pools. As such, assimilation of MODIS LAI data corrected the model’sLAI overestimation tendency, which led to a proportionally increased rootzone soil moisture and decreased plant nitrogen uptake. With these new findings, our study fills the existing knowledge gap regarding vegetation dynamics in watershed modeling and confirms that assimilation of MODIS LAI data in watershed models can effectively improve both hydrology and water quality predictions.

The Korea Land Data Assimilation System (KLDAS) has been established for agricultural drought (i.e. soil moisture deficit) monitoring in South Korea, running the Noah-MP land surface model within the NASA Land Information System (LIS) framework with the added value of local precipitation forcing dataset and soil texture maps. KLDAS soil moisture is benchmarked against three global products: the Global Land Data Assimilation System (GLDAS), the Famine Early Warning Systems Network (FEWS NET) Land Data Assimilation System (FLDAS), and the European Space Agency Climate Change Initiative (ESA CCI) satellite product. The evaluation is performed using in situ measurements for 2013–2015 and one month standardized precipitation index (SPI-1) for 1982–2016, focusing on four major river basins in South Korea. The KLDAS outperforms all benchmark products in capturing soil moisture states and variability at a basin scale. Compared to GLDAS and FLDAS products, the EAS CCI product is not feasible for long term agricultural monitoring due to lower data quality for early periods (1979–1991) of soil moisture estimates. KLDAS shows that the most recent 2015 drought event leads to highest drought areas in the Han and Geum River basins in the past 35 years. This work supports KLDAS as an effective agricultural drought monitoring system to provide continuous regional high-resolution soil moisture estimates in South Korea.

Vegetation optical depth (VOD) retrievals from passive microwave sensors provide analog estimates of above-ground canopy biomass. This study presents the development and analysis of assimilating VOD retrievals from X-, C-, and L-band passive microwave instruments within the Noah-MP land surface model over the Continental U.S. The results from this study demonstrate that the assimilation of VOD retrievals have a significant beneficial impact on the simulation of evapotranspiration and GPP, particularly over the agricultural areas of the U.S. The improvements in the water and carbon fluxes from the assimilation of VOD from X- and C-band sensors are found to be comparable to those obtained from the assimilation of vegetation indices from optical sensors.The study also quantifies the relative and joint impact of assimilating surface soil moisture and VOD from the Soil Moisture Active Passive (SMAP) mission. The utility of soil moisture assimilation for improving ET is more significant over water-limited regions, whereas VOD DA is more impactful over areas where soil moisture is not the primary controlling factor on ET. The results also indicate that the information on moisture and vegetation states from SMAP can be simultaneously exploited through the joint assimilation of surface soil moisture and VOD. Since passive microwave-based VOD retrievals are available in nearly all weather conditions, their use within data assimilation systems offers the ability to extend and improve the utility obtained from the use of optical/infrared based vegetation retrievals