The Water Residence Time distribution is the equivalent of the impulse response of a linear system allowing the propagation of water through a medium, e.g. the propagation of rain water from the top of the mountain towards the aquifers. We consider the output aquifer levels as the convolution between the input rain levels and the Water Residence Time, starting with an initial aquifer base level. The estimation of Water Residence Time is important for a better understanding of hydro-bio-geochemical processes and mixing properties of wetlands used as filters in ecological applications, as well as protecting fresh water sources for wells from pollutants. Common methods of estimating the Water Residence Time focus on cross-correlation, parameter fitting and non-parametric deconvolution methods. Here we propose a 1D full-deconvolution, regularized, non-parametric inverse problem algorithm that enforces smoothness and uses constraints of causality and positivity to estimate the Water Residence Time curve. Compared to Bayesian non-parametric deconvolution approaches, it has a fast runtime per test case; compared to the popular and fast cross-correlation method, it produces a more precise Water Residence Time curve even in the case of noisy measurements. The algorithm needs only one regularization parameter to balance between smoothness of the Water Residence Time and accuracy of the reconstruction. We propose an approach on how to automatically find a suitable value of the regularization parameter from the input data only. Tests on real data illustrate the potential of this method to analyze hydrological datasets.
The primary accomplishment of this research has been demonstrating the power of the deconvolution technique developed by Neuman and de Marsily (1976) in dealing with noisy rainfallrunoff records of short duration. Such records are encountered in arid environments where rainfall often occurs in short isolated bursts and the data are measured with a considerable margin of error. Our research work consisted of superimposing known noise on synthetic rainfallrunoff data and examining the ability of the Neuman -de Marsily deconvolution method to estimate the correct impulse response of the system when the data include only a single storm event. Approximately 50 Monte Carlo simulation runs were performed for each of three different noise models considered in our work. The results clearly demonstrated that the deconvolution model leads to reliable estimates and can be used with confidence in the presence of realistic noise levels. In addition to the Monte Carlo simulation tests and their analysis, certain improvements were introduced into the original deconvolution technique. In particular, the original version of the technique required that the hydrologist exercise subjective judgement in choosing the "best" solution for the deconvolution problem from a large number of admissible solutions. Our new method of selecting the "best" result is based on a comparative analysis of residuals and is more reliable than the earlier subjective approach. The improved method has been applied to real as well as synthetic rainfall runoff data.