Multi-objective two-stage optimization of canal system scheduling in irrigation districts incorporating groundwater dynamic variations