It is challenging to incorporate input variables' spatial distribution information when implementing long short-term memory (LSTM) models for streamflow prediction. This work presents a novel hybrid modelling approach to predict streamflow while accounting for spatial variability. We evaluated the performance against lumped LSTM predictions in 224 basins across the Great Lakes region in North America. This approach shows promise for predicting streamflow in large, ungauged basin.