In spring, rivers at middle and high latitudes in the Northern Hemisphere are prone to ice jams, which threaten the safety of hydraulic structures in rivers. Heilongjiang Province is located on the highest latitude in China, starting at 43°26′ N and reaching 53°33′ N. Rivers in Heilongjiang Province freeze in winter and break up in spring. Forecasting the break-up date of river ice accurately can provide an important reference for the command, dispatch, and decision-making of ice flood preventing and shipping. Based on the observed break-up date series of river ice from seven representative hydrological stations in Heilongjiang Province, the complete ensemble empirical mode decomposition with adaptive noise (CEEMDAN) was used to decompose the observed break-up date series of river ice into several subsequences, and the long-short term memory neural network (LSTM) was used to forecast the subsequences decomposed by CEEDMAN. Then, the forecast results of each subsequence were summed to obtain the forecasting value for the break-up date of river ice proceeded by CEEMDAN-LSTM. Compared with the LSTM, the forecast accuracy of CEEMDAN-LSTM for the break-up date of river ice had been significantly improved, with the mean absolute error reduced from 0.80–6.40 to 0.75–3.40, the qualification rate increased from 60–100% to 80–100%, the root-mean-square difference reduced from 1.37–5.97 to 0.95–1.69, the correlation coefficient increased from 0.51–0.97 to 0.97–0.98, and the Taylor skill score increased from 0.87–0.99 to 0.99. CEEMDAN-LSTM performed well in forecasting the break-up date of river ice in the Heilongjiang Province, which can provide important information for command, dispatch, and decision-making of ice flood control.

This study aimed to determine the influence of time and precipitation as covariates on the flood frequency distribution in the Songhua River tributaries under the nonstationarity assumption and to investigate the possibility of nonstationary models’ application in river management scope demarcation work. Nonstationary flood frequency analysis (NS-FFA) was conducted in three typical basins of the Songhua River (in Northeastern China) based on the generalized additive models for location, scale, and shape (GAMLSS), and stationary flood frequency analysis was used as a comparison. Under the stationarity assumption, the Pearson type Ⅲ (P-Ⅲ) distribution is the main theoretical distribution for the flood extremum at hydrological stations, followed by a lognormal (LN) distribution. Under the nonstationarity assumption, when time is considered a covariate, the optimal theoretical distribution of the flood extremum is mainly LN (with 63.75%), followed by the Weibull distribution (with 18.75%). When precipitation is considered as a covariate, the optimal theoretical distribution of the flood extremum is mainly LN (with 57.5%). We attempted to apply several FFA methods to calculate the design frequency in this study, referring to the work requirements for river management scope demarcation in three typical basins, and came to the following conclusions. From the simulation results of the p = 10% flood at the export stations of typical basins, it can be seen that time-covariate NS-FFA obtained the best simulation results. Two cases of the simulation under the stationarity assumption are positive, which will lead to a high design scale. The time-covariate GAMLSS in NS-FFA has the advantages of higher calculation accuracy and simpler processes. To better balance construction costs and disaster protection requirements, NS-FFA can be used to determine the design scale of water conservation projects; additionally, it can be used to demarcate the scope of river management. The accuracy of GAMLSS for FFA is also influenced by the complexity of the terrain, with basins with relatively simple terrain having higher calculation accuracy.