Assessing Potential Climate Change Impacts and Adaptive Measures on Rice Yields : The Case of Zhejiang Province in China

Increasing temperatures, greater carbon dioxide concentrations, and changes in related climatic variables will continue to affect the growth and yields of agricultural crops. Rice (Oryza sativa L.) is extremely vulnerable to these climatic changes. Therefore, investigating the degree to which climate changes could influence rice yields and what effective adaptive strategies could be taken to mitigate the potential adverse impacts is of vital importance. In this article, the impacts of climate change on rice yields in Zhejiang province, China, were simulated under the Representative Concentration Pathway (RCP) 4.5 and 8.5 scenarios. The impacts of climate change, with and without CO2 fertilization effects, were evaluated and the three most effective adaptive measures were examined. Compared with the yield for the baseline time of 1981–2010, the simulated average yields of all cultivars were inevitably projected to decrease under both RCPs when the CO2 fertilization effects were not considered during the three periods of the 2020s (2011–2040), 2050s (2041–2070), and 2080s (2071–2099), respectively. Declines in rice yields were able to be alleviated when the CO2 fertilization effects were accounted for, but the yields were still lower than those of the baseline. Therefore, the three adaptive measures of advancing planting dates, switching to high-temperature-tolerant cultivars, and breeding new cultivars were simulated. The results indicated that adaptive measures could effectively mitigate the adverse effects of climate change. Although the simulation had uncertainties and limitations, the results provide useful insights into the potential impacts of climate change in Zhejiang province while also proposing adaptive measures.

The recent Intergovernmental Panel on Climate Change (IPCC) special report showed that the global temperature related to human activities has increased by approximately 1.0 • C compared with preindustrial levels, and it is likely to increase to about 1.5 • C between the years 2030 and 2050 [9][10][11].
In tropical and subtropical areas, especially in developing countries, the impacts of climate change on agricultural production will be more apparent because of the tropical climate and the lack of advanced agricultural technical input [12][13][14][15].Rice is a staple food that feeds more than half of the global population and it is sensitive to even small changes in climate, which may cause a significant reduction in yield [16,17].Therefore, investigating the potential impacts of climate change on rice production and developing effective adaptive measures in subtropical areas, especially in developing countries, should be a high priority.
China is the world's largest developing country and it has a long history of rice cultivation.Rice production in China accounts for almost 30% of global rice production [18,19].Therefore, researchers in institutions and at all levels of government in China are closely observing, collecting, and documenting data on the country's crop growth and development [20,21].These high-quality data are collected and recorded by following standardized observation criteria established by agricultural technicians [22,23].
Zhejiang province, China, has a total rice cultivation area of 8.242 × 10 5 km 2 , accounting for 27.12% of all the rice cultivated in agricultural areas of China [24], thus, the rice yield in this region is critical.Furthermore, the entire area lies in a subtropical monsoon climate zone, and research evidence has shown that the advent of rising temperatures and extreme weather events will be more likely to occur here.There is an urgent need to investigate the degree to which these factors will have an impact on rice yields and to identify potential adaptive strategies to minimize their adverse effects.
So far, various methods have been adopted to investigate the underlying mechanism by which climate change affects agricultural yields.Model-based research is preferred because simulation experiments can be performed hundreds or even thousands of times in an effort to select the optimal strategies for guiding field management.The Crop Estimation through Resource and Environment Synthesis (CERES)-Rice model is a process-based real system, designed as a natural surrogate laboratory to simulate crop growth associated with factors such as weather changes, soil properties, and management practices.The model is continuously being refined and modified for simulation to ensure that credible and convincing results are obtained [25,26].The simulation results obtained with these models are the most effective methods of evaluating the potential impacts of climate change on rice [27].The effect of tillage systems has not been considered because no data are available that record real management operations after the harvest of early-mature rice [28].
Several researchers have applied various models to assess the potential impacts of climate change on rice yields in China.Erda et al. [29] simulated the potential climate change impacts on agricultural yields and yield quality with CO 2 fertilization under regional climate change model scenarios, using Providing Regional Climates for Impacts Studies (PRECIS) at the national scale in the CERES-Rice model, version 3.0.Zhang et al. [30] used the CERES-Rice, ORYZA2000, Regional Climate Model (RCM), Beta Model (BM), and Simulation Model For Rice-weather Relations (SIMRIW) to model climate change impacts on rice phenology in northeastern and southwestern China.Their results showed that the CERES-Rice and ORYZA2000 models performed well.Yao et al. [31] simulated the impacts of climate change on rice yields, including consideration of the CO 2 fertilization effects, by using data from eight representative sites in the CERES-Rice model, however, the authors proposed that no related adaptive measures were adopted for raising future rice production.Xiong et al. [32] used the decision support system for agrotechnology transfer (DSSAT), version 4.0, to assess the impacts of climate variability on rice yields, including a consideration of spatial and temporal characteristics, then predicting the possible yield increment in the principal yield regions in China at a grid resolution of 50 × 50 km.Yang et al. [33] simulated the rice yields in China with the CERES-Rice model and found that a 1 • C warmer temperature during the growing season would reduce the rice yield by roughly 4 percentage points.Their innovative findings may provide an understanding of the process and mechanisms of the impacts of climate change.However, previous studies on the impacts of climate change on rice yields have been conducted either with climate scenarios generated from a single general circulation model (GCM) or with a limited number of representative sites on a scaled-up national 0.5 • C × 0.5 • C grid level.Some studies have even assessed the global climate change impacts on agricultural production by using only specific varieties of crops, which represents a misunderstanding of the basic cropping system.Some have assessed only the impacts of climate change without considering the dynamic change in CO 2 concentration.In addition, previous studies have focused mainly on impacts, and very few have analyzed the adaptive strategies needed to deal with the potential negative effects.
In this study, the potential climate change impacts on rice yields, both with and without CO 2 fertilization effects, during three periods, the 2020s (2011-2040), 2050s (2041-2070), and 2080s (2071-2099), were evaluated based on highly resolved management data within the CERES-Rice model, version 4.7, in Zhejiang province, China.Adaptive measures to reverse the negative climate effects and cope with the challenge of climate change were also simulated.In addition, the uncertainties caused by climatic models were quantified by using climate models.The objectives of the study were as follows: (1) Assess the potential climate change impacts on rice yields, with and without the effects of atmospheric CO 2 fertilization; (2) to assess the uncertain impacts of climate change on rice yields by using five GCM models; and (3) to simulate adaptive measures, such as adjusting planting dates, switching to high-temperature-tolerant cultivars, and breeding new cultivars.

Study Area and Observation Sites
Zhejiang province (27 • 12 -31 • 31 N, 118 • 01 -123 • 10 E) belongs to the Yangtze River Basin (Figure 1) in China.It lies on the southeast coast of China and is known as "the land of fish and rice," which means rice is the grain base in China [34,35].About 57 million people inhabit in this region, and the main income here is mainly from agricultural production and booming industries.The entire area is a typical subtropical zone, thus, both the temperature and precipitation here are quite suitable for a double-cropping system, which means a rice cultivar can be divided into early-mature and late-mature.The growth period of early-mature rice is mainly from March to June and the period for late-mature rice is mainly from July to October, after the harvest of early-mature rice.Favorable soil conditions, abundant precipitation, and a scientific cropping system all make Zhejiang a suitable natural location for rice cultivation.Zhejiang is one of the main rice production areas in the world, so the effects of climate change are more obvious here than elsewhere.The rice produced here is more vulnerable to these changes than in other areas, thus, the importance of assessing the impacts of climate change on rice yields in this region cannot be understated.
Experimental records such as rice management practices, weather conditions, fertilization data, rice phenology, rice cultivars, and rice yields were used in this research.These data were observed and collected by well-trained agricultural technicians based on a standardized prescribed method [36].Not all the Chinese agrometeorological experimental stations (CAES) produced data suitable for the modeling simulations, thus, the CAES were chosen according to a strict set of procedures: (1) All the operations of the rice in the CAES were recorded from the years 1980-2010; (2) the selected CAES were good harvest sites and the typical yield exceeded 4000 kg/ha, in other words, the sites presented no obvious negative effects from pests and natural disasters; and (3) the same rice cultivars were planted for at least 3 years to ensure the stabilization of the model calibration and model validation.Generally, one year was chosen for model calibration and another two years were chosen for model validation [37].The selected CAES are shown in Figure 1 and detailed information on each site, with the corresponding cultivar, is given in Table 1.mature rice is mainly from July to October, after the harvest of early-mature rice.Favorable soil conditions, abundant precipitation, and a scientific cropping system all make Zhejiang a suitable natural location for rice cultivation.Zhejiang is one of the main rice production areas in the world, so the effects of climate change are more obvious here than elsewhere.The rice produced here is more vulnerable to these changes than in other areas, thus, the importance of assessing the impacts of climate change on rice yields in this region cannot be understated.

The CERES-Rice Model
The CERES-Rice model, which has been integrated into the DSSAT, has been in use for several decades.It has constantly been redesigned to facilitate new scientific advances and applications [38,39].The DSSAT design and approaches are used for various applications, and the model has been used in hundreds of published studies [17,40].The model is a successful crop model that has been widely used by researchers worldwide, and it is claimed to be one of the few simulation models that can simulate crop sequences, such as crop growth, development, and yield, in relation to varying levels of weather and genetics [41].The model is used for research across Asia, such as in the Philippines, Indonesia, and Thailand, as well as in areas of temperate climate, such as Japan, India, and Australia.Hence, the CERES-Rice model was expected to behave perfectly to fulfill the simulation and test the mitigation strategies [42].The CERES-Rice model performs professionally when applied to decisions on uncertain issues caused by variables in rice growth simulations.Thus, the latest CERES-Rice model, version 4.7, was adopted.
In the simulation model, crop growth could be divided into three parts according to the growing degree days (GDD): In the equation, T base , T high , and T opt represent the baseline, extremely high, and optimal temperatures during the rice growth period, respectively [43].When the temperature was between the basal temperature and the optimal temperature, the rice would grow very well and reach the highest yield.However, when the temperature was between an extreme high or low, the rice could hardly survive the environment and the result would be an obvious reduction in rice yield.

Dataset for the DSSAT
The DSSAT model requires standard input data, such as soil data, weather data, cultivar genetic coefficients, and agricultural experimental data.These are the basic input data required to guarantee rigorous calibration, validation, and evaluation of the model.Detailed information on each input data type is described in the following subsections.

Soil Data
The basic information on soil data included the soil color, average slope, potential runoff, and a fertility factor.In addition, the organic carbon, pH, cation exchange capacity, and total nitrogen of each layer were included.Detailed information on the soil was obtained from the China Soil Scientific Database (http://vdb3.soil.csdb.cn/).The organic matter in the soil and drainage conditions are crucial to rice growth.Detailed information for each site is provided in Appendix A.

Historical Weather Data
The historical weather dataset from experimental stations in this study was from the China Meteorological Data Sharing Service System (http://data.cma.cn/).The dataset includes highly resolved data, such as the daily solar radiation, daily maximum and minimum air temperatures, and daily precipitation from 8:00 am to 8:00 pm [44].The high quality of these historical weather data guaranteed reliable model calibration, validation, and evaluation.
The original weather data were converted to a standard daily solar radiation format based on the specific latitude and daily hours of sunlight at each site.Because solar radiation data were not available, the daily global solar radiation was calculated from the hours of sunlight by using the Ångström equation [45,46]: In the equation, R s represents the solar radiation (MJ/m 2 /d), n and N are the actual duration and the maximum duration of sunshine, respectively.R a is the extraterrestrial radiation (MJ/m 2 /d).Detailed information on the weather data is shown in Appendix B.

Cultivar Genetic Coefficients
Cultivar genetic coefficients are parameters that control the natural development of crops.For a given rice cultivar, the coefficients mainly control the growth stage of crops, such as the duration of flowering, the duration of maturation, and the yield [47]. Generally speaking, rice cultivar coefficients represents the nature of the genetic type, and each coefficient controls some aspects within the rice growing periods.Eight genetic coefficients for rice are defined in the CERES-Rice model, and the detailed information is given as Appendix C.

Climate Change Scenarios
Climate scenarios consisting of temperature, precipitation, and other climatic indicators are used to summarize the Representative Concentration Pathways (RCPs) for modeling in the next stages of climatic-based research [48,49].The climate scenarios used a fresh, new, coordinated parallel process, and the radiative forcing levels (W/m 2 ) by the end of the 21st century were selected to name the RCPs [50].Therefore, RCP2.6, RCP4.5, RCP6.0, and RCP8.5 were each named and used to assess potential climate change impacts.The global carbon dioxide (CO 2 ) concentration and emission for each RCP from the years 1980 to 2100, respectively, are shown in Figure 2. Conclusions from previous studies and the present research illustrated that changes in CO 2 , temperature, precipitation, and solar radiation in RCP2.6 and RCP6.0 might not be suitable for the simulations, whereas scenarios RCP4.5 and RCP8.5 could cover both medium and extreme scenarios.Therefore, the well-reasoned RCP4.5 and RCP8.5 were selected to conduct potential climate change impact simulations because they covered both medium and extreme scenarios.In this study, five GCMs named HadGEM2-ES, GFDL-ESM2M, IPSL-CM5A-LR, MIROC-ESM-CHEM, and NorESM1-M were used to generate future climate data (2011-2099) at the eight selected sites under scenarios RCP4.5 and RCP8.5.A detailed description of the selected GCMs is shown in Appendix D. In addition, the historical weather named as baseline weather data from the time scale (1980-2010) under RCP4.5 and RCP8.5 was extracted from the five GCMs to compare with weather from the future climate to prove the accuracy and stability of the model.

Calibration, Validation, and Evaluation of the CERES-Rice Model
Model calibration was done to verify that the genetic coefficients of the present rice cultivars were representative of normal levels.Model validation was conducted to test whether the confirmed genetic coefficients were suitable for evaluation.The duration of flowering, the duration of maturation, and the yield are the most commonly selected factors to confirm the genetic parameters [51].Each rice station needs to be calibrated and validated independently with a high-quality experimental record because the coefficients are cultivar specific and site specific.The generalized likelihood uncertainty estimation (GLUE) can calculate the genetic coefficients by using the observed data [52].The calibrated and validated specific rice genetic coefficients were then extracted from the CERES-Rice model to evaluate the potential rice yield under regional climate change scenarios (Figure 3).In this study, five GCMs named HadGEM2-ES, GFDL-ESM2M, IPSL-CM5A-LR, MIROC-ESM-CHEM, and NorESM1-M were used to generate future climate data (2011-2099) at the eight selected sites under scenarios RCP4.5 and RCP8.5.A detailed description of the selected GCMs is shown in Appendix D. In addition, the historical weather named as baseline weather data from the time scale (1980-2010) under RCP4.5 and RCP8.5 was extracted from the five GCMs to compare with weather from the future climate to prove the accuracy and stability of the model.

Calibration, Validation, and Evaluation of the CERES-Rice Model
Model calibration was done to verify that the genetic coefficients of the present rice cultivars were representative of normal levels.Model validation was conducted to test whether the confirmed genetic coefficients were suitable for evaluation.The duration of flowering, the duration of maturation, and the yield are the most commonly selected factors to confirm the genetic parameters [51].Each rice station needs to be calibrated and validated independently with a high-quality experimental record because the coefficients are cultivar specific and site specific.The generalized likelihood uncertainty estimation (GLUE) can calculate the genetic coefficients by using the observed data [52].The calibrated and validated specific rice genetic coefficients were then extracted from the CERES-Rice model to evaluate the potential rice yield under regional climate change scenarios (Figure 3).Indicators such as the normalized root-mean-square error (NRMSE) and the predicted deviation (PD) were selected to evaluate the precision of the simulations versus the observed data.The NRMSE and PD are each defined as follows [53]: In the equations,  represents the simulated results,  represents the observed results,  is the mean value of the observed data, n is the total number of comparisons, and  is the relative error.

Model Application
Rice grows and develops in a combination of complex, so non-countable variables should be considered for the input to perform well in the model simulation.The results were obtained with the constraint that future climate scenarios, CO2 concentrations, and adaptive measures were the only variables.The simulations were conducted to understand the impacts of climate change, the effects of CO2 as a fertilizer, and adaptive measures on rice yields in the Zhejiang province under generated weather conditions extracted from the five GCMs under RCP4.5 and RCP8.5, respectively.

Evaluation of the Impacts of Climate Change on Rice Yield
The impacts of climate change on rice phenology and yields were each simulated during the 2020s, 2050s, and 2080s.The simulation results were each compared with the average of the simulated results from the baseline period .The effects of elevated CO2 concentration on rice yields were extracted from the National Oceanic and Atmospheric Administration (Mauna Lua, Hawaii) CO2 database.
The mechanism equations for both early-mature and late-mature rice are given as follows: In the equations, , , and  represent the average temperature, precipitation, and solar radiation, respectively.In addition, P refers to the significance level and r is the correlation coefficient [54].Indicators such as the normalized root-mean-square error (NRMSE) and the predicted deviation (PD) were selected to evaluate the precision of the simulations versus the observed data.The NRMSE and PD are each defined as follows [53]: In the equations, S i represents the simulated results, O i represents the observed results, O i is the mean value of the observed data, n is the total number of comparisons, and PD i is the relative error.

Model Application
Rice grows and develops in a combination of complex, so non-countable variables should be considered for the input to perform well in the model simulation.The results were obtained with the constraint that future climate scenarios, CO 2 concentrations, and adaptive measures were the only variables.The simulations were conducted to understand the impacts of climate change, the effects of CO 2 as a fertilizer, and adaptive measures on rice yields in the Zhejiang province under generated weather conditions extracted from the five GCMs under RCP4.5 and RCP8.5, respectively.

Evaluation of the Impacts of Climate Change on Rice Yield
The impacts of climate change on rice phenology and yields were each simulated during the 2020s, 2050s, and 2080s.The simulation results were each compared with the average of the simulated results from the baseline period (1980-2010).The effects of elevated CO 2 concentration on rice yields were extracted from the National Oceanic and Atmospheric Administration (Mauna Lua, Hawaii) CO 2 database.

Adaptive Measures
The Fifth Assessment Report of the IPCC suggests several adaptation strategies to cope with potential climate change, such as altering planting dates, switching to high-temperature-tolerant cultivars, and breeding new cultivars [55,56].These measures are believed to be helpful, whereas the cropping systems and field management processes are held constant at the normal level in the future [57,58].Optimizing irrigation, controlling disease, and improving crop management practices were also suggested [59,60].
To optimize the planting dates, the simulations were conducted by advancing and delaying the planting date by 40 days at intervals of 5 days while other operations were held constant.The simulated rice yields were compared with the yields of rice when the planting dates were not changed.Switching to high-temperature-tolerant rice cultivars and breeding new rice cultivars were completed by changing the parameters in the model and then selecting the parameters that were able to achieve the highest yield, providing the most optimized rice cultivar.

Adaptive Measures
The Fifth Assessment Report of the IPCC suggests several adaptation strategies to cope with potential climate change, such as altering planting dates, switching to high-temperature-tolerant cultivars, and breeding new cultivars [55,56].These measures are believed to be helpful, whereas the cropping systems and field management processes are held constant at the normal level in the future [57,58].Optimizing irrigation, controlling disease, and improving crop management practices were also suggested [59,60].
To optimize the planting dates, the simulations were conducted by advancing and delaying the planting date by 40 days at intervals of 5 days while other operations were held constant.The simulated rice yields were compared with the yields of rice when the planting dates were not changed.Switching to high-temperature-tolerant rice cultivars and breeding new rice cultivars were completed by changing the parameters in the model and then selecting the parameters that were able to achieve the highest yield, providing the most optimized rice cultivar.The NRMSEs were all less than 15%, indicating good model performance and implying a reliable simulation [61].Details on the cultivar parameters after calibration and validation are summarized in Table 2. 1 The code is an abbreviation for the site name.The final L differentiates a late-maturation site from an early-maturation site.

Projected Climate Change
Comparisons of temperature, precipitation, solar radiation, and projected CO 2 concentration during future periods in the 2020s, 2050s, and 2080s, compared with the baseline period (1981-2010) under RCP4.5 and RCP8.5, are shown in Table 3.The mean annual average temperature relative to the baseline is expected to increase by 0.18 • C, 1.18 • C, and 1.28 • C, respectively, under RCP4.5, whereas it is projected to increase by 0.20 • C, 1.62 • C, and 1.87 • C, respectively, under RCP 8.5.The results show that the climate change variables are all in a continuously increasing trend under both RCP4.5 and RCP8.5.The results also showed that climatic variables always attained a higher level in RCP8.5 in almost all aspects.The same conclusion was reached with the projected CO 2 concentration.

Impacts of Climate Change on Rice Phenology
Figure 5 shows the average change in the flowering and maturation periods during future climate scenarios at all sites under RCP4.5 and RCP8.5, respectively.Almost all the simulations indicated that the climate change would inevitably advance rice phenology and shorten the duration of growth.More precisely, the rice flowering period from the year 2011 to the year 2099 would be shortened significantly by an average of 10 days and 12 days under RCP4.5 and RCP8.5, respectively.A maximum reduction of 20 days could possibly happen at site LYL under RCP4.5 and 18 days at site LSL under RCP8.5, respectively.
A maximum reduction of 20 days could possibly happen at site LYL under RCP4.5 and 18 days at site LSL under RCP8.5, respectively.The impacts on rice phenology attributable to climate change under RCP8.5 are more obvious than those under RCP4.5.The reason for this was that the duration of rice growth was shortened considerably by the increase in temperature.

Impacts of Climate Change on Rice Yields Without CO2 Fertilization Effects
The simulated impacts of climate change on rice yields, without accounting for the CO2 fertilization effects, under RCP4.5 and RCP8.5 during the 2020s, 2050s, and 2080s are shown in Figure 7.The rice yields appear to show a continually decreasing trend during all periods under both RCP4.5 and RCP8.5 at all sites.For all cultivars, the reductions in rice yields during the three periods were projected to be 8.79%, 14.02%, and 18.28% under RCP4.5, compared with yields from the baseline period.The decline in rice yields was projected to be 4.99%, 16.81%, and 32.58% under RCP8.5, compared with yields from the baseline period.The rice reduction at sites NH, PY, TZ, LSL, and TZL exceeded 30%, however, the most serious reduction, 48%, occurred at site PY under RCP8.5 in the 2080s.Additionally, the late-mature rice cultivars seemed to be more sensitive to future climate  A maximum reduction of 20 days could possibly happen at site LYL under RCP4.5 and 18 days at site LSL under RCP8.5, respectively.Figure 6 illustrates the maturation periods affected by projected climate change.The degree of decrease is nearly the same as that for the flowering periods.In all, the average duration of maturation would be shortened by 5.2 and 8.4 days under RCP4.5 and RCP8.5, respectively.A maximal reduction of 21 days could happen at site LSL under RCP4.5 and 24 days at site LSL under RCP8.5, indicating a decrease in the rice yield.The impacts on rice phenology attributable to climate change under RCP8.5 are more obvious than those under RCP4.5.The reason for this was that the duration of rice growth was shortened considerably by the increase in temperature.

Impacts of Climate Change on Rice Yields Without CO2 Fertilization Effects
The simulated impacts of climate change on rice yields, without accounting for the CO2 fertilization effects, under RCP4.5 and RCP8.5 during the 2020s, 2050s, and 2080s are shown in Figure 7.The rice yields appear to show a continually decreasing trend during all periods under both RCP4.5 and RCP8.5 at all sites.For all cultivars, the reductions in rice yields during the three periods were projected to be 8.79%, 14.02%, and 18.28% under RCP4.5, compared with yields from the baseline period.The decline in rice yields was projected to be 4.99%, 16.81%, and 32.58% under RCP8.5, compared with yields from the baseline period.The rice reduction at sites NH, PY, TZ, LSL, and TZL exceeded 30%, however, the most serious reduction, 48%, occurred at site PY under RCP8.5 in the 2080s.Additionally, the late-mature rice cultivars seemed to be more sensitive to future climate The impacts on rice phenology attributable to climate change under RCP8.5 are more obvious than those under RCP4.5.The reason for this was that the duration of rice growth was shortened considerably by the increase in temperature.

Impacts of Climate Change on Rice Yields without CO 2 Fertilization Effects
The simulated impacts of climate change on rice yields, without accounting for the CO 2 fertilization effects, under RCP4.5 and RCP8.5 during the 2020s, 2050s, and 2080s are shown in Figure 7.The rice yields appear to show a continually decreasing trend during all periods under both RCP4.5 and RCP8.5 at all sites.For all cultivars, the reductions in rice yields during the three periods were projected to be 8.79%, 14.02%, and 18.28% under RCP4.5, compared with yields from the baseline period.The decline in rice yields was projected to be 4.99%, 16.81%, and 32.58% under RCP8.5, compared with yields from the baseline period.The rice reduction at sites NH, PY, TZ, LSL, and TZL exceeded 30%, however, the most serious reduction, 48%, occurred at site PY under RCP8.5 in the 2080s.Additionally, the late-mature rice cultivars seemed to be more sensitive to future climate change than the early-mature ones were, presumably because late-mature rice will inevitably suffer more heat stress than early-mature rice.
change than the early-mature ones were, presumably because late-mature rice will inevitably suffer more heat stress than early-mature rice.To verify in detail, the differences in rice yield caused by meteorological data generated from different GCMs, the simulated rice yields using daily weather files were generated from five GCMs and are shown in Figure 7.The results indicate that the IPSL-CM5A-LR model always attains the highest rice yields, whereas the MIROC-ESM-CHEM model always results in the lowest yields.Despite the difference in rice yields caused by different GCMs, the difference was not entirely obvious.Thus, the results of all the models were combined together to obtain more detailed information and generate a box-picture (Figure 8).

Impacts of the CO2 Fertilization Effects on Rice Yields
Table 4 indicates that the increase in CO2 concentration is likely to mitigate the rice yields to some extent when CO2 fertilization is accounted for.The average rice yields at all experimental sites would change from −30.27% to −5.24% compared with the baseline yields.However, at most sites, the CO2 fertilization effects were unlikely to offset the total negative impacts from climate change.The reason for this was that CO2 fertilization could enhance photosynthesis by constraining photorespiration and could increase the utilization rate of water by reducing the stomata conduction and consequent transpiration.Furthermore, rice is a C3 plant and it responds to the effects of CO2 in photosynthetic carbon assimilation, thus, CO2 could contribute to rice yields [62].

Site
With CO2 Fertilization Effects (%) Without CO2 Fertilization Effects (%) To verify in detail, the differences in rice yield caused by meteorological data generated from different GCMs, the simulated rice yields using daily weather files were generated from five GCMs and are shown in Figure 7.The results indicate that the IPSL-CM5A-LR model always attains the highest rice yields, whereas the MIROC-ESM-CHEM model always results in the lowest yields.Despite the difference in rice yields caused by different GCMs, the difference was not entirely obvious.Thus, the results of all the models were combined together to obtain more detailed information and generate a box-picture (Figure 8).change than the early-mature ones were, presumably because late-mature rice will inevitably suffer more heat stress than early-mature rice.To verify in detail, the differences in rice yield caused by meteorological data generated from different GCMs, the simulated rice yields using daily weather files were generated from five GCMs and are shown in Figure 7.The results indicate that the IPSL-CM5A-LR model always attains the highest rice yields, whereas the MIROC-ESM-CHEM model always results in the lowest yields.Despite the difference in rice yields caused by different GCMs, the difference was not entirely obvious.Thus, the results of all the ned together to obtain more detailed information an

Impacts of the CO2 Fertilization Effects on Rice Yields
Table 4 indicates that the increase in CO2 concentration is likely to mitigate the rice yields to some extent when CO2 fertilization is accounted for.The average rice yields at all experimental sites would change from −30.27% to −5.24% compared with the baseline yields.However, at most sites, the CO2 fertilization effects were unlikely to offset the total negative impacts from climate change.The reason for this was that CO2 fertilization could enhance photosynthesis by constraining photorespiration and could increase the utilization rate of water by reducing the stomata conduction and consequent transpiration.Furthermore, rice is a C3 plant and it responds to the effects of CO2 in photosynthetic carbon assimilation, thus, CO2 could contribute to rice yields [62].

Impacts of the CO 2 Fertilization Effects on Rice Yields
Table 4 indicates that the increase in CO 2 concentration is likely to mitigate the rice yields to some extent when CO 2 fertilization is accounted for.The average rice yields at all experimental sites would change from −30.27% to −5.24% compared with the baseline yields.However, at most sites, the CO 2 fertilization effects were unlikely to offset the total negative impacts from climate change.The reason for this was that CO 2 fertilization could enhance photosynthesis by constraining photorespiration and could increase the utilization rate of water by reducing the stomata conduction and consequent transpiration.Furthermore, rice is a C 3 plant and it responds to the effects of CO 2 in photosynthetic carbon assimilation, thus, CO 2 could contribute to rice yields [62].
In the simulations, even when the CO 2 fertilization effects were considered, the negative impacts were still −20.643%and −24.1% at sites SX and PY, respectively, which both exceed 20%.Thus, some adaptive measures must be taken to lessen the adverse the potential impacts of climate change.

Analysis of the Impacts of Changing Climate Variables
Climate change will inevitably reduce rice yields at all sites, and the rice yields at sites PY, SX, and LYL in the 2080s, compared with the baseline period, were projected to decline by −24.1%, −20.64%, and −17.94%, respectively, under RCP8.5.From Equations ( 5) and ( 6), it could be observed that the mechanisms for early-mature rice and late-mature rice are totally different.The rice at sites PY and SX was early-mature rice, whereas that at site LYL was late-mature rice.The equations showed that the rice yields were closely connected with climatic variables during the growth of the rice.The average weather variables at sites PY, SX, and LYL are given in Table 5.The entire area is located in the Yangtze River Basin, which means that the river flows and rainfall is abundant at sites PY and SX.In addition, Table 5 clearly shows that merely the precipitation had changed from March to June during the 2020s, 2050s and 2080s, which indicates that the main influencing factors are the average solar radiation and average temperature.
For early-mature rice, the yield is positively correlated with the average solar radiation and negatively correlated with the average temperature because the yields are reduced.This result indicates that the negative impacts from the average temperature were stronger than the positive impacts of the average solar radiation.For late-mature rice, the main influence was also the negative impacts related to the average temperature.
The high temperature would certainly reduce the duration of rice growth by 11 to 15 days in the 2080s under RCP8.5.Therefore, the rice yields would be reduced by 20% to 40%.These simulation results are consistent with those by Zhang and Tao [63], who assessed climate change impacts with five models and found that the simulated growth period of rice would be shortened by 0.45 to 5.78 days.The increase in temperature during the rice growth period would reduce rice yields, and the decrease in solar radiation could reduce the crop photosynthesis rate and accumulation of biomass.The findings of this study are also similar to the results by Osborne et al. [64], who found that without the effects of CO 2 fertilization, the losses in total rice yields would be 9.7%, −1.5%, and −20.9% under scenario A2 in the 2020s, 2050s, and 2080s, respectively.The yields could reach a maximum at the optimal temperature of 25 • C and would be reduced by 10% for every 1 • C temperature increase above 25 • C until the temperature reached 35 to 36 • C, when no yields would be obtained [65].When the temperature rose above 25 • C, the rice yields would decline for the shorter grain-filling duration, thus, the maximum temperature during rice growth is a good indicator of heat stress [66].

Simulation of Adaptation Options
To reverse or lessen the potential reduction of rice yields, adaptive measures such as adjusting the planting dates and switching to high-temperature-tolerant cultivars were conducted at sites SX and PY, where the most adverse reductions in rice yields could possibly happen in the 2080s.

Adjusting the Rice Planting Dates
The rice flowering period is so sensitive to a high temperature that even an increase of 1 • C could cause a significant reduction in rice yields.The maximum temperature would increase gradually from 15 • to 40 • C from February to August at sites SX and PY.Because the temperature could be adjusted by switching planting dates, reasonably advancing the planting dates would probably help the late-mature rice survive the high temperature and extend the flowering period by several days, thereby stabilizing or increasing the rice yields.
The simulation results showed that the rice yields could be improved to a very impressive degree when the planting dates were advanced.The most suitable planting days for SX and PY are, respectively, 20 and 15 days earlier than the normal planting dates.Consequently, the rice yields could each apparently be improved by 9% and 15%, respectively (Figure 9). Figure 10 clearly shows the maximum temperature in the 2080s under the RCP8.5 scenario at the SX and PY sites.The maximum temperature increased significantly during this period at both sites.
impacts related to the average temperature.
The high temperature would certainly reduce the duration of rice growth by 11 to 15 days in the 2080s under RCP8.5.Therefore, the rice yields would be reduced by 20% to 40%.These simulation results are consistent with those by Zhang and Tao [63], who assessed climate change impacts with five models and found that the simulated growth period of rice would be shortened by 0.45 to 5.78 days.The increase in temperature during the rice growth period would reduce rice yields, and the decrease in solar radiation could reduce the crop photosynthesis rate and accumulation of biomass.The findings of this study are also similar to the results by Osborne et al. [64], who found that without the effects of CO2 fertilization, the losses in total rice yields would be 9.7%, -1.5%, and −20.9% under scenario A2 in the 2020s, 2050s, and 2080s, respectively.The yields could reach a maximum at the optimal temperature of 25 °C and would be reduced by 10% for every 1 °C temperature increase above 25 °C until the temperature reached 35 to 36 °C, when no yields would be obtained [65].When the temperature rose above 25 °C, the rice yields would decline for the shorter grain-filling duration, thus, the maximum temperature during rice growth is a good indicator of heat stress [66].

Simulation of Adaptation Options
To reverse or lessen the potential reduction of rice yields, adaptive measures such as adjusting the planting dates and switching to high-temperature-tolerant cultivars were conducted at sites SX and PY, where the most adverse reductions in rice yields could possibly happen in the 2080s.

Adjusting the Rice Planting Dates
The rice flowering period is so sensitive to a high temperature that even an increase of 1 °C could cause a significant reduction in rice yields.The maximum temperature would increase gradually from 15 °C to 40 °C from February to August at sites SX and PY.Because the temperature could be adjusted by switching planting dates, reasonably advancing the planting dates would probably help the late-mature rice survive the high temperature and extend the flowering period by several days, thereby stabilizing or increasing the rice yields.
The simulation results showed that the rice yields could be improved to a very impressive degree when the planting dates were advanced.The most suitable planting days for SX and PY are, respectively, 20 and 15 days earlier than the normal planting dates.Consequently, the rice yields could each apparently be improved by 9% and 15%, respectively (Figure 9). Figure 10 clearly shows the maximum temperature in the 2080s under the RCP8.5 scenario at the SX and PY sites.The maximum temperature increased significantly during this period at both sites.Because the high temperature would inevitably reduce the rice flowering period, a useful method would be to replace the present rice cultivar with a high-temperature-tolerant cultivar.After analyzing and considering the simulation results, the rice cultivar "Jiayu293" seemed to be more likely to survive the high temperature than the other cultivars during the flowering period.The rice cultivar "Jiayu296" was replaced with "Jiayu293" at the SX site, and the same method was performed at the PY site.
The simulation yields, when using the new cultivar at sites of SX and PY sites, were 4646 and 4914 kg/ha.When compared with the original yields of 4055 and 4126 kg/ha, they both had relative improvements of 14.57% and 19.1%.The simulation results for rice yields at the two sites where the new rice cultivars were adopted gained a significant enhancement, in contrast to the simulation results, where the rice cultivars were left unchanged.These results indicate that the flowering periods at both sites were prolonged and that rice cultivar "Jiayu296" was more tolerant to the high temperature than the original cultivar was.

Breeding New Rice Cultivars
Because genotype coefficients control the growth and development of rice cultivars, breeding new cultivars may be a very useful method of reducing the adverse impacts of climate change.There have been reports that new rice cultivars have proven very effective in increasing rice yields and coping with climate change [67].Rice coefficients are composed of several fundamental elements, such as P1, G1, and G2, and they are responsible for the GDD above 9 °C, the spikelet numbers per panicle, and the grain weight, respectively.
To determine the most suitable parameter that could raise the rice yield efficiently, the coefficients were switched within the settled range while the other conditions were held constant.Simulations of the yields were conducted with the altered genetic coefficients in the model at the SX and PY sites, where the lowest rice yield was projected to occur in the 2080s.The parameters P1, G1, and G2 were changed from the lowest to the highest within the settled range in the model: P1 was tested from 100 to 1000 at intervals of 100, G1 was tested from 20 to 290 at intervals of 30, and G2 was tested from 0.016 to 0.030 at intervals of 0.002.Each parameter was plotted with the corresponding rice yield at the SX and PY sites.Finally, the most suitable coefficients for each rice cultivar were conformed.Figure 11 clearly shows the relationships between each of the main rice coefficients and the rice yields.When parameter P1 approaches 900, the rice yield at site SX reaches the highest level,  Because the high temperature would inevitably reduce the rice flowering period, a useful method would be to replace the present rice cultivar with a high-temperature-tolerant cultivar.After analyzing and considering the simulation results, the rice cultivar "Jiayu293" seemed to be more likely to survive the high temperature than the other cultivars during the flowering period.The rice cultivar "Jiayu296" was replaced with "Jiayu293" at the SX site, and the same method was performed at the PY site.
The simulation yields, when using the new cultivar at sites of SX and PY sites, were 4646 and 4914 kg/ha.When compared with the original yields of 4055 and 4126 kg/ha, they both had relative improvements of 14.57% and 19.1%.The simulation results for rice yields at the two sites where the new rice cultivars were adopted gained a significant enhancement, in contrast to the simulation results, where the rice cultivars were left unchanged.These results indicate that the flowering periods at both sites were prolonged and that rice cultivar "Jiayu296" was more tolerant to the high temperature than the original cultivar was.

Breeding New Rice Cultivars
Because genotype coefficients control the growth and development of rice cultivars, breeding new cultivars may be a very useful method of reducing the adverse impacts of climate change.There have been reports that new rice cultivars have proven very effective in increasing rice yields and coping with climate change [67].Rice coefficients are composed of several fundamental elements, such as P1, G1, and G2, and they are responsible for the GDD above 9 • C, the spikelet numbers per panicle, and the grain weight, respectively.
To determine the most suitable parameter that could raise the rice yield efficiently, the coefficients were switched within the settled range while the other conditions were held constant.Simulations of the yields were conducted with the altered genetic coefficients in the model at the SX and PY sites, where the lowest rice yield was projected to occur in the 2080s.The parameters P1, G1, and G2 were changed from the lowest to the highest within the settled range in the model: P1 was tested from 100 to 1000 at intervals of 100, G1 was tested from 20 to 290 at intervals of 30, and G2 was tested from 0.016 to 0.030 at intervals of 0.002.Each parameter was plotted with the corresponding rice yield at the SX and PY sites.Finally, the most suitable coefficients for each rice cultivar were conformed.Figure 11 clearly shows the relationships between each of the main rice coefficients and the rice yields.When parameter P1 approaches 900, the rice yield at site SX reaches the highest level, whereas for PY, the most suitable P1 value was 546.Similar conclusions were obtained with parameters G1 and G2.
The best G1 parameters for both sites were 170 and 200, respectively, and those for G2 were 0.028 for both sites.
ainability 2019, 11, 2372 15 of 22 ameters G1 and G2.The best G1 parameters for both sites were 170 and 200, respectively, and se for G2 were 0.028 for both sites.Detailed information on the rice coefficients and their yield improvements for the selected sites hown in Table 6.The given genetic coefficients could be adopted to breed new rice cultivars.Detailed information on the rice coefficients and their yield improvements for the selected sites is shown in Table 6.The given genetic coefficients could be adopted to breed new rice cultivars.

Uncertainties
Like most simulations, this study also suffers from uncertainties and limitations, especially regarding the climate change, crop models, and agricultural techniques.Because knowing future agricultural techniques is beyond the capacity of this study, the uncertainties were mainly assessed for the climate change and crop models.
First, considering the uncertainties related to climate change, it is quite difficult to generate accurate daily weather from future climate scenarios and climate models.Comparability across climate impact studies could be compromised when little guidance is available for selecting future climate models [49].The impacts of climate change could vary significantly, so it is recommended that researchers utilize as many GCMs as possible.Although the future climate scenarios in this study were extracted from five GCMs under RCP4.5 and RCP8.5, uncertainties still remain, such as excess precipitation and heat.The uncertainties within the climate change scenarios of each model were assessed.Figure 8 shows that the maximum rice yield differences between GCMs was nearly 30% at sites SX and PY, and the minimum could be no more than 2% at site TZ.The results of this study are supported by previous studies.Li et al. [68] found that the predicted temperature contributed to approximately 59.7% of the uncertainty in rice yield projections.Krishnan et al. [69] predicted changes of −9.02%, −11.3%, and −21.35% in rice yields with future climate projections generated from the GFDL, GISS, and UKMO scenarios, respectively.Droogers et al. [70] assessed the climate change impacts on rice yields in the Volta Basin by using the SWAP and HadCM3 climate models, and the results showed that different scenarios had different impacts.
Second, the model mechanisms are critical to the accuracy of the simulation.Even though the CERES-Rice model is commonly believed to be the most reliable, the model cannot contain all the parameters or variables of the real environment.Aggarwal and Mall [71] evaluated the climate change impacts on rice yields in India by using the CERES-Rice model and the ORYZA2000 model.They found considerable differences in their assessments of the impacts of climate change on rice yields and therefore suggested using two crop models.Amiri et al. [17] studied the CERES-Rice, aqua crop, and ORYZA2000 models to examine the differences in performance when representing grain and rice yields.The CERES-Rice model represented the highest accuracy under most applications, the aqua crop model had the most precision in biological yield, and the ORYZA2000 model showed the lowest accuracy across the study.
Third, the adoption of various GCMs under RCPs and multiple simulation models could significantly reduce the uncertainties.Because more than one GCM was adopted under RCP4.5 and RCP8.5, the uncertainties from climate change were considered to have been reduced to the lowest.Even though CERES-Rice is considered the most reliable model, the uncertainty from the simulation model in this study is still unknown.The effects of extreme climatic events, such as strong winds, are not considered in the scenarios.An innovation that includes more than one simulating model should be applied in future work.

Conclusions
In this study, the potential climate change impacts and adaptation measures on rice yields were simulated by using the CERES-Rice model under the RCP4.5 and RCP8.5 scenarios, respectively.The results clearly showed a decreasing trend in the duration of flowering, the duration of maturation, and rice yields at all sites.However, in previous similar model-based studies, adaptive measures have seldom been simulated.In the present study, strategies such as advancing planting dates, exchanging rice cultivars that can tolerate high temperatures, and breeding fresh new rice varieties were evaluated at sites SX and PY, respectively.These results not only contribute to the evaluation of potential climate change impacts on rice yields, but they also provide adaptive strategies according to specific conditions.
New technologies used in agricultural and ecological applications are now being adopted more frequently than ever before.Therefore, unmanned aerial vehicle remote sensing, satellite-based remote sensing, and the modeling technique can be combined for a better prediction of agricultural yields.In addition, the model cannot handle the stresses of water, pests, and extremely high temperatures.To obtain the growth conditions of a crop quickly, sun-induced fluorescence could be adopted as an important indicator.In the near future, agricultural and ecological applications will boom.

Figure 1 .
Figure 1.Locations of selected Chinese agrometeorological experimental stations in the Zhejiang province.
forcing levels (W/m 2 ) by the end of the 21st century were selected to name the RCPs[50].Therefore, RCP2.6, RCP4.5, RCP6.0, and RCP8.5 were each named and used to assess potential climate change impacts.The global carbon dioxide (CO2) concentration and emission for each RCP from the years 1980 to 2100, respectively, are shown in Figure2.Conclusions from previous studies and the present research illustrated that changes in CO2, temperature, precipitation, and solar radiation in RCP2.6 and RCP6.0 might not be suitable for the simulations, whereas scenarios RCP4.5 and RCP8.5 could cover both medium and extreme scenarios.Therefore, the well-reasoned RCP4.5 and RCP8.5 were selected to conduct potential climate change impact simulations because they covered both medium and extreme scenarios.

Figure 3 .
Figure 3. General framework of the simulations used to determine the impacts of climate change.

Figure 3 .
Figure 3. General framework of the simulations used to determine the impacts of climate change.

Figure 4
Figure 4 clearly shows the results of the model calibration and validation for the observed and simulated duration of flowering, duration of maturation, and rice yields of the 11 sites.The NRMSEs for the duration of flowering, duration of maturation, and rice yields at all the sites were 11.12%, 10.42%, and 8.01%, respectively.

Figure 4
Figure 4 clearly shows the results of the model calibration and validation for the observed and simulated duration of flowering, duration of maturation, and rice yields of the 11 sites.The NRMSEs for the duration of flowering, duration of maturation, and rice yields at all the sites were 11.12%, 10.42%, and 8.01%, respectively.

Figure 4 .
Figure 4. Simulated and observed flowering period, maturation period, and yields at all sites.(a) Calibration results for the flowering period; (b) calibration results for the maturation period; and (c) calibration results for the rice yields.

Figure 5 .
Figure 5.Comparison of the change in flowering periods between (a) RCP4.5 and (b) RCP8.5.The numbers on the x-axis represent the rice cultivars, and those on the y-axis represent the change in flowering days.

Figure 6
Figure6illustrates the maturation periods affected by projected climate change.The degree of decrease is nearly the same as that for the flowering periods.In all, the average duration of maturation would be shortened by 5.2 and 8.4 days under RCP4.5 and RCP8.5, respectively.A maximal reduction of 21 days could happen at site LSL under RCP4.5 and 24 days at site LSL under RCP8.5, indicating a decrease in the rice yield.

Figure 6 .
Figure 6.Comparison of the change in maturation periods between (a) RCP4.5 and (b) RCP8.5.The numbers on the x-axis represent the rice cultivars, and those on the y-axis represent the change in days of maturation.

Figure 5 .
Figure 5.Comparison of the change in flowering periods between (a) RCP4.5 and (b) RCP8.5.The numbers on the x-axis represent the rice cultivars, and those on the y-axis represent the change in flowering days.

Figure 6
Figure6illustrates the maturation periods affected by projected climate change.The degree of decrease is nearly the same as that for the flowering periods.In all, the average duration of maturation would be shortened by 5.2 and 8.4 days under RCP4.5 and RCP8.5, respectively.A maximal reduction of 21 days could happen at site LSL under RCP4.5 and 24 days at site LSL under RCP8.5, indicating a decrease in the rice yield.

Figure 5 .
Figure 5.Comparison of the change in flowering periods between (a) RCP4.5 and (b) RCP8.5.The numbers on the x-axis represent the rice cultivars, and those on the y-axis represent the change in flowering days.

Figure 6 .
Figure 6.Comparison of the change in maturation periods between (a) RCP4.5 and (b) RCP8.5.The numbers on the x-axis represent the rice cultivars, and those on the y-axis represent the change in days of maturation.

Figure 6 .
Figure 6.Comparison of the change in maturation periods between (a) RCP4.5 and (b) RCP8.5.The numbers on the x-axis represent the rice cultivars, and those on the y-axis represent the change in days of maturation.

Figure 8 .
Figure 8. Future rice yields at each site, with future weather generated from a single general circulation model (GCM) for scenarios (a) RCP4.5 and (b) RCP8.5.

Figure 8 .
Figure 8. Future rice yields at each site, with future weather generated from a single general circulation model (GCM) for scenarios (a) RCP4.5 and (b) RCP8.5.

Figure 8 .
Figure 8. Future rice yields at each site, with future weather generated from a single general circulation model (GCM) for scenarios (a) RCP4.5 and (b) RCP8.5.

Figure 9 .
Figure 9. Change in rice yields in the 2080s under RCP8.5 at the (a) SX and (b) PY sites.The pluses and minuses on the horizontal ordinate represent the change in planting dates relative to the baseline date.

Figure 9 .
Figure 9. Change in rice yields in the 2080s under RCP8.5 at the (a) SX and (b) PY sites.The pluses and minuses on the horizontal ordinate represent the change in planting dates relative to the baseline date.

Figure 10 .
Figure 10.Maximum temperature in the 2080s under the RCP8.5 at sites (a) SX and (b) PY.The pluses and minuses on the horizontal ordinate represent the change in planting dates relative to the baseline date.

Figure 10 .
Figure 10.Maximum temperature in the 2080s under the RCP8.5 at sites (a) SX and (b) PY.The pluses and minuses on the horizontal ordinate represent the change in planting dates relative to the baseline date.

Table 1 .
Detailed descriptions of selected Chinese agrometeorological experimental stations (CAES) in the study area 1 .

Table 2 .
Estimated genetic coefficients for each rice cultivar 1 .

Table 3 .
Comparison of the change in climatic variables under RCP4.5 and 8.5 respectively 1 .

Table 4 .
Future (2080s) rice yields under scenario RCP8.5, with and without the CO2 fertilization effects, at all sites compared with the baseline period.

Table 4 .
Future (2080s) rice yields under scenario RCP8.5, with and without the CO2 fertilization effects, at all sites compared with the baseline period.

Table 4 .
Future (2080s) rice yields under scenario RCP8.5, with and without the CO 2 fertilization effects, at all sites compared with the baseline period.

Table 6 .
The most suitable genetic coefficients for rice cultivars at sites SX and PY.

Table 6 .
The most suitable genetic coefficients for rice cultivars at sites SX and PY.