The application of multi-mission satellite data assimilation for studying water storage changes over South America

Land Surface Model Calibration Using Satellite Remote Sensing Data

Satellite remote sensing provides a unique opportunity for calibrating land surface models due to their direct measurements of various hydrological variables as well as extensive spatial and temporal coverage. This study aims to apply terrestrial water storage (TWS) estimated from the gravity recovery and climate experiment (GRACE) mission as well as soil moisture products from advanced microwave scanning radiometer-earth observing system (AMSR-E) to calibrate a land surface model using multi-objective evolutionary algorithms. For this purpose, the non-dominated sorting genetic algorithm (NSGA) is used to improve the model's parameters. The calibration is carried out for the period of two years 2003 and 2010 (calibration period) in Australia, and the impact is further monitored over 2011 (forecasting period). A new combined objective function based on the observations' uncertainty is developed to efficiently improve the model parameters for a consistent and reliable forecasting skill. According to the evaluation of the results against independent measurements, it is found that the calibrated model parameters lead to better model simulations both in the calibration and forecasting period.

Spatiotemporal Downscaling of GRACE Total Water Storage Using Land Surface Model Outputs

High spatiotemporal resolution of terrestrial total water storage plays a key role in assessing trends and availability of water resources. This study presents a two-step method for downscaling GRACE-derived total water storage anomaly (GRACE TWSA) from its original coarse spatiotemporal resolution (monthly, 3-degree spherical cap/~300 km) to a high resolution (daily, 5 km) through combining land surface model (LSM) simulated high spatiotemporal resolution terrestrial water storage anomaly (LSM TWSA). In the first step, an iterative adjustment method based on the self-calibration variance-component model (SCVCM) is used to spatially downscale the monthly GRACE TWSA to the high spatial resolution of the LSM TWSA. In the second step, the spatially downscaled monthly GRACE TWSA is further downscaled to the daily temporal resolution. By applying the method to downscale the coarse resolution GRACE TWSA from the Jet Propulsion Laboratory (JPL) mascon solution with the daily high spatial resolution (5 km) LSM TWSA from the Ecological Assimilation of Land and Climate Observations (EALCO) model, we evaluated the benefit and effectiveness of the proposed method. The results show that the proposed method is capable to downscale GRACE TWSA spatiotemporally with reduced uncertainty. The downscaled GRACE TWSA are also evaluated through in-situ groundwater monitoring well observations and the results show a certain level agreement between the estimated and observed trends.

Multi-mission satellite remote sensing data for improving land hydrological models via data assimilation

Satellite remote sensing offers valuable tools to study Earth and hydrological processes and improve land surface models. This is essential to improve the quality of model predictions, which are affected by various factors such as erroneous input data, the uncertainty of model forcings, and parameter uncertainties. Abundant datasets from multi-mission satellite remote sensing during recent years have provided an opportunity to improve not only the model estimates but also model parameters through a parameter estimation process. This study utilises multiple datasets from satellite remote sensing including soil moisture from Soil Moisture and Ocean Salinity Mission and Advanced Microwave Scanning Radiometer Earth Observing System, terrestrial water storage from the Gravity Recovery And Climate Experiment, and leaf area index from Advanced Very-High-Resolution Radiometer to estimate model parameters. This is done using the recently proposed assimilation method, unsupervised weak constrained ensemble Kalman filter (UWCEnKF). UWCEnKF applies a dual scheme to separately update the state and parameters using two interactive EnKF filters followed by a water balance constraint enforcement. The performance of multivariate data assimilation is evaluated against various independent data over different time periods over two different basins including the Murray–Darling and Mississippi basins. Results indicate that simultaneous assimilation of multiple satellite products combined with parameter estimation strongly improves model predictions compared with single satellite products and/or state estimation alone. This improvement is achieved not only during the parameter estimation period (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim $$\end{document}∼ 32% groundwater RMSE reduction and soil moisture correlation increase from \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim $$\end{document}∼ 0.66 to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim $$\end{document}∼ 0.85) but also during the forecast period (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim $$\end{document}∼ 14% groundwater RMSE reduction and soil moisture correlation increase from \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim $$\end{document}∼ 0.69 to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim $$\end{document}∼ 0.78) due to the effective impacts of the approach on both state and parameters.

Altimetry-derived surface water data assimilation over the Nile Basin

Global hydrological models facilitate studying of water resources and their variations over time. The accuracies of these models are enhanced when combined with ever-increasing satellite remotely sensed data. Traditionally, these combinations are done via data assimilation approach, which permits the use of improved hydrological outputs to study regions with limited in-situ measurements such as the Nile Basin. This study aims at using the state-of-art satellite radar altimetry data to enhance a land-based hydrological model for studying water storage changes over the Nile Basin. Altimetry-derived surface water storage, for the first time, is assimilated into the model using the ensemble Kalman filter (EnKF) for the period of 2003 to 2016. Multiple datasets from ground measurements, as well as space observations, are used to evaluate the performance of the assimilated satellite altimetry data. Results indicate that the assimilation successfully improves model outputs, especially the surface water component. The process increases the correlation between surface water storage changes and water level variations from satellite radar altimetry by 0.44 and reduces the surface water discharge root-mean-square errors (RMSE) by approximately 33%.

Vertical difference of climate change impacts on vegetation at temporal-spatial scales in the upper stream of the Mekong River Basin

As the upper section of the Mekong River Basin, the vegetation quality of the Lancang River Basin (LRB) and the related ecological functions are critical for the whole basin. With time-series Moderate Resolution Imaging Spectroradiometer (MODIS) images from 2000 to 2015 and local daily climatic data since 1976, their vertical interaction differences were identified. The results showed that the spatial variation in Normalized difference vegetation index (NDVI) of grassland and forest were sensitive to elevation. The NDVI value in the southern area at elevations less than 3000 m was more than 0.80 and decreased to 0.30-0.60 with elevations higher than 4500 m. The general vegetation quality showed a positive trend under climate change over 16 years. The M-K test of daily precipitation and temperature from 12 local weather stations showed that the basin temperature varied more significantly than precipitation. The temporal correlation between NDVI with precipitation as well as temperature at each pixel indicated that temperature was the dominant factor affecting grassland and forest dynamics in the LRB. The interaction between vegetation and climate was more sensitive at elevations lower than 3000 m. Based on the RCP4.5 scenario, the future temperature distribution was predicted, and its impact on NDVI was simulated at the pixel scale. Under future drier and warmer climate conditions, the responded NDVI in the upper stream with higher elevation may increase soil erosion and decrease streamflow. The NDVI in the downstream area will be improved and be able to adapt to the related climate impacts. Because of the large amount of water and biomass in this basin, higher temperatures will accelerate the decomposition of forest foliar litter. Thus, more organic carbon and forest diffuse pollution will be discharged into the water, potentially affecting the water quality of the whole basin.

NCA-LDAS: Overview and Analysis of Hydrologic Trends for the National Climate Assessment

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°× 1258° 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 10.05 to 10.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.

Estimation of root zone soil moisture using passive microwave remote sensing: A case study for rice and wheat crops for three states in the Indo-Gangetic basin

This paper investigates estimation of root zone soil moisture using two passive microwave remote sensing datasets, Advanced Microwave Scattering Radiometer - 2 and Soil Moisture Active Passive satellites sensors. The study is focused on two crops, namely rice and wheat for the Indo-Gangetic basin, India, having a dynamic crop and soil type and land use land cover. A total of 21 rice crop and 23 wheat crop locations are chosen from the states of Uttar Pradesh, Madhya Pradesh and Bihar falling in the basin. The root zone soil moisture information is derived by estimating soil wetness index from surface soil moisture at 10 and 40 cm depths using a recursive exponential filter. The soil wetness index based algorithm is implementable even in the absence of ground information for a basin level study. The reference soil moisture dataset is obtained from the Global Land Database Assimilation System - NOAH at 10 and 40 cm depth. The research has also demonstrated significant potential of GLDAS-NOAH soil moisture data in the absence of ground (in-situ) soil moisture data. Of the various factors affecting surface and root zone soil moisture, this work evaluates the control of soil constituents on root zone soil moisture. The Spearman rank correlation coefficient is estimated for characteristic time delay with sand, silt and clay percentage at different locations. Coupling between and trends of surface and root zone soil moisture for rice and wheat crop locations are studied. The accuracy of estimated soil wetness index at 10 and 40 cm from two different satellite sensors at two different acquisition times (ascending and descending passes) is investigated by calculating the coefficient of determination, mean absolute error and mean biased error. This work highlights the significant difference in surface soil moisture estimation by two satellite sensors to derive root zone soil moisture for rice and wheat crops. Coefficient of determination is more (∼0.9) for SMAP derived soil wetness index whereas it is lower (∼0.65) for AMSR-2 derived soil wetness index for both crops. Characteristic time delay variation is observed at two different times and at both the depths, with characteristic time delay increasing with depth. Also, at the descending pass characteristic time delay is lower as compared to the ascending pass. A strong relationship between root zone soil moisture and soil texture is observed. For rice crop, a positive correlation with sand and clay is observed for Uttar Pradesh, Madhya Pradesh and Bihar locations having loam and sandy loam as the major soil class. And, for wheat locations, a positive correlation is observed for silt and clay for Uttar Pradesh locations and sand for Madhya Pradesh locations having loam and clay (light) soil texture. This work delivers essential information in understanding sustainable irrigation scheduling and increasing irrigation potential for rice and wheat crop locations. Having the knowledge of all the factors influencing crop cultivation and the derived root zone soil moisture, crop production can be optimized.