Parameter estimation and uncertainty analysis of the Spatial Agro Hydro Salinity Model (SAHYSMOD) in the semi-arid climate of Rechna Doab,Pakistan

Manual calibration of distributed models with many unknown parameters can result in problems of equifinality and high uncertainty. In this study, the Generalized Likelihood Uncertainty Estimation (GLUE) technique was used to address these issues through uncertainty and sensitivity analysis of a distributed watershed scale model (SAHYSMOD) for predicting changes in the groundwater levels of the Rechna Doab basin, Pakistan. The study proposes and then describes a stepwise methodology for SAHYSMOD uncertainty analysis that has not been explored in any study before. One thousand input data files created through Monte Carlo simulations were classified as behavior and non-behavior sets using threshold likelihood values. The model was calibrated (1983–1988) and validated (1998–2003) through satisfactory agreement between simulated and observed data. Acceptable values were observed in the statistical performance indices. Approximately 70% of the observed groundwater level values fell within uncertainty bounds. Groundwater pumping (Gw) and hydraulic conductivity (Kaq) were found to be highly sensitive parameters affecting groundwater recharge.     

Optimization of the shape (and topology) of the initial conditions for diffusion parameter identification

The design of an experiment, e.g., the setting of initial conditions, strongly influences the accuracy of the whole process of determining model parameters from data. We impose a sensitivity-based approach for choosing optimal experimental design variables and study the optimization of the shape (and topology) of the initial conditions for an inverse problem of a diffusion parameter identification. Our approach, although case independent, is illustrated at the FRAP (Fluorescence Recovery After Photobleaching) experimental technique. The core idea resides in the maximization of a sensitivity measure, which depends on a specific experimental setting of initial conditions. By a numerical optimization, considering radially symmetric shapes only, we find an interesting pattern of increasingly complicated (with respect to connectivity) optimal initial shapes. The proposed modification of the FRAP experimental protocol is rather surprising but entirely realistic and the resulting enhancement of the parameter estimate accuracy is significant.