Abstract

We performed bias correction in future climate change scenarios to provide better accuracy of models through adaptation to future climate change. The proposed combination of the change factor (CF) and quantile mapping (QM) methods combines the individual advantages of both methods for adjusting the bias in global circulation models (GCMs) and regional circulation models (RCMs). We selected a study site in Songwol-dong, Seoul, Republic of Korea, to test and assess our proposed method. Our results show that the combined CF + QM method delivers better performance in terms of correcting the bias in GCMs/RCMs than when both methods are applied individually. In particular, our proposed method considerably improved the bias-corrected precipitation by capturing both the high peaks and amounts of precipitation as compared to that from the CF-only and QM-only methods. Thus, our proposed method can provide high-accuracy bias-corrected precipitation data, which could prove to be highly useful in interdisciplinary studies across the world.

1. Introduction

With the increase in industrialization around the world, the quality of human life has improved in various aspects through advancement in scientific technology, professionalism in occupations, level of consumption, and so forth. However, these improvements have caused an excessive use of fossil fuels (coal, petroleum, natural gas, etc.) [1] and an increase in global temperatures triggered by the increase in greenhouse gases (carbon dioxide, nitrogen dioxide, etc.) in the troposphere [2].

In general, human activity in agriculture, industry, and commerce influences global ecosystems through factors such as aerosol effects, land use change, and deforestation. This can lead to negative impacts such as increase in air temperature and changes in precipitation patterns [3], which mainly influence climate change globally [4]. According to the 4th Assessment Report (AR4), Intergovernmental Panel on Climate Change (IPCC) [5], the number of days with precipitation and flood occurrence frequencies has steadily and gradually increased due to unusual climate changes indicating that natural disasters will occur more frequently and intensely in the future than before [5].

Republic of Korea has mountainous terrain across approximately 70% of its area and more than 70% of the country’s precipitation occurs during the summer period (July to September) [6]. Thus, significant amount of water from rainfall is usually discharged into streams as direct runoff and the topographical characteristics of the land aid this discharge. Republic of Korea has suffered from both frequent floods aided by the direct runoff in the summer period and droughts caused by water scarcity due to unbalanced rainfall during the crop growing seasons. These repetitive natural disasters have caused great economic losses and uncertainties in sustainable water resources management plans across the country every year [7].

For these reasons, the Korean government has been striving to reduce the potential risks (i.e., drought, flood, etc.) of effects from climate change and recently suggested the Four Major Rivers Restoration Project as one solution [8]. In order to verify the effects of this project at the field level, various factors, such as the dynamic interactions between meteorological and hydrological processes across the land surface, rainfall-runoff processes, seasonal weather predictions, ecosystem dynamics and alteration, and regional characteristics, need to be studied.

In Republic of Korea, Korea Meteorological Administration (KMA) provides weather predictions based on various climate change scenarios [3]. Kim et al. [9] and Lee et al. [10] assessed flood vulnerability at watershed-scales using KMA-based climate predictions. Kwon et al. [11] explored the boundary changes in a subtropical climate region (Republic of Korea) using the A1B emission scenario. Park et al. [12] evaluated the impacts of climate change on river water quality at the Chungju-dam watershed using a soil water assessment tool (SWAT) [13] model. These studies show that climate change scenarios have been used to predict and assess long-term, sustainable development and management of water resources.

The use of weather prediction has been steadily increasing, but global circulation models (GCMs)/regional scale models (RCMs) provided by the KMA still have limitations for application at local/regional scales due to the discrepancy in scale between GCMs/RCMs and field-scale (or modeling) resolutions [14]. These drawbacks may cause large uncertainties when planning for sustainable water resources management for the future. For these reasons, bias in the GCMs/RCMs needs to be corrected by comparing the weather history of local regions.

Currently, several bias-correction schemes, such as change factor (CF) [15], quantile mapping (QM) [16], and multiple linear regression [17], which have been developed and improved since a few decades ago, exist. CF is a simple downscaling method that uses the average values of observations and predictions. This method is conducted by simply scaling the average change factor to each day [18]. Due to its simplicity, this concept has been used in many climate related bias-correction applications [1923]. The QM scheme corrects GCMs/RCMs based on the cumulative distribution function (CDF) with a statistically good match [24]. The QM method has been widely employed to correct the bias of GCMs/RCMs (e.g., [18, 2527]). Both CF and QM have their main advantage in that the statistical approach to bias corrections has better computational efficiency than other physical-based approaches [18, 23]. However, bias correction using the CF method changes only the average, maxima, and minima of the climatic index in scenarios, while all the other properties, such as the number of wet/dry days and the variance of temperature, remain unchanged [28, 29]. QM has limitations in capturing extreme values beyond the range of observations (e.g., [24, 30]). Since CF and QM both have their advantages/disadvantages for bias correction, a linked CF and QM approach that can overcome both of their disadvantages may perform better in correcting GCM/RCM data to get them closer to historical observations.

Thus, we explored an enhanced bias-correction approach for improving the applicability of GCMs and RCMs to local regions. The objectives of this study were twofold: to propose an enhanced bias-correction method linking CF and QM and to assess the performance of our approach at field level. Our proposed approach can be substantially more useful for providing realistic bias-corrected GCM/RCM projections that are closer to the weather history and for creation and implementation of sustainable long-term water resources management plans for the future.

2. Experimental Section

2.1. Literature Review
2.1.1. Bias-Correction Methods

The change factor (CF) method has been widely used for bias corrections because of its simplified approach. Equation (1) describes how the CF corrects the bias between observations and predictions: where is the transformed future precipitation, is the predicted future precipitation, is the average of observed precipitation, and is the average of historical A1B precipitation.

The CF approach calculates the average values of the monthly/yearly observations and predictions. Then, the CF estimates the bias-correction coefficients indicating the differences between the observations and predictions and shifts the predicted results to the observation data based on the calculated coefficients (Figure 1).

The quantile mapping (QM) method, proposed by Panofsy and Brire [16], minimizes the differences between the observed/predicted data based on empirical probability distributions. Figure 2 shows how the QM shifts the predicted data to the observations with the same nonexceedance probability. Here, we used the one-to-one mapping scheme with the order statistics of the observed simulated data, instead of fitting the mathematical model based on the cumulative distribution functions: where is the cumulative distribution function of the observed daily data for day , is the cumulative distribution function of the simulated data from historical simulations, and and are the simulated and transformed (bias-corrected) data, respectively, for day (2).

The transformed predictions have the same probability distribution with the observations, but QM has a limitation in generating distributions on a monthly basis due to the small amount of data points available [3134].

2.1.2. Climate Change Scenarios

The Intergovernmental Panel on Climate Change (IPCC) provides past and future weather scenarios based on various meteorological, socioeconomic, and environmental data for assessing the impacts of climate changes [5]. Based on these scenarios from the IPCC, the KMA has generated global circulation model (GCM) and regional climate model (RCM) outputs. RCM predictions have relatively finer-scale resolutions (27 km × 27 km) compared to GCMs (i.e., 3.75° × 3.75°, approximately ~400 km). The Special Report on Emissions Scenarios (SRES, the 3rd Assessment Report) by the IPCC provides four different storylines (A1, A2, B1, and B2) based on different demographic, social, economic, technological, and environmental developments. The A1B (balanced emphasis on all energy sources) scenario is one of the most realistic future climate scenarios provided by the KMA.

The KMA used the A1B scenario as the input data to a GCM (ECHO-G) that included the Coupled Climate model comprised of the Atmosphere model (ECHAM4), the Thawing-Ocean model (HOPE-G), and the Atmosphere Ocean Coupled model (OASIS) for generating GCM scenarios. In addition, high-resolution climate change scenarios in Republic of Korea were provided by the KMA using the MM5 (5th Generation Mesoscale Model) [35] based regional climate model (Figure 3). The KMA climate scenarios using GCM (ECHO-G)/RCM (MM5) reflect the Korean climate features well [36]. For this reason, KMA recommends the use and integration of their scenarios as national standard science information in policies, projects, and studies related to climate change in Republic of Korea. In this regard, numerous studies have used KMA scenarios for climate change applications in Republic of Korea [3744]. Thus, we selected the KMA climate change scenarios using the A1B scenario for assessing our proposed approach.

2.2. Materials and Method
2.2.1. Enhanced Bias-Correction Approach Linking the CF and QM Methods

Both the CF and QM methods have their advantages and disadvantages in the bias-correction process. Bias-corrected (transformed) precipitation of the future will vary based on the selection of either CF or QM. In general, rainfall quantities/frequencies considerably influence the occurrence of floods/droughts. As mentioned in Section 1, heavy rainfall events/amounts, such as the 941.5 mm rainfall from July to September in 2004 [23], have been generated during the summer period indicating that small uncertainties (especially for rainfall amounts) in bias-correction methods and inherently ill-posed GCM data may lead to disastrous damages across the country. Thus, these uncertainties need to be reduced for efficient water resources management by policy/decision makers. Although the linked CF and QM method can contribute to a reduction in uncertainties in the bias-correction methods (CF or QM), the inherent uncertainty caused by drawbacks within the GCM/RCM data provided by the IPCC/KMA is unavoidable. To date, correcting daily GCM/RCM data is still a challenge globally. Therefore, we corrected the bias of monthly (precipitation) predictions using CF and then, using the QM scheme, shifted the CF-based bias-corrected predictions to the historical observations as shown in Figure 4.

2.2.2. Application of the CF, QM, and CF + QM Bias Correction

In this study, the Songwol-dong site in Seoul, Republic of Korea, shown in Figure 5, was selected to assess the bias correction of the A1B precipitation scenario using three methods (CF, QM, and CF + QM) against the observed precipitation data during a baseline period. The bias-correction coefficients (change factors for the CF method and a CDF for the QM method) obtained from the bias correction during the baseline period were applied to the A1B precipitation scenario. Twelve change factors, one for each month, were estimated, and a CDF was estimated based on the relationship between the observed precipitation data and the A1B precipitation scenario during a baseline period.

For evaluation of the bias correction by the three methods, we divided the observation period (1983–2012) into two groups: the baseline period (1983–1997) and the validation period (1998–2012). We obtained the bias-correction coefficients of CF and QM between the A1B precipitation scenario (1983–1997) and the observed precipitation data through bias correction during a baseline period (1983–1997). Then, the bias-correction coefficients were applied to the A1B precipitation scenario (1998–2012) for comparison with the observed precipitation (1998–2012) for validation.

The weather in Republic of Korea shows seasonal changes due to the monsoon climate. For this reason, we classified the monthly averaged precipitation into four categories (January–March, April–June, July–September, and October–December) and assessed the bias-corrected data using the coefficient of determination (), Nash-Sutcliffe efficiency (NSE), and root mean square error (RMSE). Then, we estimated the bias-correction coefficients of the A1B precipitation scenario relative to the observed precipitation data using another baseline period (1983 to 2012) for bias correction of future climate change scenario. We corrected the A1B precipitation projections for the future (2016–2100) based on the derived bias-correction coefficients from 1983 to 2012, and their characteristics/statistics were analyzed.

The measures the degree of collinearity between observations and simulations: where is the th observation, is the th simulation, is the mean of the observations, is the mean of the simulations, and is the total number of observations.

The NSE is a normalized statistic that gives the relative magnitude of the residual variance compared to the observed variance [45], as shown in where an NSE value of 1 indicates perfect agreement between observations and simulations.

The RMSE is one of the commonly used error index statistics for observed and simulated data [46], as shown in where an RMSE of zero indicates perfect agreement between observations and simulations.

3. Results and Discussion

3.1. Evaluation of Bias Correction Using the Three Methods

In Republic of Korea, 60–70% of the yearly precipitation amount [47] is usually generated during the summer period (July to September). To be prepared for natural disasters (drought, flood, etc.) and ensuring sustainability of water resources, efficient water management plans need to be formulated for the heavy rainfall period. The selection of precise bias-correction methods can help in obtaining reliable projected precipitation changes across the Korean peninsula for the summer period in the future, which can help in mitigating potential risks from natural disasters.

Table 1 presents the statistical analysis of the bias-corrected precipitation using the three methods. The bias-corrected precipitation results of the combined CF + QM method were considerably improved for baseline period (1983–1997) and validation period (1998–2012) as compared to those of the CF-only or QM-only methods; however, there were no significant differences between the results of the CF, QM, and CF + QM methods for the dry seasons (January to March, April to June, and October to December). The CF method overestimated the bias-corrected A1B precipitation scenarios as compared to the baseline precipitation data while the bias-corrected precipitation by the QM scheme was underestimated for baseline and validation periods. The fact that the CF + QM method performed better at correcting the A1B precipitation scenarios relative to the baseline precipitation data for the bias correction and validation of the summer period (July to September) supports the robustness of our proposed bias-correction method (combined CF + QM) for application to heavy rainfall seasons. Even though uncertainties exist in the bias correction/validation periods, the bias-corrected results of the CF + QM method are closer to the observations than those of the CF-only or QM-only methods, as shown in Figures 6 and 7.

3.2. Bias Corrections of the A1B Precipitation Projections for the Future

Since the CF + QM method showed good results in validation, we corrected the A1B precipitation scenario from 1983 to 2100 at the identical region. Note that we conducted bias correction using the three methods for the comparison using the baseline precipitation forcings (1983 to 2012), and then the bias-correction coefficients were applied to the precipitation projections from 2016 to 2100.

Figure 8 shows the comparison of the observations and transformed (bias-corrected) predictions based on the three methods. Comparison of the results of the CF, QM, and CF + QM methods shows that the monthly distribution of the transformed daily precipitation (5th/25th/50th/75th/95th percentiles and mean) from the QM is considerably different from the observations. These results imply that the QM is limited in correcting monsoon climatic predictions because intensive rainfall frequencies/amounts during heavy rainfall periods might cause over-/underestimations in correcting bias.

The results of the CF and CF + QM methods have similar distributions as shown in Figure 8. The CF method performed better in correcting the bias of the monthly precipitation data, but the CF had a limitation in that it was highly dependent on the distribution of the predicted precipitation data.

Table 2 shows the comparison of the monthly maximum precipitation data from the observations and transformed data using the CF and CF + QM methods. With respect to the monthly observation of 1,131.0 mm (the monthly maximum precipitation) in July, the CF method generated an amount of 1,540.3 mm (36.2% difference), while the CF + QM method generated 1,319.4 mm, a 16.7% lower deviation than the CF. However, the monthly averages of the observed and transformed precipitation data in July were identical at 425 mm (Figure 8). In August, the maximum precipitation using the CF method was 1,535.6 mm (24.1% difference as compared with the observed maximum precipitation of 1,237.8 mm). On the contrary, the CF + QM method generated a monthly maximum precipitation of 1,271.0 mm, a difference of only 2.7% as compared with the observation data (Table 2).

Table 3 shows the number of days with the statistics (maximum and minimum, sum, average, and standard deviation) of the weather history and the bias-corrected/uncorrected predictions. The daily maximum transformed precipitation based on the CF method is 823.1 mm (July 2004), while the maximum daily precipitation of the observation data is 332.8 mm. The total precipitation in 2004 was 1020 mm, which is approximately 20% less than the annual average precipitation in Republic of Korea. The use of a correction factor for long periods and in dry years might overestimate the daily precipitation in bias correction. These results demonstrate that the CF method can yield inaccuracies when used for dry years because correction factors may have their limitations in reflecting the variations in yearly precipitation data. The CF method corrected well the monthly/yearly precipitation amounts as well as the number of predicted precipitation occurrences when compared to the baseline precipitation data. However, this method showed high uncertainties in matching the daily outputs as shown in the statistics (maximum, average, and standard deviation, see Table 2). However, the proposed CF + QM method provided better bias-corrected predictions than the others because it reduced the limitations of both the CF and QM methods.

The bias-uncorrected predictions generated a total precipitation amount of 68,528 mm for the future period (2016–2100), but the bias-corrected results of the CF, QM, and CF + QM methods were 141,613 mm, 146,637 mm, and 140,665 mm, respectively. Although the total quantities of the bias-corrected precipitation during the entire period are similar, the temporal distributions of the transformed future precipitation data showed large differences as shown in Figures 9, 10, and 11. The yearly precipitations have extremely high frequencies from 2050 to 2070 for CF and 2070 to 2090 for QM. However, the degree of high frequencies corrected by the CF + QM method was relatively lower than those by the CF-only and QM-only methods, as shown in the bias-corrected results of bias correction and validation. This demonstrates that the CF + QM method can improve the limitations of the individual CF and QM methods. Figure 12 shows the ranges of the yearly precipitation predictions transformed by CF, QM, and CF + QM with five-year intervals. Figure 12 shows that the bias-corrected predictions based on the CF + QM are usually located near the center, between those of the CF and QM. Thus, using a single bias-correction method (CF or QM) has relatively higher uncertainties in providing accurate future precipitation projections as compared to the combined CF + QM method. Furthermore, errors incorporated within bias-correction methods might be propagated into future projections that may be used by decision makers. Therefore, appropriate bias-correction methods that reduce errors in future precipitation projections should be used. Our results show that the CF + QM method, which can improve the limitations of the CF and QM methods, can be appropriate for generating transformed precipitation data using various climate change scenarios.

4. Conclusions

Preparations for and adaptations to climate change are planned or ongoing throughout the world. This study suggests a combined bias-correction (CF + QM) method that combines the existing CF and QM methods to help prepare for the negative effects of climate change.

To assess the performance of the proposed method for bias correction for long time periods, bias corrections were conducted for the study period (1983–2012) for the Songwol-dong site in Republic of Korea. The first 15 years (1983–1997) were used as a baseline period for estimating bias-correction coefficients, and we validated three methods for the latter 15 years (1998–2012). We classified the study period into four categories (January–March, April–June, July–September, and October–December) for assessing the characteristics of seasonality when correcting projected precipitation bias. Then, we tested/compared the performance of the CF-only, QM-only, and CF + QM methods for the bias correction and validation periods. Although the CF + QM method corrected the projected precipitation bias to the baseline precipitation data slightly better than the others, there exist no/little significant differences for the dry periods (January–March, April–June, and October–December). However, the CF + QM method (R2: 0.990/0.990 and NSE: 0.920/0.916 for bias correction/validation) considerably improved the bias-corrected precipitation during the heavy rainfall season (July–September) compared to the CF-only (R2: 0.958/0.859 and NSE: 0.887/0.703 for bias correction/validation) and QM-only (R2: 0.916/0.852 and NSE: 0.860/0.700 for bias correction/validation) methods.

After comparing the three methods for bias correction, we derived the bias-correction coefficients using the CF-only, QM-only, and CF + QM methods for the entire period (1983–2012). Then, we corrected the A1B precipitation scenario based on the estimated coefficients for the future period (2016–2100). The CF + QM method gives better bias-corrected predictions with respect to statistics. In particular, the suggested CF + QM method is appropriate for monsoon climatic conditions when torrential precipitation occurs during the wet season. Furthermore, our proposed approach can contribute to a reduction in weather prediction uncertainties and the establishment of sustainable water resources management plans for mitigating natural disasters caused by future climate changes. We plan to extend the applicability and accessibility of the suggested method to the web, which will contribute considerably to the preparation for and adaptation to uncertain climate changes.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This research was supported by the Geo-Advanced Innovative Action (GAIA) Project (no. RE201402074, Surface Soil Resources Inventory & Integration: SSORII Research Group) in Republic of Korea.