Journal of Integrative Agriculture 2018, 17(0): Available online at ScienceDirect

Transcription

1 Journal of Integrative Agriculture 2018, 17(0): Available online at ScienceDirect RESEARCH ARTICLE Regionalization of wheat powdery mildew oversummering in China based on digital elevation ZOU Ya-fei 1, QIAO Hong-bo 2, CAO Xue-ren 3, Liu Wei 1, FAN Jie-ru 1, SONG Yu-li 4, WANG Bao-tong 5, ZHOU Yi-lin 1 1 State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing , P.R.China 2 College of Information and Management Science, Henan Agricultural University, Zhengzhou , P.R.China 3 Environment and Plant Protection Institute, China Academy of Tropic Agricultural Sciences/Key Laboratory of Integrated Pest Management on Tropical Crops, Ministry of Agriculture, Haikou , P.R.China 4 Institute of Plant Protection, Henan Academy of Agricultural Sciences, Zhengzhou , P.R.China 5 College of Plant Protection, Northwest A&F University, Yangling , P.R.China Abstract Blumeria graminis f. sp. tritici, the pathogen that causes wheat powdery mildew, is one of the most important diseases affecting wheat production in China, and the oversummering is the key stage of wheat powdery mildew epidemic. The more oversummering regionalization of wheat powdery mildew has played an important role in disease prediction, prevention and control. In this study, we analyzed the correlation between oversummering data of wheat powdery mildew and the meteorological factors over the past years, and determined that temperature was the key meteorological factor influencing oversummering of wheat powdery mildew. The average temperature at which wheat powdery mildew growth was terminated (26.2 C) was used as the threshold temperature to regionalize the oversummering range of wheat powdery mildew. This regionalization was done using the GIS ordinary kriging method combined with the Digital Elevation model (DEM) of China. The results showed that annual probability of oversummering region based on Model 26.2 were consistent with the actual survey of the more summer wheat powdery mildew. Wheat powdery mildew oversummering regions in China mainly cover mountainous or high-altitude areas, and these regions form a narrow north-south oversummering zone. Oversummering regions of wheat powdery mildew is mainly concentrated in the high-altitude wheat growing areas, including northern and southern Yunnan, northwestern Guizhou, northern and southern Sichuan, northern and southern Chongqing, eastern and southern Gansu, southeastern Ningxia, northern and southern Shaanxi, central Shanxi, western Hubei, western Henan, northern and western Hebei, western Liaoning, eastern Tibet, eastern Qinghai, western Xinjiang and other regions of China. Keywords: wheat powdery mildew, digital elevation model, oversummering regionalization, geographic information system Received 4 September, 2017 Accepted 29 November, 2017 ZOU Ya-fei, zouyafei@caas.cn; Correspondence ZHOU Yi-lin, Tel: , Fax: , ylzhou@ippcaas.cn 2018 CAAS. Publishing services by Elsevier B.V. All rights reserved. doi: /S (17) Introduction Wheat powdery mildew, caused by Blumeria graminis f. sp. tritici, is one of the most destructive foliar diseases of wheat world-wide, especially in areas with maritime or

2 *** et al. Journal of Integrative Agriculture 2018, 17(0): semi-continental climates (Bennett 1984). In China, with the deployment of uniform resistance cultivars with a semidwarf growth habit, improved irrigation and fertilization conditions and high-yielding traits, wheat powdery mildew epidemics have been severe since the late 1970s (Liu and Shao 1995). The most destructive epidemics occurred in 1990 and 1991, causing yield losses up to 1.44 and 0.77 million tons, respectively (Liu and Shao 1994). The disease cycle of wheat powdery mildew includes oversummering, autumn seedling infection, overwintering and spring epidemic. And oversummering is the key stage of all epidemic process. It has been shown that the pathogen can oversummer as mycelium and conidia on volunteer wheat seedlings or summer-sown wheat. Therefore, the oversummering range of this disease and the amount of pathogen are closely related with its prevalence in autumn and spring in China (Shi and Wan 1983; Li 1987). Hence, an oversummering regionalization study for this disease is of importance in its prediction, prevention and controlling. Meteorological conditions have significant effects on the prevalence rate and damage degree of wheat powdery mildew, and especially on the oversummering of the pathogen (Li 1986). A large number of studies have shown that the average temperature in the hottest continous 10-day period is the key factor influencing the oversummering of wheat powdery mildew (Lei et al. 1982; Wang et al. 1984; Wu et al. 1984; Liu et al. 1990; Duan et al. 1992; Liu et al. 1992; He et al. 1998). Researchers generally believe that temperature is the key factor influencing oversummering, but this has not been confirmed by large-scale data analysis. In a previous study, Li et al. (2013) classified regions of wheat powdery mildew oversummering in China using the temperatures of the hottest continuous 10-day period below the termination threshold (25.5 C) (model 1) and the temperatures equivalent to those terminating wheat powdery mildew growth (model 2) as major regionalization indexes. However, this study did not account for the effect of digital elevation on spatial interpolation, and the lack of spatial interpolation points resulted in low resolution maps of oversummering which possibly led to inaccurate in the oversummering regionalization of wheat powdery mildew. Another is that in the context of global warming, the temperature sensitivity of B. graminis f. sp. tritici was showing a downward trend (Wan et al. 2010; Ming et al. 2012; Caidanzhuoma et al. 2015), and thus the average termination threshold of the disease was rising accordingly, leading to changes in the oversummering range of wheat powdery mildew in China. Hence, it is necessary to further determine the oversummering range of wheat powdery mildew in China. The objectives of this study are to: 1) analyze the historical meteorological factor data and long-term oversummering survey data to pinpoint the key meteorological factor affecting oversummering of wheat powdery mildew; 2) accurately regionalize the oversummering range of wheat powdery mildew based on digital elevation and a proper threshold temperature of wheat powdery mildew. 2. Materials and methods 2.1. Data and resources Weather data Continuous data for daily average temperature, daily maximum temperature, daily average precipitation, daily average relative humidity, and sunshine duration at 752 weather stations in China (Fig. 1) from 1951 to 2013 were obtained from the National Meteorological Center of China. Sources of spatial data A map of China showing province information (1: ) was downloaded from the National Fundamental Geographic Information System of China ( The Digital Elevation model (DEM) was provided by the National Administration of Surveying, Mapping and Geoinformation of China. The resolution of the DEM was 500 m, and the elevation ranged from 159 to m. A map of wheat-growing regions in China was also scanned from the Chinese Agricultural Atlas (CCNA 1989) and digitalized by using ArcGIS. Disease data for oversummering of wheat powdery mildew Survey data of the percentage of diseased plants for oversummering of wheat powdery mildew were collected from two sources: 1) 194 survey sites from 45 counties (cities) in Anhui, Beijing, Gansu, Guizhou, Hebei, Henan, Shaanxi and Yunnan of China from 2004 to 2009; and 2) previously published oversummering survey data for 68 survey sites (Lei et al. 1982; Zhang 1983; Wang 1984; Wu et al. 1984; Li 1986; Li 1987; Liu et al. 1990; Duan et al. 1992; Liu et al. 1992; Zhang et al. 1993; Tu et al. 1999). Of the total of 262 survey sites (Fig. 2), 158 can oversummer. The occurrence of wheat powdery mildew was converted to ranked data (hereinafter referred to as oversummering data). That is, the presence of wheat powdery mildew in a surveyed site was denoted by 1, and absence was denoted by Analytical method Correlation analysis methods (1) Selection of meteorological factor data. Average temperature was mean of the hottest 10 days in July and August in , and maximum temperature was the hottest 1 day in July and August in Average precipitation, average relative humidity and sunshine duration were mean of July and August in The selected meteorological factors were shown in Table 1.

3 4 *** et al. Journal of Integrative Agriculture 2018, 17(0): N Weather station km Fig. 1 Geographical location of weather stations throughout China used in the study. N Oversummering survey sites km Fig. 2 Geographical location of the oversummering survey sites used in this study. (2) Regression analyses of temperature vs. longitude/ latitude/altitude. The regression test was carried out by fitting the regression equation with the mean temperature, latitude, longitude or altitude of 752 sites by years. Then

4 *** et al. Journal of Integrative Agriculture 2018, 17(0): Table 1 The selected meteorological factors Meteorological factor Data selection Average temperature Mean of the hottest 10 days in July and August in Maximum temperature Temperature of the hottest 1 day in July and August in Average precipitation Mean of July and August in Average relative humidity Mean of July and August in Sunshine duration Mean of July and August in the average temperature of the hottest 10 days of the surveyed sites was predicted. The same method was used to fit regression equations between the mean temperature of the hottest 1 day and latitude, longitude or altitude by year. In this study, multiple linear regression models were used for analysis. (3) Correlation analyses of oversummering data and meteorological factors. Analytical data included oversummering data, the average temperature, maximum temperature, average precipitation, average relative humidity and sunshine duration of the 262 oversummering surveyed sites. Firstly, the average and maximum temperatures of the 262 surveyed sites during the 10 hottest days in July and August were calculated respectively using the average and maximum temperatures of linear regression equations. Then data for meteorological factors (average precipitation, average relative humidity and sunshine duration) were extracted from 262 oversummering survey sites using the GIS spatial interpolation method (Li et al. 2013). Spearman s rank-order correlation (Fang et al. 2013) was used to analyze the correlations between oversummering data and meteorological factors (average temperature, maximum temperature, average precipitation, average relative humidity and sunshine duration). Oversummering regionalization of wheat powdery mildew (1) Selection of a threshold temperature. Correlation analysis revealed that temperature is a key factor influencing the oversummering of wheat powdery mildew. Hence, it is crucial to select an appropriate threshold temperature to use in oversummering regionalization. Our research team has monitored sensitivity of population of B. graminis f. sp. tritici to temperature from different regions in China for recent years. The termination threshold of the temperature of 113 single colony isolates of B. graminis f. sp. tritici collected from different regions in China in 2008 were ranged from 25.2 to 29.3 C, with an average termination threshold of 26.2 C (Wan et al. 2010). For 129 single colony isolates of B. graminis f. sp. tritici collected from different regions in China in 2009, the termination threshold of the temperature ranged from 25.6 to 28.7 C with an average termination threshold of 26.7 C (Ming et al. 2012). For 139 single spore isolates of B. graminis f. sp. tritici collected from different regions in China in 2012, the termination threshold of the temperature ranged from 25.4 to 28.7 C, with an average termination threshold of 26.2 C (Caidanzhuoma et al. 2015). These results showed average thresholds ranging from 26.2 to 26.7 C in the temperature sensitivity of B. graminis f. sp. tritici from different regions in different years. This range is higher compared with the previously reported termination threshold range (23.5 to 26.0 C) (Lei et al. 1982; Wang 1984; Wu et al. 1984; Liu et al. 1990; Duan et al. 1992; Liu et al. 1992; He et al. 1998), indicating decreased temperature sensitivity. In order to more accurately reflect the temperature sensitivity level of B. graminis f. sp. tritici in China, the median of these average termination thresholds of 2008, 2009 and 2012 (26.2 C) was selected as the threshold temperature for oversummering regionalization of wheat powdery mildew. Tests of temperature termination threshold of isolates of B. graminis f. sp. tritici were carried out in laboratory. The single colony isolates of B. graminis f. sp. tritici collected from different regions in China were tested by using detached leaf segment method at 6 different temperatures which were 18, 22, 23, 24, 25 and 26 C. On the 12th day after inoculation, the area of disease was recorded and the disease severity and inhibition rate were calculated. The equations of disease inhibition rate (Y) against temperature (X) for each isolate were constructed by linear regression analyses: Y=a+bX, when Y=100% calculated value of temperature of each equation was temperature termination threshold (ET 100 ) of individual isolate of B. graminis f. sp. tritici (Han et al. 2016). (2) Oversummering regionalization method. Using the Create Fishnet tool in ArcGIS, a fishnet with 20 km 20 km grids was created using a digital map of China as the base map, and data for longitude, latitude and altitude were extracted from all grid points, with available data points obtained. The average temperatures for each years at each of the points were calculated using the established regression equation (Table 2) and compared with the average termination thresholds (26.2 C). The annual probabilities of wheat powdery mildew oversummering at these points were then calculated as follows: Percent annual oversummering probability (%)=(The number of years when the average temperature below 26.2 C/The number of years from 1951 to 2013) 100. Based on this probability, the ordinary Kriging method (Li et al. 2013) was used for spatial interpolation and mapping oversummering regions of average termination thresholds of 26.2 C (hereinafter referred to as Model 26.2). (3) Comparison of oversummering regionalization results. Receiver operating characteristic (ROC) curves was used to determine the accuracy of oversummering regionalization maps based on the average termination thresholds of 26.2 C

5 6 *** et al. Journal of Integrative Agriculture 2018, 17(0): and the oversummering probability threshold of wheat powdery mildew. Among the 252 sites surveyed, B. graminis f. sp. tritici found to have produced disease symptoms on volunteer plants at 158 sites. Sensitivity refers to the proportion of data sets where disease can oversummer that were classified correctly, and specificity refers to the percentage of data sets where disease could not oversummer that were classified correctly. To construct ROC curves, sensitivity was plotted as the ordinate and 1 specificity (referred to as the false positive proportion) was plotted as the abscissa. A larger area under the ROC curve (AUROC) indicates higher regionalization accuracy. For each model, the cut-point where the overall error rate was the smallest was determined by selecting the point with the highest value of Youden s index, J, which is the point on the ROC curve at the greatest geometric distance from the line representing a noninformative predictor (Metz 1978). Youden s index is commonly used as a measure of overall prediction effectiveness and is calculated by J=Sensitivity+Specificity 1 (Metz 1978). The coincidence rate (%) between the model predictions and observed data for the 158 sites was then calculated as follows: Coincidence frequency (%)=[Number of oversummering sites (predicted by model)]/[number of oversummering sites (wheat powdery mildew observed)] Results 3.1. Relationship between temperature, latitude, longitude or elevation Through regression analysis, the equations for the relationship between the average of the hottest 10 days and altitude values, longitude values, latitude values during were established for each year, respectively (Table 2), and the regression equation are different in different years Correlation analysis of oversummering data and meteorological factors Spearman s rank-order correlation analysis was conducted to determine the relationship between oversummering data and meteorological factors (average temperature, maximum temperature, average precipitation, average relative humidity and sunshine duration). Based on these results (Table 3), there was only a significant correlation between oversummering and average and maximum temperature (all P<0.01) suggested that temperature is the key meteorological factor affecting oversummering of wheat powdery mildew. Table 2 Regression equations for the relationship between the average of the hottest 10 days between July and August and altitude values, longitude values, latitude values during Year Regression equation 1) r Y= X ** 2012 Y= X ** 2011 Y= X X ** 2010 Y= X ** 2009 Y= X ** 2008 Y= X ** 2007 Y= X ** 2006 Y= X X ** 2005 Y= X ** 2004 Y= X ** 2003 Y= X ** 2002 Y= X X ** 2001 Y= X ** 2000 Y= ** 1999 Y= ** 1998 Y= X ** 1997 Y= X ** 1996 Y= X X ** 1995 Y= X X ** 1994 Y= X ** 1993 Y= X ** 1992 Y= X ** 1991 Y= X ** 1990 Y= X ** 1989 Y= X ** 1988 Y= X ** 1987 Y= X X ** 1986 Y= X X ** 1985 Y= X X ** 1984 Y= X ** 1983 Y= X X ** 1982 Y= X ** 1981 Y= X ** 1980 Y= X X ** 1979 Y= X X ** 1978 Y= X ** 1977 Y= X ** 1976 Y= X ** 1975 Y= X X ** 1974 Y= X X ** 1973 Y= X X ** 1972 Y= X X ** 1971 Y= X ** 1970 Y= X ** 1969 Y= X ** 1968 Y= X ** 1967 Y= X X ** 1966 Y= X ** 1965 Y= X X ** 1964 Y= X X ** 1963 Y= X ** 1962 Y= X X ** 1961 Y= X ** 1960 Y= X ** (Continued on next page)

6 *** et al. Journal of Integrative Agriculture 2018, 17(0): Table 2 (Continued from preceding page) Year Regression equation 1) r Y= X X ** 1958 Y= X ** 1957 Y= X X ** 1956 Y= X X X ** 1955 Y= X X X ** 1954 Y= X X ** 1953 Y= X X ** 1952 Y= X X ** 1951 Y= X ** 1) Y stands for temperature; X 1 stands for latitude; X 2 stands for longitude; X 3 stands for altitude. **, correlation is significant at the 0.01 level by linear regression analysis Wheat powdery mildew oversummering regionalization Oversummering regionalization based on digital elevation and average termination threshold (26.2 C) The ordinary kriging model was used for spatial interpolation of annual oversummering probability. Using wheat-growing regions as a base map, a map of oversummering regionalization based on the average termination threshold 26.2 C of wheat powdery mildew (Model 26.2) (Fig. 3). In this map, wheat powdery mildew oversummering regions in China are located mainly in alpine or high altitude areas, forming a narrow north-south main oversummering zone. Regions with an annual oversummering probability of greater than 50% mainly cover northern and southern Yunnan, northwestern Guizhou, northern and southern Sichuan, northern and southern Chongqing, eastern and southern Gansu, southeastern Ningxia, northern and southern Shaanxi, central Shanxi, western Hubei, western Henan, northern and western Hebei, western Liaoning, eastern Tibet, eastern Qinghai, western Xinjiang and other regions. These regions are high suitable areas for oversummering of wheat powdery mildew. Regions with an annual oversummering probability of between 15 and 50% mainly cover the periphery of the Sichuan Basin, central Guizhou, central Shaanxi, southern Shanxi, central Shandong, the Liaodong Peninsula, and some mountainous areas in Hunan, Fujian and Zhejiang. These regions are at moderate suitable areas for wheat powdery mildew oversummering. Regions with an annual oversummering probability of between 5% and 15% mainly cover central Hebei, southern Beijing, Tianjin, and eastern Shandong. These regions are at low suitable areas for wheat powdery mildew oversummering. Oversummering regionalization based on digital elevation and average termination threshold (25.5 C) Li et al. (2013) made maps of oversummering regionalization based on the average termination threshold 25.5 C of wheat powdery mildew in China. In order to facilitate comparison, this study uses the same method as Model 26.2 to made a oversummering regionalization map with an average temperature termination threshold of 25.5 C (hereinafter referred to as Model 25.5) (Fig. 4). In this map, oversummering regions for wheat powdery mildew also form a narrow north-south main oversummering zone. However, unlike Model 26.2, the range of wheat powdery mildew slightly narrows, and oversummering probabilities in some regions also decline to different degrees. Regions with a great difference in oversummering probability determined using Model 25.5 and Model 26.2 mainly cover southern Yunnan, central Sichuan, central Guizhou, southeastern Chongqing, western Hubei, central Shaanxi, western Henan, southern Shanxi, mountainous areas in western Hebei, and western Liaoning. Comparison of oversummering regionalization Model Table 3 The results of correlation analysis of oversummering data of wheat powdery mildew and meteorological factors 1) Factor r and Average Maximum Average relative Average Sunshine Oversummering P 1) temperature temperature humidity precipitation duration data Average temperature r P Maximum temperature r ** P <0.001 Average relative humidity r ** * P Average precipitation r ** P <0.001 Sunshine duration r ** ** P <0.001 <0.001 Oversummering data r ** ** P < ) r stands for correlation coefficient; P stands for significance. *, correlation is significant at the 0.05 level; **, correlation is significant at the 0.01 level.

7 8 *** et al. Journal of Integrative Agriculture 2018, 17(0): N Pathogen found Annual probability (%) km Fig. 3 Regionalization of the annual probability of wheat powdery mildew oversummering based on the average termination threshold 26.2 C (Model 26.2). N Pathogen found Annual probability (%) km Fig. 4 Regionalization of the annual probability of wheat powdery mildew oversummering based on the average termination threshold 25.5 C (Model 25.5) and Model 26.2 maps The AUROC values for Model 26.2 and Model 25.5 were and , respectively (Fig. 5), indicating that Model 26.2 had higher regionalization accuracy than Model In Model 25.5,

8 *** et al. Journal of Integrative Agriculture 2018, 17(0): Youden s index is optimal for a cut-point of 15%, and the corresponding sensitivity and specificity values are 90.4 and 60.4%, respectively. we selected the cut-point in Model 25.5, 15%, to define the threshold for wheat powdery mildew oversummering. That is to say, if the local wheat powdery mildew oversummering probability reaches 15% and above, then this disease is considered able to oversummer. As shown in Table 4, of the 158 survey sites where wheat powdery mildew could oversummer, there were 144 and 152 oversummering sites located in regions with an oversummering probability of over 15% based on Model 25.5 and Model 26.2, respectively. Because Model 26.2 had a higher coincidence rate (96.20%) than Model 25.5 (91.14%), it can be concluded that predictions of oversummering based on Model 26.2 are more coincident with actual surveys, and more accurate regionalization results can be attained using this model. Sensitivity Specificity Fig. 5 Receiver operating characteristic plots for the two models, 25.5 (filled circles) and 26.2 (open circles), used to calculate the annual probability (%) of oversummering. 4. Discussion and conclusion In this study, we pinpointed temperature as the most important meteorological factor determining the probability of wheat powdery mildew oversummering. Based on this finding, we used two different termination thresholds, i.e., Model 26.2 and Model 25.5 to generate maps of wheat powdery mildew oversummering regionalization that also took into account digital elevation data. Model 26.2 had higher regionalization accuracy than Model 25.5 and was also more coincident with actual oversummering survey results. Thus Model 26.2 is a better model for generating maps of wheat powdery mildew oversummering regionalization. Generally, it takes approximately 10 days from infection of wheat by B. graminis f. sp. tritici to full development of symptoms under appropriate conditions. In the presence of the host, the fungi are able to live through the hottest 10 days in July and August of each year, i.e., to oversummer in volunteer wheat seedlings or summer-sown wheat. Therefore, average temperature was selected mean of the hottest 10 days in July and August in , and maximum temperature was the hottest 1 day in July and August in We also considered the effects of average daily precipitation, temperature and sunshine on the growth of volunteer wheat seedlings or summer-sown wheat and wheat powdery mildew in the oversummering period of disease. Accordingly, average precipitation, average relative humidity and sunshine duration were selected mean of July and August in The temperature distribution has a certain correlation with latitude, altitude or longitude in China (Fang 1992). The temperature of the survey sites can be predicted by fitting the regression relationship between temperature and latitude, longitude or altitude. Li et al. (2013) was fitted with the regression relationship between the highest temperature of 10 days in Table 4 Number of sites that Blumeria graminis f. sp. tritici can oversummer by survey and number of these sites that can oversummer using the two models used in this study in each province or city Province or City Oversummering Predicted oversummering Coincidence rate Predicted oversummering Coincidence rate site site (Model 25.5) (%) 1) site (Model 26.2) (%) Henan Shaanxi Gansu Yunnan Guizhou Sichuan Hubei Shandong Beijing Total ) Coincidence rate (%)=Number of oversummering sites (predicted bymodel)]/number of oversummering sites (powdery mildew observed) 100

9 10 *** et al. Journal of Integrative Agriculture 2018, 17(0): different sites, the elevation, latitude and longitude, so as to predict the highest temperature of 10 days of survey sites of wheat powdery mildew. However, this method ignored the difference between the annual average temperature of the hottest 10 days, and the prediction results may be inaccurate. Therefore, this study took full account of the differences in the mean temperature of the hottest 10 days of the same survey in different years. The regression test was carried out by fitting the regression equation with the mean temperature, latitude, longitude or altitude of 752 sites by years. At the same time, this study combined a high-resolution digital elevation model, the interpolation resolution was greatly improved. From the test results of oversummering sites, this study added additional 38 new oversummering sites on the basis of 121 ones selected by Li (158 sites in total). Using the temperatures of the hottest continuous 10-day period below the termination threshold (25.5 C) (model 1) and the temperatures equivalent to those terminating wheat powdery mildew growth (model 2), Li et al. (2013) obtained overall coincidence rates of and 90.91%, respectively, for oversummering regionalization. However, in this study, the coincidence rate for oversummering regionalization (Model 25.5) was 91.14% (Table 5), suggested that the method used for wheat powdery mildew oversummering regionalization in this study is more reasonable and yields more reliable results. An IPCC report predicting the future climatic conditions indicates that the global temperature will rise at an average rate of 0.3 C per decade (the predicted range is between C) (IPCC 1990). In China, the average yearly temperature has increased by 1.2 C since 1960 (Piao et al. 2010). Differences in the temperature sensitivity of B. graminis f. sp. tritici isolates collected from different provinces or cities in China in recent years and the frequency distribution of B. graminis f. sp. tritici isolates with different sensitivities to temperature do not fit a normal distribution, which indicates that global climate change may influence the evolution of B. graminis f. sp. tritici (Wan et al. 2010). Global warming and adaptive changes of B. graminis f. sp. tritici may therefore have an effect on the oversummering and overwintering of wheat powdery mildew and even on Table 5 Validation of models for predicting wheat powdery mildew oversummering in China Model 1 Model 2 Model 25.5 Oversummering site (n) Predicted oversummering site (n) Coincidence rate (%) 1) ) Coincidence rate (%)=Number of oversummering sites (predicted by model)/number of oversummering sites (wheat powdery mildew observed) 100 the whole epidemic process of this disease. Therefore, the relationship between parasitism fitness and temperature of different temperature-sensitive of B. graminis f. sp. tritici isolates should be deeply carried out, and the interaction model between different sensitive strains should be established to predict the trend of future pathogens, oversummering, overwintering and epidemic situation in climate change scenarios. 5. Conclusion In this study, we pinpointed temperature as the most important meteorological factor determining the probability of wheat powdery mildew oversummering. Based on this finding, we used two different termination thresholds i.e., Model 26.2 and Model 25.5 to generate maps of wheat powdery mildew oversummering regionalization that also took into account digital elevation data. Model 26.2 had higher regionalization accuracy than Model 25.5 and was also more coincident with actual oversummering survey results. Thus Model 26.2 is a better model for generating maps of wheat powdery mildew oversummering regionalization. According to the map Model 26.24, the main oversummering zone in China are mainly concentrated in northern and southern Yunnan, northwestern Guizhou, northern and southern Sichuan, northern and southern Chongqing, eastern and southern Gansu, southeastern Ningxia, northern and southern Shaanxi, central Shanxi, western Hubei, western Henan, northern and western Hebei, western Liaoning, eastern Tibet, eastern Qinghai and western Xinjiang and other regions. Acknowledgements This work was financially supported by the National Natural Science Foundation of China ( ), the National Key Research and Development Program of China (2016YFD ), the Special Fund for Agro-scientific Research in the Public Interest, China ( ). Thanks to the platform support of Key Laboratory of Integrated Pest Management on Crops, Ministry of Agriculture, P.R.China. References Bennett F G A Resistance to powdery mildew in wheat: a review of its use in agriculture and breeding programmes. Plant Pathology, 33, Caidanzhuoma, Gang C W, Fang J R, Zhou Y L, Duan X Y, Zou Y F, Wan Q, Na T C The sensitivity of population of Blumeria graminis f. sp. tritici isolates to temperature in Plant Protection, 41, (in Chinese) CCNA (The Compilation Committee of National Atlas)

10 *** et al. Journal of Integrative Agriculture 2018, 17(0): National Agricultural Atlas of the People s Republic of China. China SinoMaps Press, Beijing. (in Chinese) Duan S K, Wang B T, Liu C H Status and analysis of wheat powdery mildew in Guanzhong area of Shaanxi Province. Journal of Triticeae Crops, 2, (in Chinese) Fang J Q, Xu Y Y, Chen F Health Statistics. 7th ed. People s Medical Publishing House, Beijing. (in Chinese) Fang J Y Study on the geografic elements affecting temperature distribution in China. Acta Ecologica Sinca, 12, (in Chinese) Han L P, Fan J R,Zhou Y L, Zou Y F, Duan X J, Li N J, Zhu Y F Temperature sensitivity and its correlation with virulence of Blumeria graminis f. sp. tritici and isolates in some wheat growing areas in Acta Phytopathologica Sinica, 46, (in Chinese) He J B, Song Y L, Zhang Z S, He W L Wheat powdery mildew and control. I. Distribution, symptom and harm of wheat powdery mildew. Journal of Henan Agricultural Sciences, 1, (in Chinese) IPCC (Intergovernmental Panel on Climate Change) Climate Change. Cambridge Press, London. Lei T W, Feng J Y, Duan J T, Guo Y W, Cai C L, Fu S H Preliminary study on occurrence regularity of wheat powdery mildew. Journal of Henan Institute of Science and Technology, 2, (in Chinese) Li B N Zones of Oversummering and Overwintering of Wheat Powdery Mildew Based on GIS. Graduate School of Chinese Academy of Agricultural Sciences, Beijing.(in Chinese) Li B N, Cao X R, Chen L, Zhou Y L, Duan X Y, Luo Y, Fitt B D L, Xu X, Song Y, Wang B, Cao S Application of geographic information systems to identify the oversummering regions of Blumeria graminis f. sp. tritici in China. Plant Disease, 97, Li H Infection cycle of wheat powdery mildew in Sichuan Basin. Journal of Southwest University, 2, (in Chinese) Li J P Occurrence and control of wheat powdery mildew in Shangluo. Shaanxi Journal of Agricultural Sciences, 2, (in Chinese) Li T X Research progress of wheat powdery mildew forecasting method. Meteorological and Environmental Sciences, 36, (in Chinese) Liu Q J, Liu Q J, Tang Y L, Zhang B S Epidemiology and control of wheat powdery mildew. Shandong Agricultural Sciences, 6, (in Chinese) Liu W C, Shao Z R Epidemiology, occurrence and analysis of wheat powdery mildew in recent years. Plant Protection Technology and Extension, 6, (in Chinese) Liu W C, Shao Z R Analysis of the causes and trends of wheat powdery mildew in China.Plant Protection Technology and Extension, 6, 7 8. (in Chinese) Liu X K, Zhou Y L, Liang G D Investigation on the oversummering of asexual generation of Erysiphe graminis f. sp. tritici in Nanyang. Plant Protection, 1, 33. (in Chinese) Metz C E Basic principles of ROC analysis. Nuclear Medicine, 8, Ming Z Y, Wang B T, Zhou Y L, Duan X Y, Zou Y F The sensitivity of population of Blumeria graminis f. sp. tritici isolates to temperature in Acta Phytopathologica Sinica, 42, (in Chinese) Piao S, Giais P, Huang Y, Shen Z, Peng S, Li J, Zhou L, Liu H, Ma Y, Ding Y, Friedlingstein P, Liu C, Tan K, Ye Y, Zhang T, Fang J The impacts of climate change on water resources and agriculture in China. Nature, 467, Shi P G, Wan J H Research overview of Blumeria graminis f. sp. tritici. Hubei Agricultural Sciences, 1, (in Chinese) Tu J H, Li Y, Liao H M, Xu X C, Chen B Study on epidemic regularity of wheat powdery mildew in Sichuan. Southwest China Journal of Agricultural Sciences, 12, (in Chinese) Wan Q, Ding K J, Duan X Y, Zhou Y L Sensitivity of population of Blumeria graminis f. sp. tritici isolates to temperature in Acta Phytopathologica Sinica, 40, (in Chinese) Wang S T Study on the inoculum of primary infection of wheat powdery mildew. Hubei Agricultural Sciences, 11, (in Chinese) Wu M Z, Yan S H, Hu G Q Oversummering and overwintering of wheat powdery mildew and study on its biological characteristics. Guizhou Agricultural Sciences, 3, (in Chinese) Zhang Z D, Li Z Q, Zong Z F The laws of occurrence of powdery mildew of wheat in central area of Shaanxi Province. Journal of Northwest Sci-Tech University of Agriculture and Forestry, 21, (in Chinese) Zhang Z H Preliminary study on the occurrence regularity and control of wheat powdery mildew in Zibo. Zibo Agricultural Science and Technology, 1, (in Chinese) Section editor WAN Fang-hao Managing editor ZHANG Juan