Steady- and transient-state inversion in hydrogeology by successive flux estimation

A calibration method to solve the groundwater inverse problem under steady- and transient-state conditions is presented. The method compares kriged and numerical head field gradients to modify hydraulic conductivity without the use of non-linear optimization techniques. The process is repeated iteratively until a close match with piezometric data is reached. The approach includes a damping factor to avoid divergence and oscillation of the solution in areas of low hydraulic gradient and a weighting factor to account for temporal head variation in transient simulations. The efficiency of the method in terms of computing time and calibration results is demonstrated with a synthetic field. It is shown that the proposed method provides parameter fields that reproduce both hydraulic conductivity and piezometric data in few forward model solutions. Stochastic numerical experiments are conducted to evaluate the sensitivity of the method to the damping function and to the head field estimation errors.