Fast and robust solution methods for the water quality equations