Extended mild-slope equations for the propagation of small-amplitude water waves over variable bathymetry regions, recently proposed by Massel (1993) and Porter & Staziker (1995), are shown to exhibit an inconsistency concerning the sloping-bottom boundary condition, which renders them non-conservative with respect to wave energy. In the present work, a consistent coupled-mode theory is derived from a variational formulation of the complete linear problem, by representing the vertical distribution of the wave potential as a uniformly convergent series of local vertical modes at each horizontal position. This series consists of the vertical eigenfunctions associated with the propagating and all evanescent modes and, when the slope of the bottom is different from zero, an additional mode, carrying information about the bottom slope. The coupled-mode system obtained in this way contains an additional equation, as well as additional interaction terms in all other equations, and reduces to the previous extended mild-slope equations when the additional mode is neglected. Extensive numerical results demonstrate that the present model leads to the exact satisfaction of the bottom boundary condition and, thus, it is energy conservative. Moreover, it is numerically shown that the rate of decay of the modal-amplitude functions is improved from O(n−2), where n is the mode number, to O(n−4), when the additional sloping-bottom mode is included in the representation. This fact substantially accelerates the convergence of the modal series and ensures the uniform convergence of the velocity field up to and including the boundaries.