Global climate change has a significant impact on the water cycle and natural ecosystems. The quantification of the impact of climate change on precipitation, temperature, and water resources has attracted substantial attention (IPCC, 2013). Global climate models (GCMs) have been developed and applied increasingly in recent decades (Stouffer et al., 2017;IPCC, 2013;Taylor et al., 2012). However, the low spatial resolution of GCMs is inadequate for high-precision and high-resolution climate information required in regional studies. Therefore, regional climate models (RCMs) have been developed by dynamically downscaling GCMs to obtain climate data with a high temporal and spatial resolution (Giorgi, 2019;Giorgi et al., 2009;Walton et al., 2020;Xu, Han, & Yang, 2019). RCMs have significant value in areas with high altitude, complex terrain, or sparse and uneven rainfall stations. Recently, many studies have been carried out using future climate data obtained from RCMs to drive hydrological models and obtain future watershed hydrological characteristics; these studies provided valuable results (Gorguner et al., 2019;Zhou et al., 2018). Commonly used RCMs include the Weather Research and Forecast model, Regional Climate Model (RegCM), Hadley Centre Global Environmental Model Version 3 Regional Climate Model, and Providing Regional Climates for Impacts Studies; among them, the RegCM is the most widely used model (Gao & Giorgi, 2017;Gao et al., 2018). However, due to different spatial resolution of GCMs, different parameterization schemes of the RegCM, and different land-use/land cover types, substantial differences occur between the RegCM output and field observations, especially in mountainous watersheds with complex topography, such as the Upper Beijiang River Basin (UBRB) in South China (Tong et al., 2017(Tong et al., , 2020. Some studies have found that using the RegCM output in hydrological research without any bias correction amplifies the rainfall-runoff deviation, causing uncertainties (Ozturk et al., 2017;da Silva et al., 2019). Therefore, it is crucial to perform bias correction, which is also referred to as statistical downscaling of the model output, to ensure the model's applicability to the hydrological research.