How far do buckets get us? A flux tracking, topography driven conceptual lumped convolution model to rival the CAOS model
One of the objectives of the DFG (Deutsche Forschungsgemeinschaft) research group "From Catchments as Organised Systems to Models based on Dynamic Functional Units" (CAOS) is the development of the CAOS model, a physically based, flexible hydrological model for the mesoscale (Westhoff and Zehe, 2014). Aim of this substudy is the development of a more conceptual model, which might help to benchmark the CAOS model in terms of accuracy and computational performance.
Study site description
The 1.92 km^2 large Colpach catchment lies within the Attert basin in western Luxembourg directly at the Belgian border. Its elevation ranges between 305 and 538 meters above sea level. The dominating land uses are agricultural areas and grasslands on the plateaus and mainly deciduous forests on the hill slopes. Within the headwaters of the Colpach catchment lies the 47 ha large Weierbach research catchment with an elevation range between 460 and 515 m meters above sea level. The Weierbach catchment is predominantly covered by forest vegetation.
A total amount of ten discharge gauges was installed within the Colpach catchment, subdividing it into several, partly nested, subcatchments (see Tab. 1).
Table 1: study catchments, hierachy (indicates superordinate catchments in descending order) and areas
At each discharge gauge a Teledyne ISCO automatic water sampler was installed. The initial fix sampling routine of one 100 ml discharge sample every 36 h was later on changed to dynamic discharge triggered samples with a temporal resolution of up to two hours. These discharge samples were subsequently analyzed with an PICARRO cavity ring-down spectroscopy analyzer to determine their stable water isotopic compositions.
Simultaneously 23 measurement sites, measuring air temperature, relative humidity, wind speed, solar radiation, soil moisture and groundwater levels, were installed within the Colpach catchment, at open field sites as well as beneath the forest canopy.
At one of the open field sites an automatic precipitation sampler was installed an programmed to take event based precipitation samples, which were analyzed for their isotopic composition. The obtained precipitation isotope data time series could be extended to the past by fortnightly bulk sample data from the Centre de Recherche Publique Gabriel Lippman, which maintains precipitation sampling sites within the Weierbach catchment and at another site next to the Colpach catchment.
We intend to use the available data to develop a lumped catchment model that combines the flux tracking approach described by Hrachowitz et al. (2013) with the topography driven conceptual "FLEX-Topo" modeling approach described by Savenije (2010). Studies by Gharari et al. (2013) and Gao et al. (2014), showed the FLEX-Topo's approach's ability to improve simultaion results by including readily available topographic catchment data into the modeling procedure. However, as soon as a certain complexity is reached, conceptual models face the problems of over-parametrization and equifinality. A modeling approach that allows for flux tracking makes it possible to evaluate a simulation not only by the means of comparing simulated stream discharge responses to observed responses, but also by comparing simulated isotopic compositions to measured isotopic compositions. This will help to constrain the possible model parameter ranges and therefore aid in "getting the right answers for the right reasons" (Kirchner, 2006).
Furthermore, the nested catchment experimental setting creates the opportunity to assess the similarities and differences in catchment functioning and it might be suited to answer the questions of scalability and portability of the results
Methodology of a completed preliminary study
The available data was analyzed according to the "TRANSEP" lumped convolution model approach described by Weiler (2003) and applied in McGuire et al. (2005). In a first step, the effective precipitation was determined by optimizing the parameters of a non linear rainfall loss model (Jakeman and Hornberger, 1993) together with a two parallel linear reservoirs (TPLR) discharge convolution model (Weiler, 2003) to reach an optimal agreement between observed and simulated discharge amounts. The effective precipitation was then used for a mass weighted convolution of the precipitation delta18O values with either a further TPLR or a gamma distribution model (Kirchner, 2000) which was optimized to reach the best agreement between simulated delta18O values in discharge and a series of measured discharge samples. The simulated time period had a length of 12 years.
Results of the completed preliminary study
Figure 1 shows the simulated delta18O values in discharge for the two applied convolution models for the whole Colpach catchment. TPLR and GM models delivered equally appropriate simulations, reaching practically equal Kling-Gupta efficiencies (KGE: Gupta et al. 2009; Kling et al. 2012). Even though the results were equally appropriate, they still were distinctly different. The periods where simulations and observed samples agreed most differed for the two model types. Furthermore, the GM exposed a more distinct short term variability. Eventually, the mean transit times implied by the different models differed considerably.
Figure 2 shows the cumulated response time distributions (based on a TPLR) for all catchments as well as the cumulated transit time distributions (TTDs) according to the two convolution models. The KGE values plotted in Fig. 2 show, that the model failed to simulate the discharge amounts of the subcatchment U appropriately (where, as lately emerged, erroneous discharge data was used) and that
in most cases the three parameter TPLR convolution model reached a higher agreement with the observed isotopic compositions in discharge than the two parameter GM. On the other hand, the three parameters of the TPLR showed a lower identifiability. The overall comparison showed, that the TPLR model tended to imply higher mean transit times (MTTs) than the GM.
Conclusions of the preliminary study
Apparently, the length and probably also the temporal resolution of the time series of isotopic discharge composition have an considerable effect on the outcome of convolution model based transit time modeling. The comparison of two different model types reveal the high uncertainty of their results in respect to TTDs and MTTs, even though the performance of both model types might be considered acceptable for most of the modeled catchments.
Facing the gained results and their uncertainties, the limits of the convolution model approach with static transfer functions based on stable water isotope data become apparent. For a more accurate determination of MTTs and TTDs' tailings, the information content of stable water isotope data is insufficient and the inclusion of a further tracer with extended temporal scope (e.g. 3H) might be favorable (McDonnell et al. ,1999; McDonnell et al.,2010} . For a more appropriate consideration of a catchments short term behavior the assumption of static transit times has to be dropped. A more process based modeling approach which allows for dynamic transit time distributions (e.g. Hrachowitz et al., 2013) seems most promising and is currently being implemented.
Gao, H., Hrachowitz, M., Fenicia, F., Gharari, S., & Savenije, H. H. G. (2014). Testing the realism of a topography-driven model (FLEX-Topo) in the nested catchments of the Upper Heihe, China. Hydrology and Earth System Sciences, 18(5), 1895-1915. doi:10.5194/hess-18-1895-2014
Gharari, S., Hrachowitz, M., Fenicia, F., Gao, H., & Savenije, H. H. G. (2013). Using expert knowledge to increase realism in environmental system models can dramatically reduce the need for calibration. Hydrology and Earth System Sciences Discussions, 10(12), 1480-14855. doi:10.5194/hessd-10-14801-2013
Gupta, H. V., Kling, H., Yilmaz, K. K., & Martinez, G. F. (2009). Decomposition of the mean squared error and NSE performance criteria: Implications for improving hydrological modelling. Journal of Hydrology, 377(1-2), 80-91. doi:10.1016/j.jhydrol.2009.08.003
Hrachowitz, M., Savenije, H. H. G., Bogaard, T. a., Tetzlaff, D., & Soulsby, C. (2013). What can flux tracking teach us about water age distribution patterns and their temporal dynamics? Hydrology and Earth System Sciences, 17(2), 533-564. doi:10.5194/hess-17-533-2013
Kirchner, J. W. (2006). Getting the right answers for the right reasons: Linking measurements, analyses, and models to advance the science of hydrology. Water Resources Research, 42(3), 1-5. doi:10.1029/2005WR004362
Kling, H., Fuchs, M., & Paulin, M. (2012). Runoff conditions in the upper Danube basin under an ensemble of climate change scenarios. Journal of Hydrology, 424-425, 264-277. doi:10.1016/j.jhydrol.2012.01.011
McDonnell, J. J., McGuire, K. J., Aggarwal, P. K., Beven, K. J., Biondi, D., Destouni, G., Wrede, S. (2010). How old is streamwater? Open questions in catchment transit time conceptualization, modelling and analysis. Hydrological Processes, 24(12), 1745-1754. doi:10.1002/hyp.7796
McDonnell, J. J., Rowe, L., & Stewart, M. K. (1999). A combined tracer-hydrometric approach to assess the effect of catchment scale on water flow path, source and age. IAHS-AISH Publication, 265 -273.
Savenije, H. H. G. (2010). HESS Opinions "Topography driven conceptual modelling (FLEX-Topo)."Hydrology and Earth System Sciences, 14(12), 2681-2692. doi:10.5194/hess-14-2681-2010
Weiler, M. (2003). How does rainfall become runoff? A combined tracer and runoff transfer function approach. Water Resources Research, 39(11), 1315. doi:10.1029/2003WR002331