Incorporating receiving waters responses into the framework of spatial optimization of LID-BMPs in plain river network region

•A NSGA-III-based LID-BMPs spatial optimization framework was proposed.•Receiving waters responses were incorporated in the framework.•LSTM model was utilized as a surrogate model to accelerate model simulation.•Significant cost-effective optimal solutions were delivered.•Imperviousness drives the process of spatial optimization.