The prediction of water table height in unconfined layered porous media is a difficult modelling problem that typically requires numerical simulation. This paper proposes an analytical model to approximate the exact solution based on a steady-state Dupuit–Forchheimer analysis. The key contribution in relation to a similar model in the literature relies in the ability of the proposed model to consider more than two layers with different thicknesses and slopes, so that the existing model becomes a special case of the proposed model herein. In addition, a model assessment methodology based on the Bayesian inverse problem is proposed to efficiently identify the values of the physical parameters for which the proposed model is accurate when compared against a reference model given by MODFLOW-NWT, the open-source finite-difference code by the U.S. Geological Survey. Based on numerical results for a…