Next Article in Journal
Consumption-Based Carbon Emissions of Tianjin Based on Multi-Scale Input–Output Analysis
Previous Article in Journal
Consumer Perception of Remanufactured Automotive Parts and Policy Implications for Transitioning to a Circular Economy in Sweden
Previous Article in Special Issue
Evaluating Social Media Response to Urban Flood Disaster: Case Study on an East Asian City (Wuhan, China)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Global Spatial Distributions of and Trends in Rice Exposure to High Temperature

1
Faculty of Geographical Science, Beijing Normal University, Beijing 100875, China
2
The Key Laboratory of Regional Geography, Beijing Normal University, Beijing 100875, China
3
College of Biologic and Geographic Sciences, Qinghai Normal University, Xining 810008, China
4
The State Key Laboratory of Earth Surface Processes and Resource Ecology, Beijing Normal University, Beijing 100875, China
*
Author to whom correspondence should be addressed.
Sustainability 2019, 11(22), 6271; https://0-doi-org.brum.beds.ac.uk/10.3390/su11226271
Submission received: 31 July 2019 / Revised: 31 October 2019 / Accepted: 4 November 2019 / Published: 8 November 2019
(This article belongs to the Special Issue Risk Assessment and Sustainable Development in Natural Hazards)

Abstract

:
Due to the effects of global warming, extreme temperature events are posing a great threat to crop yields, especially to temperature-sensitive crops such as rice. In the context of disaster risk theory, exposure is central to disaster prevention and reduction. Thus, a comprehensive analysis of crop exposure is essential to better reduce disaster effects. By combining the maximum entropy model (MaxEnt) and a multiple-criteria decision analysis (MCDA), this paper analyzed the global distribution and change in rice exposure to high temperature. The results showed the future states of rice after exposure to high temperatures. Our results are: (1) the areas of potential rice distribution zones decreased within the representative concentration pathway (RCP) scenarios RCP2.6 to RCP8.5 in MaxEnt, where the long-term (2061–2080) decreases are greater than those seen in the medium term (2041–2060). (2) In the future, the number of high temperature hazards in potential rice distribution areas increased. In the RCP8.5 scenario, the intensities of global high temperature hazards on rice were reduced because the total area of potential rice distribution zones decreased. (3) Through the view of barycenter shift, the barycenter of the global potential rice and high temperature hazard distributions showed a trend of backward motion, which meant the global rice exposure to high temperature was in a downward trend. With the background of global change, this paper has great significance for the mitigation of high temperature risk in rice and its effect on the potential security of future global rice production. Future research is warranted to concentrate on discussing more socioeconomic factors and increasing rice exposure change from the temporal vision.

1. Introduction

The latest Intergovernmental Panel on Climate Change (IPCC) special report estimates that if the current warming rate continues, the world will reach a human-induced global warming of 1.5 °C by approximately 2040 [1]. This increase in temperature will increase the frequency of extreme weather events, especially extreme high temperature events. Rice is among the world’s most important food crops, and it is a primary source of food for over 50% of the world’s population [2,3]. Frequent and extreme high temperature events will increase the risk of rice yield losses [4], which is a serious threat to the global food supply. Furthermore, according to the latest IPCC report, temperatures in the mid-latitudes (the most important crop producing areas) may rise by approximately 3 °C [5]. The sensitivity of the mid-latitudes means any temperature changes will significantly affect the crops growing therein. As a temperature-sensitive crop [6], rice will clearly be affected by high temperature disasters. Therefore, studies on the changes in rice distribution resulting from high temperature exposure, and changes under different representative concentration pathway (RCP) scenarios, can directly reflect the severity of high temperature hazards that rice may face in the future. In addition, the study of this topic can help prevent and reduce the risk of rice yield losses, which is of great significance to world food security [7].
Current research on disaster exposure are focused on quantitative assessments of exposure, including statistics related to the spatial distributions of hazard-bearing bodies [8,9], estimations of hazard-factor intensities [10,11,12,13,14], etc. Wu used the crop choice decision model to calculate the sown areas of major crops based on the economic values, and the changes in crop habitat factors caused by climate change were not taken into consideration [15]. Liu used the maximum entropy model (MaxEnt) to simulate rice suitability zones during the historical period (1980–2010), and achieved the classification of suitability grade according to the appearance frequency from the model results [16]. Some other scholars used a layer overlay analysis of geographic information system(GIS), and took the habitat index of crop distribution as layers to obtain crop suitability zones by multiple constraints and then they classified suitability degree though the accuracy of habitat index [17,18,19,20,21,22,23,24]. These studies provide a useful means for the calculation of rice suitability zones. On the other hand, Shi investigated the spatio-temporal variation of rice post-heading heat stress in South China during the period 1981–2010 [25]. Zhang conducted a comprehensive analysis on maize exposure, vulnerability, and adaptation to extreme temperature to understand the effects on global maize production, especially in major production countries [18]. Wang examined the changes in the area exposed to heat stress and historical movements of the geographical centroid of rice exposure to heat stress (EHS) during the period 1980–2015 across southern China [17]. These studies provide a scientific and efficiency measure of rice high temperature hazard and its variation. Current research usually focuses on a single factor, hazard-bearing body or hazard. Hence, analyses of the comprehensive characteristics of many factors related to exposure are few. Moreover, most study areas are located in China and southeast Asia, and most studies usually consider only a single temporal period [23,24,25], which makes it difficult to fully reflect the trends and regional differences in crop exposure. Therefore, the research on rice exposure to high temperature needs to start with two aspects of disaster-bearing bodies and hazard factors, by taking into account the impact of climate and socioeconomic factors on the rice planting range. According to the features of high temperature hazards in multiple scenarios and multiple periods, we can obtain the overall spatial patterns of rice exposure to high temperature, then reveal the hotspots and analyze how the exposure features help in finding trend rules.
Moreover, the Representative Concentration Pathway (RCP), which is commonly used in future crop yield simulation research, is a greenhouse gas concentration trajectory adopted by the IPCC for its fifth Assessment Report (AR5) in 2014. Dias identified the yield and growth changes in rice under the global climate change scenario RCP 8.5 [26]; Dar focused on simulating the effect of climate change on crop yield with the Decision Support System for Agrotechnology Transfer (DSSAT) v 4.6.1 under RCP4.5 scenario [27]; Chun applied a multi-scale crop modeling approach to assess the impacts of climate change on future rice yields in southeast Asia under the RCP 4.5 and RCP 8.5 scenarios [28]; Van Oort calculated the impacts of climate change on rice production in Africa under four kinds of RCP scenarios [29]; Kim predicted potential epidemics of rice leaf blast and sheath blight in South Korea under the RCP 4.5 and RCP 8.5 scenarios [30]; Zhang investigated the spatio-temporal change in extreme temperature stress across China under the RCP (2.6, 4.5, 6.0, and 8.5) scenarios [31]. Therefore, for future rice research, RCP provides a reliable and easy way for calculations, which is the reason why this paper chooses RCP as the future scenarios.
In order to study the changes in rice habitats in different future scenarios, with the background of climate change, in this study, our research was devoted to the assessment of global rice exposure to high temperature in medium-term (2041–2060) and long-term (2061–2080) periods. The results are based on RCP scenarios and share socioeconomic pathway (SSP) scenario data [32], and we used the MaxEnt model and multiple-criteria decision analysis (MCDA) methods to set up cost-effective and spatially targeted protection measures in order to safeguard world food supplies in the future. (1) In Section 3.1, we analyze the characteristics of area changes and barycenter shift of potential rice distribution under different scenarios. (2) In Section 3.2, we analyze the characteristics of intensity change and barycenter shifts of high temperature hazard factors under different scenarios. (3) In Section 3.3, we analyze rice exposure to high temperature by comparing the barycenter shifts of potential rice and high temperature hazard distributions. Throughout the study, we include references on future rice exposure for governments, so that they can rationally allocate resources and formulate reasonable farming policies to mitigate the high temperature risk for rice and its effect on the potential security of future global rice production.

2. Materials and Methods

2.1. Basic Idea and Research Framework

From the view of disaster risk theory [33], disaster risks can be divided into three categories: hazard (H), vulnerability (V) and exposure (E). Among these three types of disaster risk, exposure directly determines the effect of a disaster, which can, in some scenarios, reverse the results of the risk assessment. Generally, high degree exposure means a large disaster risk [34]. Therefore, we can estimate the potential risk for a specific disaster in a region by analyzing the extent of exposure therein.
In this paper, we presented and analyzed exposure maps of crops from two aspects: hazard-bearing bodies and hazard factors. For a hazard-bearing body, based on the MaxEnt model, we screened the rice biotope data to map potential rice distribution, and classified the suitability criteria for each considered region. The hazard factor analyses were based on four different scenario data (RCP2.6 and SSP1, RCP4.5 and SSP2, RCP6.0 and SSP4, and RCP8.5 and SSP5) during three periods (Baseline: 1971–2000, medium term: 2041–2060, and long term: 2061–2080). We calculated the intensities of the high temperature hazards in different periods and overlaid them with potential rice distribution to obtain the high temperature exposure experienced by rice under the different scenarios.
We also analyzed crop exposure from two aspects: the first is a change in potential rice distribution, which is described by the change rate of an area and the shift of its barycenter. The second considers changes in high temperature hazard intensity, which is denoted by the accumulative extreme high temperature in a year (SDD). Through shifts in potential rice distribution areas and high temperature hazard intensities, we obtained changing global rice exposure during different periods, thus describing the potential risk of high temperature on rice globally. The specific process is as shown in Figure 1.

2.2. Data Collection

The data were divided into three categories: disaster-formative environment data, hazard-bearing body data, and scenario data. The disaster-formative environment data include global digital elevation model (DEM), global-slope and global-soil parameter data [35]. The hazard-bearing body data included actual rice cultivated area data from 2000 and 2005, which were used to test the accuracy of potential rice distribution estimates. The scenario data included global meteorological data and global land use data.
In order to mitigate errors among the different climate model datasets, we used three sets of CMIP5 (Coupled Model Inter-comparison Project phase 5) meteorological data from the Inter-Sectoral Impact Model Intercomparison Project (ISI-MIP) [36,37,38]: IPSL-CM5A-LR (IPSL(Institute Pierre Simon Laplace) Earth System Model for the 5th IPCC report with low resolution) [39], MIROC-ESM-CHEM (Earth System Model for Interdisciplinary Research on Climate, atmospheric chemistry coupled version) [40] and NorESM1-M (Norwegian Earth System with intermediate resolution) [41]. The meteorological data contained the maximum temperature, minimum temperature, solar radiation intensity, wind speed, and relative humidity. The bioclimatic data contained the monthly averaged data of 19 parameters. In this paper, we considered the metrological data and bioclimatic data as rice habitat suitability factors to rate the rice habitat suitability. Moreover, we used the land use data to build socioeconomic scenarios, and we took the annual C3 crop distribution zones as te potential rice distribution zones.
We normalized the spatial resolution of all the data to 0.5° × 0.5°. Table 1 shows the contents of the data used in this paper.

2.3. Methods

2.3.1. Potential Rice Distribution Estimation

The MaxEnt model [42,43] and MCDA were used to evaluate future rice potential distributions [22,23,44]. First, this paper inputted the habitat parameters, which included the global DEM data, global slope data, global soil data (20 parameters) and bioclimatic data (19 parameters). Then, this paper screened the 13 parameters to consider only those with contribution rates that were greater than 95% (area under curve (AUC) values > 0.75) as the necessary habitat variables for potential rice distribution estimates.
Second, according to the quantiles of the habitat variables, the rice habitat suitability was ranked into five grades from I to V, where I was the lowest suitability, and V was the highest suitability. Table 2 shows the suitability grade rank definitions.
Third, this paper overlaid a habitat variables layer, and eliminated areas without C3 crop distributions to obtain potential rice distribution zones. Moreover, the regions of grade II to V were regarded as suitable zones for rice growth and the region of grade I was regarded as a non-suitable zone.

2.3.2. High Temperature Hazard Intensity Calculation

According to previous research results and generally accepted standards [45,46,47,48,49,50,51], this paper took the SDD as the high temperature hazard intensity index. SDD represents the annual accumulative extreme high temperature, where a high temperature event is defined as the daily average temperature Tas ≥ 30 °C. The calculation process can be divided into two steps as follows: first, this paper calculated the daily accumulative extreme high temperature according the high temperature standard (Tas ≥ 30 °C). Secondly, this paper summed the daily accumulative extreme high temperature to obtain the SDD, as shown in Equations (1) and (2):
D D i = { T i T b a s e f o r     T i T b a s e 0 f o r     T i T b a s e
S D D = i = 1 n D D i
where Ti is the average temperature on the ith day of the year, Tbase is the threshold temperature of the high temperature event (30 °C in this study) and n is the frequency of high temperature events during the year.

2.3.3. Barycenter Shift of Potential Rice and High Temperature Hazard Distributions

This paper used the barycenter shift [52,53] to measure and analyze the spatial distributions of the rice and high temperature hazards during the different periods. The barycenter of global rice habitat suitability was used to characterize the rice distribution, where the barycenter of an SDD was used to characterize the global high temperature hazard distribution, as respectively shown in Equations (3) and (4).
(1) Barycenter shift of rice habitat suitability calculation:
{ L O N s ,   t = i = 1 n ( S D i , t × l o n i , t ) i = 1 n S D i , t L A T s ,   t = i = 1 n ( S D i , t × l a t i , t ) i = 1 n S D i , t
where LONs and LATs are the longitude and latitude of the barycenter, respectively, SDi,t is the suitability of the ith grid with a suitability grade of t, loni,t and lati,t are the respective longitude and latitude of the ith grid with a suitability grade of t, and n is the number of the suitability grade t grids.
(2) Barycenter shift of the SDD calculation:
{ L O N h = i = 1 n ( S D D i × l o n i ) i = 1 n S D D i L A T h = i = 1 n ( S D D i × l a t i ) i = 1 n S D D i
where LONh and LATh are the longitude and latitude of the SDD, respectively, SDDi is the SDD value of the ith grid, loni and lati are the longitude and latitude of the ith grid, respectively, and n is the number of grids.
Moreover, we calculated the potential rice distribution and SDD change rates to quantify the changes in high temperature exposure experienced by rice, as shown in Equations (5) and (6), respectively:
D i f f a , t , i , j = A r e a j , t A r e a i , t A r e a i , t × 100 %
where Diffa,t,i,j is the change rate of the area of the t-grade zone from the ith period to the jth period, and Areai,t is the area of the t-grade zone in the ith period.
D i f f h , k , i , j = S D D j , k S D D i , k S D D i , k × 100 %
where Diffh,k,i,j is the SDD change rate of the kth grid from the ith period to the jth period, and SDDi,k is the SDD value of the kth grid.

2.3.4. Estimation of High Temperature Exposure

Using the previously calculated barycenter shift of a hazard-bearing body and the barycenter shift of a hazard, this paper expressed the characteristics of high temperature exposure and analyzed the potential risks to rice in terms of their direction and distance. This paper then compared the moving direction and distance of the barycenter in different periods, and qualitatively divided the high temperature exposure into five grades, as shown in Figure 2.

3. Results

We calculated three sets of climate data to obtain indicators maps under four scenarios (RCP2.6 and SSP1, RCP4.5 and SSP2, RCP6.0 and SSP4, and RCP8.5 and SSP4) in the medium-term and long-term future periods. Due to the similarities of our obtained results, only the results of the IPSL-CM5A-LR model are shown in the text. The results of the MIROC-ESM-CHEM and NorESM1-M models are detailed in the appendix.

3.1. Changes in the Area and Barycenter of Rice Suitability Zones

According to the estimation method provided in Section 2.3.1, we mapped the potential rice distribution, as shown in Figure A1, Figure A2 and Figure A3. On this basis, we used the continent as the basic unit to calculate the area of rice habitat suitability zones during different periods (Figure 3), and the barycenter shift from the baseline to the long-term period (Figure 4). The specific results are as follows:
  • Compared with the baseline period, the areas of the rice habitat suitability zones in Europe, Asia, Africa and South America decreased over all periods. Among them, the suitable areas of all grades in Europe showed a continuous downward trend, which decreased by 40% to 65%. The areas of the grade II and V zones in Asia decreased by over 45%, while the areas of the grade IV zones first decreased and then increased. The areas of the grade IV and V zones in Africa decreased by 30% to 40%, respectively, while the areas of the grade II zone increased significantly. The areas of the grade III to V zones in South America showed continuously decreasing trends, which decreased by 30% to 50%, and the decrease was faster from RCP2.6 to RCP8.5.
  • Compared with the baseline period, the areas of the rice habitat suitability zones in North America and Oceania increased. Between them, the areas of the grade IV and V zones in North America increased significantly, by 40% to 60%, while the areas of the grade II zones increased by over 50% in Oceania. However, under the RCP8.5 and SSP5 scenarios, the areas of the rice moderate suitability zones decreased by over 60%.
In terms of the direction of the barycenter shift, compared with the baseline period, we obtained the following results:
  • Under the different scenarios, the rice habitat suitability zones in Africa, Asia, North America, South America and Oceania all moved to higher latitudes. The barycenters of the higher-grade suitability zones (grade IV and V zones) and the lower-grade suitability zones (grade II and III zones) in Africa moved to the southeast. The barycenters of the higher-grade suitability zones in Asia moved to the northeast, and the barycenters of the lower-grade suitability zones moved to the northwest. The barycenters of the higher-grade and the lower-grade suitability zones in South America moved to the south. The barycenters of the higher-grade suitability zones in Oceania moved to the southeast, while those of the lower-grade suitability zones moved to the north.
  • Under different scenarios, the rice habitat suitability zones of all grades in Europe moved to lower latitudes, where the barycenters of both the higher-grade (IV and V) and lower-grade (II and III) suitable zones moved to the southwest.
In terms of the distance of the barycenter shift, compared with the baseline period, we obtained the following results:
  • For the higher-grade suitability zones, Europe had the largest range of movement, up to approximately 10°, while Oceania had the smallest range of movement of only approximately 1°. The ranges of movement for the other continents were approximately 4° to 6°.
  • For lower-grade suitability zones, South America had the largest movement range, up to approximately 7°, while the ranges of movement for the other continents were approximately 3° to 5°.
In this paper, we argue against the established view where the amplitude of a temperature increase is proportional to the barycenter northward-moving distance of a rice habitat suitability zone. Instead, we found that for all continents except Europe, the barycenter of rice suitability zones shifted the most under the RCP8.5 scenario. However, under the RCP6.0 scenario, the shifts in the barycenters of the rice habitat suitability zones were usually smaller than those under the RCP4.5 scenario, and even in Oceania, they were lower than those in the RCP2.6 scenario. Therefore, we found it reasonable to include a socioeconomic factor into the analysis.

3.2. Spatial Changes in the High temperature Hazard Distributions

(1) Spatial changes in the high temperature hazard distributions
Here, we used the SDD (Tas ≥ 30 °C) as a high temperature hazard intensity index and, through our calculations, we obtained the spatial distribution of SDD in the baseline period (Figure 5), and the medium-term and long-term periods (Figure 6) in potential rice distribution areas.
In terms of the high temperature hazard distributions in the baseline period, high-value regions were concentrated in the Indian subcontinent, Asia Minor, Sudanian Savanna, and the northern plains of the Mexican Plateau. The SDD values were generally above 100, and the SDDs in most parts of the Indian subcontinent were greater than 300. The SDDs in the southwestern parts of the Great Dividing Range, the western side of the La Plata plains, Indo-China, and the southern and eastern China regions were relatively high, generally above 10–50. The remaining regions such as Europe, north China and northeast China all contained low-value regions. Overall, Asia had the highest exposure, followed by Africa, Oceania, North America and South America, and Europe had the lowest amount of high temperature exposure.
Next, for the high temperature hazard distribution on the different continents during the medium-term and long-term periods we found: (1) the high-value regions in Africa were mainly concentrated in the Sudanian Savanna, but the high-value regions gradually decreased in the area from the RCP2.6 to RCP8.5 scenarios. In particular, under the RCP8.5 and SSP5 scenarios, the high-value regions moved to southern Africa. (2) The high-value regions in Asia were mainly concentrated in southern China, Asia Minor and the Indian subcontinent. Among them, the high-value regions of the Indian subcontinent shrunk and migrated to the north under the RCP8.5 scenario. (3) The high-value regions in North America were mainly located in the northern plains of the Mexican Plateau, and the values were high in the south and low in the north. (4) The high-value regions in South America were located in the southern parts of the Brazilian Plateau and the La Plata plains, where the region gradually expanded as the RCP2.6 to RCP8.5 scenarios were considered. (5) The high-value regions in Oceania ere located on the southwestern side of the Great Dividing Range. (6) The SDD values were low in Europe, but under the RCP8.5 and SSP5 scenarios, the Mediterranean coastal transitioned into a high-value region where the SDDs were generally above 50 and even above 100 in the north. Overall, the medium-term and long-term high temperature hazard distributions are generally consistent with those of the baseline period. In North America, South America and Europe, the intensities of the high temperature hazards increased. In Asia and Africa, the intensities of the high temperature hazards decreased.
(2) Changes in the high temperature hazards
By comparing the SDD values of the baseline period with those of the medium-term and long-term periods, and those of the medium-term period with those of the long-term period, the change rates of the SDD values under different scenarios (Figure 7) were calculated to characterize the change in the intensities of the high temperature hazards.
Except for Asia Minor, the SDD values of most regions showed a continuously increasing trend from the baseline period to the long-term period. Under the different scenarios, the high-value regions were located in Chinese monsoon regions, southeastern United States, south-central South America, southeastern Australia and the Iberian Peninsula. However, the SDD values on the Indian subcontinent and the Sudanian Savanna changed only rarely.
The regional high temperature exposure also changed under the different scenarios. In central North America, southern Africa, northeast and north China, Mediterranean coastal regions, and the east coast of South America all changed from non-exposed regions to exposed regions, while the Indian subcontinent, Indo-China, the Sudanian Savanna, the Caspian Depression, the Black Sea lowlands, the Danube River Basin, and the West Indies all changed from exposed regions to non-exposed regions.
(3) Barycenter shift of high temperature hazards to rice
Using the equations provided in Section 2.3.3, we calculated the positions of each continent’s SDD barycenter in the different periods and scenarios to characterize the barycenter shift of the SDDs from the baseline period to the long-term period, as shown in Figure 8.
In terms of the barycenter shift directions, compared with the baseline period, the SDD barycenters in Africa, North America and Oceania moved to the southeast, the SDD barycenter of Asia moved to the east, the SDD barycenter of South America moved to the south, and the SDD barycenter of Europe moved to the northeast. In terms of the barycenter shift ranges, compared with the baseline period, Europe had the largest range of SDD barycenter movement, with an average distance of 4° to 6°. Oceania had the smallest range of SDD barycenter movement, with an average distance lower than 1°. The rest of the continents had average distances of 2° to 4°. Except for Europe, the range of the barycenter shift of high temperature hazards on other continents enlarged from RCP2.6 to RCP8.5.

3.3. Comprehensive Analysis of Rice Exposure to High Temperature

According to the method presented in Section 2.3.4, the severity level of future rice exposure to high temperature was obtained (Table 3). The specific calculation process is shown in Appendix A.7. In all scenarios, the high temperature exposure of the higher-grade rice habitat suitability zones in Asia were aggravated the most, which were followed by the zones in Africa. The high temperature exposure of all grades of rice habitat suitability zones in South America and Europe were alleviated, as well as the higher-grade rice habitat suitability zones in North America. Overall, the potential risks of high temperature hazards on rice seriously increased in Asia, especially on the Indian subcontinent, Indo-China and the Chinese monsoon regions. The potential risks of high temperature hazards on rice significantly reduced in southern Europe, south-central South America, and the southeastern United States.

4. Discussion

4.1. Accuracy Tests of Potential Rice Distribution Estimates

In order to test the accuracy of the estimated results of potential rice distribution areas, we compared the estimated results with existing sets of rice harvesting and seeding data (EARTHSTAT 2000 harvest range data, SPAM 2000 and 2005 harvest range and sowing range data, and MIRCA 2000 harvest range data), as shown in Table 4.
The estimated potential rice distribution zones were generally larger than the actual rice distributions (Figure 9). The reason is that although North America and other places were potential areas of rice growth, due to the local cultivation histories and eating habits, rice cultivation was lower than the predictions, which reduced the accuracy of the results given in this paper. The accuracy of the producers of each set of data was greater than 80%, indicating that the rice potential distribution zones calculated in this paper contained over 80% of the actual rice distribution areas. So, the above-mentioned potential rice distribution zone estimation method had good accuracy and can reflect the actual rice distribution cultivation in the future.

4.2. Temporal Rice Exposure to High Temperature

From the view of disaster risk theory, exposure can be defined by two aspects: spatial and temporal. The spatial aspect was analyzed in the previous sections. The temporal aspect emphasizes the occurrence of hazards during the full growth period of crops or the key stages of crop growth. In terms of the effects of high temperature hazards on rice, the high temperature hazard intensity during a whole growth period is the leading factor that affects rice loss. With global climate change, the temporal distribution of extreme high temperature events during the rice growth period will also change.
Based on the rice growth season data from SAGE (University of Wisconsin-Madison Sustainability and the Global Environment) [54], we calculated the distribution of high temperature exposure and SSD changes during the growth period of rice, and the barycenter shifts of the SDDs during the whole growth period, and then analyzed the temporal characteristics of rice exposure to high temperature. The SDD barycenter calculation results were classified into three stages of rice growth periods: vegetative phase (VE), reproductive phase (RE) and ripening phase (RI). SDD barycenter distribution during the whole rice growth period can be calculated using Equation (7):
D = i = 1 n S D D i × d i i = 1 n S D D i
where SDDi is the high temperature hazard index of day d, n is the number of days with SDD values > 0 during the whole growth period, and D is the SDD barycenter distribution during the whole growth period.
(1)
Distribution of SDDs during the rice growth period
(i)
The spatial distributions of the SDDs during the rice growth period were roughly the same as those of the SDD values determined without considering temporal exposure. The high-value regions during the growth period were in the Indian subcontinent, Mesopotamia, southern areas of the Sahara Desert, and the northern plateau of Mexico, although their scopes shrank. In particular, southern China shrunk the most because the original data used in this region only consisted of local, early-season rice growth period data.
(ii)
Changes in the SDD values without considering temporal exposure were consistent with the changes in SDDs during the rice growth period, but the change rates of the SDDs during the growth period were larger, especially on the Indian subcontinent and southern China.
(2)
Barycenter shift of the SDDs during the growth period
By comparing the SDD barycenter shifts during the rice growth periods from the baseline to long-term periods, we found that there were no significant changes in most parts of the world: the SDDs in Europe, Africa, North America and Oceania were concentrated in the RE period, and in the VE period for the Mexican plateau. The SDDs in South America were concentrated in the VE period. In Asia, there was a regional difference, where those in southern China were concentrated in the RI period, and those in northern China were concentrated in the RE period; the rest were concentrated in the VE period. However, barycenter fluctuations occurred in some areas: the SDD barycenter in the northeastern United States showed a floating change of VE to RE to VE, while the SDD barycenter in northeast China transitioned from the RE period to the RI period. The southern area of the Brazilian plateau presented RE to VE to RE floating changes. In summary, major rice-growing areas in Asia need to focus on high temperature risk prevention in the future, and according to the periods of high temperature hazards, targeted high temperature protection mitigation measures should be carried out at key growth stages [55,56].

4.3. Comparative Analysis with Related Research

According to the results of our study, with the background of the CMIP5 climate models and socioeconomic scenarios, the global rice exposure presented the differences in spatial patterns and trends in various periods. In order to ensure the rationality of the results obtained in this paper, we conducted a comparative analysis with some related research. (1) In terms of the hazard-bearing body: Wu [15] simulated the rice-growing areas for the years 2000 and 2035 using the crop choice decision model; the primary rice production regions of the results, such as the Indian subcontinent, southeast Asia, southern China, southern United States, southeastern Brazil and the Sudanian Savanna are basically consistent with the higher-grade rice suitability zones (V and IV regions) of the baseline and medium-term periods in our study. In the studies on local regions, the results showed that the area and grade of rice suitability zones in southern China were decreasing, but they were increasing in northern China. Approximately 18.2% of the new rice-growing areas are in the newly suitable areas of the northern single-cropping rice system, and the northern limit of Chinese rice production has moved northward from 47° N to 52° N. The results of our study also presented these trends above from baseline to the long-term period. (2) In terms of the hazard factor: Gourdji [11] analyzed the regional features of global crops facing extreme high temperature, and found that during the period from the 2000s to 2050s, days with over critical temperatures in the Indian subcontinent and southern China increased by over 3–5 times, while these regions of our study showed a 4–5 times increase trend. In the studies on local regions, the results showed that the high-value regions of heat stress and its variation in southern China during the period 1980–2010 were the provinces of Jiangxi, Hunan and Guangdong, while the geographical barycenters of heat stress showed an eastward or northeastward shift. The trends in our study are also consistent with theirs. The comparative analysis above can indicate to a certain extent that the spatial patterns of and trends in rice exposure to high temperature obtained in our study are rational.

4.4. Limitations of This Study

4.4.1. Regional Indicator Selection and Suitability Level Division

Habitat factors have different effects on rice production in different regions. Hence, the principle of adapting to local conditions should be adopted in the screening process. In particular, the setting of the quantile, which is used in suitability grade ranking, should also reflect any regional characteristics. For example, in the northern part of the Indian subcontinent, the annual precipitation is higher than in other major rice-growing areas, which makes the regional suitability lower than the actual situation. Therefore, in future research, crop potential distribution estimation needs to be optimized from two levels: (1) mining regional characteristic habitat variables, and (2) introducing weighted elements (such as crop yield) to set the grade quantiles.

4.4.2. Consideration of Socioeconomic Factors

In this paper, we only considered land use as socioeconomic scenario data. However, changes in population, markets and eating habits, which limit the area and distribution of crop cultivation and hence affect crop prices, demand and planting willingness [57,58,59], were only marginally considered. Therefore, in order to make the estimated results more consistent with reality, a comprehensive multi-factor exposure estimation method is imperative. In addition, the level of fortification in the considered socioeconomic factors is often neglected, and any increase in the level of fortification can greatly reduce exposure [60], which should be considered in future research.

5. Conclusions

By combining the MaxEnt model and MCDA, in this paper, we mapped the potential rice distribution under different scenarios during different periods. Based on the above-mentioned results, we calculated the characteristic variations in rice exposure to high temperature from the perspective of disaster risk theory. The specific conclusions are as follows:
(1)
The areas of potential rice distribution zones changed:
(i)
From the perspective of area changes, under different scenarios, the rice potential distribution zones continuously declined. Europe had the largest decline (approximately 50%), followed by Asia, Africa and South America (with averages of 30% to 40%). Among them, the higher-grade suitability areas in Asia and Africa dropped significantly, with a decrease of 40%–50%. There were significant declines in the areas of all suitability grades in South America, while the areas of potential distribution zones in North America and Oceania increased significantly (by approximately 30%).
(ii)
From the perspective of barycenter shifts, the global potential rice distribution zones generally moved to higher latitudes, except in Europe, which moved to lower latitudes. The movement ranges in terms of latitude were the largest in Europe, 6–10°, while Oceania had the smallest range of movement, 1–3°.
(2)
The effects of rice exposure to high temperature, and the intensities of the high temperature hazards, had the following characteristics:
(i)
High temperature exposure was concentrated on the Indian subcontinent, Asia Minor, Sudanian Savanna, the north plains of the Mexican Plateau, La Plata plains, the southwestern parts of the Great Dividing Range in Australia, and southern China, where the SDDs were high, with averages of 100–500.
(ii)
The SDD values in potential rice distribution zones showed continuous and significant increases. The SDDs increased from RCP2.6 to RCP8.5, especially on the Indian subcontinent, the monsoon region of China, southeastern United States, south-central South America, the southeastern parts of Australia and the Iberian Peninsula, where the SDD increased by over 3 to 5 times.
(iii)
High temperature exposure on the Indian subcontinent, Indo-China and the Sudanese prairie region all decreased significantly, while the Mediterranean Sea, central and southern Africa, the east coast of South America, northeast China, and central North America all increased.
(3)
The characteristics of the comprehensive changes in high temperature exposure and their effect on rice are summarized as follows. In each scenario, the degree of exposure change on each continent showed a floating variation. From the overall characteristics of high temperature exposure, changes in Asia were significantly aggravated, and South America showed significant reductions. In Africa, the performance was generally aggravated, while the rest of the continents generally had fewer or no significant changes.
In summary, according to the abovementioned changes, the areas of rice suitability zones in the main producing areas of Asia, especially the Indian subcontinent, Indo-China and the eastern monsoon region of China, were significantly reduced, and the intensities of the high temperature hazards increased significantly. The potential high temperature risks in these regions are expected to be severe in the future. To combat the expected risks, key monitoring and disaster prevention should be carried out. The suitable areas for rice crops in North America will expand, while the degree of exposure will reduce. With this prediction, it may be prudent to try to expand the existing range of rice cultivation. In view of the obvious reduction in exposure in South America, cultivation can be appropriately modified. Finally, with the background of global change, this paper has great significance for the mitigation of high temperature risk in rice and its effect on the potential security of future global rice production. For example, this paper provides a spatial distribution of global rice exposure to high temperature under RCP and SSP scenarios, which is conducive to the rational allocation of funds and resources in the disaster management of rice cultivation. What is more, the results of this study can also provide a reference for the future regional transfer of rice cultivation at the high temperature hazard level.
Further study can take the following two aspects into consideration. (1) To improve the simulation accuracy of rice suitability zones, more socioeconomic factors should be added into the constraints in the terms of reasonable and quantitative calculation (such as market supply and demand, farmers’ planting willingness and residents’ eating habits). In addition, the occurrence frequency and intensity of natural disasters (such as extreme temperatures, droughts and floods) also affects the actual crop-cultivated distribution. Therefore, the level of crop field fortification will be among the important factors in the estimation and evaluation of rice suitability zones. On the other hand, by considering the regional differences in geographical environment and rice cultivars, the sectional method of selecting habitat variables in different regions can be adopted to improve the simulation accuracy. (2) To expand the connotation of rice exposure, this study only took rice cultivation zones as the exposure subject. In the future research, we should consider the climatic and socioeconomic changes comprehensively, and use scientific methods to estimate the rice-growing areas, yields and economic values, so as to better enrich the essential content of rice exposure and use these data in the risk assessment of high temperature hazard to rice.

Author Contributions

Funding acquisition, J.W.; Conceptualization, R.W. and J.W.; Methodology, R.W. and Y.J.; Supervision, R.W. and J.W.; Validation, R.W., and Y.J. and P.S.; Formal analysis, P.S. and R.W.; Writing—original draft preparation, R.W., P.S. and Y.J.; Writing—review and editing, R.W., P.S. and J.W.

Funding

This research was funded by the National Key Research and Development Program of China (Grant No. 2016YFA0602402) and the National Natural Science Foundation of China (Grant No. 41671501).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Appendix A.1. Pattern of Rice Suitability Zones

Figure A1. Distribution of rice suitability zones with the IPSL-CM5A-LR model.
Figure A1. Distribution of rice suitability zones with the IPSL-CM5A-LR model.
Sustainability 11 06271 g0a1aSustainability 11 06271 g0a1b
Figure A2. Distribution of rice suitability zones with the MIROC-ESM-CHEM model.
Figure A2. Distribution of rice suitability zones with the MIROC-ESM-CHEM model.
Sustainability 11 06271 g0a2aSustainability 11 06271 g0a2b
Figure A3. Distribution of rice suitability zones with the NorESM1-M model.
Figure A3. Distribution of rice suitability zones with the NorESM1-M model.
Sustainability 11 06271 g0a3aSustainability 11 06271 g0a3b

Appendix A.2. Changes in the Area of the Rice Suitability Zone

Figure A4. Changes in the area of the rice suitability zone with the MIROC-ESM-CHEM model.
Figure A4. Changes in the area of the rice suitability zone with the MIROC-ESM-CHEM model.
Sustainability 11 06271 g0a4aSustainability 11 06271 g0a4b
Figure A5. Changes in the area of the rice suitability zone with the NorESM1-M model.
Figure A5. Changes in the area of the rice suitability zone with the NorESM1-M model.
Sustainability 11 06271 g0a5aSustainability 11 06271 g0a5b

Appendix A.3. Changes in the Barycenter of the Rice Suitability Zone

Figure A6. Changes in the barycenter of the rice suitability zone with the MIROC-ESM-CHEM model.
Figure A6. Changes in the barycenter of the rice suitability zone with the MIROC-ESM-CHEM model.
Sustainability 11 06271 g0a6aSustainability 11 06271 g0a6b
Figure A7. Changes in the barycenter of the rice suitability zone with the NorESM1-M model.
Figure A7. Changes in the barycenter of the rice suitability zone with the NorESM1-M model.
Sustainability 11 06271 g0a7aSustainability 11 06271 g0a7b

Appendix A.4. Pattern of Rice Exposure to High Temperature

Figure A8. Pattern of rice exposure to high temperature with the MIROC-ESM-CHEM model.
Figure A8. Pattern of rice exposure to high temperature with the MIROC-ESM-CHEM model.
Sustainability 11 06271 g0a8aSustainability 11 06271 g0a8b
Figure A9. Pattern of rice exposure to high temperature with the NorESM1-M model.
Figure A9. Pattern of rice exposure to high temperature with the NorESM1-M model.
Sustainability 11 06271 g0a9aSustainability 11 06271 g0a9b

Appendix A.5. Changes in the High Temperature Hazards (SDD)

Figure A10. Change rate (%) of SDD with the MIROC-ESM-CHEM model.
Figure A10. Change rate (%) of SDD with the MIROC-ESM-CHEM model.
Sustainability 11 06271 g0a10
Figure A11. Change rate (%) of SDD with the NorESM1-M model.
Figure A11. Change rate (%) of SDD with the NorESM1-M model.
Sustainability 11 06271 g0a11

Appendix A.6. Changes in the Barycenter of the High Temperature Hazards (SDD)

Figure A12. Changes in the barycenter of SDD with the MIROC-ESM-CHEM model.
Figure A12. Changes in the barycenter of SDD with the MIROC-ESM-CHEM model.
Sustainability 11 06271 g0a12aSustainability 11 06271 g0a12b
Figure A13. Changes in the barycenter of SDD with the NorESM1-M model.
Figure A13. Changes in the barycenter of SDD with the NorESM1-M model.
Sustainability 11 06271 g0a13aSustainability 11 06271 g0a13b

Appendix A.7. Comprehensive Characteristics of the Changes in High Temperature Exposure

Table A1. Changes in the direction of the barycenter between rice suitability zone and high temperature hazard (SDD) with the IPSL-CM5A-LR model.
Table A1. Changes in the direction of the barycenter between rice suitability zone and high temperature hazard (SDD) with the IPSL-CM5A-LR model.
ContinentsSYDRCP2.6 and SSP1RCP4.5 and SSP2RCP6.0 and SSP4RCP8.5 and SSP5Avg
BM-MLBM-MLBM-MLBM-MLBM-ML
AFHigh+++0+0+0+0
Low0+00000+0+
ASHigh0+0+-0000+
Low0+0-0+--00
EUHigh-0-0-0-+-0
Low-+-0-0-0-0
NAHigh-0-0-00+-+
Low-+-+-0-+-+
OCHigh0+0+0+0+00
Low000000000+
SAHigh-0-0-+---0
Low-----0----
Note: “+” indicates the same direction (S); “-“ indicates the opposite direction (O); “0“ indicates medium direction (between S and O).
Table A2. Changes in the distance of the barycenter between rice suitability zone and high temperature hazard (SDD) with the IPSL-CM5A-LR model.
Table A2. Changes in the distance of the barycenter between rice suitability zone and high temperature hazard (SDD) with the IPSL-CM5A-LR model.
ContinentsSYDRCP2.6& and SSP1RCP4.5& and SSP2RCP6.0& and SSP4RCP8.5& and SSP5Avg
BM-MLBM-MLBM-MLBM-MLBM-ML
AFHigh-0-0-0-0-0
Low+0+0+0+-+-
ASHigh-++-+-0++0
Low+++++0+++-
EUHigh-+---+-+-+
Low-+---+-+-0
NAHigh-0-0-0-0-0
Low0000+00000
OCHigh-0-0-0-0-0
Low+0+0+0--0-
SAHigh-0---0----
Low+0-------0
Note: “+” indicates the same direction (S); “-“ indicates the opposite direction (O); “0“ indicates medium direction (between S and O).
Table A3. Changes in the high temperature exposure with the IPSL-CM5A-LR model.
Table A3. Changes in the high temperature exposure with the IPSL-CM5A-LR model.
ContinentsSYDRCP2.6& and SSP1RCP4.5& and SSP2RCP6.0& and SSP4RCP8.5& and SSP5Avg
B-M-LB-M-LB-M-LB-M-LB-M-L
AFHighHLHLHLHLHL
LowHHNHNHHNHN
ASHighHNHNLNNHHH
LowHHLHHHLHNN
EUHighLNLLLNNNLN
LowNNLLLNLNLL
NAHighLLLLLLHLNL
LowLNNNLHNNNN
OCHighNHNHNHNHNL
LowNHNHNHNLHL
SAHighLLLLNLLLLL
LowLHLLLLLLLL
Table A4. Changes in the direction of the barycenter between rice suitability zone and high temperature hazard (SDD) with the MIROC-ESM-CHEM model.
Table A4. Changes in the direction of the barycenter between rice suitability zone and high temperature hazard (SDD) with the MIROC-ESM-CHEM model.
ContinentsSYDRCP2.6& and SSP1RCP4.5& and SSP2RCP6.0& and SSP4RCP8.5& and SSP5Avg
BM-MLBM-MLBM-MLBM-MLBM-ML
AFHigh0++-+-+++-
Low0++-+-+++-
ASHigh-0-+-0-0-0
Low-0-++0-000
EUHigh+0-0+00++0
Low-0000+00-0
NAHigh-00+-0-0-0
Low0-0+0+0000
OCHigh-+-000-000
Low-+00-0+000
SAHigh-0-0-0-0-0
Low-0-0-0-0-0
Note: “+” indicates the same direction (S); “-“ indicates the opposite direction (O); “0“ indicates medium direction (between S and O).
Table A5. Changes in the distance of the barycenter between rice suitability zone and high temperature hazard (SDD) with the MIROC-ESM-CHEM model.
Table A5. Changes in the distance of the barycenter between rice suitability zone and high temperature hazard (SDD) with the MIROC-ESM-CHEM model.
ContinentsSYDRCP2.6& and SSP1RCP4.5& and SSP2RCP6.0& and SSP4RCP8.5& and SSP5Avg
BM-MLBM-MLBM-MLBM-MLBM-ML
AFHigh-0---0---0
Low+0+0+0--0-
ASHigh++-++-+++0
Low+0+0+0++++
EUHigh+--++-0-0-
Low-+0--+-00-
NAHigh-+0+000+00
Low0++00++++0
OCHigh-0-0-0++-0
Low+0+-+0+-+0
SAHigh-+--------
Low-+--------
Note: “+” indicates the same direction (S); “-“ indicates the opposite direction (O); “0“ indicates medium direction (between S and O).
Table A6. Changes in the high temperature exposure with the MIROC-ESM-CHEM model.
Table A6. Changes in the high temperature exposure with the MIROC-ESM-CHEM model.
ContinentsSYDRCP2.6& and SSP1RCP4.5& and SSP2RCP6.0& and SSP4RCP8.5& and SSP5Avg
BM-MLBM-MLBM-MLBM-MLBM-ML
AFHighHLNLNLHLNL
LowHHNHNHHLNL
ASHighLHNNLNLHLH
LowLHNHHHLHNH
EUHighHNLNHNHLHL
LowLNNLHNNLLL
NAHighLNHHLNLHLN
LowLHHHHHNHNH
OCHighNLLLNLLHNL
LowNHNNLHHNNH
SAHighLNLLLLLLLL
LowLNLLLLLLLL
Table A7. Severity grade of the changes in the high temperature exposure with the MIROC-ESM-CHEM model.
Table A7. Severity grade of the changes in the high temperature exposure with the MIROC-ESM-CHEM model.
ContinentsSYDRCP2.6& and SSP1RCP4.5& and SSP2RCP6.0& and SSP4RCP8.5& and SSP5Avg
BM-MLBM-MLBM-MLBM-MLBM-ML
AFHighIIIIVIVIIIIV
LowIIIIIIIIIV
ASHighIIIIIIIVIIIIII
LowIIIIIIIIIII
EUHighIIIVIIIIIIII
LowIVIVIIIVV
NAHighIVIIVIIIIV
LowIIIIIIIII
OCHighIVVIVIIIIV
LowIIIIIIIIIIII
SAHighIVVVVV
LowIVVVVV
Table A8. Changes in the direction of the barycenter between rice suitability zone and high temperature hazard (SDD) with the NorESM1-M model.
Table A8. Changes in the direction of the barycenter between rice suitability zone and high temperature hazard (SDD) with the NorESM1-M model.
ContinentsSYDRCP2.6& and SSP1RCP4.5& and SSP2RCP6.0& and SSP4RCP8.5& and SSP5Avg
BM-MLBM-MLBM-MLBM-MLBM-ML
AFHigh0000000+00
Low-000-00+0-
ASHigh00-0---+00
Low00-+-0-+-0
EUHigh00-0000+-0
Low00-0000+00
NAHigh0+00-+000+
Low-0-0-0-0-+
OCHigh0000+00000
Low000+0+000+
SAHigh-0-+-----0
Low-0-+-----+
Note: “+” indicates the same direction (S); “-“ indicates the opposite direction (O); “0“ indicates medium direction (between S and O).
Table A9. Changes in the distance of the barycenter between rice suitability zone and high temperature hazard (SDD) with the NorESM1-M model.
Table A9. Changes in the distance of the barycenter between rice suitability zone and high temperature hazard (SDD) with the NorESM1-M model.
ContinentsSYDRCP2.6& and SSP1RCP4.5& and SSP2RCP6.0& and SSP4RCP8.5& and SSP5Avg
BM-MLBM-MLBM-MLBM-MLBM-ML
AFHigh-0-0------
Low----------
ASHigh+-+-+-+0+-
Low+0+0+-+++0
EUHigh---0-+-+-+
Low---0-0-0--
NAHigh-0-+-0-0-0
Low0-00000+-0
OCHigh0-000--000
Low+-+0+-+++0
SAHigh-+-0------
Low---------0
Note: “+” indicates the same direction (S); “-“ indicates the opposite direction (O); “0“ indicates medium direction (between S and O).
Table A10. Changes in the high temperature exposure with the NorESM1-M model.
Table A10. Changes in the high temperature exposure with the NorESM1-M model.
ContinentsSYDRCP2.6& and SSP1RCP4.5& and SSP2RCP6.0& and SSP4RCP8.5& and SSP5Avg
BM-MLBM-MLBM-MLBM-MLBM-ML
AFHighNLNLNLHLNL
LowLLNLLLHLLL
ASHighNNLNLNNHNN
LowNHNHLNNHLH
EUHighNLLLNNHNLN
LowNLLLNLHLNL
NAHighHLNNNLNLHL
LowLLLNLNLHNL
OCHighNLNNHLNLNN
LowNNHHHNNHHH
SAHighLNNLLLLLLL
LowLLNLLLLLNL
Table A11. Severity grade of the changes in the high temperature exposure with the NorESM1-M model.
Table A11. Severity grade of the changes in the high temperature exposure with the NorESM1-M model.
ContinentsSYDRCP2.6& and SSP1RCP4.5& and SSP2RCP6.0& and SSP4RCP8.5& and SSP5Avg
BM-MLBM-MLBM-MLBM-MLBM-ML
AFHighIVIVIVIIIIV
LowVIVVIIIV
ASHighIIIIVIVIIIII
LowIIIIIVIIIII
EUHighIVVIIIIIIV
LowIVVIVIIIIV
NAHighIIIIIIIVIVIII
LowVIVIVIIIIV
OCHighIVIIIIIIIVIII
LowIIIIIIIII
SAHighIVIVVVV
LowVIVVVIV

Appendix A.8. Pattern of High Temperature Hazard (SDD) in the Growth Period

Figure A14. Distribution of SDDs in growth period with the IPSL-CM5A-LR model.
Figure A14. Distribution of SDDs in growth period with the IPSL-CM5A-LR model.
Sustainability 11 06271 g0a14aSustainability 11 06271 g0a14b
Figure A15. Distribution of SDDs in growth period with the MIROC-ESM-CHEM model.
Figure A15. Distribution of SDDs in growth period with the MIROC-ESM-CHEM model.
Sustainability 11 06271 g0a15aSustainability 11 06271 g0a15b
Figure A16. Distribution of SDDs in growth period with the NorESM1-M model.
Figure A16. Distribution of SDDs in growth period with the NorESM1-M model.
Sustainability 11 06271 g0a16aSustainability 11 06271 g0a16b

Appendix A.9. Changes in High Temperature Hazard (SDD) in the Growth Period

Figure A17. Changes in SDDs in the growth period with the IPSL-CM5A-LR model.
Figure A17. Changes in SDDs in the growth period with the IPSL-CM5A-LR model.
Sustainability 11 06271 g0a17
Figure A18. Changes in SDDs in the growth period with the MIROC-ESM-CHEM model.
Figure A18. Changes in SDDs in the growth period with the MIROC-ESM-CHEM model.
Sustainability 11 06271 g0a18
Figure A19. Changes in SDDs in the growth period with the NorESM1-M model.
Figure A19. Changes in SDDs in the growth period with the NorESM1-M model.
Sustainability 11 06271 g0a19

Appendix A.10. Pattern of the Barycenter of High Temperature Hazard (SDD) in the Growth Period

Figure A20. Distribution of the barycenter of SDDs in the growth period with the IPSL-CM5A-LR model. Note: NS indicates unsuitable zone; NE indicates non-exposed zone.
Figure A20. Distribution of the barycenter of SDDs in the growth period with the IPSL-CM5A-LR model. Note: NS indicates unsuitable zone; NE indicates non-exposed zone.
Sustainability 11 06271 g0a20
Figure A21. Distribution of the barycenter of SDDs in the growth period with the MIROC-ESM-CHEM model.
Figure A21. Distribution of the barycenter of SDDs in the growth period with the MIROC-ESM-CHEM model.
Sustainability 11 06271 g0a21
Figure A22. Distribution of the barycenter of SDDs in the growth period with the NorESM1-M model.
Figure A22. Distribution of the barycenter of SDDs in the growth period with the NorESM1-M model.
Sustainability 11 06271 g0a22
Figure A23. Distribution of the barycenter of SDDs in the growth period (baseline).
Figure A23. Distribution of the barycenter of SDDs in the growth period (baseline).
Sustainability 11 06271 g0a23

Appendix A.11. Changes in the Barycenter of High Temperature hazard (SDD) in the Growth Period

Figure A24. Changes in the barycenter of SDDs in the growth period with the IPSL-CM5A-LR model. Note: NC indicates unchanged zone; N2E indicates non-exposed zone to exposed zone; E2N indicates exposed zone to non-exposed zone.
Figure A24. Changes in the barycenter of SDDs in the growth period with the IPSL-CM5A-LR model. Note: NC indicates unchanged zone; N2E indicates non-exposed zone to exposed zone; E2N indicates exposed zone to non-exposed zone.
Sustainability 11 06271 g0a24
Figure A25. Changes in the barycenter of SDDs in the growth period with the MIROC-ESM-CHEM model.
Figure A25. Changes in the barycenter of SDDs in the growth period with the MIROC-ESM-CHEM model.
Sustainability 11 06271 g0a25
Figure A26. Changes in the barycenter of SDDs in the growth period with the NorESM1-M model.
Figure A26. Changes in the barycenter of SDDs in the growth period with the NorESM1-M model.
Sustainability 11 06271 g0a26

References

  1. IPCC. Global Warming of 1.5 °C, Summary for Policy Makers; IPCC: Geneva, Switzerland, 2018. [Google Scholar]
  2. Khush, G.S. Origin, dispersal, cultivation and variation of rice. Plant Mol. Biol. 1997, 35, 25–34. [Google Scholar] [CrossRef]
  3. Akinbile, C.; Akinlade, G.; Abolude, A. Trend analysis in climatic variables and impacts on rice yield in Nigeria. J. Water Clim. Change 2015, 6, 534–543. [Google Scholar] [CrossRef] [Green Version]
  4. Peng, S.; Huang, J.; Sheehy, J.E.; Laza, R.C.; Visperas, R.M.; Zhong, X.; Centeno, G.S.; Khush, G.S.; Cassman, K.G. Rice yields decline with higher night temperature from global warming. Proc. Natl. Acad. Sci. USA 2004, 101, 9971–9975. [Google Scholar] [CrossRef] [Green Version]
  5. Hoegh-Guldberg, O.; Jacob, D.; Taylor, M.; Bindi, M.; Brown, S.; Camilloni, I.; Diedhiou, A.; Djalante, R.; Ebi, K.; Engelbrecht, F. Impacts of 1.5 °C global warming on natural and human systems. In Global Warming of 1.5 °C, Summary for Policy Makers; IPCC: Geneva, Switzerland, 2018. [Google Scholar]
  6. Fahad, S.; Adnan, M.; Hassan, S.; Saud, S.; Hussain, S.; Wu, C.; Wang, D.; Hakeem, K.R.; Alharby, H.F.; Turan, V. Rice responses and tolerance to high temperature. In Advances in Rice Research for Abiotic Stress Tolerance; Elsevier: Amsterdam, The Netherlands, 2019; pp. 201–224. [Google Scholar]
  7. Fraser, S.A.; Wood, N.J.; Johnston, D.; Leonard, G.S.; Greening, P.D.; Rossetto, T. Variable population exposure and distributed travel speeds in least-cost Tsunami evacuation modelling. Nat. Hazards Earth Syst. Sci. 2014, 14, 2975–2991. [Google Scholar] [CrossRef]
  8. Antwi-Agyei, P.; Fraser, E.D.; Dougill, A.J.; Stringer, L.C.; Simelton, E. Mapping the vulnerability of crop production to drought in Ghana using rainfall, yield and socioeconomic data. Appl. Geogr. 2012, 32, 324–334. [Google Scholar] [CrossRef]
  9. Gunasekera, R.; Ishizawa, O.; Aubrecht, C.; Blankespoor, B.; Murray, S.; Pomonis, A.; Daniell, J. Developing an adaptive global exposure model to support the generation of country disaster risk profiles. Earth Sci. Rev. 2015, 150, 594–608. [Google Scholar] [CrossRef] [Green Version]
  10. Gbetibouo, G.A.; Ringler, C.; Hassan, R. Vulnerability of the south African farming sector to climate change and variability: An indicator approach. Nat. Resour. Forum 2010, 34, 175–187. [Google Scholar] [CrossRef]
  11. Gourdji, S.M.; Sibley, A.M.; Lobell, D.B. Global crop exposure to critical high temperatures in the reproductive period: Historical trends and future projections. Environ. Res. Lett. 2013, 8, 024041. [Google Scholar] [CrossRef]
  12. Sun, Y.; Zhang, X.; Zwiers, F.W.; Song, L.; Wan, H.; Hu, T.; Yin, H.; Ren, G. Rapid increase in the risk of extreme summer heat in eastern China. Nat. Clim. Change 2014, 4, 1082. [Google Scholar] [CrossRef]
  13. Murthy, C.; Laxman, B.; Sai, M.S. Geospatial analysis of agricultural drought vulnerability using a composite index based on exposure, sensitivity and adaptive capacity. Int. J. Disaster Risk Reduct. 2015, 12, 163–171. [Google Scholar] [CrossRef]
  14. Anandhi, A.; Steiner, J.L.; Bailey, N. A system’s approach to assess the exposure of agricultural production to climate change and variability. Clim. Change 2016, 136, 647–659. [Google Scholar] [CrossRef]
  15. Wu, W.; Shibasaki, R.; Yang, P.; Tan, G.; Matsumura, K.-I.; Sugimoto, K. Global-scale modelling of future changes in Sown areas of major crops. Ecol. Model. 2007, 208, 378–390. [Google Scholar] [CrossRef]
  16. Liu, Z.; Yang, P.; Tang, H.; Wu, W.; Zhang, L.; Yu, Q.; Li, Z. Shifts in the extent and location of rice cropping areas match the climate change pattern in China during 1980–2010. Reg. Environ. Change 2015, 15, 919–929. [Google Scholar] [CrossRef]
  17. Wang, P.; Hu, T.; Kong, F.; Zhang, D. Changes in the spatial pattern of rice exposure to heat stress in China over recent decades. Clim. Change 2019, 154, 1–12. [Google Scholar] [CrossRef]
  18. Zhang, L.; Zhang, Z.; Chen, Y.; Wei, X.; Song, X. Exposure, vulnerability, and adaptation of major maize-growing areas to extreme temperature. Nat. Hazards 2018, 91, 1257–1272. [Google Scholar] [CrossRef]
  19. Arun-Chinnappa, K.S.; Ranawake, L.; Seneweera, S. Impacts and management of temperature and water stress in crop plants. In Abiotic Stress Management for Resilient Agriculture; Springer: Berlin, Germany, 2017; pp. 221–233. [Google Scholar]
  20. Huang, J.; Zhang, F.; Xue, Y.; Lin, J. Recent changes of rice heat stress in Jiangxi province, Southeast China. Int. J. Biometeorol. 2017, 61, 623–633. [Google Scholar] [CrossRef]
  21. Fahad, S.; Hussain, S.; Saud, S.; Hassan, S.; Ihsan, Z.; Shah, A.N.; Wu, C.; Yousaf, M.; Nasim, W.; Alharby, H. Exogenously applied plant growth regulators enhance the morpho-physiological growth and yield of rice under high temperature. Front. Plant Sci. 2016, 7, 1250. [Google Scholar] [CrossRef]
  22. Mendas, A.; Delali, A. Integration of multicriteria decision analysis in gis to develop land suitability for agriculture: Application to durum wheat cultivation in the region of Mleta in Algeria. Comput. Electron. Agric. 2012, 83, 117–126. [Google Scholar] [CrossRef]
  23. Walke, N.; Reddy, G.O.; Maji, A.; Thayalan, S. Gis-based multicriteria overlay analysis in soil-suitability evaluation for cotton (Gossypium spp.): A case study in the black soil region of Central India. Comput. Geosci. 2012, 41, 108–118. [Google Scholar] [CrossRef]
  24. You, L.; Wood, S.; Wood-Sichra, U. Generating plausible crop distribution maps for Sub-saharan Africa using a spatially disaggregated data fusion and optimization approach. Agric. Syst. 2009, 99, 126–140. [Google Scholar] [CrossRef]
  25. Shi, P.; Tang, L.; Wang, L.; Sun, T.; Liu, L.; Cao, W.; Zhu, Y. Post-heading heat stress in rice of South China during 1981–2010. PLoS ONE 2015, 10, e0130642. [Google Scholar] [CrossRef]
  26. Dias, M.; Navaratne, C.; Weerasinghe, K.; Hettiarachchi, R. Application of DSSAT crop simulation model to identify the changes of rice growth and yield in Nilwala river basin for mid-centuries under changing climatic conditions. Proc. Food Sci. 2016, 6, 159–163. [Google Scholar] [CrossRef]
  27. Dar, M.U.D.; Aggarwal, R.; Kaur, S. Effect of climate change scenarios on yield and water balance components in rice-wheat cropping system in central Punjab, India. J. Agrometeorol. 2017, 19, 226–229. [Google Scholar]
  28. Chun, J.A.; Li, S.; Wang, Q.; Lee, W.-S.; Lee, E.-J.; Horstmann, N.; Park, H.; Veasna, T.; Vanndy, L.; Pros, K. Assessing rice productivity and adaptation strategies for Southeast Asia under climate change through multi-scale crop modeling. Agric. Syst. 2016, 143, 14–21. [Google Scholar] [CrossRef]
  29. Van Oort, P.A.; Zwart, S.J. Impacts of climate change on rice production in Africa and causes of simulated yield changes. Glob. Change Biol. 2018, 24, 1029–1045. [Google Scholar] [CrossRef]
  30. Kim, K.-H.; Cho, J.; Lee, Y.H.; Lee, W.-S. Predicting potential epidemics of rice leaf blast and sheath blight in South Korea under the RCP 4.5 and RCP 8.5 climate change scenarios using a rice disease epidemiology model, EPIRICE. Agric. Forest Meteorol. 2015, 203, 191–207. [Google Scholar] [CrossRef]
  31. Zhang, Z.; Chen, Y.; Wang, C.; Wang, P.; Tao, F. Future extreme temperature and its impact on rice yield in China. Int. J. Climatol. 2017, 37, 4814–4827. [Google Scholar] [CrossRef]
  32. IPCC. Emissions Scenarios; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  33. Shi, P. Disaster Risk Science; Springer: Berlin, Germany, 2019. [Google Scholar]
  34. He, C.; Huang, Q.; Dou, Y.; Tu, W.; Liu, J. The population in China’s earthquake-prone areas has increased by over 32 million along with rapid urbanization. Environ. Res. Lett. 2016, 11, 074028. [Google Scholar] [CrossRef]
  35. Batjes, N.H. ISRIC-WISE Derived Soil Properties on a 5 by 5 Arc-Minutes Global Grid (ver. 1.2); ISRIC-World Soil Information: Wageningen, The Netherlands, 2012. [Google Scholar]
  36. Warszawski, L.; Frieler, K.; Huber, V.; Piontek, F.; Serdeczny, O.; Schewe, J. The inter-sectoral impact model intercomparison project (ISI–MIP): Project framework. Proc. Natl. Acad. Sci. USA 2014, 111, 3228–3232. [Google Scholar] [CrossRef]
  37. Wiebe, K.; Lotze-Campen, H.; Sands, R.; Tabeau, A.; van der Mensbrugghe, D.; Biewald, A.; Bodirsky, B.; Islam, S.; Kavallari, A.; Mason-D’Croz, D. Climate change impacts on agriculture in 2050 under a range of plausible socioeconomic and emissions scenarios. Environ. Res. Lett. 2015, 10, 085010. [Google Scholar] [CrossRef]
  38. Zhao, C.; Piao, S.; Wang, X.; Huang, Y.; Ciais, P.; Elliott, J.; Huang, M.; Janssens, I.A.; Li, T.; Lian, X. Plausible rice yield losses under future climate warming. Nat. Plants 2017, 3, 16202. [Google Scholar] [CrossRef]
  39. Dufresne, J.-L.; Foujols, M.-A.; Denvil, S.; Caubel, A.; Marti, O.; Aumont, O.; Balkanski, Y.; Bekki, S.; Bellenger, H.; Benshila, R. Climate change projections using the IPSL-CM5 earth system model: From CMIP3 to CMIP5. Climate Dyn. 2013, 40, 2123–2165. [Google Scholar] [CrossRef]
  40. Watanabe, S.; Hajima, T.; Sudo, K.; Nagashima, T.; Takemura, T.; Okajima, H.; Nozawa, T.; Kawase, H.; Abe, M.; Yokohata, T. MIROC-ESM 2010: Model description and basic results of CMIP5–20c3m experiments. Geosci. Model Dev. 2011, 4, 845. [Google Scholar] [CrossRef]
  41. Bentsen, M.; Bethke, I.; Debernard, J.; Iversen, T.; Kirkevåg, A.; Seland, Ø.; Drange, H.; Roelandt, C.; Seierstad, I.; Hoose, C. The Norwegian Earth system model, norESM1-M-part 1: Description and basic evaluation. Geosci. Model Dev. Discus. 2012, 5, 2843–2931. [Google Scholar] [CrossRef]
  42. Elith, J.; Graham, C.; Anderson, R.; Dudík, M.; Ferrier, S.; Guisan, A.; Hijmans, R.; Huettmann, F.; Leathwick, J.; Lehmann, A. Novel methods improve prediction of species’ distributions from occurrence data. Ecography 2006, 29, 129–151. [Google Scholar] [CrossRef]
  43. Phillips, S.J.; Anderson, R.P.; Schapire, R.E. Maximum entropy modeling of species geographic distributions. Ecol. Model. 2006, 190, 231–259. [Google Scholar] [CrossRef] [Green Version]
  44. Lane, A.; Jarvis, A. Changes in climate will modify the geography of crop suitability: Agricultural biodiversity can help with adaptation. ISRICAT 2007, 4, 1–12. [Google Scholar]
  45. Gao, S.; Wang, P. Impacts of High Temperature Stress on Rice Production in the Mid-Lower Yangtze River Valley of China; Meteorological Press: Beijing, China, 2009. (In Chinese) [Google Scholar]
  46. Lobell, D.B.; Burke, M.B. On the use of statistical models to predict crop yield responses to climate change. Agric. Forest Meteorol. 2010, 150, 1443–1452. [Google Scholar] [CrossRef]
  47. Sun, W.; Huang, Y. Global warming over the period 1961–2008 did not increase high-temperature stress but did reduce low-temperature stress in irrigated rice across China. Agric. Forest Meteorol. 2011, 151, 1193–1201. [Google Scholar] [CrossRef]
  48. Butler, E.E.; Huybers, P. Adaptation of us maize to temperature variations. Nat. Clim. Change 2013, 3, 68. [Google Scholar] [CrossRef]
  49. Zhang, Z.; Wang, P.; Chen, Y.; Song, X.; Wei, X.; Shi, P. Global warming over 1960–2009 did increase heat stress and reduce cold stress in the major rice-planting areas across China. Eur. J. Agron. 2014, 59, 49–56. [Google Scholar] [CrossRef]
  50. Xiong, W.; Feng, L.; Ju, H.; Yang, D. Possible impacts of high temperatures on China’s rice yield under climate change. Adv. Earth Sci. 2016, 31, 515–528. [Google Scholar]
  51. ISO. Animal and Vegetable Fats and Oils—Determination of Acid Value and Acidity; China Standard Press: Beijing, China, 2005. [Google Scholar]
  52. Rafiee, R.; Mahiny, A.S.; Khorasani, N.; Darvishsefat, A.A.; Danekar, A. Simulating urban growth in Mashad city, Iran through the SLEUTH model (UGM). Cities 2009, 26, 19–26. [Google Scholar] [CrossRef]
  53. Shimba, T. An estimation of center of gravity from force platform data. J. Biomech. 1984, 17, 53–60. [Google Scholar] [CrossRef]
  54. Sacks, W.J.; Deryng, D.; Foley, J.A.; Ramankutty, N. Crop planting dates: An analysis of global patterns. Glob. Ecol. Biogeogr. 2010, 19, 607–620. [Google Scholar] [CrossRef]
  55. Rezaei, E.E.; Webber, H.; Gaiser, T.; Naab, J.; Ewert, F. Heat stress in cereals: Mechanisms and modelling. Eur. J. Agron. 2015, 64, 98–113. [Google Scholar] [CrossRef]
  56. Sánchez, B.; Rasmussen, A.; Porter, J.R. Temperatures and the growth and development of maize and rice: A review. Glob. Change Biol. 2014, 20, 408–417. [Google Scholar] [CrossRef]
  57. Fischer, G.; Shah, M.; Tubiello, F.; Van Velhuizen, H. Socio-economic and climate change impacts on agriculture: An integrated assessment, 1990–2080. Philos. Trans. R. Soc. B Biol. Sci. 2005, 360, 2067–2083. [Google Scholar] [CrossRef]
  58. Yan, Y.; Zhao, J.; Deng, H.; Luo, Q. Predicting China’s cultivated land resources and supporting capacity in the twenty-first century. Int. J. Sustain. Dev. World Ecol. 2006, 13, 229–241. [Google Scholar] [CrossRef]
  59. Wei, X.; Declan, C.; Erda, L.; Yinlong, X.; Hui, J.; Jinhe, J.; Ian, H.; Yan, L. Future cereal production in China: The interaction of climate change, water availability and socio-economic scenarios. Glob. Environ. Change 2009, 19, 34–44. [Google Scholar] [CrossRef]
  60. Howden, S.M.; Soussana, J.-F.; Tubiello, F.N.; Chhetri, N.; Dunlop, M.; Meinke, H. Adapting agriculture to climate change. Proc. Natl. Acad. Sci. USA 2007, 104, 19691–19696. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Outline of the study structure.
Figure 1. Outline of the study structure.
Sustainability 11 06271 g001
Figure 2. High temperature exposure change assessment structure. Note: we defined the barycenter moving direction in the same quadrant as the direction movement, and the barycenter moving direction in a different quadrant as the different direction movement.
Figure 2. High temperature exposure change assessment structure. Note: we defined the barycenter moving direction in the same quadrant as the direction movement, and the barycenter moving direction in a different quadrant as the different direction movement.
Sustainability 11 06271 g002
Figure 3. The areas of the rice habitat suitability zones in different scenarios and periods. All changes are relative to the baseline period. M: medium term; L: Long term. AF: Africa; AS: Asia; EU: Europe, NA: North America; OC: Oceania; SA: South America.
Figure 3. The areas of the rice habitat suitability zones in different scenarios and periods. All changes are relative to the baseline period. M: medium term; L: Long term. AF: Africa; AS: Asia; EU: Europe, NA: North America; OC: Oceania; SA: South America.
Sustainability 11 06271 g003
Figure 4. Predicted barycenter shifts of the rice suitability zones on all continents in the future. Baseline-L: the baseline period of the lower-grade suitability zones; Baseline-H: the baseline period of the higher-grade suitability zones. Other items: first field (scenario); second field (period); M (medium term); L (long term); third field (grade); L (lower-grade suitability zone); H (higher-grade suitability zone).
Figure 4. Predicted barycenter shifts of the rice suitability zones on all continents in the future. Baseline-L: the baseline period of the lower-grade suitability zones; Baseline-H: the baseline period of the higher-grade suitability zones. Other items: first field (scenario); second field (period); M (medium term); L (long term); third field (grade); L (lower-grade suitability zone); H (higher-grade suitability zone).
Sustainability 11 06271 g004
Figure 5. High temperature hazard distribution during the baseline period. NS: unsuitable zone; NE: non-exposed zone.
Figure 5. High temperature hazard distribution during the baseline period. NS: unsuitable zone; NE: non-exposed zone.
Sustainability 11 06271 g005
Figure 6. High temperature hazard distribution during the medium-term and long-term periods.
Figure 6. High temperature hazard distribution during the medium-term and long-term periods.
Sustainability 11 06271 g006
Figure 7. Change rates of the high temperature hazards. NC: unchanged zone; N2E: non-exposed zone to exposed zone; E2N: exposed zone to non-exposed zone.
Figure 7. Change rates of the high temperature hazards. NC: unchanged zone; N2E: non-exposed zone to exposed zone; E2N: exposed zone to non-exposed zone.
Sustainability 11 06271 g007
Figure 8. Barycenter shift of the high temperature hazards.
Figure 8. Barycenter shift of the high temperature hazards.
Sustainability 11 06271 g008
Figure 9. Accuracy test of potential rice distribution estimates. Note: NN indicates that there is no potential rice distribution in the non-rice distribution area; NY indicates that there is potential rice distribution in the non-rice distribution area; YN indicates that there is no potential rice distribution in the rice distribution area; and YY indicates that there is potential rice distribution in the rice distribution area.
Figure 9. Accuracy test of potential rice distribution estimates. Note: NN indicates that there is no potential rice distribution in the non-rice distribution area; NY indicates that there is potential rice distribution in the non-rice distribution area; YN indicates that there is no potential rice distribution in the rice distribution area; and YY indicates that there is potential rice distribution in the rice distribution area.
Sustainability 11 06271 g009
Table 1. Datasets.
Table 1. Datasets.
Data CategoryData NameTemporal ResolutionSpatial ResolutionSource
Disaster-Formative Environment TheoryGlobal Digital Elevation Model Data GMTED2010 (DEM)20101 km × 1 kmUnited States Geological Survey (USGS) https://topotools.cr.usgs.gov/gmted_viewer/
Global-Slope Data200610 km × 10 kmInternational Institute for Applied Systems Analysis Global Agro ecological Zones (GAEZ): http://www.gaez.iiasa.ac.at
Global-Soil Parameter Data20121 km × 1 kmInternational Soil Reference and Information Centre (ISRIC): http://www.isric.org
Hazard-Bearing BodyRice Cultivation Range Data2000 or 20050.5° × 0.5°①Harvested Area and Yield for 175 Crops around the year 2000 (EARTHSTAT 2000): http://www.earthstat.org/harvested-area-yield-175-crops/
②Global monthly irrigated and rainfed crop areas around the year 2000 (MIRCA 2000): http://www.unifrankfurt.de/45218031/data_download
③Spatial Production Allocation Model (SPAM) 2000/2005: http://mapspam.info/maps/
ScenariosIPSL-CM5A-LR1971~20990.5° × 0.5°LOCEAN/IPSL: https://www.locean-ipsl.upmc.fr/smos/
MIROC-ESM-CHEM1971~20990.5° × 0.5°Japan Agency for Marine-Earth Science and Technology, Atmosphere and Ocean Research Institute (The University of Tokyo), and National Institute for Environmental Studies: http://adsabs.harvard.edu/abs/2011GMD.....4..845W
NorESM1-M1971~20990.5° × 0.5°Norwegian Earth System: https://portal.enes.org/models/earthsystem-models/ncc/noresm
Land use Data (LUH2)1971~21000.25° × 0.25°Land-Use Harmonization: http://luh.umd.edu/data.shtml
Bioclimatic Variables1971~2000
2041~2060
2061~2080
1 km × 1 kmGlobal Climate Data (WorldClim): http://worldclim.org/CMIP5v1
Table 2. Suitability classification rules.
Table 2. Suitability classification rules.
GradeVariable Valuation Rules
I>99 quantiles, <1 quantile
II97~99 quantiles, 1~3 quantiles
III95~97 quantiles, 3~5 quantiles
IV90~95 quantiles, 5~10 quantiles
V<90 quantiles, >10 quantiles
Table 3. The degree of change in high temperature exposure in rice.
Table 3. The degree of change in high temperature exposure in rice.
ContinentsSYDRCP2.6 and SSP1RCP4.5 and SSP2RCP6.0 and SSP4RCP8.5 and SSP5Avg
B-M-LB-M-LB-M-LB-M-LB-M-L
AFHighIIIIIIIIIIIIIII
LowIIIIIIIII
ASHighIIIIIIIII
LowIIIIIIIIII
EUHighIVVIVIIIIV
LowIIIVIVIVV
NAHighVVVIIIIV
LowIVIIIIIIIIIIII
OCHighIIIIIIIIIIIIIV
LowIIIIIIIIVIII
SAHighVVIVVV
LowIIIVVVV
Note: I–V represents the severity of exposure changes (from aggravated to relief), or the level of potential risk (from high to low): I indicates significant exacerbation; II indicates general exacerbation; III indicates no significant change; IV indicates general reduction; V indicates a significant reduction. SYD represents the suitability level: High indicates V and IV regions, and Low indicates III and II regions. B represents baseline; M represents medium term; L represents long term.
Table 4. The degree of change in high temperature exposure in rice.
Table 4. The degree of change in high temperature exposure in rice.
Accuracy TypeData for Verification
EARTHSTAT 2000 Harvest DataSPAM 2000 Harvest DataSPAM 2000 Sowing DataMIRCA 2000 Harvest DataSPAM 2005 Harvest DataSPAM 2005 Sowing Data
YY3.6%2.6%2.5%1.9%2.8%2.9%
NN91.6%92.6%92.5%93.0%91.2%91.4%
NY2.7%3.8%3.9%4.5%4.7%4.6%
YN2.1%1.0%1.1%0.6%1.3%1.1%
User Accuracy52.8%35.4%36.1%23.6%32.6%33.3%
Producer Accuracy80.8%83.8%84.2%80.4%81.3%81.8%
All Accuracy93.9%93.8%94.2%92.9%92.4%92.8%
Kappa Coefficient0.610.470.480.350.430.44

Share and Cite

MDPI and ACS Style

Wang, R.; Jiang, Y.; Su, P.; Wang, J. Global Spatial Distributions of and Trends in Rice Exposure to High Temperature. Sustainability 2019, 11, 6271. https://0-doi-org.brum.beds.ac.uk/10.3390/su11226271

AMA Style

Wang R, Jiang Y, Su P, Wang J. Global Spatial Distributions of and Trends in Rice Exposure to High Temperature. Sustainability. 2019; 11(22):6271. https://0-doi-org.brum.beds.ac.uk/10.3390/su11226271

Chicago/Turabian Style

Wang, Ran, Yao Jiang, Peng Su, and Jing’ai Wang. 2019. "Global Spatial Distributions of and Trends in Rice Exposure to High Temperature" Sustainability 11, no. 22: 6271. https://0-doi-org.brum.beds.ac.uk/10.3390/su11226271

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop