Next Article in Journal
Adaptation Strategies of the Hydrosocial Cycles in the Mediterranean Region
Previous Article in Journal
Energy-Recovery Pressure-Reducer in District Heating System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatio-Temporal Analysis of Meteorological Elements in the North China District of China during 1960–2015

1
College of Agronomy and Biotechnology, China Agricultural University, Yuan-Ming-Yuan West Road 2, Haidian, Beijing 100193, China
2
Key Laboratory of Farming system, Ministry of Agriculture and Rural Affairs of the People’s Republic of China, Yuan-Ming-Yuan West Road 2, Haidian, Beijing 100193, China
3
Department of Environmental Sciences, University of California, 900 University Ave, Riverside, CA 92521, USA
*
Author to whom correspondence should be addressed.
Water 2018, 10(6), 789; https://doi.org/10.3390/w10060789
Submission received: 24 May 2018 / Revised: 11 June 2018 / Accepted: 11 June 2018 / Published: 15 June 2018
(This article belongs to the Section Water Use and Scarcity)

Abstract

:
The North China District (NCD) is one of the main grain production regions in China. The double cropping system of irrigation has been leading to the groundwater table decline at the speed of 1–2 m per year. Climate change leads to uncertainty surrounding the future of the NCD agricultural system, which will have great effects on crop yields and crop water demands. In this research, the Meteorological dataset from 54 weather station sites over the period 1960–2015 were collected to quantify the long-term spatial and temporal trends of meteorological data, including daily minimum temperature (Tmin), maximum temperature (Tmax), precipitation, solar radiation, reference evapotranspiration (ET0), and aridity index (AI). The results show that the long-term wheat and maize growing season and annual average air temperatures (Tmin and Tmax) showed strong north to south increasing trends throughout the NCD. The average annual precipitation was 632.9 mm across the NCD, more than 70% of which was concentrated in the maize growing season. The regional average annual ET0 was 1026.1 mm, which was 531.2 and 497.4 mm for the wheat and maize growing season, respectively. The regional precipitation decreased from northwest to southeast in each growing season and annual timescale. The funnel areas have lower precipitation and higher ET0 than the regional average, which may lead to the mining of the groundwater funnel area. The regional average annual AI is 0.63, which lies in the humid class. For temporal analysis, the regional average trends in annual Tmin, Tmax, solar radiation, ET0, precipitation, and AI were 0.37 °C/10a, 0.15 °C/10a, −0.28 MJ/day/m2/10a, −2.98 mm/10a, −12.04 mm/10a, and 0.005/10a, respectively. The increasing trend of temperature and the decreasing trend of solar radiation may have a negative effect on the regional food security. The funnel area AI showed a significant increasing trend for the winter wheat growing season and a decreasing trend for maize, which indicated that more irrigation will be needed for the maize growing season and the winter fallow policy may lead to the increasing trend precipitation being wasted. Analyzing the growing season and the annual meteorological elements of the spatiotemporal trends can help us better understand the influence of climate change on the natural resources and agricultural development in both the past and the future, and will provide us with invaluable information for the modification of cropping patterns to protect the regional and national water and food security.

1. Introduction

Global warming, becoming a real and pressing problem, is accepted widely as the main feature of climate change. The Fifth Assessment Report of the United Nations Intergovernmental Panel on Climate Change [1] showed that the period of 1983–2012 was the warmest 30 years over the last 800 years for the average annual Northern Hemisphere temperatures. Climate change has a great effect on the water demand of plants and crops, the water consumption, and the grain yield due to increased temperatures and variable rainfall, especially in poor adaptability areas [2,3]. Climate change is expected to increasingly affect crop yields [4,5,6,7,8] and crop water requirements [9]; which impact on the regional and global food security and water security.
The North China District (NCD) is one of the most important grain producing areas in China. It produces 7.31 × 107 ton of wheat and 5.73 × 107 ton of maize with 1.17 × 107 ha and 10.00 × 107 ha area, which was 56.7% of the nation’s wheat yields and 26.1% of its maize yields [10]. Thus, NCD plays an important part in ensuring China’s food security. The winter wheat-summer maize cropping system is the main cropping system in the North China District. This double cropping system needs about 800–1100 mm water every year [11], while the average annual precipitation is around 500–600 mm at the same time [12], thus, irrigation is needed in order to meet the optimal yield. However, the water resources are limited in this area and most of the irrigation water was from underground water, which causes a dramatic decline of the water table. In addition, irrigation is related to the distribution of the precipitation crop water requirement. Thus, assessing the spatiotemporal changes of the data is beneficial for optimizing irrigation during the wheat-maize cropping system and regional food and water security.
In the context of global warming, the trend of the temperature was found to have increased for most of the world [13,14,15]. While the solar radiation [16], precipitation [17,18], ET0 [19], and AI [20,21] were found to have increased or decreased in different regions of the world. Earlier regional studies that performed a trend analysis of the meteorological elements are more focused on one climate element in an annual timescale. However, crops are affected differently by the meteorological elements in different growing seasons. Therefore, assessing meteorological element trends both spatially and temporally in different growing seasons and with an annual timescale may help us better understand how climate change impacts regional food and water security. The aims of the present study are: (1) to analyze the spatial meteorological elements in different growing seasons and with annual timescale; (2) to calculate the meteorological element trends in different growing seasons and with an annual timescale; (3) to discuss the possible impacts of climate change on the regional water and food security of NCD.

2. Materials and Methods

2.1. Study Area Description

This study area consists of five provinces and cities (Beijing, Tianjin, Hebei, Shandong, Henan), which is a generally similar area to the North China Plain (Figure 1). The region is located between the longitude of 31 and 43° N and it is a monsoon climate with medium latitudes. The annual precipitation is around an average value of 500 to 600 mm [22]. In this study, each year was divided into two periods, comprising the winter wheat growing season and the summer maize growing season. Generally, the winter wheat growing season is from October to June and the summer maize growing season is from July to September.

2.2. Input Data Sources

The daily meteorological data from 54 weather station sites across the NCD contains the maximum air temperature, minimum air temperature, average relative humidity, precipitation, bright sunshine hours, and wind speed, which were obtained from the China Meteorological Data Sharing Service System. Most selected weather stations possess data from 1960 to 2015. The Mengjin, Qingdao, and Changdao weather stations possess data from 1961 to 2015, and the Weihai station, from 1940 to 2015 (Figure 1). The groundwater contiguous groundwater funnel area is collected from the Announcement of Outlaws Area, Restricted Area, and Prohibited Area of each province.

2.3. Calculation of ET0 and Aridity Index

Reference crop evapotranspiration (ET0) expresses the evaporating power of the atmosphere at a specific location and period. ET0 is a climate parameter and is only affected by climate parameters. There are many equations to estimate ET0, such as the Hargreaves Equation [23], Priestley-Taylor Equation [24], Penman-Monteith Equation [25], and FAO Penman-Monteith Equation [26]. In this study, we use the FAO Penman-Monteith Equation, which was most widely used to estimate the ET0 around the world. The FAO Penman-Monteith Equation can be expressed as:
ET 0 = 0.408 ( R n G ) + γ 900 T + 273 μ 2 ( e s e a ) + γ ( 1 + 0.34 μ 2 )
where Rn is the net radiation at the crop surface (MJ/m2/day), which was calculated from the sunshine hours using the equation developed by Rietveld [27]; G is the soil heat flux density (MJ/day/m2); T is the air temperature at a 2 m height (°C); μ2 is the wind speed at a 2 m height (m/s); es is the saturation vapor pressure (kPa); ea is the actual vapor pressure (kPa); Δ is the slope of the saturation vapor pressure-temperature curve (kPa/°C); and γ is the psychometric constant (kPa/°C). The ET0 was calculated by the SIMETAW model, which was developed based on the FAO Penman-Monteith equation [28], and has widely used in NCD [29,30,31,32].
Aridity Index (AI) is an index for assessing drought risk in the crop growing period, which considers the total precipitation and evapotranspiration. There are many methods to evaluate an AI. In this research, the AI was calculated according to the UNESCO [33] method, which can be computed as
Aridity   index   ( AI ) = P ET 0
where P is the precipitation (mm) that is obtained from each weather station site. The annual P and ET0 of the winter-wheat and summer-maize growing season is obtained by summing the daily ET0 rasters. Any site can be classified into the five categories of hyper-arid, arid, semi-arid, dry-subhumid, and humid according to the magnitude of the AI, which is classified by the numeric 0.03, 0.2, 0.5 and 0.65.

2.4. Mann-Kendall Trend Test and Slop Estimator

The Mann-Kendall (M-K) is one of the widely used non-parametric tests for detecting significant trends and slope in time series [34,35]. There are many methods to detect the significant trend of the considered meteorological elements in time series of the selected stations, the Mann-Kendall and Spearman’s Rho test are most widely used. Many theses have discussed the difference between the Mann-Kendall and Spearman’s Rho test [36,37,38], the results show a similar performance of these two methods. In this study we used the M-K method to detect the significant trend and Sen’s method to detect trend magnitudes, which were suggested by the World Meteorological Organization [39]. The advantage of the Mann-Kendall trend test is that it is not contaminated by a small number of outliers and the sample does not need to rank within a certain distribution, so it is widely used to estimate streamflow, precipitation, and temperature trend. The M-K test statistic S is calculated as
S = i = 1 n 1 j = i + 1 n sgn ( x j x i )
where n is the number of data points, xi and xj are the data values in the time series i and j (j > i), respectively, and sgn(xjxi) is the sign function as
sgn ( x j x i ) = { + 1 ,   if   x j x i > 0 0 ,         if   x j x i = 0 1 ,   if   x j x i > 0
The variance is computed as
Var ( S ) = n ( n 1 ) ( 2 n + 5 ) i = 1 m t i ( t i 1 ) ( 2 t i + 5 ) 18
where n is the number of data points, m is the number of tied groups, and ti denotes the number of ties of extent i. A tied group is a set of sample data having the same value. In cases where the sample size n > 10, the standard normal test statistic ZS is computed as
Z S = { S 1 V a r ( S ) ,   if   S > 0 0 ,                         if   S = 0 S 1 V a r ( S ) ,   if   S < 0
Positive values of ZS indicate increasing trends while negative ZS values show decreasing trends. The testing trends are done at the specific α significance level. When |ZS| > Z(1 − α⁄2), the null hypothesis is rejected and a significant trend exists in the time series. Z(1 − α⁄2) is obtained from the standard normal distribution table. In this study, a significance level of α = 0.05 was used. At the 5% significance level, the null hypothesis of no trend is rejected if |ZS| > 1.96.
Sen’s slope estimator [40] can be used for estimating time series data trend magnitudes. The advantage of this method is that it limits the influence of the outliers on the slope in comparison with linear regression. This estimate is given by the Theil-Sen Estimator as
b = Median ( x j x z j 1 )
where xj and xz are the data values at times j and z (j > z), respectively; b denotes the annual increment under the hypothesis of a linear trend. b provides the real slope of the tendency and can vary slightly from the slope obtained from linear regression. Sen’s slope estimator is widely used to estimate the magnitudes of trends in the climatic factors.
In this study, the M-K test is used to analyze the trends of air temperature, daily temperature range, solar radiation, precipitation, ET0, and AI in the NCD. The M-K trend test and slop estimator of each weather station site is calculated by MATLAB R2017b.

3. Results

3.1. Spatial and Temporal Variation of Temperature

3.1.1. Spatial Analysis

The temperature significantly influences crop yields and crop water use. Global warming has increased the instability of agriculture and the fluctuation of crop yields. The increase in temperature has delayed autumn frosts and led to later sowing and earlier harvesting of the winter-wheat growing season in the North China District [41]. Moreover, higher daily maxim temperatures have also caused heat and drought stress during the wheat and maize growing season, which highly negatively affected the crop yields [2,42,43,44]. On the other hand, higher daily minimum temperatures have accelerated the respiration of crops, which thus, have caused yield loss. Furthermore, temperature is one of the important input data for many hydrology and crop models. Thus, it is necessary to evaluate the regional temperature trend both spatially and temporally.
During 1960–2015, the average daily minimum temperature (Tmin) was lower in the north compared to the south through the wheat and maize growing season (Figure 2a–c). The figure indicates that the average daily Tmin of the maize growing season is higher than wheat, which was above zero over the study area. The wheat and maize growing season Tmin ranged from −9.0 to 6.6 °C and from 10.3 to 22.1 °C; and from −2.6 to 11.8 °C when averaged annually. The average daily Tmin of the Henan province is the highest, which is 4.6, 20.3, and 9.8 °C for wheat, the maize growing season, and annually; the Hebei province was the lowest (Table 1). Similar to Tmin, the Hebei and Henan provinces are the highest and lowest Tmax for the three periods. The average daily maximum air temperature (Tmax) of the study region has a south-to-north decreasing trend for the wheat growing season and annually, ranging from −2.6 to 11.5 °C and 3.7 to 16.2 °C. Unlike the wheat growing stage, the Tmax for maize has shown a middle-surround decreasing trend in the NPD, which ranged from 22.3 to 30.6 °C.

3.1.2. Temporal Analysis

The trends in the winter-wheat, summer-maize growing season, and the annual timescale from 1961 to 2015 for Tmin, and Tmax are shown in Figure 3. The maps also show the statistical significance of the weather station trends during the research period. The different color circles indicate the value of the temporal trends and the positive and negative numbers represent the increasing and decreasing trend. The black points in the maps indicate that the trend quantified for the weather station was statistically significant at α = 0.05, while the red rose point represents the significance of the increasing trends, which distinguishes them from the station with insignificant trends that are assigned with a white point. Table 2 shows the trend statistics for the growing seasons and the annual timescale Tmin and Tmin on a province-wide, funnel, and regional basis. Figure 3a–c show that most of the study regions show significant increasing trends in the wheat and maize growing seasons and the annual Tmin. Only one weather station shows a decreasing trend in the winter-wheat growing season and the annual timescale. It was higher in the wheat growing season than in the maize growing season, with the highest Tmin trend province being Beijing in the wheat growing season (Table 2). The regional average increasing trends in the wheat and maize growing seasons and the annual Tmin were 0.44, 0.24, and 0.37 °C/10a, respectively. Like Tmin, the trends for both the winter-wheat growing season and the annual timescale Tmax were positive, while the central parts of the NCD show non-significant increasing trends for the wheat growing season and the annual timescale (Figure 3d–f). On the other hand, the main parts of the Henan province show decreasing trends for the maize growing season, which was different to the Tmin. Table 2 shows that the average Tmax in the Henan province for the maize growing season is the only decreasing growing season Tmax, province-wide, which was −0.08 °C/10a.

3.2. Spatial and Temporal Variation of Solar Radiation

3.2.1. Spatial Analysis

Solar radiation is the most important source of renewable and environmental energy on the planet. Because of its abundance, solar energy can play a prominent role in our future energy and reduce the dependency on fossil fuel. In addition, many researchers have discovered the importance and positive effect of solar radiation on crop yield [45,46]. During 1960 to 2015, the average daily solar radiation was higher in the north compared to the south for the wheat and maize growing seasons and the annual timescale. It varied between 12.5 and 14.9 MJ/day/m2 in the wheat growing season (Figure 4a), which was lower than the maize growing season (Figure 4b). Table 3 shows the descriptive statistics for the average daily solar radiation. The Shandong province had the highest solar radiation for the wheat growing season from the five provinces. While, for maize growing season, the Tianjin and Hebei provinces were the highest. There was little difference in the funnel area and regional average daily solar radiation.

3.2.2. Temporal Analysis

Most of the NCD daily solar radiation shows a significant decreasing trend for the wheat growing season, except for the west part of the Henan province where the Qinling Mountains exist (Figure 5a). All weather station solar radiation showed significant decreasing trends for the maize growing season, which varied between −0.87 and −0.16 MJ/day/m2/10a (Figure 5b). All provincial solar radiation showed significant decreasing trends for the wheat and maize growing seasons and the annual timescale. The tendency of the maize growing season was higher than that of the wheat growing season and the funnel area was higher than the regional scale for both the wheat and maize growing seasons (Table 4). The annual solar radiation decreased during 1960 to 2015 at a rate of −0.35 and −0.28 MJ/day/m2/10a in the funnel area and the regional basis, respectively.

3.3. Spatial and Temporal Variation of Precipitation

3.3.1. Spatial Analysis

The quantification of the precipitation spatiotemporal variability in any region is important to enhance the capacity and level of the agriculture drought monitoring, prediction, and impact assessment, especially in arid and semi-arid regions such as NCD where the groundwater table declines by 1–2 m per year [47]. During the period of 1960–2015, the average precipitation decreased from the northwest to the southeast in each of the growing seasons and the annual timescale (Figure 6). The descriptive statistics of the winter-wheat and summer-maize growing seasons and the annual precipitation for the province-wide area are presented in Table 5. Regionally, the wheat growing season precipitation varies from 83.2 to 461.6 mm, which, for the maize growing season and the annual timescale, was from 297.2 to 609.0 mm and 380.6 to 1070.0 mm, respectively. Averaged by province, the valley wheat and maize growing season precipitations were observed in Beijing and Hebei, respectively, while the maximum values were observed in Henan and Shandong. The minimum value annual precipitation was observed in Hebei too, which was 520.7 mm, while the peak annual precipitation was observed in the Henan province. The regional mean precipitation was 632.9 mm, while for the funnel area, it was 570.4 mm.

3.3.2. Temporal Analysis

The trends in the winter-wheat and summer-maize growing seasons and the annual precipitation have been computed for the period 1960–2015 using a non-parametric test (Figure 7). The winter-wheat growing season precipitation trends are positive for the study area and more than 44.4% of the weather station showed a significant increase. As shown in Figure 7a, the trend of the precipitation was higher in the east of the NCD than the west. Unlike the wheat growing season precipitation trends, the maize growing season precipitation trends show negative trends in the majority of the NCD and there were only 9 weather stations showing significant decreasing precipitation trends. It is noteworthy that Figure 7b shows that most increased weather station was concentration in the Henan province. Figure 7c shows the nature and significance of the annual precipitation trends over the study area. There were 31 weather stations showing negative trends in annual precipitation, which is a total proportion of 57.4%. Most of the weather station precipitations showed non-significant trends in the study region, with only one weather station showing dominant positive trends, and two showing dominant negative trends. The descriptive statistics for the precipitation trend when averaged for province-wide, funnel area, and regional are listed in Table 6. The magnitude of the wheat growing season precipitation trends was calculated as 6.94 mm per decade on the funnel area average basis, which was higher than the study region. Additionally, the precipitation trend was lower in the funnel area than in the regional basis for the maize growing season, which were −14.72 and −9.63 mm/10a, respectively. All the provinces showed a positive trend for the wheat growing season and a negative trend for the maize growing season, with the Henan province being the only one out of the five provinces to show a positive trend for the annual precipitation. Table 6 shows that the annual precipitation decreased during the period of 1960–2015 at a rate of −5.69 and −2.98 mm/10a in the funnel area and the regional basis, respectively.

3.4. Spatial and Temporal Variation of Reference Evapotranspiration

3.4.1. Spatial Analysis

ET0 is an important part of the hydrological cycle and occupies a more important proportion in the regional water resource balance. Different from the precipitation effects of the regional water resource supply, the ET0 mainly affects consumption, especially the crop water requirement. The results show that the mean winter-wheat growing season ET0 over a 56-year period varied from 420.3 mm to 670.0 mm across the NCD (Figure 8a). Averaged by province, the valley wheat growing season ET0 was observed in the Hebei province, which was 508.6 mm (Table 7), while the maximum values were observed in Beijing. The funnel area and regional ET0 were 545.9 and 531.2 mm for the wheat growing season. Figure 8b shows that the central part of the NCD ET0 was higher than the surrounding parts for the maize growing season. The average funnel area and regional maize growing season ET0 were 513.8 and 497.4 mm, respectively, which was lower than that of the wheat growing season. The Jinan weather station had the maximum annual ET0, which was 1240.9 mm. The annual ET0 from the lowest to the highest, province-wide, are as follows: Hebei, Henan, Shandong, Beijing, and Tianjin.

3.4.2. Temporal Analysis

The 61-year average of the ET0 trend in the research region and funnel area was about −0.53 and −3.60 mm/10a for the heat growing season (Table 8). The fluctuation range for the wheat growing season is from −15.62 to 15.61 mm/10a (Figure 9a). There are 31 weather station ET0 that show negative trends for the wheat growing season and 15 of them show significant decreasing trends. Hebei is the only one of the five provinces that shows positive trends for the winter-wheat growing season ET0. The ET0 trends for the maize growing season show northeast to southwest decreasing trends. Only five weather station ET0 show positive trends for the maize growing season, which were all located in the east part of Shandong and the north part of Hebei. More than 68% of the weather stations show significant decreasing trends in the maize growing season. All five provincial (cities) ET0 show decreasing trends in the maize growing season. The annual ET0 trend ranges from −37.70 to 19.50 mm/10a in the study area. The Tianjin, Hebei, Henan, and Shandong average ET0 trends show negative trends for the annual timescale, while Beijing shows positive trends. For the annual ET0, the funnel area and regional average trends were −17.49 and −12.04 mm/10a, respectively.

3.5. Spatial and Temporal Variation of Arid Index

3.5.1. Spatial Analysis

AI is an index that combines precipitation with ET0, which can present information about drought or wetness in each district. AI is one of the most important indexes in meteorology, ecology, and agronomy, which is a widely used measure all over the world [48,49,50]. In this research, the aridity index was calculated at point locations for the wheat and maize growing seasons and the annual timescale. These calculations are presented in Figure 10. There were southeast-to-northwest decreasing trends in the AI. The descriptive statistics for the AI when averaged for the province-wide, funnel area, and regional bases are presented in Table 9. The AI for the wheat growing season had a great difference province-wide. The maximum AI was observed in Henan, while the valley AI occurred in Beijing. The average regional and funnel area wheat growing season AI are 0.26 and 0.34, respectively, which means that the funnel area and the study area lie in the semi-arid class. For the maize growing season, the AI ranges from 0.55 to 1.30 (Figure 10b). The total study region falls into the humid class for the maize growing season, except the northwest part of the Hebei province. The regional average annual AI is 0.63, which lies in the humid class. This does not mean that the region is dominantly humid, since we can find, in Figure 10c, that more than half of the study region lies in the semi-arid and dry sub-humid classes.

3.5.2. Temporal Analysis

The trends in the wheat and maize growing season AIs and the annual AI were calculated for each weather station in the study region. Figure 11 shows the wheat and maize growing season AIs and the annual AI by the M-K test for each weather station between 1960 and 2015. The wheat growing season AI has decreased over most of the study region, with 49 weather station AI showing an increasing trend in the wheat growing season, 21 of which were significantly increasing. For the maize growing season AI, about 57% of the weather stations had negative trends and most of the weather stations are concentrated in the west part of Shandong and the north part of the Jing-Jin-Ji region. The average regional wheat and maize growing season AI trends are 0.013 and −0.004 /10a. There are 7 weather station AIs that show increasing trends for the annual timescale and all of them are in the Henan province. So, the only province with increasing trends in annual AI was observed in Henan (Table 10). The funnel area and regional annual AI show slight decreasing trends, which were 0.004 and 0.005 /10a, respectively.

4. Discussion

4.1. Perspectives and Implications

Due to the rising population, dietary changes, and increasing biofuel consumption, crop production must be doubled to meet the coming global food demand [51]. On the other hand, for most of the world, groundwater must be pumped for irrigation purposes in order to maintain a high yield, especially in semi-arid and arid regions [52]. Irrigation has led to the rapid decline of the groundwater table in highly intensified agricultural parts of the world, such as the North China Plain, the Sanjiang Plain, the Central Valley of California, and the Indus River Basin [52,53]. Water shortages have become the main factor limiting grain production in the world, as well as NCD, since the cropping system has changed from one rainfed crop to the winter wheat-summer maize system of irrigation and since the NPD’s grain production has increased and become one of the major grain producing areas in China, while the groundwater table was declined rapidly at the same time [47,54]. Meteorological elements have strong effects on crop yields and water consumption [2] and agricultural water accounts for 62.4% of the total water usage in the NCD [55]. Thus, evaluating meteorological element changes both spatially and temporally is meaningful for us to understand the possible effects of climate change on the NCD’s water and food security.

4.2. Influence of Climate Factors on Regional Water Security

The current results showed that the spatial variations of ET0 and precipitation were quite different in the wheat and maize growing seasons and the annual timescale. The average precipitation decreasing trend from the northwest to southeast and the average ET0 show decreasing trends from the central part of the NCD to the surroundings, which lead to the funnel area AI being lower than other parts of the NCD. Similar results have been found by Fan et al. [56]. Precipitation and irrigation are the only two ways to meet the crop water requirement in the NCD. Lesser precipitation and more evapotranspiration meant that more irrigation from groundwater is needed for agriculture. There was one cropping system in the north part of the Hebei province because of the low temperature in winter, while the other was a double cropping system. The double cropping system and lesser precipitation led to the agriculture of the funnel area needing more groundwater for irrigation in order to maintain a high yield level. Over pumping the region’s groundwater resources over the past three decades may have caused a cone of depression in the groundwater table (funnel area).
Our results show that the regional mean annual precipitation was 632.9 mm and more than 70% occurred during the summer maize growing season, while there was little difference for the ET0 between the wheat growing season and the maize growing season. This means that more irrigation was needed in the wheat growing season. Rainfall can meet the maize water requirement for the most part of the NCD [57]. In order to protect the groundwater resources, sustainable farmland use must be maintained and agricultural sustainable development must be forwarded. China has been developing a reasonable subsidy policy of cultivated land fallow by selecting some areas as pilot areas in the funnel area since 2016. The fallow policy changed double cropping back to one rainfed crop system with a winter fallow, which will greatly reduce the groundwater resource consumption in the funnel area. This study indicates that the precipitation for the wheat growing season showed an increasing trend in the funnel area (Table 6), while the ET0 showed a decreasing trend (Table 8). More precipitation and fewer water requirements mean that less irrigation was needed for the wheat growing season. The precipitation and the ET0 for the maize growing season both showed decreasing trends in the funnel area and the AI showed decreasing trends as well. This indicates that more groundwater must be pumped to maintain an optimal yield for the maize growing season. If the climate changing trend in the future was the same as in the past, the decreasing AI will enhance the water deficit and drought risk for the rainfed maize. Though the precipitation for wheat shows a positive trend, the winter fallow policy will waste more precipitation as weed evapotranspiration. The annual precipitation showed a slightly decreasing trend regional, while the ET0 showed a decreasing trend too, so the annual AI showed a slightly positive trend. Discarding the fallow policy, climate change has had positive impacts on the regional water security. On the other hand, the South-to-North Water Diversion Project will ultimately transfer 44.8 Gm3 water from the Yangtze River Basin to NCD annually, which will greatly release the water stress of the NCD.

4.3. Effects of Climate Factors on Regional and National Food Security

The crop yield is mainly affected by solar radiation and temperature. A lot of evidence is available, showing the importance and the positive effect of solar radiation on the grain yield [9,58,59,60]. Stanhill [61] reported a globally average solar radiation reduction of 0.43 MJ/day/m2 per decade, which was lower than the NCD average reduction. The negative trend of solar radiation has been and will be an increasing challenge for NCD’s crop production.
The observed changes in temperature in NCD were in line with other research (Fan et al., 2016; Yuan et al., 2010; Tao and Zhao, 2010). Zhao et al. [62] reported that the global production of wheat declined by 6.0% and maize by 7.4% for each 1 °C increase in the global mean temperature. Peng et al. [45] suggested that each degree-Celsius increase in the growing season minimum temperature would reduce 10% of the grain yield. The current results indicated that both Tmin showed significant increasing trends for the wheat and maize growing seasons in most of the NCD (Figure 3). Table 2 showed that the regional wheat Tmin trend was 0.44 °C/10a, which was a little higher than the value of 0.037 °C/a calculated by Fan et al. [56] for the whole nation. The trend of Tmin leads the wheat yield of the NCD to decline by about 0.44% per year by the Peng’s [39] results, discarding agronomic (variety and management practices) factor changes that occurred. When the Tmax was under 30 °C, the higher the Tmax, the more yield obtained. Once the Tmax > 30 °C, which was defined as a high temperature by Schauberger [44], the crop enzymes and tissues are damaged [63,64] and the temperature impairs flowering [65], leads to precocious maturity and senescence [66], and triggers oxidative stress [67]. So, crop yields decline rapidly due to high temperatures, especially under rainfed conditions [4,6,44]. Figure 2e showed that the average Tmax for the maize growing season was over 30 °C in the central part of the NCD, which indicates that Tmax has a negative impact on the maize yield in this area. The best temperature for wheat anthesis and grain filling ranges from 12 to 22 °C [68]. The wheat growing season Tmax was shown to have an increasing trend for all parts of the NCD, which can lead to the early maturation of wheat and the significant reduction of the wheat yield [69,70].
The increasing temperatures have led to the north boundary of the two-crop area moving north [41,71], which means that the two-crop area has increased in the NCD. This would relieve the pressures of food security of the NCD, but, on the other hand, it would increase the area that over-exploits the groundwater of the NCD. In addition, the fallow policy has reduced the agricultural groundwater consumption by reducing the winter wheat planted area. To maintain the province-wide crop production, the non-funnel area has changed from the cotton to the winter wheat-summer maize cropping system [72]. So, the wheat planted area did not change very much for the whole NCD, while the wheat planted area has increased in the non-funnel area, which will decrease the groundwater table of the non-funnel area.

5. Conclusions

The long-term wheat and maize growing season and the annual average air temperatures (Tmin, Tmax) showed strong north to south increasing trends throughout the NCD, with exceptionally lower temperatures observed in the west part of the Henan province with high elevation. The average daily solar radiation was higher in the north compared with that of the south for the wheat and maize growing seasons and the annual timescale, which had ranges of 12.5–14.9, 16.7–19.6, and 14.3–16.5 MJ/day/m2, respectively. The average precipitation decreased from the northwest to the southeast in each growing season and annual timescale. The funnel area annual precipitation, mainly located in the Hebei province, was 570.4 mm, which was lower than the regional average precipitation. While the funnel ET0 was higher than the regional average ET0. Therefore, the AI, as a result of the magnitudes of both the precipitation and the ET0, was lower in the funnel area than in the regional area. The high ET0, low rainfall, and double cropping system may cause the world’s largest mining of the groundwater funnel area. The historical spatial datasets were constructed, presented, and analyzed in this research can aid water and environmental personnel, water resource researchers, and policymaking agencies better understand the variations of the NCD climate and provide a theoretical basis for it. The trend for the wheat growing season was larger than that of the maize growing season. The solar radiation was shown to have a decreasing trend for both the wheat and maize growing seasons. The increasing trend of Tmin and the decreasing trend of the solar radiation may reduce crop production. The funnel area AI showed increasing trends for the winter wheat growing season and decreasing trends for the maize growing season. Though the winter fallow policy reduced the groundwater pumping, the increasing precipitation trend of the winter wheat growing season would cause a wastage of water by soil water evaporation. Climate change has had positive impacts on the NCD’s water security, especially for the funnel area, as shown by the increasing trend of the annual AI.

Author Contributions

Y.Y.; L.P. and Y.J. collected the data. J.T. and J.L. analyzed the data; J.T. and X.Y. wrote the paper. X.W. and F.C. edited the paper.

Funding

This research was funded by [National Key Research and Development Program of China] grant number [2016YFD0300201].

Acknowledgments

This study was financially supported by the National Key Research and Development Program of China (2016YFD0300201). We thank the anonymous reviewers and the Editor for their suggestions, which substantially improved the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Intergovernmental Panel on Climate Change (IPCC). Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK, 2013. [Google Scholar]
  2. Rosenzweig, C.; Elliott, J.; Deryng, D.; Ruane, A.C.; Müller, C.; Arneth, A.; Boote, K.J.; Folberth, C.; Glotter, M.; Khabarov, N.; et al. Assessing Agricultural Risks of Climate Change in the 21st Century in a Global Gridded Crop Model Intercomparison. Proc. Natl. Acad. Sci. USA 2014, 111, 3268–3273. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Wheeler, T.; Von Braun, J. Climate Change Impacts on Global Food Security. Science 2013, 341, 508–513. [Google Scholar] [CrossRef] [PubMed]
  4. Challinor, A.J.; Watson, J.; Lobell, D.B.; Howden, S.M.; Smith, D.R.; Chhetri, N. A Meta-Analysis of Crop Yield under Climate Change and Adaptation. Nat. Clim. Chang. 2014, 4, 287–291. [Google Scholar] [CrossRef]
  5. Kang, Y.; Khan, S.; Ma, X. Climate Change Impacts on Crop Yield, Crop Water Productivity and Food Security—A Review. Prog. Nat. Sci. 2009, 19, 1665–1674. [Google Scholar] [CrossRef]
  6. O’Leary, G.J.; Aggarwal, P.K.; Calderini, D.F.; Connor, D.J.; Craufurd, P.; Eigenbrode, S.D.; Han, X.; Hatfield, J.L. Challenges and Responses to Ongoing and Projected Climate Change for Dryland Cereal Production Systems throughout the World. Agronomy 2018, 8, 34. [Google Scholar] [CrossRef]
  7. Holzkämper, A. Adapting Agricultural Production Systems to Climate Change—What’s the Use of Models? Agriculture 2017, 7, 86. [Google Scholar] [CrossRef]
  8. Schlenker, W.; Roberts, M.J. Nonlinear Temperature Effects Indicate Severe Damages to U.S. Crop Yields under Climate Change. Proc. Natl. Acad. Sci. USA 2009, 106, 15594–15598. [Google Scholar] [CrossRef] [PubMed]
  9. Fischer, G.; Tubiello, F.N.; Van Velthuizen, H.; Wiberg, D.A. Climate Change Impacts on Irrigation Water Requirements: Effects of Mitigation, 1990–2080. Technol. Forecast. Soc. Chang. 2007, 74, 1083–1107. [Google Scholar] [CrossRef] [Green Version]
  10. National Bureau of Statistics of China. China Statistical Yearbook 2017; China Statistics Press: Beijing, China, 2017. [Google Scholar]
  11. Sun, H.; Zhang, X.; Wang, E.; Chen, S.; Shao, L. Quantifying the Impact of Irrigation on Groundwater Reserve and Crop Production—A Case Study in the North China Plain. Eur. J. Agron. 2015, 70, 48–56. [Google Scholar] [CrossRef]
  12. Yang, X.L.; Chen, Y.Q.; Steenhuis, T.S.; Pacenka, S.; Gao, W.S.; Ma, L.; Zhang, M.; Sui, P. Mitigating Groundwater Depletion in North China Plain with Cropping System That Alternate Deep and Shallow Rooted Crops. Front. Plant Sci. 2017, 8, 980. [Google Scholar] [CrossRef] [PubMed]
  13. Bale, J.S.; Masters, G.J.; Hodkinson, I.D.; Awmack, C.; Bezemer, T.M.; Brown, V.K.; Butterfield, J.; Buse, A.; Coulson, J.C.; Farrar, J.; et al. Herbivory in Global Climate Change Research: Direct Effects of Rising Temperature on Insect Herbivores. Glob. Chang. Biol. 2002, 8, 1–16. [Google Scholar] [CrossRef]
  14. Dai, A.; Fyfe, J.C.; Xie, S.P.; Dai, X. Decadal Modulation of Global Surface Temperature by Internal Climate Variability. Nat. Clim. Chang. 2015, 5, 555–559. [Google Scholar] [CrossRef]
  15. Morice, C.P.; Kennedy, J.J.; Rayner, N.A.; Jones, P.D. Quantifying Uncertainties in Global and Regional Temperature Change Using an Ensemble of Observational Estimates: The Hadcrut4 Data Set. J. Geophys. Res. Atmos. 2012, 117. [Google Scholar] [CrossRef]
  16. Liepert, B.G. Observed Reductions of Surface Solar Radiation at Sites in the United States and Worldwide from 1961 to 1990. Geophys. Res. Lett. 2002, 29, 61-1–61-4. [Google Scholar] [CrossRef]
  17. Alexander, L.V.; Zhang, X.; Peterson, T.C.; Caesar, J.; Gleason, B.; Klein Tank, A.M.G.; Haylock, M.; Collins, D.; Trewin, B.; Rahimzadeh, F.; et al. Global Observed Changes in Daily Climate Extremes of Temperature and Precipitation. J. Geophys. Res. Atmos. 2006, 111. [Google Scholar] [CrossRef]
  18. Arnell, N.W. Climate Change and Global Water Resources. Glob. Environ. Chang. 1999, 9, S31–S49. [Google Scholar] [CrossRef]
  19. Gorelick, S.M.; Zheng, C. Global Change and the Groundwater Management Challenge. Water Resour. Res. 2015, 51, 3031–3051. [Google Scholar] [CrossRef]
  20. Huang, J.; Yu, H.; Guan, X.; Wang, G.; Guo, R. Accelerated Dryland Expansion under Climate change. Nat. Clim. Chang. 2015, 6, 166–171. [Google Scholar] [CrossRef]
  21. Zarch, M.A.A.; Sivakumar, B.; Sharma, A. Assessment of Global Aridity Change. J. Hydrol. 2015, 520, 300–313. [Google Scholar] [CrossRef]
  22. Cao, G.; Han, D.; Song, X. Evaluating Actual Evapotranspiration and Impacts of Groundwater Storage Change in the North China Plain. Hydrol. Process. 2014, 28, 1797–1808. [Google Scholar] [CrossRef]
  23. Hargreaves, G.H.; Samani, Z.A. Estimating Potential Evapotranspiration. J. Irrig. Drain. Div. 1982, 108, 225–230. [Google Scholar]
  24. Priestley, C.H.B.; Taylor, R.J. On the Assessment of Surface Heat Flux and Evaporation Using Large-Scale Parameters. Mon. Weather Rev. 1972, 100, 81–92. [Google Scholar] [CrossRef] [Green Version]
  25. Penman, H.L. Natural Evaporation from Open Water, Bare Soil and Grass. Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 1948, 193, 120–145. [Google Scholar] [CrossRef]
  26. Frere, M.; Popov, G.F. Agrometeorological Crop Monitoring and Forecasting; FAO: Rome, Italy, 1979. [Google Scholar]
  27. Rietveld, M.R. A New Method for Estimating the Regression Coefficients in the Formula Relating Solar Radiation to Sunshine. Agric. Meteorol. 1978, 19, 243–252. [Google Scholar] [CrossRef]
  28. Snyder, R.L.; Orang, M.N.; Matyac, S.G.J.S.; Sarreshteh, S. Simetaw (Simulation of Evapotranspiration of Applied Water). In California Water Plan Update; Department of Water Resources: Sacramento, CA, USA, 2005. [Google Scholar]
  29. Yang, X.; Gao, W.; Shi, Q.; Chen, F.; Chu, Q. Impact of Climate Change on the Water Requirement of Summer Maize in the Huang-Huai-Hai Farming Region. Agric. Water Manag. 2013, 124, 20–27. [Google Scholar] [CrossRef]
  30. Yang, X.L.; Song, Z.W.; Wang, H.; Shi, Q.H.; Chen, F.; Chu, Q.Q. Spatio-Temporal Variations of Winter-Wheat Water Requirement and Climatic Causes in Huang-Huai-Hai Farming Region. Chin. J. Eco-Agric. 2012, 20, 356–362. [Google Scholar] [CrossRef]
  31. Song, Z.W. Construction of the Water-Saving Cropping Patterns Based on Simulation Model and Scenario Analysis: A Case of Beijing; China Agricultural University: Beijing, China, 2009. [Google Scholar]
  32. Huang, J.; Xu, C.; Ridoutt, B.G.; Chen, F. Reducing Agricultural Water Footprints at the Farm Scale: A Case Study in the Beijing Region. Water 2015, 7, 7066–7077. [Google Scholar] [CrossRef] [Green Version]
  33. United Nations Educational, Scientific and Cultural Organization (UNESCO). Map of the World Distribution of Arid Regions: Map at Scale 1:25,000,000 with Explanatory Note; MAB Technical Notes; UNESCO: Paris, France, 1979. [Google Scholar]
  34. Mann, H.B. Nonparametric Tests against Trend. Econometrica 1945, 13, 245–259. [Google Scholar] [CrossRef]
  35. Kendall, M.G.; Stuart, A. The Advanced Theory of Statistics; Griffin: London, UK, 1946; Volume 3. [Google Scholar]
  36. Shadmani, M.; Marofi, S.; Roknian, M. Trend analysis in reference evapotranspiration using Mann-Kendall and Spearman’s Rho tests in arid regions of Iran. Water Resour. Manag. 2012, 26, 211–224. [Google Scholar] [CrossRef]
  37. Yue, S.; Pilon, P.; Cavadias, G. Power of the Mann–Kendall and Spearman’s rho tests for detecting monotonic trends in hydrological series. J. Hydrol. 2002, 259, 254–271. [Google Scholar] [CrossRef]
  38. Ahmad, I.; Tang, D.; Wang, T.; Wang, M.; Wagan, B. Precipitation trends over time using Mann-Kendall and spearman’s rho tests in swat river basin, Pakistan. Adv. Meteorol. 2015, 431860. [Google Scholar] [CrossRef]
  39. Kundzewicz, Z.W.; Robson, A.J. Detecting Trend and other Changes in Hydrological Data. World Climate Programme—Water, World Climate Programme Data and Monitoring, WCDMP-45, WMO/TD no. 1013; World Meteorological Organization: Geneva, Switzerland, 2000. [Google Scholar]
  40. Sen, P.K. Estimates of the Regression Coefficient Based on Kendall’s Tau. J. Am. Stat. Assoc. 1968, 63, 1379–1389. [Google Scholar] [CrossRef]
  41. Li, W.J.; Liang, H.; Wang, P.J. Effects of Climate Warming on the Planting Boundary and Developmental Stages of Winter Wheat. J. Triticeae Crops 2013, 33, 382–388. [Google Scholar]
  42. Wang, Z.B.; Wang, M.; Yin, X.G.; Zhang, H.L.; Chu, Q.Q.; Wen, X.Y.; Chen, F. Spatiotemporal Change Characteristics of Heat and Rainfall during the Growth Period of Summer Maize in North China Plain from 1961 to 2010. Chin. J. Eco-Agric. 2015, 23, 473–481. [Google Scholar]
  43. Wang, Z.B.; Wang, M.; Yin, X.G.; Zhang, H.L.; Chu, Q.Q. Spatiotemporal Change Characteristics of Heat and Rainfall during the Growth Period of Winter Wheat in North China Plain from 1961 to 2010. J. China Agric. Univ. 2015, 20, 16–23. [Google Scholar]
  44. 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]
  45. Peng, S.; Huang, J.; Sheehy, J.E.; Laza, R.C.; Visperas, R.M.; Zhong, X.; Centeno, G.S.; Khush, G.S.; Cassman, K.G. Rice Yields Decline with Higher Night Temperature from Global Warming. Proc. Natl. Acad. Sci. USA 2004, 101, 9971–9975. [Google Scholar] [CrossRef] [PubMed]
  46. Wang, Z.; Xie, P.; Lai, C.; Chen, X.; Wu, X.; Zeng, Z.; Li, J. Spatiotemporal Variability of Reference Evapotranspiration and Contributing Climatic Factors in China during 1961–2013. J. Hydrol. 2017, 544, 97–108. [Google Scholar] [CrossRef]
  47. Kong, X.; Zhang, X.; Lal, R.; Zhang, F.; Chen, X.; Niu, Z.; Han, L.; Song, W. Groundwater Depletion by Agricultural Intensification in China’s Hhh Plains, since 1980s. In Advances in Agronomy; Academic Press: Cambridge, MA, USA, 2015. [Google Scholar]
  48. 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. 2012, 57, 275–285. [Google Scholar] [CrossRef] [PubMed]
  49. Nastos, P.T.; Politi, N.; Kapsomenakis, J. Spatial and Temporal Variability of the Aridity Index in Greece. Atmos. Res. 2013, 119, 140–152. [Google Scholar] [CrossRef]
  50. Yin, X.; Olesen, J.E.; Wang, M.; Kersebaum, K.C.; Chen, H.; Baby, S.; Öztürk, I.; Chen, F. Adapting Maize Production to Drought in the Northeast Farming Region of China. Eur. J. Agron. 2016, 77, 47–58. [Google Scholar] [CrossRef]
  51. Ray, D.K.; Mueller, N.D.; West, P.C.; Foley, J.A. Yield Trends Are Insufficient to Double Global Crop Production by 2050. PLoS ONE 2013, 8, e66428. [Google Scholar] [CrossRef] [PubMed]
  52. Aeschbach-Hertig, W.; Gleeson, T. Regional Strategies for the Accelerating Global Problem of Groundwater Depletion. Nat. Geosci. 2012, 5, 853–861. [Google Scholar] [CrossRef]
  53. Rodell, M.; Velicogna, I.; Famiglietti, J.S. Satellite-Based Estimates of Groundwater Depletion in India. Nature 2009, 460, 999–1002. [Google Scholar] [CrossRef] [PubMed]
  54. Liu, Y.; Wang, E.; Yang, X.; Wang, J. Contributions of Climatic and Crop Varietal Changes to Crop Production in the North China Plain, since 1980s. Glob. Chang. Biol. 2009, 16, 2287–2299. [Google Scholar] [CrossRef]
  55. Bureau of Provincial Water Resources. Water Resource Bulletin; China Water Power Press: Beijing, China, 2014. [Google Scholar]
  56. Fan, J.; Wu, L.; Zhang, F.; Xiang, Y.; Zheng, J. Climate Change Effects on Reference Crop Evapotranspiration across Different Climatic Zones of China during 1956–2015. J. Hydrol. 2016, 542, 923–937. [Google Scholar] [CrossRef]
  57. Zhang, X.; Chen, S.; Sun, H.; Shao, L.; Wang, Y. Changes in Evapotranspiration over Irrigated Winter Wheat and Maize in North China Plain over Three Decades. Agric. Water Manag. 2011, 98, 1097–1104. [Google Scholar] [CrossRef]
  58. Muchow, R.C.; Sinclair, T.R.; Bennett, J.M. Temperature and Solar Radiation Effects on Potential Maize Yield across Locations. Agron. J. 1990, 82, 338–343. [Google Scholar] [CrossRef]
  59. Cossu, M.; Murgia, L.; Ledda, L.; Deligios, P.A.; Sirigu, A.; Chessa, F.; Pazzona, A. Solar Radiation Distribution inside a Greenhouse with South-Oriented Photovoltaic Roofs and Effects on Crop Productivity. Appl. Energy 2014, 133, 89–100. [Google Scholar] [CrossRef]
  60. Wang, J.; Dong, J.; Yi, Y.; Lu, G.; Oyler, J.; Smith, W.K.; Zhao, M.; Liu, J.; Running, S. Decreasing Net Primary Production Due to Drought and Slight Decreases in Solar Radiation in China from 2000 to 2012. J. Geophys. Res. Biogeosci. 2017, 122, 261–278. [Google Scholar] [CrossRef]
  61. Stanhill, G.; Cohen, S. Global Dimming: A Review of the Evidence for a Widespread and Significant Reduction in Global Radiation with Discussion of Its Probable Causes and Possible Agricultural Consequences. Agric. For. Meteorol. 2001, 107, 255–278. [Google Scholar] [CrossRef]
  62. 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] [PubMed]
  63. Farooq, M.; Wahid, A.; Kobayashi, N.; Fujita, D.; Basra, S.M.A. Plant Drought Stress: Effects, Mechanisms and Management. In Sustainable Agriculture; Springer: Dordrecht, The Netherlands, 2009. [Google Scholar]
  64. Hasanuzzaman, M.; Nahar, K.; Alam, M.M.; Roychowdhury, R.; Fujita, M. Physiological, Biochemical, and Molecular Mechanisms of Heat Stress Tolerance in Plants. Int. J. Mol. Sci. 2013, 14, 9643–9684. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Boyer, J.S.; Westgate, M.E. Grain Yields with Limited Water. J. Exp. Bot. 2004, 55, 2385–2394. [Google Scholar] [CrossRef] [PubMed]
  66. Folberth, C.; Yang, H.; Gaiser, T.; Liu, J.; Wang, X.; Williams, J.; Schulin, R. Effects of Ecological and Conventional Agricultural Intensification Practices on Maize Yields in Sub-Saharan Africa under Potential Climate Change. Environ. Res. Lett. 2014, 9, 044004. [Google Scholar] [CrossRef]
  67. Kadam, N.N.; Xiao, G.; Melgar, R.J.; Bahuguna, R.N.; Quinones, C.; Tamilselvan, A.; Prasad, P.V.V.; Jagadish, K.S. Chapter Three—Agronomic and Physiological Responses to High Temperature, Drought, and Elevated CO2 Interactions in Cereals. Adv. Agron. 2014, 127, 111–156. [Google Scholar]
  68. Farooq, M.; Bramley, H.; Palta, J.A.; Siddique, K.H. Heat Stress in Wheat during Reproductive and Grain-Filling Phases. Crit. Rev. Plant Sci. 2011, 30, 491–507. [Google Scholar] [CrossRef]
  69. Tewolde, H.; Fernandez, C.J.; Erickson, C.A. Wheat Cultivars Adapted to Post-Heading High Temperature Stress. J. Agron. Crop Sci. 2006, 192, 111–120. [Google Scholar] [CrossRef]
  70. Mondal, S.; Joshi, A.K.; Huerta-Espino, J.; Singh, R.P. Early Maturity in Wheat for Adaptation to High Temperature Stress. In Advances in Wheat Genetics: From Genome to Field; Springer: Tokyo, Japan, 2015. [Google Scholar]
  71. Li, K.N.; Yang, X.G.; Liu, Z.J.; Wang, W.F.; Chen, F. Analysis of the Potential Influence of Global Climate Change on Cropping Systems in China III. The Change Characteristics of Climatic Resources in Northern China and Its Potential Influence on Cropping Systems. Sci. Agric. Sin. 2010, 43, 2088–2097. [Google Scholar]
  72. Hebei Province Bureau of Statistics. Hebei Statistical Yearbook 2016; China Statistics Press: Beijing, China, 2016. [Google Scholar]
Figure 1. Location of North China District and the weather station site.
Figure 1. Location of North China District and the weather station site.
Water 10 00789 g001
Figure 2. The spatial distribution of (ac) long-term average minimum air temperature (Tmin) and (df) maximum air temperature (Tmax) on the winter-wheat growing season, summer-maize growing season, and annual basis.
Figure 2. The spatial distribution of (ac) long-term average minimum air temperature (Tmin) and (df) maximum air temperature (Tmax) on the winter-wheat growing season, summer-maize growing season, and annual basis.
Water 10 00789 g002
Figure 3. The spatial distribution and statistical significance of the trends in (ac) Tmin and (df) Tmax on the winter-wheat growing-season, the summer-maize growing-season, and the annual basis.
Figure 3. The spatial distribution and statistical significance of the trends in (ac) Tmin and (df) Tmax on the winter-wheat growing-season, the summer-maize growing-season, and the annual basis.
Water 10 00789 g003
Figure 4. The spatial distribution of the long-term average daily solar radiation (R) for (a) the winter-wheat growing season; (b) the summer-maize growing season and (c) the annual basis.
Figure 4. The spatial distribution of the long-term average daily solar radiation (R) for (a) the winter-wheat growing season; (b) the summer-maize growing season and (c) the annual basis.
Water 10 00789 g004
Figure 5. The spatial distribution and statistical significance of the trends in the daily solar radiation (R) on (a) the winter-wheat growing season; (b) the summer-maize growing season and (c) annual basis.
Figure 5. The spatial distribution and statistical significance of the trends in the daily solar radiation (R) on (a) the winter-wheat growing season; (b) the summer-maize growing season and (c) annual basis.
Water 10 00789 g005
Figure 6. The spatial distribution of the long-term average precipitation (P) for (a) the winter-wheat growing season; (b) the summer-maize growing season and (c) the annual basis.
Figure 6. The spatial distribution of the long-term average precipitation (P) for (a) the winter-wheat growing season; (b) the summer-maize growing season and (c) the annual basis.
Water 10 00789 g006
Figure 7. The spatial distribution and statistical significance of the trends in the precipitation (P) of the (a) winter-wheat growing season; (b) the summer-maize growing season and (c) the annual basis.
Figure 7. The spatial distribution and statistical significance of the trends in the precipitation (P) of the (a) winter-wheat growing season; (b) the summer-maize growing season and (c) the annual basis.
Water 10 00789 g007
Figure 8. The spatial distribution of the long-term average reference evapotranspiration (ET0) on (a) the winter-wheat growing season; (b) the summer-maize growing season and (c) the annual basis.
Figure 8. The spatial distribution of the long-term average reference evapotranspiration (ET0) on (a) the winter-wheat growing season; (b) the summer-maize growing season and (c) the annual basis.
Water 10 00789 g008
Figure 9. The spatial distribution and statistical significance of the trends in evapotranspiration (ET0) on (a) the winter-wheat growing season, (b) the summer-maize growing season, and (c) the annual basis.
Figure 9. The spatial distribution and statistical significance of the trends in evapotranspiration (ET0) on (a) the winter-wheat growing season, (b) the summer-maize growing season, and (c) the annual basis.
Water 10 00789 g009
Figure 10. The spatial distribution of the long-term average aridity index on (a) the winter-wheat growing season; (b) the summer-maize growing season and (c) the annual basis.
Figure 10. The spatial distribution of the long-term average aridity index on (a) the winter-wheat growing season; (b) the summer-maize growing season and (c) the annual basis.
Water 10 00789 g010
Figure 11. The spatial distribution and statistical significance of the trends in the aridity index on (a) the winter-wheat growing season; (b) the summer-maize growing season and (c) the annual basis.
Figure 11. The spatial distribution and statistical significance of the trends in the aridity index on (a) the winter-wheat growing season; (b) the summer-maize growing season and (c) the annual basis.
Water 10 00789 g011
Table 1. The province-wide, funnel area and regional statistics for wheat, maize, and the annual temperature.
Table 1. The province-wide, funnel area and regional statistics for wheat, maize, and the annual temperature.
BeijingTianjinHebeiHenanShandongFunnel AreaRegional
Tmin
Wheat growing season1.32.7−1.24.62.82.41.9
Maize growing season19.320.917.420.319.319.919.0
Annual7.38.85.09.88.38.27.6
Tmax
Wheat growing season12.311.810.815.111.813.512.5
Maize growing season29.428.928.229.927.329.928.6
Annual18.017.516.620.117.018.917.8
Table 2. The province-wide, funnel area, and regional statistics for the wheat, maize, and annual temperature trends.
Table 2. The province-wide, funnel area, and regional statistics for the wheat, maize, and annual temperature trends.
BeijingTianjinHebeiHenanShandongFunnel AreaRegional
Tmin (°C/10a)
Wheat growing season0.68 *0.26 *0.53 *0.37 *0.41 *0.53 *0.44
Maize growing season0.53 *0.14 *0.32 *0.17 *0.21 *0.28 *0.24
Annual0.64 *0.22 *0.46 *0.30 *0.34 *0.42 *0.37
Tmax (°C/10a)
Wheat growing season0.25 *0.36 *0.23 *0.23 *0.24 *0.20 *0.23
Maize growing season0.19 *0.23 *0.12 *−0.080.080.020.04
Annual0.25 *0.32 *0.17 *0.10 *0.18 *0.11 *0.15
* Indicates that the temporal trend is significant (p < 0.05).
Table 3. The province-wide, funnel area, and regional statistics for wheat, maize, and the annual average daily solar radiation.
Table 3. The province-wide, funnel area, and regional statistics for wheat, maize, and the annual average daily solar radiation.
R (MJ/day/m2)BeijingTianjinHebeiHenanShandongFunnel AreaRegional
Wheat growing season13.613.413.612.914.013.413.5
Maize growing season18.018.218.217.118.117.717.8
Annual15.515.515.614.715.815.315.4
Table 4. The province-wide, funnel area, and regional statistics for the wheat, maize, and annual daily solar radiation average trends.
Table 4. The province-wide, funnel area, and regional statistics for the wheat, maize, and annual daily solar radiation average trends.
R (MJ/day/m2/10a)BeijingTianjinHebeiHenanShandongFunnel AreaRegional
Wheat growing season−0.18 *−0.28 *−0.19 *−0.20 *−0.20 *−0.25 *−0.20 *
Maize growing season−0.54 *−0.65 *−0.42 *−0.58 *−0.4 *−0.56 *−0.48 *
Annual−0.30 *−0.43 *−0.26 *−0.33 *−0.26 *−0.35 *−0.28 *
* Indicates that the temporal trend is significant (p < 0.05).
Table 5. The province-wide, funnel area, and regional statistics for the wheat, maize, and annual average precipitations.
Table 5. The province-wide, funnel area, and regional statistics for the wheat, maize, and annual average precipitations.
P (mm)BeijingTianjinHebeiHenanShandongFunnel AreaRegional
Wheat growing season99.7109.8108.3260.8185.6138.7176.4
Maize growing season445.3441.4412.5470.7507.4431.7458.0
Annual545.0551.2520.7730.8688.1570.4632.9
Table 6. The province-wide, funnel area, and regional statistics for the wheat, maize, and annual precipitation average trend.
Table 6. The province-wide, funnel area, and regional statistics for the wheat, maize, and annual precipitation average trend.
P (mm/10a)BeijingTianjinHebeiHenanShandongFunnel AreaRegional
Wheat growing season12.16 *8.46 *6.84 *4.016.696.94 *6.85 *
Maize growing season−19.38−20.53 *−17.32 *−2.43−17.65−14.72 *−9.63 *
Annual−4.95−12.59−6.936.05−9.68−5.69−2.98
* Indicates that the temporal trend is significant (p < 0.05).
Table 7. The province-wide, funnel area, and regional statistics for the wheat, maize, and annual average reference evapotranspiration.
Table 7. The province-wide, funnel area, and regional statistics for the wheat, maize, and annual average reference evapotranspiration.
ET0 (mm)BeijingTianjinHebeiHenanShandongFunnel AreaRegional
Wheat growing season571.2569.3508.6531.9552.8545.9531.2
Maize growing season515.7538.3500.2497.4487.0513.8497.4
Annual1086.91107.61008.81028.11032.01059.81026.1
Table 8. The province-wide, funnel area, and regional statistics for the wheat, maize, and annual average evapotranspiration trend.
Table 8. The province-wide, funnel area, and regional statistics for the wheat, maize, and annual average evapotranspiration trend.
ET0 (mm/10a)BeijingTianjinHebeiHenanShandongFunnel AreaRegional
Wheat growing season4.542.96−4.06 *−0.85−3.07−5.24−2.44
Maize growing season−0.01−9.58 *−6.32 *−15.71 *−7.12 *−11.35 *−8.85 *
Annual3.53−7.43−10.73 *−17.44 *−10.58 *−17.49 *−12.04 *
* Indicates that the temporal trend is significant (p < 0.05).
Table 9. The province-wide, funnel area, and regional statistics for the wheat, maize, and annual average aridity indexes.
Table 9. The province-wide, funnel area, and regional statistics for the wheat, maize, and annual average aridity indexes.
AIBeijingTianjinHebeiHenanShandongFunnel AreaRegional
Wheat growing season0.180.200.220.510.350.260.34
Maize growing season0.880.830.840.971.080.860.95
Annual0.510.500.520.730.680.550.63
Table 10. Province-wide, funnel area and regional statistics for wheat, maize, and annual average aridity index trend.
Table 10. Province-wide, funnel area and regional statistics for wheat, maize, and annual average aridity index trend.
AI (/10a)BeijingTianjinHebeiHenanShandongFunnel AreaRegional
Wheat growing season0.019 *0.013 *0.015 *0.0080.0130.016 *0.013 *
Maize growing season−0.04−0.025−0.0240.024−0.027−0.008−0.004
Annual−0.009−0.0080.0000.017−0.0050.0040.005
* Indicates that the temporal trend is significant (p < 0.05).

Share and Cite

MDPI and ACS Style

Ti, J.; Yang, Y.; Yin, X.; Liang, J.; Pu, L.; Jiang, Y.; Wen, X.; Chen, F. Spatio-Temporal Analysis of Meteorological Elements in the North China District of China during 1960–2015. Water 2018, 10, 789. https://doi.org/10.3390/w10060789

AMA Style

Ti J, Yang Y, Yin X, Liang J, Pu L, Jiang Y, Wen X, Chen F. Spatio-Temporal Analysis of Meteorological Elements in the North China District of China during 1960–2015. Water. 2018; 10(6):789. https://doi.org/10.3390/w10060789

Chicago/Turabian Style

Ti, Jinsong, Yuhao Yang, Xiaogang Yin, Jing Liang, Liangliang Pu, Yulin Jiang, Xinya Wen, and Fu Chen. 2018. "Spatio-Temporal Analysis of Meteorological Elements in the North China District of China during 1960–2015" Water 10, no. 6: 789. https://doi.org/10.3390/w10060789

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