Next Article in Journal
A Multiscale Framework for Sustainable Management of Tufa-Forming Watercourses: A Case Study of National Park “Krka”, Croatia
Previous Article in Journal
Temporal and Spatial Variations of the Biochemical Composition of Phytoplankton and Potential Food Material (FM) in Jaran Bay, South Korea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Responses of Winter Wheat Yield to Drought in the North China Plain: Spatial–Temporal Patterns and Climatic Drivers

1
State Key Laboratory of Remote Sensing Science, Beijing Normal University, Beijing 100875, China
2
Beijing Key Laboratory of Environmental Remote Sensing and Digital City, Beijing Normal University, Beijing 100875, China
3
Faculty of Geographical Science, Beijing Normal University, Beijing 100875, China
4
Center for Drought and Risk Research, Beijing Normal University, Beijing 100875, China
5
Institute of Digital Agriculture, Zhejiang Academy of Agricultural Sciences, Hangzhou 310021, China
*
Author to whom correspondence should be addressed.
Water 2020, 12(11), 3094; https://doi.org/10.3390/w12113094
Submission received: 13 September 2020 / Revised: 24 October 2020 / Accepted: 2 November 2020 / Published: 4 November 2020
(This article belongs to the Section Water, Agriculture and Aquaculture)

Abstract

:
Understanding the winter wheat yield responses to drought are the keys to minimizing drought-related winter wheat yield losses under climate change. The research goal of our study is to explore the response patterns of winter wheat yield to drought in the North China Plain (NCP) and then further to study which climatic factors drive the response patterns. For this purpose, winter wheat yield was simulated by the Environmental Policy Integrated Climate (EPIC) crop model. Drought was quantified by standardized precipitation evapotranspiration index (SPEI), and the contributions of the various climatic factors were evaluated using predictive discriminant analysis (PDA) method. The results showed that the responses of winter wheat yield to different time-scale droughts have obvious spatial differences from the north part to the south part in the NCP. Winter wheat yield is more sensitive to the medium (6–9 months) and long (9–12 months) time-scale droughts that occurred in the key growth periods (April and May). The different response patterns of winter wheat yield to the different time-scale droughts are mainly controlled by temperature and water balance (precipitation minus potential evapotranspiration) in winter in the NCP. Compared with the water balance, temperature plays a more important role in driving the response pattern characteristics. These findings can provide a reference on how to reduce drought influences on winter wheat yield in the NCP.

1. Introduction

Drought usually leads to serious agricultural damage and the negative effects of drought cannot be ignored. It is considered as the main meteorological disaster influencing crop production and food security in the world [1,2]. The agriculture of China is likely to be faced with higher drought risk because of continuous climate change [3]. Based on the record of the Chinese Academy of Agricultural Sciences (CAAS), the mean annual crop production losses in China caused by drought were up to 75.7 billion yuan from 1988 to 2004. The grain losses are up to 10 billion kg every year because of drought [4]. Understanding the influence of drought on agriculture is central to minimizing drought-related yield losses and can help farmers reduce the adverse effects of drought on agriculture [5].
As a complex natural disaster, it is not easy to establish a unitary multidisciplinary definition of drought [6,7,8]. For the purpose of reflecting the spatial–temporal distribution characteristics of drought, many drought indexes have been put forward, such as the standardized precipitation index (SPI), Palmer drought severity index (PDSI), and the standardized precipitation evapotranspiration index (SPEI). The SPEI has been widely used in drought monitoring [9,10,11]. Compared with the SPI [12], the SPEI is especially suitable for evaluating the drought influences against the backdrop of global warming [13], since SPEI takes into account the influences of the temperature. PDSI has also been widely used in drought monitoring and assessment, but the major limitation of PDSI is its incomparable parameters in different regions [14]. It should be noted that each drought index has its strengths and weaknesses. However, by comparing the performance of the above three drought indices, it also has been proven that the SPEI performed better than SPI and PDSI in estimating the drought influences on crop production in the NCP [15,16]. Additionally, by studying the drought effect on winter wheat yield of the region, Liu et al. [17] pointed out that the SPEI could explain the extensive winter wheat yield fluctuations of this region. Therefore, in our study, SPEI was selected as the drought indicator to evaluate the different response patterns of winter wheat yield to drought.
Due to the significant impact of drought on crop yield, a series of studies on the crop yield responses to drought have been carried out at different temporal and spatial scales. At the global scale, the nonlinear responses between crop yield and drought severity have been found by analyzing the responses of wheat, corn, rice, and soybeans to drought [18]. For maize, soybeans, and wheat, drought has significantly reduced their yields globally [19]. By analyzing the sensitivity of 10 main crops to drought in the United States, Lu J. et al. (2020) [20] pointed out that drought stress that occurred in the reproductive stage will severely reduce crop yield. Based on the field experiment on maize, the increased sensitivity to drought stress of maize yield has been found in the central United States [21]. Using the grid crop models to quantify the contributions of various climatic drivers to past yield variability in maize and winter wheat, Webber H. et al. (2018) found that drought stress persists as the main driver of losses for the both crops [22]. Based on the APSIM crop model simulation, the reduced rainfall was considered as one of the reasons of wheat yield decline in Australia [23]. By the empirical relationships between drought intensity and yield fluctuation which were identified by the Bayesian hierarchical model, Chen H. et al. (2018) pointed out that drought significantly reduced grain yields in 28 of 31 provinces in China and obvious spatial variability in drought sensitivity exists [24]. In the NCP, Yu, H. et al. (2018) studied the drought influences on winter wheat at different growing periods in Eastern China and found that ensuring the adequate irrigation is very important during the flowering and filling stages to guarantee normal winter wheat yield [25]. By investigating the drought effect on winter wheat yields of the region, Liu X. et al. (2018) pointed out that the SPEI could explain the extensive winter wheat yield fluctuations of this region [17].
Climate factors play an import part in driving the crop yield variations. Numerous studies have explored how climate factors drive the crop yield variability. Based on the detailed crop statistics time series, Ray D.K. et al. (2015) pointed out that climate variability led to about one third of the observed crop yield variability, including maize, rice, wheat, and soybean [26]. By crop model simulation and statistical regression methods, Qiao J. et al. (2017) [27] found that precipitation contributed more than other climatic factors in driving the crop yield variations at the semiarid region. To understand the drivers of crop production in sub-Saharan Africa, Hoffman, A.L. et al. (2017) [28] explored the relationship between crop yields and climate data with random forest, and they found that though climate variables played a comparatively small role compared technological advances, crop yields exhibited distinct responses to maximum temperature, accumulated precipitation, and vapor pressure deficit. The spatial differences in the contribution of different climatic drivers to crop yield variability also have been explored. Ceglar A. et al. (2016) [29] found the meteorological factors that dominate the variations in winter wheat and corn yield are different in different regions of France.
The above research conclusions can improve our understanding on the response patterns of crop yield to drought and how climatic factors influence crop yields. However, as for the drivers leading to the different responses of crop yield to drought, current studies are mainly carried out from the perspective of crop physiology to explain the reasons causing the different responses of winter wheat yield to drought. The differences of the winter wheat varieties are the internal cause for the different responses of winter wheat to drought. The physiological responses (such as stomatal regulation ability, osmotic adjustment capability) to drought is different among the different varieties [30]. Climatic factors are considered as the important external factors leading to the different response patterns of crop yield to drought. However, there are only few studies that have focused on which climatic factors contribute more to the different response patterns of crop yield to drought. At the global scale, Leng et al. (2019) [18] found that temperature plays an important role in determining the impact of drought by sensitivity analysis. In the United States, Pena-Gallardo et al. [31] found that the different response patterns of crop yield to drought were mainly controlled by the average climatic conditions, especially by the water availability. As mentioned above, in the NCP, previous studies mainly focused on the relationship between drought index and crop productions. Some studies explore the drivers leading to the different responses only from the perspective of vegetation physiology [30]. However, in the NCP, our understanding on which climatic factors drive the different response patterns of winter wheat yield to drought is very limited. In order to minimize drought-related winter wheat yield losses, we need to further study the role of different climatic factors in driving winter wheat drought response.
The NCP is an important winter wheat production base in China. Drought condition of the NCP has worsened in the past several decades [32]. Hence, studying the response patterns of winter wheat yield to drought, especially exploring the climatic factors causing the different response patterns, are significant for reducing the drought-related winter wheat yield losses. The main purposes of our research are (i) to identify the possible spatial patterns of winter wheat yield responses to different time-scale droughts and (ii) to explore the climatic characteristics that determine these patterns.

2. Materials and Methods

2.1. Study Region

The North China Plain (NCP, 113° E~123° E, 32° N~42° N) (Figure 1) is the main winter wheat production base and supplies 50% of China’s winter wheat production. The rainfall during the winter wheat growing season is about 100–350 mm. According to the rainfall of the winter wheat growing season, we divided the NCP into three parts: (I) the north parts, (II) the central parts, and (III) the south parts. Significant seasonal precipitation anomalies caused by the monsoon climate can be observed in this region. As a result, severe drought often occurs in the area, and drought has become one of the most prominent factors restricting the growth of winter wheat in the region.

2.2. Data

The main dataset applied in this study consists of meteorological data, observed winter wheat production records, and soil data. (1) There are 47 stations (including 31 meteorological stations and 16 agrometeorological stations) in this region that were chosen for data collection. Daily meteorological data from 1980 to 2013 were obtained from the China Meteorological Data Sharing Service System (http://data.cma.cn/), which includes the temperature (°C), precipitation (mm), relative humidity (%), sunshine duration (h), and daily average wind speed (m/s). We calculated the solar radiation (MJ/m2) using the sunshine duration by the Ångström function, and the potential evapotranspiration (PET) (mm) was estimated according to the Penman–Monteith equation [33]. (2) The soil data was collected from the Food and Agriculture Organization of the United Nations (FAO), it includes the number of soil layers and the properties for every layer. (3) The observed winter wheat yield at 16 agrometeorological stations for the crop model calibration was obtained from the National Climate Center of China.

2.3. Methods

The methodological framework of our study is shown in Figure 2. The main flowchart of our study includes (1) collect daily climate data, basic geographical data and agrometeorological data. (2) Calculate the different time-scale SPEI as the drought index. (3) Simulate the winter wheat yield by the calibrated EPIC crop model. (4) Study the response patterns of winter wheat yield to drought using the Pearson correlation analysis and principle component analysis (PCA) method. (5) Explore the climatic drivers by the predictive discriminant analysis (PDA) method.

2.3.1. Standardized Precipitation Evapotranspiration Index (SPEI)

The PET which is estimated by the Penman–Monteith equation and the precipitation were used to calculate the SPEI. SPEI is not only sensitive to the evaporation changes but also can describe the multiscale drought character. Details of the approach on how to compute the drought index have been described by Vicente-Serrano et al. [13]. In our research, we calculated the 1-month to 12-month time-scale SPEI from 1980 to 2013 in the NCP. The classifications of different levels of drought and moist by SPEI are shown in Table 1.

2.3.2. Environmental Policy Integrated Climate (EPIC) Model and its Localization

The EPIC crop model has been widely applied in the research on how drought affects crop yield due to its advantages in simulating crop growth processes; it has been applied in studies in the NCP for different topics such as irrigation management [34], potential winter wheat yield simulation [35], and crop drought risk assessments [36]. To simulate the winter wheat yield, the input data mainly consists of meteorological data, soil data, and field management data. The simulated principles and processes of the crop growth model are shown in Figure 3. Harvest index is used to convert the above-ground biomass to crop yield in the crop model. The details of the crop model can be found in Wang et al. [37].
Crop model localization helps to determine winter wheat variety parameters in the different regions of the NCP. The cultivars of winter wheat may change during the study periods. To remove the effects of changes in wheat varieties, in our research, by crop model calibration, we would obtain crop parameters that can represent the average condition of winter wheat varieties for many years in the NCP. This set of crop parameters remains unchanged throughout the simulation process, so that we can eliminate the impact of variety differences over the years. According to previous study [38,39], the following eight crop parameters (Table S1) were chosen to calibrated: (1) biomass–energy ratio (WA), (2) harvest index (HI), (3) potential heat unit (PHU), (4) fraction of growing season when leaf area declines (DLAI), (5) first point on optimal leaf area development curve (DLAP1), (6) second point on optimal leaf area development curve (DLAP2), (7) maximum potential leaf area index (DMLA) (8) leaf area index decline parameters (RLAD), which are most sensitive parameters for winter wheat yield simulation.
In this study, 16 agrometeorological stations (Figure 1) were selected to calibrate and verify the crop model based on the winter wheat plant region and the available data. A global optimization algorithm named Shuffled Complex Evolution algorithm—University of Arizona (SCE-UA) [40] was used to verify the EPIC model at these stations. Compared with traditional “trial-and-error” method, the SCE-UA is a global optimization method. There are about 180 yield records of all the stations, the observation yield records from 1999 to 2008 were used for model calibration and the observation yield records from 2009 to 2011 were used for model validation. The normalized root mean square error (N_RMSE) and Pearson correlation coefficients were used to measure the precision of the model localization.
RMSE = i = 1 n ( Yield si Yield oi ) 2
Yiel d o _ mean = 1 n i = 1 n Yiel d oi
N _ RMSE = RMSE Yield d o _ mean
where n is the number of samples, Yieldsi is the simulated winter wheat yield, Yield oi is the observed winter wheat yield and Yieldo_mean is the average observed yield. N_RMSE could express the differences between simulated productions and observed productions.
Irrigation management has influences on winter wheat in the NCP. However, it is difficult to get enough irrigation information. In our study, the averaged sowing date (Figure S1a) was used in the EPIC model, because the observed sowing date were only delayed on averaged by 1.5 days per decade in the NCP [41]. In the same method, we calculated the averaged wintering date (Figure S1b), the averaged jointing date (Figure S1c), and the averaged filling date (Figure S1d) as the reference of the irrigation date. The irrigation schedule for NCP was set up based on our field investigations and previous research conclusions [42,43,44]. In the north part of the NCP, we irrigated three times at wintering stage, jointing stage, and filling stage; in the central part of the region, we irrigated twice at jointing stage and filling stage; in the south part of the region, we only irrigated one time in the jointing stage (Table S2).

2.3.3. Statistical Methods

Three kinds of statistical methods were used in the study. The Pearson correlation analysis was used to evaluate the drought influence on winter wheat yield. The principal component analysis (PCA) method was used to identify the responses of winter wheat yield to drought. Additionally, the predictive discriminant analysis (PDA) was chosen to explore the driving climatic factors of winter wheat yield responses to different time-scale droughts.
(1)
Pearson correlation analysis
Drought influences on the winter wheat yield were assessed by Pearson correlation analysis. The correlation coefficients were computed between the winter wheat yield and the 1-month to 12-month time-scale SPEI series at each station (including meteorological stations and agrometeorological stations) in the winter wheat growing season (from October of the year to June of the following year). Thus, for each station, we obtained 108 correlation coefficients (9 months × 12 time-scales) (Figure 4).
We get the maximum correlation coefficients at each station as follows:
R i , j = cor ( Yield , SPE I i , j ) i w int erw heatgrowingseason , 1 j 12
w int erw heatgrowingseason = { 10 , 11 , 12 , 1 , 2 , 3 , 4 , 5 , 6 }
R max = max ( R i , j )
where cor is the Pearson correlation analysis, i is the month which belongs to the winter wheat growing season, j represents the time-scale of drought, ranging from 1 to 12 months; Yield is the winter wheat yield simulation series at one station from 1981 to 2013; SPEIi,j is the ith month drought index with time-scale of j; Ri,j is the Pearson correlation of Yield and SPEIi,j; Rmax is the maximum correlation coefficients at each station. As the simulation length is 33 years, the values of 0.35 are responses to 5% significant levels of those correlation coefficients.
(2)
The principal component analysis (PCA) in spatial model
The principal component analysis (PCA) in spatial model [45] was chosen to identify the response patterns of winter wheat yield to different time-scale droughts. In our study, 47 stations were chosen to simulate the winter wheat yield, and 108 correlation coefficients were obtained at each station. Therefore, we can get a correlation matrix (47 rows × 108 columns). The PCAs decompose the matrix into two parts: loading matrix and score matrix (Figure 5).
The principal component equals the one special column of the loading matrix multiplied by the corresponding row of the score matrix (Figure 5). Each principal component has different variance. The larger the variance is, the more representative the principal component is. By comparing the variance of each principal component, several representative principal components were selected, which own stronger ability to characterize the different responses of winter wheat yield to drought. The classifications of the different response patterns were according to the maximum principal component loading rules [46]. Among different stations, the closer the principal component loadings are, the more similar the response patterns of winter wheat yield to drought are. In other words, the principal component loading spatial map could help us to identify stations with similar response patterns of winter wheat yield to drought.
Materials and Methods should be described with sufficient details to allow others to replicate and build on published results. Please note that publication of your manuscript implicates that you must make all materials, data, computer code, and protocols associated with the publication available to readers. Please disclose at the submission stage any restrictions on the availability of materials or information. New methods and protocols should be described in detail while well-established methods can be briefly described and appropriately cited.
(3)
Predictive discriminant analysis
Predictive discriminant analysis (PDA) is commonly used to explain the value of a dependent categorical variable based on its relationship to one or more predictors. Given a set of independent variables, PDA attempts to identify linear combinations of those variables that best separate the groups of cases of the dependent variable. According to prior studies [31,46], the PDA method was chosen to study the driving factors of winter wheat yield to drought in the NCP. The basic problem solved by the predictive discriminant analysis (PDA) is to evaluate the relative importance of multiple independent variables to the dependent variables [47]. In this study, the PDA method was used to explore the driving climatic factors of winter wheat yield responses to different time-scale droughts. A set of climatic variables including precipitation, maximum temperature, minimum temperature, mean temperature, potential evapotranspiration, water balance (precipitation minus potential evapotranspiration) were considered as the predictive variables (independent variables), and the winter wheat yield response patterns to drought were considered as the grouping variables (dependent variables). We can get several predictive discriminate functions (Equation (7)) with different discriminate ability [48] after the PDA method. The first function which is able to separate the different response patterns as much as possible was defined automatically in the PDA, then is the second function and further functions until the number of functions reaches the maximum number according to the number of response patterns. These functions were a multiple linear combination of climatic variables which can best separate the response patterns of winter wheat yield to drought.
Func ( clim   ate   factors ) = w 1 × c f 1 + w 2 × c f 2 + + w i × c f i
where cfi represents the different climate factors, Wi is the coefficients of one specific climatic factor. If the absolute value of W1 is greater the absolute value of W2, it means cf1 has greater effect than cf2 in driving the response patterns of winter wheat yield to drought in the NCP.

3. Results

3.1. EPIC Crop Model Verification

EPIC crop model verification is the basis for the application of the model in the NCP. The cultivar parameters of the EPIC crop model were calibrated by the SCE-UA method at 16 agrometeorological stations. At each station, Firstly, the crop yield records from 1999 to 2008 were used for EPIC model calibration, the correlation coefficient between the observed and simulated winter wheat yield was 0.86 (p < 0.01), and the N_RMSE was 9.8% (Figure 6a). Secondly, the crop yield records from 2009 to 2011 were used for EPIC model validation, the correlation coefficient was 0.65 (p < 0.01), and the N_RMSE was 10.9% (Figure 6b). The N_RMSE is within the acceptable range [49], which indicated that the EPIC model after calibration was suitable for simulating winter wheat yield in the NCP [50].

3.2. Diverse and General Spatial–Temporal Patterns of Winter Wheat Yield Responses to Drought

3.2.1. Diverse Responses of Winter Wheat Yield to Drought

For each month of all 47 stations, we obtained the time-scales of SPEI (Figure 7a) and the correlation coefficients (Figure 7b), when the winter wheat yield and SPEI correlation coefficients are the highest. The mean value of the time-scales of SPEI was usually not less than 6 months (Figure 7a). Additionally, drought occurring in April and May usually have higher correlation coefficients with winter wheat yield compared with those of droughts in other months (Figure 7b), which indicates that drought occurring in April and May would have a more serious influence on winter wheat. In summary, the winter wheat yield is usually more sensitive to the medium (6–9 months) time-scale drought and long (9–12 months) time-scale drought during the critical period of winter wheat growth, such as April and May.
When the correlation coefficient reaches the highest between winter wheat yield and SPEI, we obtained the best time-scales of SPEI and the months at each station. The best time-scales of SPEI varied obviously from the north part to the south part of the NCP (Figure 8a). In the north part of the region, the long time-scale drought (9–12 months) has a more significant impact on winter wheat yield. However, in the central part of the region, winter wheat yield is more susceptible to the medium time-scale drought (6–9 months). For the south part of the region, winter wheat yield tends to be affected by short time-scale drought (1–3 months), though the influence was not significant (p > 0.05). This may be because the rainfall in the southern part of the NCP is relatively abundant and the potential evapotranspiration is smaller compared with those in the northern and central parts of the NCP. Overall, medium and long time-scale drought occurring in the critical period of winter wheat growth (heading stages, flowering stages, fill stages during April and May) has more obvious influences on the winter wheat yield (Figure 8b,c) in the NCP.

3.2.2. General Spatial Patterns of the Winter Wheat Yield Responses to Drought

The PCAs method was applied for generalizing the general spatial response patterns of winter wheat yield to different time-scale droughts. Different principal component represents the different response patterns. The larger the variance of one principal component is, the more representative of the principal component is. As Figure 9 shows, the variance of the first two principal components (the principal component 1 and the principal component 2) occupy approximately 90% of the total variance. This means a clear response pattern of the relationship between winter wheat productions and the different time-scale droughts can be summarized by the first two principal components. Therefore, the first two principal components were selected for further analysis of the characteristics and spatial distributions of these two main response patterns.
PC1 presents strong winter wheat yield-SPEI correlations at long time-scales drought (9–12 months), especially from March to May (Figure 10b). The stations that show significant responses to long time-scale drought (PC1) are mainly situated at the north region of the NCP. Conversely, stations located in the south region show a negative relationship between winter wheat yield and SPEI (Figure 10a). Additionally, higher correlations at medium time-scale (6–9 months) of SPEI during May (Figure 10d) were found at the winter wheat yield responses to drought summarized by PC2. PC2 did not exhibit significant correlations between different time-scale of SPEI and winter wheat yield in the NCP (Figure 10c). Overall, the winter wheat planted in more arid stations showed a more significant response to the SPEI. Conversely, winter wheat planted in a relatively humid environment usually show a lower sensitivity to the variability of SPEI (Figure S2). The findings were consistent with the conclusion that crop yield responses to drought are nonlinear [11,51].
To study the climatic driving factors, we used the maximum loading rule to classify the different types of winter wheat yield responses to the different time-scale droughts (Table 2). Winter wheat with responses to long time-scale drought represented 72% of the analyzed stations, but 58% showed a positive response (positive loadings) of winter wheat yield to drought, while the remaining 14% showed a negative response (negative loadings) to drought. This means we were able to obtain two different groups of winter wheat according to PC1, which, respectively, described positive (+) or negative (−) responses of winter wheat yield to long time-scale drought. Adopting the same process to the second principal component, we found that the stations with responses to medium time-scale drought (PC2) represented 28% (19% positive loadings and 9% negative loadings).

3.3. Factors Explaining the Different Responses of Winter Wheat Yield to Drought

To consider the possible influence of climate factors on the responses of winter wheat yield to different time-scale droughts, we analyzed the magnitudes of the mean air temperature, maximum air temperature, minimum air temperature, precipitation, PET, water balance (precipitation minus PET), relative humidity, wind speed, and sun radiation during the winter wheat growth season.
The stations presenting a negative response to long time-scale drought (PC1−) are generally situated in the south region of the NCP, in which the climate conditions are more humid and warmer than the winter wheat stations that showed a positive response to long time-scale drought (PC1+) (Figure 11). Thus, the negative relationship between winter wheat yield and drought are probably caused by higher temperature and precipitation conditions. The main climatic factors that influence the positive responses to the medium time-scale drought (PC2+) and the negative responses to the medium time-scale drought (PC2−) are complex and include wind speed, relative humidity, and sun radiation in different seasons (Figures S3–S6). Overall, winter wheat stations represented by PC2+ are usually with higher relative humidity, wind speed, and sun radiation than those represented by PC2− (Figure 11). When the winter wheat showed positive responses to drought, the winter wheat plant region with a positive response to long time-scale drought (PC1+) were usually much drier than the region with a positive response to medium time-scale drought (PC2+). However, in terms of winter wheat corresponding to negative responses to drought, the region where winter wheat showed higher correlations with long time-scale drought are usually with a higher water balance and higher temperature compared with the areas where winter wheat was more sensitive to medium time-scale drought (Figure 11).
The PDA function quantified the relative importance of the various climatic factors for driving the winter wheat yield response patterns to drought (Table 3). The first PDA function (PDA1) accounts for 62.2% of the total variance and allows the four patterns (PC1+, PC2+, PC1−, PC2−) of winter wheat yield responses to the different time-scale droughts to be differentiated. According to the coefficient of each climatic factor (Table 3), we can find that the PDA1 mainly represents the latitude, winter climatic conditions (including mean temperature, precipitation, and water balance), and mean temperature and max temperature in autumn at each station. In other words, the PDA1 showed that the winter climatic conditions and temperature have an important impact on winter wheat yield response to drought. The second PDA function (PDA2) could explain 30% of the variance and shows significant correlations with the water balance, especially the autumn water balance. The PDA2 revealed that the water balance also plays an important part in influencing the response characteristics. The third PDA function (PDA3, 7.8% of variance) mainly represents the influence of PET on the winter wheat yield responses to drought. In summary, the above three PDA functions highlight the importance of the different climatic factors to discriminate the different patterns of winter wheat yield responses to drought. Additionally, the climate conditions in winter play a more important part in discriminating the different patterns of winter wheat yield responses to drought. Additionally, compared with the water balance, the temperature plays a more important role in driving the response pattern characteristics.
PDA1 discriminates between winter wheat with a positive response to long time-scale drought (PC1+) and a negative response to long time-scale drought (PC1−) well (Figure 12). This means that whether the winter wheat yield show positive responses to long time-scale drought (PC1+) or negative responses to long time-scale drought (PC1−) is mainly controlled by latitude and winter climate conditions (including mean temperature, precipitation, and water balance). This conclusion also applies to the positive response to medium time-scale drought (PC2+) and the negative response to medium time-scale drought (PC2−). It also can be found that the centroids of PC1+ (winter wheat yield with positive response to long time-scale drought) and PC2− (winter wheat yield with negative response to medium time-scale drought) are close to each other, which means that the different responses to drought characterized by PC1+ and PC2− are not mainly controlled by the above climate conditions. However, PDA2, mainly representing water balance, was able to distinguish between wheat with PC1+ and PC2−, which means the differences between the positive response to long time-scale drought (PC1+) and the negative response to medium time-scale drought (PC2−) are more susceptible to water balance. Though PDA3 only accounts for 7.8% of the total variance, it can be considered to represent the impacts of PET, especially the summer PET and spring PET (Table 3). PDA3 can clearly differentiate between winter wheat that has a positive response to medium time-scale drought (PC2+, humid but with low PET) and winter wheat that has a negative response to long time-scale drought (PC1−, humid with high PET).

4. Discussion

The results of our study clearly showed that the correlation between the winter wheat production and drought is spatially variable and relied on the drought time scales quantified by the SPEI. However, there is still a uniform pattern in the responses of winter wheat yield to drought: the medium time-scale (6–9 months) drought and long time-scale (9–12 months) drought occurring in the key growth periods (April and May) usually have the greatest impact on the yield in the north and central region of the NCP, though winter wheat planted in the south part of the NCP were more sensitive to the short time-scale (1–3 months) drought. As for the other question, which climatic factors lead to this kind of response pattern, we found climate conditions in winter play a more important part in discriminating the different patterns of winter wheat yield responses to drought. Additionally, compared with the water balance, the temperature plays a more important role in driving the response pattern characteristics of winter wheat yield to drought.

4.1. The Mechanisms of the Different Responses of Winter Wheat Yields to Drought

Winter wheat yields are easily influenced by severe droughts that occur in the key growth stages from April to May compared with other stages, just as Xu et al. (2018) has found [52]. In addition, a series of studies which aimed at different crop types and were carried out in different regions all indicated that the crop yields were more easily influenced by drought that occurred at the critical growth stages [11,31,53]. There are several reasons leading to the phenomenon. Firstly, the different growth stages own different resistances to drought [54] due to the different physiological characteristics and the different field management [55] of each growth period. When drought occurred in the earlier stage, winter wheat may have enough time to remedy itself with a positive irrigation management and the influence of drought on winter wheat will be alleviated [52]. Secondly, droughts that occurred at different growth stages have different effects on yield formation. During the flowering stage, drought mainly reduced the success rate of pollination and increased the floret death and which leads to a reduction of spikes number. Drought that occurred at the wheat filling periods affects the transfer of substances from other organs to wheat ears, resulting a thousand-grain weight reduction [56].
Different climatic conditions from the north region to the south region of the NCP caused the different sensitivity of winter wheat yield to the time-scale of drought. In the relative dry region (including the north region and the central region of the NCP), winter is the key season for the soil moisture supply, because of the low atmospheric evaporative demand caused by the low temperature [57,58,59]. The medium and long time-scale SPEI of April and May have considered the water supply condition in winter, so they show higher relationship with the winter wheat yield in the north and central part of the NCP. However, in the south part of the NCP, the precipitation and temperature conditions are more suitable for winter wheat growth than those at the north and central region of the NCP. The more suitable climate conditions lead to a short time-scale SPEI that can directly represent the soil water content, so in the south region of the NCP, short time-scale SPEI has high correlations with winter wheat yield. Meanwhile, the correlations between winter wheat yield and SPEI in the southern region are not significant, which indicates that the main restricted factor of winter wheat production in this region is not drought stress, as some previous studies have found [60]. By counting the number of winter wheat growing season with extreme moist degree at each station based on the SPEI, we also found that the stations where extreme moist occurred are mainly distributed in the southern part of the NCP (Figure S7). It indicates that excessive precipitation in some years may adversely affect wheat production in the south part of the NCP.

4.2. Global Warming Enhances the Role of Temperature in Driving the Response Patterns of Winter Wheat Yield to Drought Compared with Other Climatic Factors

Studies focusing on the driving factors of natural vegetation activities and crop yield responses to the different time-scales of drought have revealed that climate characteristics may be the main controller of the different responses [31,46,61]. Those studies suggested that the water balance (precipitation minus PET) and precipitation are usually the most influential of the different climatic variables. Generally, the different response patterns to the different time-scale SPEI could be explained by the different resistance to water deficiency between the different vegetation [62] and the varied vegetation strategies for coping with drought [63]. As in previous studies, the importance of the spatial differences in water balance, precipitation, and PET to the winter wheat yield responses to drought was also found in this study. However, compared with above three climatic factors, temperature plays a more import role in influencing the different spatial response patterns of winter wheat yield to drought in the NCP.
Some prior studies have also emphasized the importance of temperature in influencing crop yields. Lobell et al. (2011) pointed out that the majority climate influences on crop yields were driven by temperature trends rather than the rainfall trends in most regions of the world [64]. Schauberger B. et al. (2017) [65] research found that water stress caused by high temperature is the main driving force of the yield decline. By sensitivity analysis, Leng G. et al. (2019) [18] also revealed that temperature plays an important role in determining the impacts of drought under global warming. A series of studies in the NCP have found that the temperature of the winter wheat growth season has risen obviously, while precipitation has not changed significantly [66,67]. The length of the winter wheat growth period is susceptible to temperature [68]. The reproductive growth duration of winter wheat has extended because of the shift in the flowering date caused by significant warming trends [67,69]. Thus, the increasing temperature has increased the probability of winter wheat suffering from drought during the critical growth stage in the NCP. On the other hand, atmospheric water demand increases with temperature as an immediate effect [65]. The above studies may be used to explain why temperature has become increasingly important in affecting the response of winter wheat yield to drought. Although a large number of studies on how temperature affects crop production and food security have been carried out and have reached different conclusions [65,70], how temperature rise affects crop yield responses to drought requires further research and analysis against the backdrop of global warming.

4.3. The Contributions and Limitations of Our Study

Effective drought monitoring and management is becoming increasingly essential in the context of climate change. Selecting the appropriate time scale for a drought indicator (SPEI, etc.) is essential for a reasonable assessment of the influence of drought on crop yield. Our results contribute to our understanding on the response patterns of winter wheat yield to drought in the NCP. Such understanding can help us to choose a suitable time scale SPEI to evaluate the drought influence on winter wheat yield. In addition, our findings emphasize the importance of temperature in driving the response patterns of winter wheat yield to drought, this could provide references for choosing a reasonable sowing time and suitable wheat varieties.
It is important to note there are several limitations in our study. Firstly, though SPEI has been widely used in exploring the drought impacts on crop production in the NCP according to previous studies [16,17]. However, there are still limitations for the SPEI to monitor the drought condition in the irrigated area. Some drought indices which consider the irrigation influences, have been proposed [69]. The further study should collect more irrigation data and use these indices to evaluate the drought response characteristics. Secondly, the EPIC crop model was used in our study. However, there are only 16 agrometeorological stations with enough data to verify the crop model, for other stations with incomplete data, we use the crop model parameters of the site closest to it, which may reduce the simulation accuracy of winter wheat yield.

5. Conclusions

In this research, we investigated the responses of winter wheat yield to different time-scale droughts and then explored the climatic factors which may determine the different response patterns. Our results suggested that (i) the responses of winter wheat yield to different time-scale droughts have obvious spatial differences from the north part to the south part in the NCP. In the north part of the NCP, the long time-scale (9–12 months) droughts in the key growth periods (April and May) usually have the greatest impact on the winter wheat yield. In the central part of the NCP, the medium time-scale (6–9 months) droughts are more likely to affect wheat yield. However, winter wheat planted in the south part of the NCP were more sensitive to the short time-scale (1–3 months) droughts. (ii) The different patterns of winter wheat yield responses to different time-scale droughts are mainly controlled by the climate conditions (including mean temperature, precipitation, and water balance) in winter in the NCP. Global warming enhances the role of temperature in driving the response patterns of winter wheat yield to drought. Compared with the water balance, the temperature plays a more important role in driving the response pattern characteristics. (iii) The water availability and PET also affect the responses of winter wheat yield to drought. Winter wheat planted in the humid environments are easily influenced by shorter time-scale drought, while winter wheat grown under dry conditions are more sensitive to longer time-scale drought. In the context of global warming, our study can provide a reference on how to reduce drought influence on winter wheat yield and the results are also effective at sustaining long-term winter wheat productivity in the NCP.

Supplementary Materials

The following are available on line at https://www.mdpi.com/2073-4441/12/11/3094/s1, Figure S1: Long-term mean sowing dates (a), wintering dates (b), jointing dates (c) and filling dates (d), Figure S2: Distribution of average precipitation during the winter wheat growing season from 1981 to 2013, Figure S3: Mean spring values of the main climatic factors characterized by the first two principal components of the winter wheat responses to different time-scale SPEI obtained in a PCA method. Error bars represent 1 standard error of the mean value, Figure S4: Same as Figure S3, but for summer, Figure S5: Same as Figure S3, but for autumn, Figure S6: Same as Figure S3, but for winter, Figure S7: The count of winter wheat growing season with different moist degree at each station from 1980 to 2013, Table S1: Calibrated parameters for winter wheat, its default value and calibration range, Table S2: The irrigation schedule.

Author Contributions

Methodology, J.Y.; validation, J.Y., H.Z.; formal analysis, J.Y, L.L.; data curation, J.Y.; A.G.; X.H.; W.Z.; writing—original draft preparation, J.Y.; visualization, J.Y.; L.L.; supervision, J.W.; L.L.; funding acquisition, J.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (2017YFC1502402); the National Natural Science Foundation of China (No. 41671424).

Acknowledgments

We thank Qianfeng Wang for his advice on EPIC crop model localizations, we thank Cunjie Zhang and Juqi Duan for the data collection.

Conflicts of Interest

The authors declare that we do not have any commercial or associative interest that represents a conflict of interest in connection with the work submitted.

References

  1. Zampieri:, M.; Ceglar, A.; Dentener, F.; Toreti, A. Wheat yield loss attributable to heat waves, drought and water excess at the global, national and subnational scales. Environ. Res. Lett. 2017, 12, 064008. [Google Scholar] [CrossRef]
  2. Daryanto, S.; Wang, L.; Jacinthe, P.-A. Global Synthesis of Drought Effects on Maize and Wheat Production. PLoS ONE 2016, 11, e0156362. [Google Scholar] [CrossRef] [PubMed]
  3. Hao, L.; Zhang, X.; Liu, S. Risk assessment to China’s agricultural drought disaster in county unit. Nat. Hazards 2012, 61, 785–801. [Google Scholar] [CrossRef]
  4. He, B.; Lü, A.; Wu, J.; Zhao, L.; Liu, M. Drought hazard assessment and spatial characteristics analysis in China. J. Geogr. Sci. 2011, 21, 235–249. [Google Scholar] [CrossRef]
  5. Madadgar, S.; AghaKouchak, A.; Farahmand, A.; Davis, S.J. Probabilistic estimates of drought impacts on agricultural production. Geophys. Res. Lett. 2017, 44, 7799–7807. [Google Scholar] [CrossRef]
  6. Wilhite, D.A.; Glantz, M.H. Understanding: The Drought Phenomenon: The Role of Definitions. Water Int. 1985, 10, 111–120. [Google Scholar] [CrossRef] [Green Version]
  7. Paulo, A.A.; Pereira, L.S. Drought Concepts and Characterization. Water Int. 2006, 31, 37–49. [Google Scholar] [CrossRef]
  8. Lloyd-Hughes, B. The impracticality of a universal drought definition. Theor. Appl. Clim. 2014, 117, 607–611. [Google Scholar] [CrossRef] [Green Version]
  9. Zeng, Z.; Li, Y.; Wu, W.; Zhou, Y.; Wang, X.; Huang, H.; Li, Z. Spatio-Temporal Variation of Drought within the Vegetation Growing Season in North Hemisphere (1982–2015). Water 2020, 12, 2146. [Google Scholar] [CrossRef]
  10. Beguería, S.; Vicente-Serrano, S.M.; Reig, F.; Latorre, B. Standardized precipitation evapotranspiration index (SPEI) revisited: Parameter fitting, evapotranspiration models, tools, datasets and drought monitoring. Int. J. Climatol. 2014, 34, 3001–3023. [Google Scholar] [CrossRef] [Green Version]
  11. Zipper, S.C.; Qiu, J.; Kucharik, C.J. Drought effects on US maize and soybean production: Spatiotemporal patterns and historical changes. Environ. Res. Lett. 2016, 11, 094021. [Google Scholar] [CrossRef]
  12. Guttman, N.B. Comparing the Palmer Drought Index and the Standardized Precipitation Index. J. Am. Water Resour. Assoc. 1998, 34, 113–121. [Google Scholar] [CrossRef]
  13. Vicente-Serrano, S.M.; Beguería, S.; López-Moreno, J.I. A Multiscalar Drought Index Sensitive to Global Warming: The Standardized Precipitation Evapotranspiration Index. J. Clim. 2010, 23, 1696–1718. [Google Scholar] [CrossRef] [Green Version]
  14. Alley, W.M. The Palmer Drought Severity Index: Limitations and Assumptions. J. Climate Appl. Meteor. 1984, 23, 1100–1109. [Google Scholar] [CrossRef] [Green Version]
  15. Wang, H.; Vicente-Serrano, S.M.; Tao, F.; Zhang, X.; Wang, P.; Zhang, C.; Chen, Y.; Zhu, D.; Kenawy, A.E. Monitoring winter wheat drought threat in Northern China using multiple climate-based drought indices and soil moisture during 2000–2013. Agric. For. Meteorol. 2016, 228–229, 1–12. [Google Scholar] [CrossRef]
  16. Liu, X.; Zhu, X.; Pan, Y.; Bai, J.; Li, S. Performance of different drought indices for agriculture drought in the North China Plain. J. Arid. Land 2018, 10, 507–516. [Google Scholar] [CrossRef] [Green Version]
  17. Liu, X.; Pan, Y.; Zhu, X.; Yang, T.; Bai, J.; Sun, Z. Drought evolution and its impact on the crop yield in the North China Plain. J. Hydrol. 2018, 564, 984–996. [Google Scholar] [CrossRef]
  18. Leng, G.; Hall, J. Crop yield sensitivity of global major agricultural countries to droughts and the projected changes in the future. Sci. Total Environ. 2019, 654, 811–821. [Google Scholar] [CrossRef]
  19. Matiu, M.; Ankerst, D.P.; Menzel, A. Interactions between temperature and drought in global and regional crop yield variability during 1961-2014. PLoS ONE 2017, 12, e0178339. [Google Scholar] [CrossRef]
  20. Lu, J.; Carbone, G.J.; Huang, X.; Lackstrom, K.; Gao, P. Mapping the sensitivity of agriculture to drought and estimating the effect of irrigation in the United States, 1950–2016. Agric. For. Meteorol. 2020, 292–293, 108124. [Google Scholar] [CrossRef]
  21. Lobell, D.B.; Roberts, M.J.; Schlenker, W.; Braun, N.; Little, B.B.; Rejesus, R.M.; Hammer, G.L. Greater sensitivity to drought accompanies maize yield increase in the U.S. Midwest. Science 2014, 344, 516–519. [Google Scholar] [CrossRef]
  22. Webber, H.; Ewert, F.; Olesen, J.E.; Müller, C.; Fronzek, S.; Ruane, A.C.; Bourgault, M.; Martre, P.; Ababaei, B.; Bindi, M.; et al. Diverging importance of drought stress for maize and winter wheat in Europe. Nat. Commun. 2018, 9, 4249. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Hochman, Z.; Gobbett, D.L.; Horan, H. Climate trends account for stalled wheat yields in Australia since 1990. Glob. Chang. Biol. 2017, 23, 2071–2081. [Google Scholar] [CrossRef] [PubMed]
  24. Chen, H.; Liang, Z.; Liu, Y.; Jiang, Q.; Xie, S. Effects of drought and flood on crop production in China across 1949–2015: Spatial heterogeneity analysis with Bayesian hierarchical modeling. Nat. Hazards 2018, 92, 525–541. [Google Scholar] [CrossRef]
  25. Yu, H.; Zhang, Q.; Sun, P.; Song, C. Impact of Droughts on Winter Wheat Yield in Different Growth Stages during 2001–2016 in Eastern China. Int. J. Disaster Risk Sci. 2018, 9, 376–391. [Google Scholar] [CrossRef] [Green Version]
  26. Ray, D.K.; Gerber, J.S.; MacDonald, G.K.; West, P.C. Climate variation explains a third of global crop yield variability. Nat. Commun. 2015, 6, 5989. [Google Scholar] [CrossRef] [Green Version]
  27. Qiao, J.; Yu, D.; Liu, Y. Quantifying the impacts of climatic trend and fluctuation on crop yields in northern China. Environ. Monit. Assess. 2017, 189, 532. [Google Scholar] [CrossRef]
  28. Hoffman, A.L.; Kemanian, A.R.; Forest, C.E. Analysis of climate signals in the crop yield record of sub-Saharan Africa. Glob. Chang. Biol. 2018, 24, 143–157. [Google Scholar] [CrossRef]
  29. Ceglar, A.; Toreti, A.; Lecerf, R.; van der Velde, M.; Dentener, F. Impact of meteorological drivers on regional inter-annual crop yield variability in France. Agric. For. Meteorol. 2016, 216, 58–67. [Google Scholar] [CrossRef]
  30. Feng, F.; Huang, G.; Wang, Y.; Zhang, X.; Wang, T.; Zhong, X. Cultivar Difference in Yield and Physiological Response of Winter Wheat after Exposed to Drought Stress at Jointing Stage. Chin. J. Agrometerol. 2020, 41, 43–50. [Google Scholar]
  31. Peña-Gallardo, M.; Vicente-Serrano, S.M.; Quiring, S.; Svoboda, M.; Hannaford, J.; Tomas-Burguera, M.; Martín-Hernández, N.; Domínguez-Castro, F.; Kenawy, A.E. Response of crop yield to different time-scales of drought in the United States: Spatio-temporal patterns and climatic and environmental drivers. Agric. For. Meteorol. 2019, 264, 40–55. [Google Scholar] [CrossRef] [Green Version]
  32. Dai, A. Drought under global warming: A review. Wiley Interdiscip. Rev. Clim. Chang. 2011, 2, 45–65. [Google Scholar] [CrossRef] [Green Version]
  33. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration. In Guidelines for Computing Crop Water Requirements; FAO: Roma, Italy, 1998. [Google Scholar]
  34. Liu, M.; Wu, J.; Lv, A.; Zhao, L. Water stress of winter wheat and irrigation strategy in typical region of Huang-Huai-Hai Plain. Trans. Chin. Soc. Agric. Eng. 2010, 26, 40–44. [Google Scholar]
  35. Lu, C.; Fan, L. Winter wheat yield potentials and yield gaps in the North China Plain. Field Crops Res. 2013, 143, 98–105. [Google Scholar] [CrossRef] [Green Version]
  36. Yin, Y.; Zhang, X.; Lin, D.; Yu, H.; Wang, J.A.; Shi, P. GEPIC-V-R model: A GIS-based tool for regional crop drought risk assessment. Agric. Water Manag. 2014, 144, 107–119. [Google Scholar] [CrossRef]
  37. Wang, Q.; Wu, J.; Li, X.; Zhou, H.; Yang, J.; Geng, G.; An, X.; Liu, L.; Tang, Z. A comprehensively quantitative method of evaluating the impact of drought on crop yield using daily multi-scale SPEI and crop growth process model. Int. J. Biometeorol. 2017, 61, 685–699. [Google Scholar] [CrossRef]
  38. Yu, F. Yield Estimation of Winter Wheat with EPIC Model in Huang-Huai-Hai Plain-A Case Study in Luancheng, Hebei. Master’s Thesis, Chinese Academy of Agricultural Science, Beijing, China, 2007. [Google Scholar]
  39. Yue, Y.; Li, J.; Ye, X.; Wang, Z.; Zhu, A.-X.; Wang, J.-A. An EPIC model-based vulnerability assessment of wheat subject to drought. Nat. Hazards 2015, 78, 1629–1652. [Google Scholar] [CrossRef]
  40. Duan, Q.; Sorooshian, S.; Gupta, V. Effective and efficient global optimization for conceptual rainfall-runoff models. Water Resour. Res. 1992, 28, 1015–1031. [Google Scholar] [CrossRef]
  41. Xiao, D.; Tao, F.; Liu, Y.; Shi, W.; Wang, M.; Liu, F.; Zhang, S.; Zhu, Z. Observed changes in winter wheat phenology in the North China Plain for 1981-2009. Int. J. Biometeorol. 2013, 57, 275–285. [Google Scholar] [CrossRef]
  42. Qin, X.; Zhou, G.; Jv, H.; Li, X.; Liu, Q. Simulation of Different Irrigation Strategy on Wheat Yield in Huang-Huai-Hai Plain under the RCP8.5 Scenario. Chin. J. Agrometeorol. 2018, 39, 220–232. [Google Scholar]
  43. Peng, Z.; Zhang, B.; Liu, Y.; Cai, J.; Wang, L.; Du, L. Simulation for Regional Water Consumption and Optimization of Irrigation Schedule for Winter Wheat in North China. Trans. Chin. Soc. Agric. Mach. 2017, 48, 238–246. [Google Scholar]
  44. Fang, W.; Liu, R.; Zhu, Z.; Ma, Z.; Li, S.; Xu, P. Influencing factors and characteristics of winter wheat irrigation requirement in Huanghuai Plain. Chin. J. Ecol. 2009, 28, 2177–2182. [Google Scholar]
  45. Demšar, U.; Harris, P.; Brunsdon, C.; Fotheringham, A.S.; McLoone, S. Principal Component Analysis on Spatial Data: An Overview. Ann. Assoc. Am. Geogr. 2013, 103, 106–128. [Google Scholar] [CrossRef]
  46. Vicente-Serrano, S.M.; Camarero, J.J.; Azorin-Molina, C. Diverse responses of forest growth to drought time-scales in the Northern Hemisphere. Glob. Ecol. Biogeogr. 2014, 23, 1019–1030. [Google Scholar] [CrossRef] [Green Version]
  47. Huberty, C.J. Applied Discriminant Analysis; Wiley: New York, NY, USA, 1994. [Google Scholar]
  48. Hair, J.F.; Black, W.C.; Babin, B.J.; Anderson, R.E. Multivariate Data Analysis, 7th ed.; Pearson New International Edition; Pearson: Harlow, UK, 2014. [Google Scholar]
  49. Jamieson, P.D.; Porter, J.R.; Wilson, D.R. A test of the computer simulation model ARCWHEAT1 on wheat crops grown in New Zealand. Field Crops Res. 1991, 27, 337–350. [Google Scholar] [CrossRef]
  50. Wen, X.; Chen, F. Simulation of climatic change impacts on yield potential of typical wheat varieties based on DSSAT model. Trans. Chin. Soc. Agric. Eng. 2011, 27, 74–79. [Google Scholar] [CrossRef]
  51. Meyer, S.J.; Hubbard, K.G.; Wilhite, D.A. The relationship of climatic indices and variables to corn (maize) yields: A principal components analysis. Agric. For. Meteorol. 1991, 55, 59–84. [Google Scholar] [CrossRef]
  52. Xu, X.; Gao, P.; Zhu, X.; Guo, W.; Ding, J.; Li, C. Estimating the responses of winter wheat yields to moisture variations in the past 35 years in Jiangsu Province of China. PLoS ONE 2018, 13, e0191217. [Google Scholar] [CrossRef]
  53. Potopová, V.; Štěpánek, P.; Možný, M.; Türkott, L.; Soukup, J. Performance of the standardised precipitation evapotranspiration index at various lags for agricultural drought risk assessment in the Czech Republic. Agric. For. Meteorol. 2015, 202, 26–38. [Google Scholar] [CrossRef]
  54. NeSmith, D.S.; Ritchie, J.T. Maize (Zea mays L.) response to a severe soil water-deficit during grain-filling. Field Crops Res. 1992, 29, 23–35. [Google Scholar] [CrossRef]
  55. Wu, H.; Hubbard, K.G.; Wilhite, D.A. An agricultural drought risk-assessment model for corn and soybeans. Int. J. Climatol. 2004, 24, 723–741. [Google Scholar] [CrossRef] [Green Version]
  56. Zhao, S.; Guan, X.; Wu, J. Effects of Water Deficits on Yield and WUE in Winter Wheat. Irrig. Drain. 2001, 20, 56–69. [Google Scholar]
  57. Austin, R.B.; Cantero-Martínez, C.; Arrúe, J.L.; Playán, E.; Cano-Marcellán, P. Yield–rainfall relationships in cereal cropping systems in the Ebro river valley of Spain. Eur. J. Agron. 1998, 8, 239–248. [Google Scholar] [CrossRef] [Green Version]
  58. Liu, C.; Zhang, X.; Zhang, Y. Determination of daily evaporation and evapotranspiration of winter wheat and maize by large-scale weighing lysimeter and micro-lysimeter. Agric. For. Meteorol. 2002, 111, 109–120. [Google Scholar] [CrossRef]
  59. Burba, G.G.; Verma, S.B. Seasonal and interannual variability in evapotranspiration of native tallgrass prairie and cultivated wheat ecosystems. Agric. For. Meteorol. 2005, 135, 190–201. [Google Scholar] [CrossRef]
  60. Zhang, H.; Ma, X.; Peng, N.; Xv, Y.; Zhang, H. Study on assessment of waterlogging disaster losses of winter wheat in Huaihe river basin. J. Environ. Manag. 2015, 31, 123–129. [Google Scholar]
  61. Pasho, E.; Camarero, J.J.; Luis, M.D.; Vicente-Serrano, S.M. Impacts of drought at different time scales on forest growth across a wide climatic gradient in north-eastern Spain. Agric. For. Meteorol. 2011, 151, 1800–1811. [Google Scholar] [CrossRef]
  62. Chaves, M.M.; Maroco, J.P.; Pereira, J.S. Understanding plant responses to drought—From genes to the whole plant. Funct. Plant Biol. 2003, 30, 239. [Google Scholar] [CrossRef]
  63. McDowell, N.; Pockman, W.T.; Allen, C.D.; Breshears, D.D.; Cobb, N.; Kolb, T.; Plaut, J.; Sperry, J.; West, A.; Williams, D.G.; et al. Mechanisms of plant survival and mortality during drought: Why do some plants survive while others succumb to drought? New Phytol. 2008, 178, 719–739. [Google Scholar] [CrossRef]
  64. Lobell, D.B.; Schlenker, W.; Costa-Roberts, J. Climate trends and global crop production since 1980. Science 2011, 333, 616–620. [Google Scholar] [CrossRef] [Green Version]
  65. Schauberger, B.; Archontoulis, S.; Arneth, A.; Balkovic, J.; Ciais, P.; Deryng, D.; Elliott, J.; Folberth, C.; Khabarov, N.; Müller, C.; et al. Consistent negative response of US crops to high temperatures in observations and crop models. Nat. Commun. 2017, 8, 13931. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Tao, F.; Xiao, D.; Zhang, S.; Zhang, Z.; Rötter, R.P. Wheat yield benefited from increases in minimum temperature in the Huang-Huai-Hai Plain of China in the past three decades. Agric. For. Meteorol. 2017, 239, 1–14. [Google Scholar] [CrossRef] [Green Version]
  67. Lv, Z.; Liu, X.; Cao, W.; Zhu, Y. Climate change impacts on regional winter wheat production in main wheat production regions of China. Agric. For. Meteorol. 2013, 171–172, 234–248. [Google Scholar] [CrossRef]
  68. Tao, F.; Zhang, Z.; Xiao, D.; Zhang, S.; Rötter, R.P.; Shi, W.; Liu, Y.; Wang, M.; Liu, F.; Zhang, H. Responses of wheat growth and yield to climate change in different climate zones of China, 1981–2009. Agric. For. Meteorol. 2014, 189–190, 91–104. [Google Scholar] [CrossRef]
  69. Yu, H.; Zhang, Q.; Xu, C.-Y.; Du, J.; Sun, P.; Hu, P. Modified Palmer Drought Severity Index: Model improvement and application. Environ. Int. 2019, 130, 104951. [Google Scholar] [CrossRef]
  70. Zhao, C.; Liu, B.; Piao, S.; Wang, X.; Lobell, D.B.; Huang, Y.; Huang, M.; Yao, Y.; Bassu, S.; Ciais, P.; et al. Temperature increase reduces global yields of major crops in four independent estimates. Proc. Natl. Acad. Sci. USA 2017, 114, 9326–9331. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Study area and the distribution of meteorological stations and agrometeorological stations in the NCP.
Figure 1. Study area and the distribution of meteorological stations and agrometeorological stations in the NCP.
Water 12 03094 g001
Figure 2. Methodological framework of this study.
Figure 2. Methodological framework of this study.
Water 12 03094 g002
Figure 3. The simulated principle and process of the EPIC crop model.
Figure 3. The simulated principle and process of the EPIC crop model.
Water 12 03094 g003
Figure 4. The correlation analysis in one special station. i is the month which belongs to the winter wheat growing season, j represents the time-scale of drought, ranging from 1 to 12 months; Ri,j is the Pearson correlation between Yield and the ith month SPEI with time-scale of j.
Figure 4. The correlation analysis in one special station. i is the month which belongs to the winter wheat growing season, j represents the time-scale of drought, ranging from 1 to 12 months; Ri,j is the Pearson correlation between Yield and the ith month SPEI with time-scale of j.
Water 12 03094 g004
Figure 5. The follow chart of the PCAs method in our study. “Sta” is the abbreviation of “Station”, “PC” is the abbreviation of “principal component”.
Figure 5. The follow chart of the PCAs method in our study. “Sta” is the abbreviation of “Station”, “PC” is the abbreviation of “principal component”.
Water 12 03094 g005
Figure 6. The EPIC crop model verification result at all 16 agrometeorological stations. (a) The calibration result based on the yield record from 1999 to 2008, (b) the validation result based on the yield records from 2009 to 2011.
Figure 6. The EPIC crop model verification result at all 16 agrometeorological stations. (a) The calibration result based on the yield record from 1999 to 2008, (b) the validation result based on the yield records from 2009 to 2011.
Water 12 03094 g006
Figure 7. Box plot: (a) the time-scales of the SPEI which have the highest correlation coefficients with the winter wheat yield from October to June at each station. The red dotted line represents the 6-month time-scale SPEI, the red point in figure a represents the mean value of the time scales of drought. (b) The highest correlation coefficients from October to June at each station between winter wheat yield and SPEI. As the simulation length is 33 years, the values of 0.35 are responses to 5% significant levels of those correlation coefficients (the red dotted line). The red point in (b) represents the mean value of the correlation coefficients.
Figure 7. Box plot: (a) the time-scales of the SPEI which have the highest correlation coefficients with the winter wheat yield from October to June at each station. The red dotted line represents the 6-month time-scale SPEI, the red point in figure a represents the mean value of the time scales of drought. (b) The highest correlation coefficients from October to June at each station between winter wheat yield and SPEI. As the simulation length is 33 years, the values of 0.35 are responses to 5% significant levels of those correlation coefficients (the red dotted line). The red point in (b) represents the mean value of the correlation coefficients.
Water 12 03094 g007
Figure 8. (a) Spatial map of best SPEI timescale and month at each station. The best timescale of the SPEI refers to the time-scale when the correlation coefficient reaches the highest between winter wheat yield and SPEI, the best month of the SPEI refers to the month when the correlation coefficient reaches the highest between winter wheat yield and SPEI. (b) Histograms of months with the best relationship between the SPEI and winter wheat yield across all stations (c) Distributions of the timescales of SPEI with the best relationship between the SPEI and winter wheat yield across all stations.
Figure 8. (a) Spatial map of best SPEI timescale and month at each station. The best timescale of the SPEI refers to the time-scale when the correlation coefficient reaches the highest between winter wheat yield and SPEI, the best month of the SPEI refers to the month when the correlation coefficient reaches the highest between winter wheat yield and SPEI. (b) Histograms of months with the best relationship between the SPEI and winter wheat yield across all stations (c) Distributions of the timescales of SPEI with the best relationship between the SPEI and winter wheat yield across all stations.
Water 12 03094 g008
Figure 9. Scree plots displaying the explanatory capability of each principal component for the response patterns of winter wheat yield to drought by the variance percentage. The greater the variance is, the stronger the explanatory capability of the principal components owns.
Figure 9. Scree plots displaying the explanatory capability of each principal component for the response patterns of winter wheat yield to drought by the variance percentage. The greater the variance is, the stronger the explanatory capability of the principal components owns.
Water 12 03094 g009
Figure 10. (a) The PC-loading map of the first principal components. (b) The PC-scores that that represent the first principal components. (c) The PC-loading map of the second principal components. (d) The PC-scores that that represent the second principal components. We obtained the PC-loading spatial map by the Kriging interpolation method.
Figure 10. (a) The PC-loading map of the first principal components. (b) The PC-scores that that represent the first principal components. (c) The PC-loading map of the second principal components. (d) The PC-scores that that represent the second principal components. We obtained the PC-loading spatial map by the Kriging interpolation method.
Water 12 03094 g010
Figure 11. Mean annual values of the main climatic factors characterized by the first two principal components of the winter wheat responses to different time-scale droughts obtained in a PCA method. Error bars represent 1 standard error of the mean value.
Figure 11. Mean annual values of the main climatic factors characterized by the first two principal components of the winter wheat responses to different time-scale droughts obtained in a PCA method. Error bars represent 1 standard error of the mean value.
Water 12 03094 g011
Figure 12. Centroids of each groups got by the PCAs method, corresponding to the first three PDA functions. The farther the vertical projection distance between different groups on a certain PDA function is, the better the PDA function can distinguish them.
Figure 12. Centroids of each groups got by the PCAs method, corresponding to the first three PDA functions. The farther the vertical projection distance between different groups on a certain PDA function is, the better the PDA function can distinguish them.
Water 12 03094 g012
Table 1. Drought/moisture level based on SPEI.
Table 1. Drought/moisture level based on SPEI.
Drought/Moisture LevelSPEI Values
Mild moistureSPEI ≥ 2.0
Moderate moisture1.5 < SPEI <2.0
Severe moisture1.0 < SPEI ≤ 1,5
Extreme moisture0.5 < SPEI ≤ 1.0
Normal −0.5 < SPEI ≤ 0.5
Mild drought−1.0 < SPEI < −0.5
Moderate drought−1.5 < SPEI ≤ −1.0
Severe drought−2.0 < SPEI ≤ −1.5
Extreme droughtSPEI ≤ −2.0
Table 2. Percentages of the first two principal components (PC1 and PC2) summarizing the winter wheat yield responses to the SPEI.
Table 2. Percentages of the first two principal components (PC1 and PC2) summarizing the winter wheat yield responses to the SPEI.
Group (Effect)PC1PC2
PC1+PC1−PC2+PC2−
Percentage58%14%19%9%
Table 3. Structure matrix of the three PDA functions for winter wheat (the explained variance by each function is indicated in parentheses; the greater the variance is, the stronger the explanatory power of the discriminant function owns.). The correlation values computed for every environmental factor are displayed as well.
Table 3. Structure matrix of the three PDA functions for winter wheat (the explained variance by each function is indicated in parentheses; the greater the variance is, the stronger the explanatory power of the discriminant function owns.). The correlation values computed for every environmental factor are displayed as well.
VariablesPDA Functions (Percentage of Variance)
PDA1 (62.2%)PDA2 (30%)PDA3 (7.8%)
Pre(winter)−0.372 *−0.126−0.008
Latitude0.366 *0.109−0.044
Tmean (winter)−0.334 *−0.064−0.029
Water balance (winter)−0.330 *−0.2200.020
Tmean (autumn)−0.324 *0.016−0.083
Tmax(autumn)−0.324 *−0.069−0.128
Water balance (spring)−0.322 *−0.2340.277
Rhu (winter)−0.317 *−0.0450.061
Pre (spring)−0.315 *−0.1870.046
Tmin (winter)−0.315 *−0.028−0.025
Tmax (winter)−0.311 *−0.116−0.039
Tmin (annual)−0.311 *−0.014−0.151
Rhu (annual)−0.304 *−0.0930.242
Tmean (annual)−0.300 *−0.057−0.232
Pre (autumn)−0.299 *−0.2080.092
Water balance (autumn)−0.294 *−0.2600.262
Rhu (spring)−0.293 *−0.0990.233
Pre (annual)−0.284 *−0.1770.085
Tmin (autumn)−0.276 *0.035−0.062
Rhu (autumn)−0.243 *−0.1240.182
Rs (summer)0.228 *0.197−0.194
Pre (summer)−0.175 *−0.1390.102
Water balance (autumn)−0.254−0.364 *0.145
Rs (autumn)0.2400.267 *−0.009
Rs (annual)0.1430.227 *−0.009
PET (autumn)−0.0320.223 *−0.078
Wind (summer)−0.0210.219 *0.027
Wind (spring)0.0540.216 *0.021
Wind (annual)0.0140.202 *0.055
Wind (autumn)0.0100.201 *0.054
Rs (autumn)0.0440.199 *0.140
Wind (winter)0.0160.167 *0.110
Rs (winter)−0.0670.163 *0.123
PET (winter)−0.0840.103 *−0.034
PET (summer)0.1900.205−0.608 *
PET (spring)0.2470.230−0.488 *
PET (annual)0.1420.283−0.487 *
Tmean (summer)−0.116−0.044−0.433 *
Rhu (summer)−0.211−0.0640.416 *
Tmin (summer)−0.271−0.023−0.381 *
Tmax (summer)−0.024−0.057−0.366 *
Tmean (spring)−0.151−0.076−0.320 *
Tmax (spring)−0.063−0.081−0.312 *
Tmin (spring)−0.260−0.044−0.285 *
Water balance (summer)−0.203−0.1770.263 *
Tmax (annual)−0.188−0.097−0.252 *
Altitude0.0020.0100.199 *
Longitude0.0100.1070.146 *
The environmental factors are marked with asterisks (*) while p-value < 0.05. “Pre” is “precipitation”; “Tmean” is the “mean temperature”; “Tmax” is “maximum temperature”, “Rhu” is “relative humidity”, “PET” is “potential evapotranspiration”, “Water balance” is “precipitation minus potential evapotranspiration”, “Wind” is “mean wind speed”, “Rs” is “sun radiation”.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yang, J.; Wu, J.; Liu, L.; Zhou, H.; Gong, A.; Han, X.; Zhao, W. Responses of Winter Wheat Yield to Drought in the North China Plain: Spatial–Temporal Patterns and Climatic Drivers. Water 2020, 12, 3094. https://doi.org/10.3390/w12113094

AMA Style

Yang J, Wu J, Liu L, Zhou H, Gong A, Han X, Zhao W. Responses of Winter Wheat Yield to Drought in the North China Plain: Spatial–Temporal Patterns and Climatic Drivers. Water. 2020; 12(11):3094. https://doi.org/10.3390/w12113094

Chicago/Turabian Style

Yang, Jianhua, Jianjun Wu, Leizhen Liu, Hongkui Zhou, Adu Gong, Xinyi Han, and Wenhui Zhao. 2020. "Responses of Winter Wheat Yield to Drought in the North China Plain: Spatial–Temporal Patterns and Climatic Drivers" Water 12, no. 11: 3094. https://doi.org/10.3390/w12113094

APA Style

Yang, J., Wu, J., Liu, L., Zhou, H., Gong, A., Han, X., & Zhao, W. (2020). Responses of Winter Wheat Yield to Drought in the North China Plain: Spatial–Temporal Patterns and Climatic Drivers. Water, 12(11), 3094. https://doi.org/10.3390/w12113094

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