Long-term changes in Central European river discharge 1869-2016: impact of changing snow covers, reservoir constructions and an intensified hydrological cycle

Abstract. Recent climatic changes have the potential to severely alter river runoff, particularly in snow-dominated river basins. Effects of changing snow covers superimpose with changes in precipitation and anthropogenic modifications of the watershed and river network. In the attempt to identify and disentangle long-term effects of different mechanisms, we employ a set of analytical tools to extract long-term changes in river runoff in high resolution. We combine quantile sampling with moving average trend statistics and empirical mode decomposition and apply these tools to discharge data recorded along rivers with 5 nival, pluvial and mixed flow regimes as well as temperature and precipitation data covering the time frame 1869-2016. With a focus on Central Europe, we analyze the long-term impact of snow cover and precipitation changes along with their interaction with reservoir constructions. Our results show that runoff seasonality of snow-dominated rivers decreases. Runoff increases in winter and spring, while discharge decreases in summer and beginning of autumn. We attribute this redistribution of annual flow mainly to reservoir 10 constructions in the alpine ridge. During the course of the last century, large fractions of the alpine rivers have been dammed to produce hydropower. In recent decades, runoff changes induced by reservoir constructions seem to overlap with changes in snow cover. We suggest that alpine signals propagate downstream and affect runoff far outside the alpine area in river segments with mixed flow regimes. Furthermore, our results hint at more (intense) rainfall in recent decades. Detected increases in high discharge can be traced back to corresponding changes in precipitation. 15


Introduction
In many regions of the world, rivers constitute essential lifelines and form the basis of human livelihood. However, recent climate changes may severely affect the hydrological cycle and jeopardize the functional diversity of river systems. Most severe changes are expected to occur in snowdominated river basins. In a warmer world, snow cover characteristics and snowmelt contribution to river runoff will change fundamentally. Rising temperatures are expected to cause less winter precipitation to fall as snow and existing snow covers to melt earlier in spring (Barnett et al., 2005;Simpkins, 2018;Kormann et al., 2015;Birsan et al., 2005). Recent studies suggest that rainfall amount and the number of extreme rainfall events increase due to warmer air holding more water along with enhanced evaporation (Lehmann et al., 2015;Coumou and Rahmstorf, 2012;Mueller and Pfister, 2011). Investigating changes in features of snowpack and snowmelt for key mountain regions, Stewart (2009) summarizes "that both temperature and precipitation increases to date have impacted mountain snowpacks" already. For the Rhine River, one of the most important rivers in Europe, Stahl et al. (2016) indicate that "the influence of climate change is visible particularly in the temporal shifts of sea-sonal minima and maxima of the hydrological regimes of snow and glacier melt dominated alpine headwater catchments." In addition to changes in snowpacks and precipitation, anthropogenic modifications of land surface, subsurface properties and the river network alter river runoff. During the 20th century, more than 45 000 large dams were constructed around the world (Word Commision on Dams, 2000). Also in the Rhine River basin, human activities change runoff with regard to amount, its temporal distribution as well as water quality (Wildenhahn and Klaholz, 1996;Belz et al., 2007;Wildi et al., 2004).
The current knowledge of how climatic changes and changing watershed properties impact river runoff comes largely from instrumental records of hydro-climatic variables, particularly temperature, precipitation and runoff. Birsan et al. (2005) state that "as a spatially integrated variable streamflow is more appealing for detecting regional trends than point measurements of precipitation which is highly variable in space and time", but also point out that watershed properties and their changes over time constitute an "obvious complication in interpreting trends in streamflow data." In addition, quality and length of recorded time series often are insufficient to identify and disentangle effects of the various mechanisms. A sufficient length of the time series inter alia is crucial to be able to distinguish between natural climate variability and signals of climate change. Variability of largescale atmospheric flow on annual to multi-decadal scales, for example, can cause variations in hydro-climatic data, which can either counterbalance or reinforce signals of long-term changes (Hanson et al., 2006;Frei et al., 2000;Kerr, 2000;Scherrer et al., 2016). Studies preparing and investigating long time series of high quality are of great importance and form the basis of our current understanding of features and magnitudes of recent climatic changes (e.g. Vincent et al., 2002;Begert et al., 2005;Schmidli and Frei, 2005;Moberg et al., 2006;Scherrer et al., 2016). Often, simple linear regression approaches are applied to assess characteristics of climatic changes. One frequently used analytical tool in this regard is the robust non-parametric Mann-Kendall trend test (Kendall, 1975;Theil, 1950;Sen, 1968). However, restricting the assessment to linear trends only is hard to justify. The potential of more detailed analyses with regard to seasons, moving time windows (e.g. Kormann et al., 2015) or quantiles of the target variable have hardly been tapped. To further consolidate and extend findings obtained so far, new sets of analytical tools to extract information stored in this time series need to be developed, tested and applied to climatological and hydrological records.
Our study aims at a better understanding of long-term changes in river runoff and identifying potential underlying driving mechanisms, by analysing daily resolution hydroclimatic time series recorded in central Europe between 1869 and 2016. We assess long-term changes in a highly resolved manner by combining quantile sampling, moving average trend statistics and empirical mode decomposition. The two main research question we want to address are the following.
-What is the long-term impact of changes in snow cover on river runoff?
-How do runoff changes induced by changes in snow cover compare with changes caused by reservoir constructions and changes in precipitation?
2 Study area and data We investigate discharge time series from four gauging stations ( Fig. 1 and Table 1). The depicted gauges stand out by the exceptional length of their records and represent different types of flow regimes: nival, pluvial and complex. Gauge Wasserburg is located at the Inn River in Upper Bavaria, Germany. The Inn River is a right tributary of the Danube. The river's source is located in the Swiss Alps and most of its drainage area (1.20 × 10 4 km 2 until gauge Wasserburg) possesses a high Alpine character. The other three gauges investigated, namely Basel, Würzburg and Cologne, are located in the Rhine River basin. The Rhine River is one of the largest rivers in Europe. It is a heavily used waterway and livelihood for the region. At gauge Basel, river runoff is dominated by snowmelt and rainfall runoff from the Alps. Gauge Würzburg is located at the Main River in northern Bavaria, Germany. The Main River is a right tributary of the Rhine River. The catchment area until gauge Würzburg is 1.40 × 10 4 km 2 . The city of Cologne is the largest city along the Rhine River and located in the Lower Rhine region after the confluences with all major tributaries. Until Cologne, the Rhine River drains an area of 1.44 × 10 5 km 2 . For all selected gauges, discharge data at daily resolution have been available since at least 1869. For gauge Basel, statistical tests on daily runoff means conducted by Pfister et al. (2006) show that measured discharge has been homogeneous since 1869 (digitally available part of the time series); i.e. the values are free of anthropogenic effects such as change in instrumentation, change in daily recording frequency or lowering of the river bed (Pfister et al., 2006). Other gauging stations investigated are part of the hydrometric observation network of the water authorities in Germany. Recordings are regularly checked to ensure high quality and reliability. Discharge time series were obtained from the Global Runoff Data Centre (GRDC). Data from the GRDC were used as-is without any further treatment. Elevation distributions and monthly Pardé coefficients for investigated river basins are presented in the Appendix (Fig. A1). Furthermore, we analyse daily resolution temperature and precipitation data provided by the Federal Office of Meteorology and Climatology of Switzerland (MeteoSwiss). At MeteoSwiss, a standardized homogenization procedure is applied to a set of monthly temperature and precipitation time series (Begert et al., 2005). During this homogenization pro-  cedure attained monthly correction values are also applied to daily resolution data. The homogenization of long climatological time series is necessary to correct for non-climatic factors influencing the data. Currently, homogenized daily temperature/precipitation data are available for 28/73 stations. In the following, we focus on meteorological stations where both temperature and precipitation data have been available since at least 1869 and there is no gap in the data longer than 60 d. In total, nine stations fulfil these criteria. Results of the three most prominent stations are displayed and discussed in the main paper ( Fig. 1 and Table 1). Information and results of the remaining stations are given in the Appendix (Table A1 Fig. D1). These strict selection criteria with regard to data length and quality strongly limit the number of recordings suitable for analysis. However, even if a smaller database can reduce the significance of attained results, this is a trade-off we need to accept. Only with recordings having sufficient length and quality can we ensure that the advantages of the proposed analytical tools described in the following chapter can take effect.

Methods
To detect long-term changes in the investigated hydroclimatic data, we combine quantile sampling with moving average trend statistics and empirical mode decomposition (EMD). The selected analytical tools and their combined application to daily time series enable a highly resolved investigation of changes throughout the investigated time frame. The analysis is divided into four steps. Each analysis step complements and extends the information of the previous one, so that step by step, a comprehensive picture of longterm changes takes shape (Fig. 3). A list of all abbreviations and acronyms used can be found in the Appendix (Table B1).

Seasonality of river runoff
To investigate the seasonality of river runoff, we estimate quantiles on a daily basis (QDAY). For every day of the year (DOY), we take all available measurements (i.e. 148 daily values for the period 1869-2016) and calculate QDAYs empirically for probabilities ranging from 0.01 to 0.99. In the framework of this study, quantiles are calculated as the k−1/3 n+1/3 plotting position, with n as the sample size and k = 1, . . .n as the rank (e.g. Hyndman and Fan, 1996). This approach corresponds to type 8 of the "quantile" function in the R environment (R Core Team, 2018).

Changes in seasonality
In order to get a first insight into changes in runoff seasonality, we estimate quantiles from and within a 30 d moving window (QMOV). QMOV operates on discharge data and is independent of previously computed QDAYs. Quantiles are calculated for probabilities between 0.01 and 0.99. The results are continuous quantile time series for each discharge series (Fig. 3). To assess the temporal evolution of these values over the observation period, we employ trend analysis. We calculate trend magnitudes of QMOV using the robust Theil-Sen trend estimator (TST) on a daily basis for all quantiles. Since the computation of the trend uses the values of the same DOY of successive years, auto-correlation should not be of any concern. Within the linear regression approach of the TST, trend magnitudes are estimated as the median slope of ranked data values (Theil, 1950;Sen, 1968;Bronaugh and Werner, 2013).

Onset and evolution of changes
The use of linear trends to quantify the temporal evolution of hydro-climatic variables often lacks physical justification. The respective signals are likely to be non-linear ( Fig. 2c and d). Even when using parametric functions for capturing the non-linear behaviour, e.g. exponential or power-law functions, it is not guaranteed that results will reflect the actual characteristics of underlying processes in the data. An adaptive approach, which does not require a predetermined basis function, is required to get a more flexible characterization of the trend. We employ EMD for this purpose. EMD is an empirical, direct and adaptive method to analyse non-linear trends. It decomposes the signal into oscillatory modes and provides a powerful tool to separate short timescale signals from a general trend (Wu et al., 2007;Huang et al., 1998Huang et al., , 1999Luukko et al., 2016). To avoid mode mixing issues, we performed EMD on an ensemble of the initial data signal: ensemble EMD (EEMD) (Wu and Huang, 2009). Each ensemble member is perturbed by low-amplitude white noise and the results are averaged at the end of the computations.
To keep the characteristics of a complete decomposition, i.e. all extracted intrinsic mode functions (IMFs) sum up to the original signal, the averaging process is carried out separately for each IMF component (Torres et al., 2011). This extension results in a complete EEMD with additive noise (CEEM-DAN). We use an ensemble of 10 000 members, a noise strength of 0.5 times the standard deviation of the input signal, and R package "Rlibeemd" (Luukko et al., 2016) to perform CEEMDAN. The residual of CEEMDAN "can be used to represent the intrinsic trend of the data" (Luukko et al., 2016). We assess these residuals for discharge, temperature and precipitation on a daily basis after calculating moving average values within a window with a width of 30 d for discharge and temperature and a width of 90 d for precipitation. Testing different window sizes, commonly used monthly (30) and seasonal (90) values proved to be a good compromise between robustness of the signal and preservation of signal variability. To make results of different days comparable, we centre each residual by subtracting its mean. To enable the comparison between CEEMDAN residuals and more commonly used linear approaches, we assess whether the nonparametric Mann-Kendall trend test (MK) detects statistical significant monotonic trends in the data CEEMDAN was applied to (α = 0.05) (Mann, 1945;Kendall, 1975). Days with significant monotonic changes are marked with points on top of the respective plot panels (third column in Fig. 4 and columns one and three in Fig. 5).

Changes in quantiles
Furthermore, we investigate changes in quantile magnitudes over time. Therefore, quantiles are estimated on an annual level (QYEA) (Fig. 3). The temporal evolution of QYEAs over the investigated time frame is assessed by applying CEEMDAN. In the case of precipitation, we only use values from "rainy days" (i.e. precipitation > 1 mm). The MK test serves to assess the significance of the trends (marked with points on top of the panel; see Sect. 3.3).

Seasonality of river runoff
Runoff recorded at gauges Wasserburg and Basel is highly seasonal, with high/low runoff during summer/winter ( Fig. 4a1 and b1). Compared to gauge Wasserburg, more runoff is recorded at Basel during winter; i.e. the contrast between summer and winter is less pronounced. At gauge Wasserburg, very high discharge values are almost solely recorded between the months of May and September. Conversely, at gauge Basel, days with very high runoff show up throughout the year. Downstream from gauge Basel, runoff from rain-dominated tributaries such as the Neckar, Main and Mosel blend with Alpine runoff (see gauge Cologne, Fig. 4c1). Rainfall-runoff-dominated basins are characterized by high discharge during winter and the beginning of   spring and low discharge in summer, as seen for Würzburg (Fig. 4d1).

Changes in seasonality
At gauges Wasserburg and Basel, runoff increases during winter and spring for all quantiles, while it decreases during summer and at the beginning of autumn ( Fig. 4a2 and b2). This corresponds to a reduction in runoff seasonality. A very similar overall pattern of changes in runoff can be detected at gauge Cologne: runoff increases during winter and spring and decreases during summer and autumn (Fig. 4c2). In contrast, at gauge Würzburg, discharge quantiles increase throughout, except for high levels at the end of February and March (Fig. 4d2). Similarly to gauges Basel and Cologne, the strongest increases occur during winter.

Onset and evolution of changes
At gauge Wasserburg, pronounced changes in seasonality started in the second half of the 20th century during the 1960s (Fig. 4a3). In contrast, changes at gauge Basel seem to be more gradual and starting earlier in the investigated time period already (Fig. 4b3). At gauge Würzburg, a clear onset of change cannot be detected (Fig. 4d3); however, increases seem to be more uniform and enhanced in recent decades. Patterns of change from snowmelt and rainfall-runoffdominated tributaries overlap at gauge Cologne (Fig. 4c3).
Looking at the respective evolution of potential drivers, temperatures continuously increased throughout the year (Fig. 5a1, b1, and c1). Similar amplitude and interannual patterns are apparent in the three time series. The amount of precipitation increases in recent decades, particularly during winter (Fig. 5a3, b3, and c3). The MK trend test detects significant monotonic increases/decrease in runoff during winter/summer for gauges Wasserburg and Basel ( Fig. 4a3  and b3). For temperature, the MK detects significant increases throughout the year (Fig. 5a1, b1, and c1). Precipitation increases significantly during winter (Fig. 5a3, b3, and c3).

Changes in quantiles
Since the 1960s, QYEAs have strongly increased/decreased at levels below/above 0.6 at gauge Wasserburg (Fig. 4a4). These changes in QYEAs correspond to the strong decrease in seasonality in recent decades (see Sect. 4.2): runoff diminishes in summer and increases in winter. Likewise, at gauge Basel, lower QYEAs (levels < 0.6) increase and higher QYEAs (levels 0.6-0.8) decrease (Fig. 4b4). However, the onset of changes is earlier and changes are smoother compared to detected signals at gauge Wasserburg. Particularly changes in low QYEAs start to increase at the beginning of the investigated time frame already. In contrast to results from gauge Wasserburg, QYEAs at the highest levels (> 0.8) have been increasing at gauge Basel since the 1960s (Fig. 4b4). QYEAs from gauge Würzburg increase over the entire range investigated (Fig. 4d4). Changes in QYEAs below a level of approximately 0.6 occur earlier and are smoother than for higher levels. There, the increases have been enhanced in recent decades. At gauge Cologne, high QYEAs have increased in recent decades (Fig. 4c4), making it similar to findings from gauges Basel and Würzburg. Also, the lower QYEAs have experienced an increase. This increase, however, is not a gradual one over the entire time frame, but is rather a U-shaped process (decline until the 1940s, then increase).
For precipitation, similar pattern seem to show up. Increasing QYEAs hint at more (intense) rainfall in recent decades (Fig. 5a4, b4 and c4). Increases in QYEAs in temperature seem to occur earlier and seem to be enhanced at lower temperatures (Fig. 5a2, b2 and c2). Changes in quantiles for individual seasons are given in the Appendix (Figs. B1 and C1). At gauge Wasserburg, changes in QYEAs are significant according to the MK (Fig. 4a4). The more a trend pattern deviates from a monotonic increase and more U-shaped signals emerge, the more often the MK results in non-significant p values. The main results depicted in Fig. 4 are summarized in Table 2.

Seasonality of river runoff
Runoff at gauge Wasserburg is dominated by the accumulation and depletion of a seasonal snowpack. The intra-annual variability of runoff is very high, and high flows mainly occur during the snowmelt season and during summer, when higher temperatures enable liquid precipitation in large fractions of the catchment (Fig. 4a1). There are no bigger lakes that could attenuate flood or low-flow events generated in the basin. In comparison, large lakes constitute an important element of the Rhine River basin until gauge Basel. Furthermore, large parts of the basin are sub-Alpine terrain. As a result, liquid rainfall is an important streamflow component throughout the year and runoff less seasonal compared to gauge Wasserburg (Stahl et al., 2016). A detailed overview of hydrological regimes in Switzerland and their characteristics can inter alia be found in Weingartner and Aschwanden (1992) or Speich et al. (2015). Reconstructing the largest flood events in the High Rhine basin since 1268, Wetter et al. (2011) indicate that about half of all major floods occur during summer. Flood events during summer usually are the result of high baseflow due to a melting Alpine snow cover superimposing with heavy rainfall (Wetter et al., 2011). Extreme flood events during autumn, winter and spring often are caused by long-lasting precipitation events coinciding with strong snowmelt due to rain on snow (RoS) and/or a temporary temperature increase (Wetter et al., 2011;Schmocker-Fackel and Naef, 2010). For higherelevated river basins, RoS events play an important role in runoff formation (Sui and Koehler, 2001;Merz and Blöschl, 2003). The RoS flood that occurred in the Bernese Alps, Switzerland, in October 2011 showed how damaging these kinds of events can be (Rössler et al., 2014). Another example is the RoS events from January 2011, where rainfall released vast amounts of water stored in a temporary snow cover and caused RoS-driven flood events in the whole of central Europe (Freudiger et al., 2014). However, the importance of different flood-generating mechanisms is changing with recent climatic changes (see e.g. Blöschl et al., 2017Blöschl et al., , 2019Berghuijs et al., 2019).
Even though about one-third of the runoff in the Main River originates from snowmelt (Stahl et al., 2016), there is only little impact of snow accumulation and snowmelt on the seasonal distribution of discharge. It seems that low temperatures rarely prevail long enough to enable the accumulation and preservation of snow over a longer period. Runoff is dominated by large-scale rainfall events occurring in winter and increased evapotranspiration during summer (Fig. 4d1). At gauge Cologne, we have the situation of superimposing nival and pluvial runoff components (Fig. 4c1). This overlap results in a more uniform seasonal distribution of discharge. High QDAYs are higher during winter, whereas low QDAYs are higher during summer (Fig. 4c1). This reversal in the seasonal distribution hints at the importance of different flow components for different flow situations. Runoff due to large-scale rainfall events over the middle and lower parts of the catchment are important for high discharge values, particularly during winter. During summer, snowmelt and glacier melt from the Alpine part of the basin play an important role in the sustenance of runoff in the lower reaches of the Rhine River (Stahl et al., 2016).

Changes in seasonality
In the snow-dominated river basins Wasserburg and Basel, the seasonality of river runoff decreases over the investigated time frame. For the increasing runoff values during winter and early spring, several mechanisms have to be taken into account. First of all, changes in the Alpine snow cover have to be considered. In recent decades, rising temperatures have caused less snow accumulation during winter (Laternser and Schneebeli, 2003;Marty, 2008;Scherrer et al., 2004;Wielke et al., 2004;Marty et al., 2017). Thus, a greater fraction of total precipitation is liquid and reaches the river system without being stored in snowpacks. In addition, the frequency of days with temperatures above 0 • C increases, causing parts of any existing snow cover to melt (Scheifinger et al., 2003;Kreyling and Henry, 2011;Zubler et al., 2014;Schädler and Weingartner, 2010). Rising temperatures also result in shorter snow duration, where "shorter snow duration is mainly caused by earlier snow melting in spring than by later first snowfalls in autumn" (Laternser and Schneebeli, 2003). The earlier onset of snowmelt in spring represents a much-noticed effect of rising temperatures on Alpine river runoff (e.g. Kormann et al., 2015;Birsan et al., 2005;Stewart, 2009).
Less snow accumulation during the preceding winter results in lower discharges during the following melting period, i.e. late spring and early summer. Furthermore, recent studies suggest that rising temperatures might lead to a reduction in snowmelt rates (Musselman et al., 2017;Wu et al., 2018). "Slower snowmelt in a warmer world may decrease the likelihood that wetness thresholds that permit hydrologic connectivity will be exceeded, leading to spring and summer stream- Gradual + U-shape flow declines and lower runoff efficiency" (Musselman et al., 2017). Changes in the liquid/solid fraction of precipitation overlap with changes in the total amount of rainfall. We observed increased rainfall during winter for all stations investigated. Likewise, numerous other studies point to a recent increase in precipitation, particularly during winter (e.g. Begert et al., 2005;Scherrer et al., 2016;Frei and Schär, 2001). However, increasing catchment evaporation due to increasing radiation, air temperature and vegetation activity might at least partly compensate detected changes in precipitation (Duethmann and Blöschl, 2018;Schädler and Weingartner, 2010;Norris and Wild, 2007;Wild et al., 2007).
At gauge Cologne, we also detect a decrease in discharge during summer and autumn (Fig. 4c2). We hypothesize that this decrease is the result of a downstream propagation of the Alpine signal, possibly overlapping with increasing evaporation rates in the basin. The decrease in summer discharge in the Lower Rhine cannot be attributed to reduced ice melt contributions from the Alpine glaciers. Assessing the snow and glacier melt components of streamflow of the Rhine River for the time frame 1901-2006, Stahl et al. (2016 showed that "despite the glacier retreat the modelled ice melt component of the streamflow in the Rhine does not show a strong long-term trend over the entire study period, i.e. a systematic decline or increase of this component. The detailed results of the modelling suggests that an increased ice melt due to increased temperature may have been compensated by the reduction in glacier area". Gauge Würzburg, with its discharge increasing throughout the entire year, does not show any detectable changes in seasonality.

Onset and evolution of changes
Investigating long-term snow trends of the Swiss Alps, Laternser and Schneebeli (2003) suggest that "mean snow depth, the duration of continuous snow cover and the number of snowfall days in the Swiss Alps all show very similar trends during the observation period 1933-1999: a gradual increase until the early 1980s (with significant interruptions during the late 1950s and early 1970s) followed by a statistically significant decrease towards the end of the century".
At gauge Basel, these changes in Alpine snowpacks seem to be insufficient to explain the decrease in runoff seasonality detected. In particular, winter discharge (low QYEAs) increases already from the beginning of the investigated time frame on ( Fig. 4b3 and b4). Instead, we suspect anthropogenic alterations of the river network, particularly reservoir constructions, to be an important driver. These might have caused the redistribution of water from summer to winter earlier in the investigated time frame already. Large fractions of the Swiss and Austrian Alpine river systems have been dammed to produce hydropower. The two Alpine countries have the highest specific hydroelectric production per surface area globally (Truffer et al., 2001). The first hydropower station in Switzerland was constructed in 1899 (Verbunt et al., 2005). Dam constructions in the Alpine Rhine and along other Alpine rivers, such as the Aare, Limmat and Reuss, gained momentum in the 1920s, and most of the large storage lakes were constructed between 1950 and 1970 (Meile et al., 2011;Wildenhahn and Klaholz, 1996;Wagner et al., 2015;Bosshard et al., 2013) (Fig. E1). The total storage volume of large storage lakes (river weirs not included) of the High Rhine/entire Rhine basin is estimated to amount to 1.86/3.12 × 10 9 m 3 (Wildenhahn and Klaholz, 1996). In order to ensure full functional capability of highhead storage hydropower stations, reservoirs need to have sufficient water volume stored at all times. Therefore, reservoirs tend to be filled during summer when discharge is high. Conversely, storages are depleted during low flow in winter (Belz et al., 2007;Meile et al., 2011;Farinotti et al., 2016;Wesemann et al., 2018). A rough estimation supports this notion: assuming the 1.86 × 10 9 m 3 of storage is emptied between December and April (and filled between June and October), mean runoff would increase/decrease by approximately 10 m 3 s −1 decade −1 in these months during the investigated time frame, which corresponds to the trend magnitudes depicted in Fig. 4b2. In addition to reservoir constructions, regulations of lake levels and routing of rivers through lakes, e.g. the diversion of the Aare River into Lake Biel in 1887 (part of the First Jura-Waters Corrections), need to be considered (Wetter et al., 2011).
At gauge Wasserburg, pronounced changes in runoff seasonality do not show up until the second half of the 20th century ( Fig. 4a3 and a4). In the Inn basin, the constructions of key reservoirs, such as the Gepatsch reservoir (Tyrol, Austria), the Lago di Livigno reservoir (Grisons, Switzerland, and Lombardy, Italy) and the Lai da Ova Spin compensation reservoir (Grisons, Switzerland), were not completed until the 1960s. The construction of those big reservoirs coincides with the detected onset of changes in river runoff. We suspect that also in the Inn River basin, the construction and management of reservoirs for hydropower might be an important factor changing seasonality of river runoff. In addition to changes in seasonality, the operation of high-head hydropower stations causes unnatural fluctuation on (sub-)daily timescales (hydropeaking) (Meile et al., 2011;Pérez Ciria et al., 2019) (see also Fig. F1). Effects of reservoirs possibly overlap (with) changes induced by changes in snow cover.
Also, rainfall-runoff-dominated rivers, such as the Main River at Würzburg, are strongly affected by hydroengineering installations. One large-scale project inaugurated in 1992 after numerous decades of constructions is the Rhine-Main-Danube waterway. In order to raise low water discharge in the Main River, about 1.55 × 10 8 m 3 (3.50 × 10 8 m 3 ) of water is transferred on average per year (in a dry year) from the Danube into the Main River basin via the Main-Danube Canal (Maniak, 2016). The connection from the Rhine River until Würzburg with constructions of weirs to regulate the river´s water level was completed in the 1940s (Wirth, 1995). This onset of water-level regulations in the 1940s coincides with increasing low QYEAs (< 0.6) at gauge Würzburg (Fig. 4d4). We suspected that anthropogenic alterations strongly impacted the discharge of the Main River, particularly during low-discharge periods. However, they seem to be insufficient to explain changes in higher QYEAs.

Changes in quantiles
We detect increasing high QYEAs at gauges Basel, Cologne and Würzburg (Fig. 4b4, c4, and d4). Possible driving mechanisms might be changes in precipitation: our results hint at more (intense) rainfall in recent decades (Fig. 5). In the following, we discuss possible underlying forcing mechanisms of detected signals. These include changes in largescale circulation patterns, solar dimming/brightening and temperature-moisture feedbacks.
Long-term changes in the occurrence frequencies and/or characteristics of circulation patterns are known to have a strong impact on local climate. More frequent zonal circulation in winter since the 1970s, for example, might be responsible for "more frequent mild and humid winters in Central Europe" (Bárdossy and Caspary, 1990). This increase in zonal circulation follows upon several decades with in-creased numbers of blocking days during winter (Häkkinen et al., 2011). Blocking in the Atlantic region is anti-correlated with phases of the North Atlantic Oscillation (NAO) (Scherrer et al., 2006;Stein, 2000;Pavan et al., 2000). Negative values of the NAO index "indicate periods of reduced north-south pressure gradient, reduced westerly winds and weaker advection of warm oceanic air onto the cold European landmass" (Parker et al., 2007). Wintertime NAO on the other hand is influenced by the Atlantic Multidecadal Oscillation (AMO): a "positive phase of the AMO results in more frequent negative NAO" (Peings and Magnusdottir, 2014). The AMO depicts multi-decadal (60-70 years) variations in sea surface temperatures in the North Atlantic basin (Peings and Magnusdottir, 2014;Kerr, 2000). After several decades of warm anomalies of sea surface temperatures (positive phase of AMO, which coincided with more frequent blocking days), the North Atlantic started to cool down and to transition into a negative AMO phase in the 1960s (Peings and Magnusdottir, 2014;Häkkinen et al., 2011). This transition into a negative AMO phase and fewer blocking days coincides with more frequent zonal circulation in winter (Bárdossy and Caspary, 1990), more (intense) rainfall (Fig. 5a3, b3, c3, a4, b4, and c4) and an increase in discharge, whereas changes seem to overlap with changes induced by anthropogenic alteration of the river network and changes in Alpine snowpacks (Fig. 4b4, c4, and d4).
Generally, detected patterns in temperature time series investigated in the framework of this study are consistent across all stations (Figs. 5 and D1). In the case of precipitation, overall patterns are similar; however, stronger variations among stations show up. Precipitation is subject to stronger local and regional variability than temperature. This evidently limits the informative value of precipitation recorded at individual points for discussions on catchment scale. However, variations of local meteorological variables are strongly influenced by large-scale flow and regional-scale weather patterns (Scherrer et al., 2016;Murawski et al., 2018;Weusthoff, 2011). Weather patterns/types represent specific synoptic conditions and lead to certain meteorological conditions in a region. For western Germany, precipitation correlates well on scales of hundreds of kilometres, particularly during winter (Schönwiese and Rapp, 1997). During summer, local convective storms are an important source of rainfall, and strong differences in amount and intensity over very short distances are possible (Sodemann and Zubler, 2010;Lavers et al., 2013). However, such convective storms "are hardly of any relevance for the formation of floods in the large river basins of Central Europe, because the extent of convective rainstorms is restricted to local occurrence" (Bronstert et al., 2007). Rainfall-runoff processes on larger scales are dominated by advective precipitation. The main moisture source then is the Atlantic Ocean (Sodemann and Zubler, 2010). Following the above-mentioned aspects, we hypothesize that even if superimposed by local variability due to smaller-scale processes and regional variations due to general precipitation gradients, long-term signals in precipitation detected on point scale can provide important information for discussions on catchment scale. However, caution has to be exercised, results of available stations compared and findings not transferred to places outside the region of influence. Marty (2008) relate detected shifts in snow days in Switzerland to an enhanced temperature increase due to changes in circulation patterns coinciding with "the full magnitude of the greenhouse effect, which is no longer masked by solar dimming". After a multi-decadal decrease from about the 1950s to the 1980s (solar dimming), recent decades saw an increase in regional solar irradiance (solar brightening) due to decreasing amounts of anthropogenic aerosols in the atmosphere (Ruckstuhl et al., 2008;Norris and Wild, 2007;Ruckstuhl and Norris, 2009).
Rapidly rising temperatures have the capacity to affect the entire hydrological cycle. A feedback mechanism of major importance in this respect is the temperature-moisture feedback: rising temperatures result in increasing evaporation and precipitation, which in turn leads to an intensification of the entire hydrological cycle (Huntington, 2006;Held and Soden, 2000). Lehmann et al. (2015) indicate that a "thermally driven moisture increase has significantly contributed to the intensification of extreme rainfalls since the 1980s".
Against the background of recent changes in temperature, precipitation, and frequencies in zonal circulation, and following Labat et al. (2004), we suggest that rapid increases in temperatures in recent decades have resulted in an increased sea-land transport of moisture and increases in precipitation and runoff. Signals possibly overlap with changes in moisture transport due to varying frequencies in zonal circulation and blocking days.
At gauge Wasserburg, in contrast to other gauges investigated, high QYEAs do not increase. We suspect that in this case, high discharges are also controlled by snowmelt processes rather than by liquid rainfall. Furthermore, investigated rain gauges might not depict changes in precipitation in the complex Alpine topography of the catchment. In the High Rhine basin up to gauge Basel, large fractions of the basin are located outside the Alpine ridge and liquid precipitation plays an important role throughout the year. Therefore, impacts of changes in snow cover, reservoir constructions and effects of more (intense) rainfall all seem to be detectable in measured discharge (Fig. 4b4). Changes in snow cover and river regulations, which decrease runoff seasonality, seem to primarily affect QYEAs below a level of 0.85, and more intense rainfall events increase the magnitude of higher QYEAs.

Conclusions
We investigate daily observational data from key river gauges and meteorological stations located in central Europe cover-ing the time frame 1869-2016. Investigated time series stand out by the exceptional length and quality of their continuous recordings. A cascading sequence of analytical tools is used to extract high-resolution signals of long-term changes. In order to acquire a comprehensive picture of long-term changes, we combine quantile sampling with moving average trend statistics and empirical mode decomposition. Given that the recordings have sufficient length and quality, presented tools enable investigations of high resolution and provide detailed insights into underlying trend patterns. A very high quality of the time series is required to prevent non-climatic factors, such as changes in observation practices or site relocation, from affecting the determination of trends (Begert et al., 2005;Scherrer et al., 2016;Begert and Frei, 2018). A sufficient length of the time series is vital to be able to distinguish between natural climate variability and signals of climate change. When adopting strict criteria regarding data length and data quality, the number of stations suitable for analysis strongly decreases. A small database limits the significance of attained results, though. This is a trade-off we need to accept in order to ensure full functionality of presented analytical tools. Consistent results for stations and parameters investigated make us confident that even if the number of stations is limited, we will attain meaningful results that are worth discussing.
The seasonality of the analysed snow-dominated rivers decreases. We suspect river regulations, particularly reservoir constructions, to be the main driver of detected changes. Reservoirs are filled during summer when discharge is high and storages depleted during low flow in winter (e.g. Belz et al., 2007;Meile et al., 2011). In recent decades, runoff changes induced by reservoir constructions seem to overlap with changes in Alpine snowpacks (Laternser and Schneebeli, 2003;Scherrer et al., 2004). Rising temperature reduces seasonal snow covers and the seasonal redistribution of runoff from winter to summer. An exact separation of effects of reservoirs and changes in snow cover and investigations of possible counterbalancing interactions are still pending and are the focus of future research. Furthermore, we suspect that detected decreases in discharge during summer and autumn in the Lower Rhine region at gauge Cologne are the result of a downstream propagation of the Alpine signal, possibly further overlapping with increasing evaporation rates in the basin.
In addition, our results hint at more (intense) rainfall in recent decades, particularly during winter. Detected changes in precipitation seem to intensify high discharges. The detected increase in precipitation (intensity) is not a gradual one over the entire time frame, but rather follows a U-shape (decline until the 1940s, then increase). Further research is necessary to pin down underlying mechanisms of detected changes in precipitation and runoff. We suspect that detected signals might be due to an increase in sea-land transport of moisture, particularly during winter, being part of a recent intensification of the entire hydrological cycle (Huntington, 2006;Held and Soden, 2000;Lehmann et al., 2015;Labat et al., 2004). Temperature-driven increases in moisture and precipitation possibly overlap with natural multi-decadal variations in sea-land-moisture transport (Parker et al., 2007;Kerr, 2000;Häkkinen et al., 2011;Peings and Magnusdottir, 2014;Bárdossy and Caspary, 1990;Scherrer et al., 2006;Pavan et al., 2000).
Over recent decades, hydrological regimes have been changing at a very fast pace. Some progress has been made in extracting long-term signals of change in hydroclimatic data. However, further studies investigating longterm changes in river runoff focusing on the detection of underlying mechanisms and the disentanglement of their effects are of great urgency and importance.
Appendix A Figure A1. Elevation distribution (raster cells at 500 m resolution calculated based on EU-DEM v.1.1 by the EU Copernicus Programme) and Pardé coefficients (mean monthly discharge divided by the mean annual discharge) (Pardé, 1933;Spreafico and Weingartner, 2005) for investigated river basins Wasserburg, Basel, Cologne and Würzburg.     Appendix F Figure F1. Raster hydrograph for gauge Wasserburg. In recent decades, hydropeaking (weekly pattern) due to the operation of high-head storage hydropower stations imprinted. Time frame displayed: 1869-2016.
Data availability. Climatological data used in this study were obtained from the Federal Office of Meteorology and Climatology of Switzerland, CH-8058 Zurich Airport (MeteoSwiss). Discharge data analysed were provided by the Global Runoff Data Centre, 56068 Koblenz, Germany (GRDC). Data analysis was carried out using the R statistical software (https://www.r-project.org/, last access: 7 April 2020; R Core Team, 2018).
Author contributions. ER conducted the analysis and wrote the manuscript. TF, GB and AB provided guidance in the process of data analysis and preparation of the manuscript.
Competing interests. The authors declare that they have no conflict of interest.