In drainage basin evolution models the implementation of sediment transport and sediment balance is constrained by numerical instabilities. In order to overcome these difficulties an implicit finite difference method is proposed. The method simultaneously adjusts the elevation at each node of the numerical grid, is unconditionally stable, and significantly reduces the computational time. The performance of the method is tested and briefly discussed with a numerical example.