This study explores the benefits of assimilating SMOS soil moisture retrievals for hydrologic modeling, with a focus on soil moisture and streamflow simulations in the Murray Darling Basin, Australia. In this basin, floods occur relatively frequently and initial catchment storage is known to be key to runoff generation. The land surface model is the Variable Infiltration Capacity (VIC) model. The model is calibrated using the available streamflow records of 169 gauge stations across the Murray Darling Basin. The VIC soil moisture forecast is sequentially updated with observations from the SMOS Level 3 CATDS (Centre Aval de Traitement des Données SMOS) soil moisture product using the Ensemble Kalman filter. The assimilation algorithm accounts for the spatial mismatch between the model (0.125°) and the SMOS observation (25 km) grids. Three widely-used methods for removing bias between model simulations and satellite observations of soil moisture are evaluated. These methods match the first, second and higher order moments of the soil moisture distributions, respectively. In this study, the first order bias correction, i.e. the rescaling of the long term mean, is the recommended method. Preserving the observational variability of the SMOS soil moisture data leads to improved soil moisture updates, particularly for dry and wet conditions, and enhances initial conditions for runoff generation. Second or higher order bias correction, which includes a rescaling of the variance, decreases the temporal variability of the assimilation results. In comparison with in situ measurements of OzNet, the assimilation with mean bias correction reduces the root mean square error (RMSE) of the modeled soil moisture from 0.058 m 3 /m 3 to 0.046 m 3 /m 3 and increases the correlation from 0.564 to 0.714. These improvements in antecedent wetness conditions further translate into improved predictions of associated water fluxes, particularly runoff peaks. In conclusion, the results of this study clearly demonstrate the merit of SMOS data assimilation for soil moisture and streamflow predictions at the large scale.

No abstract

The current study presents an algorithm to retrieve surface Soil Moisture (SM) from multi-temporal Synthetic Aperture Radar (SAR) data. The developed algorithm is based on the Cumulative Density Function (CDF) transformation of multi-temporal RADARSAT-2 backscatter coefficient (BC) to obtain relative SM values, and then converts relative SM values into absolute SM values using soil information. The algorithm is tested in a semi-arid tropical region in South India using 30 satellite images of RADARSAT-2, SMOS L2 SM products, and 1262 SM field measurements in 50 plots spanning over 4 years. The validation with the field data showed the ability of the developed algorithm to retrieve SM with RMSE ranging from 0.02 to 0. and a correlation coefficient of approximately 0.9. The developed model is compared and found to be better than the change detection and delta index model. The approach does not require calibration of any parameter to obtain relative SM and hence can easily be extended to any region having time series of SAR data available.Remote Sens. 2015, 7 8129

Soil moisture retrievals, delivered as a CATDS (Centre Aval de Traitement des Données SMOS) Level-3 product of the Soil Moisture and Ocean Salinity (SMOS) mission, form an important information source, particularly for updating land surface models. However, the coarse resolution of the SMOS product requires additional treatment if it is to be used in applications at higher resolutions. Furthermore, the remotely sensed soil moisture often does not reflect the climatology of the soil moisture predictions, and the bias between model predictions and observations needs to be removed. In this paper, a statistical framework is presented that allows for the downscaling of the coarse-scale SMOS soil moisture product to a finer resolution. This framework describes the interscale relationship between SMOS observations and model-predicted soil moisture values, in this case, using the variable infiltration capacity (VIC) model, using a copula. Through conditioning, the copula to a SMOS observation, a probability distribution function is obtained that reflects the expected distribution function of VIC soil moisture for the given SMOS observation. This distribution function is then used in a cumulative distribution function matching procedure to obtain an unbiased fine-scale soil moisture map that can be assimilated into VIC. The methodology is applied to SMOS observations over the Upper Mississippi River basin. Although the focus in this paper is on data assimilation apcations, the framework developed could also be used for other purposes where downscaling of coarse-scale observations is required.

Abstract:Groundwater accounts for half of Indian urban water use. However, little is known about its sustainability, because of inadequate monitoring and evaluation. We deployed a dense monitoring network in 154 locations in Bengaluru, India between 2015 and 2017. Groundwater levels collected at these locations were analyzed to understand the behavior of the city's groundwater system. At a local scale, groundwater behavior is non-classical, with valleys showing deeper groundwater than ridge-tops. We hypothesize that this is due to relatively less pumping compared to artificial recharge from leaking pipes and wastewater in the higher, city core areas, than in the rapidly growing, lower peripheral areas, where the converse is true. In the drought year of 2016, groundwater depletion was estimated at 27 mm, or 19 Mm 3 over the study area. The data show that rainfall has the potential to replenish the aquifer. High rainfall during August-September 2017 led to a mean recharge of 67 mm, or 47 Mm 3 for the study area. A rainfall recharge factor of 13.5% was estimated from the data for 2016. Sustainable groundwater management in Bengaluru must account for substantial spatial socio-hydrological heterogeneity. Continuous monitoring at high spatial density will be needed to inform evidence-based policy.

⃝ Massively parallel computation Transfers in soils Porous media a b s t r a c tIn this paper we present a massively parallel open source solver for Richards equation, named the RichardsFOAM solver. This solver has been developed in the framework of the open source generalist computational fluid dynamics tool box OpenFOAM R ⃝ and is capable to deal with large scale problems in both space and time. The source code for RichardsFOAM may be downloaded from the CPC program library website.It exhibits good parallel performances (up to ∼90% parallel efficiency with 1024 processors both in strong and weak scaling), and the conditions required for obtaining such performances are analysed and discussed. These performances enable the mechanistic modelling of water fluxes at the scale of experimental watersheds (up to few square kilometres of surface area), and on time scales of decades to a century. Such a solver can be useful in various applications, such as environmental engineering for long term transport of pollutants in soils, water engineering for assessing the impact of land settlement on water resources, or in the study of weathering processes on the watersheds. Program summaryProgram title: RichardsFOAM Catalogue identifier: AEUF_v1_0 Program summary URL: L. Orgogozo et al. / Computer Physics Communications ( ) -External routines: OpenFOAM R ⃝ (version 2.0.1 or later) Nature of problem:This software solves the non-linear three-dimensional transient Richards equation, which is a very popular model for water transfer in variably saturated porous media (e.g.: soils). It is designed to take advantage of the massively parallel computing performance of OpenFOAM R ⃝ . The goal is to be able to model natural hydrosystems on large temporal and spatial scales. Solution method:A mixed implicit (FVM in the object oriented OpenFOAM R ⃝ framework) and explicit (FVC in the object oriented OpenFOAM R ⃝ framework) discretization of the equation with a backward time scheme is coupled with a linearization method (Picard algorithm). Due to the linearization loop the final solution of each time step tends towards a fully implicit solution. The implementation has been carried out with a concern for robustness and parallel efficiency. Restrictions:The choice of the maximum and initial time steps must be made carefully in order to avoid stability problems. A careful convergence study of mesh cell size, linear solver precision and linearization method precision must be undertaken for each considered problem, depending on the precision required for the expected results, the spatial and temporal scales at stake, and so on. Finally, the solver in its current version only handles meshes with a constant cell volume (a crash will not necessarily occur with an irregular mesh but some problems may arise with the convergence criterion of the linearization method). Running time:Highly variable, depending on the mesh size and the number and nature of cores involved. The test run provided requires less than 2 s on a 64 bit machine with Intel R ⃝ CoreT...

The main objective of this study is to analyze the potential use of L-band radar data for the estimation of soil moisture over tropical agricultural areas under dense vegetation cover conditions. Ten radar images were acquired using the Phased Array Synthetic Aperture Radar/Advanced Land Observing Satellite (PALSAR/ALOS)-2 sensor over the Berambadi watershed (south India), between June and October of 2018. Simultaneous ground measurements of soil moisture, soil roughness, and leaf area index (LAI) were also recorded. The sensitivity of PALSAR observations to variations in soil moisture has been reported by several authors, and is confirmed in the present study, even for the case of very dense crops. The radar signals are simulated using five different radar backscattering models (physical and semi-empirical), over bare soil, and over areas with various types of crop cover (turmeric, marigold, and sorghum). When the semi-empirical water cloud model (WCM) is parameterized as a function of the LAI, to account for the vegetation’s contribution to the backscattered signal, it can provide relatively accurate estimations of soil moisture in turmeric and marigold fields, but has certain limitations when applied to sorghum fields. Observed limitations highlight the need to expand the analysis beyond the LAI by including additional vegetation parameters in order to take into account volume scattering in the L-band backscattered radar signal for accurate soil moisture estimation.

Characterization of the soil water reservoir is critical for understanding the interactions between crops and their environment and the impacts of land use and environmental changes on the hydrology of agricultural catchments especially in tropical context. Recent studies have shown that inversion of crop models is a powerful tool for retrieving information on root zone properties. Increasing availability of remotely sensed soil and vegetation observations makes it well suited for large scale applications. The potential of this methodology has however never been properly evaluated on extensive experimental datasets and previous studies suggested that the quality of estimation of soil hydraulic properties may vary depending on agro-environmental situations. The objective of this study was to evaluate this approach on an extensive field experiment. The dataset covered four crops (sunflower, sorghum, turmeric, maize) grown on different soils and several years in South India. The components of AWC (available water capacity) namely soil water content at field capacity and wilting point, and soil depth of two-layered soils were estimated by inversion of the crop model STICS with the GLUE (generalized likelihood uncertainty estimation) approach using observations of surface soil moisture (SSM; typically from 0 to 10 cm deep) and leaf area index (LAI), which are attainable from radar remote sensing in tropical regions with frequent cloudy conditions. The results showed that the quality of parameter estimation largely depends on the hydric regime and its interaction with crop type. A mean relative absolute error of 5% for field capacity of surface layer, 10% for field capacity of root zone, 15% for wilting point of surface layer and root zone, and 20% for soil depth can be obtained in favorable conditions. A few observations of SSM (during wet and dry soil moisture periods) and LAI (within water stress periods) were sufficient to significantly improve the estimation of AWC components. These results show the potential of crop model inversion for estimating the AWC components of two-layered soils and may guide the sampling of representative years and fields to use this technique for mapping soil properties that are relevant for distributed hydrological modelling