Abstract
The inverse problem is formulated as a linear least squares problem. Lower and upper bounds on the parameters are added to ensure hydrologic plausibility. The proposed algorithm utilizes the special structure of the problem to enable efficient handling of models which are composed of a large number of cells. An example of calibrating a 241‐cell model of the coastal aquifer of Israel with the proposed algorithm is presented.