2.1. SWAT Hydrological Model
The Soil and Water Assessment Tool (SWAT) is a distributed hydrological model developed by the USDA-ARS, based on physical mechanisms [
22,
23]. It incorporates various basin characteristics, such as climate, topography, geology, soil, land use types, and management practices, to effectively simulate the hydrological processes of translocation and transformation. The SWAT model has been widely adopted both domestically and internationally in watershed studies, demonstrating commendable simulation performance [
24,
25,
26,
27]. Given its proven track record, the SWAT model is selected to simulate the water cycle processes within the basin in this study.
The SWAT distributed hydrological model consists of three primary sub-modules: hydrological process, soil erosion, and water quality simulation. The selection of specific sub-modules depends on the research objectives. In this study, the hydrological simulation sub-module is employed to simulate the water cycle processes within the basin. The hydrological process sub-module of SWAT can be further divided into two components: the land surface part, which encompasses runoff generation and slope confluence, and the water surface part, which focuses on river confluence within the watershed. The former governs the input of water, sediment, nutrients, and chemicals into the main river channel within each sub-basin, while the latter determines the transport of water, sediment, and other substances from the river network to the basin outlet.
The watershed is initially divided into multiple sub-basins, based on the actual river network. These sub-basins are then further subdivided into smaller hydrological response units, using clustering techniques that consider the similarities in land use, soil type, and surface slope within each sub-basin. Through this process, the water yield within each sub-basin is obtained by effectively simulating the water cycle processes, including rainfall, evaporation, surface runoff, soil water dynamics, and groundwater interactions. The total runoff at the basin outlet is subsequently obtained by aggregating the individual runoff contributions from each sub-basin through the river network.
Within the hydrological simulation sub-module of the SWAT model, evapotranspiration is calculated using the Penman-Monteith method, which accounts for various climatic and surface characteristics. Surface runoff is simulated using the Soil Conservation Service (SCS) curve method, which provides a physically-based approach for runoff estimation considering different land use and soil types. The simulation of soil water dynamics employs a dynamic storage model, which incorporates factors such as slope, hydraulic conductivity, and temporal and spatial variations in soil water content. The simulation of groundwater includes two components: shallow groundwater and deep groundwater. Shallow groundwater represents water within the shallow saturation zone, which eventually contributes to river runoff as base flow. Deep groundwater refers to water within the deep pressure saturation zone.
(1) The water balance equation of shallow groundwater is as below:
Where, aqsh,j and aqdp,i-1 are the water storage in shallow water layer on day i and day i-1, respectively (unit, mm). wrchrg,sh is the recharge from shallow aquifer on day i (unit, mm). Qgw is the groundwater flowing into the main river on day i (unit, mm). wrecap is the amount of water entering the soil zone due to insufficient soil water on day i (unit, mm). wpump,sh is the groundwater extraction from shallow aquifer on day i (unit, mm).
(2) The water balance equation of deep groundwater is as below:
Where, aqdp,i and aqdp,i-1 are the water storage in deep aquifer on day i and day i -1, respectively (unit, mm). wdeep is the water infiltration from shallow aquifers into deep aquifers on day i (unit, mm). wpump,dp is the groundwater extraction from deep aquifer on day i (unit, mm).
The calculation formula of water yield in the hydrological response unit is as below:
Where, WYLD is the water yield (unit, mm). SURQ is the surface runoff (unit, mm). LATQ is the lateral flow (unit, mm). is the groundwater (unit, mm). TLOSS is the riverbed transmission water loss (unit, mm). PA is the water retention in ponds (unit, mm).
To evaluate the performance of the SWAT model in simulating monthly runoff, three statistical metrics were utilized: the Nash-Sutcliffe efficiency coefficient (Ns) [
28], the relative error (Re), and the coefficient of determination (R
2) [
29]. The simulated and observed runoff time series were judged to have an acceptable fit if the following thresholds were met concurrently on a monthly scale: Ns > 0.5, |Re| < 25%, and R
2 > 0.6 [
30].
2.3. Drought Propagation Time
Drought propagation time refers to the lag between the onset of meteorological drought and agricultural drought, or between agricultural drought and hydrological drought. The strongest correlation coefficient method has been widely utilized in recent studies to quantify propagation time [
40,
41,
42,
43]. Taking meteorological to agricultural drought as an example, the approach is illustrated in
Figure 2. Since crop growth cycles are generally less than one year, the cumulative drought scale is defined from 1-12 months. As shown in
Figure 2, let SPEI-j and SSMI-1 denote the j-month meteorological drought and 1-month agricultural drought in month i, respectively, where i = 1, 2,...,12. When the correlation between SPEI-j and SSMI-1 is maximized, the meteorological to agricultural drought propagation time in month i is j months. Similarly, the lag between SSMI-j and NHDI-1 represents the agricultural to hydrological drought time.
Seasonal mean propagation times are computed by averaging the monthly values - for instance, springtime lag is the mean of March, April, and May. Overall, this correlation-based approach quantitatively elucidates the cascading effects of drought across interlinked meteorological, agricultural, and hydrological systems.