Dynamic Linkages among Climate Change, Mechanization and Agricultural Carbon Emissions in Rural China

Climate change has become a major environmental issue facing all countries, having a significant effect on all aspects of agricultural production, such as the agricultural mechanization process and fertilizer use. Greenhouse gases produced by agricultural machinery and fertilizers during agricultural production are an important cause of climate change. On the basis of the above facts, researching the connection between agricultural mechanization, climate change, and agricultural carbon emissions is crucial for the development of low-carbon agriculture and for addressing climate change. We used a variety of econometric models and methods to analyze data from China’s multiple provinces (cities) covering the years 2000 through 2019, in order to meet the research objectives. Furthermore, we utilized rainfall and sunlight as variables to assess climate change and adopted Granger tests to establish the link between rainfall, sunlight, agricultural mechanization, and carbon emissions in farming. The findings indicate a bidirectional causality relationship between rainfall, sunlight, agricultural mechanization, and carbon emissions in farming. Rainfall and sunlight are Granger causes of agricultural mechanization. Furthermore, agricultural mechanization has favorable effects on carbon emissions of agriculture, and climate change has long-term implications on agricultural mechanization and carbon emissions of agriculture. Finally, this paper investigated the green path suitable for the low-carbon development of Chinese agriculture, arguing that the government should formulate low-carbon agricultural policies by region and actively promote the upgrading of agricultural machinery.


Introduction
Climate change has become a common problem faced by all countries. During the 20th century, the average global temperature rose about 1 degree Celsius [1], which led to sea level rise [2], species extinction [3], and frequent climate extremes [4]. China's carbon dioxide emissions hit 11.9 billion tons in 2021, making up 33 percent of the world's total and far outpacing those of other nations and regions, according to the International Energy Agency (IEA). As a consequence, controlling China's carbon emissions is vital for mitigating climate change. The Chinese government has also shown its determination to tackle global warming. China has pledged to peak its carbon emissions by 2030 and become carbon neutral by 2060, which is referred to as the "dual carbon" target. In China, about 20 percent of carbon emissions are attributable to agricultural production and land use [5], and agriculture-related carbon emissions make a contribution of about a quarter of the globe's total [6]. Due to the extensiveness and universality of China's agricultural production activities and the dispersion of agricultural production entities, coupled with the wide range, randomness, concealment, difficulty in monitoring, and quantification of agricultural carbon emissions, China's agricultural carbon emissions are difficult to control. Meanwhile, China's agricultural production is large and growing rapidly. If there are no corresponding reduction measures, agricultural greenhouse gas emissions will continue to rise rapidly. Therefore, in order to achieve the "double carbon" goal, China must find a way to reduce agricultural carbon emissions as soon as possible, as well as realize the transformation of agricultural production to green, low-carbon, and sustainable development.
Most agricultural GHG (greenhouse gas) emissions currently come from fertilizer, feed, pesticides, and energy consumption, as revealed in the research conducted by World Food and Agriculture Organization [7]. Numerous studies have shown that climate change increases the volatility of agricultural production and has a differential impact on different regions, but climate change damages agricultural production in general [8]. Farmers and the government have also taken corresponding actions to comply with climate variation. Examples include cultivating arable land and increasing the use of chemical fertilizers to cope with agricultural production reduction brought by bad weather [9], increasing the use of pesticides to address pests and diseases caused by temperature rise [10], and investing in more mechanical equipment and changing irrigation methods to deal with extreme droughts and floods [11]. Each country should use agricultural machinery on the basis of climatic conditions [12] and guide the agricultural mechanization following its requirements. Therefore, the direction and pace of agricultural mechanization will be influenced by climate change.
From 1980 to 2015, China's agricultural mechanization continued to deepen, which enhanced agricultural productivity [13]. However, due to the uncentralized distribution of agricultural land in China, the use efficiency of agricultural machinery was hampered. Giant energy use was caused by the inefficient operation of equipment, which aggravates carbon emissions to a large extent [14]. In the last few years, energy consumption brought about by agricultural mechanization has gradually become a crucial source of carbon emissions in farming production [15]. Therefore, investigating the connection between agricultural carbon emissions, climate change, and agricultural mechanization in the backdrop of "dual carbon" today can make up for the vacancy in research in related fields and provide suggestions on measures to promote the growth of low-carbon agriculture for the government.
Accordingly, China Statistical Yearbook and China Rural Statistical Yearbook covering 2000 to 2019 are the data sources for the analysis of this study. We chose rainfall and sunlight as variables to measure climate change, and then used ARDL (autoregressive distributed lag model) and PVAR (panel data vector auto regression) to look into the connection between climate change, agricultural mechanization, and agricultural carbon emissions.
The key contributions of this paper in comparison to earlier material are listed below. First, climate change and agricultural carbon emissions have been hot issues in recent years, but there are few empirical studies on the interaction between them. At the same time, most of the literature adopts only one primary explanatory factor, while the empirical study in this paper ingeniously integrates multiple variables including agricultural mechanization, climate change, and agricultural carbon emissions, demonstrating the interaction and transmission mechanism between them. Second, the dynamic test framework established in this paper considers the cross-sectional dependence, co-integration relationship, hysteresis effect, and different mutual inspection methods. The purpose is to guarantee the precision and reliability of the research results and provide a reference for the subsequent research of other scholars. Finally, this paper puts forward recommendations for advancing energysaving and low-carbon agriculture, as well as establishing low-carbon agricultural policies.

Literature Review
A substantial amount of carbon emissions are now produced by China's agriculture, drawing the interest of numerous academics. They have discussed the elements that influence carbon emissions in China's farming sector, believing that human activities such as pesticides, fertilization, irrigation, mechanization, and mulch coverage significantly contribute to carbon emissions in farming production [16]. In recent years, due to extreme droughts and floods caused by climate change, mechanization has become an indispensable factor in agricultural carbon emissions [17]. Mechanized equipment has significantly improved labor productivity and resource utilization, further contributing to the agricultural technological progress and output growth [18]. However, agricultural machinery is mainly driven by fossil fuels, and the gas generated by fossil fuel combustion significantly harms the environment [19]. Unfortunately, there has been little reference to climate change in the discussion of the factors that contribute to agricultural carbon emissions.
Comprehensive research has already been conducted on how climate change would affect the agricultural industry. By and large, climate change has uncertain impacts on agriculture in different regions, but mainly adverse effects [20]. Grain security is seriously threatened by climate variation since agricultural production depends on relatively consistent weather patterns [21]. Changes in rainfall patterns, floods, droughts, temperatures, and pests and diseases have hindered the development of agriculture, adversely affecting crop production and yields [22][23][24]. Additionally, climate change affects the direction and process of agricultural mechanization. The continuous severe weather is unfavorable to the input and use of machinery and equipment [25]. According to studies, climate change in northern Norway has made spring sowing windows even more constrained, affecting soil working days and farmers' management of agricultural mechanization [26]. Therefore, countries should use agricultural machinery based on climate conditions and guide mechanization following its requirements. Although the existing literature has confirmed the impact of climate change on the agriculture department, few scholars have looked into the connection between climate change and agricultural carbon emissions.
However, many scholars have found that climate change has affected the factors contributing to agricultural carbon emissions. We can roughly divide the factors that lead to the connection between climate change and agricultural carbon emissions into natural and artificial factors. In terms of natural factors, climate change considerably influences the carbon sequestration efficiency of farm production. Although crop cover can slow down the erosion rate of heavy rainfall on soil, the increase in rainfall poses a threat to the extent of soil erosion in farmland [27]. On the one hand, soil erosion can destroy soil's physical and chemical properties, resulting in the loss of soil nutrients and the decline of carbon sequestration. On the other hand, soil and water loss caused by heavy rainfall leads to the accumulation of sediments carrying large amounts of chemical fertilizers and pesticide residues in rivers and reservoirs, which adversely affects the carbon cycle in the atmosphere [28]. From the perspective of human factors, farmers have taken corresponding agricultural production measures according to the changing climate conditions [29], which may indirectly affect carbon emissions. Soil moisture deficiency caused by reduced rainfall can reduce the biological function of crops, which is more likely to induce pests and diseases [24]. Increased pesticide use can reduce agricultural output losses caused by illnesses and insects, but it can be damaging to the environment [30]. In addition, researchers found that adverse climatic conditions affected soil processing conditions in certain parts of the UK, leading to a corresponding reduction in the use of farm machinery [31]. Moreover, farmers may extract more groundwater for agricultural irrigation to ease the pressure on water use from changing temperature and rainfall patterns. However, groundwater extraction consumes energy and exacerbates agricultural carbon emissions [32]. Farmers may also use drip irrigation, sprinkler irrigation, and other water-saving irrigation techniques in agricultural production to cope with climate change. However, sprinkler irrigation technology requires high water pressure, and the application of irrigation equipment may increase GHG emissions [33]. Therefore, climate change is closely related to agricultural carbon emissions through human adaptation behavior, especially the input of mechanized equipment. Overall, we should not ignore the link between carbon emissions and climate change.
Through the literature, we note that the existing research focuses on the effects of agricultural mechanization on agricultural carbon emissions and the effects of climate change on agricultural mechanization. Although a great deal of research has confirmed the association between climate change and agricultural carbon emissions, few scholars have taken both climate change and mechanization into account when studying the factors that contribute to carbon emissions. Additionally, previous studies mostly used single provincial, municipal or regional data when examining the effects of mechanization or climate change on agricultural departments, but few studies chose data from multiple provinces or regions of a country. Therefore, this paper thoroughly examined the data covering China's multiple provinces (cities); innovatively took rainfall and sunlight as variables to measure climate change; and incorporated agricultural mechanization, climate change, and agricultural carbon emissions into a homogeneous measurement framework. The purpose was to look into the long-term relationship between agricultural mechanization, climate change, and carbon emissions in agriculture, especially the relationship between the latter two, in order to fill the gap in existing research and propose countermeasures and suggestions for sustainable growth of ecological agriculture.

Calculation of Agricultural Carbon Emission
Agricultural carbon emissions (ACE) chiefly include fossil fuels, deforestation, burning, fertilizer use, arable land, and animal digestion and excretion [34]. This study focused on carbon emissions from agriculture diesel inputs, selecting fertilizers, pesticides, agricultural irrigation, agricultural electricity, agricultural farming, etc. The collected data were then multiplied by the usage of carbon sources and carbon emission factors using the carbon emission factor method to derive the total agricultural carbon emission of each province (municipality). At the same time, in order to avoid the length of the article, the empirical methods and results of this paper used ACE to represent agricultural carbon emissions.
where E represents the total agricultural carbon emissions, E i represents the agricultural carbon emissions of different carbon emission sources, T i represents the carbon source consumption, and δ i represents the carbon emission factor. Table 1 lists the emission coefficients of each carbon source for the readers' reference.

Data and Variables
In this paper, agricultural mechanization intensity (total power of agricultural machinery/cultivated land area) was used to measure the degree of agricultural mechanization. Climate change will have many impacts on the Earth's environment, including rainfall and the duration of sunlight [37], which play an important role in agriculture [38]. The direct effect of these two factors on agricultural carbon emissions is relatively rare in previous studies, and the correlation between the two is small (compared to other climate variables such as temperature). Therefore, this paper chose rainfall and sunlight as variables representing climate change for further empirical study. This study adopted the annual data of 30 provinces (municipalities) from the China Rural Statistical Yearbook (2000-2019) and the China Statistical Yearbook (2000-2019) for empirical analysis. To carry out the research smoothly, this paper set the following variables on the basis of previous literature: agricultural carbon emission (PERAGCARBON), rainfall (RAIN), sunlight (SUNLIGHTS), and agricultural mechanization index (MACHINE). To alleviate heteroscedasticity, we analyzed the data by logarithmic processing. In addition, we provide readers with the definitions of the above variables in Table 2. Table 2. Definition of variables.

Variables Definition
Carbon emissions from agricultural production (PERAGCARBON) Average carbon emissions from agricultural production

Rainfall (RAIN)
The average total rainfall of each province (municipality) Sunlight intensity (SUNSUNLIGHTS) The average total sunlight days of each province (city)

Index of mechanization (MACHINE)
Average agricultural mechanization intensity of each province (municipality) Table 3 shows descriptive statistics for each variable from 2004 to 2019. The average carbon emission was 0.313025, and the interval was 1.055501, which indicates that the data had volatility. At the same time, the average rainfall and sunlight were also different from the interval. Thus, the climate conditions in the 30 provinces (municipalities) were diverse. At the same time, there were differences in agricultural carbon emissions and agricultural mechanization in 30 provinces (cities), which could have a symphony for farmers' income and planting area. At the same time, the standard deviation of the original data was too large, so we logged the data to reduce the impact of spurious regression and heteroscedasticity The standard deviation after the logarithm of the data was larger, indicating that the data were more dispersed.

Test for Cross-Sectional Dependence
One of the characteristics of panel data is that there is a varying degree of correlation between cross-sections, leading to the correlation between cross-sectional heterogeneity or regression errors, which affects the unbiasedness, consistent, and valid estimation of standard panel data [39]. Therefore, before panel testing, cross-sectional correlation tests should be performed in order to address cross-section correlation [40].
The Breusch-Pagan Lagrange multiplier (LM) test was proposed by Breusch and Pagan in 1980 to test cross-sectional correlation. Pesaran improved the disadvantages of the Breusch-Pagan LM test and proposed the Pesaran cross-sectional dependence (CD) and Pesaran LM tests. Subsequently, some scholars have proposed the bias-corrected LM test and extended it to the dynamic panel data model [41]. We provide the calculation formula of the section test method in Appendix A for the readers' reference.

Panel Unit Root Tests
When common least squares are used, the presence of horizontal unit roots in panel data can result in pseudo-regression, leading to false estimates. Therefore, the variable unit root tests should be performed first. In addition, the unit root of a non-stationary sequence can be eliminated by difference [42]. On the basis of previous studies, this paper used the Levin-Lin-Chu test (LLC), the Im-Pesaran-Shin test (IPS), Fisher's augmented Dickey-Fuller (ADF) test, the Breitung t-stat method, and Fisher's Phillip and Perron (PP) test in order to examine the panel unit root. We provide the specific process of these five methods in Appendix B for the readers' reference.

Panel Cointegration Test
The paper used the Kao test [43] to conduct a co-integration test on panel data. According to the test results, we can judge if there is a long-term relationship between variables. This study lists two main steps of the test: Step 1: Set each section to have the same coefficient and different intercept: In the above equation, α i represents a single intercept parameter, and e it represents the residual.
Step 2: Perform the unit root test on the residual sequence of step 1. The null hypothesis is H 0 : ρ = 1, and the following statistics are constructed:

Granger Causality Test
After the above research, we focused on the causal relationship between the variables. Granger [44] first proposed a method to analyze the causality of time series data. On this basis, Dumitrescu and Hurlin [45] improved it. This paper adopted the improved Granger causality test of Dumitrescu and Herlin. This paper lists the formulas for the causality test: where x i,t and y i,t are the observed value of individual i in period t, respectively, and k represents the number of lags of the individual.

Autoregressive Distributed Lag
To better understand the dynamic relationships between the variables studied, an autoregressive distributed lag model (ARDL) proposed by Pesaran [46] was used. The ARDL model has many advantages, such as dealing with different delay-order variables and analyzing statistical regression and other common models. The typical ARDL models are as follows: In the above equation, P represents the order of lag y t and q i represents the order of lag of the i-th independent variable x it . L is the lag operator, which can be defined as follows: Ly t = y t−1 , and w t is the determination vector of S row 1 column.

FMOLS and DOLS
The paper used two static models, FMOLS (full modified ordinary least squares) [47] and DOLS (dynamic ordinary least squares) [48], to ensure the robustness of the results in Section 3.7. Compared with the OLS (orthogonal least square method) model, these two models construct a second-order bias consisting of endogenous bias and non-central bias.
The panel FMOLS estimator β is given by The DOLS is written as follows: where the estimated coefficient β is given by where Z it = X it − X l , ∆X it−j , . . . , ∆X it+k is 2 (K + 1) vectors of regressors.

Impulse Response Approach and Variance Decomposition
To further study the dynamic influence of each variable, this paper adopted impulse response and variance decomposition analysis to grasp the mutual influence and action of each variable dynamically. The impulse response function refers to the impact trajectory brought about by one standard deviation of random disturbance terms on the current and future values of itself and other endogenous variables, which can describe the dynamic influence among variables more intuitively. The variance decomposition can determine the influence of explanatory variables and other variables on endogenous variables by the degree of contribution of each shock to endogenous variables. Lanne [49] proposes that the estimated impulse response exists in the vector autoregressive (VAR) model as follows: (12) is the impulse response function, which can be transformed into infinite moving regression estimation by the following equation: A j is a coefficient matrix that transforms VAR into an infinite vector moving average, k is the optimal lag term, and ε t represents the error term.
Equation (15) represents the contribution of variable N to the prediction error variance of variable M in period h.

Cross-Sectional Dependence Tests Results
The results of cross-sectional correlation tests in Tables 4 and 5 show that the null hypothesis of no cross-sectional correlation at the significance level of 1‰ was rejected by all four tests. Therefore, the cross-section of panel data was correlated.

Unit Root Test Results
LNRAIN and LNSUNLIGHT both rejected the null hypothesis at a significance level of 1% in each of the five trials, but the other two variables only rejected the null hypothesis in the LLC trial. Therefore, we performed a first-order differential and tested the data, wherein the final result rejected the zero hypothesis for the unit root of each variable at the critical 1% level. However, this suggests that there may be a spurious regression and therefore we need to use the KAO test for cointegration. We show the results of the five unit root test methods in Appendix C for reference. For details see Table A1 The results of panel unit root tests. Tables 6 and 7 show the co-integration test results of rainfall, agricultural mechanization, and ACE, and the co-integration test results of sunlight, agricultural mechanization, and ACE, respectively. Both groups rejected the original hypothesis, indicating that the panel data existed in a cointegration relationship. The results confirm a long-run equilibrium cause-and-effect relationship between these two sets of variables and help us further study the effects of climate change and agricultural mechanization on ACE.  Tables 8 and 9 show the estimation results of FMOLS and DOLS panel models. From the level of explicitness of parameters, DOLS had a better fitting effect. Thus, we know that when the agricultural mechanization index increases by 1%, ACE will increase by 0.24%. When rainfall increases by 1%, ACE will increase by 0.88%. When sunlight increases by 1%, ACE will decrease by 0.61%, indicating that rainfall and agricultural mechanization have a certain positive correlation with agricultural carbon emissions, while sunlight is negatively correlated with ACE.

Results of ARDL
We used the Akaike information criteria (AIC) to determine the optimal lag length of the model. Finally, the model for agricultural carbon emissions, agricultural mechanization, and rainfall was ARDL (2, 3, 3), while the model for agricultural carbon emissions, agricultural mechanization, and sunlight was ARDL (1, 1, 1).
The long-term relationship of variables obtained from the ARDL model is shown in Table 10. In the long term, agricultural mechanization and rainfall will increase ACE, consistent with the results shown in Section 4.4. In the short term, looking at the first-order difference data, the results showed that agricultural mechanization will promote ACE. However, the impact of rainfall on ACE in the short term is not strong, indicating that the impact of rainfall on ACE may be a long-term process. The effect of agricultural mechanization on ACE in Table 11 is similar to the conclusion in Table 11. Sunlight can inhibit ACE in the long and short term, but the impact of shortterm sunlight on ACE is not obvious, and thus it is speculated that the impact of sunlight on ACE is a long-term process, such as with rainfall.

Granger Causality Tests
Granger causality exists in the previous cointegration, but we cannot guarantee that the causal relationship between variables was identified. Therefore, we further conducted the Granger causality test. Tables 12 and 13 show the results of the causality test among several groups of variables. The null hypothesis is that there is no Granger causality between variables. We generated Figure 1 to facilitate readers' understanding of the causal relationship between variables.
We found bidirectional causality between agricultural mechanization and ACE, bidirectional causality between rainfall and ACE, and bidirectional causality between sunlight and ACE at the 5% significance level. In addition, sunlight and rainfall were unidirectional Granger causes of agricultural mechanization at the 5% level of significance.  We found bidirectional causality between agricultural mechanization and ACE, bidirectional causality between rainfall and ACE, and bidirectional causality between sunlight and ACE at the 5% significance level. In addition, sunlight and rainfall were unidirectional Granger causes of agricultural mechanization at the 5% level of significance.

Analysis Results of Rainfall, Mechanization, and Agricultural Carbon Emissions
The optimal lag order for agricultural mechanization, rainfall, and agricultural carbon emissions should be determined before VAR systems are used to analyze their pulse effect and variance decomposition as endogenous variables. In this paper, LR test statistic (LR), final prediction error (FPE), Akaike information criterion (AIC), Schwarz information criterion (SIC), and Hannan-Quinn information criterion (HQ) were the five methods for comprehensive judgment. As shown in Table 14, it can be seen that lag order 2 was selected as the optimal lag term. According to this sequence, Figure 2 was produced. It can be seen from Figure 2 that each root was within the unit circle, so this VAR model had the conditions for variance decomposition analysis and impulse response analysis.

Analysis Results of Rainfall, Mechanization, and Agricultural Carbon Emissions
The optimal lag order for agricultural mechanization, rainfall, and agricultural carbon emissions should be determined before VAR systems are used to analyze their pulse effect and variance decomposition as endogenous variables. In this paper, LR test statistic (LR), final prediction error (FPE), Akaike information criterion (AIC), Schwarz information criterion (SIC), and Hannan-Quinn information criterion (HQ) were the five methods for comprehensive judgment. As shown in Table 14, it can be seen that lag order 2 was selected as the optimal lag term. According to this sequence, Figure 2 was produced. It can be seen from Figure 2 that each root was within the unit circle, so this VAR model had the conditions for variance decomposition analysis and impulse response analysis.  Note: *** indicates that this column has 10% visibility. The impulse response function can visually describe the VAR model of a standard deviation of the random disturbance impact on other variables' trajectories and the effect of the influence of current and future values. Therefore, this paper further analyzed the impact of the changes in ACE, agricultural mechanization, and precipitation on the other two variables through an impulse response function diagram. The response time set in this paper was 15 years. In Figure 3, longitudinal coordinates indicate the degree of response of endogenous variables to an impact, the abscissa indicates the lag length of the impact, and dotted lines on either side of the solid lines indicate the range of possible impulse response. The results in Figure 3 are as follows: (1) The impulse response to the impacts of ACE was positive and large, attenuated slightly in the first phase and then was stabilized and sustained. The impulse response of ACE to rainfall and agricultural mechanization was positive in the initial stage, and it turned negative in the fourth stage and converged to 0 in the long run. This indicates that the increase in both will promote the increase in ACE in a short amount of time. (2) The impulse response of agricultural mechanization to ACE was positive in the initial stage, and then tended towards zero after reaching the maximum in the second stage.
The impulse response of agricultural mechanization to rainfall was negative at the beginning, turned positive after several periods, and then approached zero, indicating that the increase in ACE and the decrease in rainfall will promote agricultural mechanization to some extent. The impulse response function can visually describe the VAR model of a standard deviation of the random disturbance impact on other variables' trajectories and the effect of the influence of current and future values. Therefore, this paper further analyzed the impact of the changes in ACE, agricultural mechanization, and precipitation on the other two variables through an impulse response function diagram. The response time set in this paper was 15 years. In Figure 3, longitudinal coordinates indicate the degree of response of endogenous variables to an impact, the abscissa indicates the lag length of the impact, and dotted lines on either side of the solid lines indicate the range of possible impulse response. The results in Figure 3 are as follows: (1) The impulse response to the impacts of ACE was positive and large, attenuated slightly in the first phase and then was stabilized and sustained. The impulse response of ACE to rainfall and agricultural mechanization was positive in the initial stage, and it turned negative in the fourth stage and converged to 0 in the long run. This indicates that the increase in both will promote the increase in ACE in a short amount of time. (2) The impulse response of agricultural mechanization to ACE was positive in the initial stage, and then tended towards zero after reaching the maximum in the second stage. The impulse response of agricultural mechanization to rainfall was negative at the beginning, turned positive after several periods, and then approached zero, indicating that the increase in ACE and the decrease in rainfall will promote agricultural mechanization to some extent. (3) The initial impulse response of rainfall to agricultural mechanization and ACE was weak, but in the long run, rainfall had a negative response to ACE. This suggests that rainfall is almost unaffected by agricultural mechanization and that ACE has a lagged impact on rainfall.
(3) The initial impulse response of rainfall to agricultural mechanization and ACE was weak, but in the long run, rainfall had a negative response to ACE. This suggests that rainfall is almost unaffected by agricultural mechanization and that ACE has a lagged impact on rainfall.  On the basis of the above analysis, the importance of each structural impact on endogenous variables can be assessed by variance decomposition in VAR models. The calculation results of variance decomposition of the VAR model of ACE, rainfall, and agricultural mechanization are shown in Table 15.  On the basis of the above analysis, the importance of each structural impact on endogenous variables can be assessed by variance decomposition in VAR models. The calculation results of variance decomposition of the VAR model of ACE, rainfall, and agricultural mechanization are shown in Table 15. Agricultural carbon emissions are mainly affected by self-impact, followed by agricultural mechanization. With a lag of 15 cycles, the variance contribution rates of rainfall and mechanization were 0.44% and 1.88%, respectively. Agricultural mechanization was mainly affected by its impact and gradually weakened. The influence of rainfall gradually became stronger, indicating that the influence of rainfall on it was a long-term process. Rainfall was mainly affected by its impact, followed by agricultural carbon emissions, but not obviously.

Analysis Results of Sunlight, Mechanization, and Agricultural Carbon Emissions
The optimal lag order for agricultural mechanization, sunlight, and agricultural carbon emissions should be determined before VAR systems can be used to analyze their pulse effect and variance decomposition as endogenous variables. The same five methods used in Section 4.7.1 were used here. As shown in Table 16, it can be seen that the third lag order was selected as the optimal lag term. It can be seen from Figure 4 that each root was within the unit circle, and therefore this VAR model had the conditions for variance decomposition analysis and impulse response analysis.

Analysis Results of Sunlight, Mechanization, and Agricultural Carbon Emissions
The optimal lag order for agricultural mechanization, sunlight, and agricultural carbon emissions should be determined before VAR systems can be used to analyze their pulse effect and variance decomposition as endogenous variables. The same five methods used in Section 4.7.1 were used here. As shown in Table 16, it can be seen that the third lag order was selected as the optimal lag term. It can be seen from Figure 4 that each root was within the unit circle, and therefore this VAR model had the conditions for variance decomposition analysis and impulse response analysis.   The assumptions in Figure 5 are similar to those in Figure 3 in Section 4.7.2 and are not detailed here. The results in Figure 4 are as follows: The impulse response of ACE to its impact was positive and large, and the impulse response of ACE to agricultural mechanization was consistent with the analysis in Section 4.7.2. However, the pulsing response of ACE to sunlight was negative at first, positive over several periods, and then tended towards zero. This was consistent with the results of the FMOLS and DOLS tests, indicating that sunlight can reduce ACE. The impulse response of sunlight to ACE was always negative, which was strong in the first few periods, and then decreased. These results indicate that ACE would have a certain weakening effect on sunlight. The impulse response of sunlight duration to agricultural mechanization was positive at the beginning, turned negative after several periods, and then approached zero, indicating that sunlight duration has a promoting effect on agricultural mechanization. Table 17 shows the calculation results of variance decomposition of the VAR model for agricultural carbon emission, sunlight, and agricultural mechanization.  The impulse response of ACE to its impact was positive and large, and the impulse response of ACE to agricultural mechanization was consistent with the analysis in Section 4.7.2. However, the pulsing response of ACE to sunlight was negative at first, positive over several periods, and then tended towards zero. This was consistent with the results of the FMOLS and DOLS tests, indicating that sunlight can reduce ACE.
The impulse response of sunlight to ACE was always negative, which was strong in the first few periods, and then decreased. These results indicate that ACE would have a certain weakening effect on sunlight. The impulse response of sunlight duration to agricultural mechanization was positive at the beginning, turned negative after several periods, and then approached zero, indicating that sunlight duration has a promoting effect on agricultural mechanization. Table 17 shows the calculation results of variance decomposition of the VAR model for agricultural carbon emission, sunlight, and agricultural mechanization. ACE was mainly affected by self-impact, followed by agricultural mechanization. This was highly consistent with the results in Section 4.7.2. Agricultural mechanization was mainly affected by its impact and gradually weakened. Then, it was influenced by carbon emissions from agriculture, but the impact was not significant. Sunlight was mainly affected by its impact, followed by ACE. In phase 15, the contribution of ACE variance reached 3.26%.

Discussion
We use cross-sectional correlation tests to verify the correlation between variables. In this article, we tested the stability of a unit root of panel data using the IPS test, ADF test, PP test, and LLC test. The results show that the variable after the first order difference was stable, indicating that each variable was an integrated sequence of the same order and can be used in the PVAR model. In addition, we validated the long-term cointegration relationship between variables using the Kao test. The results show that there was a long-term integration relationship between these two groups of variables.
Then, ARDL-and VAR-based impulse response methods were used to conduct an empirical study on the relationship between variables. The results show that the impulse response function more intuitively reflected the dynamic influence between the variables studied in this paper. Results from FOMLS and DOLS tests confirmed the robustness of long-term results. In addition, the Granger causality test was used to study the causal relationship between variables, and the findings can be summarized as follows: First, there was bidirectional causality between rainfall, sunlight, agricultural mechanization, and agricultural carbon emissions. Rainfall and sunlight were Grange causes of agricultural mechanization. Second, agricultural mechanization and rainfall both increased carbon emissions, while sunlight reduced them. In the short term, rainfall reduced the level of agricultural mechanization, and sunlight increased it.
We can speculate the reasons for the interaction between variables by focusing on the empirical results. The results of this paper show that agricultural mechanization will promote agricultural carbon emissions, which is similar to the results of Liu [50] and Feng [51]. The consumption of fossil fuels brought about by agricultural mechanization will undoubtedly increase agricultural carbon emissions. However, agricultural mechanization will increase agricultural productivity, reduce straw burning, and enhance fertilizer utilization to a certain extent. In the long run, technological innovation will reduce mechanized energy consumption and agricultural carbon emissions [52]. At the same time, we found that rainfall has a positive impact on agricultural mechanization and a negative impact on agricultural carbon emissions. This is close to the findings of Barberis [53] and Liu et al. [54]. This phenomenon should be due to the multiple effects of climate change on agricultural production [55]. Climate change has led to unusual droughts in some areas, and low rainfall means increased irrigation equipment [56], which increases energy consumption and further increases agricultural carbon emissions. At the same time, reduced rainfall will also induce abnormal pests in some areas, resulting in increased use of pesticides, leading to more agricultural carbon emissions [57].

Conclusions and Policy Implications
Global warming has become a major environmental problem facing all countries of the world, one that is linked to the sustainable development of future human life. Moreover, agricultural carbon emissions are significantly connected to climate change. Therefore, we used panel data from 30 provinces (cities) of China between 2000 and 2019, constructing an empirical research framework on the impact of climate change and agricultural mechanization on agricultural carbon emissions, with empirical results. The results can be summarized as follows. Firstly, there was a two-way causal relationship between rainfall, sunlight, agricultural mechanization, and agricultural carbon emissions. Rainfall and sunshine were the main reasons for agricultural mechanization. Second, agricultural mechanization and rainfall will increase carbon emissions, while sunlight will reduce carbon emissions. In the short term, rainfall reduces agricultural mechanization, while sunlight increases it.
On the basis of our findings, we draw the following policy implications: First, with improving resource utilization efficiency and eco-environmental protection at the core, we should increase financial input, promote the renewal and iteration of agricultural machinery equipment, improve agricultural production efficiency, and encourage research and use of renewable energy. At the same time, a large number of low-carbon agricultural technologies, corresponding equipment, and human resources should be introduced to guide and train relevant universities, research institutions, and social organizations. Secondly, each region should optimize the agricultural industrial structure according to local climatic conditions, adjust measures to local conditions, and find its own positioning. Moreover, the application of renewable energy in agriculture should be explored. At the same time, we should build resource-saving and climate-smart agriculture. Furthermore, we should build resourceefficient and climate-smart agriculture. Third, we should strengthen publicity so that farmers and agricultural subjects will be fully aware of the seriousness of climate change. Only in this way can we deeply understand that agriculture is a double-edged sword in climate change and build awareness of low-carbon economic development.
This study is innovative on the basis of predecessors, but there are some limitations: First of all, this paper only selected two climate variables to study, and climate change on agricultural carbon emissions is not systematic. In addition, the direction of agricultural development in China's provinces is different, and the impact of agricultural mechanization on agricultural carbon emissions in different regions is also different. This difference between regions has not been reflected in this article. It is suggested that future research should focus on regionally differentiated and multi-level low-carbon policies and government behaviors.

Conflicts of Interest:
The authors declare no conflict of interest.

Breush-Pagan
Equation (A1) is suitable for small sample numbers and time series, Equations (A2) and (A4) are suitable for large sample and dynamic panel data, and Equation (A3) is suitable for large samples. In Equation (A5), ε ij and ε jt refer to the standard error, and µ ij refers to the correlation coefficient of the residual.
Substitute the above calculation results into the expression of Z ρ and Z t statistics, obtain the value of statistics, then check the critical value for comparison, and finally make a judgment.

Appendix B.3. LLC and Breitung
The LLC test uses the proxy variables that are removed from the correlation and the deterministic effect, and the principle of the LLC test is in the form of the ADF test. The following is the basis for the LLC test: It is also assumed that all individuals in an LLC are homogeneous under the original and alternative assumptions, as follows: H 0 : ρ 1 = ρ 2 = · · · = ρ N = 0,H 1 : ρ 1 = ρ 2 = · · · = ρ N < 0.
Using combined data under the above assumptions to calculateρ and se(ρ), then t ρ =ρ se(ρ) and modified t * ρ = where T = T − P − 1, P = 1 N ∑ N i=1 P i . The Breitung test is similar to the LLC test. Dynamic terms are first removed, ∆y it−j from ∆y it and y it , and then normalized to proxy variables: Finally, ADF regressionε it = ρ ε it−1 + v it is used to inspect unit root.

Appendix B.4. IPS
Previous studies [58] pointed out that the stability of the whole panel data would be affected by the differences of some individuals. IPS testing recognizes overall data stability, but individual data instability, overcoming the drawback that LLC testing is too rigorous and more consistent with economic data.
The IPS test models are as follows: ∆y i,t = α t + ρ i y i,t−1 + ε i,t , i = 1, 2, . . . , N, t = 1, 2, . . . , T The null and alternative hypotheses of the IPS test are as follows: H 0 : ρ i = 0, ∀ i ∈ N H 1 : ρ i < 0, i = 1, 2, . . . , N 1 ρ i = 0, i = N 1 + 1, N 2 + 2, . . . , T In addition, IPS requires that the temporal data of stable individuals be partially non-zero-sum lim The t-statistic of the IPS test is calculated as follows: t ρ i is the ADF statistic of the i-th sequence; the IPS test statistics are as follows: Appendix C Note: *** indicates that this column has 1% visibility.