The SWAT-T model is a further development of the SWAT ecohydrological model . The key enhancement of SWAT-T is the modification of the plant growth module to better represent the plant phenology in tropical regions . A detailed explanation of the advances as well as a comprehensive case study of SWAT-T is presented in . The SWAT-T model leads to more accurate predictions of AET in regions with perennial plant growth . Aside from the plant module modification, SWAT-T remains identical to SWAT. The original SWAT model has been applied for river basins worldwide , as well as regionally in West Africa and for different catchments in Benin . In Benin, most studies focused on the discharge simulation for the Ouémé River basin and as its subcatchments . Remotely sensed AET has also been a primary model optimization target to predict streamflow . To the best of our knowledge, the SWAT-T model has only been applied in Benin for the assessment of the LAI–AET interaction , which demonstrates the AET prediction through LAI on the footprint scale.

The SWAT/SWAT-T model is a process-oriented ecohydrological model for the catchment-scale water balance estimation. It partitions the computation of the water balance into five storages: snow, canopy storage, soil profiles, and a shallow and a deep aquifer. The water balance in SWAT/SWAT-T is described as: 1ΔS=∑i=1NP-Qtotal-E-wlosses, where ΔS is the change in storage; N is the time in days or months; P is the precipitation; Qtotal is the total runoff (in SWAT “water yield”) as a sum of surface runoff, lateral flow, and base flow; E is AET; and wlosses are deep groundwater losses. We use the Soil Conservation Service (SCS) curve number (CN) method to model the surface runoff. The SWAT/SWAT-T model is spatially discretized into subbasins and further subdivided into hydrological response units (HRUs). Generally, three methods are available to compute potential ET (E0) in SWAT/SWAT-T: the Hargreaves method , the Priestley-Taylor method , and the Penman-Monteith method . Once E0 is calculated, it is partitioned into potential plant transpiration (Tplant) and potential soil evaporation (Esoil). We use the Penman-Monteith method since it integrates temperature- and energy-based assumptions and vegetation. demonstrated that this method outperforms the other two options for AET prediction for characteristic West African land cover types (forest and shrubland), which are also most dominant in the Bétérou Catchment. E0 after Penman-Monteith in SWAT/SWAT-T is computed as: 2E0=Δ⋅(Rnet-G)+ρair⋅cp⋅(ez0-ez)/raλ⋅Δ+γ⋅(1+rc/ra), where λ is the latent heat of vaporization; Rnet is the net radiation; Δ is the slope of the saturation vapor pressure-temperature curve; ρair is the air density; cp is the specific heat at constant pressure; ez0 is the saturation vapor pressure of air at height z; ez is the water vapor pressure of air at height z; ra is the aerodynamic resistance; and rc is the plant canopy resistance. In SWAT/SWAT-T with E0 after Penman-Monteith, ra and rc are attributed to an alfalfa grass reference with constant values . For the computation of Tplant, the Penman-Monteith equation is also used, and ra and rc are derived from an actual plant, where the plant's modeled canopy height and LAI are considered. Esoil is then determined as Esoil=E0-Tplant. Actual plant transpiration and soil evaporation are computed based on water availability and biophysical factors, such as LAI, root depth, and soil properties such as field capacity. The total actual evapotranspiration (AET) is the sum of actual plant transpiration and soil evaporation. A detailed description of the plant growth modeling with SWAT-T can be found in the Supplement to this article.

A total of 26 parameters have been investigated in this study (Table ). The list of 26 parameters is gathered with a literature review, where selected SWAT studies in Benin and SWAT-T studies have been screened. The focus was on peer-review studies, where the parameter choice and calibration strategy is transparent. All 26 parameters in Table  are evaluated in the sensitivity analysis. Based on the sensitivity ranking, the 15 most influential parameters are considered for the model optimization. The plant parameters FRGRW₂, LAIMX₂, DLAI, T_BASE, PHU, and GSI are specified in SWAT-T for each land cover type. The Bétérou Catchment is primarily covered by forest, shrubland, and cropland. To avoid an excessive number of parameters (e.g., 5 plant parameters × 3 land cover types), we transfer the parameter values for FRGRW₂, LAIMX₂, DLAI, T_BASE, PHU, and GSI for each land cover type from , who investigated specific sites within the Bétérou Catchment. In the model optimization, each plant parameter is changed through a multiplier from 0.75 to 1.25 (see Table ).

Insensitive parameters are not considered in the optimization but kept constant (Table ). The parameters ALAI_MIN and BLAI are derived based on the maximal and minimal values of the GLASS-LAI time series, respectively. The multipliers for the plant parameters and constant values for insensitive parameters are used in the optimization to reduce the parameter space and address parameter equifinality.

The Bétérou catchment is the Northern most headwater catchment of the Ouémé River Basin in Benin (Fig. ). It covers an area of 10 100 km². The catchment is relatively flat with an altitude range from 248 to 598 m above sea level. The land cover is a typical West African combination mainly consisting of gallery forests (52 %), savannah (33 %), and agriculture (15 %), where the landscape is marked by scattered mosaics of woody savannah with grassy understories. The region is in a natural state with a small share of urban areas (< 1 %) .

The agriculture in this region is combination of cash crops, such as cotton, and subsistence-driven smallholders. The cultivated crops are usually crops like maize, yam, sorghum, bean, millet, or cassava . The cropping calender in this region follows the rainy season and can generally be dated from May to November – with crop specific variations in planting and harvesting . The distinct dry and wet streamflow seasonality is typical for catchments along the Guinean Savanna. The maximum daily average streamflow occurs in September with an estimated value of 290 m³ s⁻¹. The climate in the Bétérou catchment is typical for Sub-Saharan, sub-humid Africa. The annual precipitation ranges from 1100 to 1500 mm . The precipitation pattern is unimodal, with a rainy season between April and October, whereas from November to March, the dry season occurs. The annual mean daily temperature is 25 °C . The soils in the Bétérou catchment consist of ferric soils with loamy sand present in the upper soil horizons . Generally, the AET data follows the seasonal pattern of LAI (Fig. c). In addition, a decrease in AET in the wet season can be observed for this region. During this season, the atmospheric water demand is reduced due high air humidity resulting to low water pressure deficits as shown by .

For the model set up, we subdivided the Bétérou Catchment into 14 subbasins and 126 HRUs. The dominant land cover types forest, shrubland, and cropland are assigned to the SWAT/SWAT-T data base of “FRSD”, “RNGE”, and “AGRL”, respectively. The share of cropland cover is relatively small (15 %) and data for detailed crop modeling is scarce. Thus, we apply a conceptual approach to integrate crops. We assign the generic plant data base “AGRL”. For cropland, the plant growth is not automatically triggered by the soil moisture index but its planting and harvesting is scheduled in the management operations to follow the rainy season . We validate the conceptual approach by comparing the simulated LAI of croplands to the reference of GLASS-LAI. The simulations are run for daily timesteps.

The SWAT-T model is set up using the data listed in Table . We force the models with daily data for precipitation, temperature, relative humidity, solar radiation, and wind speed from two data sets. The first set of forcing data is the ISIMIP3b baseline data set from ISIMIP3b (W5E5 data; ). In ISIMIP3b, the baseline (historical) forcing is not a raw GCM historical simulation but the observation-based W5E5 dataset, to which all future GCM projections are bias-adjusted. This ensures that the statistical properties (mean, variability, and seasonality) of the future projections are consistent with the historical reference climate. For the Bétérou Catchment, the W5E5 baseline captures the main precipitation and temperature patterns reasonably well (see Fig.  in Appendix). Generally, the W5E5-precipitation tends to smooth individual extreme rainfall events and the W5E5-temperature reflects higher maximal daily temperatures than the monitored data.

Secondly, we use observed data (AMMA-CATCH and precipitation data from KASS-D ). The precipitation point observations were regionalized to a 0.05° grid (5 km resolution) using a simple nearest-neighbor interpolation approach . The resulting gridded rainfall fields were then spatially averaged over each sub-catchment to produce 14 sub-catchment-specific rainfall time series, which served as precipitation inputs in the model setup. Both data sets (W5E5 and observed data) are used for a separate calibration.

As a reference to SWAT-T's simulated LAI, we use the GLASS-LAI data, because it is reasonable performance in applications worldwide . Typically, raw LAI data from satellite-based products (e.g., MODIS) encounter uncertainties and noise in tropical regions . The GLASS-LAI algorithm provides robust LAI data through tailoring MODIS and CYCLOPES products with general regression neural networks . In SWAT-T, LAI is simulated for each HRU. The raster cell size of the GLASS-LAI data does not align with the HRUs. We compute the mean GLASS-LAI value for each subbasin by averaging the values of the raster cells that cover it. We then aggregate the simulated LAI at the subbasin scale to match the GLASS-LAI reference. For each subbasin, the subbasin-averaged LAI value is calculated based on the corresponding HRUs, with the LAI weighted by the area of each HRU within the subbasin. In the model evaluation, we compare the weighted-average subbasin values of the GLASS-LAI to simulated LAI. For AET, we use FLUXCOM-AET , which combines MODIS, meteorological and eddy covariance tower data to provide a global data set of energy fluxes.

We have compared the FLUXCOM-AET product and with eddy-covariance (EC) observations from the Bellefoungou flux tower (latitude 9.791° N, longitude 1.718° E) which is part of the AMMA-Catch network in the western part of the Bétérou Catchment in Benin (see Fig. ). This validation shows a reasonable fit of FLUXCOM-AET to the EC data (see Fig.  and Table  in Appendix). The EC data represents a localized area and does not fully capture the spatial heterogeneity within the corresponding FLUXCOM-AET grid cell (0.083°) or the entire Bétérou catchment. Consequently, some discrepancies in both magnitude and short-term dynamics between FLUXCOM and tower-based AET can be expected. Nevertheless, FLUXCOM-AET captures the seasonal timing and variability of monitored AET reasonably well. Therefore, we apply FLUXCOM-AET as a promising validation dataset for catchment-scale AET patterns.

Similarly to LAI, the AET data is aggregated at the subbasin scale, with FLUXCOM raster cells weighted by the area of each subbasin. In SWAT-T, AET is simulated at the HRU level and averaged at the subbasin scale. We compare the weighted-averaged AET data from both FLUXCOM and SWAT-T at the subbasin scale. We use the Kling-Gupta efficiency (EKGE) after , percent bias (EPBIAS), and the coefficient of determination (ER2) (see Appendix for additional equations) to quantify the SWAT-T performance to model daily discharge, LAI, and AET. Observed streamflow, GLASS-LAI, and FLUXCOM-AET data is respectively available from January 1996 to December 2020, January 2001 to December 2015, and January 2001 to December 2014. Observed precipitation data is available from Januaray 1981 to December 2017. The meteorological data from the AMMA-Catch network is available from January 2006 to December 2020.

We use the daily meteorological data of five GCMs and three scenarios from CMIP6 for the climate impact assessment. The selected scenarios refer to the low, medium-to-high, and high narrative SSP scenarios 1–2.6, 3–7.0, and 5–8.5, respectively, where number before the hyphen represents the shared socio-economic pathway narrative and the number after indicates the radiative forcing by 2100 in W m⁻². We use the climate forcing data from the ISIMIP3b initiative with five GCMs (GFDL-ESM4, IPSL-CM6A-LR, MPI-ESM1-2-HR, MRI-ESM2-0, UKESM1-0-LL), since it provides bias-adjusted, downscaled, and standardized climate data. It has been widely used by the hydrological impact modeling community, enabling cross-sectoral comparisons. For instance, report the applicability of ISIMIP3b data for West Africa, and assess uncertainties and impacts on agriculture for the region.

Table  lists 26 parameters that have been used to assess the vegetation, AET, and discharge simulation with SWAT/SWAT-T. We perform a sensitivity analysis to understand the parameter response complexity on the catchment-scale. Based on a sensitivity ranking, the parameter space can be reduced to the most important parameters. We apply the elementary effects (or Morris) method . The application of the Morris method allows the quantification of the parameter sensitivity through statistical means (μ∗-values) of the elementary effects . While the elementary effects are a local sensitivity measure, the computation of multiple elementary effects for varying parameter subsets enables a quasi-global sensitivity analysis . Compared to other sensitivity approaches, it is moreover computationally efficient due to its one-at-a-time sampling strategy . A brief description on the mathematical background of the Morris method can be found in Sect. . We apply Latin hypercube sampling to generate a distributed parameter space. We define r=500 and q=26 parameters with which the total sample size is N=13500. The model sensitivity is evaluated with KGE. We quantify the influence of the SWAT-T parameters on streamflow, LAI, and AET separately, denoted as EKGE,Q, EKGE,LAI, and EKGE,AET, respectively. Additionally, we assess the parameter sensitivity on the equally-weighted combination of all three objectives, EKGE,all: 3EKGE,all=13⋅EKGE,Q+13⋅EKGE,LAI+13⋅EKGE,AET. The combination of separated and composed sensitivity analysis objective enhances the understanding of the parameter response complexity. For the application of the Morris method, we integrated the approaches of , , and into a set of MATLAB scripts.

Table  summarizes the calibration strategies with their level of complexity. The approaches consider (i) only streamflow at the catchment outlet (“Q only”); (ii) streamflow at the outlet and LAI at the subcatchment scale (“Q + LAI”); and (iii) streamflow at the outlet, along with LAI and AET at the subcatchment scale (“Q + LAI + AET”). For each calibration strategy, ROPE is applied and a separate objective function EKGE,eff is defined with an equally weighted contribution EKGE,v of the considered variable v (Table ). We integrate the modeling of LAI as proxy of hydrological process representation spatially distributed in the entire catchment. AET data for model performance evaluation is often limited or the necessary data to compute AET, e.g., with Penman-Monteith, is lacking. Given the LAI–AET linkage in this region , we test whether a calibration approach with streamflow and LAI can be used to predict AET.

The robust parameter estimation (ROPE) algorithm after is used for the optimization of the three modeling approaches. ROPE is a depth function-based algorithm to iteratively find parameter sets that give good model predictions and are robust. The ROPE algorithm does not necessarily give the best parameter sets (global optimum), but rather an ensemble of different near-optimum parameters sets . We apply ROPE to derive an ensemble of good parameter sets for climate impact assessment. In addition, the limitations and dependency on a single parameter set (equifinality) are overcome by incorporating multiple well-performing sets. We define “well-performing” if the parameter set gives satisfactory KGE values (EKGE≥0.6) according to .

With ROPE, parameter robustness is quantified using the halfspace depth function introduced by . Generally, data depth is a quantitative approach to define how central a sampled point is in a multivariate distribution. The more central and, thus, deep a point is located, the less the corresponding objective value differs if the point is subject to a small change, e.g., parameter changes . The calculated parameter set depth value is denoted as L in the following.

To apply the ROPE algorithm, random samples XN of N parameters and their bounds d are defined. The model is run for each sample and its performance is quantified. A subset XN∗ of the best performing parameters is defined, e.g., the best 10 % of XN. The performance is evaluated with an objective function (here EKGE,eff, Table ). Based on XN∗, M new random samples (labeled YM) are generated according to the corresponding depth values L. A new parameter vector of YM is only considered if it lays within the convex hull given by the depth values. The boundary of the convex hull is defined as L≥1. Each new parameter vector inside the hull (L≥1) is considered for the sampling of the next iteration. The iterations of the ROPE algorithm stop as soon as a convergence criteria is reached. We figured that after 12 (20) iterations for the “Q only” and “Q + LAI” (“Q + LAI + AET”) the model performance plateaus. From one iteration to the next, the 10 % best parameter values are used. After the final iteration, we consider the 20 best parameter sets for the climate impact assessment.

We calibrate against streamflow, GLASS-LAI, and FLUXCOM-AET. The calibration (validation) is defined from January 2008 to December 2011 (January 2012 to December 2015). We use eight years to enable a computational cost-efficient calibration for the total of six applied model optimizations. The period from 2006 to December 2015 is an adequate overlap of W5E5 and observed data to enable a fair and comprehensive comparison of the forcing data.

The model performance to predict streamflow, LAI, and AET is separately evaluated with EKGE and for each forcing data set. For model validation, we evaluate the efficiencies EKGE, percent bias EPBIAS, and the coefficient of determination ER2 of streamflow, AET, and LAI (see Appendix for supplementary equations).

We assess the reliability of the derived parameter sets through a benchmark evaluation. For hydrological models, benchmarks facilitate an improved interpretation of the model performance and KGE metrics. . In this study, we define a lower level of expected model performance based on the interannual mean for every calendar day (“seasonal benchmark”) as suggested by . For this purpose, we compute the daily mean seasonal values for discharge from the observed data of the period 2001 to 2015. For LAI and AET, the mean seasonal values are based on the 8-daily temporal resolution for according to the GLASS-LAI and FLUXCOM-AET data source. Each seasonality is duplicated to be as long as the evaluation period (2008 to 2015), and, thus, represents the identical temporal resolution of the model simulation. We compare the duplicated seasonality to the observed time series and compute EKGE,Q, EKGE,LAI, and EKGE,AET for the entire evaluation period. We use the duplicated seasonality (“seasonal benchmark”) as a synthetic predictor to benchmark SWAT-T's modeling performance.

The applicability of the parameter sets is evaluated with a cross validation. The parameter sets derived with observed data input are run with model forcing of ISIMIP3b reference data W5E5; and vice-versa the W5E5-parameters with observed data forcing. The cross validation is conducted for the evaluation period (January 2008 to December 2015) and assessed through EKGE,Q, EKGE,LAI, and EKGE,AET. We focus the cross validation on the most comprehensive calibration approach “Q + LAI + AET”.

The ISIMIP3b project provides meteorological forcing data from five climate models for three future scenarios, and the corresponding baseline, or historical data. The baseline data set represents the observed climate characteristics, usually for the period from 1901 to 2015, and is used as a reference point to evaluate the future climate scenarios. For meteorological changes, we calculate the mean annual precipitation for each baseline data of the climate models. For each climate model and scenario, we compute the annual precipitation for each year of the near (2031 to 2050) and far (2070 to 2099) future periods. For each period, the mean annual precipitation is calculated. The resulting mean annual precipitation per period is compared to the mean annual precipitation of the baseline data, and the increase or decrease (compared to the baseline data) is determined. For the temperature change, we compute the mean annual temperature for each climate model/scenario and its baseline data. The multi-model mean is calculated from all five GCMs. For hydrological changes, we evaluate the target variables AET and discharge. As baseline for AET (discharge), we use the simulated mean annual AET (discharge) from the period January 2001 to December 2015 from the W5E5 forcing.

The sensitivity of the parameters is quantified with the elementary effects method after whereby a ranking is possible. The higher the μ∗ value of a parameter is, the higher is the parameter sensitivity . Basically, the elementary effects are calculated according to the model responses, represented as objective function KGE, to the parameter changes. Four different objectives EKGE,Q, EKGE,LAI, EKGE,AET, and EKGE,all are defined to evaluate the parameter sensitivity when modeling discharge, LAI, AET, and the combined performance, respectively, with SWAT-T. Figure  shows the values μ∗. The sensitivity analysis considers the entire evaluation period from 2001 to 2015.

The parameter sensitivity for the equally-weighted objective EKGE,all is displayed in Fig. a. The ranking of the first 8 parameters is similar to the discharge-only analysis in Fig. c indicating a strong influence of the discharge performance EKGE,Q in EKGE,all. The ranking 9 to 15 is assigned to relevant plant, soil, and AET parameters that are sensitive for EKGE,LAI and EKGE,AET.

Figure b summarizes the sensitivity ranking with respect to LAI. The LAI parameters (DLAI, T_BASE, FRGRW2, PHU, LAIMX2, LAIMX1; green circles) are generally the most sensitive, together with the SOL_AWC (soil; in purple) and GSI (AET; in red) parameters. The sensitivity ranking confirms the linkage of LAI and AET in the region. The plant growth depends on water and temperature stresses in SWAT/SWAT-T. The water stress is a function of AET and water availability, in which the parameters SOL_AWC and GSI are pivotal. The other soil and AET parameters are less influential, while groundwater and runoff parameters are irrelevant for predicting LAI.

The sensitivity for the discharge modeling is summarized in Fig. c. It can be observed that the most sensitive parameters are the groundwater (GWQMN, RCRG_DP, GW_REVAP, REVAPMN) and soil parameters (SOL_AWC, SOL_BD), where GWQMN is most influential. It governs the water depth for which baseflow can occur and partitions the flow into surface and baseflow contributions. Large GWQMN values allow more surface runoff and less baseflow – and vice versa for small values. The streamflow at the Bétérou gauge is seasonal and connected to the dry and wet precipitation periods. The groundwater parameters, and GWQMN in particular, are essential to model the low flows in the dry season. Figure c shows that the groundwater parameters have more impact on the discharge modeling than the runoff (CN_2, CH_K2, HRU_SLP, SURLAG, CH_N2) or plant parameters. Similar findings of the strong groundwater relevance in tropical regions are also reported in the SWAT-T application of .

The sensitivity ranking in Fig. d highlights the soil (SOL_AWC, SOL_BD) and AET (ESCO, GSI, EPCO, CAN_MX) parameters to be most sensitive for AET modeling. The AET parameters directly govern the calculation of AET. The soil parameters determine the soil water and availability and, thus, indirectly contribute to the actual amount of AET. The plant parameters are comparably less influential, yet are relevant. The contribution of the plant parameters affirm the significance of LAI for AET in this region. Groundwater and runoff parameters are not substantial for AET prediction.

The sensitivity analysis is applied to identify relevant parameters and reduce the parameter space for a more efficient optimization. Overall, 15 parameters are selected that are substantial for the prediction of discharge, LAI, and AET. The final parameters are listed in Table .

The models are forced with observed and W5E5 data. The derived parameter sets from the W5E5 forcing are used for the climate impact assessment. The purposes of the optimization with observed input data are to demonstrate the general performance potential of the ROPE algorithm and to enable a parameter cross validation for the models with W5E5 parameters. We select the 20 best runs from the ROPE optimization for the subsequent model applications. In the following analysis, Fig.  shows the comparison of simulated and observed discharge, LAI, and AET. Table  summarizes the quantitative model performance as KGE values for the calibration and validation periods.

Overall, the mean annual patterns of discharge, LAI, and AET are moderately well represented in all optimization approaches and with both forcing datasets (Fig. a to c). For discharge, the performance of the daily simulation is tied to the forcing data, as spotlighted by the scatter plots in Fig. d and g. The discharge simulation with observed data align well with the 1 : 1 line in Fig. d, i.e., it reproduces the recorded streamflow accurately. In contrary, the discharge simulation with W5E5 shows a spread in the scatter, where single high flow events are recorded, but not modeled; and vice-versa for single simulated high flow events. The performance quantification according to KGE is given in Table . Model optimization with observed data (EKGE,Q = 0.49 to 0.92, Table , calibration and validation) outperforms the W5E5 model optimization (EKGE,Q = 0.10 to 0.64, calibration and validation). The variation in performance can be explained with differences in the input precipitation data sets. For instance, the precipitation intensity in 2010 is too low in W5E5 compared to the observation data (see Fig. ). In this period, the simulated discharge is underestimated with W5E5.

Figure b shows the mean annual LAI pattern from all 20 parameter sets for each optimization approach. Simulated LAI is compared to GLASS-LAI (8-daily) in Fig. e and h for model forcing with observed and W5E5 data, respectively. The performance to model LAI is particularly good in the modeling approaches with LAI as target variable (“Q + LAI” and “Q + LAI + AET”) for both forcing data sets. In the growing phase, a small underestimation compared to GLASS-LAI is generally present. However, the plant maturity and leaf senescence are modeled well with “Q + LAI” and “Q + LAI + AET”. The performance metrics EKGE,LAI ranges from 0.75 (validation) to 0.88 (calibration) for these approaches (Table ). For “Q only”, the LAI modeling is less adequate, particularly for the W5E5 forcing, where the models fail to simulate plant maturity and leaf senescence. The LAI calibration varies from EKGE,LAI = 0.25 to EKGE,LAI = 0.86 (Table ) for discharge-only approaches (“Q-o” and “Q-w”).

Similar to LAI, Fig. c shows the modeled mean annual AET pattern as an average of the 20 parameter sets. Overall, the simulated mean AET closely reproduces the annual pattern of FLUXCOM reference. The SWAT-T model can represent the two “AET peaks” in the beginning and end of the rainy season, as well as the “AET dip” in the middle of the rainy season. However, the FLUXCOM-AET values are continuously underestimated. It can be observed that “Q only” particularly underestimates the reference AET and even fails to mimic the second “AET peak” due to the lack of LAI representation in this period. The scatter plots in Fig. f and i confirm SWAT-T's capability to predict AET, but also the simulated underestimation to FLUXCOM-AET. Generally, the optimization approaches with AET (“Q + LAI + AET”) result in very good performance values for both forcing data sets (EKGE,AET = 0.75 to EKGE,AET = 0.90, Table ). If the model is calibrated with respect to discharge and LAI, adequate AET predictions are possible. For the “Q + LAI”, EKGE,AET varies between 0.71 to 0.90 (Table ).

The investigation of multiple parameter sets confirms that the modeling performance of SWAT-T is tied to the forcing data, particularly for discharge. Figure  shows the empirical cumulative distribution functions for the 20 best runs from the ROPE optimization for the calibration period. In the figure, the performance of the seasonal benchmark of discharge, LAI, and AET is integrated.

The performance distribution of the parameter sets to model discharge is evaluated in Fig. a. Two distinct clusters can be identified: the distribution of the W5E5-based parameters (in blue, with EKGE,Q from 0.54 to 0.64) and based on the observed data set (in red, with EKGE,Q from 0.81 to 0.91). The W5E5 modeling results in a performance as good as the seasonal discharge benchmark (pink line, KGE = 0.63). Notably, the model outperforms the seasonal discharge benchmark, if the observed data set is applied. Figure b shows the distribution with respect to EKGE,LAI. The values of EKGE,LAI from discharge-only (Q-o, Q-w) decline rapidly, while the other approaches yield equally good parameter sets for LAI. Figure c displays the distribution functions of EKGE,AET. All parameter sets for the “Q + LAI” and “Q + LAI + AET” approaches result in values of EKGE,AET ≥ 0.68. In contrary, only a small subset of the “Q-only” approach yields EKGE,AET ≥ 0.68. Figure  highlights that the ROPE algorithm is applicable to derive multiple equally-good parameter sets for LAI and AET modeling. The analysis in Fig.  further indicates that the integration of LAI alone does not necessarily improve the streamflow modeling with SWAT-T. As shown in Fig. , the influence of LAI parameters on streamflow prediction is less sensitive. Yet, the SWAT-T model can't beat the benchmark for LAI or AET. GLASS-LAI and FLUXCOM-AET are reanalysis outputs based on a data assimilation of observations and satellite-based data, such as MODIS data. Both data sets inherit observed seasonal patterns and are processed to mimic reality . Their benchmark is thus hard to outperform for a physical-describing model like SWAT-T.

We assess the applicability of the W5E5-parameters with a parameter cross validation. Figure  shows the evaluation of the cross validation for “Q + LAI + AET”. If the model is run with W5E5 parameters and observed data forcing (red dashed lines in Fig. ), the resulting metrics EKGE,Q and EKGE,AET are mainly ≥ 0.6, which aligns with a satisfactory to good modeling performance according to . The median EKGE,Q for the calibration and validation period is 0.70 and 0.69, respectively, which also outperforms the seasonal discharge benchmark (KGE = 0.63). For the prediction of LAI, EKGE,LAI varies for the W5E5-parameters (red dashed lines in Fig. b, e), where the median of EKGE,LAI is 0.48 (calibration) and 0.42 (validation). The variation for W5E5 parameters and observed forcing can be explained through the final parameter values (see Table ). For W5E5 optimization, the T_BASE and PHU parameter values are smaller than for observed data forcing. Lower T_BASE and PHU values trigger an earlier start of the plant growth along with a shorter period of maturity. If the input data is replaced, even small temperature changes can influence the growth start and plant maturity. Given the dependency of LAI on temperature in SWAT-T, its robustness for climate impact assessment is crucial. For climate impact assessment, multiple LAI parameter sets can be used to avoid the reliance on single LAI parameter sets and increase the LAI parameter robustness.

Figure  compares the future changes of precipitation and temperature of the climate models to their baseline data. For all climate models and scenarios, the average temperature is projected to increase between 0.5 to 5.8 °C in the future. All climate models indicate a temperature increase of more than 3 °C for the higher emission scenarios (SSP3–7.0 and SSP5–8.5) in the far future (filled orange and purple symbols), with highest change values of 5.8 and 5.4 °C modeled by UKESM1-0-LL and IPSL-CM6A-LR, respectively. The temperature increase for the low emission scenario (SSP1–2.6) is comparably smaller (0.9 to 2.0 °C) in the near and far future. The multi-model mean temperature change for this region ranges from 1.2 °C (SSP1–2.6, far future) to 4.2 °C (SSP5–8.5, far future). The mean temperature is expected to get warmer especially in the dry period of the year (January to June and October to December), where a prolongation of the dry and hot period (historically from February to May) to January to June is projected (Fig. ).

The precipitation projection is mixed with increases and decreases in the future, where the multi-model precipitation mean varies between increases of 1.2 % (SSP1–2.6) and 7.2 % (SSP5–8.5). For the SSP1–2.6 scenarios, the precipitation changes are ±10 % in the near and far future, except for the near future in the MRI-ESM2-0 model. Positive precipitation changes for selected scenarios are higher than 15 % in UKESM1-0-LL, IPSL-CM6A-LR, MPI-ESM1-2-HR, and MRI-ESM2-0 GCMs. The GFDL-ESM4 model (circles) mainly projects precipitation decreases. The change in precipitation is most substantial in the rainy season. From all GCMs, an increase in the amount of monthly precipitation is projected for August, September, and October (see Fig. ). The change values are summarized in Table . The future changes of solar radiation, relative humidity, and wind speed for the near and far future are displayed in Figs.  to .

The focus of this comparative study is the analysis of the implications of different modeling approaches on the prediction of future hydrological changes, particularly of AET and discharge. We evaluate the future discharge and AET response to climatic changes for the near (2031–2050) and far (2070–2099) future period. Generally, the climate change signal for AET is similar across all calibration strategies, where the multi-model mean values of AET indicate future increases. Yet, the future annual AET can vary dependent on the calibration strategy, the underlying GCM, the scenario, and the future period. Figure  shows the annual AET for the applied GCM-SSP sets and the multi-model mean of all GCMs for the near and far future period. The analysis proves that the discharge only strategy (“Q only”, blue boxcharts) predicts lower annual AET rates compared to the integrated optimization approaches for each applied GCM-SSP set and for the multi-model mean values for annual AET. The simulated future annual AET values are lower for “Q only” (846 to 889 mm yr⁻¹, depending on the GCM-SSP scenario), while the climate sensitivity of “Q + LAI” (886 to 918 mm yr⁻¹) and “Q + LAI + AET” (904 to 932 mm yr⁻¹) is similar (Table ). The lower rates of AET for “Q only” is furthermore evident if the future AET values are compared to the respective AET baseline. The change rates of annual AET from “Q only” to “Q + LAI + AET” are close for the near future period. While “Q only” predicts an increase of annual AET of 5.5 %, 5.7 %, and 4.2 %, “Q + LAI + AET” results in increases of 7.4 %, 7.2 %, and 5.8 % for the SSP1–2.6, SSP 3–7.0, and SSP5–8.5 scenarios (near period, multi-model mean values, Table ), respectively. For the far future period, the change rates are generally higher for “Q + LAI + AET” than “Q only”. However, for the SSP5–8.5 scenario, the future change of AET for “Q only” is 9.2 % compared to an increase of 8.4 % for “Q + LAI + AET”. Figure  moreover indicates that “Q + LAI” results in predictions of annual AET which are close to “Q + LAI + AET”. Thus, “Q + LAI” can approximate AET rates on the catchment scale even if AET is not integrated in the parameter estimation. For the given setup of SWAT-T modeling in the Bétérou Catchment, the LAI computation can be potentially used as a proxy for AET on the catchment scale.

The future annual AET rates also depend on the GCM model. For the near future (Fig. ), the range of predicted annual AET is smaller (750 to 1000 mm) than compared to the far future GCM model applications (710 to 1120 mm, Fig. ). The modeling with GFDL-ESM4 and MPI-ESM-1-2-HR result in comparably smaller annual AET rates, while the simulations forced by IPSL-CM6A-LR, MRI-ESM2-0, and UKESM1-0-LL indicate higher annual AET values for the near future. For the far future period, high annual AET rates are particularly predicted with IPSL-CM6A-LR. The GCM forecasts for the far future an increase in precipitation and temperature (Fig. ), which translates to available water and energy for AET. The simulation with forcing data from the MPI-ESM2-0 model shows distinct differences between the SSP scenarios, where the also future decreases of annual AET rates can occur. For the far future, the MPI-ESM2-0 model simulates the low temperature changes (ΔT ≤ 2 °C, Fig. ), and, thus, the increase in AET is the lowest with MPI-ESM2-0 compared to the other GCM applications.

The differences in predicting annual AET are explained with variations in monthly AET dynamics. Figures  and show the monthly mean AET for the calibration approaches for the near and far future period, respectively. The figures also display the monthly mean AET values from the each model optimization forced with W5E5 data (black cross). The analysis shows for which months the calibration approaches predict different AET rates. For the near future period (Fig. ), the monthly AET rates predicted with “Q only” are lower in the period from August to December. This period describes the end of the wet and transition into the dry period. Given the discharge-centered optimization, the model does not represent the plant and soil water processes in this period. The smaller AET rates for “Q only” are evident for all three SSP scenarios in the near future period. Concurrently, the predicted mean seasonal values from “Q + LAI” are close to the “Q + LAI + AET” where only small differences per month are computed.

For the far future period (Fig. ), the pattern of monthly mean AET differences is similar to the predictions of the near future period. The monthly AET rates of “Q only” are underestimated from August to December. The analysis for the far future period moreover shows large variations within each calibration strategy itself. The box charts are wider in the far period, than in the near period evaluations, particularly for SSP5-8.5 (see Fig. c). The chart sizes indicate that within each calibration strategy, the range of monthly AET varies, and, thus, the forcing from the GCMs is varying, too. The comparison of future monthly mean AET values to the monthly AET rates from the W5E5 optimization points out differences particularly in April for all calibration approaches and SSP scenarios. The contrasting high rates of future monthly AET in April is related to the strong increase of temperature in the future spring months (see Fig. ). In April, the warmer climate overlaps with the onset of the rainy season, which eventually catalyzes the AET generation.

Generally, the future rate of AET is computed to increase for the Bétérou Catchment based on the GCMs integrated in the ISIMIP3 project. Depending on the calibration strategy and the SSP scenario, the multi-model mean changes range from 4.0 % (“Q only”, SSP1–2.6, far period) to 9.2 % (“Q + LAI + AET”, SSP3–7.0, far period; and “Q only”, SSP5–8.5, far period) as indicated from the multi-model mean values in Fig.  and Table . The projected changes of AET (multi-model mean) are lower for “Q only” than “Q + LAI + AET”. The overview of future changes in Fig.  and Table  indicates high AET increases particularly for the IPSL-CM6A-LR, which projects high precipitation and temperature increase. Negative AET changes are computed only partially and mainly with forcing from MPI-ESM-1-2-HR (Fig. a).

For discharge, the future changes can be positive and negative depending on the underlying GCM forcing and calibration strategy. The application of the MPI-ESM1-2-HR model data results in a strong increase of discharge (up to 101.1 %, SSP5–8.5, “Q only”). The forcing data of GFDL-ESM4 yields mostly decreases of future discharge due to projected precipitation decreases and temperature increases of > 3 °C. For the near future, the discharge changes from “Q only” range from -3.3 % to 3.2 % (multi-model mean, Table ). For “Q + LAI” and “Q + LAI + AET”, the multi-model mean changes of discharges are simulated from -7.9 % to -4.0 % and from -10.8 % to -6.7 %, respectively. For the far-future period, the multi-model mean changes of discharge show a substantial decrease independent of the calibration strategy, where the change rates for “Q + LAI + AET” are the lowest due to the higher simulated evaporative demand. The changes are -18.9 % to -1.6 % for “Q only”, -21.7 % to -5.1 % for “Q + LAI”, and -22.9 % to -7.0 % for “Q + LAI + AET”.

This study primarily focus on the variability of future AET depending on the calibration strategy in hydrological modeling. For this purpose, three calibration methods with different levels of complexity and information indices (“Q only”, “Q + LAI”, and “Q + LAI + AET”) are evaluated. Similar to previous comparative studies , we investigate simple to comprehensive calibration strategies, yet with a particular focus on AET. The different calibration complexities from “Q only” to “Q + LAI + AET” also serve as a proxy for a benchmark test with lower and upper arguments, as suggested by . The integration of an “AET only” approach may add another layer of “upper benchmark”. Yet, we also assess the calibration strategy impacts on discharge changes for the region. Disregarding discharge in the hydrological modeling can imply large discharge variability . We have thus avoided to include “AET only” to minimize discharge uncertainties. The investigation of a Q + AET calibration could also serve as an upper benchmark and is expected to reproduce AET rates well. It is however not considered in the present study, because it would potentially neglect vegetation dynamics to satisfy the atmospheric water demands. Isolating Q + AET from LAI adds a further layer in understanding the role of vegetation processes in AET-focused calibrations. Looking ahead, future work could explicitly explore Q + AET as a further upper limit together with Q + LAI + AET.

Model evaluation in this study is based on in-situ monitoring and satellite-derived reanalysis datasets (GLASS-LAI and FLUXCOM-AET). GLASS-LAI and FLUXCOM-AET have been validated by for two distinct locations within the Bétérou Catchment. GLASS-LAI and FLUXCOM-AET rely on models and assumptions, and are provided at a coarser resolution than the hydrological model. To address scale differences, the SWAT-T outputs, GLASS-LAI and FLUXCOM-AET data were spatially projected onto the same units. While this approach allows for a spatially explicit evaluation that reflects land cover variability within the Bétérou catchment, it also introduces some information loss, e.g., from the aggregation of AET heterogeneity. This is particularly important in the process representation of the LAI-AET interaction for localized land-atmosphere feedbacks. The subbasin-level evaluation offers a practical compromise, where spatial context is preserved at a scale appropriate for hydrological modeling. Smaller scales (e.g., HRU levels) and dynamic land use changes have not been considered given the spatial mismatch of modeled and simulated LAI and AET. As demonstrated in , dynamic land use modeling can further enhance the process representation of surface runoff and AET in the region. Future research can consider a coupling of detailed LAI, AET, and dynamic land use modeling to further enhance robust hydrological modeling for the region.

The application of satellite-based GLASS-LAI and FLUXCOM-AET data for a tropical catchment like the Bétérou Catchment carries uncertainties. Generally, GLASS-LAI and FLUXCOM-AET are both widely used and are considered high quality. Yet, satellite-based datasets in tropical regions can be subject to cloud contamination and reflectance noise or the lack of observation networks for validation . For the Bétérou region, the GLASS-LAI dataset shows promising applicability when compared to monitored data . FLUXCOM-AET has been extensively validated against eddy covariance (EC) measurements across a wide range of climatic conditions, supporting its general reliability . This study relies on a pointwise validation of FLUXCOM-AET and GLASS-LAI. While this approach does not fully capture spatial heterogeneity or all sources of uncertainty, it provides a consistent and pragmatic basis for model evaluation in data-scarce tropical regions. This approach ensures that the validation framework remains robust and fit for purpose.

The hydrological cycle in the Western African region is driven by precipitation, AET, and discharge . Although soil moisture is a key internal state variable in hydrological modeling, we have not included soil moisture products for the model validation. Accurate simulations of the dynamics of AET and discharge has been shown to reliably represent the hydrological behavior in the Ouémé catchment . We thus focused the model evaluation on discharge, AET, and LAI, and excluded soil moisture to minimize additional uncertainty and potential conflicts between observational variable constraints.

We integrate LAI as a proxy for vegetation in the hydrological modeling. The relevance of LAI on the AET and discharge estimation in sub-humid regions has been highlighted prior, e.g., in , , , and . Previous studies for the region, like or , include AET in the hydrological modeling of the Ouémé River in Benin, while others have focused on “Q only” approaches . We demonstrate that a detailed LAI modeling can predict adequate AET outputs.

The evaluation of modeled discharge from the ROPE application yields different KGE ranges depending on the forcing data. For calibration with data from the observation network, the maximum EKGE,Q (validation) is 0.91 (0.92), 0.85 (0.88), and 0.90 (0.90) for “Q only”, “Q + LAI”, and “Q + LAI + AET”, respectively. These results align with findings from other studies; for instance, report 0.83 (0.77) with HBV model, or validate the Parflow-CLM model with a mean KGE value of 0.82. The strong KGE performance confirms that the ROPE algorithm and the SWAT-T are well-suited for discharge calibration in the Bétérou Catchment.

However, it is still essential to calibrate the hydrological impact model with the W5E5 forcing dataset if the bias-adjusted ISIMIP3b GCM outputs are used. In this study, the maximum EKGE,Q from calibration (validation) are 0.64 (0.50), 0.64 (0.49), and 0.63 (0.49) for “Q only”, “Q + LAI”, and “Q + LAI + AET”, respectively. We assume the KGE values are smaller compared to the observed data due to the W5E5 data forcing, particularly precipitation. The W5E5 data set represents the seasonality accurately well, but it lacks the depiction on sub-seasonal variations and single rainfall events for the Bétérou Catchment in West Africa. Similar results with CMIP6 are reported by for the Congo River in Africa who figured challenges in discharge representation through CMIP6 climate data. However, we benchmark the discharge performance with the seasonal benchmark after , which shows an adequate behavioral modeling of discharge. Combined with the parameter cross validation, the W5E5-models are applicable and can still represent the mean seasonal patterns of the hydrological cycle in Bétérou. We thus recommend benchmarking and cross validation with observed data if the ISIMIP3b is used for future studies in West Africa.

The benchmarking moreover contains a Julian Day method for FLUXCOM-AET and GLASS-LAI, where a relatively strong performance for both datasets can be observed. The remotely sensed FLUXCOM-AET and GLASS-LAI products are subject to temporal smoothing and gap-filling, which can regularizes the signal and enhances the skill of a climatological benchmark. Therefore, the benchmark performance should be cautiously interpreted as it is not an evidence of superior model process representation. It is possible that satellite-based datasets do not fully capture the dynamic year-to-year variability in vegetation and evapotranspiration, particularly during cloudy periods or in highly dynamic land surface conditions . The retrieval algorithms and data assimilation for both FLUXCOM-AET and GLASS-LA rely on temporal compositing and interpolation. Thus, the processing influences and smooths the real interannual fluctuations and make the reference signal appear more repeatable from year to year. It is likely that this regularization of the observational data contributes to the high skill of the Julian Day benchmark for FLUXCOM-AET and GLASS-LAI.

The model optimization and selection of the 20 parameter sets for the climate impact assessment are based on the performance metrics KGE, which is a robust criteria to asses the runoff generated within a catchment . Single metrics can however overlook pitfalls of the model structure and defer the model and parameter choice . We apply multiple parameter sets, the benchmarking after , and parameter cross validation to increase the reliability of the single-metrics based model choice for climate impact assessments. Looking ahead, the model optimization can include further performance criteria to assess specific features of hydrological process in the region, such as the Nash-Sutcliff Efficiency (NSE) or the logarithmic NSE, among many others, to better include high or low flows, respectively.

The evaluation of LAI as proxy for AET is site- and model-specific. Although the Bétérou Catchment is representative of parts of West Africa, it captures only a subset of the region's diverse hydrological, phenological, and climatic conditions. The ET dynamics in the region can be generally considered water-limited, since the strong seasonal water deficit governs ecosystem functioning . From a modeling perspective, LAI can be expected to serve as a proxy for AET primarily in physical-oriented hydrological models that use LAI as a vegetation attribute in the framework of ET computation, such as the Penman-Monteith method. In contrast, the relation of LAI and AET is weaker in hydrological models that estimate ET using more simplified or empirical approaches, such as the Hargreaves method. In such models, the suitability of LAI as a proxy for AET is limited. Consequently, the conclusions regarding the suitability of LAI as a proxy for AET may be specific to water-limited environments and should be evaluated in regions with different hydroclimatic regimes as well as the hydrological model applied.

We used five GCMs with three SSP scenarios from ISIMIP3b (based on CMIP6) equally weighted, where we followed the conventional model democracy approach. Model democracy in climate impact studies can be crucial . However, the use of multiple GCM-scenario combinations helps to reduce structural GCM model uncertainty, while enhancing transparency and comparability. For example, the application of CMIP6 data in this study for the Bétérou Catchment in Benin enables comparisons with similar assessments in this region, such as the findings of for the Tougou Basin in Burkina Faso. Both studies indicate substantial increases in projected discharge under the MPI-ESM1-2-HR model, with changes up to 101.1 % in this study and up to 117 % in .

However, assuming all GCM models perform reasonably well globally is critical, since poorly performing GCMs can influence ensemble-based hydrological impact statements . highlight that GCMs are major source of uncertainty in the climate-hydrological impact modeling. Consequently, the influence of the calibration strategy on future AET is likely to be smaller compared to the GCM model uncertainty. Nevertheless, our study demonstrates that the future changes of AET strongly depend on the calibration strategy. Thus, disregarding AET during model calibration may introduce an additional layer of uncertainty other than the GCM uncertainty.

The application of multiple ISIMIP3b GCMs offers a transparent and standard way to understand climate changes and hydrological responses across studies, and to avoid uncertainties from individual downscaling and bias corrections. demonstrate the potential benefits of the ISIMIP3b climate data since its multi-model mean ensembles can even outperform individual models when compared to observed precipitation in Western Africa. The seasonal patterns of temperature and precipitation in the GCM historical data closely align with those observed in the monitoring network. The monthly values are similar, and qualitatively meet the criteria suggested by . The meteorological changes for the Bétérou Catchment in Benin indicate a clear projection for temperature, but a divers picture for precipitation. For temperature, all GCMs forecast increases of 0.5 to 5.8 °C, which generally aligns with the changes projected with CMIP6 for the West African Region .

The projected precipitation derived from the ISIMIP3b data set indicates a diversified scatter of changes. The multi-model mean values indicate small (1.2 %, SSP1–2.6) to large (7.2 %, SSP5–8.5) precipitation increases. Yet, the variability between single models is distinct. Substantially wetter periods are indicated by MPI-ESM1-2-HR, IPSL-CM6A-LR, and UKESM1-0-LL, while GFDL-ESM4 and MRI-ESM2-0 project drier future rainfall. The ensemble spread in CMIP6 precipitation projections can be attributed to differences in how models simulate the West African Monsoon . Similar precipitation variability for West Africa has also been reported in other climate forcing datasets. For example, found substantial spreads in CORDEX precipitation projections that are also linked to the complex monsoon dynamics.

The annual AET is generally predicted to increase independent of the GCM, the time period (near-, far future), or the calibration strategy (“Q only”, “Q + LAI”, “Q + LAI + AET”). For the low emission scenario (SSP1–2.6), the multi-model mean changes range from 4.0 % to 7.4 % across all calibration strategies. The high emission scenario (SSP5–8.5) predicts multi-model mean changes of 4.2 % to 9.2 %. Previous impact studies for the Ouémé region report increases, but also decreases in future AET for the region . However, these studies relied on models calibrated against discharge, without explicitly including AET in the calibration. As our results demonstrate, disregarding AET from calibration can defer the simulated amount of future AET, where “Q only” simulates lower future annual AET changes than “Q + LAI” or “Q + LAI + AET”. These findings align with , who demonstrate the influence of the AET calibration strategy on the climate impact assessment for the Upper Godavari in India.

Differences in projected AET across the calibration approaches are most pronounced between August and December. The calibration strategy has a smaller impact on AET simulation during the early part of the year, which is the shift from dry to wet period. The seasonal AET pattern in Bétérou is characterized by two peaks (April to June and September to November) with a notable dip during the peak rainy season (June to September). This mid-season decline is primarily caused by cloud cover and rainfall, which limit solar radiation and thus reduce AET. SWAT-T can simulate AET dynamics in the beginning of the season (January to July), which is largely independent of the calibration strategy. However, from August onward, the AET simulation becomes more sensitive to the calibration strategy, demonstrating the need for carefully adjusted model parameters. The soil (e.g., SOL_AWC, SOL_BD), AET (e.g., GSI), and plant growth (DLAI, T_BASE, FRGRW2) parameters can be essential to realistically represent post-monsoon and dry-season AET in Bétérou.

The analysis of future discharge changes shows increases and decreases depending on the GCM, emission scenario, and evaluation period. The different projections are largely attributed to the complexity of how GCMs simulate precipitation dynamics and projections of the West African Monsoon . Consequently, the model forcing with MPI-ESM-1-2-HR outputs can also result in high future discharges (change of 101.1 % for “Q only”, SSP5–8.5, far future). Similar substantial discharge changes with MPI-ESM-1-2-HR are also reported by for a catchment in Burkina Faso. This GCM projects substantial precipitation increases with minor temperature changes for the West African region. Overall, the analysis of the multi-model mean for future discharge changes essentially indicates a decrease of discharge for the Bétérou region. The multi-model mean changes vary from 22.9 % to 3.2 %. Discharge decreases for the Ouémé region have also been found in other studies . The decreases can potentially be attributed to the increase in temperature, which intensify the evaporative demand.

In this study, the evaluation of hydrological changes is focused on the primary components of the hydrological cycle AET and discharge in the Bétérou Catchment. These variables provide a robust basis for interpretation, since they are reasonably well constrained by observations. Other water balance components, such as the partitioning into surface runoff, lateral flow, and baseflow, are not specifically investigated at this stage. Soil moisture dynamics are not analyzed in detail, because no particular calibration or validation has been conducted. Looking ahead, a comprehensive analysis of the water balance components will be conducted with a special focus on soil moisture, the partitioning of the flow components, and hydrological variables.