Next Article in Journal
Wood Vinegar Impact on the Growth and Low-Temperature Tolerance of Rapeseed Seedlings
Next Article in Special Issue
Factors Driving Autumn Quercus Flowering in a Thermo-Mediterranean Area
Previous Article in Journal
Research on Dynamic Monitoring of Grain Filling Process of Winter Wheat from Time-Series Planet Imageries
Previous Article in Special Issue
A Novel Crop Shortlisting Method for Sustainable Agricultural Diversification across Italy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Climate-Smart Tillage Practices with Straw Return to Sustain Crop Productivity

1
State Key Laboratory of Crop Biology, Key Laboratory of Crop Water Physiology and Drought-Tolerance Germplasm Improvement, Ministry of Agriculture and Rural Affairs, College of Agronomy, Shandong Agricultural University, Tai’an 271018, China
2
Shandong Key Laboratory of Eco-Environmental Science for Yellow River Delta, Binzhou University, Binzhou 256603, China
*
Author to whom correspondence should be addressed.
Agronomy 2022, 12(10), 2452; https://doi.org/10.3390/agronomy12102452
Submission received: 1 September 2022 / Revised: 27 September 2022 / Accepted: 30 September 2022 / Published: 10 October 2022
(This article belongs to the Special Issue Adaptations to Climate Change in Agricultural Systems)

Abstract

:
Climate change seriously threatens global crop production. However, there are few reports on field crop yield and yield components based on long-term different climate conditions. The objectives of this study were to identify and compare the differences in crop yield and yield components in long-term tillage and straw returning under different climate regions. Conventional tillage (CT) and rotary tillage (RT) in combination with no straw return and whole straw return (S) were conducted under a wheat (Triticum aestivum L.)–maize (Zea mays L.) cropping system in cool-wet and warm-dry regions from 2010 to 2019. We hypothesized that long-term suitable tillage under warm-dry or cool-wet regions can increase the yield and components of wheat and maize, and temperature and precipitation had significant effects on crop yield and yield components. Conventional tillage with straw return (CTS) in the warm-dry region and rotary tillage with straw return (RTS) in the cool-wet region can increase the yield and yield components of wheat and maize, respectively, compared with CT. The yield stability of wheat was higher than that of maize under the two climate conditions. Compared with tillage practices, the effects of experimental sites and straw return on crop yield and yield components were more remarkable. The combination of mean temperature, annual precipitation, and yield components explained 75% and 100% of the variance in the wheat yield and maize yield, respectively. The thousand-kernel weight was the key factor in regulating wheat yield, and kernel number was the key factor in regulating maize yield. In conclusion, the combination of rotary tillage in cool-wet regions or conventional tillage in warm-dry regions with straw return is a good technique for increasing crop security.

1. Introduction

In the past 100 years, the temperature of the Earth’s surface has increased by 0.74 °C [1], and it will increase by 1.0–1.5 °C before 2100 [2]. Due to global warming, the number of regions with relatively warm and extremely dry weather may increase by the end of this century [3]. Climate change is very unstable and difficult to control by people through farmland management measures [4]. Climate change is altering water supplies worldwide and further threatening food productivity [5]. In the coming decades or even centuries, climate change will have a significant impact on the agroecological food system [6]. Due to climate change and global population growth, achieving food security on limited arable land is a major challenge of the 21st century [7]. Climatic vulnerabilities, especially heat waves, drought, or rainstorms occurring in key growth stages, such as anthesis and filling, will lead to serious yield loss [8]. For instance, relatively high temperatures during grain filling can decrease the grain filling rate [9], shorten the grain filling duration [10], and ultimately lead to a significant grain yield reduction. Globally, there is a negative correlation between climate warming and some crop productivity, such as wheat, maize, and barley [11]. Future agriculture will be challenged by changes in the spatial and temporal distribution of precipitation [12]. Precipitation was identified as the most important determinant of many crop yields [13]. The interannual variation in wheat yield is very large due to the unstable seasonal climate pattern, especially the irregular precipitation distribution and high temperature during the grain filling period [14]. In the past 50 years, temperature and precipitation have increased in Northeast China, with reduced snow and increased spring drought [15]. Wheat production in northern Kazakhstan has decreased due to the combination of warming and drought since 1980 [16]. However, the results are often different or even opposite due to the differences in the duration of the experiment, soil type, management measures, and climatic conditions. As a result, the effectiveness of cropping systems is highly geographically specific, and the impact of yield constraints can vary greatly depending on environmental conditions and their interactions with management practices [17]. Crops need to adapt to changing climate conditions to guarantee food security. Understanding the link between climate change and crop productivity is therefore crucial to assessing the adaptability of agricultural systems to future climate changes [18].
Further improvement in crop yields is needed for future food security [19]. Winter wheat–summer maize rotation is the most typical cropping system on the North China Plain [20]. The main objective of the winter wheat–summer maize cropping system is to produce sufficient crop yields to sustain the population [21]. The precipitation in the wheat growing season accounts for 20–30% of the annual precipitation [22]. Low precipitation and high annual variation are the main reasons for the low and unstable wheat yield in this region [23]. Crop yield is determined by yield components, including ear density, kernel number, and thousand-kernel weight [24]. The genetic process of crop yield is mainly related to kernel number and kernel weight [25,26]. Yield and yield components are determined by genotype, cropping system, agronomic management, climatic conditions, and plant–plant interactions [27,28,29]. Therefore, reasonable climate-smart crop management is very important to sustain crop yield.
Conservation tillage, such as minimum or no-tillage combined with mulching and crop rotation, can reduce soil runoff, evaporation, erosion, and degradation [30]. A study found that minimum tillage combined with straw mulching is effective in reducing soil temperature and increasing soil moisture and crop yield [31]. Moreover, the combination of minimum or no-tillage and straw return can effectively promote soil macrocohesion, change the abuse of conventional tillage systems to consume soil nutrients, and increase soil fertility and crop yield [32]. Straw mulching works worldwide by absorbing the kinetic energy of raindrops to reduce soil compaction and erosion caused by precipitation [33]. Straw mulching and returning can conserve water by reducing surface loss and increasing infiltration, thereby improving the availability of soil water, thus increasing crop yield and water use efficiency [34,35]. Conservation tillage under well-drained soils can increase maize yield, while minimum tillage does not appear to significantly affect yield stability [13].
However, the effects of long-term tillage and straw return under different climate conditions on field crop yield and yield components are still unclear. The objectives of this study were: (1) to identify the differences in wheat and maize yields and yield components in long-term tillage and straw returning; (2) to compare the yield and yield component stability of wheat and maize under tillage and straw returning practices; (3) to clarify the correlation between crop yield, yield components, and mean temperature and precipitation. It was hypothesized that long-term suitable tillage under warm-dry or cool-wet regions can increase the yield and components of wheat and maize, and that temperature and precipitation have a significant effect on crop yield and yield components.

2. Materials and Methods

2.1. Experimental Site

One warm-dry region (Tai’an) and one cool-wet (Longkou) region of Shandong Province in northern China were selected in this study (Figure 1). This former and the latter have a temperate continental monsoon climate, with mean temperatures of 14.4 °C and 13 °C and mean annual precipitation values of 580.6 mm and 650 mm, respectively, during 2010–2019 (Figure 2). The experiment was first established in Tai’an (36°09′ N, 117°09′ E) in 2002. Considering the climatic factors, another experiment with the same treatments was conducted in Longkou (37° 38′ 30″ N, 120° 24′ 9″ E) in 2006. The soils in the two regions were the same and were classified as Cambisols (Food and Agriculture Organization). Only the temperature and precipitation were different. The major initial (2002) properties in the 0–20 cm soil layer in Tai’an were: 40% sand, 44% silt, and 16% clay; the pH was 7.09, and the soil contained 7.2 g kg−1 soil organic carbon, 1.3 g kg−1 total nitrogen, 0.8 mg kg−1 available phosphorus, and 41.3 mg kg−1 exchangeable potassium. The major initial (2006) properties within the 0–20 cm soil depth in Longkou were: 47.0% sand, 37.0% silt, and 15.0% clay; the pH was 6.4, and the soil contained 7.6 g kg−1 soil organic carbon, 115.0 mg kg−1 total nitrogen, 36.7 mg kg−1 available phosphorus, and 91.2 mg kg−1 exchangeable potassium.

2.2. Experimental Design

A typical cropping system of winter wheat–summer maize was chosen at the Tai’an and Longkou experimental sites. Winter wheat was generally sown in mid-October and harvested in early June of the next year. Summer maize was sown in mid-June and harvested in early October. Fertilization application was carried out according to the local high-yield field. A 10-year (2010–2019) long-term experiment was conducted in both regions. Two tillage practices included conventional tillage (CT, 30 cm depth) and rotary tillage (RT, 10 cm depth). Four treatments were designed by two maize and wheat straw return methods, with no straw return or whole maize and wheat straw return (S). Three replicates were set for each treatment, and the size of each plot was 15 × 4 m2 at Tai’an and 15 × 15 m2 at Longkou. The operational procedures for conventional tillage were: maize was harvested mechanically, and straw was crushed (3–5 cm long) and returned to the field at the same time for straw return treatment (or after harvesting the maize manually, cut the straw against the ground and move it out of the farmland for no straw return treatment)–fertilizer application–stubble plowing with notched harrow–plowing with moldboard plow (25 cm)–rotary tillage with rotavator (10 cm)–wheat sowing–irrigation–herbicide application. The operational procedures for rotary tillage were: maize was harvested mechanically, and straw was crushed (3–5 cm long) and returned to the field at the same time for straw return treatment (or after harvesting the maize manually, cut the straw against the ground and move it out of the farmland for no straw return treatment)–fertilizer application–stubble plowing with notched harrow–rotary tillage with rotavator (10 cm)–wheat sowing–irrigation–herbicide application. Rotary cultivator is an agricultural machine that breaks the soil by rotating the blade teeth. Tractors are often used as power sources for rotary cultivators to provide power for driving and operation.

2.3. Determination of Crop Yield and Components

The wheat grains within a 1 m2 area were collected in each treatment with triplicates at maturity. The maize cobs of 5 m double rows were harvested from the center of each treatment with three replicates. The samples were air dried and threshed, and then ear density, kernel number, and thousand-kernel weight were measured to calculate the yield.
Wheat grain yield (t ha−1) = Ear density (ears ha−1) × kernel number per ear
× thousand-kernel weight (g 1000 kernels−1)/109
× (1 − sample moisture content%)/(1–13%)
Maize grain yield (t ha−1) = Ear density (ears ha−1) × kernel number per ear
× thousand-kernel weight (g 1000 kernels−1)/109
× (1 − sample moisture content%)/(1–14%)

2.4. Determination of Equivalent Ratios (Ers)

The relative values of the wheat and maize yield and components were standardized based on those in the CT, which was considered to be 1.

2.5. Determination of the Stability Index of Yield and Yield Components

The stability index (SI) is an important parameter for measuring the sustainability and stability of a system. The formula for calculating the SI is as follows [36]:
SI = (Y − σ)/Ymax
where Y represents the average yield or components per unit area during the experiment (t ha−1), σ is the standard deviation, and Ymax is the largest crop yield or components among all experimental years (t ha−1). The SI can be any value between 0 and 1; the greater the SI, the more stable the yield [37].

2.6. Statistical Analyses

Illustrations were prepared via Sigmaplot 12.5. The significance of differences between the treatments was analyzed using Duncan’s test and independent-samples t tests at p < 0.05 and p < 0.01 via SPSS 18.0 Statistical Analysis System software [38]. Structural equation modeling (SEM) was used to evaluate the direct and indirect relationships between wheat and maize yield and mean temperature, annual precipitation, and yield components. SEM analysis was performed using Amos 21.0 (Amos Development Corporation, Chicago, IL, USA).

3. Results

3.1. Wheat Yield, Yield Components, and Ers Analysis

After the ten-year experiment, the mean ear density and kernel number ranged from 651.8 (RTS)–676.7 (104 ha−1) (CTS) and 39.2 (RT)–41.4 (CT), respectively, at the warm-dry site (Figure 3A). The mean thousand-kernel weight and wheat yield of RTS increased by 10.2% and 6.8%, respectively; RT decreased the mean thousand-kernel weight and wheat yield by 1.3% and 8.5%, respectively, compared with CT (Table 1). The wheat yield and yield components, except for the thousand-kernel weight of rotary tillage, were lower than those of conventional tillage. Straw return increased the wheat yield and yield components compared with no straw application.
At the cool-wet site, compared with CT, the mean kernel number of CTS increased by 2.4%; the mean thousand-kernel weight of RT increased by 2.1%; and the ear density and wheat yield of RTS increased by 3.4% and 7.3%, respectively (Figure 3B; Table 1). The wheat yield and components except for kernel number under rotary tillage were higher than those under conventional tillage. Straw return increased the wheat yield and components compared with no straw application.
Box and whisker plots of the Ers of wheat yield and components were constructed for tillage and straw management at the warm-dry (Figure 4A) and cool-wet sites (Figure 4B). The median Ers of the ear density in the CTS was the largest at the warm-dry site (Figure 4A). The median Ers of the kernel number, thousand-kernel weight, and wheat yield in the RTS were the largest. The median Ers of the thousand-kernel weight and wheat yield in RT were the lowest. In addition to the Ers value of the kernel number in the CTS, the Ers value of the thousand-kernel weight in the RT, and the Ers value of the wheat yield in the RTS, the box plots for the other tillage and straw management practices were noticeably asymmetric. The comparative shortness of the box plots for the CTS and RT indicated that the Er values of the ear density and wheat yield were distributed similarly. The comparative tallness of the box plots for the RTS indicated that the Er values of the wheat yield and components showed wider distributions.
The median Ers of the ear density, kernel number, thousand-kernel weight, and wheat yield in the RTS were the largest at the cool-wet site (Figure 4B). The median Ers of the thousand-kernel weight and wheat yield in CT and RT were lower. There was no significant difference in these indexes among the four treatments. In addition to the Ers value of the wheat yield in the CTS, the box plots for the other tillage and straw management practices were noticeably asymmetric. The comparative shortness of the box plots for the CTS and RT indicated that the Er values of the ear density were distributed similarly. The comparative tallness of the box plots for the RTS indicated that the Er values of the kernel number showed wider distributions.

3.2. Maize Yield, Yield Components, and Ers Analysis

At the warm-dry site, compared with CT, the mean ear density, kernel number, thousand-kernel weight, and maize yield of CTS increased by 0.8%, 2.9%, 5.7%, and 10.2%, respectively; the mean ear density, kernel number, and maize yield of RT decreased by 4.3%, 0.6%, and 2.8%, respectively (Figure 5A; Table 1). The maize yield and components except for the thousand-kernel weight of rotary tillage were lower than those of conventional tillage. Straw return increased the maize yield and components compared with no straw application. At the cool-wet site, compared with CT, the mean ear density, kernel number, thousand-kernel weight, and maize yield of RTS increased by 1.5%, 5.4%, 3.8%, and 10.9%, respectively; the mean ear density of RT decreased by 1.0% (Figure 5B; Table 1). The maize yield and components of rotary tillage were higher than those of conventional tillage. Straw return increased maize yield and components compared with no straw application.
Box and whisker plots of the Ers of the maize yield and components were constructed for tillage and straw management at the warm-dry (Figure 6A) and cool-wet (Figure 6B) sites. The median Ers of the ear density in the CT, kernel number in the RTS, thousand-kernel weight, and maize yield in the CTS were the largest at the warm-dry (Figure 6A) site. There were no significant differences in the median Ers of the ear density, kernel number, and thousand-kernel weight between these four treatments. The maize yield of CTS was significantly higher than that of CT and RT. The box plots for tillage and straw management were noticeably asymmetric. The comparative shortness of the box plots for the RT and RTS indicated that the Er values of the maize yield were distributed similarly. The comparative tallness of the box plots for the RT and CTS indicated that the Er values of the ear density and maize yield showed wider distributions, respectively.
The median Ers of the ear density, kernel number, thousand-kernel weight, and maize yield in the RTS were the largest at the cool-wet site (Figure 6B). There were no significant differences in the median Ers of the ear density, kernel number, and thousand-kernel weight between these four treatments. The maize yield of RTS was significantly higher than that of CT. In addition to the Ers value of the kernel number in the RT and thousand-kernel weight in the CTS, the box plots for the other tillage and straw management were noticeably asymmetric. The comparative shortness of the box plots for the CTS, RT, and RTS indicated that the Er values of the kernel number were distributed similarly. The comparative tallness of the box plots for the RT and CTS indicated that the Er values of the ear density and maize yield showed wider distributions, respectively.
From 2010 to 2015, there was a significant difference in wheat yield between straw return and no straw return, and there was no significant difference from 2016 to 2019, and there was a significant difference in maize yield only from 2015 to 2017 at the warm-dry site (Figure 7A). Wheat yield at the cool-wet site was significantly different only in 2010, and maize yield was significantly different only in 2014 and 2015 between straw return and no straw return (Figure 7B).
In addition, 2010 and 2014 were the eighth years of the warm-dry and cool-wet trials, respectively, and the comparison was conducted from 2010 to 2015 at the warm-dry site and from 2014 to 2019 at the cool-wet site. At the warm-dry site, in addition to RT, the wheat yield decreased at first (2010–2014), and then increased gradually, and the maize yield also showed the same trend under CT and CTS. At the cool-wet site, the wheat yield increased gradually, but it fluctuated greatly. The maize yield first increased (2014–2017), and then decreased (2018–2019). The changes in crop yield at the two sites were opposite.

3.3. Yield and Yield Component Stability of Wheat and Maize

The ear density and kernel number stability of maize were greater than those of wheat at the warm-dry site (Figure 8A). The thousand-kernel weight stability of wheat and maize was similar. The yield stability of wheat increased by 10.4% compared with that of maize. The ear density stability of wheat under CT was obviously higher than that under RTS, and the other yield components and yields of wheat and maize showed no obvious distinction under tillage and straw management. The ear density stability in wheat and maize and kernel number stability in maize under rotary tillage were lower than those under conventional tillage. The kernel number stability in wheat and yield stability in wheat and maize under rotary tillage were greater than those under conventional tillage.
The ear density and kernel number stability of wheat and maize were similar at the cool-wet site (Figure 8B). The kernel number and yield stability of wheat increased by 4.4% and 4.8%, respectively, compared with those of maize. The thousand-kernel weight stability of wheat was lower than that of maize. The yield components and yield of wheat and maize were not obviously different under tillage and straw management. The ear density and yield stability of maize under rotary tillage were greater than those under conventional tillage. The kernel number stability in wheat and thousand-kernel weight stability in maize under rotary tillage were lower than those under conventional tillage.

3.4. Correlation Analysis between Crop Yield and Components

There was an apparent negative correlation between ear density and wheat yield at the 0.01 level at the warm-dry and cool-wet sites (Table 2). There was a positive correlation between kernel number and wheat yield at the 0.05 level, thousand-kernel weight, and wheat yield at the 0.01 level. There were positive correlations between ear density, kernel number, thousand-kernel weight, and yield of maize at the 0.01 level.
Under the conditions of long-term tillage practice, straw management, the experiment site, and their interaction had significant effects on the thousand-kernel weight of wheat (Table 3). Site and straw management had significant effects on wheat and maize yields. The site had significant effects on the ear density of wheat and the kernel number of maize. The interaction of site and tillage practice had a significant effect on maize yield. Tillage practice, the interaction of tillage practice and straw management, and the interaction of site, tillage practice, and straw management had no significant effects on crop yield or components.

3.5. The Correlation Analysis between Climate and Crop Yield

At the warm-dry site, there was a positive correlation between the ear density of wheat and the mean temperature in the wheat season at the 0.05 level (Figure 9). The difference between the wheat yield and mean temperature in the wheat season was nearly significant at the 0.05 level. The maize yield was positively correlated with the mean temperature in the maize season at the 0.01 level (Figure 10). At the cool-wet site, there was a positive correlation between the kernel number of wheat and the mean temperature in the wheat season at the 0.05 level (Figure 9). The difference between the wheat yield and mean temperature in the wheat season was nearly significant at the 0.05 level. The kernel number, thousand-kernel weight of maize, and maize yield were positively correlated with the mean temperature in the maize season at the 0.05 level (Figure 10). Other indices had no obvious divergence with mean temperature and precipitation.

3.6. The Relationship between Climate and Crop Yield Trends

Climate and crop yield trends were the values of the following year minus the previous year. At the warm-dry site, the maize yield trend was obviously correlated with the mean temperature trend in the maize season (Figure 11). At the cool-wet site, the maize yield trend was obviously correlated with the mean temperature trend at the 0.01 level in the maize season. The wheat yield trend had no obvious correlation with the mean temperature and precipitation trends, and the maize yield trend had no apparent divergence from the mean precipitation trend.

3.7. Structural Equation Modeling Analysis

Structural equation modeling further revealed the direct and indirect effects of mean temperature, annual precipitation, and yield components on crop yield (Figure 12). The combination of mean temperature, annual precipitation, and yield components explained 75% and 100% of the variance in the wheat yield and maize yield, respectively. The mean temperature and annual precipitation had weak direct effects but stronger indirect effects on crop yield through their associations with yield components. The mean temperature had a significant positive correlation with the ear density and thousand-kernel weight of wheat and with the ear density, kernel number, and thousand-kernel weight of maize. The annual precipitation had no significant correlation with the yield components, which was inconsistent with our hypothesis. The ear density of wheat had a strong negative effect on thousand-kernel weight. The thousand-kernel weight played a major role in regulating wheat yield, and the kernel number played a major role in regulating maize yield.

4. Discussion

4.1. Effects of Tillage Practices and Straw Management on Crop Yield and Components

The ear density and kernel number of wheat distinguish the sink stage from the sowing and anthesis stages. Minimum-tillage or no-tillage conditions were shown to increase the ear density and kernel number [39]. However, there is research that shows the opposite result [40]. In our study, the yields of wheat and maize under rotary tillage at the warm-dry site were lower than those under conventional tillage, mainly because of the higher ear density under conventional tillage, which played a major role. CTS had the best effect on crop yield because the thousand-kernel weight in CTS was higher than that of CT, while RT had the worst effect. At the cool-wet site, the yields of wheat and maize under rotary tillage were higher than those under conventional tillage. Rotary tillage mainly increased the ear density of wheat, and the kernel number played a major role in the maize yield. RTS had the best effect on crop yield and yield components because the kernel number in RTS was higher than that in RT, while CT had the worst effect. However, no significant differences in yield components were observed between RTS and CTS or between RT and CT. In warm-dry and cool-wet conditions, tillage methods showed opposite effects on crop yield and yield components. The yield stability of wheat was higher than that of maize because the stability of the kernel number was better in long-term planting. Moreover, the difference in yield and yield components between treatments in the cool-wet treatment was less than that in the warm-dry treatment. However, there was no significant difference in yield stability between wheat and maize in different treatments. This may be because tillage and straw returning can increase soil organic carbon, which is conducive to the growth of plant roots, increasing the ability to absorb water and nutrients, thus stabilizing and improving crop growth [41]. The improvement of agricultural productivity depends on the local environment and soil conditions [42]. In Europe, the implementation of conservation tillage may indeed reduce yield, which is related to tillage technology, rotation, crop type, and soil texture [43].
Different tillage methods are important for increasing crop yield. Minimum tillage or no-tillage combined with straw return can increase crop yield and reduce environmental pollution and soil degradation [44,45]. Compared with conventional tillage, the wheat and maize yields increased by 12.3–16.9% under minimum tillage combined with straw return [46]. Reduced or no-tillage alone cannot improve crop yield and yield stability [47]. After long-term wheat and maize cropping, rotary tillage alone cannot produce good effects. This is because long-term rotary tillage formed a plow bottom at a depth of 10 cm, the tractor pulls the rotary cultivator to cultivate the land, and the subsequent process of sowing, fertilization, and harvesting makes the soil below 10 cm mechanically compacted, and hindering root elongation and reducing crop yield. Crops under conventional tillage can grow more fine roots [48], thus enhancing yield growth. Straw return increased crop yield and yield components compared with no straw return. A global meta-analysis showed that straw returning can increase soil organic matter by 12.8% [49]. A previous study has also shown that crop yields increase as soil organic matter increases [50]. However, in the long-term experiment, straw return did not continuously show a significantly higher effect on crop yield than no straw return. Therefore, there was no need to return straw to the field every year after several years, and the excess straw can be used for feed, energy, and other purposes [51].
In addition, the production costs and energy consumption of different tillage practices and straw use are also very important to evaluate the yield and components of wheat and maize. In this study, production inputs include mechanical oil consumption such as tilling (28.9 USD ha−1 for CT and 17.5 USD ha1 for RT), sowing (35.4 USD ha1), stubble control (36.6 USD ha1), pesticide spraying (12.8 USD ha1), and harvesting (59.8 USD ha1), and other costs such as seed (86.6 USD ha1), irrigation (0.3 USD ha1), fertilization (571.4 USD ha1), pesticide (103 USD ha1), and labor (107.1 USD Person1 ha1). The results showed that the production cost of rotary tillage was 1.1% lower than that of conventional tillage.

4.2. Effects of Temperature and Precipitation on Crop Yield and Components

The ear density and kernel number of wheat in warm-dry conditions were 16.1% and 0.9% higher than those in cool-wet conditions, and the thousand-kernel weight and yield were 17.9% and 16.4% lower, respectively. The thousand-kernel weight of maize in warm-dry conditions was 1.9% lower than those in cool-wet conditions, and the ear density, kernel number, and yield of maize were 3.1%, 9.5%, and 9.8% higher, respectively. This result indicated that there was a certain regional effect of tillage that was due to the different climatic conditions. In the wheat season, the temperature at the warm-dry site was higher than that at the cool-wet site, and precipitation at the cool-wet site was greater than that at the warm-dry site. In the maize season, cool-wet maize had wider temperature and precipitation range values than warm-dry maize. The temperature and precipitation trend range values at the cool-wet site were greater than those at the warm-dry site.
Wheat yield is greatly influenced by climatic factors, such as temperature and precipitation. In rainy years, the harvest may be good, but there will be serious yield loss because either the crops cannot be harvested or the grains cannot be dried [52]. The wheat yield decreases with higher temperature but increases with more precipitation [53]. In this study, temperature and precipitation had weak positive and negative effects on wheat yield, respectively (Figure 10). From heading to anthesis, if the precipitation is less, the ear density of wheat will be seriously reduced [54]. This study showed that compared with 229.3 mm precipitation, 76 mm precipitation reduced the ear density, kernel number, and thousand-kernel weight and then reduced the wheat yield by 65%. Compared with other years, the lowest precipitation occurred at the anthesis stage of wheat in 2019, which led to a sharp decrease in ear density at the warm-dry site. For summer maize, the precipitation decreased by 14.1%, the ear density, kernel number, and thousand-kernel weight decreased, and the yield decreased by 45.2%. The kernel number was the main reason for the decrease in summer maize yield [46], which is consistent with the results of this study (Figure 10). At the warm-dry site, the highest kernel number and yield of maize were achieved in 2016 due to relatively higher precipitation during the anthesis stage. In 2018 and 2019, however, although the thousand-kernel weight was higher, the ear density and kernel number were the lowest because of heavy precipitation, which resulted in the lowest yield. At the cool-wet site, the low maize yield and yield components in 2011 may be attributed to the low temperature during the filling stage. Thus, temperature was a major limiting factor in cool-wet crop yields.
An assessment of the impact of seasonal climate change on interannual wheat and maize yields in 92 regions of France revealed regional differences in the most important variables affecting crop yields [18]. A temperature increase of 0.5 °C in winter may reduce India’s wheat production by 0.45 t ha−1 [55]. A 1 °C rise in temperature could reduce global wheat production by 6% [56]. The global maize yields decreased by 7.4% for each degree increase in the global mean temperature [57]. In this study, an approximately 38.1% maize trend was indicated by the mean temperature. When the temperature increased by 1 °C, the maize yield increased by 1.5 t ha−1 at the warm-dry site. Approximately 71.7% of maize trends can be indicated by the mean temperature. When the temperature increased by 1 °C, the maize yield increased by 1.8 t ha−1 at the cool-wet site. Therefore, drought-related temperature anomalies led to a significant decline in maize yield. All crops have a temperature threshold, and the yield will decrease with increasing temperature when it exceeds the threshold [58]. A study reported that maize and soybean yields declined by 17% for a 1 °C increase in temperature [59].

5. Conclusions

After ten years (2010–2019) of tillage and straw return to the field, CTS and RTS were more effective in increasing wheat and maize yields than RT and CT at the warm-dry and cool-wet sites. The yield of rotary tillage at the warm-dry site was lower than that in conventional tillage because of the lower ear density of wheat and maize. However, the yield of rotary tillage at the cool-wet site was higher than that in conventional tillage because of the higher ear density of wheat and higher kernel number of maize. The effects of straw return and site on crop yield and yield components were more obvious than tillage practices. Because the stability of the kernel number was better, the yield stability of wheat was higher than that of maize at the two sites. In the long-term experiment, straw return did not continuously show a higher effect on crop yield than no straw return. Therefore, it can be considered that there was no need to return straw to the field year after year.
In terms of long-term results, wheat and maize yields and yield components in both regions were not always significantly correlated with temperature and precipitation. The explanatory degree of the combination of mean temperature, annual precipitation, and yield components on maize yield was 25% higher than that of wheat yield. The mean temperature and annual precipitation had the highest explanatory degree for the thousand-kernel weight of wheat and kernel number of maize. Moreover, the thousand-kernel weight and kernel number also played the most important roles in wheat yield and maize yield, respectively. Approximately 38.1% and 71.7% of maize trends can be indicated by the mean temperature; when the temperature increased by 1 °C, the maize yield increased by 1.5 and 1.8 t ha−1 at the warm-dry and cool-wet sites, respectively.

Author Contributions

Conceptualization, T.N.; methodology, Z.L.; software, Z.L., N.W. and L.W.; validation, N.W., J.L. and L.W.; formal analysis, Z.L., G.L. and T.N.; investigation, Z.L., N.W. and J.L.; resources, G.L. and T.N.; data curation, Z.L.; writing—original draft preparation, Z.L.; writing—review and editing, Z.L.; visualization, Z.L. and L.W.; supervision, T.N.; project administration, Z.L. and T.N.; funding acquisition, Z.L. and T.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 32101853), the Natural Science Foundation of Shandong Province (Grant No. ZR2021QC189), the Major Science and Technology Innovation Projects of Shandong Province (Grant No. 2019YQ014 and 2021CXGC010804-05), and the Funds of Shandong “Double Tops” Program.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

CT: conventional tillage without straw returning; CTS, conventional tillage with straw returning; RT, rotary tillage without straw returning; RTS, rotary tillage with straw returning; SYI, the yield stability index; L, experimental site; T, tillage practice; S, straw management.

References

  1. Solomon, S.; Qin, D.; Manning, M.; Chen, Z.; Marouis, M.; Averyt, K.B.; Tignor, M.; Miller, H.L. (Eds.) Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: New York, NY, USA, 2007. [Google Scholar]
  2. Ju, H.; van der Velde, M.; Lin, E.D.; Xiong, W.; Li, Y.C. The impacts of climate change on agricultural production systems in China. Clim. Chang. 2013, 120, 313–324. [Google Scholar]
  3. Kogan, F.; Adamenko, T.; Guo, W. Global and regional drought dynamics in the climate warming era. Remote Sens. Lett. 2013, 4, 364–372. [Google Scholar]
  4. Schmidt, L.; Odening, M.; Schlansteina, J.; Ritter, M. Exploring the weather-yield nexus with artificial neural networks. Agric. Syst. 2022, 196, 103345. [Google Scholar]
  5. Mitchell, D.; Heaviside, C.; Vardoulakis, S.; Huntingford, C.; Masato, G.; Guillod, P.B.; Frumhoff, P.; Bowery, A.; Wallom, D.; Allen, M. Attributing human mortality during extreme heat waves to anthropogenic climate change. Environ. Res. Lett. 2016, 11, 74006. [Google Scholar]
  6. Schattman, R.E.; Méndez, V.E.; Merrill, S.C.; Zia, A. Mixed methods approach to understanding farmer and agricultural advisor perceptions of climate change and adaptation in Vermont, United States. Agron. Sustain. Dev. 2018, 42, 121–148. [Google Scholar]
  7. Gong, F.; Wu, X.; Zhang, H.; Chen, Y.; Wang, W. Making better maize plants for sustainable grain production in a changing climate. Front. Plant Sci. 2015, 6, 835. [Google Scholar]
  8. Lesk, C.; Rowhani, P.; Ramankutty, N. Influence of extreme weather disasters on global crop production. Nature 2016, 529, 84–87. [Google Scholar]
  9. Hunt, L.A.; Poorten, G.V.D.; Pararajasingham, S. Postanthesis temperature effects on duration and rate of grain filling in some winter and spring wheats. Can. J. Plant Sci. 1991, 71, 609–617. [Google Scholar]
  10. Sodield, I.; Evans, L.T.; Cook, M.G.; Wardlaw, I.F. Factors influencing the rate and duration of grain filling in wheat. Aust. J. Plant Physiol. 1977, 4, 785–797. [Google Scholar]
  11. Lobell, D.B.; Field, C.B. Global scale climate-crop yield relationships and the impacts of recent warming. Environ. Res. Lett. 2007, 2, 014002. [Google Scholar]
  12. Boisier, J.P.; Ciais, P.; Ducharne, A.; Guimberteau, M. Projected strengthening of Amazonian dry season by constrained climate model simulations. Nat. Clim. Chang. 2015, 5, 656–660. [Google Scholar]
  13. Rusinamhodzi, L.; Corbeels, M.; van Wijk, M.T.; Rufino, M.C.; Nyamangara, J.; Giller, K.E. A meta-analysis of long-term effects of conservation agriculture on maize grain yield under rain-fed conditions. Agron. Sustain. Dev. 2011, 31, 657–673. [Google Scholar]
  14. Lopez-Bellido, L.; Fuentes, M.; Castillo, J.E.; Lopez-Garrido, F.J.; Fernandez, E.J. Long-term tillage, crop rotation, and nitrogen fertilizer effects on wheat yield under Mediterranean conditions. Agron. J. 1996, 88, 783–791. [Google Scholar]
  15. Tao, F.; Yokozawa, M.; Hayashi, Y.; Lin, E. A perspective on water resources in China: Interactions between climate change and soil degradation. Clim. Chang. 2005, 68, 169–197. [Google Scholar]
  16. Schierhorn, F.; Hofmann, M.; Adrian, I.; Bobojonov, I.; Muller, D. Spatially varying impacts of climate change on wheat and barley yields in Kazakhstan. J. Arid Environ. 2020, 178, 104164. [Google Scholar]
  17. Subedi, K.D.; Ma, B.L. Assessment of some major yield-limiting factors on maize production in a humid temperate environment. Field Crops Res. 2009, 110, 21–26. [Google Scholar]
  18. Ceglar, A.; Toreti, A.; Lecerf, R.; Van der Velde, M.; Dentener, F. Impact of meteorological drivers on regional inter-annual crop yield variability in France. Agric. For. Meteorol. 2016, 216, 58–67. [Google Scholar]
  19. Spiertz, H. Challenges for crop production research in improving land use, productivity and sustainability. Sustainability 2013, 5, 1632–1644. [Google Scholar]
  20. Huang, T.; Yang, H.; Huang, C.; Ju, X. Effect of fertilizer N rates and straw management on yield-scaled nitrous oxide emissions in a maize-wheat double cropping system. Field Crops Res. 2017, 204, 1–11. [Google Scholar]
  21. Guo, L.; Li, J.; Li, Y.; Xu, D. Balancing the nitrogen derived from sewage effluent and fertilizers applied with drip irrigation. Water Air Soil Poll. 2017, 228, 12. [Google Scholar]
  22. Meng, Q.F.; Sun, Q.P.; Chen, X.P.; Cui, Z.L.; Yue, S.C.; Zhang, F.S.; Römheld, V. Alternative cropping systems for sustainable water and nitrogen use in the North China Plain. Agric. Ecosyst. Environ. 2012, 146, 93–102. [Google Scholar]
  23. Wang, H.G.; Yu, Z.W.; Zhang, Y.L.; Shi, Y.; Wang, D. Effects of tillage practices on water consumption and dry matter accumulation in dryland wheat. Acta Agron. Sin. 2012, 38, 675–682. (In Chinese) [Google Scholar]
  24. Sharma, R.C.; Dubin, H.J. Effect of wheat cultivar mixtures on spot blotch (Bipolaris sorokiniana) and grain yield. Field Crops Res. 1996, 48, 95–101. [Google Scholar]
  25. Beche, E.; Benin, G.; Silva, C.L.D.; Munaro, L.B.; Marchese, J.A. Genetic gain in yield and changes associated with physiological traits in Brazilian wheat during the 20th century. Eur. J. Agron. 2014, 61, 49–59. [Google Scholar]
  26. Aisawi, K.A.B.; Reynolds, M.P.; Singh, R.P.; Foulkes, M.J. The physiological basis of the genetic progress in yield potential of CIMMYT spring wheat cultivars from 1966 to 2009. Crop Sci. 2015, 55, 1749–1764. [Google Scholar]
  27. Liu, R.J.; Sheng, P.P.; Hui, H.B.; Lin, Q.; Chen, Y.L. Integrating irrigation management for improved grain yield of winter wheat and rhizosphere AM fungal diversity in a semi-arid cropping system. Agric. Syst. 2015, 132, 167–173. [Google Scholar]
  28. Lu, D.; Lu, F.; Pan, J.; Cui, Z.; Zou, C.; Chen, X.; He, M.; Wang, Z. The effects of cultivar and nitrogen management on wheat yield and nitrogen use efficiency in the North China Plain. Field Crops Res. 2015, 171, 157–164. [Google Scholar]
  29. Sun, S.; Yang, X.G.; Lin, X.M.; Sassenrath, G.F.; Li, K.N. Climate-smart management can further improve winter wheat yield in China. Agric. Syst. 2018, 162, 10–18. [Google Scholar]
  30. Bossio, D.; Geheb, K.; Critchley, W. Managing water by managing land: Addressing land degradation to improve water productivity and rural livelihoods. Agric. Water Manag. 2010, 97, 536–542. [Google Scholar]
  31. Bhatt, R.; Khera, K.L. Effect of tillage and mode of straw mulch application on soil erosion in the submontaneous tract of Punjab, Indian. Soil Tillage Res. 2006, 88, 107–115. [Google Scholar]
  32. Pittelkow, C.M.; Liang, X.; Linquist, B.A.; van Groenigen, K.J.; Lee, J.; Lundy, M.E.; van Gestel, N.; Six, J.; Venterea, R.T.; van Kessel, C. Productivity limits and potentials of the principles of conservation agriculture. Nature 2014, 517, 365–368. [Google Scholar]
  33. Nyssen, J.; Poesen, J.; Moeyersons, J.; Haile, M.; Deckers, J. Dynamics of soil erosion rates and controlling factors in the Northern Ethiopian Highlands–towards a sediment budget. Earth Surf. Proc. Landf. 2008, 33, 695–711. [Google Scholar]
  34. Wang, Y.; Xie, Z.; Malhi, S.S.; Vera, C.L.; Zhang, Y.; Wang, J. Effects of rainfall harvesting and mulching technologies on water use efficiency and crop yield in the semi-arid Loess Plateau, China. Agric. Water Manag. 2009, 96, 374–382. [Google Scholar]
  35. Sui, P.; Tian, P.; Lian, H.; Wang, Z.; Ma, Z.; Qi, H.; Mei, N.; Sun, Y.; Wang, Y.; Su, Y.; et al. Straw Incorporation Management Affects Maize Grain Yield through Regulating Nitrogen Uptake, Water Use Efficiency, and Root Distribution. Agronomy 2020, 10, 324. [Google Scholar]
  36. Singh, R.P.; Das, S.K.; Rao, V.M.B.; Reddy, M.N. Towards Sustainable Dryland Agricultural Practices; CRIDA: Hyderabad, India, 1990. [Google Scholar]
  37. Sharma, K.L.; Mandal, U.K.; Srinivas, K.; Vittal, K.P.R.; Mandal, B.; Grace, J.K.; Ramesh, V. Long-term soil management effects on crop yields and soil quality in a dryland Alfisol. Soil Tillage Res. 2005, 83, 246–259. [Google Scholar]
  38. SPSS. PASW Statistics, Version 18.0, PASW Statistics for Windows; SPSS Inc.: Chicago, IL, USA, 2009. [Google Scholar]
  39. Liu, Y.; Sui, Y.W.; Gu, D.D.; Wen, X.X.; Chen, Y.; Li, C.J.; Liao, Y.C. Effects of conservation tillage on grain filling and hormonal changes in wheat under simulated rainfall conditions. Field Crops Res. 2013, 144, 43–51. [Google Scholar]
  40. Latifmanesh, H.; Deng, A.X.; Nawaz, M.M.; Li, L.; Chen, Z.J.; Zheng, Y.T.; Wang, P.; Song, Z.W.; Zhang, J.; Zheng, C.Y.; et al. Integrative impacts of rotational tillage on wheat yield and dry matter accumulation under corn–wheat cropping system. Field Crops Res. 2018, 184, 100–108. [Google Scholar]
  41. Zhang, X.; Sun, N.; Wu, L.; Xu, M.; Bingham, I.J.; Li, Z. Effects of enhancing soil organic carbon sequestration in the topsoil by fertilization on crop productivity and stability: Evidence from long-term experiments with wheat-maize cropping systems in China. Sci. Total Environ. 2016, 562, 247–259. [Google Scholar]
  42. Mekuria, W.; Noble, A.; Sengtaheuanghoung, O.; Hoanh, C.T.; Bossio, D.; Sipaseuth, N.; McCartney, M.; Langan, S. Organic and clay-based soil amendments increase maize yield, total nutrient uptake, and soil properties in Lao PDR. Agroecol. Sustain. Food 2014, 38, 936–961. [Google Scholar]
  43. Van den Putte, A.; Govers, G.; Diels, J.; Gillijns, K.; Demuzere, M. Assessing the effect of soil tillage on crop growth: A meta-regression analysis on European crop yields under conservation agriculture. Eur. J. Agron. 2010, 33, 231–241. [Google Scholar]
  44. Wang, X.B.; Cai, D.X.; Hoogmoed, W.B.; Oenema, O.; Perdok, U.D. Developments in conservation tillage in rainfed regions of North China. Soil Tillage Res. 2007, 93, 239–250. [Google Scholar]
  45. Peng, Z.; Wang, L.; Xie, J.; Li, L.; Coulter, J.A.; Zhang, R.; Luo, Z.; Kholova, J.; Choudhary, S. Conservation Tillage Increases Water Use Efficiency of Spring Wheat by Optimizing Water Transfer in a Semi-Arid Environment. Agronomy 2019, 9, 583. [Google Scholar]
  46. Shao, Y.H.; Xie, Y.X.; Wang, C.Y.; Yue, J.Q.; Yao, Y.Q.; Li, X.D.; Liu, W.X.; Zhu, Y.J.; Guo, T.C. Effects of different soil conservation tillage approaches on soil nutrients, water use and wheat-maize yield in rainfed dry-land regions of North China. Eur. J. Agron. 2016, 81, 37–45. [Google Scholar]
  47. Thierfelder, C.; Chivenge, P.; Mupangwa, W.; Rosenstock, T.S.; Lamanna, C.; Eyre, J.X. How climate-smart is conservation agriculture (CA)?—Its potential to deliver on adaptation, mitigation and productivity on smallholder farms in southern Africa. Food Secur. 2017, 9, 537–560. [Google Scholar]
  48. Qin, R.J.; Stamp, P.; Richner, W. Impact of tillage on maize rooting in a cambisol and luvisol in Switzerland. Soil Tillage Res. 2006, 85, 50–61. [Google Scholar]
  49. Liu, C.; Lu, M.; Cui, J.; Li, B.; Fang, C. Effects of straw carbon input on carbon dynamics in agricultural soils: A meta-analysis. Glob. Chang. Biol. 2014, 20, 1366–1381. [Google Scholar]
  50. Stewart, C.E.; Paustian, K.; Conant, R.T.; Plante, A.F.; Six, J. Soil carbon saturation: Concept, evidence and evaluation. Biogeochemistry 2007, 86, 19–31. [Google Scholar]
  51. Liu, Z.; Gao, T.P.; Liu, W.T.; Sun, K.; Xin, Y.N.; Liu, H.J.; Wang, S.Z.; Li, G.; Han, H.F.; Li, Z.J.; et al. Effects of part and whole straw returning on soil carbon sequestration in C3–C4 rotation cropland. J. Plant Nutr. Soil Sci. 2019, 182, 429–440. [Google Scholar]
  52. Autio, A.; Johansson, T.; Motaroki, L.; Minoia, P.; Pellikka, P. Constraints for adopting climate-smart agricultural practices among smallholder farmers in Southeast Kenya. Agric. Syst. 2021, 194, 103284. [Google Scholar]
  53. Qu, C.H.; Li, X.X.; Ju, H.; Liu, Q. The impacts of climate change on wheat yield in the Huang–Huai–Hai Plain of China using DSSAT–CERES–Wheat model under different climate scenarios. J. Integr. Agric. 2019, 18, 1379–1391. [Google Scholar]
  54. Elhani, S.; Martos, V.; Rharrabti, Y.; Royo, C.; Garcia del Moral, L.F. Contribution of main stem and tillers to durum wheat (Triticum turgidum L. var. durum) grain yield and its components grown in Mediterranean environments. Field Crops Res. 2007, 103, 25–35. [Google Scholar]
  55. Easterling, W.E.; Aggarwal, P.K.; Batima, P.; Brander, K.M.; Erda, L.; Howden, S.M.; Kirilenko, A.; Morton, J.; Soussana, J.F.; Schmidhuber, J.; et al. Food, fibre and forest products. Clim. Chang. 2007, 2007, 273–313. [Google Scholar]
  56. Asseng, S.; Ewert, F.; Martre, P.; Rötter, R.P.; Lobell, D.B.; Cammarano, D.; Kimball, B.A.; Ottman, M.J.; Wall, G.W.; White, J.W.; et al. Rising temperatures reduce global wheat production. Nat. Clim. Chang. 2015, 5, 143–147. [Google Scholar]
  57. Zhao, C.; Liu, B.; Piao, S.; Wang, X.; Lobell, D.B.; Huang, Y.; Huang, M.T.; Yao, Y.T.; 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]
  58. 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]
  59. Lobell, D.B.; Asner, G.P. Climate and management contributions to recent trends in U.S. agricultural yields. Science 2003, 299, 1032. [Google Scholar]
Figure 1. The locations of the Tai’an and Longkou experimental sites.
Figure 1. The locations of the Tai’an and Longkou experimental sites.
Agronomy 12 02452 g001
Figure 2. Temperature and precipitation of the Tai’an and Longkou experimental sites.
Figure 2. Temperature and precipitation of the Tai’an and Longkou experimental sites.
Agronomy 12 02452 g002
Figure 3. Wheat yield and yield components over time under tillage and straw returning practices at the Tai’an (A) and Longkou (B) experimental sites. CT: conventional tillage without straw returning; CTS: conventional tillage with straw returning; RT: rotary tillage without straw returning; RTS: rotary tillage with straw returning. T-like lines: standard deviation. Different letters in each column indicate significant differences between different treatments (p < 0.05; Duncan’s test).
Figure 3. Wheat yield and yield components over time under tillage and straw returning practices at the Tai’an (A) and Longkou (B) experimental sites. CT: conventional tillage without straw returning; CTS: conventional tillage with straw returning; RT: rotary tillage without straw returning; RTS: rotary tillage with straw returning. T-like lines: standard deviation. Different letters in each column indicate significant differences between different treatments (p < 0.05; Duncan’s test).
Agronomy 12 02452 g003
Figure 4. Box plots of equivalent ratios of the wheat yield and yield components under tillage and straw returning practices at the Tai’an (A) and Longkou (B) experimental sites. The horizontal line inside the boxes denotes the median value. The error bars show the upper/lower quartiles (1.5 interquartile distance). The line extending vertically from the box represents the variability beyond the range of the upper and lower quartiles, and the values represent the 10th and 90th percentiles, respectively; the points that lie outside of this range are shown as solid circles. CT: conventional tillage without straw returning; CTS: conventional tillage with straw returning; RT: rotary tillage without straw returning; RTS: rotary tillage with straw returning. Different letters in each group indicate significant differences between different treatments (p < 0.05; Duncan’s test).
Figure 4. Box plots of equivalent ratios of the wheat yield and yield components under tillage and straw returning practices at the Tai’an (A) and Longkou (B) experimental sites. The horizontal line inside the boxes denotes the median value. The error bars show the upper/lower quartiles (1.5 interquartile distance). The line extending vertically from the box represents the variability beyond the range of the upper and lower quartiles, and the values represent the 10th and 90th percentiles, respectively; the points that lie outside of this range are shown as solid circles. CT: conventional tillage without straw returning; CTS: conventional tillage with straw returning; RT: rotary tillage without straw returning; RTS: rotary tillage with straw returning. Different letters in each group indicate significant differences between different treatments (p < 0.05; Duncan’s test).
Agronomy 12 02452 g004
Figure 5. Maize yield and yield components over time under tillage and straw returning practices at the Tai’an (A) and Longkou (B) experimental sites. CT: conventional tillage without straw returning; CTS: conventional tillage with straw returning; RT: rotary tillage without straw returning; RTS: rotary tillage with straw returning. T-like lines: standard deviation. Different letters in each column indicate significant differences between different treatments (p < 0.05; Duncan’s test).
Figure 5. Maize yield and yield components over time under tillage and straw returning practices at the Tai’an (A) and Longkou (B) experimental sites. CT: conventional tillage without straw returning; CTS: conventional tillage with straw returning; RT: rotary tillage without straw returning; RTS: rotary tillage with straw returning. T-like lines: standard deviation. Different letters in each column indicate significant differences between different treatments (p < 0.05; Duncan’s test).
Agronomy 12 02452 g005
Figure 6. Box plots of equivalent ratios of the maize yield and yield components under tillage and straw returning practices at the Tai’an (A) and Longkou (B) experimental sites. The horizontal line inside the boxes denotes the median value. The error bars show the upper/lower quartiles (1.5 interquartile distance). The line extending vertically from the box represents the variability beyond the range of the upper and lower quartiles, and the values represent the 10th and 90th percentiles, respectively; the points that lie outside of this range are shown as solid circles. CT: conventional tillage without straw returning; CTS: conventional tillage with straw returning; RT: rotary tillage without straw returning; RTS: rotary tillage with straw returning. Different letters in each group indicate significant differences between different treatments (p < 0.05; Duncan’s test).
Figure 6. Box plots of equivalent ratios of the maize yield and yield components under tillage and straw returning practices at the Tai’an (A) and Longkou (B) experimental sites. The horizontal line inside the boxes denotes the median value. The error bars show the upper/lower quartiles (1.5 interquartile distance). The line extending vertically from the box represents the variability beyond the range of the upper and lower quartiles, and the values represent the 10th and 90th percentiles, respectively; the points that lie outside of this range are shown as solid circles. CT: conventional tillage without straw returning; CTS: conventional tillage with straw returning; RT: rotary tillage without straw returning; RTS: rotary tillage with straw returning. Different letters in each group indicate significant differences between different treatments (p < 0.05; Duncan’s test).
Agronomy 12 02452 g006
Figure 7. Wheat and maize yields over time under no straw return and straw return practices at the Tai’an (A) and Longkou (B) experimental sites. ** significant difference at 0.01 level; * significant difference at 0.05 level; ns: not significant difference.
Figure 7. Wheat and maize yields over time under no straw return and straw return practices at the Tai’an (A) and Longkou (B) experimental sites. ** significant difference at 0.01 level; * significant difference at 0.05 level; ns: not significant difference.
Agronomy 12 02452 g007
Figure 8. Yield and yield component stability of wheat and maize under tillage and straw returning practices at the Tai’an (A) and Longkou (B) experimental sites. CT: conventional tillage without straw returning; CTS: conventional tillage with straw returning; RT: rotary tillage without straw returning; RTS: rotary tillage with straw returning. T-like lines: standard deviation. Different letters in each group indicate significant differences between different treatments (p < 0.05; Duncan’s test).
Figure 8. Yield and yield component stability of wheat and maize under tillage and straw returning practices at the Tai’an (A) and Longkou (B) experimental sites. CT: conventional tillage without straw returning; CTS: conventional tillage with straw returning; RT: rotary tillage without straw returning; RTS: rotary tillage with straw returning. T-like lines: standard deviation. Different letters in each group indicate significant differences between different treatments (p < 0.05; Duncan’s test).
Agronomy 12 02452 g008
Figure 9. Correlation between wheat yield, yield components and mean temperature and precipitation in the wheat season. The blue circles, lines, and formulas represent the data of the Tai’an experimental site, and the pink circles and lines represent the data of the Longkou experimental site. r: correlation coefficient; p: significance level; n: sample size; *: Significant at p < 0.05. ns: not significant.
Figure 9. Correlation between wheat yield, yield components and mean temperature and precipitation in the wheat season. The blue circles, lines, and formulas represent the data of the Tai’an experimental site, and the pink circles and lines represent the data of the Longkou experimental site. r: correlation coefficient; p: significance level; n: sample size; *: Significant at p < 0.05. ns: not significant.
Agronomy 12 02452 g009
Figure 10. Correlation between maize yield, yield components and mean temperature and precipitation in the maize season. The blue circles, lines, and formulas represent the data of the Tai’an experimental site, and the pink circles and lines represent the data of the Longkou experimental site. r: correlation coefficient; p: significance level; n: sample size; **: significant at p < 0.01; *: significant at p < 0.05. ns: not significant.
Figure 10. Correlation between maize yield, yield components and mean temperature and precipitation in the maize season. The blue circles, lines, and formulas represent the data of the Tai’an experimental site, and the pink circles and lines represent the data of the Longkou experimental site. r: correlation coefficient; p: significance level; n: sample size; **: significant at p < 0.01; *: significant at p < 0.05. ns: not significant.
Agronomy 12 02452 g010
Figure 11. Correlation between wheat and maize yield and mean temperature trend and precipitation trend. The trend represents the difference of the following year minus the previous year. The blue circles and lines represent the data of the Tai’an experimental site, and the pink circles, lines, and formulas represent the data of the Longkou experimental site. r: correlation coefficient; p: significance level; n: sample size; **: significant at p < 0.01; ns: not significant.
Figure 11. Correlation between wheat and maize yield and mean temperature trend and precipitation trend. The trend represents the difference of the following year minus the previous year. The blue circles and lines represent the data of the Tai’an experimental site, and the pink circles, lines, and formulas represent the data of the Longkou experimental site. r: correlation coefficient; p: significance level; n: sample size; **: significant at p < 0.01; ns: not significant.
Agronomy 12 02452 g011
Figure 12. Direct and indirect effects of mean temperature, annual precipitation, and yield components on crop yield. Structural equation modeling was conducted for wheat yield (A) and maize yield (B). Single-headed arrows indicate the hypothesized direction of causation. Red and blue arrows indicate positive and negative relationships, respectively. The arrow width was proportional to the strength of the relationship. The numbers adjacent to the arrows are the standardized path coefficients. Significant differences are indicated: *, p<0.05; **, p<0.01; ***, p<0.001.
Figure 12. Direct and indirect effects of mean temperature, annual precipitation, and yield components on crop yield. Structural equation modeling was conducted for wheat yield (A) and maize yield (B). Single-headed arrows indicate the hypothesized direction of causation. Red and blue arrows indicate positive and negative relationships, respectively. The arrow width was proportional to the strength of the relationship. The numbers adjacent to the arrows are the standardized path coefficients. Significant differences are indicated: *, p<0.05; **, p<0.01; ***, p<0.001.
Agronomy 12 02452 g012
Table 1. Mean yield of wheat and maize under different tillage and straw returning from 2010 to 2019.
Table 1. Mean yield of wheat and maize under different tillage and straw returning from 2010 to 2019.
Warm-Dry Site Cool-Wet Site
Wheat Yield (t ha−1)Maize Yield (t ha−1)Wheat Yield (t ha−1)Maize Yield (t ha−1)
CT7.7 a10.0 b8.1 a8.4 b
CTS7.9 a11.0 a8.6 a8.9 ab
RT7.0 b9.7 b8.3 a8.7 ab
RTS8.2 a10.3 ab8.7 a9.3 a
CT: conventional tillage without straw returning; CTS: conventional tillage with straw returning; RT: rotary tillage without straw returning; RTS: rotary tillage with straw returning. Different letters in each column indicate significant differences between different treatments (p < 0.05; Duncan’s test).
Table 2. Correlation between crop yield and yield components.
Table 2. Correlation between crop yield and yield components.
Wheat YieldMaize Yield
Ear density–0.301 **0.738 **
Kernel number0.244 *0.867 **
Thousand kernel weight0.724 **0.460 **
*: Significant at p < 0.05; **: Significant at p < 0.01.
Table 3. Three-way analysis of variance of the effects of long-term experiments on crop yield and yield components.
Table 3. Three-way analysis of variance of the effects of long-term experiments on crop yield and yield components.
FactorWheatMaize
Ear DensityKernel NumberThousand-Kernel WeightYieldEar DensityKernel NumberThousand-Kernel WeightYield
DFpDFpDFpDFpDFpDFpDFpDFp
L10.00010.41710.00010.00010.15310.00010.12110.000
T10.43110.77310.38510.98910.21510.55510.58010.727
S10.31110.50010.02310.00110.71610.09510.07410.001
L × T10.28510.82810.80010.18610.11110.40410.77910.011
L × S10.58910.51710.02710.29210.90810.83510.63010.376
T × S10.85810.18710.58010.10610.70710.93710.42710.516
L × T × S10.71110.16310.44010.08910.54310.84410.84210.470
L: experimental site; T: tillage practice; S: straw management; DF: degree of freedom; p: significance level.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, Z.; Wang, N.; Lü, J.; Wang, L.; Li, G.; Ning, T. Climate-Smart Tillage Practices with Straw Return to Sustain Crop Productivity. Agronomy 2022, 12, 2452. https://doi.org/10.3390/agronomy12102452

AMA Style

Liu Z, Wang N, Lü J, Wang L, Li G, Ning T. Climate-Smart Tillage Practices with Straw Return to Sustain Crop Productivity. Agronomy. 2022; 12(10):2452. https://doi.org/10.3390/agronomy12102452

Chicago/Turabian Style

Liu, Zhen, Ning Wang, Jinling Lü, Lan Wang, Geng Li, and Tangyuan Ning. 2022. "Climate-Smart Tillage Practices with Straw Return to Sustain Crop Productivity" Agronomy 12, no. 10: 2452. https://doi.org/10.3390/agronomy12102452

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