“…Establishing an additional buffer distance for the static attributes, including the grid cells outside each model grid, may allow us to consider such effects (Xu et al, 2023). Also, future improvements on the CNN-based integrated model can involve novel structures like ConvLSTM (Li et al, 2021;Shi et al, 2015) that can enhance the utilization of the spatio-temporally varying predictors (e.g., dynamic landscape features), adopting more advanced loss functions like Kling-Gupta efficiency (Gupta et al, 2009), and applying multi-objective optimization (e.g., soil moisture, terrestrial water storage, and runoff) to further avoid the issue of parameter identifiability (Beven, 2006;Yi & Park, 2021). Our approach can be flexibly transferred to other regions, improving the regionalization of the model and increasing our understanding of the natural environment.…”