This study investigates the spatiotemporal variability of minimum river discharge in Northern Kazakhstan and its sensitivity to climatic factors under continental and semi-arid conditions. Long-term meteorological and hydrological data from key regional rivers, including Tobyl, Toguzak, Torgay, Karatorgay, and Ayat, were analyzed to detect structural shifts, trends, and stationarity. Structural changes were identified using Pettitt’s test and supported by difference-integral and cumulative-integral curves, while linear trend analysis with 10-year moving averages highlighted long-term tendencies. Correlation analysis assessed the influence of temperature and precipitation, distinguishing low- and high-flow years. Stationarity and homogeneity were evaluated with ADF, KPSS, and complementary parametric and non-parametric tests to ensure robust statistical interpretation. Results reveal significant structural shifts during 1980–1990, with rising winter air temperatures and marked heterogeneity of minimum discharge, especially in cold seasons. Winter flows are predominantly controlled by temperature and groundwater recharge, while summer flows are strongly linked to precipitation dynamics and evapotranspiration processes. Differences between low- and high-flow years confirm the dynamic non-stationarity of the system. These findings enhance understanding of regional hydro-climatic interactions and provide a scientific basis for adaptive water resource management, reservoir regulation, and long-term strategies of climate change adaptation in the steppe and semi-arid regions of Central Asia.
INTRODUCTION
Within a framework for global climate change, there is an observed increase in air temperatures combined with shifts in precipitation patterns, which together impact hydrological cycles across different spatial and temporal scales (Trenberth et al., 2014; IPCC, 2021). River catchments in arid and semi-arid regions are particularly sensitive to climate variations, making low flow a key indicator of both water supply sustainability and ecological health in catchments within a river basin (Kundzewicz et al., 2008; Shiklomanov & Lammers, 2020; Birimbayeva et al., 2024; Makhmudova et al., 2024; Abdykadyrov et al., 2025).
Analysis of low-flow discharge allows for a determination of long-term trends in either deterioration or recovery of the water regime, as well as for an estimation of climatic sensitivity of hydrological catchments. In north Kazakhstan and in adjacent Soviet Urals regions marked by a prevalence of steppe landscapes with a sharply continental climate, minimum discharge values, primarily those for summer and winter periods with low flow, are dependent upon an interaction between thermobaric regimes, freezing and moisture-accumulation regimes, and infiltration- and groundwater contributions (Dyachenko & Maznev, 2017; Belyaev et al., 2019; Makhmudova et al., 2024; Smagulov et al., 2025).
Despite a substantial amount of empirical data, the spatial variation in how climate impacts low flow remains poorly understood. This lack of knowledge applies both to differences between winter and summer regimes and to how much these relationships are influenced by year-to-year changes in water availability. A major research effort should focus on examining structural breaks in climatic and hydrological time series and on understanding the relative roles of precipitation and temperature in driving extremely low water discharge values (Zhang et al., 2020; Ivanov & Shlykov, 2022).
This study aims to provide a comprehensive assessment of climatic factors, namely precipitation and air temperature, and their effect on low-flow discharge at hydrological measurement stations in northern Kazakhstan. The analysis includes consideration of seasonal features, between-year variations in wet and dry periods, alongside structural shifts and longer-term trends.
Research includes a comprehensive overview of the geographical location of measurement stations and climate zones, sources and duration of observations, along with employed analytical procedures, including statistical analysis, trend analysis, and correlation analysis. Special focus is given to breakpoint selection criteria in addition to year classification for water availability.
MATERIALS AND METHODS
The research scope comprises hydrological gauging stations located in small to medium catchments in north Kazakhstan, specifically the Toguzak River, the Tobyl River (in proximity to Grishenka village and Kostanay city), the Ayat River (reach near Varvarinka village), the Karatorgay River, and the Torgay River. These gauging stations encompass lowland as well as foot-slope catchments within the belted basin of Esil and some parts of the basin of Tobyl (Figure 1). Low flow discharge was estimated independently for winter and summer low flow periods, thereby allowing consideration of variations in feeding mechanisms like groundwater, infiltration, and precipitation.
|
|
|
Figure 1. Tobyl–Torgay water management basin (Nugmanov et al., 2025) |
Data for discharge cover the period from the 1930s up to 2020 and were taken from RSE "Kazhydromet" archival collections and processing reports. Also, auxiliary information about air temperature and precipitation was taken from the Kostanay meteorological station, providing a fine approximation for the climatic conditions of the study area. The observation period includes 1935–2021.
The analysis employed a range of statistical and graphical forms of analysis commonly used in hydrology and climatology. To identify structural breakpoints in the time series for temperature, precipitation, and low-flow discharge, use was made of the Pettitt test (Pettitt, 1979), so as to make it possible for a shift in a dataset's median structure with an indeterminate transition duration to be identified at a specific point. Analysis of trend was performed using linear regression, calculation of the slope for the time series, and assessment of significance using the Student's t-test, following procedures outlined in previous literature (Kundzewicz & Robson, 2000). Tests for stationarity were performed using the Augmented Dickey-Fuller (ADF) test and the KPSS test, which allow for an estimation of the presence or absence of a unit root or a component of trend (Said & Dickey, 1984, Kwiatkowski et al., 1992). To test the correlation between climatic variables and discharge at low flow, the Pearson pairwise correlation coefficient was used. Multivariate analysis was conducted on the full set of data as a whole and separately for wet years (Qmax ≥ 0.25 percentile) and dry years (Qmin ≤ 0.75 percentile). Such an approach allowed for an estimation of modified climatic impacts on river discharge for varying hydrological conditions.
Further, the residual mass curve as well as the cumulative sums technique was employed, a standard technique for representing extended variations within hydrological settings (Burn & Hag Elnur, 2002; Maznev, 2009).
RESULTS AND DISCUSSION
Breakpoints and trends in climatic parameters
The analysis of long-term series of air temperature and precipitation from the Kostanay meteorological station for the period 1935–2021 revealed the presence of structural breakpoints and persistent positive trends, particularly pronounced in temperature variables. Application of the non-parametric Pettitt test (Pettitt, 1979), designed to detect unknown points of structural change in time series, indicated statistically significant breakpoints at the following moments:
In all cases, the detected trends are positive and statistically significant (p < 0.01), with the highest rate of increase observed for cold-season temperature (+0.0491 °C per year, R² = 0.61). The analysis of long-term climatic variables—air temperature and precipitation—thus demonstrates persistent trend shifts and structural non-stationarity, particularly evident in temperature characteristics.
It should be noted that in all of the above cases, temperature series display stable positive dynamics, confirmed by linear trend analysis. The most pronounced increase was recorded during the cold season, where the rate of warming reached +0.0491 °C per year, with a coefficient of determination (R²) of 0.61, indicating a high level of explained variance and robustness of the trend.
By contrast, the analysis of precipitation revealed substantially lower variability. Weakly positive trends were observed, but they were not always statistically significant. Precipitation exhibited lower variability and a weakly positive trend, with significant breakpoints detected only for the cold season; the structural shift in precipitation sums was identified at the boundary of significance (p = 0.063).
Thus, it can be concluded that temperature regimes exhibit a more pronounced and persistent dynamic than precipitation regimes, which is consistent with global and regional climatic trends documented in comparable studies (Kundzewicz et al., 2008, IPCC, 2021).
Hydrological restructuring and changes in low-flow discharge
The low-flow regime of the rivers analyzed in this investigation exhibits distinct indications of structural alterations in hydrological dynamics, evident in both the statistical framework of time series and the trajectory and magnitude of long-term trends. The analysis indicated that all monitoring stations within the region underwent structural shifts during the interval from 1976 to 1992, as determined through the application of the Pettitt test (Pettitt, 1979) (Table 1). In various instances, these shifts were associated with modifications in trend slope—from stable or positive values to negative ones—most notably observed for the Toguzak River and the Tobyl River (Grishenka station) during the winter low-flow season.
Table 1. Structural shifts of minimum flow in rivers of the Tobyl–Torgay water management basin
|
No. |
Observation Point |
Year of Structural Shift |
Pettitt Statistic |
p-value |
Note |
|
Summer period |
|||||
|
1 |
Toguzak – s. Togyzak |
1989 |
615.0 |
0.09203 |
– |
|
2 |
Tobyl – s. Grishenka |
1985 |
1044.0 |
0.00028 |
Climate humidification, flow stabilization |
|
3 |
Tobyl – g. Kostanay |
1989 |
1129.0 |
0.00006 |
Increase in minimum flow after regulation (Karatomar Reservoir) |
|
4 |
Kara-Torgay – s. Urpek |
1981 |
1077.0 |
0.000005 |
Sharp zeroing of flow after the 1990s |
|
Winter period |
|||||
|
1 |
Toguzak – s. Togyzak |
1992 |
1364.0 |
0.00000006 |
Transition from zero values to stable flow after 1992 |
|
2 |
Tobyl – s. Grishenka |
1976 |
1021.0 |
0.000085 |
No flow before 1976, followed by a stable increase |
|
3 |
Tobyl – g. Kostanay |
1988 |
1335.0 |
0.000001 |
Sharp rise in winter flow after 1988 |
|
4 |
Kara-Torgay – s. Urpek |
– |
– |
– |
The river freezes to the bottom in winter, with the flow absent. |
The structural breakpoints identified by the Pettitt test were further confirmed by residual mass and cumulative sum curves (Figure 2).
The patterns that were observed underscore the seasonal specificity of river regime responses to alterations in climatic variables. In the winter low-flow period, the minimum discharges exhibit the greatest sensitivity to cold-season temperatures, indicating the impact of soil freezing depth, snow cover conditions, and infiltration feeding. Additionally, a consistent positive correlation was identified with mean annual air temperature, which serves as an indicator of an integrated thermal effect.
For low flows in summer, autumn-winter precipitation is typically the dominant driver. However, where evaporation loss is strong in arid climatic regimes, there may be a reverse relationship with temperature: high summer temperature leads not to an increase but a decrease in discharge due to greater evaporative loss in addition to a decrease in recharge to moisture in soils.
Therefore, the results confirm the occurrence of several climate-hydrological processes in controlling low-flow discharge in both winter and summer seasons. Such delineations are thus essential for developing adaptive measures for water resource management in a situation involving an increase in climate variability.
|
|
|
|
a) |
b) |
|
|
|
|
c) |
d) |
|
|
|
|
e) |
f) |
|
|
|
|
g) |
h) |
|
Figure 2. Determination of the breakpoint in stationarity for the case of increase, stability, and decrease of minimum summer discharges (a, b – Tobyl River, Grishenka station; c, d – Tobyl River, Kostanay station; e, f – Toguzak River, Toguzak station; g, h – Karatorgay River, Urpek station). Left: breakpoint identified using the Pettitt test; right: by cumulative ΣQ and residual mass Σ(Ki-1)/Cv curves. |
|
Impact of temperature and precipitation on low-flow discharge
A comparison of long-term low-flow series with climatic parameters (temperature and precipitation) revealed a strong dependence of discharge on climate factors, with pronounced seasonal and spatial variability (Table 2).
Table 2. Correlation analysis between climatic parameters and minimum river flow
|
№ |
Observation point |
Correlation coefficient between temperature and minimum flow |
Correlation coefficient between precipitation and minimum flow |
||||
|
Year |
Cold period |
Warm period |
Year |
Cold period |
Warm period |
||
|
1 |
Toguzak – s. Togyzak |
0,37/-0,1 |
0,41/0,06 |
0,16/0,32 |
0,33/0,43 |
0,27/0,1 |
0,21/0,43 |
|
2 |
Tobyl – s.Grishenka |
0,28/0,68 |
0,31/0,60 |
0,13/0,55 |
0,34/0,28 |
0,27/0,34 |
0,23/0,12 |
|
3 |
Tobyl – Kostanay |
0,31/0,68 |
0,31/0,60 |
0,19/0,57 |
0,11/0,27 |
0,20/0,34 |
0,014/0,10 |
|
4 |
Ayat – s. Varvarinka |
-/0,10 |
-/0,23 |
-/-0,18 |
-/0,34 |
-/0,12 |
//0,31 |
|
5 |
Kara-Torgay – s. Urpek |
-/0,27 |
-/-0,22 |
-/-0,25 |
-/-0,16 |
-/-0,24 |
-/-0,04 |
|
6 |
Torgay – Peski Tusum |
-/-0,28 |
-/-0,16 |
-/-0,28 |
-/-0,04 |
-/-0,17 |
-/0,06 |
Correlation analysis showed that, on average across all gauging stations in the region, air temperature exerts a more consistent and pronounced influence on low-flow discharge than precipitation. This is most evident in the winter period, where correlation coefficients between air temperature and discharge reach R = 0.40–0.68, indicating a moderate and statistically significant relationship. Particularly important is the influence of cold-season temperature, reflecting soil freezing conditions, infiltration feeding, and the persistence of groundwater storage (Nasiyev et al., 2015; Shayakhmetova et al., 2023).
However, a number of gauging stations exhibited local deviations from the general pattern, which in turn highlights the complexity of interactions between climatic and hydrological factors. For example, at the Tobyl River – Kostanay station, a positive correlation was found between warm-season temperature and minimum summer discharge. This can be explained by the specific features of runoff formation at this site: earlier snowmelt, infiltration recharge, and a significant contribution of groundwater inflow, driven by the geomorphological and hydro-engineering characteristics of the catchment.
Therefore, summer discharge is primarily determined by precipitation and its interaction with evaporation and infiltration processes, whereas winter discharge is mostly influenced by the temperature regime. The identified heterogeneity of climatic sensitivity therefore requires a spatially differentiated approach to the analysis of river regimes under changing climate conditions.
Correlation analysis also showed that in the basins of the Karatorgay and Torgay rivers, a negative relationship is observed between total precipitation and minimum river discharge, both in summer and in winter. For the Karatorgay River during the summer season, correlation coefficients between precipitation (annual, cold season, and warm season) and minimum discharge ranged from R = –0.16 to –0.24, indicating a weak but persistent negative relationship. For the Torgay River, the relationship is weaker, approaching zero (R ≈ –0.036), which likewise points to the absence of a direct dependence between precipitation volume and the formation of minimum discharge.
During the winter period, river flow in these catchments often ceases completely due to deep channel freezing. Nevertheless, even under conditions of strongly limited water movement, a negative correlation was recorded between winter precipitation and discharge, which requires further detailed investigation (Ahn et al., 2025; Castellano-Rioja, 2025; Cavero & Ferraz, 2025; Hart & Reed, 2025; Khalil & Nassar, 2025; Rani & Gehrke, 2025; Rivas & Molina, 2025; Rojas & Paredes, 2025; Shen & Bao, 2025; Silik, 2025).
The basins of the Kara-Torgai and Torgay rivers are located in an arid climate zone, characterized by low annual precipitation and high evaporation. The geological structure is dominated by sandy loams, sands, and carbonate deposits with high permeability. Under such conditions, most atmospheric precipitation infiltrates into deep soil horizons, accumulates in local depressions, and rapidly evaporates under high summer air temperatures. As a result, even heavy rainfall does not lead to an increase in channel runoff; on the contrary, it reduces its volume because the water balance shifts toward irreversible losses (Andersson et al., 2001; Scanlon et al., 2006). Minimum river discharge, particularly in arid basins, should be regarded as an inertial component of the water regime, shaped by long-term accumulation, infiltration, and groundwater recharge, rather than by an immediate response to current precipitation. As noted in a number of international studies (Shanafield, 2021; Birimbayeva et al., 2024), under arid climate conditions, there is often a time lag between precipitation events and river discharge. In summer, intense rainfall is accompanied by high evaporation and deep infiltration, so actual runoff may not increase and in some cases may even decline due to falling groundwater levels and downward drainage of soil moisture (Goodrich, 2004).
In continental climates with prolonged sub-zero temperatures in winter, small rivers often freeze to the channel bed, effectively halting water movement. Precipitation falling as snow accumulates in the snowpack but does not contribute to current discharge formation, which makes the correlation between winter precipitation and runoff not only weak but, in some cases, negative. Moreover, snow accumulation can aggravate channel blockages by ice crusts and icings, further restricting already limited winter flow (Woo, 2008).
Anthropogenic influences cannot be ruled out either. During low-water years, precipitation events may trigger intensified water use for agriculture (e.g., filling ponds, irrigation), leading to redistribution of water resources and a reduction in river discharge available downstream (Falkenmark & Rockström, 2006).
The negative relationship between precipitation and minimum discharge in the Karatorgay and Torgay basins can thus be explained by a combination of factors:
These findings highlight the importance of an integrated approach to analyzing river flow response to climatic signals, especially in arid and semi-arid regions where the classical “precipitation–runoff” relationship does not always operate in a straightforward positive direction (Wheater & Evans, 2009).
Climatic contrasts in dry and wet years
A comparative analysis of climate–hydrological relationships depending on interannual water availability revealed qualitatively different regimes of low-flow formation. Separate assessments of climatic controls in wet years (P ≤ 0.25 quantile) and dry years (P ≥ 0.75 quantile) demonstrated significant differences in the governing factors.
In dry years, the role of precipitation decreases, while the influence of temperature becomes more pronounced, especially the temperatures of the preceding warm season. For example, at the Toguzak River gauging station (winter flow), a positive correlation was identified between the previous warm-season temperature and low-flow discharge (R ≈ 0.5). This may reflect prolonged recharge of groundwater under high evaporative demand and subsequent transformation of the water balance.
In wet years, the role of precipitation mainly during the warm season increases. This reflects a more direct response of river flow to meteorological conditions under:
– sufficient soil moisture saturation;
– the presence of active surface drainage;
– reduced infiltration and evaporation losses.
However, for several rivers, such as the Torgai and Kara-Torgai, earlier studies assumed a strong positive correlation between summer precipitation and minimum flow (up to R ≈ 0.6). A later re-examination of the data, however, indicated that actual correlation coefficients remain at weak to moderate levels (R < 0.4), and in some cases even show a negative sign. This points to the nonlinearity and instability of precipitation–low flow relationships, even in water-abundant years, which may be linked to infiltration capacity of the catchment and landscape-related water storage processes (Scanlon et al., 2006).
In dry years, the role of precipitation decreases markedly, while temperature becomes the dominant factor. Of particular importance is the temperature of the preceding warm season, which affects soil moisture losses through evaporation, groundwater dynamics, and the degree of watershed desiccation.
For example, at the Toguzak hydrological station (winter discharge), a moderate positive correlation (R ≈ 0.5) was found between the mean temperature of the previous warm season and the minimum winter flow. This likely reflects enhanced infiltration and groundwater storage under high evaporative demand, which in turn supports a more stable subsurface contribution to winter discharge despite weak surface inflow (Shanafield, 2021).
Overall, in wet years, minimum discharge responds more directly to precipitation, especially summer rainfall if surface drainage remains active. In contrast, in dry years, climatic control shifts towards temperature, which governs potential moisture losses and redistribution within the catchment. This suggests that modeling and forecasting of low-flow conditions should differentiate between wet and dry years, accounting for seasonal water balance and the phase state of available resources.
These contrasts illustrate a dynamic nonstationarity of flow-controlling factors: in water-deficient years, the thermal regime and antecedent moisture storage are critical, while in wetter years, the flow regime is driven mainly by the amount and timing of rainfall.
The findings highlight the need to incorporate climate-related variability into predictive models and adaptation strategies for water management. They are consistent with results reported in a number of international studies (Kundzewicz et al., 2008; Belyaev et al., 2019), which show an increasing influence of air temperature on river discharge under global warming, particularly where precipitation is uneven, and evaporative demand is high. Similar trends have been documented in other regions (Kundzewicz et al., 2008; Shiklomanov & Lammers, 2020; Birimbayeva et al., 2024; Makhmudova et al., 2024; Smagulov et al., 2025), where changes in the thermal regime are driving a reorganization of river flow (Agrawal et al., 2024; Bona & Lozano, 2024; Khan et al., 2024; Ha et al., 2025; Yilmaz & Erkol, 2025).
This study contributes to the regional understanding of climate control on river discharge under continental climate and steppe landscapes. A key finding is the spatial–seasonal variability of climate sensitivity, where the dominant factors, such as temperature or precipitation, change not only vary with the season but also with geographic location and catchment water availability (Danchin et al., 2024, Essah et al., 2024; Hillman, 2024; Qiao et al., 2024; Rivera & Carter, 2024; Snodin & McCrossen, 2024; Souza et al., 2024; Braun et al., 2025; Cuenca-Martínez et al., 2025; Eid et al., 2025).
Particularly noteworthy are the contrasts between dry and wet years, which are rarely addressed in international studies focusing mainly on averaged conditions. This approach allows for a more detailed understanding of hydrological responses to extreme climatic events and improves the realism of regional hydrological forecasts. The results provide an empirically based complement to global climate–hydrology models, highlighting the specific characteristics of steppe river systems in Central Asia and Eastern Europe.
The homogeneity of low-flow time series is critical for accurate statistical modeling, hydrological forecasts, and water resource assessment. Structural breaks may reflect changes in climate, hydrotechnical interventions, or anthropogenic impacts. Steppe and semi-arid regions of Central and Northern Kazakhstan are particularly sensitive, as runoff is episodic and strongly dependent on meteorological conditions (Pettitt, 1979; IPCC, 2021).
To assess homogeneity, the following methods were applied: parametric tests—Student’s t-test (for mean differences) and Fisher’s F-test (for variance differences); non-parametric tests—Mann–Whitney U-test (for differences in distributions). All calculations were performed using Python with the SciPy Stats package. The significance level was set at α = 0.05. The results indicate substantial differences in runoff structure before and after identified breakpoints:
The results confirm runoff heterogeneity across almost all studied rivers. Possible contributing factors include:
Similar findings regarding the impact of climatic fluctuations on runoff homogeneity are reported in previous studies (IPCC, 2021; Kundzewicz et al., 2008). In Central and Northern Kazakhstan, the homogeneity of river discharge series has been disrupted by the combined effects of climatic and anthropogenic factors. Statistical tests (Student’s t-test, Fisher’s F-test, Mann–Whitney U-test) confirmed the presence of structural breaks, particularly for rivers with seasonal or episodic flows. Future hydrological models should consider phase-based segmentation of time series according to breakpoints, modeling each regime separately, and taking into account winter low-flow degradation as an indicator of climate vulnerability.
Stationarity analysis of low-flow time series
Low-flow time series were tested for stationarity using two widely applied methods:
The two tests complement each other: ADF detects trends or drifts, while KPSS evaluates the stability of the level. The results are summarized in Table 3.
The summer low-flow series of the Toguzak River was found to be stationary according to both tests, allowing the use of standard modeling approaches without additional transformation. In contrast, the winter low-flow series of the same river was non-stationary by both tests, indicating a persistent trend and a possible structural shift.
In some cases, discrepancies between the tests were observed. For example, at the Tobyl River – Grishenka station (winter period), the series was stationary according to the ADF test but non-stationary according to KPSS, suggesting a phase shift in mean values while the autocorrelation structure remains largely intact.
Low-flow series for the Ayat, Karatorgay, and Torgay rivers also exhibited unstable behavior, with ADF failing to reject the null hypothesis of non-stationarity while KPSS confirmed it. This highlights the need for caution when analyzing trends. Overall, most series, particularly for the winter period, show signs of non-stationarity, which should be taken into account when developing predictive models, assessing water availability, or calculating normative flows. Under such conditions, the use of adaptive models or prior transformation of the series (e.g., differencing) is justified.
Table 3. Stationarity of minimum flow time series in the Tobyl–Torgay river basin
|
№ |
River – Observation point |
Season |
ADF statistic |
p-ADF |
KPSS statistic |
p-KPSS |
Interpretation |
|
1 |
Toguzak – s. Togyzak
|
Summer |
–3.57 |
0.0063 |
0.1726 |
0.1 |
The time series is stationary by both tests. |
|
2 |
Winter |
–1.36 |
0.6014 |
0.7504 |
0.01 |
The time series is non-stationary by both tests. |
|
|
3 |
Tobyl – s. Grishenka
|
Summer |
–2.94 |
0.0403 |
0.2556 |
0.1 |
The time series is stationary by both tests. |
|
4 |
Winter |
–4.11 |
0.00094 |
0.5579 |
0.0286 |
Stationary by ADF, non-stationary by KPSS |
|
|
5 |
Tobyl – Kostanay |
Summer |
–2.65 |
0.0837 |
0.5388 |
0.0329 |
The time series is non-stationary by both tests. |
|
6 |
Winter |
–7.17 |
0.0000003 |
0.7224 |
0.01 |
The time series is non-stationary by both tests. |
|
|
7 |
Ayat – s. Varvarinka |
Summer |
–4.21 |
0.00063 |
0.2831 |
0.1 |
The time series is non-stationary by both tests. |
|
8 |
Kara-Torgay – s. Urpek |
Summer |
–4.14 |
0.00081 |
0.6851 |
0.015 |
The time series is non-stationary by both tests. |
|
9 |
Torgay – Peski Tusum |
Summer |
–1.94 |
0.312 |
0.7078 |
0.0128 |
The time series is non-stationary by both tests. |
To evaluate the direction of changes in low flows, a linear trend analysis based on 10-year moving averages was performed (Figure 3). This approach allows the identification of persistent trends by smoothing short-term fluctuations and highlighting long-term changes in the hydrological regime. The analysis was conducted separately for summer and winter periods.
For the summer low flow of the Toguzak River – Togyzak station, the 10-year moving average analysis revealed a slight positive change. The average growth rate was approximately 0.0003 m³/s per year. However, the coefficient of determination (R² = 0.026) indicates very low explanatory power and statistical reliability. Therefore, the observed trend is not statistically robust and likely reflects short-term fluctuations rather than a systematic change in river regime.
For the summer period, a more pronounced positive trend was detected, with an increase of about 0.0035 m³/s per year. Despite the relatively low R² (~0.07), this indicates a consistent increase in discharge, possibly linked to the recurrence of warm-season precipitation and local catchment characteristics.
At this site, a moderate increase in summer low flow was recorded, with an average growth rate of 0.04 m³/s per year. The coefficient of determination exceeds 31 %, confirming the presence of a stable trend. The observed behavior is likely influenced by reservoir regulation, which helps maintain minimum flows during summer by smoothing seasonal fluctuations.
|
|
|
|
(а) Toguzak – s. Togyzak |
|
|
|
|
|
(b) Tobyl – s. Grishenka |
|
|
|
|
|
(с) Tobyl – Kostanay |
|
|
Figure 3. Dynamics of minimum flow changes in the rivers of the Tobyl basin |
|
The winter low flow of the Toguzak River – Togyzak station shows a pronounced and consistent increase, with an average growth rate of 0.063 m³/s per year. The coefficient of determination is nearly 40 %, indicating a strong statistical basis. These changes can be interpreted as resulting from climate warming, reduced riverbed freezing, and enhanced groundwater contributions.
Analysis of the Tobyl River – Grishenka station revealed a steady increase in winter low flow at a rate of 0.046 m³/s per year. The coefficient of determination is 0.345, meaning that approximately 35 % of the variability in minimum flows can be explained by the linear trend, indicating moderate stability. The increase in winter low flows is likely associated with gradual warming, manifested as higher average air temperatures, reduced soil and riverbed freezing, and enhanced groundwater supply at the end of winter. Such processes are possible even under harsh climatic conditions due to the accumulated thermal background and extended transitional temperature periods.
The highest trend stability among all studied sites is observed at the Tobyl River – Kostanay section. Winter low flow increases on average by 0.048 m³/s per year, with a coefficient of determination approaching 50 %, confirming a systemic shift to a new hydrological regime. In addition to the climatic factor, reservoir regulation significantly affects the flow regime, smoothing seasonal variability and maintaining minimum flows even during winter months. The results of this study complement existing scientific knowledge, emphasizing the need to consider regional features of climate–hydrology interactions when developing adaptive water management strategies under changing climate conditions.
CONCLUSION
Through an extensive examination of longitudinal data concerning minimum river flows alongside climate observations gathered from the Kostanay meteorological station, significant alterations in the thermal and hydrological patterns of northern Kazakhstan were discerned, particularly throughout the 1980s and 1990s. The most notable changes were detected in air temperatures during the cold season, which exhibit a consistent rise at a rate reaching +0.049 °C annually, corroborated by elevated R² values and trends that are statistically significant (p<0.01).
The analysis of climate–hydrology relations identified minimum flows in winter as mainly controlled by air temperature because it affects soil freezing depth, infiltration for meltwater, and recharge for groundwater. Precipitation is involved secondarily when it comes into play under moisture-accumulation-sufficient conditions, like wet winters.
For warm-period conditions in the summertime, warm-period atmospheric precipitation is the dominant controller of minimum flow, highlighting the importance of rainfall at the immediate surface for sustaining low flows. Temperature effects are less certain: positive (via increased snowmelt and recharge for groundwater stores) or negative (via increased evaporation and less efficient runoff). Interannual variations in climatic responsiveness in arid and wet years were identified as requiring consideration in integrated water resources management. A comparative analysis across interannual flow conditions indicated that the climatic influence on flow framework is characterized by dynamic non-stationarity: in arid years, temperature is dominant (explaining about 50–60%), while in wet years, precipitation and stored water reserves become increasingly important. That is why hydrological assessments need to consider spatiotemporal differentiation.
ACKNOWLEDGMENTS: None
CONFLICT OF INTEREST: None
FINANCIAL SUPPORT: This research was funded by the Committee of Science of the Ministry of Science and Higher Education of the Republic of Kazakhstan (Grant No. BR24992785 "Organization and implementation of comprehensive research to ensure sustainable development of the agro-industrial complex of the Kostanay region with the creation of a research and technology center").
ETHICS STATEMENT: None
Abdykadyrov, A., Marxuly, S., Mechsheryakova, T., Tastemirova, B., Tolen, G., Kulshikova, E., Duissenov, N., & Uzak, M. (2025). Scientific foundations of industrial engineering and geoecological utilization of water resources. Water Conservation and Management, 9(4), 707–713.
Agrawal, M., Shrivastava, S., Khare, R. L., Jaiswal, S., Singh, P., & Hishikar, R. (2024). Comparative analysis of nephrotoxicity in patients on tenofovir-containing and non-tenofovir antiretroviral therapy. Annals of Pharmacy Practice and Pharmacotherapy, 4, 24–33. doi:10.51847/ZTJxzvNZt4
Ahn, S., Kim, S., Zhang, H., Dobalian, A., & Slavich, G. M. (2025). The long-term effects of early-life adversity on health, lifestyle, and cognition. Bulletin of Pioneer Research in Medical and Clinical Sciences, 5(1), 13–27. doi:10.51847/6ZhsvRWk26
Andersson, L., Rosberg, J., & Pers, C. (2001). Hydrological impact of climate change in a small basin in Sweden. Climatic Change, 49(4), 349–366.
Belyaev, V. R., Gusarov, A. V., & Motovilov, Y. G. (2019). Climate-induced changes in the hydrological regime of rivers in Eastern Europe. Russian Meteorology and Hydrology, 44(8), 567–574.
Birimbayeva, L., Makhmudova, L., Alimkulov, S., Tursunova, A., Mussina, A., Tigkas, D., Beksultanova, Z., Rodrigo-Clavero, M.-E., & Rodrigo-Ilarri, J. (2024). Analysis of the spatiotemporal variability of hydrological drought regimes in the lowland rivers of Kazakhstan. Water, 16(16), 2316. doi:10.3390/w16162316
Bona, C., & Lozano, R. (2024). Delayed hematologic toxicities following axicabtagene ciloleucel and tisagenlecleucel CAR-T cell therapy: A retrospective analysis of lymphocytopenia, neutropenia, and thrombocytopenia. Asian Journal of Current Research in Clinical Cancer, 4(2), 1–4. doi:10.51847/13fIAq9Er0
Braun, D., Eckert, J., Lang, T., Schneider, L., & Wirth, F. (2025). High-resolution HLA-B allele and haplotype frequencies in the Qatari population: Implications for pre-emptive pharmacogenetic screening of drug hypersensitivity. Special Journal of Pharmacognosy, Phytochemistry and Biotechnology, 5, 77–87. doi:10.51847/43hsE3qngg
Burn, D. H., & Hag Elnur, M. A. (2002). Detection of hydrologic trends and variability. Journal of Hydrology, 255(1–4), 107–122. doi:10.1016/S0022-1694(01)00514-5
Castellano-Rioja, C. (2025). Investigating the effect of providing required training to mothers of children with surgery and its effect on mothers' anxiety. Journal of Integrative Nursing and Palliative Care, 6, 7–11. doi:10.51847/m0J08PS92
Cavero, M., & Ferraz, J. P. (2025). Sedentary behavior determinants in pregnant women during the second and third trimesters: Prospective findings from the large-scale Japan Environment and Children’s Study. Bulletin of Pioneer Research in Medical and Clinical Sciences, 5(1), 1–12. doi:10.51847/bdDIEuY7cl
Cuenca-Martínez, F., Herranz-Gómez, A., Madroñero-Miguel, B., Reina-Varona, Á., Touche, R. L., Angulo-Díaz-Parreño, S., Pardo-Montero, J., del Corral, T., & López-de-Uralde-Villanueva, I. (2025). A systematic review of the literature on the connection between cervical spine abnormalities and internal disorders of the temporomandibular joint. Journal of Current Research in Oral Surgery, 5, 1–10. doi:10.51847/e4CoCM6iSZ
Danchin, A., Ng, T. W., & Turinici, G. (2024). Transmission pathways and mitigation strategies for COVID-19. Interdisciplinary Research in Medical Sciences (Special Issue), 4(1), 1–10. doi:10.51847/p0YhQPxvkW
Dyachenko, N. A., & Maznev, N. P. (2017). Assessment of the impact of climatic factors on the minimum flow of rivers in Kazakhstan. Bulletin of Kazakh National University. Geographical Series, 2(63), 56–62.
Eid, M., Hassan, A., Abdelrahman, Y., & Saeed, O. (2025). Association of UGT1A1*28 and CYP2B6 c.516G>T variants with adverse effects in HIV patients receiving antiretroviral therapy in Chile. Special Journal of Pharmacognosy, Phytochemistry and Biotechnology, 5, 88–95. doi:10.51847/gHbHJPAf9Q
Essah, A., Igboemeka, C., & Hailemeskel, B. (2024). Exploring gabapentin as a treatment for pruritus: A survey of student perspectives. Annals of Pharmacy Education, Safety, Public Health and Advocacy, 4, 1–6. doi:10.51847/h8xgEJE3NE
Falkenmark, M., & Rockström, J. (2006). The new blue and green water paradigm: Breaking new ground for water resources planning and management. Journal of Water Resources Planning and Management, 132(3), 129–132. doi:10.1061/(ASCE)0733-9496(2006)132:3(129)
Goodrich, D. C. (2004). Variation in runoff and sediment yield due to vegetation change. Water Resources Research, 40(5), W05201.
Ha, M. N., Thanh, T. L., & Thanh, V. P. T. (2025). Factors influencing retail customer satisfaction in mobile banking services: A case study of Sacombank - Hanoi Branch. Annals of Organizational Culture, Leadership, and External Engagement Journal, 6, 41–50. doi:10.51847/0718r25rSv
Hart, O. B., & Reed, C. J. (2025). Anti-neuroinflammatory effects of 50% ethanolic Curcuma longa extract in LPS-stimulated BV2 microglia: Involvement of Nrf2/HO-1 signaling and inhibition of NF-κB/MAPK. Journal of Medical Sciences: Interdisciplinary Research, 5(1), 1–12. doi:10.51847/bMiGiIU1tk
Hillman, T. (2024). A predictive approach for determining the optimal anti-CCR5 monoclonal antibody. Interdisciplinary Research in Medical Sciences (Special Issue), 4(1), 11–22. doi:10.51847/jVJ3mRs3yR
Intergovernmental Panel on Climate Change. (2021). Climate change 2021: The physical science basis. Cambridge University Press. https://www.ipcc.ch/report/ar6/wg1/
Ivanov, A. A., & Shlykov, A. A. (2022). Changes in the minimum flow of rivers in southern Siberia under the influence of climate. Geography and Natural Resources, 3, 91–101.
Khalil, Y. H., & Nassar, O. S. (2025). Small-molecule inhibition of SOST suppresses breast cancer bone metastasis. Archive International Journal of Cancer and Allied Sciences, 5(1), 1–17. doi:10.51847/XR3tpAzZri
Khan, T. M., Tahir, H., Adil, Q., Baig, M. R., Jaber, A. A. S., Khaliel, A. M., & Mohammed, Z. M. (2024). A three-decade overview of female-specific cancers in Malaysia: A thorough examination. Asian Journal of Current Research in Clinical Cancer, 4(2), 5–18. doi:10.51847/LIdazW7afN
Kundzewicz, Z. W., & Robson, A. (2000). Detecting trend and other changes in hydrological data (WCDMP-45; WMO-TD No. 1013). World Meteorological Organization.
Kundzewicz, Z. W., Mata, L. J., Arnell, N. W., Döll, P., Jimenez, B., Miller, K., Oki, T., Şen, Z., & Shiklomanov, I. (2008). The implications of projected climate change for freshwater resources and their management. Hydrological Sciences Journal, 53(1), 3–10. doi:10.1623/hysj.53.1.3
Kwiatkowski, D., Phillips, P. C. B., Schmidt, P., & Shin, Y. (1992). Testing the null hypothesis of stationarity against the alternative of a unit root. Journal of Econometrics, 54(1–3), 159–178. doi:10.1016/0304-4076(92)90104-Y
Makhmudova, L., Moldakhmetov, M., Mussina, A., Kurmangazy, E., Kambarbekov, G., Zharylkassyn, A., & Zhulkainarova, M. (2024). Study of water stress in plains rivers: Climate and human influence. Evergreen Joint Journal of Novel Carbon Resource Sciences & Green Asia Strategy, 11(3), 1530–1543. doi:10.5109/7236809
Maznev, N. P. (2009). Variability of river flow and methods of its assessment. Nauka.
Nasiyev, B., Tulegenova, D., Zhanatalapov, N., Bekkaliev, A., & Shamsutdinov, Z. (2015). Studying the impact of grazing on the current state of grassland in the semi-desert zone. Biosciences Biotechnology Research Asia, 12(2), 1735–1742. doi:10.13005/bbra/1837
Nugmanov, A., Moldakhmetov, M., Makhmudova, L., Yskak, A., Chashkov, V., Kuanyshbayev, S., & Joldassov, A. (2025). Methodological foundations for assessing the water management balance of water bodies in the Tobyl River Basin. International Journal of Agriculture and Biosciences, 14(1), 31–39. doi:10.47278/journal.ijab/2024.1918
Pettitt, A. N. (1979). A non-parametric approach to the change-point problem. Journal of the Royal Statistical Society: Series C (Applied Statistics), 28(2), 126–135. doi:10.2307/2346729
Qiao, J., Luo, B., Ming, J., Zhou, S., Chen, Y., & Zhang, X. (2024). Prevalence and implications of non-prescription antibiotic dispensing in Baghdad community pharmacies. Annals of Pharmacy Practice and Pharmacotherapy, 4, 34–41. doi:10.51847/5SuGTfpren
Rani, N., & Gehrke, P. (2025). Promoting intercultural competence in German medical students via innovative medical ethics education focused on Muslim patients: A pilot study. Asian Journal of Ethics in Health and Medicine, 5, 1–12. doi:10.51847/0foncaeXr1
Rivas, D. S., & Molina, P. A. (2025). UBE2J1 suppresses colorectal cancer progression via ubiquitin-mediated degradation of RPS3. Archive International Journal of Cancer and Allied Sciences, 5(1), 18–39. doi:10.51847/7UbioAJLAO
Rivera, J., & Carter, E. (2024). Selective targeting of WNK kinases: High-throughput screening and structural insights into isoform specificity. Pharmaceutical Science and Drug Design, 4, 51–62. doi:10.51847/ilB2Ak8vDC
Rojas, V. M., & Paredes, D. E. (2025). Exploring acceptability, implementation, and contextual factors in a school-based physical activity intervention: Qualitative process evaluation from a feasibility trial in Yangzhou, China. Journal of Medical Sciences: Interdisciplinary Research, 5(1), 13–28. doi:10.51847/duZGj20T8V
Said, S. E., & Dickey, D. A. (1984). Testing for unit roots in autoregressive-moving average models of unknown order. Biometrika, 71(3), 599–607. doi:10.1093/biomet/71.3.599
Scanlon, B. R., Keese, K. E., Flint, A. L., Flint, L. E., Gaye, C. B., Edmunds, W. M., & Simmers, I. (2006). Global synthesis of groundwater recharge in semiarid and arid regions. Hydrological Processes, 20(15), 3335–3370. doi:10.1002/hyp.6335
Shanafield, M. (2021). Understanding episodic recharge and its impact on arid zone hydrology. Hydrogeology Journal, 29, 49–66.
Shayakhmetova, S., Savenkova, I., Nasiyev, B., Akhmetov, A., Useinov, A., Taskulova, A., & Temirbulatova, S. (2023). Agrotechnology for feed cultivation and creation of hayfields and pastures in the forest and steppe zone of Northern Kazakhstan. Sabrao Journal of Breeding and Genetics, 55(4), 1245–1258. doi:10.54910/sabrao2023.55.4.18
Shen, F., & Bao, L. (2025). Studying the effects of music on the time to gain independent oral feeding in premature infants. Journal of Integrative Nursing and Palliative Care, 6, 1–6. doi:10.51847/xBTC4CiH10
Shiklomanov, A. I., & Lammers, R. B. (2020). Changing discharge patterns of the Russian Arctic rivers. Hydrological Processes, 34(23), 4317–4330. doi:10.1029/2006JG000352
Silik, O. E. (2025). Integrating artificial intelligence in oncology: Turkish medical oncologists’ knowledge, attitudes, and ethical perspectives. Asian Journal of Ethics in Health and Medicine, 5, 13–24. doi:10.51847/ObxojlDDEH
Smagulov, Z., Alimkulov, S., Makhmudova, L., Talipova, E., Zagidullina, A., Birimbayeva, L., & Sailaubek, A. (2025). Transformation of the seasonal distribution of river flow in the Zhaiyk–Caspian water basin under changing climate conditions. Journal of Water and Climate Change, 16(2), 400–429. doi:10.2166/wcc.2024.537
Snodin, D. J., & McCrossen, S. D. (2024). Regulatory considerations of pharmaceutical impurities with emphasis on genotoxic impurities. Pharmaceutical Science and Drug Design, 4, 1–15. doi:10.51847/ck2yogXhAS
Souza, J. S., Reis, E. A., Godman, B., Campbell, S. M., Meyer, J. C., Sena, L. W. P., & Godói, I. P. D. (2024). Designing a healthcare utilization index to enable worldwide patient comparisons: A cross-sectional study. Annals of Pharmacy Education, Safety, Public Health and Advocacy, 4, 7–15. doi:10.51847/EeWKtBkVgK
Trenberth, K. E., Dai, A., van der Schrier, G., Jones, P. D., Barichivich, J., Briffa, K. R., & Sheffield, J. (2014). Global warming and changes in drought. Nature Climate Change, 4(1), 17–22. doi:10.1038/nclimate2067
Wheater, H., & Evans, E. (2009). Land use, water management, and future flood risk. Land Use Policy, 26, S251–S264. doi:10.1016/j.landusepol.2009.08.019
Woo, M. K. (2008). Cold region hydrology. In Cold region atmospheric and hydrologic studies. Springer.
Yilmaz, D., & Erkol, C. (2025). Understanding and reducing alienation among assembly line workers: Insights from leadership, compensation, and work formalisation. Annals of Organizational Culture, Leadership, and External Engagement Journal, 6, 56–72. doi:10.51847/FegO2CUoQF
Zhang, Y., Singh, V. P., & Wang, D. (2020). Impacts of climate change on streamflow regimes: A global assessment. Journal of Hydrology, 586, 124818.
This work is licensed under a Creative Commons Attribution 4.0 International License.