Impact of Snowmelt Conditions on the Isotopic Composition of the Surface Waters of the Upper Ob River during the Flood Period

: For many of the Siberian rivers, and the Upper Ob in particular, 70–80% of the volume of the annual water runoff is formed during the spring ﬂood. Thus, factors inﬂuencing the formation of water runoff during the spring ﬂood are paramount. We explain changes in the isotopic composition of the Upper Ob surface waters by changing different components’ contribution to the runoff water discharge over the spring ﬂood period. We suggest estimating the time of meltwater ﬂow from the Upper Ob watershed to the outlet section using the difference between the date of the complete melting of the snow cover in the catchment area and the date of the maximum light isotope composition of water in the outlet section. We show that a sharp short-term weighting of the isotopic composition of water in the river at the end of the ﬁrst phase of the ﬂood may be associated with the inﬂux of autumn soil moisture, displaced from the soils by snowmelt waters.


Introduction
A significant number of papers on the use of stable water isotopes (deuterium and oxygen-18) in the study of hydrological processes have recently been published [1][2][3][4][5][6][7][8][9][10]. Data on the water isotopic composition allow us to determine the genetic relationship between the surface water and atmospheric precipitation [10][11][12][13], rivers and groundwaters [11,[13][14][15]; to quantify the formation of sea-ice [16]; to evaluate the impact of water evaporation and the rate of water exchange [17,18]; and to study the mechanisms of formation of the chemical composition of water masses coming from tributaries [15,19]. In contrast to the analysis of the isotopic composition of atmospheric precipitation, the study of the isotopic composition of river water has a great advantage, since river water reflects the isotopic composition of the entire watershed rather than a single sampling site [20]. At the same time, the formation of the isotopic composition of surface waters is associated not only with the main hydrological cycle and geomorphological and climatic characteristics of the catchment area, but, according to many researchers, it largely depends on the isotopic composition of precipitation falling on the surface of the catchment area [10,21,22].
Since 2002, the International Atomic Energy Agency (IAEA) at the United Nations (UN) has been conducting research under the project "Isotope tracing of hydrological processes in a large river basin", aimed at studying the isotope composition of water, water balance, and related processes, as well as monitoring environmental changes in the watershed. For this purpose, a global network of data on the isotopic composition of rivers was created-the "Global Network for Isotopes in Rivers". This project is supported by leading research institutes and universities of the world [6,23]. The use of water isotope data is a common tool in hydrological models as it helps in understanding the processes of the water cycle and identifying sources of water masses [8,[24][25][26][27]. Models of hydrological processes using stable isotopes of water as indicators are built mainly for watersheds of temperate and warm climates and are poorly represented for watersheds with long-term seasonal snow cover. Therefore, such studies are especially important for the territory of Siberia, where the rivers are completely covered with ice for a significant part of the year, and their catchment area is covered with a thick layer of snow in winter.
Among the Siberian rivers, the Ob River has the largest catchment area, occupying almost the entire territory of Western Siberia. In this regard, the assessment of the effect of the snowmelt stage on the water balance of the Upper Ob system (south of Western Siberia) during the open-water period provides important information for understanding the features of hydrological processes occurring at the initial stage of the formation of the Ob River runoff at the northern border of the Central Asian watershed.
In this work, we used the stable isotopes of oxygen and hydrogen in precipitation, snow cover, and river water to study the sources of water runoff formation in the territory of the Upper Ob watershed (outlet section in Barnaul) and to assess their contribution to the water balance of the river during the snowmelt and flood period.

Study Area
The studied catchment area of the Upper Ob ( Figure 1a) is typically flat. The water regime in this section of the river is formed by glacial and snow feeding, groundwater, and summer precipitation [28]. The climate in the catchment area is continental, characterized by hot summers and cold winters, and is formed as a result of frequent changes in air masses coming from the Arctic, Central Asia, and the Atlantic. Average annual temperatures are positive at approximately 0.5-2.1 • C. Average maximum temperatures in July are 26-28 • C, and average minimum temperatures in January are (−20)-(−24) • C. The average annual precipitation is 433 mm, about 65% of which falls during the warm period of the year [29]. The snow cover is established mainly in mid-November and collapses in early April; its average depth is 40-60 cm [30]. The relief of the study area is determined by the Ob plateau and the valleys of the Ob River and its tributaries. The left-bank part of the catchment area is located at an altitude of 185-251 m. It is composed of loams in the upper layers and heavy loams in the lower layers; clay outcrops are also noted everywhere [31]. The right bank of the river is floodplain everywhere, up to 10 km wide, with numerous lakes, oxbow lakes, and channels. The altitude of the floodplain is 3-4 m [32,33]. Alluvial formations of the first-third floodplain terraces represented by sands of various compositions and loams are distributed here [34].
The width of the river at the sampling site is about 500 m (Figure 1b), and the depth varies on average from 2-3.5 m during the low-water to 6-8 m during the high-water period. Freezing is usually established by mid-November and lasts for five months. The thickness of the ice layer on the river by the end of the winter period reaches 0.9-1.1 m. In mid-April (or the first ten days of April, at the earliest), the Ob River opens up and the ice drift begins. The periods for the maximum flow and the beginning and the end of the flood period fluctuate significantly from year to year. The hydrograph of the Ob River is characterized by an extended spring-summer flood, in which two phases (waves) can be distinguished, and stable autumn and winter low-water periods. Up to 70-80% of the river's annual runoff passes during the flood period [30,33]. The first phase of the spring flood, associated with the melting of snow from the catchment area of the river, occurs from the end of April. This snowmelt phase includes 15-20% of the total runoff during the entire flood. The waters of this phase largely accumulate in the Ob channel, rarely reaching the floodplain. The maximum levels and floodplain flooding for up to 2-3 months are observed during the second main flood phase, which occurs in May-June and is associated with melting of the snow cover and glaciers in the sources of the Ob River in the mountains [35]. The width of the river at the sampling site is about 500 m (Figure 1b), and the depth varies on average from 2-3.5 m during the low-water to 6-8 m during the high-water period. Freezing is usually established by mid-November and lasts for five months. The thickness of the ice layer on the river by the end of the winter period reaches 0.9-1.1 m. In mid-April (or the first ten days of April, at the earliest), the Ob River opens up and the ice drift begins. The periods for the maximum flow and the beginning and the end of the flood period fluctuate significantly from year to year. The hydrograph of the Ob River is characterized by an extended spring-summer flood, in which two phases (waves) can be distinguished, and stable autumn and winter low-water periods. Up to 70-80% of the river's annual runoff passes during the flood period [30,33]. The first phase of the spring flood, associated with the melting of snow from the catchment area of the river, occurs from the end of April. This snowmelt phase includes 15-20% of the total runoff during the entire flood. The waters of this phase largely accumulate in the Ob channel, rarely reaching the floodplain. The maximum levels and floodplain flooding for up to 2-3 months are observed during the second main flood phase, which occurs in May-June and is associated with melting of the snow cover and glaciers in the sources of the Ob River in the mountains [35].

Sampling
Sampling of the surface water was carried out at the outlet section of the Ob River before the city of Barnaul at three verticals (left bank, middle, right bank) every two weeks from early-March to mid-July (from the end of the winter low-water period to the end of the flood) in 2020-2021. A schematic map of the water sampling sites is shown in Figure

Sampling
Sampling of the surface water was carried out at the outlet section of the Ob River before the city of Barnaul at three verticals (left bank, middle, right bank) every two weeks from early-March to mid-July (from the end of the winter low-water period to the end of the flood) in 2020-2021. A schematic map of the water sampling sites is shown in Figure 1b. Sampling was carried out with a bathometer at a depth of 20 cm from the water surface. Sampling was carried out through drilled holes in the ice during the freeze-up period (for safety reasons, before the ice drift, sampling was carried out only at two verticals near the riverbank), and a boat was used for sampling during the open-water period (after the ice drift). The collected samples were placed in a hermetically sealed plastic container and delivered to the laboratory, where they were filtered through a membrane filter with a 0.45 µm pore diameter using sterile syringes and Minisart NML Plus syringe nozzles. Three to five parallel samples were taken from the resulting filtrate and placed in sealed test tubes, which were stored in a refrigerator until isotopic analysis.
Rain samples were collected during August-October and then during April-June immediately after they fell on the experimental site located on the roof of the building of the Institute for Water and Environmental Problems of the Siberian Branch of the Russian Academy of Sciences (IWEP SB RAS), at a height of 25 m from the ground. Samples were passed through a funnel into a plastic flask-receiver. After the end of rainfall, the flaskreceiver was disconnected from the funnel, and the volume of the sample was measured and filtered through the membrane filter with a 0.45 µm pore diameter (using sterile syringes and Minisart NML Plus syringe nozzles). Three to five parallel samples were taken from the resulting filtrate and placed in sealed test tubes, which were stored in a refrigerator until isotopic analysis.
It is known that the snow cover of temperate and polar latitudes as a whole determines the integral seasonal (for the winter period) characteristics of the isotopic composition of moisture entering the study area [36,37]. For the watersheds of a cold climate, the contribution of winter precipitation to the water runoff occurs only from the beginning of the snowmelt period; therefore, the bulk snowpack sampled just before the snowmelt is a convenient object for studying the distribution of winter snowfalls over the catchment area. For this purpose, the bulk snowpack samples were taken at the sites of the network located in the catchment area above the studied outlet section (Figure 1a) at the time of maximum snow accumulation (during 10-14 March in 2020 and 15-18 March in 2021). All sites were located in a field on flat territory free from trees and bushes. Sampling was performed using the envelope method (10 m × 10 m). The composite sample consisted of 5 snow pits collected with a plastic pipe of 4.5 cm inner diameter. After collection, composite samples of the snowpack were placed in clean tight-closing plastic bags and stored frozen until analysis. Before instrumental analysis, snow samples were transferred into specially prepared closed plastic containers [36] and melted at room temperature. Then, by analogy with river and rain waters, they were filtered through a membrane filter with 0.45 µm pore diameter; three to five parallel samples were taken from the resulting filtrate and placed in sealed test tubes.
Over the 2 years, 147 samples (55 surface water, 24 snowpack, 24 rainfalls in the spring-summer of the current year, and 44 rainfalls in the autumn of the previous year) were taken and analyzed.

Analytical Methods
The main experimental and theoretical positions of the isotopic systematics of deuterium and oxygen-18 were developed in the works of Craig and Dansgaard [38,39]: where R sample and R std are the ratios of 2 H/ 1 H and 18 O/ 16 O in the measured sample and in the standard, respectively. The relationship between the ratio of stable isotopes of oxygen (δ 18 O) and hydrogen (δD or δ 2 H) in precipitation is described by an empirical relationship called the global meteoric water line (GMWL) [38,40]: The ratio of δ 18 O and δD in precipitation at a given location can be described by its local meteoric water line (LMWL), whose deviation from GMWL can be used to estimate the isotope fractionation of local precipitation. In isotope studies, the calculation criterion d-excess proposed by Dansgaard [39] is used for this purpose, together with measurements of δ 18 O and δD. The d-excess parameter is associated with the kinetic processes of isotope fractionation, which characterize the processes of evaporation or freezing. Seasonal variations in d-excess can provide information about moisture sources and their hydrological cycle, as well as the influence of local climatic conditions [41]. The deuterium excess (d-excess) is calculated as: The stable isotope (δ 18 O, δD) analysis of samples was carried out at the Chemical Analytical Center of IWEP SB RAS. The analysis was performed by laser absorption IR spectrometry on a PICARRO L2130-i instrument (WS-CRDS). The measurement accuracy of δD and δ 18 O (1σ, n = 5) was ±0.4‰ and ±0.1‰, respectively. International standards GRESP and USGS-47 were used for instrumental calibration.
Precipitation-weighted mean values of stable isotope compositions (δD, δ 18 O, and d-excess) in precipitation for a certain period were calculated by the formula: where Cvwm is the precipitation-weighted mean values for the calculated period, ‰; Cj is the value of δD, δ 18 O, or d-excess in the j-th precipitation sample, ‰; Qj is the amount of the j-th precipitation sample, mm weq.; and Q is the total precipitation amount for the calculated period, mm weq.

Results
The study of the isotopic composition of surface water in the Ob River near Barnaul City showed that over the period from March 15 to July 15, its isotopic composition varied: To analyze the temporal variability of the isotope composition of water in the Ob River during the spring flood, curves of the change in water discharge and water isotopic composition (δ 18 O) were plotted for three verticals of the river outlet section from 15 March to 15 July 2020 and 2021 ( Figure 2). The water-discharge curve was constructed according to data from the daily discharges measured at the hydro-meteorological station located 300 m downstream from the sampling site [42]. This curve makes it possible to identify the phases of the water regime in accordance with the change in the values of the water discharge in the river section. For the Upper Ob, it often happens that the first and the second flood phases overlap each other and are not identified on the discharge curve [30,33]. In our case, especially for 2021, the flood phases can be distinguished.
In 2020 and 2021, the lightest isotopic compositions of river water were at the end of the first and second flood phases (Figure 2), while in each year, at the beginning of the second flood stage, a short-term weighting of the isotopic composition of water (by 0.5-0.9‰ in terms of δ 18 O) was observed. Temporal changes in the isotopic composition of water at different verticals of the river for both years also have a similar character (with small local differences). Table 1 presents the average values of the isotopic composition of water in the Ob outlet section near Barnaul in the different seasons of the hydrological regime of the river in 2020 and 2021: the winter low-water (under-ice period), the first phase of the spring flood (snowmelt period, opening of the river from ice), and the second phase of the spring flood (melting of snow and glaciers in the mountains at the source of the Ob River). For comparison, the table contains the precipitation-weighted mean values of the isotope composition of rainfall during the autumn period of the previous year (from August to October) and during the current flood period, as well as the depth-weighted mean values of the snow cover of the past winter. In 2020 and 2021, the lightest isotopic compositions of river water were at the end of the first and second flood phases (Figure 2), while in each year, at the beginning of the second flood stage, a short-term weighting of the isotopic composition of water (by 0.5-0.9‰ in terms of δ 18 O) was observed. Temporal changes in the isotopic composition of water at different verticals of the river for both years also have a similar character (with small local differences).  Figure 1c) and water discharge in the Ob River during the flood period in 2020 and 2021. I-the end of the winter low-water period (under-ice period); II-the first phase of the spring flood; III-the second phase of the spring flood; IV-the beginning of the summer low-water period; the vertical arrow line indicates the date of the ice drift. Table 1. Average values of the isotopic composition of the surface water of the Ob outlet section in different periods of the hydrological regime in 2020 (1) and 2021 (2); precipitation-weighted mean of autumn (August-October 2019 (1) and 2020 (2)) and spring-summer (first and second phases of the floods in 2020 (1) and 2021 (2)) rainfalls; depth-weighted mean of snow cover in winter 2019-2020 (1) and 2020-2021 (2). According to the results of the analysis of the snow cover, the average value of the isotopic composition of atmospheric precipitation was −20.3‰ for δ 18 O and −157.8‰ for δD in winter 2019-2020, and it was −19.8‰ for δ 18 O and −151.3‰ for δD in winter 2020-2021 (Table 1).
The data in Table 1 indicate similar average values of the isotopic composition of surface water in the Ob outlet section in the different hydrological seasons of 2020 and 2021. For all studied hydrological seasons, the difference in values between 2020 and 2021 was 0.2-0.3‰ for δ 18 O and 0.1-0.8‰ for δD. At the same time, the heaviest isotopic composition of water was observed in the winter low-water period, and in the first phase of the flood, it was lighter relative to the winter low-water period by 0.6-0.7‰ and 4.4-5.1‰, and in the second phase by 0.6‰ and 5.2-5.4‰ for δ 18 O and δD, respectively. The value of the deuterium excess (d-excess) in all studied hydrological periods indicates that the river water was affected by evaporative fractionation. Among all objects of study, snowmelt waters had the lightest isotopic composition, and the difference between the average values between years was quite small (0.5‰ for δ 18 O and 6.5‰ for δD). The spring-summer rainfall had the heaviest isotopic composition, and the autumn precipitation significantly affected by evaporative fractionation was the most depleted in deuterium atoms (δD) relative to δ 18 O.

Discussion
The local meteoric water lines (LMWLs) for the Ob River water during the spring-summer flood period in 2020 and 2021, as well as for the snowpack melting water in the catchment area in 2020 and 2021 and autumn rainfalls in previous 2019 and 2020, are shown in Figure 3. For comparison, the figure also shows the global meteoric water line (GMWL). It follows from Figure 3 that for both years, the values of the slope coefficients in the LMWL equations for rainfalls and river water are somewhat lower than the slope of the global meteor water line (GMWL) described by the equation δD = 8 × δ 18 O + 10 [38], while the slope of LMWL for the snow cover did not differ from the GMWL one. The slopes of LMWL equations for rainfalls and river water indicate the process of evaporative fractionation of their initial isotopic composition of water or moisture, and the values of d-excess (Table 1) indicate maximum depletion in deuterium atoms relative to δ 18 O for autumn precipitation.
in the LMWL equations for rainfalls and river water are somewhat lower than the slope of the global meteor water line (GMWL) described by the equation δD = 8×δ 18 O + 10 [38], while the slope of LMWL for the snow cover did not differ from the GMWL one. The slopes of LMWL equations for rainfalls and river water indicate the process of evaporative fractionation of their initial isotopic composition of water or moisture, and the values of d-excess (Table 1) indicate maximum depletion in deuterium atoms relative to δ 18 O for autumn precipitation.    Using the data from Table 1, and taking into account that during the flood period, the contribution of groundwater runoff and spring precipitation to the water flow in the river can be neglected, the contribution of the autumn rain and snow precipitation to the water flow in the river at the outlet section was estimated for the snowmelting period (the first phase of the flood). For a preliminary assessment of the autumn rain and snow contributions, the following equations were solved: where X is the share of snowmelt water in the river runoff, (1 − X) is the share of autumn precipitation in the river runoff; and the numbers are the average values of the isotopic composition of the surface water, autumn precipitation, and snow cover (Table 1). Solving these equations showed that during the snowy winter of 2019-2020 (where snow cover reached up to 170 mm weq.) and lower autumn precipitation in 2019, the contribution of melted snow to the water flow in the outlet section during the first phase of the flood in 2020 reached 52%, while during winter 2020-2021 (snow cover reached only 115 mm weq.) and a larger amount of autumn rainfall in 2020, it dropped to 43% in 2021. Figure 2 shows that the isotopic composition of the water in the Ob River begins to become lighter due to snowmelt more than a week before the ice drift, and the sharpest change in the isotopic composition was observed in the spring of 2021, although a rise in the water level and an increase in water discharge were noted earlier in 2020 [37,43,44]. In our opinion, these differences can be associated with the conditions of snow accumulation, the main characteristics of which are given in Table 2 for the winters of 2019-2020 and 2020-2021. Table 2. Conditions of snow accumulation and snowmelting in winter-spring 2019-2020 and 2020-2021 in the study area near the city of Barnaul according to [45].

2019-2020 2020-2021 Standard Average [34]
The day of snow cover formation 13 Reference information on meteorological parameters (daily data on air temperature, precipitation, snow depth) of the study area (Table 2) was taken from the site "Specialized data sets for climate studies" [45], where, among other data, data from the weather station located in the suburbs of Barnaul were published. Based on these daily data, we calculated the average monthly daily air temperatures during the winters 2019-2020 and 2020-2021. In addition, we determined the day of the beginning of the formation and the end of the melting of the snow cover (according to daily data on snow depth from 0 cm in early November to 0 cm by early April) and the time of the onset of the climatic spring (the beginning of stable positive average daily temperatures), as well as the depth of the snow cover at the first low winter temperatures (<−25 • C) lasting more than three days. The data presented in the "standard average" column in Table 2 were taken from [34].
The authors [46], when studying the features of floods on the tributaries of the Upper Ob, concluded that if the height of the snow cover in the watershed before the onset of significant frosts is more than 20-25 cm, then the soil under the snow does not freeze and melt water intensively penetrates and is retained in the soil in the spring. According to Table 2, stable snow cover in the study catchment area was established on 13 November for both years; however, the depth of snow cover before the first severe frosts lasting more than three days was 20 cm at the beginning of winter 2019-2020 and only 10 cm at the beginning of winter 2020-2021. Thus, in the winter of 2020-2021, deep freezing of the soil could have occurred, and melted snow water of light isotope composition (Table 1) could have directly flowed to the riverbed along the frozen soil surface during the snowmelt from 5 to 13 April 2021 (Table 2), which could explain the sharp change in the isotopic composition of river water from the very beginning of the snowmelt (Figure 3).
The climatic spring in 2020 began a week earlier than in 2021, but since the soil under the snow was not frozen, the first portions of meltwater were absorbed and retained by the upper soil layer, and a sharp change in the isotopic composition of a river water began only by the end of the snowmelt (by 09.04.2020, Table 2 and Figure 2).
The lowest values of the isotope composition of river water in the first phase of the spring flood indicate the maximum contribution of snowmelt water to the water discharge in the riverbed. Therefore, using the difference between the date of the complete melting of the snow cover in the catchment area and the date of its maximum contribution to the river discharge (the lightest isotope composition of water in the outlet section), it is possible to estimate the time it takes for meltwater to flow from the surface of the Upper Ob catchment area (the section of the river from Biysk to Barnaul) to the outlet section. In 2020, the end of the snowmelt fell on 9 April, and in 2021 on 13 April (Table 2), while the maximum contribution of meltwater was on 20 April in 2020, and on 23 April in 2021 ( Figure 2). From here, the time for the melted snow runoff to reach the outlet section in 2020 was 11 days, and in 2021, 10 days. It should be noted that during the first phase of the spring flood, as a rule, a small amount of precipitation falls in the region under study (in 2020 and 2021, the amount of rainfall was 5 and 4 mm, respectively), which, in our opinion, does not warrant significant adjustments in the calculation of the time for melted snow water to reach the outlet section.
Both in 2020 and 2021, during the first phase of the flood, the weighting of the isotopic composition of water in the river outlet section begins with a decrease in the peak of water discharge (Figure 2). We associate the initial stage of the weighting of the isotope composition of water with the possible discharge of groundwater due to the expulsion of autumn moisture by snowmelt waters. During this period, the groundwater flow is a mixture of last year's autumn rainfall and melted snow water of the current year.
In the second phase of the spring flood, changes in the isotopic composition of river water are determined by the mixing of relatively isotopically heavy spring-summer rainfalls and maximally isotopically lighter melt snow and glacial waters in the high mountains at the source of the Ob River. At this time, due to the maximum high-water level in the river (Figure 2), the contribution of groundwater to the Ob runoff is significantly reduced, so the inflow of rainwater is mainly associated with surface flat runoff. As in the case of snowmelt, by the end of the second phase of the flood, the lowest values of the isotopic composition of river water are noted. This situation characterizes the maximum contribution of melted snow and glacier waters from the high mountains of the Ob source to the water flow in the river channel at the outlet section.

Conclusions
This paper determined the features of the formation of water runoff in the outlet section of the Ob River near Barnaul during the spring floods in 2020-2021 using the stable isotope composition of these runoff components, namely, the snow cover of the catchment area, as well as river water and rain precipitation in different phases of the hydrological regime of the river. For the Siberian rivers, and for the Upper Ob in particular, the study of the factors influencing the formation of water runoff during the spring flood is very important, because 70-80% their annual water runoff occurs during the spring-summer flood. In the present study, we paid special attention to the first phase of the flood (the snowmelt period). The results can be summarized as follows:

•
The lightest isotopic composition of the Upper Ob River water occurred at the end of the first and the second flood phases due to the influx of melted snow from the flat catchment area (first phase), and then due to melted snow and glacial waters from the high mountains (second phase).
• In our opinion, a short-term weighting of the water isotope composition by 0.5-0.9‰ in terms of δ 18 O at the beginning of the second stage of the flood is attributed to the discharge of groundwater due to the expulsion of autumn moisture by snowmelt waters. During this period, the groundwater flow is a mixture of last year's autumn rainfalls and the melted snow water of the current year.

•
The main factors controlling the rate and volume of snowmelt water flow from the catchment area to the Upper Ob channel are the initial conditions of the snow cover formation (autumn moisture content and freezing of the soil in the catchment area at the beginning of winter). Therefore, monitoring these conditions will allow the prediction of the magnitude and features of the first phase of the flood. • An assessment of the rain and snow contributions to the Upper Ob River runoff during the snowmelt period was carried out using the average values of the isotopic composition of the surface water, the precipitation-weighted mean of autumn rainfalls, and the depth-weighted mean of snow cover. The calculations showed that the contributions of the melted winter snow and autumn rainfalls from the previous year to the river water flow were, respectively, 52% and 48% in spring 2020 and 43% and 57% in spring 2021.

•
The difference between the date of complete melting of the snow cover in the catchment area and the date of the maximum light isotopic composition of water in the outlet section was proposed for estimating the time for meltwater runoff from the catchment area to the outlet section. This substantiated approach allowed us to estimate the time for the melted snow runoff to reach the outlet section of the Upper Ob River, being 11 and 10 days in 2020 and 2021, respectively.