the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Using normalised difference infrared index patterns to constrain semidistributed rainfall–runoff models in tropical nested catchments
Nutchanart Sriwongsitanon
Wasana Jandang
James Williams
Thienchart Suwawong
Ekkarin Maekan
Hubert H. G. Savenije
A parsimonious semidistributed rainfall–runoff model has been developed for flow prediction. In distribution, attention is paid to both the timing of the runoff and the heterogeneity of moisture storage capacities within subcatchments. This model is based on the lumped FLEXL model structure, which has proven its value in a wide range of catchments. To test the value of distribution, the gauged upper Ping catchment in Thailand has been divided into 32 subcatchments, which can be grouped into five gauged subcatchments at which internal performance is evaluated. To test the effect of timing, first the excess rainfall was calculated for each subcatchment, using the model structure of FLEXL. The excess rainfall was then routed to its outlet using the lag time from the storm to peak flow (TlagF) and the lag time of recharge from the root zone to the groundwater (TlagS), as a function of catchment size. Subsequently, the Muskingum equation was used to route subcatchment runoff to the downstream subcatchment, with the delay time parameter of the Muskingum equation being a function of channel length. Other model parameters of this semidistributed FLEXSD model were kept the same as in the calibrated FLEXL model of the entire upper Ping River basin (UPRB), controlled by station P.1 located at the centre of Chiang Mai province. The outcome of FLEXSD was compared to the (1) observations at the internal stations, (2) calibrated FLEXL model, and (3) the semidistributed URBS model – another established semidistributed rainfall–runoff model. FLEXSD showed better or similar performance during calibration and especially in validation. Subsequently, we tried to distribute the moisture storage capacity by constraining FLEXSD on patterns of the NDII (normalised difference infrared index). The readily available NDII appears to be a good proxy for moisture stress in the root zone during dry periods. The maximum moistureholding capacity in the root zone is assumed to be a function of the maximum seasonal range of NDII values and the annual average NDII values to construct two alternative models, namely FLEXSDNDII_{MaxMin} and FLEXSDNDII_{Avg}. The additional constraint on the moistureholding capacity (Sumax) by NDII, particularly in FLEXSDNDII_{Avg}, improved both the model performance and the realism of its distribution across the UPRB, which corresponds linearly to the percentage of evergreen forests (R^{2}=0.69). To check how well the models represents simulated root zone soil moisture (Su_{i}), the performance of the FLEXSDNDII models was compared to the time series of the soil water index (SWI). The correlation between the Su_{i} and the daily SWI appeared to be very good and was even better than the correlation with the NDII, which does not provide good estimates during wet periods. The SWI, which is modelbased, was not used for calibration but appeared to be an appropriate index for validation.
Runoff is one of the most important components of the hydrological cycle and can be monitored by the installation of a gauging station. Unfortunately, only a limited number of highquality gauging stations are available due to topographic, financial, and human resources limitations. A wide variety of rainfall–runoff models has been developed in gauged and ungauged catchments in different parts of the globe. Most rainfall–runoff models are categorised as lumped models, which provide runoff estimates only at the site of calibration. These models include FLEXL, FLEXTopo (Euser et al., 2015; Gao et al., 2014), NAM (Bao et al., 2011; Tingsanchali and Gautam, 2000; Vaitiekuniene, 2005; Yew Gan et al., 1997), SCS (Hawkins, 1993; Lewis et al., 2000; Mishra et al., 2005; Suresh Babu and Mishra, 2011; Yahya et al., 2010), and others. Among the wide range of existing lumped rainfall–runoff models, FLEXL has proven to be an adequate model for runoff estimation in a wide range of catchments (Fenicia et al., 2008, 2011; Gao et al., 2014; Kavetski and Fenicia, 2011; Tekleab et al., 2015). This model was further developed by Gharari et al. (2011) and Gao et al. (2016) to account for the spatial variability in landscape characteristics (FLEXTOPO), which is useful for prediction in ungauged basins (Savenije, 2010).
However, the generation of runoff exhibits the spatial and temporal variability in nature, which is not effectively accounted for by the conceptualisation of lumped models. Hence, the distribution of key routing and storage parameters is implemented in semidistributed frameworks to better understand their effects on the partitioning of moisture and water balance. The URBS model accounts for the spatiotemporal variability in rainfall by separating the catchment of interest into a series of subcatchments (Mapiam and Sriwongsitanon, 2009). This framework allows runoff to be estimated at any required upstream location (Carroll, 2004; Malone, 1999), thus providing useful application for realtime flood forecasting in a variety of catchments in Australia and globally (Malone, 2006; Malone et al., 2003; Mapiam and Sriwongsitanon, 2009; Mapiam et al., 2014; Rodriguez et al., 2005; Sriwongsitanon, 2010). But accounting for distributed routing and storage parameters alone does not address the variability in moisture storage capacities across heterogeneous subcatchments, which is a key parameter for runoff generation. For this reason, room for improvement with respect to the realism of conceptual hydrological models lies within their effectiveness at better capturing the variability in moisture stresses.
The fact that remote sensing (RS) proxies for moisture storage are available in ungauged basins makes RS an essential additional data source for distributed modelling, even though such proxies themselves have intrinsic modelbased uncertainties. Remote sensing observation techniques have been demonstrated in several studies to account for the spatial patterns of different vegetation types and moisture states, such that they could be valuable in constraining semidistributed hydrological models (e.g. Savenije and Hrachowitz, 2017).
The normalised difference infrared index (NDII) is an index that detects canopy water content (Hardisky et al., 1983), which has recently been investigated for monitoring drought conditions (Moricz et al., 2018; Xulu et al., 2018). The index is indicative of differences in moisture capacities and has been shown to correspond with root zone soil moisture (RZSM) dynamics. Sriwongsitanon et al. (2016) used NDII as a proxy for RZSM and showed its effectiveness in eight subcatchments of the upper Ping River basin in Thailand. This agrees with the study carried out by Castelli et al. (2019), who found reasonable correlations between Landsat 7 NDII values and the measured RZSM contents of rainfed olive trees growing in the arid regions of southeastern Tunisia. Mao and Liu (2019) found RZSM signatures to be wellcorrelated to NDII in most regions, except in river basins with high forest coverage, in addition to those with low moisture stress or those with trees intercepting deep groundwater (e.g. Eucalyptus). As such, the above studies reveal the worthwhileness of incorporating NDII to constrain semidistributed models to enhance model realism, as indicated by the achieved accuracy improvements to runoff estimation. This study utilises the fundamental model structure of FLEXL, which includes the same distributed time lags and channel routing as used in URBS and includes distributed root zone soil moisture capacity per subcatchment to create a new parsimonious semidistributed FLEX model for flood and flow monitoring within the (ungauged) subcatchments of the gauged upper Ping River basin. (1) To introduce the effect of runoff timing in a catchment with multiple subcatchments, travel times to the outfall of each individual subcatchment are computed on the basis of topographical indicators, and the routing of the discharge from the subcatchment outfall to stations further downstream is computed using the Muskingum method. These time lags are then applied in the FLEXSD model system and in the wellestablished URBS model for the purpose of comparison. These two semidistributed models only account for timing but not for the distribution of the moisture storage capacity, which is a crucial parameter in runoff generation. The distribution of time lags is expected to properly simulate the hydrograph shape, particularly the timing and shape of the peaks of all subcatchments within calibrated gauging station, but it does not affect the partitioning of the hydrological fluxes or the water balance and the accuracy of runoff estimates. (2) Subsequently, the effect of the distribution of the root zone moisture storage is studied in the FLEXSD model, making use of the spatial distribution pattern of the maximum and minimum range of NDII values and the annual average NDII values to construct two alternative models, namely FLEXSDNDII_{MaxMin} and FLEXSDNDII_{Avg}. (3) As a validation of the realism of the hydrological models tested, as indicated by their ability to capture the variability in internal moisture states across subcatchments, the derived simulated root zone moisture storages are compared to the independent time series of the modelbased soil water index (SWI) and the RSbased NDII.
2.1 Study area
The upper Ping River basin (UPRB) is situated between 17^{∘}14^{′}30^{′′} to 19^{∘}47^{′}52^{′′} N and 98^{∘}4^{′}30^{′′} to 99^{∘}22^{′}30^{′′} E in the provinces of Chiang Mai and Lamphun. The catchment area of the basin is approximately 25 370 km^{2}. The basin is dominated by wellforested, steep mountains, generally in a north–south alignment (Sriwongsitanon and Taesombat, 2011). The areal average annual rainfall and runoff of the basin from 2001–2016 is 1224 and 235 mm yr^{−1}, respectively. The land use for the UPRB in 2013 can be classified into six main classes comprising forest, irrigated agriculture, rainfed agriculture, bare land, waterbody, and others, which cover approximately 77.40 %, 3.11 %, 12.54 %, 1.99 %, 1.23 %, and 3.73 % of the catchment area, respectively (Land Development Department, LDD). The landform of the UPRB varies from an undulating to a rolling terrain, with steep hills at elevations of 1500–2000 m, and valleys of 330–500 m (Mapiam and Sriwongsitanon, 2009; Sriwongsitanon, 2010). The Chiang Dao district, north of Chiang Mai, is the origin of the Ping River, which flows downstream to the south to become the inflow of the Bhumibol Dam – a large dam with an active storage capacity of about 9.7 × 10^{9} m^{3} (Sriwongsitanon, 2010). The climate of the basin is dominated by tropical monsoons. The southwestern monsoon causes a rainy season between May and October, and the northeastern monsoon brings dry weather and low temperatures between November and April. Only 6142 km^{2} of the total area controlled by the runoff station P.1 (situated at the centre of Chiang Mai) is selected for this study (Fig. 1). The catchment area of the station P.1 is divided into 32 subcatchments (Fig. 1), which is where the semidistributed rainfall–runoff models are tested.
2.2 Rainfall data
Daily rainfall data from 48 nonautomatic rain gauge stations located within the UPRB and its surroundings from 2001–2016 were used in this study. These data are owned and operated by the Thai Meteorological Department and the Royal Irrigation Department. These data have been validated for their accuracy on a monthly basis using double mass curve, and some inaccurate data were removed from the time series before spatially averaging using an inverse distance square (IDS) to be applied as the forcing data of URBS, FLEXL, and FLEXSD. Mean areal rainfall depth for each of 32 subcatchments varies between 1100 (S17) and 1402 mm yr^{−1} (S11), as shown in Fig. 1b, while the average rainfall depth of P.1 is approximately 1224 mm yr^{−1}.
2.3 Runoff data
The Royal Irrigation Department (RID) operated seven daily runoff stations in the study area between 2001 and 2016, as shown in Fig. 1. Catchment P.56A was rejected from the study because it is located upstream of the Mae Ngat reservoir. Outflow data from the reservoir were used as input data in model calibration. Runoff data at the remaining six stations were used for the study, since they are not affected by large reservoirs. The data have been checked for their accuracy by comparing them with average rainfall data covering their catchment areas at the same periods. Table 1 presents the catchment characteristics and hydrological data for these six gauging stations in the UPRB. In this study, the catchments of these six stations were divided into 32 subcatchments (see Fig. 1), with areas ranging from 57 to 230 km^{2}. High variation in catchment size is due to the proximity between the locations of these runoff stations and the outlets of the tributaries. Runoff data have been checked for their accuracy by comparing the annual runoff coefficient between all stations. The comparison revealed that the runoff coefficients at P.20 in 2006 and 2011 are overestimated, while the runoff coefficient at P.21 in 2004 is underestimated, and in 2007 and 2009, the runoff coefficients are overestimated due to the incorrect rating curves (see Fig. 2). These inaccurate data would affect the results of the model calibration.
2.4 NDII data
The normalised difference infrared index (NDII) is a ratio of the nearinfrared (NIR) and shortwave infrared (SWIR) bands, centred at 859 and 1640 nm, respectively, as shown in Eq. (1). In this study, the NDII was calculated using the MODIS level 3 surface reflectance product (MOD09A1), which is available at 500 m resolution in an 8 d composite of the gridded level 2 surface reflectance products. Atmospheric correction has been carried out to improve the accuracy and can be downloaded from ftp://e4ftl01.cr. usgs.gov/MOLT (last access: 14 November 2020; Vermote et al., 2011). The 8 d NDII values between 2002–2016 were averaged over each of 31 subcatchments of the UPRB to be used for estimating the model parameter within subcatchment and for it to be compared to the 8 d average Su (root zone storage) values extracted from the model results at each station.
2.5 SWI data
The nearrealtime SWI is derived from the reprocessed surface soil moisture (SSM) data derived from the Advanced SCATterometer (ASCAT) sensor (Brocca et al., 2011; Paulik et al., 2014), which is a Cband scatterometer measuring at a frequency of 5.255 GHz in vertical transmit and receive polarisation (VV) (Paulik et al., 2014). The product makes use of a twolayer water balance model to describe the time series relationship between surface and profile soil moisture. This dataset of moisture conditions is available on a daily basis for eight characteristic time windows (1, 5, 10, 15, 20, 40, 60, and 100 d). The globalscale SWI dataset is available at 0.1^{∘}, which is about 10 km resolution, within 3 d after observation and can be downloaded from the Copernicus Global Land Service website. The dataset is available from January 2007 onwards. Since the SWI dataset is not complete in 2007, only the data between 2008 and 2016 have been used in this study. The characteristic time length is the only parameter in the SWI procedure. Bouaziz et al. (2020) specified the optimal T values (T_{opt}) by matching the modelled time series of RZSM from a calibrated FLEXL model to several SWI products in 16 contrasting catchments in the Meuse river basin. They concluded that the characteristic time lengths are differentiated amongst land cover (percent agriculture), soil properties (percent silt), and runoff signatures (flashiness index). In Sect. 5.4.2, the appropriate timescale for the Ping basin will be determined in 40 d.
3.1 FLEXL model
FLEXL is a lumped hydrological model comprising five reservoirs, namely a snow reservoir (Sw), an interception reservoir (Si), an unsaturated soil reservoir (Su), a fastresponse reservoir (Sf), and a slowresponse reservoir (Ss; Gao et al., 2014). Excess rainfall from a snow reservoir, an interception reservoir, and an unsaturated soil reservoir is divided and routed into a fastresponse reservoir and a slowresponse reservoir using two lag functions. It includes the lag time from storm to peak flow (TlagF) and the lag time of recharge from the root zone to the groundwater (TlagS). Each reservoir has process equations that connect the fluxes entering or leaving the storage compartment to the storage in the reservoirs (socalled constitutive functions; Sriwongsitanon et al., 2016). The water balance equations and constitutive equations for each conceptual reservoir are summarised in Fig. 3 and Table 2. The total number of model parameters is 11. Forcing data include the daily average rainfall and potential evaporation derived by the Penman–Monteith equation.
3.1.1 Snow reservoir
The snow routine, which is not very relevant in Thailand, can play an important role in areas with snow. When there is snow cover and the temperature (T_{i}) is above Tt, then the effective precipitation is equal to the sum of the rainfall (P_{i}) and snowmelt (M_{i}). The snowmelt (M_{i}) is calculated by the melted water per day per degree Celsius above Tt (F_{DD}; Eq. 2). The snow reservoir uses the water balance equation, Eq. (3), where Sw_{i} (mm) is the storage of the snow reservoir.
3.1.2 Interception reservoir
Interception is more important in summer and autumn. The interception evaporation (Ei_{i}) was calculated by the potential evaporation (Ep_{i}) and the storage in the interception reservoir (Si_{i}), with a daily maximum storage capacity (Imax; Eqs. 4 and 5). The interception reservoir uses the water balance equation, Eq. (6), as presented in Table 2.
3.1.3 Root zone reservoir
The root zone routine, which is the core of the hydrological models, determines the amount of runoff generation. In this study, we applied the widely used beta function of the Xinanjiang model (RenJun, 1992) to compute the runoff coefficient for each time step as a function of the relative soil moisture. In Eq. (7), Cr_{i} indicates the runoff coefficient, Su_{i} is the storage in the root zone reservoir, Sumax is the maximum moistureholding capacity in the root zone, and β is the parameter describing the spatial process heterogeneity of the runoff threshold in the catchment. In Eq. (8), Pe_{i} indicates the effective rainfall and snowmelt into the root zone routine, and Ru_{i} represents the generated flow during rainfall events. In Eq. (9), Su_{i}, Sumax, and potential evaporation (Ep_{i}) were used to determine the actual evaporation from the root zone Ea_{i}; Ce indicates the fraction of Sumax above which the actual evaporation is equal to potential evaporation, which is set here to 0.5, as previously suggested by Savenije (1997). Otherwise, Ea_{i} is constrained by the water available in Su_{i}. The unsaturated soil reservoir uses the water balance equation, Eq. (10), as presented in Table 2.
3.1.4 Fastresponse reservoir
In Eq. (11), Rf_{i} indicates the flow into the fastresponse routine, and D is a splitter to separate the recharge from preferential flow. Equations (12) and (13) were used to describe the lag time between storm and peak flow. Rf${}_{ti+\mathrm{1}}$ is the generated fast runoff in the unsaturated zone at time $ti+\mathrm{1}$, TlagF is a parameter which represents the time lag between storm and fastrunoff generation, c_{lagF}(i) is the weight of the flow in i−1 days before, and Rfl_{i} is the discharge into the fastresponse reservoir after convolution.
A linearresponse reservoir, representing a linear relationship between storage and release, was applied to conceptualise the discharge from the surface runoff reservoir, fastresponse reservoirs and slowresponse reservoirs. In Eq. (14), Qff_{i} is the surface runoff, with timescale Kff, which is active when the storage of the fastresponse reservoir exceeds the threshold Sfmax. In Eq. (15), Qf_{i} represents the fast runoff, Sf_{i} represents the storage state of the fastresponse reservoirs, and Kf is the timescales of the fast runoff. The fastresponse reservoir uses the water balance equation, Eq. (16), as presented in Table 2.
3.1.5 Slowresponse reservoir
In Eq. (17), Rs_{i} indicates the recharge of the groundwater reservoir. Equations (18) and (19) were used to describe the lag time of recharge from the root zone to the groundwater. Rs${}_{ti+\mathrm{1}}$ is the generated slow runoff in the groundwater zone at time $ti+\mathrm{1}$, TlagS is a parameter which represents the lag time of recharge from the root zone to the groundwater, c_{lagS}(i) is the weight of the flow in i−1 days before, and Rsl_{i} is the discharge into the slowresponse reservoir after convolution. In Eq. (20), Qs_{i} represents the slow runoff, Ss_{i} represents the storage state of the groundwater reservoir, and Ks is the timescale of the slow runoff. The slowresponse reservoir uses the water balance equation, Eq. (21), as presented in Table 2.
3.2 URBS model
URBS was developed by the Queensland Department of Natural Resources and Mines in 1990, based on the structures of RORB (Laurenson and Mein, 1990) and WBNM (Boyd et al., 1987). URBS is a semidistributed rainfall–runoff model that can provide runoff estimates not only at the calibrated station but also at the outlet of every subcatchment at any required location upstream. The calibrated catchment area needs to be divided into subcatchments to obtain different areal rainfall and different catchment and channel travelling time.
Table 3 presents five main processes used in URBS, comprising the calculation of the initial loss, proportional loss, excess rainfall, catchment routing, and channel routing. Excess rainfall is calculated separately between the pervious and impervious areas. For the pervious area, URBS assumes that there is a maximum initial loss rate (IL_{max}) to be reached before any rainfall can become the effective rainfall (${R}_{i}^{\mathrm{eff}}$). The initial loss (IL_{i}) can be recovered when the rainfall rate (R_{i}) is less than the recovering loss rate (rlr) per time interval (δt; see Eq. 22).
Excess rainfall for each time step is calculated using Eq. (23), by weighting the excess rainfall between pervious and impervious area using a ratio of the cumulative infiltration (F_{i}) and the maximum infiltration capacity (F_{max}). The recovering rate is included by simply reducing the amount infiltrated after every time step, using the reduction coefficient (k_{δt}), as shown in Eq. (24), and the pervious excess rainfall (${R}_{i}^{\mathrm{per}}$) is calculated using the Eq. (25), where pr is the proportional runoff coefficient. The remaining water $(\mathrm{1}\mathrm{pr}){R}_{i}^{\mathrm{eff}}$ will infiltrate to the root zone storage (dF_{i}; see Eq. 26). Excess rainfall is then routed to the centroid of any subcatchment using a nonlinear reservoir relationship (${S}_{i}=K{Q}_{i}^{m}$). The parameter m is the catchment nonlinearity, and K is the catchment travel time, which can be calculated for different subcatchment using the multiplication between the catchment lag time coefficient (β) and square root of each subcatchment area (A; see Eq. 27).
Thereafter, the outflow at the centroid of each subcatchment is routed along a reach downstream of each subcatchment, using the Muskingum equation (${S}_{i}^{\mathrm{ch}}={K}_{\mathrm{ch}}(X{I}_{i}+(\mathrm{1}X\left){Q}_{i}\right)$). The parameter X is the Muskingum coefficient, and K_{ch} is the channel travel time, which can be calculated for different subcatchment using the multiplication between the channel lag coefficient (α) and the reach length (L) between the closest location in the channel to the centroid and the outlet of each subcatchment (see Eq. 28).
4.1 Development of the semidistributed FLEX models
The first step in distribution is to account for the timing of floods and the rooting of flood waves as a function of topographical factors. The resulting semidistributed FLEXSD model is therefore expected to better represent the shape of hydrographs, although it would not affect the partitioning of fluxes or the water balance. The root zone storage capacity is a strong control on partitioning, affecting both runoff generation and evaporation. Therefore, the distribution of this parameter would potentially affect overall model performance more strongly than merely affecting the timing of the peaks. Therefore, in a second step, the NDII, as a proxy for moisture storage, is used to assess the distribution of moisture storage among subcatchments.
4.1.1 Accounting for distributed timing and channelrouting
FLEXSD is set up by applying lumped models for each subcatchment, which adds up to a semidistributed model for a downstream calibration site. Therefore, the catchment area of any gauging station needs to be divided into subcatchments. Runoff estimates at each subcatchment can be simulated using the structure of the original FLEXL, by calculating different excess rainfall for each subcatchment. The excess rainfall of each subcatchment is routed to its outlet using the lag time from rainfall to surface runoff (TlagF) and the lag time of recharge from the root zone to the groundwater (TlagS). In this study, TlagF and TlagS are calculated in hours instead of days to increase the model performance. The lag time is distributed among subcatchments using the following equations.
where TlagF and TlagS are lag time parameters for the entire catchment of a calibrated gauging station. The lag time of each subcatchment (TlagS_{sub}) is scaled by the square root of each subcatchment area divided by the overall catchment area (A).
Runoff estimates from an upstream subcatchment are later routed from their outlet to the outlet of a downstream subcatchment using the Muskingum method (Eq. 31) before adding to the runoff estimates of the downstream subcatchment.
where α and X are the delay time parameter and the channel routing parameter for the entire catchment, respectively. The delay time parameter of each subcatchment (K_{sub}) can be calculated by the multiplication between α and the main channel length of each subcatchment, as shown in Eq. (32).
4.1.2 Accounting for distributed root zone storage at subcatchment scale using the maximum and minimum values of NDII (FLEXSDNDII_{MaxMin} model)
The normalised difference infrared index (NDII) was used to estimate the root zone storage capacity for each subcatchment. The NDII values, which are available at 8 d intervals, were found to correlate well with the 8 d average root zone moisture content (Su) simulated by FLEXL during the dry period in eight subcatchments in the UPRB (Sriwongsitanon et al., 2016). The relation between NDII and Su can be described by an exponential function of the type ae^{b(NDII)}+c, with c being close to zero. The maximum value that Su can achieve is Sumax, which is the storage capacity of the root zone. The hypothesis is that the ecosystem creates sufficient storage to overcome a critical period of drought (Gao et al., 2014; Savenije and Hrachowitz, 2017). Every year has a maximum range of storage variation. If a sufficiently long NDII record is available, then the maximum of the annual ranges of the NDII should provide an estimate of the root zone storage capacity Sumax. By calibrating the hydrological FLEX model to discharge observations at the gauging stations, a Sumax value can be calibrated for each gauged catchment. This is a representative Sumax value for a particular gauging station, consisting of n subareas, as indicated by Sumax_{n}.
By using the NDII as proxy for root zone storage, we have developed the following equation for the proxy root zone storage capacity Sumax${}_{i}^{\prime}$ for a subarea within a river basin consisting of 31 subcatchments:
where Sumax${}_{i}^{\prime}$ is a scaled proxy for the root zone storage capacity of each subcatchment, and b is the remaining calibration parameter because the constant c and the factor a of the exponential function drop out. NDII_{i,max} and NDII_{i,min} represent the maximum and minimum values of NDII for each year of each subcatchment, while NDII_{n,max} and NDII_{n,min} indicate the maximum and minimum values of NDII for each year in the reference basin, which, in this case, is the entire upper Ping basin controlled by station P.1. The unscaled root zone storage capacity per subcatchment then becomes the following:
where Sumax_{n} is the calibrated value of the root zone storage capacity of the gauged catchment, and Sumax${}_{n}^{\prime}$is the areaweighted proxy for the root zone storage capacity.
4.1.3 Accounting for distributed root zone storage at subcatchment scale using average value of NDII (FLEXSDNDII_{Avg})
Instead of applying the maximum and minimum of the annual ranges of the NDII to distribute root zone storage at a subcatchment scale, we tested the annual average NDII value of each subarea to calculate Sumax${}_{i}^{\prime}$ as presented in the following equation.
where NDII_{i} represents the annual average NDII value of each subcatchment, while $({e}^{b\times {\mathrm{NDII}}_{i\to n}}{)}_{\mathrm{max}}$ and $({e}^{b\times {\mathrm{NDII}}_{i\to n}}{)}_{\mathrm{min}}$ indicate the maximum and minimum values of exponential function produced by the annual average NDII value within 32 subcatchments. The parameters b and R can be determined by model calibration. The parameter R is suggested to vary between 0.2 and 0.8 to force a scaled factor (Sumax${}_{i}^{\prime}$) to be more than 0 and less than 1. The average NDII value is supposed to reflect the maximum moisture storage capacity as well, since a high maximum value also leads to a higher average but is much easier to calculate. However, this method requires the introduction of the additional calibration parameter R.
4.2 Model applications
First, FLEXL and the four semidistributed models, URBS, FLEXSD, FLEXSDNDII_{MaxMin}, and FLEXSDNDII_{Avg} were individually calibrated (2001–2011) and validated (2012–2016) at the stations P.4A, P.20, P.21, P.75, P.67, and P.1 to provide baselines for comparison of model performances. Furthermore, additional runoff estimates at the outlet of 31 subcatchments, inclusive of the internal gauging stations upstream of P.1, were extracted from the four semidistributed models. The runoff estimates at the above stations were compared to the observed data and also for their comparability to the performance of FLEXL at individual stations.
The model parameters of the calibrated models were determined using the MOSCEMUA (MultiObjective Shuffled Complex Evolution Metropolis–University of Arizona) algorithm (Vrugt et al., 2003), by finding the Paretooptimal solutions defined by three objective functions of the Kling–Gupta efficiencies for high flows, low flows, and the flow duration (KGE_{E}, KGE_{L} and KGE_{F}), respectively. KGE_{E} is analysed using the following equations, where $\overline{X}$ is the average observed discharge, $\stackrel{\mathrm{\u203e}}{Y}$ is the average simulated discharge, S_{X} is the standard deviation of observed discharge, S_{Y} is the standard deviation of simulated discharge, and r is the linear correlation between observations and simulations. KGE_{L} can be calculated using the logarithm of flows to emphasise low flows. The Nash–Sutcliffe efficiency (NSE) is an independent statistical indicator, which is not utilised in the objective function but merely used to summarise model performance. The model calculates at daily time steps, but this is disaggregated to hourly to take the time lags into account. The output is again aggregated to daily time steps.
The results of this section are presented in Sect. 5.1 and 5.2.
Additionally, to test the efficacy of the distribution of moisture capacities by NDII, we investigated the relationship between the proportion of evergreen forests and the Sumax values produced by FLEXSDNDII_{Avg} and FLEXSDNDII_{MaxMin} in each subcatchment. The results are presented in Sect. 5.3.
4.3 Estimation of uncertainty in FLEXSDbased models
In hydrological models, inherent uncertainties caused by imperfect model structures and model parameters are unavoidable (Solomatine and Shrestha, 2009). To identify uncertainty in developed models, FLEXSD, FLEXSDNDII_{MaxMin}, and FLEXSDNDII_{Avg} were calibrated (2001–2011) and validated (2012–2016) at P.1 station using 50 000 random parameter sets. The 5 % bestperforming parameter sets were identified as feasible (Hulsman et al., 2020) and utilised to evaluate the uncertainties. The results of KGE_{E}, KGE_{L}, and KGE_{F} at the calibrated station (P.1) and at five upstream stations (P.20, P.4A, P.21, P.75, and P.67) were assessed.
4.4 Relationship between the average root zone soil moisture storage (Su_{i}) and the average NDII and SWI
Sriwongsitanon et al. (2016) suggested that NDII can be used as a proxy for soil moisture storage in hydrology. Here, 8 d average NDII values were compared to 8 d average root zone moisture storage (Su_{i}) derived at the six gauging stations, as calculated by FLEXL, FLEXSD, FLEXSDNDII_{MaxMin}, and FLEXSDNDII_{Avg}. The Su estimates were compared to the SWI product at the daily timescale. The Su–SWI relationship was examined for all characteristic time lengths to deduce the optimal time length which best represents RZSM dynamics in the UPRB. The results of Su–NDII and Su–SWI were aggregated to present these relationships at the seasonal basis. The coefficient of determination (R^{2}) and NSE were used as objective functions. Subsequently, for the FLEXSDbased models, their Su time series were extracted at the outlets of 31 subcatchments to be compared to time series of NDII and the optimal time length of SWI. These derived relationships serve as an indicator of the realism for the FLEXSDbased hydrological models examined.
5.1 Model performances
The assessment of model performances during the calibration and validation periods are presented in two subsections. Section 5.1.1 presents the results at internal gauging stations through the calibration/validation at all stations and from all models. Section 5.1.2 shows the results at internal gauging stations from the semidistributed models through calibration/validation at P.1, which were compared to the corresponding performance of FLEXL through calibration/validation at all stations. A discussion of the model parameters is provided in Sect. 5.1.3. Figures A1–A3 provide a comparison between performances undertaken at all stations and at P.1 only, in terms of accumulated flows, hydrographs, and duration curves, respectively. In the following sections, the term “average” refers to the average across all gauging stations (P.4A, P.20, P.67, P.75, P.21, and P.1).
5.1.1 Performance of all models by calibration/validation at all stations
For the calibration period (Fig. 4a), similar overall accuracy was produced by all models, with URBS and FLEXL showing NSE of 0.69 and 0.73, respectively, and the three FLEXSDbased models (FLEXSDNDII_{Avg}, FLEXSD, and FLEXSDNDII_{MaxMin}) yielding NSE of 0.76. However, FLEXL provided a slightly lower average KGE_{L} at some stations. For the validation period (Fig. 4b), a lower overall accuracy was achieved by FLEXL (NSE = 0.53) and URBS (NSE = 0.59) and with notably lower KGE_{L} at some stations. Meanwhile, the FLEXSDbased models showed NSE of 0.65–0.66. All models provided similar KGE_{E} and KGE_{F} during both periods.
In conclusion, all developed FLEXSDbased models can simulate runoff with similar accuracy or perform even better than the lumped model, FLEXL, and the established semidistributed model, URBS.
5.1.2 Performance of semidistributed models by calibration/validation at P.1
As previously mentioned, the performance of FLEXL at each gauging station provides a baseline comparison for the results of the semidistributed models in the predictive mode. Over the calibration period (Fig. 5a), the FLEXSDbased models produced similar accuracies (average NSE = 0.72–0.75), while URBS performed slightly worse (NSE = 0.68). This is comparable to, if not slightly better than, FLEXL (NSE = 0.73). However, in terms of the KGE_{E} indicator, FLEXL (KGE_{E}=0.86) performs notably better than the semidistributed models (KGE_{E}=0.78–0.82). For the validation period (Fig. 5b), the FLEXSDbased models (NSE = 0.67–0.70) and URBS (NSE = 0.65) performed notably better than FLEXL (NSE = 0.53).
As shown in Fig. A1, the semidistributed models are not capable of closing the water balance in four stations, except at the most downstream stations (P.67 and P.1). Additionally, the calibration/validation performed at all stations also closes the water balance better than achieved in predictive mode, although this may be due to overfitting. Nonetheless, the fact that the validation mode of all semidistributed models obtains more accurate results than the lumped and calibrated FLEXL model indicates a higher predictive capacity of the semidistributed models.
Needless to say, issues such as flow regulation and water withdrawal pose challenges to these semidistributed models. This is apparent at P.75, where outflows from Mae Ngat reservoir were potentially abstracted for agricultural demands but were unaccounted for, thereby causing notable overestimations of observed flows. This is reinforced in Fig. A3, which shows the lowest observed flows to be mostly below the modelled flows. P.4A also revealed overestimated flows, which drain a mountainous catchment with evergreen forest. In contrast, flow underestimations were seen in P.21 and P.20, which are intensively used catchments with rating problems (see Fig. 2). The lumped models are apparently not yet able to distinguish well between different landscapes. A landscapebased model, as suggested by Gharari et al. (2011) and Savenije (2010), could be the next step for improvement.
5.1.3 Discussion of model parameters
The model parameters used in FLEXL, FLEXSD, and FLEXSDNDII_{MaxMin} and FLEXSDNDII_{Avg} are summarised in Table A1. The FLEXSDbased models provide different values for TlagF (the time lag between storm and fastrunoff generation) and TlagS (the lag time of recharge from the root zone to the groundwater); the other parameters are kept the same as the calibrated values for P.1. Since TlagF and TlagS were designed to be related to the catchment area, the parameter values for each station are more reasonable compared to the values given by FLEXL. It can be noted that the values of TlagF obtained by FLEXSDNDII_{Avg} are much closer to the ones presented by FLEXSD compared to the values obtained by FLEXSDNDII_{MaxMin}.
In general, FLEXSD provides lower Sumax estimates than the other models, constraining evaporation in the dry season but compensating for this reduction with a smaller β value so as to limit excessive flood generation. Since these parameters jointly control Eq. (7), they can compensate for each other, leading to equifinality. If one of the parameters is constrained by additional information, as is the case here using the NDII, then this is no longer possible. The performance with respect to bestfit parameters may reduce in the process, but the model has gained realism and hence predictive power.
We see that the FLEXLSDNDII models show the highest realism (illustrated in Fig. A2) but not a very good performance in the subcatchment P.20, although they are still better than the other semidistributed models. P.20 remains a difficult subcatchment to predict, due to its flow regulation and water consumption. Also, we see that adding constraints to model calibration does not always improve bestfit performance, as compared to free calibration, but that realism can be improved.
5.2 Relationship between Sumax and percentage of evergreen forests at the subcatchment scale
The relationship between modelled Sumax values and the proportion of evergreen forests in each subcatchment has been presented in Fig. 6. R^{2} of 0.69 and 0.01 were yielded by FLEXSDNDII_{Avg} and FLEXSDNDII_{MaxMin}, respectively, indicating that Sumax values from the former model increase with greater vegetation coverage. With this unforeseen distinguishing power between subcatchments of varying land use compositions, it is worth investigating the merit of implementing FLEXSDNDII_{Avg} in the ungauged basins. To further test the realism of the models, the outputs of the models were compared to observations of NDII and the globalscale SWI dataset for verification, as described in Sect. 5.3.
5.3 Uncertainty in runoff estimation using FLEXSD, FLEXSDNDII_{MaxMin}, and FLEXSDNDII_{Avg}
Figure 7 displays the 2500 parameter sets selected to demonstrate the uncertainty (KGE_{E}, KGE_{F}, and KGE_{L}) at the six gauging stations. The results of all indicators are similar at P.67 and P.1. At Station P.75, KGE_{L} values are similar, while KGE_{E} and KGE_{F} provided by FLEXSDNDII_{MaxMin} are slightly higher than the others. At P.21, FLEXSD performed marginally better than the others for all indicators. The most notable observations are at P.4A and P.20, where FLEXSDNDII_{Avg} is shown to outperform the others in terms of KGE_{E} and KGE_{F}. This is reinforced by the observed and calculated hydrographs in Fig. A4, where this model shows the narrowest uncertainty band at these stations. Moreover, the uncertainty bands in the flow duration curves of Fig. A5 show that FLEXSDNDII_{Avg} produced the narrowest bands at P.4A and P.21, but FLEXSD showed the best performance at P.20 and P.67. In summary, all models show similar uncertainties; however, FLEXSDNDII_{Avg} reveals significantly better performance for upstream subcatchments and not much difference for other areas compared to the other two models.
5.4 Relationship between the average root zone soil moisture storage (Su_{i}) and the average NDII and SWI
5.4.1 Su–NDII relationships
Figure 8 shows the R^{2} and NSE during the wet and dry seasons for all six stations generated by three FLEXSD models. Figure 8a shows that the time series of NDII correlates well with Su values during the dry season by giving an average R^{2} value of 0.75–0.79. The average NSE values given by these models are 0.500.58, respectively. During the wet season, these correlations are lower (average R^{2}=0.41–0.46; NSE = 0.44–0.49).
The same procedure was also carried out for all 31 subcatchments and the results shown in Fig. 8b. During the dry season, the average R^{2}=0.71–0.74, but FLEXSDNDII_{Avg} provided notably higher NSE (0.52) than FLEXSD (0.41) and FLEXSDNDII_{MaxMin} (0.45). The relationships are significantly lower in the wet season (average R^{2}=0.36–0.41; average NSE = 0.36–0.40). All in all, FLEXSDNDII_{Avg} provided higher R^{2} and NSE than FLEXSDNDII_{MaxMin} in both seasons.
It is worth noting that FLEXSDNDII_{Avg} was able to differentiate the signatures in catchments of various soil moisture capacities. For a number of subcatchments, particularly those with more evergreen forest (e.g. S14 = 72.2 %, S23 = 58.8 %, and S29 = 44.7 %), the estimates of Su from FLEXSDNDII_{Avg} provide notably higher NSE than others. Figure 9 presents the time series of simulated root zone moisture storage (Su), NDII, and SWI (discussed in following section) for a selection of contrasting subcatchments. NDII signatures are well correlated with Su in subcatchments with low percentages of evergreen forest (S4, S6, S9, S15, and S32). Low correlations were found in evergreenforestrich subcatchments (S14, S23, and S29). The evergreen forest probably experiences less moisture stress compared to other land use land cover (LULC), so the situation NDII does not relate as well to simulated root zone soil moisture.
5.4.2 Su–SWI40 relationships
Through testing all available characteristic time lengths (TLs) of SWI, the TL of 40 d produced the optimal Su–SWI relationship across the six gauging stations. This was thereby used for subsequent investigation and referred to as SWI40. Figure 10 shows the R^{2} and NSE at all six stations simulated by the FLEXSDbased models. Figure 10a shows that SWI40 correlates well with Su values during the dry season (R^{2}=0.86–0.89 and NSE = 0.76–0.81). During the wet season, these correlations are of the same order of magnitude as in the dry season (R^{2}=0.87–0.89 and NSE = 0.80–0.84). The results for all 31 subcatchments are shown in Fig. 10b, which show average R^{2} values of 0.86–0.87 and NSE of 0.78–0.79 during the dry season. During the wet season, R^{2} values of 0.87 and NSE of 0.79–0.81 were yielded.
Despite the overall consistency in the average performance amongst these models and limited seasonality, it is evident that FLEXSDNDII_{Avg} shows the most variability across the 31 subcatchments. In addition to the sensitivity of NDII to the proportion of evergreen forests, Fig. 9 also reveals the sensitivity of SWI40 to it. This is reflected in the notable time lag of the Su time series from SWI40 in the productive subcatchments (e.g. S14, S23, and S29). This observation was not apparent with Su derived from FLEXSD or FLEXSDNDII_{MaxMin}.
Inspired by the significance of Sumax in governing the hydrological cycle, we investigated the utilisation of NDII to constrain our formulated FLEXSD model with the objective of distributing the parameter of interest across the UPRB's 31 subcatchments. The distribution of Sumax, particularly as achieved by FLEXSDNDII_{Avg}, helped produce betterinformed runoff estimates for the gauged (especially for stations P.4A and P.20, which are densely forested areas) and ungauged subcatchments than FLEXSD, FLEXSDNDII_{MaxMin}, and URBS. The analysis of model uncertainties reinforces the improvements to flow estimation in ungauged subcatchments. As explained below, this study has raised a few points worthy of discussion.

With the better representation of Su, the Su–NDII relationship yielded has become more informative about the underlying degree of aridity; that is, arid and productive subcatchments exhibit greater differences than presented by FLEXSD and FLEXSDNDII_{MaxMin}. The increased variation challenges the preconception from the study by Sriwongsitanon et al. (2016) that NDII is a suitable proxy for RZSM dynamics in the dry season, as it was hereby shown to barely correspond with Su_{i} in evergreen forests under any circumstances.

Upon testing the accuracy of derived Su_{i} against SWI, we realise the potential of implementing SWI to indicate moisture states of river basins. However, we should be reminded of its simplicity, such as using a single parameter, T, to emulate RZSM signatures. Furthermore, the model lacks the ability to account for the effects of local characteristics (e.g. LULC), which likely explains the variable Su–SWI correlation across the 31 subcatchments. Nonetheless, in future work, we are convinced that appropriately manipulating the T parameter to specific catchments, instead of using a T value of 40 d to represent the RZSM characteristics of the entire catchment as done in this study, could widen the potential of implementing SWI to indicate moisture states in river basins of contrasting characteristics.

It is worth noting the increased variability in Su–SWI across the 31 subcatchments (as seen in Fig. 10), in spite of the systematic increase in Su–NDII, which, counterintuitively, could be attributed to the improved realism of the FLEXSDNDII_{Avg}. That is, by using NDII to distribute Sumax, and thereby improving accuracy of estimated runoff, this accounts for the greater variability in Su signatures across the 31 subcatchments than was produced by FLEXSD.
Most lumped rainfall–runoff models are controlled by a gauging station at the outfall on which it is calibrated. Runoff estimation at any location upstream requires indirect approaches, such as model parameter transfer from gauged stations to ungauged locations, or applying relationships between model parameters and catchment characteristics to the ungauged locations. By using any of these approaches, the uncertainty in runoff estimation for ungauged catchments is unavoidable. A semidistributed hydrological model could offer a better alternative. Besides considering lag times and flood routing (as in FLEXSD), it has been shown that it is required to account for the spatial variation in the moistureholding capacity of the root zone. Therefore, the model was constrained by using NDII patterns (particularly average NDII, so as to produce the FLEXSDNDII_{Avg}) as a proxy for the spatial variation in root zone moisture leading to distributed Sumax values among subcatchments. The model parameters provided by the semidistributed FLEX models are more realistic compared to the original FLEXL, since they are distributed according to catchment characteristics comprising catchment area, reach length, and the NDII.
With the inclusion of NDII, the estimated catchmentscale root zone soil moisture (Su_{i}) has been shown to be increasingly sensitive to the underlying degree of aridity, which is arguably more representative of the hydrological responses across heterogeneous subcatchments. Such knowledge should be used in further studies to explore the opportunities of implementing the modelbased SWI to estimate soil moisture in different land use land cover.
The code is available at https://bit.ly/3ONm29s (Google Drive, 2023a).
The dataset is available at https://bit.ly/3N3sS9n (Google Drive, 2023b).
NS planned the research and wrote the manuscript draft; WJ simulated FLEX models and all analyses; JW edited the manuscript and prepared the figures; TS simulated URBS model; EM supported and analysed land use data; HHGS advised the research and edited the manuscript.
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.
The authors would like to express our sincere gratitude to the Royal Irrigation Department and Thai Meteorology Department for providing the hydrological data. We would also like to thank the reviewers for their invaluable feedback that led to the revision of this paper.
This research has been supported by the Master's Degree Research Assistantship Program, Faculty of Engineering, Kasetsart University, and Remote Sensing Research Centre for Water Resources Management (SENSWAT) under project number 62/03/WE.
This paper was edited by Roger Moussa and reviewed by Federico GómezDelgado and one anonymous referee.
Bao, A. M., Liu, H. L., Chen, X., and Pan, X. l.: The effect of estimating areal rainfall using selfsimilarity topography method on the simulation accuracy of runoff prediction, Hydrol. Process., 25, 3506–3512, https://doi.org/10.1002/hyp.8078, 2011.
Bouaziz, L. J. E., SteeleDunne, S. C., Schellekens, J., Weerts, A. H., Stam, J., Sprokkereef, E., Winsemius, H. H. C., Savenije, H. H. G., and Hrachowitz, M.: Improved understanding of the linkbetween catchmentscale vegetation accessible storage and satellitederivedSoil Water Index, Water Resour. Res., 56, e2019WR026365, https://doi.org/10.1029/2019WR026365, 2020.
Boyd, M. J., Bates, B. C., Pilgrim, D. H., and Cordery, I.: WBNM: A General Runoff Routing Model Computer Programs and User Guide, Water Research Laboratory, The University of New South Wales, https://doi.org/10.4225/53/57996b382f17b, 1987.
Brocca, L., Hasenauer, S., Lacava, T., Melone, F., Moramarco, T., Wagner, W., Dorigo, W., Matgen, P., MartínezFernandez, J., and Llorens, P.: Soil moisture estimation through ASCAT and AMSRE sensors: An intercomparison and validation study across Europe, Remote Sens. Environ., 115, 3390–3408, https://doi.org/10.1016/j.rse.2011.08.003, 2011.
Carroll, D.: URBS a Rainfall Runoff Routing Model for flood forecasting and design version 4.00, https://www.scribd.com/document/93746264/URBSManualV440 (last access: 15 January 2020), 2004.
Castelli, G., Oliveira, L. A. A., Abdelli, F., Dhaou, H., Bresci, E., and Ouessar, M.: Effect of traditional check dams (jessour) on soil and olive trees water status in Tunisia, Sci. Total Environ., 690, 226–236, https://doi.org/10.1016/j.scitotenv.2019.06.514, 2019.
Euser, T., Hrachowitz, M., Winsemius, H. C., and Savenije, H. H. G.: The effect of forcing and landscape distribution on performance and consistency of model structures, Hydrol. Process., 29, 3727–3743, https://doi.org/10.1002/hyp.10445, 2015.
Fenicia, F., Savenije, H. H. G., Matgen, P., and Pfister, L.: Understanding catchment behavior through stepwise model concept improvement, Water Resour. Res., 44, W01402, https://doi.org/10.1029/2006WR005563, 2008.
Fenicia, F., Kavetski, D., and Savenije, H. H. G.: Elements of a flexible approach for conceptual hydrological modeling: 1. Motivation and theoretical development, Water Resour. Res., 47, W11510, https://doi.org/10.1029/2010WR010174, 2011.
Gao, H., Hrachowitz, M., Schymanski, S. J., Fenicia, F., Sriwongsitanon, N., and Savenije, H. H. G.: Climate controls how ecosystems size the root zone storage capacity at catchment scale, Geophys. Res. Lett., 41, 7916–7923, https://doi.org/10.1002/2014GL061668, 2014.
Gao, H., Hrachowitz, M., Sriwongsitanon, N., Fenicia, F., Gharari, S., and Savenije, H. H. G.: Accounting for the influence of vegetation and landscape improves model transferability in a tropical savannah region, Water Resour. Res., 52, 7999–8022, https://doi.org/10.1002/2016WR019574, 2016.
Gharari, S., Hrachowitz, M., Fenicia, F., and Savenije, H. H. G.: Hydrological landscape classification: investigating the performance of HAND based landscape classifications in a central European mesoscale catchment, Hydrol. Earth Syst. Sci., 15, 3275–3291, https://doi.org/10.5194/hess1532752011, 2011.
Google Drive: Code, https://bit.ly/3ONm29s (last access: 3 June 2023), 2023a.
Google Drive:: Data, https://bit.ly/3N3sS9n (last access: 3 June 2023), 2023b.
Hardisky, M. A., Klemas, V., and Smart, R. M.: The influence of soil salinity, growth form, and leaf moisture on the spectral radiance of Spartina alterniflora canopies, Photogram. Eng. Remote Sens., 49, 77–83, 1983.
Hawkins, R. H.: Asymptotic Determination of Curve Numbers from Rainfall–Runoff Data, J. Irrig. Drain. Eng., 119, 67–76, https://doi.org/10.1061/(ASCE)07339437(1993)119:2(334), 1993.
Hulsman, P., Winsemius, H. C., Michailovsky, C. I., Savenije, H. H. G., and Hrachowitz, M.: Using altimetry observations combined with GRACE to select parameter sets of a hydrological model in a datascarce region, Hydrol. Earth Syst. Sci., 24, 3331–3359, https://doi.org/10.5194/hess2433312020, 2020.
Kavetski, D. and Fenicia, F.: Elements of a flexible approach for conceptual hydrological modeling: 2. Application and experimental insights, Water Resour. Res., 47, W11511, https://doi.org/10.1029/2011WR010748, 2011.
Laurenson, E. M. and Mein, R. G.: Version 4 Runoff Routing Program User Mannual, Department of Civil Engineering, Monash University, Autralia, https://library.deakin.edu.au/record=b1580310~S1 (last access: 20 January 2020), 1990.
Lewis, D., Singer, M. J., and Tate, K. W.: Applicability of SCS curve number method for a California oak woodlands watershed, J. Soil Water Conserv., 55, 226–230, 2000.
Malone, T.: Using URBS for Real Time Flood Modelling: in: 25th Hydrology & Water Resources Symposium, 2nd International Conference on Water Resources & Environment Research, Handbook and Proceedings, 1 January 1999, Institution of Engineers, Australia, Barton, ACT, 603–608, https://doi.org/10.3316/informit.ENG_70A, 1999.
Malone, T.: Roadmap mission for the development of a flood forecasting system for the Lower Mekong River, Technical Component – Main Report, Mekong River Commission Flood Management and Mitigation Programme, 152–156, http://archive.iwlearn.net/mrcmekong.org/download/free_download/AFF4/session4/Road_Map_Mission.pdf (last access: 18 Febuary 2020), 2006.
Malone, T., Johnston, A., Perkins, J., and Sooriyakumaran, S.: HYMODEL – a Real Time Flood Forecasting System, in: 28th International Hydrology and Water Resources Symposium: About Water, Symposium Proceedings, 1 January 2003, Institution of Engineers, Australia, Barton, ACT, 1.121–1.126, https://doi.org/10.3316/informit.350947309072083, 2003.
Mao, G. and Liu, J.: WAYS v1: a hydrological model for root zone water storage simulation on a global scale, Geosci. Model Dev., 12, 5267–5289, https://doi.org/10.5194/gmd1252672019, 2019.
Mapiam, P. and Sriwongsitanon, N.: Estimation of the URBS model parameters for flood estimation of ungauged catchments in the upper Ping river basin, Thailand, Science Asia, 35, 49–56, https://doi.org/10.2306/scienceasia15131874.2009.35.049, 2009.
Mapiam, P. P., Sharma, A., and Sriwongsitanon, N.: Defining the Z–R relationship using gauge rainfall with coarse temporal resolution: implications for flood forecasting, J. Hydrol. Eng., 19, 04014003, https://doi.org/10.1061/(ASCE)HE.19435584.0000616, 2014.
Mishra, S. K., Jain, M. K., Bhunya, P. K., and Singh, V. P.: Field applicability of the SCSCNbased Mishra–Singh general model and its variants, Water Resour. Manage., 19, 37–62, https://doi.org/10.1007/s1126900510763, 2005.
Moricz, N., Garamszegi, B., Rasztovits, E., Bidlo, A., Horvath, A., Jagicza, A., Illes, G., Vekerdy, Z., Somogyi, Z., and Galos, B.: Recent DroughtInduced Vitality Decline of Black Pine (Pinus nigra Arn.) in SouthWest Hungary – Is This DroughtResistant Species under Threat by Climate Change?, Forests, 9, 414, https://doi.org/10.3390/f9070414, 2018.
Paulik, C., Dorigo, W., Wagner, W., and Kidd, R.: Validation of the ASCAT Soil Water Index using in situ data from the International Soil Moisture Network, Int. J. Appl. Earth Obs. Geoinf., 30, 1–8, https://doi.org/10.1016/j.jag.2014.01.007, 2014.
RenJun, Z.: The Xinanjiang model applied in China, J. Hydrol., 135, 371–381, 1992.
Rodriguez, F., Morena, F., and Andrieu, H.: Development of a distributed hydrological model based on urban databanks – production processes of URBS, Water Sci. Technol., 52, 241–248, https://doi.org/10.2166/wst.2005.0139, 2005.
Savenije, H. H. G.: Determination of evaporation from a catchment water balance at a monthly time scale, Hydrol. Earth Syst. Sci., 1, 93–100, https://doi.org/10.5194/hess1931997, 1997.
Savenije, H. H. G.: HESS Opinions “Topography driven conceptual modelling (FLEXTopo)”, Hydrol. Earth Syst. Sci., 14, 2681–2692, https://doi.org/10.5194/hess1426812010, 2010.
Savenije, H. H. G. and Hrachowitz, M.: HESS Opinions “Catchments as metaorganisms — a new blueprint for hydrological modelling”, Hydrol. Earth Syst. Sci., 21, 1107–1116, https://doi.org/10.5194/hess2111072017, 2017.
Solomatine, D. P. and Shrestha, D. L.: A novel method to estimate model uncertainty using machine learningtechniques, Water Resour. Res., 45, W00B11, https://doi.org/10.1029/2008WR006839, 2009.
Sriwongsitanon, N.: Flood Forecasting System Development for the Upper Ping River Basin, Kasetsart J., Nat. Sci., 44, 717–731, 2010.
Sriwongsitanon, N. and Taesombat, W.: Effects of land cover on runoff coefficient, J. Hydrol., 410, 226–238, https://doi.org/10.1016/j.jhydrol.2011.09.021, 2011.
Sriwongsitanon, N., Gao, H., Savenije, H. H. G., Maekan, E., Saengsawang, S., and Thianpopirug, S.: Comparing the Normalized Difference Infrared Index (NDII) with root zone storage in a lumped conceptual model, Hydrol. Earth Syst. Sci., 20, 3361–3377, https://doi.org/10.5194/hess2033612016, 2016.
Suresh Babu, P. and Mishra, S. K.: Improved SCSCN–inspired model, J. Hydrol. Eng., 17, 1164–1172, https://doi.org/10.1061/(ASCE)HE.19435584.0000435, 2011.
Tekleab, S., Uhlenbrook, S., Savenije, H. H. G., Mohamed, Y., and Wenninger, J.: Modelling rainfall–runoff processes of the Chemoga and Jedeb mesoscale catchments in the Abay/Upper Blue Nile basin, Ethiopia, Hydrolog. Sci. J., 60, 2029–2046, https://doi.org/10.1080/02626667.2015.1032292, 2015.
Tingsanchali, T. and Gautam, M. R.: Application of tank, NAM, ARMA and neural network models to flood forecasting, Hydrol. Process., 14, 2473–2487, https://doi.org/10.1002/10991085(20001015)14:14<2473::AIDHYP109>3.0.CO;2J, 2000.
Vaitiekuniene, J.: Application of rainfallrunoff model to set up the water balance for Lithuanian river basins, Environ. Res. Eng. Manage., 1, 34–44, 2005.
Vermote, E. F., Kotchenova, S. Y., and Ray, J. P.: MODIS Surface Reflectance user's guide, version 1.3, MODIS Land Surface Reflectance Science Computing Facility, https://lpdaac.usgs.gov/documents/445/MOD09_User_Guide_V5.pdf (last access: 22 Febuary 2020), 2011.
Vrugt, J. A., Gupta, H. V., Bastidas, L. A., Bouten, W., and Sorooshian, S.: Effective and efficient algorithm for ultiobjective optimization of hydrologic models, Water Resour. Res., 39, 1214, https://doi.org/10.1029/2002WR001746, 2003.
Xulu, S., Peerbhay, K., Gebreslasie, M., and Ismail, R.: Drought Influence on Forest Plantations in Zululand, South Africa, Using MODIS Time Series and Climate Data, Forests, 9, 528, https://doi.org/10.3390/f9090528, 2018.
Yahya, B. M., Devi, N. M., and Umrikar, B.: Flood hazard mapping by integrated GISSCS model, Int. J. Geomat. Geosci., 1, 489–500, 2010.
Yew Gan, T., Dlamini, E. M., and Biftu, G. F.: Effects of model complexity and structure, data quality, and objective functions on hydrologic modeling, J. Hydrol., 192, 81–103, https://doi.org/10.1016/S00221694(96)031149, 1997.