Accurate prediction of river discharge is critical for a wide range of sectors, from human activities to environmental hazard management, especially in the face of increasing demand for water resources and climate change. To address this need, a multivariate model that incorporates both local and global data sources, including river and piezometer gauges, sea level, and climate parameters. By employing phase shift analysis, the model optimizes correlations between the target discharge and 12 parameters related to hydrologic and climatic systems, all sampled daily. In addition, a stacked LSTM - a more complex neural network architecture - is used to improve information extraction ability. Exploring river dynamics in the Loire-Bretagne basin and its surroundings, the investigation delves into predictions in daily time steps for one, three, and six months ahead. The resulting forecast features high accuracy and efficiency in predicting river discharge fluctuations, showcasing superior performance in forecasting drought periods over flood peaks. A detailed examination on data used highlights the significance of both local and global datasets in predicting river discharge, where the former dictates short-term predictions, while the latter drives long-range forecasts. Seasonally extended forecasting confirms a strong connection between the forecast leading time and the shift in data correlation, with lower correlation at a lag of 3 months due to seasonal changes affecting forecast quality, compensated by a higher correlation at a longer lag of 6 months. Such mutual effect in this multi-time-step forecasting improves the predictive quality of a six-month horizon, thus encourages progress in long-term prediction to a seasonal scale. The research establishes a practical foundation for effectively utilizing big data to leverage long-term forecasting of environmental dynamics.
