A massively parallel OpenFOAM solver for Richards Equation: towards mechanistic modelling of transfers at the watershed scale

We focus on the development and test of a massively parallel OpenFOAM solver for Richards Equation, towards mechanistic modelling of flow and transport phenomena at the watershed scale

An open source massively parallel solver for Richards equation: Mechanistic modelling of water fluxes at the watershed scale

In 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® 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

The Indian COSMOS Network (ICON): validating L-band remote sensing and modelled soil moisture data products

Availability of global satellite based Soil Moisture (SM) data has promoted the emergence of many applications in climate studies, agricultural water resource management and hydrology. In this context, validation of the global data set is of substance. Remote sensing measurements which are representative of an area covering 100 m2 to tens of km2 rarely match with in situ SM measurements at point scale due to scale difference. In this paper we present the new Indian Cosmic Ray Network (ICON) and compare it’s data with remotely sensed SM at different depths. ICON is the first network in India of the kind. It is operational since 2016 and consist of seven sites equipped with the COSMOS instrument. This instrument is based on the Cosmic Ray Neutron Probe (CRNP) technique which uses non-invasive neutron counts as a measure of soil moisture. It provides in situ measurements over an area with a radius of 150–250 m. This intermediate scale soil moisture is of interest for the validation of satellite SM. We compare the COSMOS derived soil moisture to surface soil moisture (SSM) and root zone soil moisture (RZSM) derived from SMOS, SMAP and GLDAS_Noah. The comparison with surface soil moisture products yield that the SMAP_L4_SSM showed best performance over all the sites with correlation (R) values ranging from 0.76 to 0.90. RZSM on the other hand from all products showed lesser performances. RZSM for GLDAS and SMAP_L4 products show that the results are better for the top layer R = 0.75 to 0.89 and 0.75 to 0.90 respectively than the deeper layers R = 0.26 to 0.92 and 0.6 to 0.8 respectively in all sites in India. The ICON network will be a useful tool for the calibration and validation activities for future SM missions like the NASA-ISRO Synthetic Aperture Radar (NISAR)

Modeling the vulnerability of an urban groundwater system due to the combined impacts of climate change and management scenarios

Climate change impact on a groundwater-dependent small urban town has been investigated in the semiarid hard rock aquifer in southern India. A distributed groundwater model was used to simulate the groundwater levels in the study region for the projected future rainfall (2012-32) obtained from a general circulation model (GCM) to estimate the impacts of climate change and management practices on groundwater system. Management practices were based on the human-induced changes on the urban infrastructure such as reduced recharge from the lakes, reduced recharge from water and wastewater utility due to an operational and functioning underground drainage system, and additional water extracted by the water utility for domestic purposes. An assessment of impacts on the groundwater levels was carried out by calibrating a groundwater model using comprehensive data gathered during the period 2008-11 and then simulating the future groundwater level changes using rainfall from six GCMs Institute of Numerical Mathematics Coupled Model, version 3.0 (INM-CM. 3.0); L'Institut Pierre-Simon Laplace Coupled Model, version 4 (IPSL-CM4); Model for Interdisciplinary Research on Climate, version 3.2 (MIROC3.2); ECHAM and the global Hamburg Ocean Primitive Equation (ECHO-G); Hadley Centre Coupled Model, version 3 (HadCM3); and Hadley Centre Global Environment Model, version 1 (HadGEM1)] that were found to show good correlation to the historical rainfall in the study area. The model results for the present condition indicate that the annual average discharge (sum of pumping and natural groundwater outflow) was marginally or moderately higher at various locations than the recharge and further the recharge is aided from the recharge from the lakes. Model simulations showed that groundwater levels were vulnerable to the GCM rainfall and a scenario of moderate reduction in recharge from lakes. Hence, it is important to sustain the induced recharge from lakes by ensuring that sufficient runoff water flows to these lakes

Groundwater Level Dynamics in Bengaluru City, India

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

Uncertainty in estimation of potential recharge in tropical and temperate catchments using a crop model and microwave remote sensing.

International audienc

Retrieval and Multi-scale Validation of Soil Moisture from Multi-temporal SAR Data in a Semi-Arid Tropical Region

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.06 m(3)/m(3) for the majority of plots. Comparison with the SMOS SM showed a good temporal behaviour with RMSE of approximately 0.05 m(3)/m(3) 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

Copula-based downscaling of coarse-scale soil moisture observations with implicit bias correction

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