Effects of Bias-Correcting Climate Model Data on the Projection of Future Changes in High Flows

Bias-correction methods are commonly applied to climate model data in hydrological climate impact studies. This is due to the often large deviations between simulated and observed climate variables. These biases may cause unrealistic simulation results when directly using the climate model data as input for hydrological models. Our analysis of the EURO-CORDEX (Coordinated Downscaling Experiment for Europe) data for the Northwestern part of Germany showed substantial biases for all climatological input variables needed by the hydrological model PANTA RHEI. The sensitivity for climatological input data demonstrated that changes in only one climate variable significantly affect the simulated average discharge and mean annual peak flow. The application of bias correction methods of different complexity on the climate model data improved the plausibility of hydrological modeling results for the historical period 1971–2000. The projections for the future period 2069–2099 for high flows indicate on average small changes for representative concentration pathway (RCP) 4.5 and an increase of approximately 10% for RCP8.5 when applying non-bias corrected climate model data. These values significantly differed when applying bias correction. The bias correction methods were evaluated in terms of their ability to (a) maintain the change signal for precipitation and (b) the goodness of fit for hydrological parameters for the historical period. Our results for this evaluation indicated that no bias correction method can explicitly be preferred over the others.


Introduction
An intensification of the hydrological cycle due to climate change is likely, but varies in dependence of the geographic region and spatial scale [1]. This implies precipitation, flood intensities and frequencies as well as the occurrence and magnitude of droughts and low water situations. Changes in intensities and frequencies of precipitation related to global warming investigated both on the global [2] and regional scale [3][4][5] compared climate impact projections of river floods in Europe under climate change. They concluded that flood risk will consistently increase in most of Central and Western Europe. But in addition to this, investigations on smaller spatial scales are needed to identify flood endangered areas due to regional varying effects of climate change [6]. These studies are necessary in order to develop specific adaptation measures of public authorities and local governments [7].
When assessing the impacts of global change on mesoscale catchments, the general procedure is the application of a chain of several models/steps: (1) general circulation models (GCM), (2) regional climate models (RCMs) or statistical downscaling, (2b) bias-correction (BC) and (3) hydrological models. Due to the coarse spatial resolution of GCMs of 1-2 • [1], the direct application of these data for hydrological models is only useful in large-scale catchments. At mesoscale level, step (2) is necessary in • Sensitivity analysis for assessing the necessity of bias correcting various climate variables when using climate scenario data as input for hydrological models. • Projecting future changes in high flow parameters in six medium-scale catchments in the Northwestern part of Germany.

•
Evaluating the performance of BC methods of different complexity.

Observational Data
The observed meteorological input data (precipitation, mean, minimum and maximum temperature, relative humidity, global radiation and wind speed) was available on a daily time step for the period 1951-2015 from the Lower Saxony Water Management, Coastal Defence and Nature Conservation Agency (NLWKN). Data from 771 stations for precipitation and 123 stations for the other meteorological variables was used to generate 1 × 1 km gridded data fields for Lower Saxony by means of geo-statistical analyses. This was executed by Institute of Water Resources Management, Hydrology and Agriculture (WAWI) of the Leibniz University of Hannover in the framework of the research project KliBiW. These fine-resolution grids were used as input data for the calibration and validation of the hydrological model.
For the comparison with the RCM data for the historical period 1971-2000, the observational 1 × 1 km grid was aggregated to a 10 × 10 km grid. In the following, this grid data is referred to as observational data X obs .
Observed discharge data, Q obs , was available on a daily basis for at least the year 1977 up to 2011 from the NLWKN.

Coordinated Downscaling Experiment for Europe (CORDEX) Data
The RCM data were obtained from the CORDEX (Coordinated Downscaling Experiment for Europe) project [17]. The 17 GCM/RCM combinations, hereafter referred to as climate scenarios, were composed of six different GCMs and six different RCMs (Table 1). Table 1. General circulation model/regional climate model (GCM/RCM) combination matrix. S7 and S8, as well as S16 and S17 are combinations of the same GCM and RCM. They are different in terms of the GCM run with which the RCM was driven. CNRM_CM5 S1 S12 EC-Earth S2 S7, S8 S10 S13 IPSL-CM5A-MR S3 S6 HadGEM2-ES S4 S9 S14 MPI-ESM-LR S5 S15 S16, S17 NorESM1-M S11

Institution
In this study, two RCP scenarios were considered: RCP4.5 and RCP8.5. The name of the RCP refers to the projected additional radiative forcing in W/m 2 [18]. Additionally, historical runs for the climate scenarios were also available. All climate scenarios have a spatial resolution of 0.11 • (about 12.5 km) and a temporal resolution of one day.
Most of the climate variables needed by the hydrological model were available from the CORDEX project. Relative humidity was not provided by all climate scenarios. It was generally calculated in dependence of specific humidity and daily mean temperature [19] as the ratio of actual water vapor pressure to saturated water vapor pressure. Before further processing, a reprojection of the CORDEX data was necessary in order to comply with the coordinate system and resolution of X obs (see Figure 1). This was done using the inverse distance method with the consideration of three neighbours.

Study Area
The study area comprises six mesoscale catchments with spatial extends of 309 to 1731 km 2 in the Northwestern part of Germany (see Figure 1). The region is mostly flat with a mean elevation of the catchments between 45 and 210 m above sea level.
The annual precipitation for the period 1971-2000 varies between 620 and 820 mm for the six catchments. The annual mean of temperature is similar in all catchments (8.6-9.6 • C).
The texture of the soil in the catchments Emlichheim and Hellwege is sandy whereas in the catchments of Derneburg and Wunstorf, the soil texture is dominated by silt. In Harxbüttel and Bramsche, loams with both fine and coarse texture size can be found.
The catchment land use is mainly agricultural with winter crops and corn being the predominant crop types. The annual precipitation for the period 1971-2000 varies between 620 and 820 mm for the six catchments. The annual mean of temperature is similar in all catchments (8.6-9.6 °C).
The texture of the soil in the catchments Emlichheim and Hellwege is sandy whereas in the catchments of Derneburg and Wunstorf, the soil texture is dominated by silt. In Harxbüttel and Bramsche, loams with both fine and coarse texture size can be found.
The catchment land use is mainly agricultural with winter crops and corn being the predominant crop types.

General Procedure
Projections of 17 climate scenarios and two RCP scenarios for the future period 2070-2099, Xfut, as well as corresponding data for the historical period 1971-2000, Xhist, were used as input for the hydrological model PANTA RHEI in order to identify hydrological climate change signals with a focus on high flow parameters. High precipitation values being a crucial factor for high flows are also investigated in terms of climate change signals. All time periods are based on the hydrological year and comprise 30 years (e.g., 11/1970-10/2000). The influence of BC was assessed by correcting climate variables with different methods and comparing the resulting parameters for the historical period and the climate change signals with those obtained by non-corrected climate scenario data. Non-corrected climate scenario data are referred to as Xhist,orig/Xfut,orig, bias-corrected data are Xhist,bc/Xfut,bc.
For the computation of climatological parameters and the corresponding climate change signals, all grid cells with the center within the catchment boundaries or within a distance of maximum 5 km from the catchment boundaries were considered. The number of grid cells per catchment varied between 6 and 31. The total number of grid cells for all catchments was 93.
First, the necessity of the application of BC for the different climate variables was identified by comparing Xhist,orig to Xobs in terms of annual mean values and monthly deviations. The number of wet days and parameters for high precipitation events were also considered. In accordance with the majority of studies reviewed [20,21] a threshold of 1 mm/d was applied to define a wet day (WD1). In terms of high precipitation events, the 95% percentile of (a) daily precipitation on wet days (P95) and (b) the 3-day precipitation sum above 3 mm (P95d3) were compared. The overall result of the comparison between Xhist,orig and Xobs is calculated as the average of all 93 grid cells. Simulation results of the calibrated hydrological model forced with both Xobs and Xhist,orig were analysed

General Procedure
Projections of 17 climate scenarios and two RCP scenarios for the future period 2070-2099, X fut , as well as corresponding data for the historical period 1971-2000, X hist , were used as input for the hydrological model PANTA RHEI in order to identify hydrological climate change signals with a focus on high flow parameters. High precipitation values being a crucial factor for high flows are also investigated in terms of climate change signals. All time periods are based on the hydrological year and comprise 30 years (e.g., 11/1970-10/2000). The influence of BC was assessed by correcting climate variables with different methods and comparing the resulting parameters for the historical period and the climate change signals with those obtained by non-corrected climate scenario data. Non-corrected climate scenario data are referred to as X hist,orig/ X fut,orig , bias-corrected data are X hist,bc/ X fut,bc .
For the computation of climatological parameters and the corresponding climate change signals, all grid cells with the center within the catchment boundaries or within a distance of maximum 5 km from the catchment boundaries were considered. The number of grid cells per catchment varied between 6 and 31. The total number of grid cells for all catchments was 93.
First, the necessity of the application of BC for the different climate variables was identified by comparing X hist,orig to X obs in terms of annual mean values and monthly deviations. The number of wet days and parameters for high precipitation events were also considered. In accordance with the majority of studies reviewed [20,21] a threshold of 1 mm/d was applied to define a wet day (WD1). In terms of high precipitation events, the 95% percentile of (a) daily precipitation on wet days (P95) and (b) the 3-day precipitation sum above 3 mm (P95d3) were compared. The overall result of the comparison between X hist,orig and X obs is calculated as the average of all 93 grid cells. Simulation results of the calibrated hydrological model forced with both X obs and X hist,orig were analysed regarding the deviations in average discharge (MQ) and the mean of annual maximum discharge (MHQ). Hereafter, the simulated discharge applying X obs /X hist,orig/ X hist,bc is referred to as Q ref /Q hist,orig/ Q hist,bc . Discharge simulations for the future period are correspondingly abbreviated with Q fut,orig and Q fut,bc .
The BC effects on the simulation of high flow parameters during the historical period were assessed by comparing the parameters MHQ and flood events with a 100-year return period (HQ100) of Q ref and Q hist,bc . The parameter HQ100 is estimated using the Gumbel distribution and L-Moments for the parameter estimation of the function. In terms of the climate change signals, the parameters P95d3 and MHQ were compared. For P95d3, the mean value of all grid cells was used, for MHQ, it is the area-weighted mean of all six catchments.

Hydrological Model PANTA RHEI
The hydrological model PANTA RHEI has been developed by the Department of Hydrology, Water Management and Water Protection, Leichtweiß Institute for Hydraulic Engineering and Water Resources, University of Braunschweig in cooperation with the "Institut für Wassermanagement IfW GmbH", Braunschweig [22]. It has been successfully applied for scientific questions in numerous national and international projects [7,23,24]. Furthermore, PANTA RHEI is applied in the operational flood forecast of the Federal State of Lower Saxony, Germany [25].
PANTA RHEI is a deterministic, semi-distributed, physically based hydrological model with a three level spatial disaggregation scheme: catchment, subcatchments and hydrotopes. Within a subcatchment, hydrotopes are homogeneous regarding the soil type and land use type.
The internal temporal discretization of PANTA RHEI is one hour. However, the time step of climatological input data can be different. For this study, discharge data were evaluated on a daily basis. A detailed description of PANTA RHEI is given by [26].

Calibration and Validation
The six catchments of the study area were calibrated and validated with observational climate and discharge data for the period 1971-2011. The first 20 years were used for calibration and the last 20 years for validation. Due to data availability, for the gauges Derneburg and Harxbüttel, a slightly shorter period was used.
The calibration was done semi-automatically applying the lexicographic strategy developed by [26] using four criteria: the percentage deviation of the total discharge volume, the mean monthly deviation of discharge, the Nash-Sutcliffe efficiency [27] and the root mean squared error (RMSE) of the maximum annual discharge (HQa). The mean monthly deviation is calculated by the absolute values of percentage deviations.
The Nash-Sutcliffe efficiencies for both periods range between 0.79 and 0.91 with a mean of 0.83 for calibration as well as for validation, which can be denoted as values of high quality. The deviation in total discharge volume amounts to a maximum of 8.6% during both periods. The mean monthly deviation of discharge varies between 3% and 16% and is slightly higher for the validation period for all catchments.

Sensitivity Analysis for Climate Variables
The sensitivity of discharge in dependence of climate variables was carried out for all climate variables that serve as input for the hydrological model. The variables of X obs were altered in order to assess the impact of changes in one variable on the simulation of discharge. The alteration of X obs was carried out either by multiplication with a factor (wind speed, humidity, radiation, precipitation) or by addition of a value (temperature). The level of change is derived from the deviations between X hist,orig and X obs for each variable during the historical period. Therefore, three sensitivity levels (SL) were defined. The SL take into account the mean deviations and the standard deviation of the deviations between X hist,orig and X obs (see Equations (1)-(6)) for the two examples of precipitation and temperature).
where P is precipitation, T is temperature, d is the day, µ is the mean of observational data, µ s and σ s are the mean and the standard deviation of climate variables of all climate scenarios considered. The adjustment of the variables wind speed, humidity and radiation was undertaken analogous to that for precipitation. In terms of relative humidity, values that amounted to >100% were set to 100%. The daily minimum and maximum temperature were adjusted with the same deltas as mean temperature, so that the temperature amplitudes of the X hist,orig were not changed.
One climate variable was altered according to one sensitivity level for each simulation run of the hydrological model, while X obs was used for the other input variables. In addition to these local sensitivity studies, one global sensitivity simulation run was carried out with alterations of all variables according to SL1.
The sensitivity of discharge in dependence of climate variables was evaluated for the two parameters MQ and MHQ.

Bias Correction for Hydrological Modeling
The main focus when comparing different BC methods was on precipitation as it is the most critical factor for the deviations between Qref and Qhist,orig (see Section 4.2). Various methods with different complexity were applied to Phist,orig. All other climate variables were corrected with the linear scaling method.
The following short description of BC methods therefore focuses on the adjustment of precipitation although some methods are also applicable to other climate variables. As a part of some methods, the wet-day frequency is adjusted, as it is overestimated by most of the climate scenarios. In this study, this is generally done by applying the approach described by [21] as part of the local intensity scaling (LOCI) method with a 1mm threshold for a wet day. In contrast to a 0 mm threshold, which is used by [12], this also allows the adjustment in case that the wet-day frequency is underestimated by climate scenarios. In our study, this method is part of the BC methods LOCI, power transformation and distribution mapping. For all BC methods applies that Xfut,orig, is corrected the same way as the corresponding data Xhist,orig and all corrections are carried out on a monthly basis.

Linear Scaling
Linear scaling is a well-established BC method. The adjustment is made based on the difference in monthly mean values between X obs and X hist,orig . A correction factor calculated out of these differences is applied to X hist,orig and X fut,orig either by summation (temperature, Equation (7)) or multiplication (all other variables, see Equation (8) for the example of precipitation): where "bc" indicates the bias corrected data, µ is the mean and i is the month (1 . . . 12). The LOCI correction, is described in detail e.g., by [21] and implies the adjustment of the wet-day frequency. This is done by defining a wet-day threshold (WDT) from P hist,orig in such a way that the number of days in that time series matches the number of days in P obs exceeding 1mm. In a second and third step, rain intensities on wet days is calculated (Equation (9)) and the time series is adjusted (Equation (10)).
where s is the scaling factor and WDT is the wet-day threshold.

Power Transformation
Additional to the adjustment of wet days and rain-intensities, the power transformation, described by [28] also addresses the variance of precipitation data. As proposed by [12,29], the wet-day frequency was adjusted in advance. Then, P hist,orig was corrected applying the parameters a and b: First, the parameter b is determined using a root-finding algorithm so that the coefficient of variance of the adjusted P hist matches the coefficient of variance of P obs . The parameter a is determined afterwards in order to adjust the mean.

Distribution Mapping
The main principle is the identification of a transfer function in order to adjust the distribution function of P hist . For precipitation, the Gamma distribution is suggested [30]: where α is the shape parameter and β is the scale parameter. The application of the distribution mapping approach for adjusting climate scenario data is described in detail by [12]. The limitation of this method is that the precipitation data can be well approximated with the chosen distribution function. This should be checked in advance.

Inter-Sectoral Impact-Model Intercomparison Project (ISI-MIP) Approach
This bias correction method was applied within the Inter-Sectoral Impact-Model Intercomparison Project (ISI-MIP) and is described by [31]. Within this method, the daily variability of P hist,orig is adjusted, in addition to the correction of long-term monthly mean values and wet-day frequency. The ISI-MIP approach preserves the relative changes for monthly mean precipitation. Two versions of this BC method are described: the FAST-TRACK approach and the extended approach. In this study, the latter one is applied. This implicates the adjustment of the daily variability by using a nonlinear transfer function. The linear option is chosen, if the non-linear adjustment does not provide better results than a linear fitting. The correction procedure is carried out as follows: In a first step, the long-term monthly mean values are adjusted in the same way it was done for the linear scaling method.
In the second step, the frequency of wet days is adjusted. A wet-day threshold is defined from P hist,orig analogous to the procedure of the LOCI method, but with the difference that all values below the threshold are set to zero which can reduce monthly precipitation. In order to compensate for that, the amount from dry days of a specific month in a specific year is redistributed uniformly on the wet days.
The steps three and four imply the correction of the precipitation intensity on wet days. Therefore, in step three, the values of a specific month are normalised by the mean of this month: where P** is the normalised time series, P* indicates the time series after the adjustments up to the second step and j is the year (1 . . . 30). The redistribution of precipitation on dry days as well as the normalisation is also done for P obs .
In step four of the extended approach a nonlinear transfer function is derived from the normalised values originating from both P obs and P hist,orig : where P*** hist is the rank ordered and normalised data originating from P hist,orig ; a is the offset and b is the slope of the linear part of the function, τ is a parameter of the exponential part. All three parameters have to be fitted. Finally, the correction is completed by multiplying the mean of a specific month with the value from the transfer function for a specific day: The above presented procedure of the ISI-MIP method is limited by the requirements that (1) there are not too few rain days of a month in the whole period and (2) it is possible to identify a parameter set of a, b and τ of Equation (14) by the fitting algorithm. If one of these requirements is not fulfilled, simplifications of the method are possible (see [31]).

Delta Change
The delta change method, [32] differs from the previously described bias-correction methods. It is not P hist,orig that is adjusted for a better representation of the dynamics of P obs . Instead, a climate signal, obtained from P hist,orig and P fut,orig is used to modify P obs in order to create a synthetic climatological time series for the future, P fut,bc . The climate signal is calculated on the basis of monthly mean values: By definition of the method, the 'corrected' data P hist,bc equals X obs and the climate change signals on a monthly basis do not differ from those of the original climate scenario data. Because of the comparatively low performance in terms of the climate change signal for high precipitation events (see Section 4.3), the delta change approach was not considered for the simulation with the hydrological model.

Evaluation of Bias-Correction Methods
The performances of the different BC methods were evaluated in terms of two objectives: (1) the goodness of fit between Q hist,bc and Q ref and (2) the maintenance of climate change signals for precipitation. For the historical period the different methods were evaluated in terms of the respective ability to reproduce discharge parameters obtained by Q ref . The two parameters MMQ D (mean monthly deviation of discharge) and MHQ D (deviation of MHQ) were evaluated. Within the historical period, there is no deviation between X obs and X hist,bc for the delta change method which implies a perfect agreement between Q ref and Q hist,bc . Therefore, this method is not considered for the evaluation of discharge.
Evaluation of BC regarding discharge parameters is not feasible for future projections as there is no 'correct' climate change signal for discharge on which the evaluation could be based on. For the evaluation regarding future projections, the ability to maintain the precipitation change signal was considered. This is done under the assumption that the climate signal of the bias-corrected climate scenario data should be as close as possible to the signal of the original climate scenario data. Hempel, S. et al. [31] argued that the projected trend, which is a result of the climate models sensitivity, should not be affected by the correction. Cannon, A.J. et al. [33] pointed out that the maintenance of the change signal ensures maintaining the physical relationships of climate variable, e.g., the relationship between projected changes in precipitation and temperature as indicated by the Clausius-Clapeyron equation. It indicates that there is an increase of ca. 7% of column water vapor per K of increase in temperature [34]. While the correction methods linear scaling, delta change and the ISI-MIP approach preserve the relative monthly trend, trend-preserving of extremes is not given for any BC method.
For the evaluation of the change signals for high precipitation values, the parameter P95d3 was chosen as a high correlation between the change signal of this parameter and the change signal for MHQ was identified. The mean of the correlation coefficient for all BC methods, both RCP scenarios and all climate scenarios amounts to 0.77. Furthermore, the maintenance of the monthly change signal is also evaluated. Regarding this criterion, only the methods LOCI, distribution mapping and power transformation were taken into account, as the others, by definition, do not change the monthly signal. Hence the two criteria for the maintenance of the precipitation change signal are the deviation of the P95d3 change signal, P95d3 SigD , and the deviation of the monthly signal, MMP SigD .
The performance of BC methods was calculated for each criterion of evaluation by applying a similar approach as [35] used to evaluate the performance of climate scenarios. The method compares the RMSE of a single criterion for one BC method with the median of all methods (Equation (18)). In terms of the signal for P953d, for example, the criterion is the deviation between the change signal of the original climate scenario data and the signal of the bias-corrected data. The performance for discharge parameters is calculated based on the agreement between Q ref and the Q hist,bc .
where Y is one of the criteria MMQ D , MHQ D , MMP SigD or P95d3 SigD and med is the median. The parameters MMQ D , MHQ D , MMP SigD or P95d3 SigD each comprise 17 values, one for every climate scenario. The computation of these parameters is given in Equations (19)- (22).
abs Sig P,bc (r, g, i) − Sig P,orig (r, g, i) Sig P,orig (r, g, i) P95d3 SigD = 1 2 * nr r=RCP45 1 ng * ng g=1 abs Sig P95d3,bc (r, g) − Sig P95d3,orig (r, g) Sig P95d3,orig (r, g) with Sig P (i) = µ P f ut,bc (i) − µ P hist,bc (i) µ P hist,bc (i) * 100 (23) and Sig P95d3 = P95d3 f ut,bc − P95d3 hist,bc P95d3 hist,bc * 100 (24) where s is the total area of all catchments, c is the catchment (1 . . . 6), s c is the catchment size, r is the RCP scenario (RCP4.5 or RCP8.5), g is the grid cell (1 . . . 93), abs is the absolute value, Sig P and Sig P95d3 are the change signals for mean precipitation and P95d3. The BC method with the lowest performance value for MMQ D or MHQ D is the one that represents best the monthly discharge dynamics or the MHQ simulated with X obs , respectively. The method with the lowest value for performance value for MMP SigD or P95d3 SigD preserves best the climate signal that is projected by the non-bias corrected climate scenario data. Values < 1 indicate that the method performs better than the median of all methods. A value of 0.9, for example, means that the performance of the method is 10% better than the median of all methods.

Climatological and Hydrological Bias in the Historical Period
The comparison of the historical, non-bias-corrected climate scenario data X hist,orig and the observational data X obs for the period 1971-2000 as well as the deviations in corresponding hydrological data Q ref and Q hist,orig, are displayed in Table 2. Table 2. Deviations of climate variables and simulated discharge parameters between climate scenario data and observational data for the period 1971-2000. Climate variables: temperature (T), precipitation (P), 95% percentile of precipitation/3-day precipitation sum (P95/P95d3), number of wet days >1 mm (WD1), wind Speed (Wind), global radiation (Rad) and relative humidity (Hum); hydrological parameters: mean discharge (MQ) and mean annual peak flow (MHQ).

Scenario
∆T The degree of agreement between X obs and X hist,orig as well as between Q ref and Q hist,orig varies considerably among the different climate scenarios. The mean of all climate scenarios shows a cold bias for mean temperature. The mean precipitation and the amount of wet days as well as global radiation are mostly overestimated. However, high precipitation parameters show a good agreement with observations. Wind speed is massively overestimated by all scenarios (approx. +100%). The simulated discharge parameters show a positive bias of 52% for MQ and 43% for MHQ on average. The best agreement for the discharge parameters shows the climate scenario S13, the highest deviations (>100% for MQ or MHQ) are simulated with S3, S5 and S11.
The long-term monthly mean values of the climate variables X obs and X hist,orig and discharge Q ref and Q hist,orig are displayed in Figure 2 for the catchment Wunstorf. Generally, the seasonality of simulated and observed climate variables shows a good agreement, whereby the absolute values often show large deviations. Discharge is overestimated by most scenarios for most of the year. S10 shows high positive deviations of up to 300% in winter. The analyses of monthly data for the other five catchments led to similar results. Overestimations of discharge in such magnitudes may indicate hydrological characteristics of a catchment that are very different from realistic conditions.

Sensitivity of Simulated Discharge against Biases in CORDEX Data
The sensitivity analysis was carried out for all climatological variables that were used as input for the hydrological model. The results are shown in Figure 3. The influence of the climate variables on the simulation of the hydrological parameters MQ and MHQ is similar. The applied alterations of precipitation have the greatest influence on the simulated discharge. However, the other variables have an influence on the hydrological simulations that should not be neglected. The deviation of discharge when changing all variables simultaneously is comparable to that of changing precipitation only. First of all, this can be explained by the by far greatest influence of precipitation. But secondly, the deviations of the other climate variables may have opposing effects on discharge and, therefore, they might neutralize each other to a certain extent. The overestimation of radiation and wind speed, for example, leads to an underestimation of discharge, while the overestimation of humidity leads to more discharge.

Sensitivity of Simulated Discharge against Biases in CORDEX Data
The sensitivity analysis was carried out for all climatological variables that were used as input for the hydrological model. The results are shown in Figure 3. The influence of the climate variables on the simulation of the hydrological parameters MQ and MHQ is similar. The applied alterations of precipitation have the greatest influence on the simulated discharge. However, the other variables have an influence on the hydrological simulations that should not be neglected. The deviation of discharge when changing all variables simultaneously is comparable to that of changing precipitation only. First of all, this can be explained by the by far greatest influence of precipitation. But secondly, the deviations of the other climate variables may have opposing effects on discharge and, therefore, they might neutralize each other to a certain extent. The overestimation of radiation and wind speed, for example, leads to an underestimation of discharge, while the overestimation of humidity leads to more discharge.

Effects of Bias Correction on High Flow Simulation during the Historical Period
The simulated high flow parameters obtained by Qref and Qhist,bc were compared for all six catchments. The deviations for MHQ and HQ100 are displayed in Figure 4.
The agreement between Qhist,bc and Qref for all BC methods is much better than without the application of BC (Section 4.1). The deviations show similar patterns for all catchments. The median of deviations for linear scaling and the ISI-MIP approach shows mostly an underestimation of high flow values, whereas for LOCI, power transformation and distribution mapping, the deviations are smaller. Further, it can be concluded, that the range among the different climate scenarios is high for all catchments and methods. The range is even larger for HQ100 than for MHQ.

Effects of Bias Correction on High Flow Simulation during the Historical Period
The simulated high flow parameters obtained by Q ref and Q hist,bc were compared for all six catchments. The deviations for MHQ and HQ100 are displayed in Figure 4.

Effects of Bias Correction on High Flow Simulation during the Historical Period
The simulated high flow parameters obtained by Qref and Qhist,bc were compared for all six catchments. The deviations for MHQ and HQ100 are displayed in Figure 4.
The agreement between Qhist,bc and Qref for all BC methods is much better than without the application of BC (Section 4.1). The deviations show similar patterns for all catchments. The median of deviations for linear scaling and the ISI-MIP approach shows mostly an underestimation of high flow values, whereas for LOCI, power transformation and distribution mapping, the deviations are smaller. Further, it can be concluded, that the range among the different climate scenarios is high for all catchments and methods. The range is even larger for HQ100 than for MHQ.  The agreement between Q hist,bc and Q ref for all BC methods is much better than without the application of BC (Section 4.1). The deviations show similar patterns for all catchments. The median of deviations for linear scaling and the ISI-MIP approach shows mostly an underestimation of high flow values, whereas for LOCI, power transformation and distribution mapping, the deviations are smaller. Further, it can be concluded, that the range among the different climate scenarios is high for all catchments and methods. The range is even larger for HQ100 than for MHQ. Figure 5 shows the climate change signals for P95d3 and for MHQ for bias-corrected and non-corrected climate scenario data. In terms of P95d3, an increase was simulated by all climate scenarios and both RCP scenarios. These climate change signals are generally higher for RCP8.5 than for RCP4.5 and mostly they are higher for the bias-corrected data. The deviations between the signals of corrected and non-corrected data are higher for the BC methods power transformation, distribution mapping or LOCI than for linear scaling and the ISI-MIP approach.  Figure 5 shows the climate change signals for P95d3 and for MHQ for bias-corrected and non-corrected climate scenario data. In terms of P95d3, an increase was simulated by all climate scenarios and both RCP scenarios. These climate change signals are generally higher for RCP8.5 than for RCP4.5 and mostly they are higher for the bias-corrected data. The deviations between the signals of corrected and non-corrected data are higher for the BC methods power transformation, distribution mapping or LOCI than for linear scaling and the ISI-MIP approach.

Impact of Bias Correction on the Future Projections of Precipitation and Discharge
In terms of MHQ, the climate change signals projected by non-corrected data show on average only small changes for RCP4.5 and an increase of ca. 10% for RCP8.5. BC mostly leads to higher change signals in case the signal of non-corrected data was positive. In case of a negative signal obtained by non-corrected data, the BC led to either a smaller negative or even positive signal. Typically, for MHQ the deviations of the signals among the different BC methods are smaller than the difference between the signals of one BC method and non-corrected data for one climate scenario and one RCP scenario.
The deviations in climate change signals are higher for MHQ than for P95d3.  In terms of MHQ, the climate change signals projected by non-corrected data show on average only small changes for RCP4.5 and an increase of ca. 10% for RCP8.5. BC mostly leads to higher change signals in case the signal of non-corrected data was positive. In case of a negative signal obtained by non-corrected data, the BC led to either a smaller negative or even positive signal. Typically, for MHQ the deviations of the signals among the different BC methods are smaller than the difference between the signals of one BC method and non-corrected data for one climate scenario and one RCP scenario.
The deviations in climate change signals are higher for MHQ than for P95d3.

Evaluation of Climate Scenarios
The results for the evaluation of BC methods for both objectives, the agreement of Q hist,bc with Q ref and the maintenance of the climate change signal obtained by non-corrected data P orig by the corrected data P bc are shown in Figure 6.
The results for the evaluation of BC methods for both objectives, the agreement of Qhist,bc with Qref and the maintenance of the climate change signal obtained by non-corrected data Porig by the corrected data Pbc are shown in Figure 6.
The methods do not show great differences regarding monthly mean discharge, MMQD, but there are differences regarding the performance for the simulated MHQ. Power transformation, distribution mapping and LOCI perform better than linear scaling and the ISI-MIP approach.
For the evaluation of MMPSigD, the methods delta change, linear scaling and ISI-MIP are not included, because of the perfect agreement with the signal obtained from non-bias-corrected data (see Section 3.5). LOCI performs better than distribution mapping and power transformation. For P95d3SidD, the performance of the ISI-MIP approach and linear scaling is significantly better than the other methods, whereas delta change shows a much lower performance.

Discussion and Conclusions
The comparison of observed and simulated climate variables showed a significant bias in climate scenario data. When applying this data to hydrological models, significant impacts on hydrological parameters are the consequence. Substantial differences in the discharge volume indicate fundamentally changed hydrological processes, e.g., an altered ratio of discharge and infiltration due to much higher water saturation of the soil or a significantly higher groundwater recharge. Therefore, applying bias-correction on climate model data prior to the use as input data for hydrological models seems necessary as it leads to simulations of hydrological processes that correspond much better to observed conditions.
The sensitivity analysis showed that deviation in precipitation of the climate scenario data are by far the most critical factor with respect to the simulated discharge. Therefore, the variable precipitation should be bias-corrected for a plausible simulation of discharge. Besides, temperature The methods do not show great differences regarding monthly mean discharge, MMQ D , but there are differences regarding the performance for the simulated MHQ. Power transformation, distribution mapping and LOCI perform better than linear scaling and the ISI-MIP approach.
For the evaluation of MMP SigD , the methods delta change, linear scaling and ISI-MIP are not included, because of the perfect agreement with the signal obtained from non-bias-corrected data (see Section 3.5). LOCI performs better than distribution mapping and power transformation. For P95d3 SidD , the performance of the ISI-MIP approach and linear scaling is significantly better than the other methods, whereas delta change shows a much lower performance.

Discussion and Conclusions
The comparison of observed and simulated climate variables showed a significant bias in climate scenario data. When applying this data to hydrological models, significant impacts on hydrological parameters are the consequence. Substantial differences in the discharge volume indicate fundamentally changed hydrological processes, e.g., an altered ratio of discharge and infiltration due to much higher water saturation of the soil or a significantly higher groundwater recharge. Therefore, applying bias-correction on climate model data prior to the use as input data for hydrological models seems necessary as it leads to simulations of hydrological processes that correspond much better to observed conditions. The sensitivity analysis showed that deviation in precipitation of the climate scenario data are by far the most critical factor with respect to the simulated discharge. Therefore, the variable precipitation should be bias-corrected for a plausible simulation of discharge. Besides, temperature is often bias-corrected in hydrological climate impact studies. Our findings moreover suggest investigating biases in all climate variables that are used as input for the hydrological model as they can cause substantial deviations in discharge (mean flows and high flows) as well.
Higher climate change signals for high precipitation parameters were calculated for bias-corrected climate scenario data than for non-corrected data. The deviations for the climate change signals for high flow parameters are even larger than those for high precipitation parameters when using the corrected/non-corrected climate scenario data as input for the hydrological model. In most cases when the signal for the high flow parameter MHQ obtained by non-corrected climate scenario data was positive, the application of BC led to higher signals. In case the signal with non-corrected climate scenario data was negative, the signal with bias-corrected data was lower or even positive.
It could be argued that in climate change studies the most relevant result is the climate change signal and not absolute values, e.g., of discharge. However, it should be considered that signals that are based on simulations are supposed to represent the process dynamics of a certain region. If these processes significantly differ from realistic conditions, because of high deviations in input data, the magnitude of climate change signals obtained by these simulations have to be critically questioned. Therefore, BC should be applied when climate scenario data significantly differ from observations.
The evaluation results also showed that no BC method is clearly better than the others in terms of all evaluated criteria. While power transformation, LOCI and distribution mapping led to a more sophisticated representation of high flows during the reference period, linear scaling and ISI-MIP were better capable of maintaining the projected changes of the non-corrected data in terms of monthly precipitation and high precipitation values.
An assumption of BC in general is that the deviations found for the reference period remain constant over time. This can be doubted, as it was shown, for example by [14], that BC procedure may be sensitive to the choice of the reference period. Lafon, T. et al. [14] found less complex methods to be more robust and less vulnerable to over-tuning. From this study covering the historic period it cannot be concluded that the less complex methods lead to more realistic corrections of the simulations of future climate as variations between slightly modified reference periods can be expected to be smaller than the variations between future climate, especially at the end of the 21 st century, and the reference period. Due to the projected changes in climate variables, the reliability of the calibrated BC procedure for all BC methods may decline in the course of the century.
When investigating climate change impacts on hydrological parameters on a regional scale, the additional source of uncertainty of BC to the entire model chain uncertainty should not be neglected. In case only one BC method can be applied, the performance of the selected method related to the study objective should be evaluated for the study area in advance.