High mountain areas are prone to extreme hydrological events, and their study is especially important in the context of ongoing intensive deglaciation. In this research, a model “chain” consisting of a hydrodynamic model and a runoff formation model is adopted to simulate a glacier lake outburst flood (GLOF) from Bashkara Lake (the Central Caucasus, Russia) and its effect on downstream. In addition to an actual GLOF event that occurred on 1 September 2017 and led to casualties and significant destruction in the Adylsu and Baksan Rivers valleys, possible scenarios for the re-outburst of the lake are considered. The hydrographs of the outburst and the downstream movement of the flood wave along the Adylsu River valley are estimated using STREAM_2D two-dimensional hydrodynamic model. The water discharges in the entire river network of the Baksan River are assessed using the ECOMAG (ECOlogical Model for Applied Geophysics) runoff formation model. The output flood hydrograph from the hydrodynamic model is set as additional input into the Baksan River runoff formation model in the upper reaches of the Adylsu River. As a result of the simulations, estimates for the contribution of GLOFs and precipitation to an increase in peak discharge along the Baksan River were obtained. The actual outburst flood contributed 45% and precipitation 30% to the peak flow in the Baksan River at the mouth of the Adylsu River (10 km from the outburst site). In Tyrnyauz (40 km from the outburst site), the contributions of the outburst flood and precipitation were equal and, in Zayukovo (70 km from the outburst site), the outburst flood contributed only 20% to the peak flow, whereas precipitation contributed 44%. Similar calculations were made for future potential re-outburst flood, taking into account climatic changes with an increase in air temperatures of 2 °C, an increase in precipitation of 10% in winter and a decrease of 10% in summer. The maximum discharge of the re-outburst flood in the Adylsu River mouth, according to model estimations, will be approximately three times less than the discharge of the actual outburst on 1 September 2017 and can contribute up to 18% of the peak discharge in the Baksan River at the confluence.

This study presents a dataset on long-term multidisciplinary glaciological, hydrological, and meteorological observations and isotope sampling in a sparsely monitored alpine zone of the North Caucasus in the Djankuat research basin. The Djankuat glacier, which is the largest in the basin, was chosen as representative of the central North Caucasus during the International Hydrological Decade and is one of 30 “reference” glaciers in the world that have annual mass balance series longer than 50 years (Zemp et al., 2009). The dataset features a comprehensive set of observations from 2007 to 2017 and contains yearly measurements of snow depth and density; measurements of dynamics of snow and ice melting; measurements of water runoff, conductivity, turbidity, temperature, δ¹⁸O, δD at the main gauging station (844 samples in total) with an hourly or sub-daily time step depending on the parameter; data on δ¹⁸O and δ²H sampling of liquid precipitation, snow, ice, firn, and groundwater in different parts of the watershed taken regularly during melting season (485 samples in total); measurements of precipitation amount, air temperature, relative humidity, shortwave incoming and reflected radiation, longwave downward and upward radiation, atmospheric pressure, and wind speed and direction – measured at several automatic weather stations within the basin with 15 min to 1 h time steps; gradient meteorological measurements to estimate turbulent fluxes of heat and moisture, measuring three components of wind speed at a frequency of 10 Hz to estimate the impulse of turbulent fluxes of sensible and latent heat over the glacier surface by the eddy covariance method. Data were collected during the ablation period (June–September). The observations were halted in winter. The dataset is available from PANGAEA (https://doi.org/10.1594/PANGAEA.894807, Rets et al., 2018a) and will be further updated. The dataset can be useful for developing and verifying hydrological, glaciological, and meteorological models for alpine areas, to study the impact of climate change on hydrology of mountain regions using isotopic and hydrochemical approaches in hydrology. As the dataset includes the measurements of hydrometeorological and glaciological variables during the catastrophic proglacial lake outburst in the neighboring Bashkara valley in September 2017, it is a valuable contribution to study lake outbursts.

STREAM_2D software package was applied to retrospective and predictive simulations of the Lena River near city Yakutsk hydraulic and channel changes during ice-free period. The modelling results indicate significant correspondence of simulated water discharges distribution and water levels with observed ones for the period 2001–2016. Model has captured main erosion and depositional zones observed in the 2009–2016 years. Combination of typical monthly average hydrographs of 1 %, 10 % and 50 % return intervals was used as an initial parameter for 10-year channel changes forecast. According to the simulation results, degradation of the Adamovskaya branch will be continued, which is the most negative impact for the Yakutsk city infrastructure maintenance. At the end of the forecast period the equal distribution between Adamovskaya and Buorylarskaya channels is possible, which is a positive trend from the point of view of the city water supply system.

The water resources of the North Caucasus depend mostly on the state of glaciers, which have been intensely losing their mass in the recent decades against the background of climate changes. The deglaciation leads not only to a decrease in the glacier runoff of mountain rivers, but also to changes in the annual distribution of runoff. The focus of this study is the adaptation of ECOMAG software complex to simulating river runoff in the Baksan River basin based on data on the relief and underlying surface of the drainage basin (soil, vegetation) and daily data on the surface air temperature, air saturation deficit, and precipitation. The calibration and validation of the model and the statistical estimate of calculation efficiency were based on the data on water discharges in the Baksan River over 2000–2017. The developed model of runoff formation in the Baksan River basin was used to carry out numerical experiments for assessing the sensitivity of runoff characteristics to glacier area variations. Depending on the rate of deglaciation process, the runoff of the Baksan River can drop by 10–30% because of a decrease in its glacial component, and the maximal water discharges can drop by 10–15%.

In this study, we adapted the ECOMAG model of the runoff formation for analysis of the Terek River basin using comprehensive hydrometeorological information as well as data on soils, landscape, and glaciation. To take account of regional characteristics of the glaciation, the additional ice module was used with the model. This improvement has resulted in a satisfactory agreement between the modeled runoff hydrographs and the observed ones. In our simulations we used the updated glacier cover predictions from the- global glaciological model GloGEMflowdebris together with regional climate projections from the CORDEX experiment to determine possible future changes in the Terek River flow in the 21st century. The results show that the runoff will change between −2% and +5% according to the RCP2.6 scenario, and from −8% to +14% in the RCP8.5 scenario. The directedness of the runoff changes in particular subbasins of the River will essentially depend on the altitude position of the snow and glacier feeding zones, that is responsible for the intensity of their degradation. Thus, in the RCP8.5 scenario, the flow of the Chegem River will begin to decrease significantly in the second half of the 21st century. In contrast, the predicted increasing of the runoff in Malka and Baksan rivers, which are primarily fed by meltwater from glaciers and snow on Elbrus and other high-mountain zones, is expected to be continued until the end of the century. But this increase may be caused only by a growth of a part of the snowmelt feeding due to greater winter precipitation. The model estimates confirm the present-day observed trends within the intra-annual runoff distribution, demonstrating the earlier start of the spring flood, a decrease in summer runoff volumes and then its increase in the autumn months. The results of the research may be used for more efficient management of water resources in the North Caucasus in the future, including electricity generation and water supply.