Abstract

The effects of agricultural development on observed changes in near-surface atmospheric temperature and moisture from 1960 to 2014 over Northeast China are evaluated using data from 109 meteorological stations. Cultivated land fraction (CF) within a 3 km radius of the meteorological station is used as a quantitative indicator of agricultural intensity. Stations with large CFs experience a less significant increase in air temperature, especially in daily maximum temperature () and a more rapid increase in vapor pressure () and relative humidity (RH) than stations with small CFs, especially during the main growing season (from May to September). Compared with the reference station group with , cooling effects during May to September in terms of daily mean, maximum, and minimum temperature by −0.067°C, −0.081°C, and −0.069°C per decade and wetting effects of May to September regarding by 0.075 hPa and RH by 0.56% per decade exist for the station group with . The cooling and wetting effects can be attributed to the agricultural development and thus should be considered when analyzing the near-surface atmospheric temperature and moisture records in Northeast China.

1. Introduction

Agricultural activities, including the conversion of rainfed to irrigated land and the enhancement of irrigation or fertilization, may increase vegetation activities [1] and alter the energy partitioning between sensible heat and latent heat [2]. Agricultural development-induced increase in evapotranspiration would lower the near-surface air temperature [1, 35] and lead to additional water vapor into the near-surface atmosphere [6, 7]. The effects of agricultural development on climate have been investigated over almost all the main agricultural areas of the world [5, 710], especially about the irrigation effects on the land surface fluxes [1114], temperature [15, 16], atmospheric moisture content [17], cloud [18, 19] and precipitation [20, 21], and so forth.

Previous studies, for both simulations and observations, have shown that the agricultural development, especially irrigation expansion in cropland regions, exerts significant cooling and wetting effects on near-surface atmosphere regardless of their climate regime [5, 2226]. Observational studies typically rely on pairwise comparisons of the changes in temperature, water vapor pressure, relative humidity, or dewpoint temperature [5, 7, 17] between agricultural and nonagricultural locations. As an example, a study in California has revealed a highly significant effect of irrigation on summer average maximum temperature, with substantial 5.0°C cooling for 100% irrigation cover [5, 25]. Comparatively, the agricultural development-induced cooling effect is more obvious on the maximum temperature [17, 27]. As discussed by Gaffen and Ross [28], the considered evapotranspiration from irrigation may be one of the possible causes for the increase in observed near-surface humidity over the United States. The cooling and wetting effects always accompany agricultural development, especially in boreal summer with crop growth. A study in the Great Plains of North America has revealed that irrigated agriculture development has modified the near-surface atmospheric temperature and moisture records in irrigated locations and resulted in a decrease in the mean maximum temperatures during the growing season [29] and an increase in near-surface atmospheric moisture content [17, 30]. Observational studies also have indicated that irrigation had cooling effects on local temperature [7] and wetting effects on atmospheric moisture content [31] in the growing season over Xinjiang, Northwest of China.

Benefiting from the expanded irrigation, increased fertilizers using, and improved crop management, China experienced remarkable agricultural development since the 1960s. The irrigated area in China increased from  ha in 1961 to  ha in 2013, with more rapidly expansion in the northern region. Meanwhile, the effects of irrigation on climate have been investigated over Northern China: from the arid Northwest region [7, 31] and the Yellow River Basin [32, 33] to the Huang-Huai-Hai Plain [34, 35]. As a main agricultural area in China, Northeast China has undergone rapid cropland extension in the last century, where the cropland area increased from 9.8% of Northeast China area in 1908 to 29.6% in 2000 [36]. Studies in this area have paid significant attention to the effect of climate change on food production [37] and agriculture expansion [38]. The changes of near-surface temperature and moisture would contribute to changes of evaporation demand [39], which require comprehensive understanding. In the past 50 years, significant warming trends [40] had been accompanied by land cover change [38, 41], which had created favorable conditions for agriculture development. However, a significant warming and a slight decreasing in precipitation had also been accompanied in this region [42], which might result in increasing irrigation water requirement. A large irrigation project aiming to increase food production is conducting in Northeast China since 2012, in which the irrigated area will be expanded by  ha after the project. Observational evidence of cooling effects with irrigation expansion was found according to the comparison of two rainfed and two irrigated sites in this region [43].

Although the impressions that agricultural development would affect observed air temperature and atmospheric moisture have been acknowledged, the extent of their contributions on temperature and moisture records over Northeast China with a boreal climate remains elusive. The objective of this study is to quantitatively investigate the effects of agricultural development on observed temperature and atmospheric moisture changes in Northeast China.

2. Study Area, Data, and Methods

2.1. Study Area

Northeast China (38°42′N to 53°36′N, 115°24′E to 135°12′E) covers Heilongjiang, Jilin, and Liaoning provinces, as well as the eastern part of the Inner Mongolia Autonomous Region (Figure 1), with an area of  km2 (12.8% of China in its entirety), of which nearly 30% is arable [36, 41]. The topography is dominated by a few mountain ranges, such as the Greater Khingan and the Lesser Khingan. The low-lying terraces formed by the Heilong River and the Ussuri River are the most food productive land, where most of the farmland are transferred from wetlands [38]. The mean annual precipitation decreases from southeastern to northwestern [44], ranging from 1000 mm in the southeastern Changbai Mountains to 300 mm in the Inner Mongolia Plateau [45]. Most precipitation occurred during summer (65.7% of annual precipitation). The main crops in this region are rice, soybean, spring maize, and wheat. The dry spring results in water stress in the initial growing season, and supplemental irrigation is essential to ensure the grain product [46, 47]. At present, the gross irrigation amount per unit area is  m3/ha. Surface water is the main source for irrigation. However, the increasing groundwater withdrawal for irrigation has resulted in serious groundwater crisis in regions with primarily groundwater-fed irrigation.

As one of the important food bases in China, the study area plays an extremely important and irreplaceable role in the national food security and had presented a rapid agricultural development from 1960 to 2014. Taking the Heilongjiang, Jilin, and Liaoning provinces into account, the total grain production also exhibited a rapid rise under agricultural development, from  tons in 1960 to  tons in 2014 (Figure 2). The agricultural cultivated area increased by  ha from 1960 to 2014. The area equipped for irrigation had also experienced a significant increase from  ha in 1960 to  ha in 2014, especially a sharp increase since the late 1990s.

2.2. Data and Methods

The data used in this study includes the daily average temperature (), maximum temperature (), minimum temperature (), diurnal temperature range (DTR) (), and relative humidity (RH) from 1960 to 2014 of stations across Northeast China. The water vapor pressure (, hPa) was calculated according to , , and RH, as referred to by Allen et al. [48], as follows:The meteorological data were provided and quality-tested by the National Meteorological Information Center of the China Meteorological Administration http://data.cma.cn/. The stations with time series of less than 30 years were excluded. After the data were filtered for missing observations, a total of 109 stations remained (Figure 1).

Land use data for 2000, in which the area proportions of cultivated land, forest, residential, industrial, and traffic lands were provided at grid intervals of 1 km, were obtained from the Data Center for Resources and Environmental Sciences of the Chinese Academy of Sciences [49]. The area fractions of the cultivated lands in 2000 are shown in Figure 1. The 109 stations are located in landscapes with markedly different cultivated land uses. The cultivated land fractions (CFs) within radii of 1 km to 10 km, 15 km, and 30 km, centered at each station, were used to evaluate the intensity of cultivated land uses. According to the CFs within a 3 km radius of each station, the 109 stations (with CFs ranging from 0% to 91.4%) could be classified into three groups in this study; namely, 23 were characterized by ( for 9 of these stations), 37 by , and 49 by . The stations with (average ) are mainly located in the north and west mountainous regions with an insignificant agricultural influence. The stations with intensive cultivated land uses (, average ) are mainly located in the central region with abundant water supply for irrigation. The mean annual precipitation for stations with is only 539.3 mm (Table 1) (453.9 mm during the growing season). Crop water requirement cannot be supplied by precipitation during the growth period. Taking the paddy rice as an example, the water requirement ranges from 250 mm to 750 mm, with irrigation requirement from 80 mm to 450 mm [50]. Therefore, supplemental irrigation is needed to ensure grain production [46, 47].

Both parametric and nonparametric methods (e.g., the linear regression method and the Mann–Kendall test) are widely used to identify trends in the data. The nonparametric tests are more suitable for identifying the trends in the field of hydrology and meteorological study, where the nonnormally distributed, censored data, including missing values, are frequently encountered in these time series [51]. By contrast, the time series are required to be independent and normally distributed using parametric methods. The nonparametric Mann–Kendall test with a trend-free, prewhitening method [52] was widely used to identify temperature and atmosphere moisture trends [31, 43, 53]. Therefore it is used in this study.

The dependence of the trends in observed surface air temperature and atmospheric moisture content on CFs was evaluated on the basis of regression slopes and linear correlation coefficients according to the linear regressions of temperature and atmospheric moisture content trends on the CFs, following the method of Han et al. [7, 54]. The statistical significance of the correlation coefficient and linear regression was evaluated through a -test. The CFs within radii of 1 km to 10 km, 15 km, and 30 km, centered at each station, were all tested. A certain radius generating the most significant correlations would be regarded as optimal and be selected as the index representing the intensity of agricultural development.

3. Results

3.1. Observed Temperature and Atmospheric Moisture Trends

The trends of yearly mean temperature (, , , and DTR) and atmospheric moisture (ea, RH) in the growing season (from May to September), nongrowing season (from October to April), and all 12 months of the 109 stations from 1960 to 2014 are evaluated. The average annual trend slopes of , , , DTR, ea, and RH are 0.30°C·decade−1, 0.20°C·decade−1, 0.43°C·decade−1, −0.23°C·decade−1, 0.08 hPa·decade−1, and −0.33%·decade−1, respectively (Table 2). Owing to the more significant increase in , DTR presents a decreasing trend, especially in the middle and south regions with large CFs. The air temperature in the growing season presents a less obvious increase than that in the nongrowing season, whereas the vapor pressure in the growing season presents a more significant increase than that in the nongrowing season (Table 2). For all stations, 106, 90, 108, and 82 are significant at 95% confidence level for the growing season in terms of , , , and DTR, respectively. The trends in ea and RH in the growing season are significant at 95% confidence level for 66 and 54 of all stations, respectively.

The temperature trends during the growing season of stations in the north, mainly with mountainous terrain and small CF, are generally higher than those in middle and south regions with large CFs, especially for (Figure 3). On the contrary, ea exhibits a more obvious increase, and RH indicates a less obvious decrease in the middle and southern regions. The near-surface air temperature and atmospheric moisture trends of the 109 stations are significantly correlated with the CFs within radii of 1 km to 10 km, 15 km, and 30 km, and the correlations are significant during the growing season (Table 3). Stations with large CFs are likely to experience less significant near-surface atmospheric warming and more significant wetting than stations with small CFs. Although the correlations are essentially the same for different radii, CFs within a 3 km radius of each station are used to denote the extent of agricultural land use thereafter because they are the most significant.

The trends of , , , and DTR decreased with increasing CFs, especially during the growing season (Table 4). The correlation between trends and CFs is the most significant (, ), and and DTR trends are only significant in the growing season with 90% confidence level. In October to April, no significant correlation is found, except for a weak correlation in trends (, ). The trends of , , , and DTR in the growing season are expected to decrease by 0.10°C·decade−1, 0.14°C·decade−1, 0.10°C·decade−1, and 0.08°C·decade−1, respectively, with full-covered cultivated land use (i.e., 100% CF) (Figure 4). The ea and RH trends of all the 109 stations have high positive correlations with the CFs (at 99% confidence level, except the nongrowing season ea, which is at 95% confidence level) (Table 4). With the full-covered cultivated land use (i.e., 100% CF), the ea trends in annual mean, growing season, and nongrowing season are expected to increase by 0.07, 0.13, and 0.02 hPa·decade−1, respectively. Meanwhile, the trends in annual mean, growing season, and nongrowing season RH are expected to increase by 0.72, 0.98, and 0.52%·decade−1, respectively.

3.2. Cooling Effects on Stations with High Levels of Cultivated Land Use

, , and . The 49 stations with are taken as with significant influences of agricultural development. As shown in Figure 4, the plots of temperature and moisture trends against CFs keep relatively stable after 0.5, which was usually taken as a break value in other studies [5, 7]. In order to isolate the impact of agricultural development from other forcing, a reference station group containing enough stations with minimal human activities is needed. Limited by the fact that the stations with minimal influences of agricultural development are less (only 10 stations are characterized with ), the station-averaged trend slopes of the station group with are close to those of stations with , and they are widely distributed in the study area. Besides, no significant statistical correlation exists between temperature trends and CFs as well as the urban land fractions. Therefore, the 23 stations with can provide a reference of background climate change, and the cooling and wetting effects of agricultural development can be detected through the comparison between the two groups in spite of their different number. The average values of May to September in terms of , , and trends of the station group with are 0.225°C·decade−1, 0.152°C·decade−1, and 0.336°C·decade−1, respectively, which are considerably smaller than those of the reference station group (0.301°C·decade−1, 0.253°C·decade−1, and 0.410°C·decade−1, resp.) (Table 5). Influenced by the more significant cooling effect on during the daytime because of enhanced evapotranspiration by agricultural development, the DTR trend of the station group with shows more significant decreasing (−0.194°C·decade−1) than that with (−0.145°C·decade−1). However, the differences during the period from October to April are fairly smaller than those in the growing season.

The variations in the station-averaged growing season in terms of , , , and DTR from 1960 to 2014 of the two groups are compared according to their anomalies from the means of 1960 to 1965 (Figure 5). The two station groups are characterized by similar variations; namely, a slight increase in , , and occurred from 1960 to the early 1990s, followed by a more obvious increase in 2000; and the DTR showed a slight decrease in the early 1980s. However, the station group with exhibited a less significant increase in , , and and a more significant decrease in DTR than the station group with , which can be easily detected from their differences. The difference has begun to increase rapidly since the late 1990s, when Northeast China began to experience a rapid agricultural development. Similar to the urban-minus-rural method used in detecting urban heat island effects [7], the cooling effects on the station group with are determined by evaluating the differences in temperature trends with the reference station group with CF < 0.2. , , , and DTR in the growing season from 1960 to 2014 of the station group with had slowly declined by approximately −0.067°C·decade−1, −0.081°C·decade−1, −0.069°C·decade−1, and −0.052°C·decade−1 with respect to the reference station group (Table 5).

3.3. Wetting Effects on Stations with High Levels of Cultivated Land Use

The atmospheric moisture for the stations with is apparently higher than that for the reference group . The trends in the growing season in terms of ea and RH of the station group with are 0.141 hPa·decade−1 and −0.196%·decade−1, respectively, whereas those of the reference group are only 0.077 hPa·decade−1 and −0.801%·decade−1, respectively (Table 5). The differences in the station-averaged data from October to April in terms of ea and RH between the two groups (0.011 hPa·decade−1 and 0.341%·decade−1) are fairly smaller than those in the growing season (0.064 hPa·decade−1 and 0.605%·decade−1).

The anomalies in the station-averaged near-surface ea and RH during the growing season of the two groups are compared in Figure 5. The averaged ea increases while RH decreases during 1960–2014, which ought to be caused by the significant warming over the study area. According to (1), the significant warming has enlarged the saturation vapor pressure, as well as ea, although the increasing amount is smaller than the saturation vapor pressure with decreasing RH. The averaged ea of stations with showed a minor upward trend overall, with a large decrease from 1960 to 1969, followed by a slight increase from 1970s to 2000, and then a small decrease until 2009. Compared with the reference group, the averaged ea of stations with presented an apparent increase after the 1980s, reflecting a more significant atmospheric wetting. The near-surface RH of the two groups had shown a slight decrease from 1960 to 2014. However, the averaged RH of stations with exhibited a less significant decrease than that of stations with , the difference of which had enlarged since the late 1990s. According to the difference in the trends in station-averaged data regarding ea and RH, the wetting effects on the station group with resulted in an increase on ea with 0.075 hPa·decade−1 and an increase on RH with 0.560%·decade−1 from 1960 to 2014 in Northeast China (Table 5).

4. Discussion

In Northeast China, the temperature trends during the growing season are significantly negatively correlated with the CFs of meteorological stations, whereas the near-surface atmospheric content trends are significantly positively correlated with the CFs. Nevertheless, most of the cultivated land distributes on the middle region with low elevation, making the CFs of stations negatively correlate with altitude and leaving a possibility that the dependence of temperature trends and atmospheric content trends on the CFs comes from the geographic locations of stations. No significant correlations exist between the CFs and the longitude of stations (the correlation coefficient is 0.13), and the correlation between the CFs and the latitude of stations is weaker than that between the CFs and the temperature and moisture content trends. Although the CFs of stations are significantly correlated with altitude (correlation coefficient is −0.51), the trends in temperature and atmospheric moisture content at stations with for all 78 stations with elevation between 90 and 700 m (39 stations, with the average growing season trend slope of 0.241, 0.163, and 0.353°C·decade−1 for , , and , resp., as well as 0.129 hPa·decade−1 and −0.291%·decade−1 for ea and RH) remain more significant than those at stations with (13 stations, with the average growing season trend slope of 0.305, 0.248, and 0.435°C·decade−1 for , , and , resp., as well as 0.104 hPa·decade−1 and −0.729%·decade−1 for ea and RH). However, the difference in elevation (average values are 247 and 340 m) is negligible. Besides, the 23 reference stations with are widely distributed in the study area (116°7′E to 129°36′E, 38°54′N to 50°29′N, and an altitude from 91.5 to 1003 m), and no correlation exists between the trends and the locations of these stations. These findings suggest that the geographic locations of stations have a minimal effect or no effect on the relationship among the trends in temperature, ea, RH, and CFs during the growing season. Other topographic factors, such as the terrain slope and aspect, were not considered in the analysis as most stations in the study area are located in the plain.

From 1960 to 2014, Northeast China had experienced a rapid agriculture development (Figure 2). Accordingly, the land uses surrounding the stations with large CFs would be significantly affected in the growing season, which can be detected from the normalized difference vegetation index (NDVI) of the pixel where the stations are located; it is calculated from the biweekly NDVI data of 1982 to 2011 obtained from the third-generation Global Inventory Monitoring and Modeling Studies-NDVI3g database [55], with a spatial resolution of 8 km. The differences in the average yearly mean NDVI at the pixel where the stations are located between the station groups and from 1982 to 2011 are compared in Figure 6. The land surrounding the stations with had experienced relatively obvious increasing vegetation activities in the growing season compared with the stations with . By contrast, a weak declining trend in the difference in the average yearly mean NDVI in the nongrowing season existed. The difference in NDVI in the growing season is highly negatively correlated with the differences in , , and (the correlation coefficients are −0.50, −0.40, and −0.39, resp., with ) but is positively correlated with the difference in and RH (the correlation coefficients are 0.18 and 0.43, resp., with ). Meanwhile, no significant correlation exists between precipitation trends and CFs, and no obvious difference is identified in the variation of annual precipitation between the two groups. Therefore, the cooling and wetting effects on the stations with in Northeast China would be potentially attributed to the agricultural development of the areas surrounding these stations. However, the mechanism of agricultural development-induced cooling and wetting effects is usually complicated; whether it is dominated by direct evaporation cooling or indirect effects, such as increase in cloud cover in Northeast China, should be further detected.

According to pairwise comparisons, the agricultural development from 1960 to 2014 is associated with cooling effect on from May to September for stations with by −0.081°C·decade−1, the magnitude of which is 35.4% of the trend slope of the reference station group (0.229°C·decade−1), indicating that the warming for the station group with has been reduced by 35.4%, compared with the warming of the reference group when regarding it as natural climate variability. The cooling effects on and during the growing season (−0.067°C and −0.069°C·decade−1, the magnitude of which are 22.8% and 17.4% of the trend slope of the reference station group, resp.) are smaller than that of .

Accompanied with cooling, agricultural development is associated with the wetting of from May to September and RH by approximately 0.064 hPa·decade−1 and 0.605%·decade−1, respectively. Our result on cooling effects agrees with the results reported by Zhu et al. [43], in which the difference in the trends of from 1956 to 2008 between two highly and lightly irrigated sites in the Jilin province is 0.085°C·decade−1. However, the cooling effects are lower than those in Northwest and North Plain of China with intensive irrigation. In Xinjiang, Northwest China, the irrigation expansion from 1959 to 2006 is associated with the cooling of , , and from May to September by −0.148, −0.119, and −0.138°C·decade−1, respectively, in the station group with extensive irrigation [7]. In the Huang-Huai-Hai Plain of China, irrigation indicated a cooling effect of 0.12°C·decade−1 on summertime daily maximum temperature [34]. In California, USA, the irrigation expansion from June to August leads to the fact that the current irrigation level exceeds 50%, and the summer decreased by −0.14°C to −0.25°C·decade−1 from 1915 to 1979 owing to irrigation [25]. Irrigation expansion is an important aspect of agricultural development that is at the origin of cooling and wetting effects [5]. The lower irrigation intensity in Northeast China (the irrigated area is only 35.9% of the total sown area in 2014) may be one of the reasons for its lower cooling effect than above regions. However, the present irrigation project with water-saving technologies in Northeast China, by which the newly increase of irrigated area is approximately 30% of that in 2013, would result in further cooling and wetting effects in these agricultural regions, which could have positive effects on agriculture. However, given that the irrigation amount per unit area with water-saving technologies is less than that with flood irrigation, the cooling and wetting effects of irrigation will therefore likely be different in the future than in the past century. However, further studies are needed.

Specific uncertainties associated with this study should also be noted. Only the land use data for year 2000 are used in this study because of data shortage, which may lead to uncertainties illustrated in Figure 4. According to Lobell and Bonfils [25], the spatial distributions of cultivated lands, types of crops grown, or irrigation methods used could cause differences among stations with similar CFs. Urban land uses may also affect the atmospheric moisture, clouds, and precipitation. These factors are not considered in the present study, although stations with an urban land fraction larger than 25% within a 3 km radius are excluded to eliminate the influence of urban land uses. More observations and temporal information are needed for further understanding of the cooling and wetting effects of agricultural irrigation in Northeast China. Furthermore, we should pay attention to how land cover dynamics, types of crops growth, and especially irrigation interact with climate, as well as the magnitude of this interaction.

5. Conclusions

The preceding analysis demonstrates that observations of near-surface atmospheric temperature and moisture in Northeast China were disturbed by agricultural development. Stations surrounded with large fractions of cultivated land experienced less rapid increases in air temperature, more rapid increases in water vapor pressure, and weak decreases in RH from 1960 to 2014. Compared with the average , , and , as well as and RH, of 23 stations with , those values of 49 stations with exhibited cooling of −0.067, −0.081, and −0.069°C·decade−1 and wetting of 0.075 hPa·decade−1 and 0.56%·decade−1. The difference was associated with agricultural development. In light of this scenario, the agricultural development introduces excessive cooling and wetting signals to the observations in Northeast China. The agricultural development around meteorological stations should be seriously considered for a better understanding of observed near-surface atmospheric temperature and moisture. However, special uncertainties, such as the spatial distributions of cultivated lands, types of crops grown, or irrigation methods, still need to be further noted.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This research was partially sponsored by the National Key Technology R&D Program (2014BAD12B05), the National Natural Science Foundation of China (no. 51579249), the Research Fund (no. 2016ZY06) and Open Research Fund (no. IWHR-SKL-201612) of State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin, China Institute of Water Resources and Hydropower Research, and the National Nonprofit Institute Research Grant from CAFINT (Grant no. CAFINT2015K06).