Next Article in Journal
Direct and Indirect Effects of Urease and Nitrification Inhibitors on N2O-N Losses from Urea Fertilization to Winter Wheat in Southern Germany
Next Article in Special Issue
Achieving Carbon Neutrality for A Future Large Greenhouse Gas Emitter in Quebec, Canada: A Case Study
Previous Article in Journal
Traceable Determination of Atmospheric Mercury Using Iodinated Activated Carbon Traps
Previous Article in Special Issue
The Relationship of Energy and CO2 Emissions with GDP per Capita in Colombia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Energy Balances and Greenhouse Gas Emissions of Agriculture in the Shihezi Oasis of China

1
State Key Laboratory of Grassland Agro-Ecosystems, Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture, Lanzhou 730020, Gansu, China
2
College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou 730000, Gansu, China
3
College of Information and Science Technology, Gansu Agricultural University, Lanzhou 730070, Gansu, China
4
Xinjiang Western Animal Husbandry Limited Company, Shihezi 832000, Xinjiang, China
*
Author to whom correspondence should be addressed.
Atmosphere 2020, 11(8), 781; https://doi.org/10.3390/atmos11080781
Submission received: 1 June 2020 / Revised: 10 July 2020 / Accepted: 12 July 2020 / Published: 24 July 2020
(This article belongs to the Special Issue Greenhouse Gas Emission Mitigation: Feasibility and Economics)

Abstract

:
The objective of this study was to evaluate the difference of crop and livestock products regarding energy balances, greenhouse gas (GHG) emissions, carbon economic efficiency, and water use efficiency using a life cycle assessment (LCA) methodology on farms in three sub-oases within the Shihezi Oasis of China. The three sub-oases were selected within the Gobi Desert, at Shizongchang (SZC), Xiayedi (XYD), and Mosuowan (MSW), to represent the various local oasis types: i. Oasis; ii. overlapping oasis-desert; and iii. Gobi oasis. The results indicated that crop production in XYD Oasis had higher energy balances (221.47 GJ/ha), and a net energy ratio (5.39), than in the other two oases (p < 0.01). The production of 1 kg CW of sheep in XYD Oasis resulted in significantly higher energy balances (18.31 MJ/kg CW), and an energy ratio (2.21), than in the other two oases (p < 0.01). The water use efficiency of crop production in the SZC Oasis was lower than that of the XYD and MSW oases (p < 0.05). Alfalfa production generated the lowest CO2-eq emissions (8.09 Mg CO2-eq/ha. year) and had the highest water use efficiency (45.82 MJ/m3). Alfalfa (1.18 ¥/kg CO2-eq) and maize (1.14 ¥/kg CO2-eq) had a higher carbon economic efficiency than other crops (p < 0.01). The main sources of GHG emissions for crop production were fertilizer and irrigation. The structural equation modelling (SEM) of agricultural systems in the Shihezi Oasis showed that the livestock category significantly influenced the economic income, energy, and carbon balances.

1. Introduction

The typical Shihezi meta-ecosystem of mountains, oases, and desert in northwest China consists of Tianshan Mountain, Shihezi Oasis, and the Gurbantunggut Desert. The areas of the mountain, oasis, and desert in this ecosystem are 2541 km2, 7681 km2, and 10,996 km2, respectively [1]. The Manasi River is the lifeblood of the Shihezi Oasis. It originates from Tianshan Mountain and runs dry in the Gurbantunggut Desert. Agricultural production in the Shihezi Oasis is located in three sub-oases: Shizongchang (SZC), Xiayedi (XYD), and Mosuowan (MSW). More than 95% of people, farm produce, and energy production are concentrated in these sub-oases. At present, many problems such as secondary salinization of cropland, desertification, chemical pollution, and climate change have seriously threatened the sustainability of oasis agriculture through a lack of knowledge and understanding of oasis agricultural management [2]. There exists a close relationship between agricultural production and energy use [3]. Agricultural production requires high human-applied energy inputs, of which a large proportion are imported (i.e., fertilizer, diesel, electricity, pesticide, etc.), the remaining being produced on farms as bio-energy (i.e., seed, manure, and animate energy provided by living plants in nature). Crop yields can be improved (i.e., at economic maximum) by increasing the net energy (outputs/inputs) of crops [4].
Since 1950, the average rate of increase in global temperatures has been 0.17 °C per decade due to agricultural practices [5]. Agriculture has been considered as a major contributor to GHG emissions [6]. Crop yields in developing countries are projected to decrease due to global climate change [7]. As elsewhere, there are enormous challenges of alleviating greenhouse gas (GHG) emissions from agricultural production due to the heterogeneous nature and biophysical complexity of farming systems in China [8].
Characterizing the energy balances and GHG emissions of agricultural practices in the Shihezi Oasis offer key information to mitigate carbon emissions and ensure food security [6]. Energy consumption and GHG emissions involve on-farm and off-farm inputs, which are carbon-based operations and products [8]. In this study, the system boundary and scope of calculating energy balances and GHG emissions only included farm agricultural production practices using the life cycle assessment (LCA) methodology. The objectives of this study were to calculate the differences of energy balances and GHG emissions from various oasis ecotypes (SZC vs. XYD vs. MSW) in the Shihezi Oasis associated with producing per unit of crop and livestock products using the LCA technique. Moreover, finding the difference in energy balance and GHG emissions between six crops and water use efficiency of crop production in the Shihezi Oasis was the other aim of this study.

2. Materials and Methods

2.1. Study Site

Shihezi Oasis, which is divided into the three sub-oases, the SZC, XYD, and MSW oases, is located in the center of the northern foothills of Tianshan Mountain in XinJiang, China (84°58′–86°64′ E, 43°26′–45°20′ N, Figure 1 and Figure S1). Agricultural production in the Shihezi Oasis is divided geographically into three contrasting systems (Figure 1 and Figure S1):
(a)
SZC sub-oasis, with a total area of 37.58 × 103 ha, of which 20.5 × 103 ha is cultivated, sits on an inclined piedmont plain. The Manasi River winds its way through the northeast of the sub-oasis; other surface and groundwater sources are abundant, and in a rural zone of Shihezi City, a freshwater spring overflow is used for flood (border dyke) irrigation;
(b)
XYD sub-oasis, with a total area of 300 × 103 ha, of which 61 × 103 ha is cultivated, lies where the Manasi River downstream oasis system and Gurban Tunggut Desert overlap;
(c)
MSW sub-oasis, with a total area of 137 × 103 ha, of which 44 × 103 ha is cultivated, lies where the Manasi River allows the sub-oasis to cut through the Gurban Tunggut Desert for a distance of 60 km.

2.2. Data Collection

Data were collected from official records, published literature, and farm survey data. The farm surveys were carried out from 2015 to 2016 with a selection of a random 354 farmers in this study (Tables S1 and S2). Inputs of crop production comprised of labor, seed, fertilizer, pesticide, fuel and electricity consumption, plastic film, and the life and working hours of farm machines. Outputs of crop production accounted as crop products including grain, straw, and roots. Inputs of livestock production were comprised of labor, feed (intake, source, and processing), and fuel and electricity consumption. Outputs of livestock production only included livestock products (carcass weight, milk, and wool). To quantify the GHG emissions from livestock production, data from livestock production included categories, livestock numbers, age, feed sources, and feed usage. Structural equation modelling (SEM) can evaluate direct and indirect effects between variables by expected statistical relationships. An SEM analysis program named AMOS 19.0 (IBM Corporation, New York, USA) was used to evaluate the direct and indirect effects of predictor variables on net income, energy balances, and carbon balances. A good fit for that model is 0 ≤ χ2/df ≤ 2 and 0.05 < p ≤ 1 [9].

2.3. Goal and Scope Definition

Based on the ISO standard, the overall goal and scope of this study was to evaluate energy balances and GHG emissions on farms [10]. The specific aim was to quantify the sustainability impacts associated with energy balances and GHG emissions per ha of cropland and kg of livestock products between three sub-oases in Shihezi of China. The target audience included the farm and its stakeholders, consumers, and the public.

2.4. Functional Unit

The functional unit was defined for the purposes of this study, with GJ/ha for carbon balances of crop production and MJ/kg of carcass weight (CW) and milk for carbon balances of livestock. Based on information from the global warming potential for a 100-year period, the data of CH4 and N2O emissions were converted into CO2 equivalents, with 34 for CH4 and 298 for N2O. The units of GHG emissions were kg CO2-eq/ha for cropland, kg CO2-eq/kg for crop products, and kg CO2-eq/kg CW and milk for livestock products.

2.5. Calculation of Energy Balances

The human-applied energy inputs of crop production are calculated as Equation (1).
E I c r o p = C I l × E C l + C I s × E C s + C I f × E C f + C I p × E C p + C I i e × E C i e + C I p m × E C p m + C I d c × E C d c + C I m d × E C m d
where EIcrop is the unit area of total energy inputs for crop production (MJ/ha), CI is the unit area of energy input for crop production, and EC is the energy coefficient (Table 1). In regards to indexing, l, s, f, p, ie, pm, dc, and md are inputs of human labor (h), seed (kg), fertilizers (kg), pesticides (kg), electricity consumption for irrigation (kwh), and plastic films (kg), respectively. The energy outputs for each category of crop including the grain, straw, and roots are calculated as Equation (2).
E O c r o p = C Y g r a i n × E C g r a i n + C Y s t r a w × E C s t r a w + C Y r o o t × E C r o o t
where EOcrop is the unit area of total energy outputs from crop products (MJ/ha) and CY is the unit area of yield. In regards to indexing, grain, straw, and root are grain (kg), straw (kg), and roots (kg) of crop products, respectively.
The energy inputs of livestock production are calculated as Equation (3).
E I l i v e s t o c k = ( j = 1 m ( L I f e e d , j × E C f e e d , j ) + L I d r u g × E C d r u g + L I l a b o r × E C l a b o r + L I e l e c × E C e l e c + L I c o a l × E C c o a l ) / C W l i v e s t o c k
where EIlivestock is the unit carcass (milk) weight of total energy inputs for livestock production (MJ/kg CW and milk), LI is the unit livestock of energy input for livestock production, and CWlivestock is the unit livestock of carcass weight (i.e., the weight of meat with bone except for fur, viscera, head, hooves, and blood). In regards to indexing, feed, drug, labor, elec, coal, j, and m are feed inputs (kg), veterinary drugs (kg), human labor (h), electricity consumption for lighting of livestock housing (kwh), coal consumption for heating of livestock housing in winter (kg), feed classified j, and feed types, respectively.
The energy outputs of livestock products are calculated as Equation (4).
E O l i v e s t o c k = ( L Y c a r c a s s × E C c a r c a s s + L Y m i l k × E C m i l k + L Y w o o l × E C w o o l ) / C W l i v e s t o c k
where EOlivestock is the unit carcass (milk) weight of total energy outputs from livestock products (MJ/kg CW and milk) and LY is the unit livestock of yield. In regards to indexing, carcass, milk, and wool are carcass weight (kg), milk (kg), and wool of livestock (kg), respectively.
The energy indices of balances and ratios are calculated as follows:
E B c r o p & l i v e s t o c k = E O c r o p & l i v e s t o c k E I c r o p & l i v e s t o c k
N E R c r o p & l i v e s t o c k = E O c r o p & l i v e s t o c k E I c r o p & l i v e s t o c k
where EBcrop&livestock represents the energy balances of crops (MJ/ha) or livestock (MJ/CW and milk). NERcrop&livestock represents the net energy ratio (output/input) of crop or livestock production. EOcrop&livestock and EIcrop&livestock represent the same parameters of the above equation, respectively.

2.6. Calculation of GHG Emissions

The GHG emissions for each category of crop production input are calculated as Equation (7).
C E c r o p = C I s × E F s + C I f × E F f + C I p × E F p + C I i e × E F i e + C I p m × E F p m + A I d c × E F d c + A I m d × E F m d + S O I L r e s
where CEcrop is the unit area of GHG emissions (kg CO2-eq/ha) and EF is the emission factor (Table 2). SOILres is the unit area of GHG emissions from soil respiration (kg CO2-eq/ha) using the following, Equation (8) [17].
S O I L r e s = 1.55 × e 0.031 × T × P × S O C ( P + 0.68 ) × ( P + 2.23 )
where T, P, and SOC are the annual mean temperature (°C), the annual rainfall (m), and soil organic carbon between a 0 and 20 cm depth (kg C.m−2), respectively [10].
Carbon stocks of crop production are calculated using Equation (9) [28].
C S c r o p = C S g r a i n + C S s t e m + C S r o o t
where CScrop, CSgain, CSstem, and CSroot represent the unit area of carbon values accumulated in crops (kg CO2-eq/ha), grain (kg CO2-eq/ha), stems (kg CO2-eq/ha), and roots (kg CO2-eq/ha), respectively.
Carbon balances of crop production are calculated as Equation (10).
C B c r o p = C S c r o p C E c r o p
where CBcrop represents the unit area of carbon balances (kg CO2-eq/ha).
The yearly GHG emissions from each livestock species are calculated using Equation (11).
C E l i v e s t o c k = ( T C O 2 f e e d + T C O 2 d r u g + T C O 2 e l e c + T C O 2 c o a l + T C H 4 e n t e r i c + T C H 4 m a n u r e + T N 2 O m a n u r e ) / C W l i v e s t o c k
where CElivestock is the unit carcass (milk) weight of GHG emissions from livestock (kg CO2-eq/kg CW and milk), TCO2 is the emissions of carbon dioxide, TCH4 is the emissions of methane, and TN2O is the emissions of nitrous oxide. In regards to indexing, enteric and manure are ruminant enteric fermentation (kg CO2-eq/kg CW and milk) and manure management (kg CO2-eq/kg CW and milk), respectively.
The carbon stock and carbon balances of livestock production are calculated as the following, Equations (12) and (13) [29].
C S l i v e s t o c k = ( L W × 0.2 ) / C W l i v e s t o c k
C B l i v e s t o c k = C S l i v e s t o c k C E l i v e s t o c k
where CSlivestock is the carbon stock of livestock (kg CO2-eq/kg CW and milk) and LW is the live weight of livestock. CBlivestock is the carbon balance (kg CO2-eq/kg CW and milk).

2.7. Calculation of Carbon Economic Efficiency

Based on emissions per one kilogram of carbon dioxide equivalence from agricultural products, the total carbon economic efficiency associated with the mean market price of these products in 2015 and 2016 (Table S3) is calculated using Equation (14) [28].
C E E = Y P × P R I C E C E
where CEE is the carbon economic efficiency (¥/kg CO2-eq), YP is the yield of agricultural products (kg), PRICE is the price of agricultural products (¥), and CE is the GHG emissions from agricultural production (kg CO2-eq).

2.8. Calculation of Water Use Efficiency

The analysis of energy balance per 1 cubic meter of water can find the relationship between irrigation and the net energy of crop production in the Shihezi Oasis. The water use efficiency of crop production is calculated using Equation (15):
W U E c r o p = E B c r o p W U c r o p
where WUEcrop, EBcrop, and WUcrop represent the water use efficiency of crops grown (MJ/m3), energy balances (MJ/ha), and the water use of crops grown (m3/ha), respectively.

2.9. Statistical Analyses

We used the statistical program named Genstat17.0 (VSNI Corporation, Hemel Hempstead, UK) to analyze differences in energy indices (energy balances and net energy ratio), carbon indices (GHG emissions, carbon stock, carbon balances, and carbon economic efficiency), water use, and water use efficiency using inline linear models.

3. Results

3.1. Yield and Water Use of Crop Production

The results of the yield and water use of crop production in three sub-oases of the Shihezi Oasis are presented in Table 3. The yields (kg DM/ha) of maize and cotton in SZC Oasis were significantly higher than in the other two oases (p < 0.05). The water use for wheat (spring) and maize production in SZC Oasis was higher than in the other two oases (p < 0.05).

3.2. Energy Balances, Net Energy Ratio, and Water Use Efficiency of Agricultural Production

The results of energy balances, the net energy ratio, and water use efficiency in three sub-oases of Shihezi Oasis are presented in Table 4 and Table S4. For crop production, the output energy (308.69 GJ/ha), energy balances (221.47 GJ/ha), and net energy ratio (5.39) in the XYD Oasis were the highest among the three production systems; the corresponding values of MSW Oasis were higher than in the SZC Oasis, respectively (Table 4). However, the input energy (42.58 GJ/ha) in the MSW Oasis was lower than in the other two oases (p < 0.01). The net energy ratio (6.45) and energy balances (309.46 MJ/ha) of maize in the XYD Oasis were significantly higher than that of the SZC and MSW oases (Table 4). Alfalfa (12.71) and maize (6.27) had higher net energy ratios than other crops (p < 0.01) (Table S4). For livestock production, the production of 1 kg CW of sheep in XYD Oasis resulted in significantly higher energy balances than in the other two oases. Similarly, the net energy ratio (2.21) of sheep in the XYD Oasis was significantly higher than that of the SZC and MSW oases (Table 4).

3.3. GHG Emissions and Carbon Economic Efficiency of Agriculture Production

Carbon indices to evaluate agriculture production are presented in Table 5 and Table S4. Carbon balances (−1.12 Mg CO2-eq/ha) of crop production in SZC Oasis were lower than in the other two oases (p < 0.05) (Table 5). GHG emissions (17.72 Mg CO2-eq/ha) from cotton production were significantly higher than for the other five crops (p < 0.01) (Table S4). Alfalfa had higher carbon stocks (23.76 Mg CO2-eq/ha) and carbon balances (15.97 Mg CO2-eq/ha) than other crops (p < 0.01) (Table S4). According to the price of agricultural products (Table S3) and carbon input for farm production, SZC Oasis had a lower carbon economic efficiency (0.34 ¥/kg CO2-eq) of crop production than the other two oases (p < 0.05) (Table 5). The carbon economic efficiency of maize (1.14 ¥/kg CO2-eq) and alfalfa (1.18 ¥/kg CO2-eq) were significantly higher than other crops (p < 0.01) (Table S4).

3.4. Water Use Efficiency Based on Energy

Water use and the water use efficiency of crop production in the Shihezi Oasis are presented in Table 3 and Table 4. Among the three sub-oases in the Shihezi Oasis, the water use of crop production in the MSW Oasis is lower than the corresponding value of the SZC and XYD oases (p < 0.05). However, the water use efficiency of crop production in the SZC Oasis is lower than that of the XYD and MSW oases (p < 0.05). The water use (4590 m3/ha) for alfalfa production in the Shihezi Oasis is the lowest, with the highest water use efficiency (45.82 MJ/m3) (Table S4).

3.5. Contribution of Carbon Emissions

To further illustrate the relationship between GHG emissions and other factors, Figure 2 shows the contribution of GHG emissions from inputs of crop production. The major GHG emission sources within the crop production system are fertilizer and irrigation.

3.6. Effect Analysis of Structural Equation Modelling (SEM)

The effect analysis of SEM is presented in Table S5. The path modes show that the direct and total effects of the livestock breeding structure on predicted variables (farm net income, energy, and carbon balances) are much stronger than other dependent variables (Figure 3a, total effects = 0.647; Figure 3b, total effects = 0.898; Figure 3c, total effects = 1.091; Figure 3d, total effects = 0.980; Figure 3e, total effects = 0.898; and Figure 3f, total effects = 0.1.091). Similarly, the path modes show that the indirect effects of water use efficiency on economy income, energy, and carbon balances are much stronger than those of other variables (Figure 3a, indirect effects = 0.196; Figure 3b, indirect effects = 0.297; and Figure 3c, indirect effects = 0.430).

4. Discussion

4.1. Energy Balances and Net Energy Ratio

The present energy balances of agricultural production in the Shihezi Oasis are comparable to those published elsewhere in the world using the similar methodology of calculating the input and output energy values on farms. For example, the input energy of cotton production (58.80 GJ/ha) in the Shihezi Oasis is much higher than that (31.237 GJ/ha) in Iran [30]. This difference refers to the low energy input of weed controlling and harvesting operations using human labor instead of applying machinery in Iran. Nevertheless, Shihezi’s input energy for maize production (49.71 GJ/ha) is close to that reported in Iran (50.458 GJ/ha) due to similar intensive and high input crop production systems [31]. However, our output energy including grain, straw, and roots of maize production (341.21 GJ/ha) is much higher than that only related to grain estimated in Iran (134.946 GJ/ha) [30]. The net energy ratios of wheat and maize production in the research area are much higher than that (3.95 vs. 2.13, 6.27 vs. 2.67, respectively) in Iran [31,32]. The reason for the lower output energy and net energy ratio of crop production estimated in Iran is that the outputs included only grain. However, the net energy ratio of wheat production in this study is higher than that (3.95 vs. 1.59) in Pakistan due to less wheat yield as output energy [33]. The corresponding value of maize production in the research area is also higher than that (6.87 vs. 5.52) in Turkey due to less crop yields [34].
The flow of energy and matter is the driving force of agricultural development [35]. Agricultural production in the Shihezi Oasis is a high inputs and outputs system. With the rapid development of agriculture, high inorganic energy inputs can satisfy people’s expectations of living standards, but at the cost of environmental sustainability. From data in this research study (cf. Table S1), a more sustainable approach for Shihezi Oasis farmers would be to increase both alfalfa acreage and sheep numbers, in tandem.

4.2. GHG Emissions from Agricultural Systems

Similarly, GHG emissions from agricultural systems in the Shihezi Oasis are comparable to those evaluated in other places. For example, GHG emissions for maize production (12.1 MgCO2-eq/ha) in this study are similar to those reported in Iran (12.9 Mg CO2-eq/ha) due to similar energy inputs of maize production [31]. As the effective value of economic output produced by the unit carbon input, the signals of carbon economic efficiency can explain the benefits of carbon cost to society. The present carbon economic efficiency for wheat production ($0.027/kg CO2-eq) is higher than that ($0.023/kg CO2-eq) in the USA [36]. The most significant reason for this difference is GHG emissions from wheat production, ranging from 0.14–0.38 CO2-eq per kg of wheat produced on farms in the USA. Apart from this paper, there is no other related research on carbon balances in China which is significant enough to suggest adjusting the balance of crops with livestock production. As already stated, in the Shihezi Oasis, high inputs such as fertilizer, mulch, and machinery resulted in high outputs in crop production, but also in high GHG emissions. In addition, if GHG emissions are to be related to climate policy framework in the future, similar to legislation introduced in some other countries, it will be essential to know the impact of those polices on crop and livestock production costs.

4.3. Balance between Livestock and Forage Crops

There exists an imbalance between livestock and forage crop energy inputs vs. outputs in the Shihezi Oasis (Figure 4). Livestock numbers in the Shihezi Oasis are far greater than can be supported by the forage crop supply. For example, the current forage crop supply in XYD Oasis could not meet the feed demand of livestock, which had the largest number of livestock among the three sub-oases in the Shihezi Oasis. Livestock feedstuff prices in the Shihezi Oasis, whether bought or sold, are the same (Table S2), indicating that the cost of growing feedstuff “on farm” is less than the cost of buying it from outside. This is good reason for oasis forage maize and alfalfa acreage to be greatly increased, and thus bringing the demand and supply into balance.
Similarly, the revenue (13,976 ¥/ha) of alfalfa production in 2014 was higher than that (10,275 ¥/ha) of cotton production in the oasis of the Manasi River. The differences of net income per 100 sheep in 2014 were ¥29,492 between proposed crop-livestock production compared to the current practices of Nileke County of the Xinjiang Autonomous Region in China [37]. The revenue per hectare of maize production in Dingxi City of Gansu Province in China was much higher than that (¥14,070 vs. ¥3315) of wheat production in 2018. The net income of integrated crop-livestock production in 2018 was 3.19 times that of intensive crop production [38].

4.4. Uncertainty of Energy Balances and GHG Emissions Assessment

As a basic assessment of energy balances and GHG emissions from agricultural systems in the Shihezi Oasis, there still existed many factors of uncertainty. Firstly, there were differences in data of arable land area, livestock numbers, and inputs of agricultural production within each sub-oasis. Secondly, the statistical method of partial data in China’s literature is different from that in developed countries. Thirdly, the coefficient for calculation of carbon balances and GHG emissions in the Shihezi Oasis were estimated using overseas data reported by Zeng et al. [39], and Cheng et al. [6]. In addition, the other factors such as energy consumption and GHG emissions from plastic mulch were only associated with plastic production. In this study, we only used the Tier 1 method of IPCC 2013 to evaluate GHG emissions from livestock production [11]. Nevertheless, the basic estimation of energy balances and GHG emissions in the Shihezi Oasis may offer fundamental information for the Chinese government to develop long-term agricultural policies on food security, energy conservation, and GHG emissions reduction in northwest China.

5. Conclusions

The models of energy balances, carbon balances, carbon economic efficiency, and water use efficiency developed in this study were used to calculate the differences of energy balances and GHG emissions from the three various local oasis ecotypes in the Shihezi Oasis. The evaluation indicated that crop production in the XYD Oasis had higher energy balances (221.47 GJ/ha) and an energy ratio (5.39) than in the SZC and MSW oases. The sheep production per kg of CW in the XYD Oasis had higher energy balances (18.31 MJ/kg CW) and an energy ratio (2.21) than in the other two oases. The water use efficiency of crop production in the SZC Oasis was lower than that of the other two oases (p < 0.05). These evaluation models in the present study were also used to calculate the differences of energy balances, GHG emissions, and carbon economic efficiency from local main crops in the Shihezi Oasis. We found one hectare of forage crops (i.e., alfalfa and maize) generated fewer emissions than any of the other four crops (wheat, cotton, grapes, and tomatoes). Alfalfa production in the Shihezi Oasis had the lowest water use (4594.0 m3/ha) and highest water use efficiency (45.82 MJ/m3) compared with all other five crops (wheat, maize, cotton, grapes, and tomatoes). The carbon economic efficiency of maize, relative to its market value, was significantly higher than that for each of the other five crops. Fertilizer and irrigation were the two main GHG emission sources from crop production in all three sub-oases together. Analysis of SEM showed that the livestock breeding structure and water use efficiency significantly influenced farm income, and energy and carbon balances in the Shihezi Oasis. The agro-system of Shihezi Oasis was dominated by crop production derived by high energy inputs (e.g., irrigation, fertilizer, and pesticide). The present production model of Shihezi Oasis will bring high risk to the environment and market. However, integrated crop-livestock production will increase energy use efficiency and decrease GHG emissions from the agriculture system in the Shihezi Oasis. A more sustainable approach for the agricultural development of Shihezi Oasis would be to increase both forage crop (i.e., alfalfa and maize) acreage and sheep numbers, in tandem.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4433/11/8/781/s1, Figure S1: Longitudinal section of Shihezi study site to show Mountain-Oasis-Desert coupling ecological system (84°58′–86°64′ E), Table S1: Average data of climate, crop, and livestock of the three sub-oases in the Shihezi Oasis (2015–2016), Table S2: The structured questionnaire of farm survey, Table S3: Average market price of inputs and outputs for agricultural production (2015–2016), Table S4: Energy balances, carbon balances, carbon economic efficiency, water use, and water use efficiency of crop grown in the Shihezi Oasis, Table S5: The standardized direct, indirect and total effects between dependent variables and predict variables.

Author Contributions

F.H. (Fujiang Hou) and Z.Y. conceived and designed the experiments; F.H. (Fuqin Hou) and Z.Y. performed the experiments and farm survey; Z.Y. and F.H. (Fuqin Hou) analyzed the data; Z.Y. wrote the main manuscript text; and F.H. (Fujiang Hou) reviewed the manuscript. 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 (No. 31660347 and 31172249), the National Key Project of Scientific and Technical Supporting Programs (No. 2014CB138706), and Programme for Changjiang Scholars and University Innovation and Research Teams (No. IRT13019).

Acknowledgments

The authors thank Roger Davies (NZCFS) for his support in fruitful discussion.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ren, J.Z.; He, D.H.; Wang, N.; Zhu, X.Y.; Li, Z.Q. Models of coupling agro-grassland systems in desert-oasis region. Acta Prataculturae Sin. 1995, 2, 11–19. (In Chinese) [Google Scholar]
  2. Kramer, K.J.; Moll, H.C.; Nonhebel, S. Total greenhouse gas emissions related to the Dutch crop production system. Agric. Ecosyst. Environ. 1999, 72, 9–16. [Google Scholar] [CrossRef]
  3. Ozkan, B.; Akcaoz, H.; Fert, C. Energy input–output analysis in Turkish agriculture. Renew. Energy 2004, 29, 39–51. [Google Scholar] [CrossRef]
  4. Khoshnevisan, B.; Rafiee, S.; Omid, M.; Mousazadeh, H.; Rajaeifar, M.A. Application of artificial neural networks for prediction of output energy and GHG emissions in potato production in Iran. Agric. Syst. 2014, 123, 120–127. [Google Scholar] [CrossRef]
  5. Lal, R. Carbon emission from farm operations. Environ. Int. 2004, 30, 981–990. [Google Scholar] [CrossRef] [PubMed]
  6. Cheng, K.; Pan, G.; Smith, P.; Luo, T.; Li, L.; Zheng, J.; Zhang, X.; Han, X.; Yan, M. Carbon footprint of China’s crop production—An estimation using agro-statistics data over 1993–2007. Agric. Ecosyst. Environ. 2011, 142, 231–237. [Google Scholar] [CrossRef]
  7. Broad, K.; Agrawala, S. The ethiopia food crisis-uses and limits of climate forecasts. Science 2000, 289, 1693–1694. [Google Scholar] [CrossRef]
  8. Marland, G.; West, T.; Schlamadinger, B.; Canella, L. Managing soil organic carbon in agriculture: The net effect on greenhouse gas emissions. Tellus B 2003, 613–621. [Google Scholar] [CrossRef]
  9. Yan, Z.G.; Li, W.; Yan, T.H.; Chang, S.H.; Hou, F.J. Evaluation of energy balances and greenhouse gas emissions from different agricultural production systems in Minqin Oasis, China. PeerJ 2019, 7, e6890. [Google Scholar] [CrossRef] [Green Version]
  10. International Organization for Standardization. ISO14044: Environmental Management—Life Cycle Assessment—Requirements and Guidelines; ISO: Geneva, Switzerland, 2006; pp. 1–16. [Google Scholar]
  11. Intergovernmental Panel on Climate Change (IPCC). Climate Change 2014: Mitigation of Climate Change. Available online: http://www.ipcc.ch/report/ar5/wg3/ (accessed on 21 July 2020).
  12. Wen, D.; Pimentel, D. Energy flow through an organic agroecosystem in China. Agric. Ecosyst. Environ. 1984, 11, 145–160. [Google Scholar] [CrossRef]
  13. Pimentel, D. Handbook of Energy Utilization in Agriculture; Bender, M., Ed.; CRC Press: Boca Raton, FL, USA, 1980; pp. 155–161. [Google Scholar]
  14. Huang, Z.W.; Yang, D.G.; Li, X.P.; Refu, K.T. Analysis on the energy flow of the farmer households and the characteristics of the eco-economic fractals in the middle and lower reaches of the Tarim river. Arid Zone Res. 2004, 21, 308–312. (In Chinese) [Google Scholar]
  15. Lu, F.B. Energy flow in the agroecosystem of farming-livestock-fruit. Ecol. Agric. Res. 1994, 2, 40–46. (In Chinese) [Google Scholar]
  16. Steve, M. How much carbon dioxide does the human breathe out every day in the world? China Sci. Tech. Panor. Mag. 2005, 6, 176–177. (In Chinese) [Google Scholar]
  17. Chen, S.T.; Huang, Y.; Zou, J.; Shen, Q.; Hu, Z.; Qin, Y.; Chen, H.; Pan, G. Modeling interannual variability of global soil respiration from climate and soil properties. Agric. For. Meteorol. 2010, 150, 590–605. [Google Scholar] [CrossRef]
  18. West, T.O.; Marland, G.A. Synthesis of carbon sequestration, carbon emissions, and net carbon flux in agriculture: Comparing tillage practices in the United States. Agric. Ecosyst. Environ. 2002, 91, 217–232. [Google Scholar] [CrossRef]
  19. Shi, L.G.; Chen, F.; Kong, F.L.; Fan, S.C. The carbon footprint of winter wheat-summer maize cropping pattern on north China plain. China Popul. Resour. Environ. 2011, 21, 93–98. (In Chinese) [Google Scholar]
  20. Blook, H.; Kool, A.; Luske, B.; Ponsioen, T.; Scholten, J. Methodology for Assessing Carbon Footprints of Horticultural Products. Available online: https://www.researchgate.net/publication/285690888 (accessed on 21 July 2020).
  21. Lu, F.; Wang, X.K.; Han, B. Assessment on the availability of nitrogen fertilization in improving carbon sequestration potential of China’s cropland soil. Chin. J. App. Ecol. 2008, 19, 2239–2250. (In Chinese) [Google Scholar]
  22. Dubey, A.; Lal, R. Carbon footprint and sustainability of agricultural production systems in Punjab, India, and Ohio, USA. J. Crop Improv. 2009, 23, 332–350. [Google Scholar] [CrossRef]
  23. Adom, F.; Maes, A.; Workman, C.; Clayton-Nierderman, Z.; Thoma, G.; Shonnard, D. Regional carbon footprint analysis of dairy feeds for milk production in the USA. Int. J. Life Cycle Assess. 2012, 17, 520–534. [Google Scholar] [CrossRef]
  24. Lal, R. Enhancing crop yields in the developing countries through restoration of the soil organic carbon pool in agricultural lands. Land Degrad. Dev. 2010, 17, 197–209. [Google Scholar] [CrossRef]
  25. Wang, J.Q. Agricultural Standards—Feeding Standard of Meat Producting Sheep and Goats (NY/T 816-2004); Ministry of Agriculture of China, Ed.; China Agriculture Press: Beijing, China, 2004; pp. 17–20.
  26. Dyer, J.A.; Desjardins, R.L. Carbon dioxide emissions associated with the manufacturing of tractors and farm machinery in Canada. Biosyst. Eng. 2006, 93, 107–118. [Google Scholar] [CrossRef]
  27. Meng, X.H.; Cheng, G.Q.; Zhang, J.B.; Wang, Y.B.; Zhou, H.C. Analyze on the spatialtemporal characteristics of GHG estimation of livestock’s by life cycle assessment in China. China Environ. Sci. 2014, 34, 2167–2176. (In Chinese) [Google Scholar]
  28. Shi, L.G.; Fan, S.G.; Kong, F.L.; Chen, F. Preliminary study on the carbon efficiency of main crops production in north China plain. Acta Agron. Sin. 2011, 37, 1485–1490. (In Chinese) [Google Scholar] [CrossRef]
  29. Wu, C.C.; Gao, X.Y.; Hou, F.J. Carbon balance of household production system in the transition zone from the loess plateau to the Qinghai-Tibet Plateau, China. Chin. J. Appl. Ecol. 2017, 28, 3341–3350. (In Chinese) [Google Scholar]
  30. Pishgarkomleh, S.H.; Sefeedpari, P.; Ghahderijani, M. Exploring energy consumption and CO2 emission of cotton production in Iran. J. Renew. Sustain. Energy 2012, 4, 427–438. [Google Scholar]
  31. Yousefi, M.; Damghani, A.M.; Khoramivafa, M. Energy consumption, greenhouse gas emissions and assessment of sustainability index in corn agroecosystems of Iran. Sci. Total Environ. 2014, 493, 330–335. [Google Scholar] [CrossRef] [PubMed]
  32. Khoshroo, A. Energy use pattern and greenhouse gas emission of wheat production: A case study in Iran. Agric. Commun. 2014, 2, 9–14. Available online: https://www.researchgate.net/publication/262105641 (accessed on 21 July 2020).
  33. Abbas, A.; Yang, M.L.; Ahmad, R.; Yousaf, K.; Iqbal, T. Energy use efficiency in wheat production, a case study of Paunjab Pakistan. Fresenius Environ. Bull. 2017, 26, 6773–6779. Available online: https://www.prt-parlar.de/download_feb_2017/ (accessed on 21 July 2020).
  34. Baran, M.F.; Gokdogan, O. Comparison of energy use efficiency of different tillage methods on the secondary crop corn silage production. Fresenius Environ. Bull. 2016, 25, 3808–3814. Available online: https://www.prt-parlar.de/download_afs_2016/ (accessed on 21 July 2020).
  35. Sere, C.; Steinfeld, H.; Groenewold, J. Food and Agricultural Organization of the United Nations (FAOUN). Available online: http:/www.fao.org/WAIRDOCS/LEAD/X6101E/X6101E00.HTM (accessed on 21 July 2020).
  36. Soussana, J.F.; Tallec, T.; Blanfort, V. Mitigating the greenhouse gas balance of ruminant production systems through carbon sequestration in grasslands. Animal 2010, 4, 334–350. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Zhang, Y.H. The Research of Grassland Agricultural Development Model in Xinjiang. Ph.D. Thesis, Shihezi University, Xinjiang, China, 2015. (In Chinese). [Google Scholar]
  38. Zhang, P.P. Comparison of Pastoral Agriculture and Traditional Agricultural Production Pattern in Dingxi City. Master’s Thesis, Lanzhou University, Gansu, China, 2019. (In Chinese). [Google Scholar]
  39. Zeng, X.F.; Zhao, S.W.; Li, X.X.; Li, T.; Liu, J. Main crops carbon footprint in pingluo county of the Ningxia Hui autonomous region. Bull. Soil Water Conserv. 2012, 32, 61–65. (In Chinese) [Google Scholar]
Figure 1. Satellite map of the study site in Xinjiang, China.
Figure 1. Satellite map of the study site in Xinjiang, China.
Atmosphere 11 00781 g001
Figure 2. Total GHG emissions from crop production and the contribution of emission sources.
Figure 2. Total GHG emissions from crop production and the contribution of emission sources.
Atmosphere 11 00781 g002
Figure 3. Path models showing the direct and indirect effects of predictor variables on farm net income, energy balances, and carbon balances. The path models with significant correlations are presented as solid lines. The values on solid lines represent standardized regression weights. Interrupted lines indicate no significant correlation between two variables. Black arrows indicate positive effects. For each endogenous variable, the relative amount of explained variance is given. Shading indicates the greatest positive direct effect, indirect effect, and total effect between dependent and independent variables. OtoD: The distance from the oasis to the desert (km); SPD: Soil particle diameter (μm); PS: Planting structure (planting crop type); WUE: Water use efficiency (MJ/m3); BS: Livestock breeding structure (breeding livestock category); ECO: Net income (1000 ¥/ha); EB: Energy balance (GJ/ha); CB: The difference value of carbon stock minus GHG emissions from agricultural production inputs (Mg CO2-eq/ha); and OtoM: The distance from the oasis to the desert (km). χ2: Chi-square. p: Probability level. df: Degrees of freedom. n: Sample size.
Figure 3. Path models showing the direct and indirect effects of predictor variables on farm net income, energy balances, and carbon balances. The path models with significant correlations are presented as solid lines. The values on solid lines represent standardized regression weights. Interrupted lines indicate no significant correlation between two variables. Black arrows indicate positive effects. For each endogenous variable, the relative amount of explained variance is given. Shading indicates the greatest positive direct effect, indirect effect, and total effect between dependent and independent variables. OtoD: The distance from the oasis to the desert (km); SPD: Soil particle diameter (μm); PS: Planting structure (planting crop type); WUE: Water use efficiency (MJ/m3); BS: Livestock breeding structure (breeding livestock category); ECO: Net income (1000 ¥/ha); EB: Energy balance (GJ/ha); CB: The difference value of carbon stock minus GHG emissions from agricultural production inputs (Mg CO2-eq/ha); and OtoM: The distance from the oasis to the desert (km). χ2: Chi-square. p: Probability level. df: Degrees of freedom. n: Sample size.
Atmosphere 11 00781 g003
Figure 4. Feed demand compared to supplemental feed in the research site.
Figure 4. Feed demand compared to supplemental feed in the research site.
Atmosphere 11 00781 g004
Table 1. Factors used for the calculation of energy inputs and energy outputs.
Table 1. Factors used for the calculation of energy inputs and energy outputs.
Energy Factors of Agricultural Production Inputs
Seed
(MJ/kg seed)
Wheat (spring)17.9[11]
Maize104.65[12]
Cotton22.024[13]
Alfalfa108.82[11]
Grape15.16[13]
Tomato16.33[14]
Fertilizer
(MJ/kg fertilizer)
N78.1[12]
P17.4[12]
K13.7[12]
Farmyard manure
(MJ/kg manure)
Animal manure14.63[12]
Pesticide
(MJ/kg pesticide)
Herbicides278[12]
Insecticides233[12]
Fungicides121[12]
Mulch (MJ/kg mulch)Plastic mulch51.9[6]
Fuel (MJ/kg fuel)Diesel47.78[6]
Electricity (MJ/kwh electricity)Electricity
for irrigation
12
Transportation
(MJ/kg truck)
Truck8.8[11]
Maintenance of
machinery (MJ/kg tractor)
Tractor5.21[15]
Human labor (MJ/h)Male0.68[12]
Female0.52[12]
Forage feed
(MJ/kg feed)
Wheat hay15.05[16]
Maize hay15.22[16]
Alfalfa hay18.8[16]
Concentrate feed
(MJ/kg feed)
Maize18.26[16]
Soybean18.83[16]
Wheat husk13.72[16]
Energy Factors of Agricultural Products
Grain (MJ/kg grain)Wheat (spring)12.56[16]
Maize18.26[16]
Cotton22.024[16]
Grape2.341[16]
Tomato1.258[13]
Hay (MJ/kg hay)Wheat (spring)15.05[16]
Maize15.22[16]
Alfalfa18.8[16]
Cotton17.37[16]
Livestock products (MJ/kg product)Lamb12.877[13]
Beef13.88[13]
Milk2.889[13]
Wool23.41[11]
Table 2. Factors used for the calculation of GHG emissions.
Table 2. Factors used for the calculation of GHG emissions.
ItemSub-ItemFactorsReferences
Emission Factors of GHG for Agricultural Production
Seed 1
(kg CO2-eq/kg seed)
Wheat (spring)0.477[18]
Maize3.85[19]
Cotton2.383[18]
Alfalfa9.643[18]
Grape2.35[2]
Tomato1.63[20]
Fertilizer
(kg CO2-eq/kg fertilizer)
N6.38[21]
P0.733[22]
K0.55[22]
Soil emissions CO2 after N application0.633[11]
Soil emissions N2O after N application6.205[23]
Pesticide
(kg CO2-eq/kg pesticide)
Herbicides23.1[5]
Insecticides18.7[5]
Fungicides13.933[24]
Mulch (kg CO2-eq/kg mulch)Plastic mulch18.993[6]
Electricity (tCO2-eq/kwh electricity)Electricity
for irrigation
0.917[25]
Fuel (kg CO2-eq/L fuel)Diesel2.629[6]
Tractor depreciation
(kg CO2-eq/year)
Tractor 781014.07[26]
Tractor 55/600.49[26]
Tractor 1002/12021.32[26]
Tractor 2500.16[26]
Harvester 12000.66[26]
Harvester 1541.34[26]
Feed processing
(kg CO2-eq/kg feed)
Maize0.0102[27]
Soybean0.1013[27]
Wheat0.0319[27]
CH4 emissions from enteric fermentation
(kg CO2-eq/head/year)
Sheep125[11]
Beef cattle1175[11]
Dairy cattle1525[11]
CH4 emissions from
manure management
(kg CO2-eq/head/year)
Sheep2.75[11]
Beef cattle25[11]
Dairy cattle250[11]
N2O emissions from
manure management
kg CO2-eq/head/year)
Sheep62.3[11]
Beef cattle120.4[11]
Dairy cattle106.7[11]
1 Including seed production, cleaning, and packaging; 2 Including exhaled carbon dioxide of adult labor.
Table 3. Yield and water use of crop production in the Shihezi sub-oases.
Table 3. Yield and water use of crop production in the Shihezi sub-oases.
SZC 1XYD 2MSW 3SED 4p-Value
Crop products (Mg DM/ha)
Wheat (spring)15.65 a16.44 a10.67 b0.379<0.05
Maize60.00 a45.45 b48.89 b0.218<0.05
Cotton13.42 a15.68 b15.53 c0.157<0.05
Alfalfa-14.8113.62--
Grape-30.3430.67--
Tomato8175----
Water use (1000 m3/ha)
Wheat (spring)6.77 a5.59 c6.39 b0.176<0.001
Maize9.06 a7.59 b7.623 b0.244<0.05
Cotton8.32 a8.23 a6.30 b0.33<0.05
Alfalfa-4.175.06--
Grape-7.395.96--
Tomato6.56----
Farmland7.67 a7.53 a6.501 b0.202<0.05
1 SZC: Shizongchang; 2 SZC: Xiayedi; 3 SZC: Mosuowan; 4 SED: Standard error of differences; a, b, and c represent means, with different letters in a row differing significantly (p < 0.05).
Table 4. Energy balances, net energy ratio, and water use efficiency of agricultural production in the Shihezi sub-oases.
Table 4. Energy balances, net energy ratio, and water use efficiency of agricultural production in the Shihezi sub-oases.
SZC 3XYD 4MSW 5SED 6p-Value
Energy input
Crop production (GJ/ha)
Wheat (spring)63.19 a60.01 a55.39 b3.239<0.05
Maize56.93 a42.17 c50.03 b1.913<0.001
Cotton76.91 a58.37 b41.12 c5.166<0.001
Alfalfa-19.0923.61--
Grape-35.5837.54--
Tomato49.55----
Farmland69.33 a58.47 b42.58 c4.416<0.001
Livestock production (MJ/kg CW and milk)
Sheep14.96 a10.96 b11.97 b1.225<0.05
Beef cattle38.84 a30.32 b30.01 b5.230<0.05
Dairy cattle1.44 a1.30 b1.32 b0.040<0.05
Energy output (GJ/ha)
Crop production (GJ/ha)
Wheat (spring)255.82 b268.81 a174.36 c14.779<0.001
Maize313.63360.01350.247.0480.317
Cotton99.72 b116.49 a115.41 a2.709<0.05
Alfalfa-278.48255.97-
Grape-101.28101.9--
Tomato204.15----
Farmland230.69 c308.69 a247.64 b16.517<0.001
Livestock production (MJ/kg CW and milk)
Sheep29.2329.2629.290.0100.35
Beef cattle87.0487.0187.020.0020.561
Dairy cattle2.562.542.540.0010.570
Energy balances (GJ/ha)
Crop production (GJ/ha)
Wheat (spring)195.63 b205.8 a118.97 c13.699<0.001
Maize257.07 b309.46 a269.97 b8.519<0.05
Cotton22.82 c58.12 b74.29 a7.598<0.001
Alfalfa-255.39236.36--
Grape-265.7264.36--
Tomato154.6----
Farmland161.12 c221.47 a207.95 b11.452<0.001
Livestock production (MJ/kg CW and Milk)
Sheep17.11 b18.31 a17.18 b1.035<0.05
Beef cattle48.1349.0148.220.9240.125
Dairy cattle1.101.121.080.0510.214
Net energy ratio 1
Crop production
Wheat4.05 a4.27 a3.15 b0.185<0.05
Maize6.09 c6.45 a6.20 b0.225<0.001
Cotton1.30 c2.03 b2.81 a0.219<0.001
Alfalfa-12.0613.05--
Grape-8.478.04--
Tomato4.12----
Farmland3.32 c5.39 a4.63 b0.625<0.001
Livestock production
Sheep1.95 b2.21 a1.97 b0.167<0.05
Beef cattle2.29 b2.43 a2.31 b0.256<0.05
Dairy cattle1.781.801.770.0040.524
Water use efficiency 2
Crop production (MJ/m3)
Wheat (spring)28.92 a,b36.84 a18.61 b3.196<0.05
Maize7.57 c14.42 a11.80 b1.102<0.001
Cotton2.73 c6.99 b11.79 a1.321<0.001
Alfalfa-27.8523.50--
Grape-25.5724.38--
Tomato23.59----
Farmland9.90 b13.99 a14.15 a0.812<0.05
1 net energy ratio = output energy/input energy; 2 water use efficiency: Water use based on energy balances; 3 SZC: Shizongchang; 4 SZC: Xiayedi; 5 SZC: Mosuowan; 6 SED: Standard error of differences; a, b, and c represent means, with different letters in a row differing significantly (p < 0.05).
Table 5. GHG emissions, carbon stocks, carbon balances, and the carbon economic efficiency of agricultural production in the Shihezi sub-oases.
Table 5. GHG emissions, carbon stocks, carbon balances, and the carbon economic efficiency of agricultural production in the Shihezi sub-oases.
SZC 2XYD 3MSW 4SED 5p-Value
GHG emissions
Crop production (Mg CO2-eq/ha. year)
Wheat (spring)8.58 b8.60 b8.99 a0.065<0.05
Maize12.45 a12.12 b12.09 b0.057<0.05
Cotton17.7217.7517.690.1370.981
Alfalfa-8.098.10--
Grape-12.2612.19--
Tomato17.02----
Farmland13.2212.9412.270.3770.534
Livestock production (kg CO2-eq/kg CW and milk)
Sheep9.23 a8.35 a,b7.62 b0.3050.072
Beef cattle22.95 b24.19 a22.87 b0.237<0.05
Dairy cattle0.670.700.710.0130.639
Carbon stocks
Crop production (Mg CO2-eq/ha year)
Wheat (spring)10.44 b10.40 b10.86 a0.0750.069
Maize23.52 a22.84 b22.83 b0.107<0.05
Cotton13.1313.1013.110.1010.995
Alfalfa-23.8323.74--
Grape-10.9811.00--
Tomato10.35----
Farmland12.10 b15.94 a17.33 a0.582<0.001
Livestock production (kg CO2-eq/kg CW and milk)
Sheep1.06 a1.76 a,b1.93 b0.0610.610
Beef cattle3.573.223.540.0970.287
Dairy cattle0.100.110.110.0030.593
Carbon balances
Crop production (Mg CO2-eq/ha year)
Wheat (spring)1.851.801.870.0120.160
Maize11.06 a10.72 b10.75 b0.050<0.05
Cotton−4.59−4.64−4.560.0370.601
Alfalfa-15.7415.63--
Grape-−1.28−1.18--
Tomato−6.68----
Farmland−1.12 b2.99 a5.06 a0.798<0.05
Livestock production (kg CO2-eq/kg CW and milk)
Sheep−7.63 b−6.60 a,b−5.69 a0.330<0.05
Beef cattle−19.39 a−20.97 b−19.32 a0.298<0.05
Dairy cattle−0.53−0.59−0.570.0110.653
Carbon economic efficiency (¥/kg CO2-eq)
Crop production
Wheat0.170.180.170.0010.185
Maize1.161.171.140.0160.802
Cotton0.700.720.700.0150.872
Alfalfa-1.181.17--
Grape-0.410.42--
Tomato0.32----
Farmland0.34 b0.75 a0.82 a0.042<0.001
Livestock production
Sheep0.24 a0.22 a,b0.20 b0.0080.063
Beef cattle0.38 b0.40 a0.38 b0.004<0.05
Dairy cattle0.170.180.180.6090.524
1 net energy ratio = output energy/input energy; 2 SZC: Shizongchang; 3 SZC: Xiayedi; 4 SZC: Mosuowan; 5 SED: Standard error of differences; a, b, and c represent means, with different letters in a row differing significantly (p < 0.05).

Share and Cite

MDPI and ACS Style

Yan, Z.; Hou, F.; Hou, F. Energy Balances and Greenhouse Gas Emissions of Agriculture in the Shihezi Oasis of China. Atmosphere 2020, 11, 781. https://doi.org/10.3390/atmos11080781

AMA Style

Yan Z, Hou F, Hou F. Energy Balances and Greenhouse Gas Emissions of Agriculture in the Shihezi Oasis of China. Atmosphere. 2020; 11(8):781. https://doi.org/10.3390/atmos11080781

Chicago/Turabian Style

Yan, Zhengang, Fuqin Hou, and Fujiang Hou. 2020. "Energy Balances and Greenhouse Gas Emissions of Agriculture in the Shihezi Oasis of China" Atmosphere 11, no. 8: 781. https://doi.org/10.3390/atmos11080781

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