A regional-scale land surface parameterization based on areally-averaged hydrological conservation equations