Abstract. General circulation models (GCMs) are the primary tools to evaluate the possible impacts of climate change; however, their results are coarse in temporal and spatial dimensions. In addition, they often show systematic biases compared to observations. Downscaling and bias correction of climate model outputs is thus required for local applications. Besides the computationally intensive strategy of dynamical downscaling, statistical downscaling offers a relatively straightforward solution by establishing relationships between small and large scale variables. This study compares four statistical downscaling methods (SDMs) of bias correction (BC), change factor of mean (CFM), quantile perturbation (QP) and event based weather generator (EBWG) to assess climate change impact on drought by the end of the 21st century (2071–2100) relative to a baseline period of 1971–2000. A set of drought related aspects is analysed: dry day frequency, dry spell duration and total precipitation. The downscaling is applied to a 14-member ensemble of CMIP6 GCMs, each powered by four future scenarios of SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5. A 25-member ensemble of CanESM5 GCM is also used to assess the significance of the climate change signals in comparison to the internal variability of the climate. While all methods show a good agreement on downscaling total precipitation, the CFM method fails to downscale dry day frequency well. The QP method outperforms the others in downscaling dry spells. Using this method, dry day frequency is projected to increase significantly in the summer months, with relative changes of up to 20.4 % in the worst-case climate change scenario. At the same time, total precipitation is projected to decrease significantly by up to 41.9 % in these months. Lastly, extreme dry spells are projected to increase in length by up to 7.4 %.