Forecasting the occurrence of hydrological drought according to a forecasting system is an important disaster reduction strategy. In this paper, a new drought prediction model adapted to changing environments was constructed. Taking the Luanhe River basin in China as an example, first, nonstationarity analysis of hydrological sequences in the basin was carried out. Then, conditional distribution models with the human activity factor as an exogenous variable were constructed to forecast hydrological drought based on meteorological drought, and the results were compared with the traditional normal distribution model and conditional distribution model. Finally, a scoring mechanism was applied to evaluate the performance of the three drought forecasting models. The results showed that the runoff series of the Luanhe River basin from 1961 to 2010 were nonstationary; moreover, when human activities were not considered, the hydrological drought class tended to be the same as the meteorological drought class. The calculation results of the models involving HI as an exogenous variable were significantly different from the models that did not consider human activities. When the current drought class tended towards less severe or normal, the meteorological drought tended to turn into more severe hydrological drought with the increase in human index values. According to the scores of the three drought forecasting models, the conditional distribution models involving the human index can further improve the forecasting accuracy of drought in the Luanhe River basin.

Typically, meteorological drought is regarded as the beginning of a drought event; after the occurrence of meteorological drought, other drought phenomena occur, such as hydrological drought (Fendeková et al., 2018; Fuentes et al., 2022; Wang et al., 2021). However, there is a delay period from meteorological drought to hydrological drought (Ding et al., 2021; Xu et al., 2019; Onyutha, 2017; Cammalleri and Vogt, 2018). Therefore, the occurrence of hydrological drought can be forecasted according to meteorological drought monitoring. Accurate hydrological forecast information is beneficial to reduce the losses caused by hydrological drought (Behzad and Hamid, 2019; Oertel et al., 2018 Dixit et al., 2022; Jehanzaib et al., 2020).

To identify the drought characteristics of the region, scholars have developed drought indices. For example, the standardized precipitation index (SPI) is typically used to identify and capture the characteristics of meteorological drought (McKee et al., 1993). Considering the influence of precipitation and temperature, Vicente-Serrano et al. (2010) proposed the standardized precipitation evapotranspiration index (SPEI) to characterize meteorological drought. The standardized runoff index (SRI), which focuses on the surface runoff of catchments, is typically used to indicate hydrological drought (Shukla and Wood, 2008). Aghelpour and Varshavian (2021) proposed the multivariate standardized precipitation index (MSPI) to forecast hydrological drought in Iran.

Statistical technology is an effective prediction method that has been widely used in drought forecasting in recent years (Alquraish et al., 2021; Abbasi et al., 2021; Dehghani et al., 2019; Mallya et al., 2013). For instance, neural network models have been proposed to combine multiple data for drought prediction (Rezaeianzadeh et al., 2016; Mokhtarzad et al., 2017; Ahnadi et al., 2011), and time series models can be used to analyse the variation in time series such as rainfall and runoff to achieve drought prediction (Moghimi et al., 2020; Natsagdorj et al., 2021; Stojković et al., 2020). The conditional probability model was proposed by Cancelliere et al. (2007) and developed for drought forecasting by Bonaccorso et al. (2015). Bonaccorso et al. (2015) showed that the conditional probability model can calculate the transition probabilities from the current drought index values to the future drought classes, and this is a more robust method that can be used to forecast drought than the traditional probability prediction models (such as the multivariate normal distribution model and Markov method).

A change in the environment may lead to the nonstationarity of the relationship between hydrological series (for example, precipitation and runoff series), which also occurs in the Luanhe River basin (Wang et al., 2016, 2018, 2020; Li et al., 2015). Traditional drought prediction methods need to be further improved to adapt to nonstationary conditions (Wang et al., 2022; Zhao et al., 2018; Chen et al., 2021). Ren et al. (2017) found that the conditional distribution model using large-scale climatic indices as covariates can improve the accuracy of meteorological drought forecasting in the Luanhe River basin. Although some progress has been made in the study of drought forecasting, there are relatively few studies considering the impact of the changing environment.

In this paper, to analyse the impact of human activities on hydrological drought, we constructed the human activity index (HI) based on the restoration method. Subsequently, conditional distribution models with the HI as the exogenous variable were developed to forecast hydrological drought based on meteorological drought, and then the results were compared with the traditional normal distribution model and conditional distribution model; as a result, the impact of human activities on transition probabilities was illustrated. A scoring mechanism was applied to the evaluation of the three probability models.

In addition to the introduction, this paper contains the following sections. Section 2 introduces the study area and data. Section 3 briefly describes the methods used in the research. Section 4 introduces the model construction and calculation results and analyses the results. Section 5 presents the prospects.

The Luanhe River basin, located in the subtropical monsoon region, covers an area of approximately 33 700 km

The geographical location of the Luanhe River basin.

Influenced by topography, meteorology, hydrology and hydrogeological conditions, the spatial distribution of groundwater resources in the Luanhe River basin is quite different. The recharge and storage conditions of shallow groundwater in plain areas and intermountain basins are relatively superior, and the content of groundwater in mountainous areas is relatively small (the area of mountainous areas in the Luanhe River basin accounts for 98.2 %). Therefore, the total amount of water resources in the Luanhe River basin is mainly considered to be affected by the amount of surface water resources.

In this paper, the monthly rainfall data from 26 stations in the Luanhe River basin from 1961 to 2010 were provided by the Hebei Provincial Hydrology and Water Resources Investigation Bureau. The average monthly rainfall data of the area were obtained by the inverse distance weighting interpolation method. The runoff data from 1961 to 2010 came from the inflow runoff series of the Panjiakou Reservoir. The SPI and SRI can be calculated for 1-month, 3-month, 6-month, and 12-month timescales to characterize meteorological drought and hydrological drought based on these data.

In the case of environmental changes, nonstationarity may occur in
hydrological series. The Pettitt test, as one of the important methods to
test whether there is nonstationarity in time series, can identify whether
there are change points in the sample series (Malede et al., 2022). Assuming
that the sample sequence is

The sliding

The Spearman correlation test can be applied to test the trend of time series, and the specific description refers to the article of Bishara and Hittner (2012).

The rainfall and runoff series of the watershed are usually strongly correlated. However, under the interference of human activities, the relationship between rainfall and runoff changes.

The double cumulative curve method can test the nonstationarity of the bivariate correlation between rainfall series and runoff series, and the point where the underlying surface is significantly altered by human activities can be determined according to the position of the slope change of the curve. Due to the short data series before and after the change point (20 years before the change point and 30 years after the change point), a linear equation was used to fit the relationship between precipitation and runoff.

The linear regression relationship of the cumulative rainfall and runoff
series can be calculated according to the following formula:

Human activities are the main reason for the nonstationarity of the runoff
series in the watershed, so the HI can be constructed to quantify the impact of human activities on runoff. Based on the linear regression relationship established between the accumulated precipitation and the accumulated runoff before the change point, the theoretical runoff sequence during the human activity period can be calculated from the measured precipitation sequence. SRI

The SPI is one of the important indicators for evaluating meteorological drought in the basin, and the SRI is an important indicator for evaluating hydrological drought in the basin. According to the rainfall data and runoff data in the basin, the SPI and SRI can be calculated at different timescales. Table 1 provides the drought class classification and corresponding SPI values and SRI values (Kolachian and Saghafian, 2021).

Drought class classification and corresponding SPI values and SRI values.

As a traditional drought class forecasting model, the multivariate normal
distribution model (Model 1) can forecast the future SRI class according to
the current SPI class. Assuming that the current SPI and SRI series both
satisfy a multivariable normal distribution, the joint probability density
can be expressed as follows (Chang et al., 2022):

Furthermore, according to the joint probability density function of the SPI
value

The conditional distribution model (Model 2) proposed by Bonaccorso et al. (2015) is described as follows: when one group of sample data

When sequence

Then, the probability of the current SPI value transitioning to the future
SRI drought class can be deduced as follows (Ren et al., 2017):

The calculation of the correlation coefficient

According to the above conditional probability model, when considering the HI as an exogenous variable, the model (Model 3) can be extended as follows:

A scoring mechanism was applied to evaluate the performance of the drought
forecasting models. In this method, the monthly drought transition
probability is summed to evaluate the model (Chen et al., 2013), where

In this paper, the area average monthly rainfall data of the Luanhe River
basin from 1961 to 2010 were obtained by spatial interpolation. The runoff
data came from the inflow runoff series of the Panjiakou Reservoir. Given
the significance level

The change points of the runoff series.

Figure 2a shows that the years of possible runoff change were 1979, 1996,
1997, 1998 and 1999. The

The results of the Spearman correlation test (Table 2) indicate that the runoff series showed an upwards trend before 1979, but the trend was not significant. However, there was a significant downwards trend in the series after 1979. In general, the runoff series showed a significant downwards trend.

Spearman correlation test results of runoff series trend.

In addition, according to historical records, local human activities (such as land use change and reservoir construction) are regarded as the main factors influencing runoff (Yan et al., 2018; Chen et al., 2021). Synthesizing the above analysis, 1979 was determined as the change point for the runoff sequence in the basin, and this conclusion was consistent with Li et al. (2015) and Wang et al. (2015).

According to the normality test results of rainfall and runoff series, it was reasonable to apply the conditional distribution model. To analyse the influence of different timescales of the SPI on the transition probabilities, using the forecast period as 1 month and the SPI timescales at 1, 3, 6 and 12 months as examples, the probabilities of converting SPI values to SRI classes were calculated (Fig. 3).

Influence of the SPI timescale on transition probabilities (

As shown in Fig. 3, when meteorological drought is categorized as extreme drought, the probabilities of maintaining the SRI class in extreme drought increased with the increasing SPI timescale. While the SPI had a short timescale, the response of the future SRI class to rainfall was fast, so the hydrological drought was more likely to tend to a normal status. This situation also occurred when the current meteorological drought was in another status.

In addition, the transition probabilities of drought were distinct for
different forecast periods. As seen in Fig. 4, when the forecast periods
were short (

Influence of forecast period on transition probabilities (

The effects of human activities are complex. To quantify the impact of human activities, the change point was identified, and then it was believed that the difference in the relationship between precipitation and runoff before and after the change point was caused by human activities. Moreover, the HI is easy to calculate and can approximately replace the influence of human activities. According to historical records, local human activities (such as land use change and reservoir construction) were regarded as the main factors influencing runoff (Yan et al., 2018; Chen et al., 2021). According to the above nonstationarity test results, 1979 was the change point, and the linear regression relationship of the cumulative rainfall and runoff series before and after the change point was established. The calculation results are shown in Table 3.

Linear regression relationship between cumulative precipitation
(

The HI results for different timescales are shown in Fig. 5.

Different average periods of the HI (HI-1: HI with a 1-month timescale; HI-3: HI with a 3-month timescale; HI-6: HI with a 6-month timescale; HI-12: HI with a 12-month timescale).

As shown in Fig. 5, the HI at all monthly scales generally ranged upwards, which means that human activities have intensified the occurrence of hydrological drought. According to historical statistics, many water conservancy projects were built in the basin from 1980 to 2000, and the construction and operation of large reservoirs in the mid-1990s may be the main reason for the serious negative values of the HI.

The HIs of different monthly scales were standardized, taking the 12-month timescale as an example, and the results were calculated as shown in Table 4.

HI-12 monthly mean and standard deviation (SD).

Furthermore, the drought transition probabilities involving the HI can be calculated from Eq. (23). Using the forecast period of 1 month from December and the SPI timescale of 12 months as an example, the drought transition probabilities from the current SPI values to the future SRI classes were calculated (Fig. 6). To analyse the effect of human activities on the drought transition probability more clearly, the calculation results of the three models are compared here separately. The horizontal coordinate indicates the drought classes corresponding to the SRI for the coming month, and the vertical coordinate is the drought transition probability.

Drought transition probability under the influence of human
activities (

In Fig. 6a, when the initial

From the initial

While the initial meteorological drought class was severe drought (Fig. 6c), the probabilities of the future SRI drought class being in the normal class became larger as the HI decreased. When the effect of human activities was not considered, the probability that the current SPI value transitioned to the SRI class under severe drought in the future month was the highest, and the probability of being in the normal class was the lowest. For Model 2, the probability of the SRI classes transitioning to severe drought was higher than the result of Model 1.

It was noteworthy that when the initial

In general, for the evaluation of drought transition probabilities in the
future month, hydrological drought classes tended to be the same as
meteorological drought when human activities were not considered, and this
situation was more significant in Model 2 than in Model 1. The calculation
results of the model involving the HI as an exogenous variable were
significantly different from those of the models that did not consider human
activities. The calculation results of Model 1 and Model 2 showed that the
future hydrological drought classes were more likely to be the same as the
meteorological drought classes in the current period, and they were more
significant in Model 2. In addition, it was obvious that the drought
transition probabilities of Model 3 were significantly different from those
of Model 1 and Model 2. Taking Fig. 6b as an example, when

To quantitatively evaluate the prediction accuracy of Model 1, Model 2 and
Model 3, the study period was divided into a correction period (1961–2003)
and a verification period (2004–2010), and then the drought transition
probability from the SPI value or class to the SRI class in the future

Model evaluation (Model 1: multivariate normal distribution model; Model 2: conditional distribution model; Model 3: conditional distribution model with the HI).

With the same timescale of the SPI, the model scores of Model 1 and Model 2
decreased as the forecast period

Many studies have noted that human activities have a significant impact on watershed runoff in the Luanhe River basin. In this paper, three probability models were constructed to calculate the transition probabilities from the current SPI classes or values to the future SRI classes; then, a scoring mechanism was applied to evaluate the performance of the models.

The calculation results of Model 1 and Model 2 showed that the future
hydrological drought classes were more likely to be the same as the
meteorological drought classes in the current period, and they were more
significant in Model 2. In addition, it was obvious that the drought
transition probabilities of Model 3 were significantly different from those
of Model 1 and Model 2. Under the condition of considering the HI, the results of the drought transition probability showed that when HI

Finally, a scoring mechanism was applied to the evaluation of the models, and the forecast results of the three models were evaluated. The results demonstrate that when the SPI timescale was the same, the scores of Model 1 and Model 2 decreased as the forecast period lengthened. In most cases, Model 2 performed better than Model 1, and the performance of Model 3 was the most stable of the three models and had the highest score. The conditional probability model considering the HI was more suitable for the Luanhe River basin, where human activities have a high influence.

Although this study has made some progress in the forecasting of hydrological drought in a changing environment, only the HI was considered as the exogenous variable in this paper, and human activities were generalized. In future studies, the HI can be analysed specifically. For example, the impact of land use and socioeconomics on drought prediction can be specifically analysed. In addition, climate factors can be further considered in future research.

The data and materials of the research are available from the Hydrology and Water Resource Survey Bureau of Hebei Province.

ML (first author and corresponding author): conceptualization, methodology, software, investigation, formal analysis, writing – original draft; MZ: data curation, writing – original draft; RC: visualization, investigation; YS: resources, supervision; XD: visualization, writing – review and editing. ML, MZ, RC, YS and XD all agreed to participate in the research for the article. ML, MZ, RC, YS and XD all agreed to publish this article.

The contact author has declared that none of the authors has any competing interests.

Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

We are grateful to the Hydrology and Water Resource Survey Bureau of Hebei Province for providing runoff data. This work was supported by the State Key Laboratory of Hydraulic Engineering Simulation and Safety (Tianjin University) program.

This research has been supported by the State Key Laboratory of Hydraulic Engineering Simulation and Safety (Tianjin University) (grant nos. HESS-2206 and HESS-2222).

This paper was edited by Brunella Bonaccorso and reviewed by Fawen Li and two anonymous referees.