These files contain the model results underlying the paper "Reduction of vegetation-accessible water storage capacity after deforestation affects catchment travel time distributions and increases young water fractions in a headwater catchment" (Hrachowitz et al., 2021; HESS):
(1) Model results calibrated for the 2009-2013 pre-deforestation period ("1_Predeforestation_2009_2013")
(2) Model results from the application of the model calibrated to the pre-deforestation period and without re-calibration applied to the 2013-2016 post deforestation period ("2_Predeforestation_model_applied_to_postdeforestation_2013_2016")
(3) Model results calibrated for the 2013-2016 post-deforestation period ("3_Postdeforestation_2013_2016")

each of the 3 zip files contains:
(a) Modelled daily fluxes for all feasible parameter sets as defined in the paper (Figure 4):
Water fluxes [mm/d] in the following files: "PR" = rain; "PM" = melt; "EP" = potential evaporation; "PEH" = hillslope throuhgfall; "PER" = riparian zone throughfall; "QS" = groundwater flow; "QR" = shallow lateral flow from riparian zone; "Q" = total flow; "EA" = total actual evaporation 
Tracer [%o] in the following files: "O18_Q" = O-18 in stream flow; "O18_QS" = O-18 in groundwater flow; "O18_QR" = O-18 in lateral flow from riparian zone; "O18_EA" = O-18 in total actual evaporation 
Note: first column is the date in Excel format. Second column are the results of the best performing parameter set.

(b) Modelled daily states for all feasible parameter sets as defined in the paper (Figure 4):
Water [mm]: "Ssnow" = snowwater equivalent storage; "SIH" = hillslope interception storage; "SIR" = riparian interception storage; "SUH" = hillslope unsaturated rootzone storage; "SUR" = riparian unsaturated rootzone storage; "SSa" = hydrologically active groundwater storage; "SFR" = fast responding riparian zone storage  
Note: first column is the date in Excel format. Second column are the results of the best performing parameter set.

(c) Feasible Parameter sets as defined in the paper (Table 2) in file "Par":
Column 1: melt factor FM [mm d^-1 C^-1]
Column 2: upwelling groundwater factor fQS [-]
Column 3: hillslope interception capacity ImaxH [mm]
Column 4: riparian interception capacity ImaxR [mm]
Column 5: storage coefficent of fast responding riparian storage kR [d^-1]
Column 6: storage coefficent of groundwater storage kS [d^-1]
Column 7: transpiration water stress factor Lp [-]
Column 8: maximum percolation rate Rs,max [mm d^-1]
Column 9: hillslope unsaturated root zone storage capacity SUmaxH [mm]
Column 10: riparian unsaturated root zone storage capacity SUmaxR [mm]
Column 11: Threshold temperature TT [C]
Column 12: shape factor gamma [-]
Column 13: shape factor for SAS function alpha0 [-]
Column 14: hydrologically passive groundwater reservoir Ssp [mm]
Note: first row is the best performing parameter set. 

(d) Performance metrics for feasible Parameter sets as defined in the paper (Table 3) in file "obj_fun":
Column 1: ENS,Q
Column 2: ENS,log(Q)
Column 3: EV,Q
Column 4: ENS,FDC
Column 5: ENS,FDCh
Column 6: ENS,PD
Column 7: ER,RLD
Column 8: ER,DLD
Column 9: ENS,AC
Column 10: ER,AC1
Column 11: ER,AC1l
Column 12: ER,CR
Column 13: ENS,d18O
Column 14: ER,RD
Note: first row are the performance metrics for the best performing parameter set. 

(e) Modelled hydrological signatures using the feasible parameter sets as defined in the paper (Table 3) in the following files: 
"FDC": flow duration curve (Column 1: exceedance probability [-], Columns 2-last: Q [mm/d])
"FDCh": high flow duration curve (Column 1: exceedance probability [-], Columns 2-last: Q [mm/d])
"PD": peak distribution (Colums 1-last: Q [mm/d])
"RLD": rising limb density [-]
"DLD": declining limb density [-]
"AC": autocorrelation function (Column 1: lag times [d]; Column 2-last: correlation coefficient r [-])
"AC1": lag-1 autocorrelation [-]
"AC1l": low flow lag-1 autocorrelation [-]
"CR": Runoff ratio [-]


