Refine
Year of publication
Document Type
- Article (32)
- Postprint (11)
- Doctoral Thesis (1)
- Review (1)
Language
- English (45)
Is part of the Bibliography
- yes (45)
Keywords
- Nordeste (4)
- Precipitation (4)
- time-series (4)
- variability (4)
- Connectivity (3)
- Erosion (3)
- seasonality (3)
- snow (3)
- Climate (2)
- Forecasting Framework (2)
Water shortage is a serious threat for many societies worldwide. In drylands, water management measures like the construction of reservoirs are affected by eroded sediments transported in the rivers. Thus, the capability of assessing water and sediment fluxes at the river basin scale is of vital importance to support management decisions and policy making. This subject was addressed by the DFG-funded SESAM-project (Sediment Export from large Semi-Arid catchments: Measurements and Modelling). As a part of this project, this thesis focuses on (1) the development and implementation of an erosion module for a meso-scale catchment model, (2) the development of upscaling and generalization methods for the parameterization of such model, (3) the execution of measurements to obtain data required for the modelling and (4) the application of the model to different study areas and its evaluation. The research was carried out in two meso-scale dryland catchments in NE-Spain: Ribera Salada (200 km²) and Isábena (450 km²). Adressing objective 1, WASA-SED, a spatially semi-distributed model for water and sediment transport at the meso-scale was developed. The model simulates runoff and erosion processes at the hillslope scale, transport processes of suspended and bedload fluxes in the river reaches, and retention and remobilisation processes of sediments in reservoirs. This thesis introduces the model concept, presents current model applications and discusses its capabilities and limitations. Modelling at larger scales faces the dilemma of describing relevant processes while maintaining a manageable demand for input data and computation time. WASA-SED addresses this challenge by employing an innovative catena-based upscaling approach: the landscape is represented by characteristic toposequences. For deriving these toposequences with regard to multiple attributes (eg. topography, soils, vegetation) the LUMP-algorithm (Landscape Unit Mapping Program) was developed and related to objective 2. It incorporates an algorithm to retrieve representative catenas and their attributes, based on a Digital Elevation Model and supplemental spatial data. These catenas are classified to provide the discretization for the WASA-SED model. For objective 3, water and sediment fluxes were monitored at the catchment outlet of the Isábena and some of its sub-catchments. For sediment yield estimation, the intermittent measurements of suspended sediment concentration (SSC) had to be interpolated. This thesis presents a comparison of traditional sediment rating curves (SRCs), generalized linear models (GLMs) and non-parametric regression using Random Forests (RF) and Quantile Regression Forests (QRF). The observed SSCs are highly variable and range over six orders of magnitude. For these data, traditional SRCs performed poorly, as did GLMs, despite including other relevant process variables (e.g. rainfall intensities, discharge characteristics). RF and QRF proved to be very robust and performed favourably for reproducing sediment dynamics. QRF additionally excels in providing estimates on the accuracy of the predictions. Subsequent analysis showed that most of the sediment was exported during intense storms of late summer. Later floods yielded successively less sediment. Comparing sediment generation to yield at the outlet suggested considerable storage effects within the river channel. Addressing objective 4, the WASA-SED model was parameterized for the two study areas in NE Spain and applied with different foci. For Ribera Salada, the uncalibrated model yielded reasonable results for runoff and sediment. It provided quantitative measures of the change in runoff and sediment yield for different land-uses. Additional land management scenarios were presented and compared to impacts caused by climate change projections. In contrast, the application for the Isábena focussed on exploring the full potential of the model's predictive capabilities. The calibrated model achieved an acceptable performance for the validation period in terms of water and sediment fluxes. The inadequate representation of the lower sub-catchments inflicted considerable reductions on model performance, while results for the headwater catchments showed good agreement despite stark contrasts in sediment yield. In summary, the application of WASA-SED to three catchments proved the model framework to be a practicable multi-scale approach. It successfully links the hillslope to the catchment scale and integrates the three components hillslope, river and reservoir in one model. Thus, it provides a feasible approach for tackling issues of water and sediment yield at the meso-scale. The crucial role of processes like transmission losses and sediment storage in the river has been identified. Further advances can be expected when the representation of connectivity of water and sediment fluxes (intra-hillslope, hillslope-river, intra-river) is refined and input data improves.
The Southern Pre-Pyrenees experienced a substantial land-use change over the second half of the 20th century owing to the reduction of agricultural activities towards the formation of a more natural forest landscape. The land-use change over the last 50 years with subsequent effects on water and sediment export was modelled with the process-based, spatially semi-distributed WASA-SED model for the meso-scale Canalda catchment in Catalonia, Spain. It was forwarded that the model yielded plausible results for runoff and sediment yield dynamics without the need of calibration, although the model failed to reproduce the shape of the hydrograph and the total discharge of several individual rainstorm events, hence the simulation capabilities are not yet considered sufficient for decision-making purposes for land management. As there are only a very limited amount of measured data available on sediment budgets with altered land-use and climate change settings, the WASA-SED model was used to obtain qualitative estimates on the effects of past and future change scenarios to derive a baseline for hypothesis building and future discussion on the evolution of sediment budgets in such a dryland setting. Simulating the effects of the past land-use change, the model scenarios resulted in a decrease of up to 75% of the annual sediment yield. whereas modelled runoff remained almost constant over the last 50 years. The relative importance of environmental change was evaluated by comparing the impact on sediment export of land-use change, that are driven by socio-economic factors, with climate change projections for changes in the rainfall regime. The modelling results suggest that a 20% decrease in annual rainfall results in a decrease in runoff and sediment yield, thus an ecosystem stabilisation in regard to sediment export which can only be achieved by a substantial land-use change equivalent to a complete afforestation. At the same time, a 20% increase in rainfall causes a large export of water and sediment resources out of the catchment, equivalent to an intensive agricultural use of 100% of the catchment area. For wet years, the effects of agricultural intensification are more pronounced, so that in this case the intensive land-use change has a significantly larger impact on sediment generation than climate change. The WASA-SED model proved capable in quantifying the impacts of actual and potential environmental change, but the reliability of the simulation results is still circumscribed by considerable parameterisation and model uncertainties.
Suspended sediment transport in a highly erodible catchment : the River Isábena (Southern Pyrenees)
(2009)
Understanding and quantifying sediment load is important in catchments draining highly erodible materials that eventually contribute to siltation of downstream reservoirs. Within this context, the suspended sediment transport and its temporal dynamics have been studied in the River Isabena (445 km(2), south-central Pyrenees, Ebro basin) by means of direct sampling and turbidity recording during a 3-year dry period. The average flood-suspended sediment concentration was 8 g l(-1). with maximum instantaneous values above 350 g l(-1). The high scatter between discharge and suspended sediment concentrations (up to five orders of magnitude) has not permitted the use of rating curve methods to estimate the total load. Interpolation techniques yielded a mean annual sediment load of 184,253 t y(-1) for the study period, with a specific yield of 414 t km(-2) y(-1). This value resembles those reported for small torrents in nearby mountainous environments and is the result of the high connectivity between the badland source areas and stream courses, a fact that maximises sediment conveyance through the catchment. Floods dominated the sediment transport and yield. However, sediment transport was more constant through time than that observed in Mediterranean counterparts; this can be attributed to the role of base flows that entrain fine sediment temporarily stored in the channel and force the river to carry high sediment concentrations (i.e., generally in the order of 0.5 g l(-1)), even under minimum flow conditions.
Forests seem to represent low-erosion systems, according to most, but not all, studies of suspended-sediment yield. We surmised that this impression reflects an accidental bias in the selection of monitoring sites towards those with prevailing vertical hydrological flowpaths, rather than a tight causal link between vegetation cover and erosion alone. To evaluate this conjecture, we monitored, over a 2-year period, a 3.3 ha old-growth rainforest catchment prone to frequent and widespread overland flow. We sampled stream flow at two and overland flow at three sites in a nested arrangement on a within-event basis, and monitored the spatial and temporal frequency of overland flow. Suspended-sediment concentrations were modeled with Random Forest and Quantile Regression Forest to be able to estimate the annual yields for the 2 years, which amounted to 1 t ha(-1) and 2 t ha(-1) in a year with below-average and with average precipitation, respectively. These estimates place our monitoring site near the high end of reported suspended-sediment yields and lend credence to the notion that low yields reflect primarily the dominance of vertical flowpaths and not necessarily and exclusively the kind of vegetative cover. Undisturbed forest and surface erosion are certainly no contradiction in terms even in the absence of mass movements.
In this paper, we analyse the effectiveness of flood management measures based on the concept known as "retaining water in the landscape". The investigated measures include afforestation, micro-ponds and small-reservoirs. A comparative and model-based methodological approach has been developed and applied for three meso-scale catchments located in different European hydro-climatological regions: Poyo (184 km(2)) in the Spanish Mediterranean, Upper Iller (954 km(2)) in the German Alps and Kamp (621 km(2)) in Northeast-Austria representing the Continental hydro-climate. This comparative analysis has found general similarities in spite of the particular differences among studied areas. In general terms, the flood reduction through the concept of "retaining water in the landscape" depends on the following factors: the storage capacity increase in the catchment resulting from such measures, the characteristics of the rainfall event, the antecedent soil moisture condition and the spatial distribution of such flood management measures in the catchment. In general, our study has shown that, this concept is effective for small and medium events, but almost negligible for the largest and less frequent floods: this holds true for all different hydro-climatic regions, and with different land-use, soils and morphological settings.
This study focuses on the prediction of event-based runoff coefficients (an important descriptor of flood events) for nested catchments up to an area of 50?km(2) in the Eastern Ore Mountains. The four main objectives of the study are (i) the prediction of runoff coefficients with the statistical method of generalized linear models, (ii) the comparison of the results of the linear models with estimates of a distributed conceptual model, (iii) the comparison of the dynamics of observed soil moisture and simulated saturation deficit of the hydrological model and (iv) the analysis of the relationship between runoff coefficient and observed and simulated wetness. Different predictor variables were selected to describe the runoff coefficient and were differentiated into variables describing the catchment’s antecedent wetness and meteorological forcing. The best statistical model was estimated in a stepwise approach on the basis of hierarchical partitioning, an exhaustive search algorithm and model validation with jackknifing. We then applied the rainfall runoff model WaSiM ETH to predict the runoff processes for the two larger catchments. Locally measured small-scale soil moisture (acquired at a scale of four to five magnitudes smaller than the catchment) was identified as one of the key predictor variables for the estimation of the runoff coefficient with the general linear model. It was found that the relationship betweenobserved and simulated (using WaSiM ETH) wetness is strongly hysteretic. The runoff coefficients derived from the rainfall runoff simulations systematically underestimate the observed values. Copyright (C) 2012 John Wiley & Sons, Ltd.
The paper presents the sediment budget of the Isabena basin, a highly dynamic 445-km(2) catchment located in the Central Pyrenees that is patched by highly erodible areas (i.e., badlands). The budget for the period 2007-2009 is constructed following a methodology that allows the interpolation of intermittent measurements of suspended sediment concentrations and enables a subsequent calculation of sediment loads. Data allow specification of the contribution of each subbasin to the water and sediment yield in the catchment outlet. Mean annual sediment load was 235,000 t y(-1). Specific sediment yield reached 2000 t km(-2) y(-1), a value that indicates very high sedimentary activity, especially in the case of Villacarli and Lascuarre subcatchments, were most badlands are located. The specific sediment yield obtained for the entire Isabena is 527 t km(-2) y(-1), a high value for such a mesoscale basin. Results show that a small part of the area (i.e., 1%) controls most of the catchment's gross sediment contribution. Sediment delivery ratio (ratio between sediment input from primary sources and basin export) has been estimated at around 90%, while in-channel storage represents the 5% of the annual load on average. The high connectivity between sediment sources (i.e., badlands) and transfer paths (i.e., streamcourses) exacerbates the influence of the local sediment production on the catchment's sediment yield, a quite unusual fact for a basin of this scale.
Saturated hydraulic conductivity (K-s) is an important soil characteristic affecting soil water storage, runoff generation and erosion processes. In some areas where high-intensity rainfall coincides with low K-s values at shallow soil depths, frequent overland flow entails dense drainage networks. Consequently, linear structures such as flowlines alternate with inter-flowline areas. So far, investigations of the spatial variability of K-s mainly relied on isotropic covariance models which are unsuitable to reveal patterns resulting from linear structures. In the present study, we applied two sampling approaches so as to adequately characterize K-s spatial variability in a tropical forest catchment that features a high density of flowlines: A classical nested sampling survey and a purposive sampling strategy adapted to the presence of flowlines. The nested sampling approach revealed the dominance of small-scale variability, which is in line with previous findings. Our purposive sampling, however, detected a strong spatial gradient: surface K-s increased substantially as a function of distance to flowline; 10 m off flowlines, values were similar to the spatial mean of K-s. This deterministic trend can be included as a fixed effect in a linear mixed modeling framework to obtain realistic spatial fields of K-s. In a next step we used probability maps based on those fields and prevailing rainfall intensities to assess the hydrological relevance of the detected pattern. This approach suggests a particularly good agreement between the probability statements of K-s exceedance and observed overland flow occurrence during wet stages of the rainy season.
Soils in various places of the Panama Canal Watershed feature a low saturated hydraulic conductivity (K-s) at shallow depth, which promotes overland-flow generation and associated flashy catchment responses. In undisturbed forests of these areas, overland flow is concentrated in flow lines that extend the channel network and provide hydrological connectivity between hillslopes and streams. To understand the dynamics of overland-flow connectivity, as well as the impact of connectivity on catchment response, we studied an undisturbed headwater catchment by monitoring overland-flow occurrence in all flow lines and discharge, suspended sediment, and total phosphorus at the catchment outlet. We find that connectivity is strongly influenced by seasonal variation in antecedent wetness and can develop even under light rainfall conditions. Connectivity increased rapidly as rainfall frequency increased, eventually leading to full connectivity and surficial drainage of entire hillslopes. Connectivity was nonlinearly related to catchment response. However, additional information on factors such as overland-flow volume would be required to constrain relationships between connectivity, stormflow, and the export of suspended sediment and phosphorus. The effort to monitor those factors would be substantial, so we advocate applying the established links between rain event characteristics, drainage network expansion by flow lines, and catchment response for predictive modeling and catchment classification in forests of the Panama Canal Watershed and in similar regions elsewhere.
Deforestation is a prominent anthropogenic cause of erosive overland flow and slope instability, boosting rates of soil erosion and concomitant sediment flux. Conventional methods of gauging or estimating post-logging sediment flux often focus on annual timescales but overlook potentially important process response on shorter intervals immediately following timber harvest. We resolve such dynamics with non-parametric quantile regression forests (QRF) based on high-frequency (3 min) discharge measurements and sediment concentration data sampled every 30-60 min in similar-sized (similar to 0.1 km(2)) forested Chilean catchments that were logged during either the rainy or the dry season. The method of QRF builds on the random forest algorithm, and combines quantile regression with repeated random sub-sampling of both cases and predictors. The algorithm belongs to the family of decision-tree classifiers, which allow quantifying relevant predictors in high-dimensional parameter space. We find that, where no logging occurred, similar to 80% of the total sediment load was transported during extremely variable runoff events during only 5% of the monitoring period. In particular, dry-season logging dampened the relative role of these rare, extreme sediment-transport events by increasing load efficiency during more efficient moderate events. We show that QRFs outperform traditional sediment rating curves (SRCs) in terms of accurately simulating short-term dynamics of sediment flux, and conclude that QRF may reliably support forest management recommendations by providing robust simulations of post-logging response of water and sediment fluxes at high temporal resolution.
The 933 km(2) Bengue catchment in northeastern Brazil is characterized by distinct rainy and dry seasons. Precipitation is stored in variously sized reservoirs, which is essential for the local population. In this study, we used TerraSAR-X SM(HH) data for an one-year monitoring of seasonal changes in the reservoir areas from July 2011 to July 2012. The monitoring was based on acquisitions in the ascending pass direction, complemented by occasional descending-pass images. To detect water surface areas, a histogram analysis followed by a global threshold classification was performed, and the results were validated using in situ GPS data. Distinguishing between small reservoirs and similar looking dark areas was difficult. Therefore, we tested several approaches for identifying misclassified areas. An analysis of the surface area dynamics of the reservoirs indicated high spatial and temporal heterogeneities and a large decrease in the total water surface area of the reservoirs in the catchment by approximately 30% within one year.
In a study from 2008, Lariviere and colleagues showed, for the field of natural sciences and engineering, that the median age of cited references is increasing over time. This result was considered counterintuitive: with the advent of electronic search engines, online journal issues and open access publications, one could have expected that cited literature is becoming younger. That study has motivated us to take a closer look at the changes in the age distribution of references that have been cited in water resources journals since 1965. Not only could we confirm the findings of Lariviere and colleagues. We were also able to show that the aging is mainly happening in the oldest 10-25% of an average reference list. This is consistent with our analysis of top-cited papers in the field of water resources. Rankings based on total citations since 1965 consistently show the dominance of old literature, including text books and research papers in equal shares. For most top-cited old-timers, citations are still growing exponentially. There is strong evidence that most citations are attracted by publications that introduced methods which meanwhile belong to the standard toolset of researchers and practitioners in the field of water resources. Although we think that this trend should not be overinterpreted as a sign of stagnancy, there might be cause for concern regarding how authors select their references. We question the increasing citation of textbook knowledge as it holds the risk that reference lists become overcrowded, and that the readability of papers deteriorates.
To support scientifically sound water management in dryland environments a modelling system has been developed for the quantitative assessment of water and sediment fluxes in catchments, transport in the river system, and retention in reservoirs. The spatial scale of interest is the mesoscale because this is the scale most relevant for management of water and land resources.
This modelling system comprises process-oriented hydrological components tailored for dryland characteristics coupled with components comprising hillslope erosion, sediment transport and reservoir deposition processes. The spatial discretization is hierarchically designed according to a multi-scale concept to account for particular relevant process scales. The non-linear and partly intermittent run-off generation and sediment dynamics are dealt with by accounting for connectivity phenomena at the intersections of landscape compartments. The modelling system has been developed by means of data from nested research catchments in NE-Spain and in NE-Brazil.
In the semi-arid NE of Brazil sediment retention along the topography is the main process for sediment retention at all scales, i.e. the sediment delivery is transport limited. This kind of deposition retains roughly 50 to 60 % of eroded sediment, maintaining a similar deposition proportion in all spatial scales investigated. On the other hand, the sediment retained in reservoirs is clearly related to the scale, increasing with catchment area. With increasing area, there are more reservoirs, increasing the possibility of deposition. Furthermore, the area increase also promotes an increase in flow volume, favouring the construction of larger reservoirs, which generally overflow less frequently and retain higher sediment fractions. The second example comprises a highly dynamic Mediterranean catchment in NE-Spain with nested sub-catchments and reveals the full dynamics of hydrological, erosion and deposition features. The run-off modelling performed well with only some overestimation during low-flow periods due to the neglect of water losses along the river. The simulated peaks in sediment flux are reproduced well, while low-flow sediment transport is less well captured, due to the disregard of sediment remobilization in the riverbed during low flow.
This combined observation and modelling study deepened the understanding of hydro-sedimentological systems characterized by flashy run-off generation and by erosion and sediment transport pulses through the different landscape compartments. The connectivity between the different landscape compartments plays a very relevant role, regarding both the total mass of water and sediment transport and the transport time through the catchment.
The Isabena catchment (445 km(2)), Spain, features highly diverse spatial heterogeneity in land use, lithology and rainfall. Consequently, the relative contribution in terms of water and sediment yield varies immensely between its subcatchments, and also temporally. This study presents the synthesis of similar to 2.5 years of monitoring rainfall, discharge and suspended sediment concentration (SSC) in the five main subcatchments of the Isabena and its outlet.
Continuous discharge at the subcatchment outlets, nine tipping bucket rainfall and automatic SSC samplers (complemented by manual samples), were collected from June 2011 until November 2013. The water stage records were converted to discharge using a rating curve derived with Bayesian regression. For reconstructing sediment yields, the data from the intermittent SSC sampling needed to be interpolated. We employed non-parametric multivariate regression (Quantile Regression Forests, QRF) using the discharge and rainfall data plus different aggregation levels of these as ancillary predictors. The subsequent Monte Carlo simulations allowed the determination of monthly sediment yields and their uncertainty.
The Isabena catchment shows high erosion dynamics with great variability in space and time, with stark contrasts even between adjacent subcatchments. The natural conditions make water and sediment monitoring and instrumentation very challenging; the measurement of discharge is particularly prone to considerable uncertainties. The QRF method employed for reconstructing sedigraphs and monthly yields proved well suited for the task.
The results of streamflow trend studies are often characterized by mostly insignificant trends and inexplicable spatial patterns. In our study region, Western Austria, this applies especially for trends of annually averaged runoff. However, analysing the altitudinal aspect, we found that there is a trend gradient from higher-altitude to lower-altitude stations, i.e. a pattern of mostly positive annual trends at higher stations and negative ones at lower stations. At mid-altitudes, the trends are mostly insignificant. Here we hypothesize that the streamflow trends are caused by the following two main processes: on the one hand, melting glaciers produce excess runoff at higher-altitude watersheds. On the other hand, rising temperatures potentially alter hydrological conditions in terms of less snowfall, higher infiltration, enhanced evapotranspiration, etc., which in turn results in decreasing streamflow trends at lower-altitude watersheds. However, these patterns are masked at mid-altitudes because the resulting positive and negative trends balance each other. To support these hypotheses, we attempted to attribute the detected trends to specific causes. For this purpose, we analysed trends of filtered daily streamflow data, as the causes for these changes might be restricted to a smaller temporal scale than the annual one. This allowed for the explicit determination of the exact days of year (DOYs) when certain streamflow trends emerge, which were then linked with the corresponding DOYs of the trends and characteristic dates of other observed variables, e.g. the average DOY when temperature crosses the freezing point in spring. Based on these analyses, an empirical statistical model was derived that was able to simulate daily streamflow trends sufficiently well. Analyses of subdaily streamflow changes provided additional insights. Finally, the present study supports many modelling approaches in the literature which found out that the main drivers of alpine streamflow changes are increased glacial melt, earlier snowmelt and lower snow accumulation in wintertime.
Owing to average temperature increases of at least twice the global mean, climate change is expected to have strong impacts on local hydrology and climatology in the Alps. Nevertheless, trend analyses of hydro-climatic station data rarely reveal clear patterns concerning climate change signals except in temperature observations. However, trend research has thus far mostly been based on analysing trends of averaged data such as yearly, seasonal or monthly averages and has therefore often not been able to detect the finer temporal dynamics. For this reason, we derived 30-day moving average trends, providing a daily resolution of the timing and magnitude of trends within the seasons. Results are validated by including different time periods. We studied daily observations of mean temperature, liquid and solid precipitation, snow height and runoff in the relatively dry central Alpine region in Tyrol, Austria. Our results indicate that the vast majority of changes are observed throughout spring to early summer, most likely triggered by the strong temperature increase during this season. Temperature, streamflow and snow trends have clearly amplified during recent decades. The overall results are consistent over the entire investigation area and different time periods.
The results of streamflow trend studies are often characterized by mostly insignificant trends and inexplicable spatial patterns. In our study region, Western Austria, this applies especially for trends of annually averaged runoff. However, analysing the altitudinal aspect, we found that there is a trend gradient from higher-altitude to lower-altitude stations, i.e. a pattern of mostly positive annual trends at higher stations and negative ones at lower stations. At midaltitudes, the trends are mostly insignificant. Here we hypothesize that the streamflow trends are caused by the following two main processes: on the one hand, melting glaciers produce excess runoff at higher-altitude watersheds. On the other hand, rising temperatures potentially alter hydrological conditions in terms of less snowfall, higher infiltration, enhanced evapotranspiration, etc., which in turn results in decreasing streamflow trends at lower-altitude watersheds. However, these patterns are masked at mid-altitudes because the resulting positive and negative trends balance each other. To support these hypotheses, we attempted to attribute the detected trends to specific causes. For this purpose, we analysed trends of filtered daily streamflow data, as the causes for these changes might be restricted to a smaller temporal scale than the annual one. This allowed for the explicit determination of the exact days of year (DOYs) when certain streamflow trends emerge, which were then linked with the corresponding DOYs of the trends and characteristic dates of other observed variables, e.g. the average DOY when temperature crosses the freezing point in spring. Based on these analyses, an empirical statistical model was derived that was able to simulate daily streamflow trends sufficiently well. Analyses of subdaily streamflow changes provided additional insights. Finally, the present study supports many modelling approaches in the literature which found out that the main drivers of alpine streamflow changes are increased glacial melt, earlier snowmelt and lower snow accumulation in wintertime.
The results of streamflow trend studies are often characterized by mostly insignificant trends and inexplicable spatial patterns. In our study region, Western Austria, this applies especially for trends of annually averaged runoff. However, analysing the altitudinal aspect, we found that there is a trend gradient from higher-altitude to lower-altitude stations, i.e. a pattern of mostly positive annual trends at higher stations and negative ones at lower stations. At midaltitudes, the trends are mostly insignificant. Here we hypothesize that the streamflow trends are caused by the following two main processes: on the one hand, melting glaciers produce excess runoff at higher-altitude watersheds. On the other hand, rising temperatures potentially alter hydrological conditions in terms of less snowfall, higher infiltration, enhanced evapotranspiration, etc., which in turn results in decreasing streamflow trends at lower-altitude watersheds. However, these patterns are masked at mid-altitudes because the resulting positive and negative trends balance each other. To support these hypotheses, we attempted to attribute the detected trends to specific causes. For this purpose, we analysed trends of filtered daily streamflow data, as the causes for these changes might be restricted to a smaller temporal scale than the annual one. This allowed for the explicit determination of the exact days of year (DOYs) when certain streamflow trends emerge, which were then linked with the corresponding DOYs of the trends and characteristic dates of other observed variables, e.g. the average DOY when temperature crosses the freezing point in spring. Based on these analyses, an empirical statistical model was derived that was able to simulate daily streamflow trends sufficiently well. Analyses of subdaily streamflow changes provided additional insights. Finally, the present study supports many modelling approaches in the literature which found out that the main drivers of alpine streamflow changes are increased glacial melt, earlier snowmelt and lower snow accumulation in wintertime.
The characteristics of a landscape pose essential factors for hydrological processes. Therefore, an adequate representation of the landscape of a catchment in hydrological models is vital. However, many of such models exist differing, amongst others, in spatial concept and discretisation. The latter constitutes an essential pre-processing step, for which many different algorithms along with numerous software implementations exist. In that context, existing solutions are often model specific, commercial, or depend on commercial back-end software, and allow only a limited or no workflow automation at all. Consequently, a new package for the scientific software and scripting environment R, called lumpR, was developed. lumpR employs an algorithm for hillslope-based landscape discretisation directed to large-scale application via a hierarchical multi-scale approach. The package addresses existing limitations as it is free and open source, easily extendible to other hydrological models, and the workflow can be fully automated. Moreover, it is user-friendly as the direct coupling to a GIS allows for immediate visual inspection and manual adjustment. Sufficient control is furthermore retained via parameter specification and the option to include expert knowledge. Conversely, completely automatic operation also allows for extensive analysis of aspects related to landscape discretisation. In a case study, the application of the package is presented. A sensitivity analysis of the most important discretisation parameters demonstrates its efficient workflow automation. Considering multiple streamflow metrics, the employed model proved reasonably robust to the discretisation parameters. However, parameters determining the sizes of subbasins and hillslopes proved to be more important than the others, including the number of representative hillslopes, the number of attributes employed for the lumping algorithm, and the number of sub-discretisations of the representative hillslopes.
The characteristics of a landscape pose essential factors for hydrological processes. Therefore, an adequate representation of the landscape of a catchment in hydrological models is vital. However, many of such models exist differing, amongst others, in spatial concept and discretisation. The latter constitutes an essential pre-processing step, for which many different algorithms along with numerous software implementations exist. In that context, existing solutions are often model specific, commercial, or depend on commercial back-end software, and allow only a limited or no workflow automation at all.
Consequently, a new package for the scientific software and scripting environment R, called lumpR, was developed. lumpR employs an algorithm for hillslope-based landscape discretisation directed to large-scale application via a hierarchical multi-scale approach. The package addresses existing limitations as it is free and open source, easily extendible to other hydrological models, and the workflow can be fully automated. Moreover, it is user-friendly as the direct coupling to a GIS allows for immediate visual inspection and manual adjustment. Sufficient control is furthermore retained via parameter specification and the option to include expert knowledge. Conversely, completely automatic operation also allows for extensive analysis of aspects related to landscape discretisation.
In a case study, the application of the package is presented. A sensitivity analysis of the most important discretisation parameters demonstrates its efficient workflow automation. Considering multiple streamflow metrics, the employed model proved reasonably robust to the discretisation parameters. However, parameters determining the sizes of subbasins and hillslopes proved to be more important than the others, including the number of representative hillslopes, the number of attributes employed for the lumping algorithm, and the number of sub-discretisations of the representative hillslopes.