Abstract
It is estimated that farm practices are responsible for as much as 25–30% of greenhouse gas (GHG) emissions. The scale of such emissions from agriculture depends on the model to which that sector conforms. The aim of this study was to investigate how the attributes of particular models of agriculture (i.e. labour, capital, land productivity, machinery, chemical, energy and water inputs) affect GHG emissions at the sectoral level. We test the hypothesis that the productivity of capital and labour as well as the intensity of inputs affect GHG emissions differently in countries with varying levels of economic development. The results indicate that in the poorest countries, a significant barrier to sustainable development is mechanization. In moderately wealthy countries there is a lack of stimuli to implement energy-saving technologies, while in the richest countries it is desirable that mechanization and chemicalization be reduced in favour of energy-saving technologies.
Introduction
Scientists make efforts to conceptualize the term ‘sustainable development’ and to give sustainability a more practical and applicative dimension (Von Wirén-Lehr 2001; Van Cauwenbergh 2007). Definitions of sustainable development usually refer to the need to consider the economic, social and environmental aspects of development. This last element is taking on even greater importance in relation to agriculture, which is the chief component of the natural resources sector (NRS). Agriculture is a sector with particularly strong links to the natural environment, which also means that agriculture faces a broad range of challenges in the context of environmental well-being (Walls, 2006). The difficulty in assessing the impact of agriculture on the environment results primarily from the limited availability of data and the multiplicity of variables which might be used but also from the links between agriculture and other components of the NRS (e.g. some farms are also engaged in fishing or forestry). Environmental impact may be evaluated through the lens of greenhouse gas (GHG) emissions, maintenance of biodiversity or water and soil quality. Nonetheless, the most perceptible global effects currently arise from the first of these factors.
GHG emissions are among the most significant external effects of agricultural production having a global impact on the environment. It is estimated that farm practices are responsible for as much as 25–30% of GHGs, with particularly high emissions in the case of livestock production (Hamuda and Patko, 2010). Agriculture does not emit large quantities of carbon dioxide (CO2), on average 15% of total GHG emission, since agricultural land stores up significant quantities of carbon, leading to a reduction in the presence of CO2 in the atmosphere. The chief gases produced by agriculture are methane (CH4) and nitrous oxide (N2O). These gases have much greater potential to produce a greenhouse effect than CO2. A single kilogram of methane has 25 times the global warming effect of a kilogram of CO2, and a kilogram of nitrous oxide has 298 times the global warming effect of CO2 (Eurostat, 2016; Lesschen et al., 2011). N2O is emitted into the atmosphere from agricultural land chiefly as a result of microbiological processing of nitrogen fertilizers in the soil. Direct emissions of N2O from the soil account for 25% and indirect emissions from the soil represent 16% of GHGs in European Union (EU) agriculture (EEA, 2009). Methane emissions are chiefly a result of the digestive processes of ruminants (mainly cattle and sheep). According to EEA (2009), CH4 accounts for 59% of GHG emissions in EU livestock farming. The global proportions of N2O and CH4 in agriculture are similar (40/60).
The unclear relationships between farm practices and their impact on GHG emissions arise from the large number of indicators which might be used. These may be classified as effect based and means based (van der Werf and Petit, 2002). These efforts, however, are focused on the microeconomic level and most of the existing indicators are of that type. They show how farm practices (such as use of fertilizers, pesticides) at the farm level impact various aspects of the natural environment (Galan et al., 2007; Torrellas et al., 2013). The question is whether the microeconomic perspective is sufficient? Although there have been references to the need for a broader view on the NRS, including the agricultural sector (Payraudeau and van der Werf, 2005), as far as we are aware insufficient attempt has been made to explore the dependence of a sectoral model of agriculture in different parts of the world and the environmental impact produced by that sector. There are reports indicating macroeconomic determinants of GHG emissions, but these refer to the economy as a whole (EEA, 2016). According to the literature, the reduction in GHG emissions in the EU can be attributed to several factors (using the Kaya decomposition identity). Of these, the lower energy intensity of gross domestic product (GDP) plays the largest role, followed by the carbon intensity of energy, which reflects the lower use of very carbon-intensive fuels such as coal, and of fossil fuels in general. These two factors driving emissions down are partially offset by higher population and GDP per capita as well as by a higher share of non-energy-related (non-combustion) emissions, including agricultural production. Studies of this type relating specifically to the agriculture are nonetheless lacking.
The sectoral perspective is however of great importance for the shaping of agricultural policy towards small farms dominating the agricultural sector in many countries. If one considers a small farm as a holding that sells less than 50% of its agricultural output and consumes the remainder within the farm household, then there are over 6 million such units in the EU27 (FAO and OECD, 2012). Hence, the aim of this study was to investigate how sectoral characteristics (models) of agriculture affect GHG emissions across a sample of 40 countries over a long-term period. To obtain robust results, we divided the sample into clusters of countries where conditions for agricultural production at national level (i.e. level of support, GDP per capita, average farm size) were relatively similar.
Methodology
Our approach differs to a significant extent from those adopted in previous studies on the environmental impact of agriculture (Czyżewski and Kryszak, 2017). We evaluate environmental effects based on a scale of GHG emissions from the NRS (agriculture, hunting, forestry and fishing) in individual countries. We study the NRS as a whole and adopt an input–output approach based on a macroeconomic model of the economy, as pioneered by Leontief (1966). In the macroeconomic input–output table, the category ‘agriculture, hunting, forestry and fishing’ is viewed as a ‘proxy’ for agriculture. The source of the data is the World Input–Output Database
1
and the associated Socio Economic Accounts and Environmental Accounts (WIOD, 2016) as well as figures from Eurostat, Food and Agriculture Organization of the United Nations (FAOSTAT) and the World Bank (concerning number of farms, GDP and support for agriculture measured by the nominal rate of assistance (NRA) index). The study covers the period 1995–2009 in 40 countries. We divided our analysis into two stages. In the first, we attempt to address the effect of structural factors that constitute the conditions for agriculture development (Czyżewski and Kryszak, 2017): GDP per capita (IMF, 2016); Average size of a farm in a given country in hectares (source: Lowder et al., 2014); Level of support to agriculture measured by the NRA in 2009 (from the World Bank’s Distortions to Agricultural Incentives database, Anderson, Nelgen 2013).
For this purpose, we performed cluster analysis (by Ward’s method) (Figure 1) using standardized values of the three listed structural variables. The disjointness of the clusters was verified using the Silhouette index S(i) as recommended by Gatnar and Walesiak (2004). The importance of the above-listed structural variables for GHG emissions has been confirmed by EEA (2016). We know that higher wealth of a country and its inhabitants is often correlated with higher per capita emissions (Aldy, 2006), while on the other hand it enables the implementation of modern, more environmental-friendly technologies. The introduction of consumption-based rather than production-based accounting of CO2 emission changes to some extent popular beliefs about GHG emission. It turns out that countries such as China perform better, whereas Western Europe, the United States and Japan have even higher emissions (Davis and Caldeira, 2010). Aldy (2006) discusses an example of environmental Kuznets curve showing that firstly when per capita income grows the emission is also growing but when the level of income is sufficiently high, the emission may decrease. There is also uncertainty as to whether an increase in per capita income in itself leads to an increase in environmental awareness. Assistance in this regard may come from institutional stimuli, including agricultural and environmental policies, which stimulate pro-environmental methods of production. It is clear that richer countries have greater possibilities in this regard, but an increase in GDP per capita is not synonymous with a pro-environmental orientation of agricultural policy. There is also the issue of the agrarian structure, shaped both historically and geographically (i.e. depending on natural conditions), which we express by means of the average farm size.

Clusters of countries according to the structural variables: GDP per capita, average size of farms and support to agriculture (NRA).
In the second stage, we computed panel regression models for each of the clusters A, B and C, beginning with the least wealthy countries, as in equation (1) (variables expressed in US dollars based on fixed prices from 1995):
where,
i denotes a country; t denotes a year; GHG denotes total methane, nitrous oxide and carbon dioxide in tonnes of CO2 equivalent per hectare; W prod denotes labour productivity (own labour compensation, ‘000 US$ per 1 engaged); C prod denotes capital productivity (‘000 US$ per Gross Fixed Capital Formation (GFCF) in ‘000 US$); L prod denotes land productivity (agricultural output in ‘000 US$ per ha); Chem input denotes chemical input (‘000 US$ per ha); Mach input denotes machinery input (‘000 US$ per ha); Elec input denotes electricity, gas and water supply (‘000 US$ per ha); Chem input + Mach input + Elec input = ‘intensive input’; β′ is a vector of dummy variables for countries; λ′ is a vector of dummy variables for years; and u is a random error.
A linear form of the function was chosen, as this achieved the best fit. We believe that the above specification of variables reflects well the differentiation between agriculture models at sectoral level. The variables reflecting productivity (W prod , C prod , L prod) enable the capture of a combination of labour intensity, capital intensity and land fertility. High labour intensity is a distinctive feature of subsistence agriculture and of the organic and conservation models (in the last case, however, in combination with higher capital profitability) (Aliber and Hart, 2009; FAO, 2016; OWF, 2016). High capital productivity is a feature of industrial agriculture and precision farming (Kansas State University, 2016, UCSUSA, 2016). Land productivity depends largely on natural factors (climate, soil quality). We treat intermediate consumption in the areas of mechanization, chemicalization, and consumption of energy and fuels as a ‘proxy’ for the diffusion of technological progress and knowledge in these areas. We assumed that the more modern a solution is used, the more likely it is environmental friendly. We selected these three groups of providers by means of production, because they make the greatest contributions to intermediate consumption in the agriculture (Eurostat use table). GHG emission stands for the dependent variable covering methane, nitrous oxide and carbon dioxide in kilogram CO2 equivalent per hectare (Brander and Davis, 2012) produced by the NRS in a given country. To eliminate the effect of exchange rate fluctuations, all variables are expressed in ‘000 US dollars at 1995 prices. It is assumed that the differences in rates of inflation between countries are reflected in the dollar exchange rates of national currencies.
In the first step, the multicollinearity of the variables from model 1 was evaluated on the basis of variance inflation factors (VIFs). Variables not exceeding VIF = 5 were retained in the model (Chatterjee and Hadi, 2006). Next we computed ordinary least squares (OLS) base models. In case of rejection of the hypothesis of applicability of this approach (based on the Breusch–Pagan test), we computed panel models with fixed (FE) and random effects (RE). In these models, we introduced the variables step by step, checking whether the model was stable and whether the addition of a further variable caused changes in the signs of the other regression coefficients. The effect of the time factor was also shown in the panel models (as the linear trend variable), provided that it was significant in capturing the efforts of different countries to reduce GHG emissions. 2 The final models were computed taking into account the Beck–Katz robust standard errors (PCSE). In the results, we present, for each cluster, a final panel model with robust standard errors. We present the models in order from the poorest to the wealthiest countries.
Results and discussion
Three clusters of countries were identified (Figure 1 and Table 1). Cluster A includes 11 countries (Cyprus, Greece, Japan, South Korea, Malta, Poland, Russia, Slovenia, Taiwan, Hungary and Italy). Most of the countries in this cluster have moderate GDP values per capita purchasing power parity (PPP) (outside Japan and South Korea) (mean 21,774), while the average farm size is the lowest (5.24 ha), and the level of support (NRA) is high (mean 0.41). Cluster B consists of countries (Romania, Indonesia, India, China, Turkey, Latvia, Lithuania, Slovakia, Estonia, Mexico, Bulgaria, Brazil) that are usually relatively poor (mean GDP per capita 10,300), while the average level of support to agriculture is the lowest (mean NRA 0.18). In the countries of this cluster, the average farm size is usually small or moderate (mean, 14.4 ha). There are certain exceptions, however, such as Brazil, where large farms dominate, although its agricultural model is assigned as a whole to cluster B. The largest of the clusters is cluster C. This contains 17 relatively wealthy countries, with a mean GDP per capita (PPP) of 32,983 (Portugal, Spain, Germany, Finland, Ireland, Belgium, the Netherlands, Austria, Luxembourg, the United Kingdom, Sweden, France, Canada, Denmark, the Czech Republic, the United States and Australia). The mean value of support (NRA) is 0.26, and compared with the other clusters this value is relatively uniform between countries (low standard deviation). The average farm size for countries in this cluster is relatively high (see Table 1 for cluster characteristics). In each case, the FE model was found to be appropriate. All the models fit fairly well: the least square dummy variable (LSDV) R 2 value ranges from 0.96 to 0.99, and within R 2 from 0.17 to 0.51. Some variables were excluded due to excessive collinearity. All variables in the three models are statistically significant (p values do not exceed 0.1, except for one variable in cluster C for which p = 0.12).
Descriptive statistics of clusters A, B and C.
GDP: gross domestic product; PPP: purchasing power parity; NRA: nominal rate of assistance.
Source: Own study based on data from the World Bank, the IMF, and Lowder et al (2014).
Cluster B (Table 3) contained, relatively speaking, the poorest countries, with fragmented and highly polarized agriculture (farms of moderate size; Table 1), receiving limited state support. The appropriate model was found to be the panel model with FE. Due to strong collinearity, the variables ‘output per ha’ and ‘capital productivity’ were excluded. The model explains more than 96% of the variation in GHG emissions, taking account of individual country effects (LSDV R 2 = 0.9657) which are constant in time (but vary in space). These are contained in the free term of the model but are expressed in detail by the dummy variables in the OLS model. The within R 2 explains 17% of the intragroup (within) variation attributing variation in GHG emissions to the explanatory variables, which vary in time. Their variation in time seems to have significance. There is evidence for the lack of the pro-ecological evolution of the NRS model, in these countries. This is presented in Table 3, which shows that in 2002–2009 compared with 1995–2001, upward change in GHG emissions was visible.
Model for cluster A.a
GHG: green house gas; LSDV: least square dummy variable.
aFixed effects, using 158 observations. Dependent variable (Y): GHG emissions in kilogram CO2 equivalent per hectare fixed prices in 1995 US$.
Source: Own computations using Gretl software based on the database of WIOD (2016).
Model for cluster B.a
GHG: greenhouse gas; LSDV: least square dummy variable.
aFixed effects, using 180 observations. Dependent variable (Y): GHG emissions in kilogram CO2 equivalent per hectare; fixed prices in 1995 US$.
Source: Own computations using Gretl software based on the database of WIOD (2016).
A surprising finding is that an increase in chemical inputs, in particular, leads to lower GHG emissions. Published data shows that technological progress in fertilizer production, and agro-technical progress relating to fertiliser use, is relatively speaking the fastest and leads to significant reduction in GHGs. The effectiveness of nitrogen application in highly developed countries has increased from 30% in the mid-1980s to an average of 60% at present (EEA, 2009). Many analyses have previously shown that when managed properly, fertilizer use fosters carbon capture. ‘When fertilizers are used properly, they help plants to produce more energy than is consumed during the production, transport and application of fertilizers. They also encourage the conversion of carbon dioxide to biomass through photosynthesis’ (IFA, 2009). The positive impact of agro-technology and fertilizer quality on emissions has been confirmed by many authors (Henault et al., 1998; Kaiser et al., 1998; Zieliński, 2014). Recent advances in fertilizer technology undoubtedly help to reduce N2O emissions. ‘Timed-release fertilizers and fertilizers with nitrification inhibitors provide a gradual supply of nitrogen to the crop, synchronous with plant demand for nitrogen’ (Parton et al., 2011). We should also bear in mind that in cluster B the initial level of fertilizer use was low, and so it is undoubtedly subject to low-baseline effects. In spite of the increase in fertilizer use, the level of soil nitrogen concentration at which N2O emissions increase was not attained. Conversely, in this cluster an increase in mechanization inputs leads to an increase in GHG emissions (by 1192 kg CO2 equivalent per hectare for each US$1000). An increase in own labour compensation also produces higher emissions (in simplified terms, wealthy farmers do more harm to the environment). In highly simplified terms, then, we may adopt the interpretation that wealthier farmers (or fishermen) in these countries invest chiefly in machinery, but this involves investment in older devices which are less environmentally friendly.
Paradoxically, although these are the poorest countries in terms of GDP per capita, the mean unit compensation for own labour is higher than in cluster A (US$19,000 vs. US$16,000 per annum), although it is only half of the value for cluster C. Average use of chemicals per hectare is 50% lower than in cluster A but interestingly is 60% higher than in cluster C. This shows that progress in fertilizer usage is most visible in the wealthiest countries, probably because they use more effective fertilizers in smaller quantities per hectare. In view of the high standard deviations, however, it is necessary to treat the mean values of variables here with caution.
For the countries of cluster A (Table 2) (medium GDP per capita, small farms, fragmented agriculture, medium support), the FE model was again found to be appropriate. All variables in this model were statistically significant at α = 0.05. The model explains more than 98% of the variation in GHG emissions, taking account of the individual effects of countries (LSDV R 2 = 0.9897), which are constant in time (but vary in space). The variables ‘output per ha’ and ‘machinery input’ were excluded due to collinearity (VIF > 5). The within-group R 2 explains 33% of the variation within the cluster. This indicates that the evolution of the NRS model in these countries is faster than in cluster B. This is confirmed by analysis of the time influence (linear trend), given in Table 2, which shows that there is a clear downward trend in GHG emissions throughout the period (151 kg (CO2 equivalent)) per hectare per year). Pro-environmental adaptations of the NRS were most visible in this group of countries, in the whole analysis period. An increase in chemical inputs leads to lower GHG emissions (by 1174 kg (CO2 equivalent)) per hectare for every US$1000. The same interpretation as cluster B may be given, although in this case the diffusion of technological and agro-technical progress seems to be faster (the marginal effect is stronger). One of the reasons is probably the higher level of support given to agriculture.
At the same time, increases in labour compensation and in return on capital lead to higher GHG emissions, which is in line with environmental Kuznets curve theory (Aldy, 2006). We may therefore adopt the assumption, in simplified terms, that wealthier farmers in these countries are investing not in modern fertilizers but in energy-intensive technologies that are less friendly to the environment (probably technologies relating to livestock production, which generates more GHG than crop production and substitutes manure and slurry for mineral fertilizer) (Gołębiewska, 2010). For example, in a research conducted in Poland, the authors showed that in the beginning of the period (1996–2000) Polish farmers invested mainly in machinery and other technical devices. However, in the future they are likely to invest in buildings and their improvement. Both these investments are connected with high energy use.
The variable ‘electricity, gas, water supply’ has a negative effect on GHG emissions in this group. This may be an effect of the large share of livestock production in a whole agricultural sector. In 1995, the value of livestock production in total agricultural production in this cluster (excluding Taiwan) was 37%. In 2009, it increased to 39% (own calculations based on FAOSTAT data). Interestingly, the average return on capital in this group was negative, which supports the hypothesis of investment in obsolescent technologies. The intensity of chemical application is the highest (US$400 vs. US$130/ha/annum in the cluster containing the wealthiest countries). The high standard deviations again indicate the need for caution when considering the mean values of the variables.
In the cluster C countries (Table 4) (high GDP per capita, large farms, homogeneous agriculture, high support), the FE model again proved appropriate. It is noteworthy that in this model the variables ‘own labour compensation’ and ‘energy consumption’ take negative signs. This indicates that it is possible to combine environmental and socioeconomic goals (a win-win situation). Similar conclusions were reached by German et al (2017). Cluster C is dominated by the wealthy countries of the EU. It may be said that the European Agricultural Model, in spite of criticism and claims that its pro-environmental measures are merely ‘for show’ (Erjavec and Erjavec, 2015), has the greatest chance of realizing the concept of sustainable development. Our model explains 99% of the variation in GHG emissions, taking account of the individual effects of countries (LSDV R 2 = 0.9956), which are constant in time (but vary in space). The variables ‘output per ha’ and ‘capital productivity’ were excluded due to collinearity. The within-group R 2 explains 50% of the variation within the cluster. This indicates rapid implementation of technical progress and evolution of the NRS model compared with the other clusters. Nonetheless, the time factor was found to be statistically less significant (p = 0.15) with the lower pro-ecological effect compared to cluster A.
Only in this model does an increase in own labour compensation lead to lower GHG emissions, that is to say, wealthier farmers do less harm to the environment. In simplified terms, we may thus adopt the interpretation that wealthier farmers (or fishermen) in these countries invest in energy-saving, environmental-friendly technologies (OECD, 2001), which is the evidence of the increasing environmental awareness and agro-technical knowledge as well as rapid diffusion of technical progress in that area.
The marginal effect of the variable ‘electricity, gas, water supply’ is exceptionally strong: an increase in spending by US$1000 leads to a decrease in GHG emissions by 2784 kg (CO2 equivalent per hectare). This is the evidence of rapid progress in the field of energy-saving technologies. This view is shared by other authors who point out that a key element in pro-environmental farm practices that reduce GHG emissions is effective irrigation (IFA, 2009). The variables representing machinery and chemical inputs do nonetheless have positive signs, which mean that in this group of countries opportunities for further pro-environmental progress in reducing CO2 emissions in these areas are limited. This result has also been corroborated by other studies (Thirtle et al., 2004). As regards the descriptive statistics, own labour compensation was the highest in this group, while the value of chemical inputs was relatively speaking the lowest.
Table 5 summarizes the direction and relative strength of the influence of particular variables (ceteris paribus) on GHG emissions (an upward-pointing arrow indicates a positive effect). The downward arrows in the ‘intensive input’ (grey shading) indicate problematic areas in agriculture development. An interesting evolutionary path can be observed in this area by analysing in turn the poorest, the moderately wealthy and then the wealthiest countries: Initially it is mechanization that is most problematic for the environment, then consumption of energy (electricity, gas and water), while at a higher level of development the most problematic factors are chemicalization as well as mechanization.
Model for cluster Ca.
GHG: greenhouse gas; LSDV: least square dummy variable.
aFixed effects using 245 observations. Dependent variable (Y): GHG emissions in kilogram CO2 equivalent per hectare fixed prices in 1995 US$.
Source: Own computations using Gretl software based on the database of WIOD (2016).
Comparison of marginal effects for models A, B and C.
GHG: greenhouse gas; GDP: gross domestic product.
Concluding comments
This analysis provides evidence that the conflict between the attainment of environmental and socioeconomic goals in the NRS is one that can be reduced. An example is provided by cluster C, which consists largely of countries of the ‘old’ EU15. Our models indicate a positive correlation in these countries between the level of labour compensation and the attainment of environmental goals in the form of reduced GHG emissions. This provides a partial response to the doubts expressed in the cited report (OECD, 2001). Higher wealth of farmers (fishermen, foresters), evaluated in terms of their income, favours the attainment of environmental goals on condition that they make use of energy-saving technologies. Further research will attempt to determine to what degree this is caused by increased awareness among farmers and by the requirements of agricultural and environmental policy or simply by the possibility of investing in more environmental-friendly technologies.
A second general finding is the indication of problem areas in the pro-ecological development of the NRS in the analysed clusters, which ought in principle to be the focus of agricultural and environmental policy. The aim is to indicate a general direction for action and not to lay down detailed guidelines. Among the poorest countries, a significant barrier to the sustainable development of the NRS is mechanization. In the moderately wealthy countries, there is a lack of market conditions and stimuli for the implementation of energy-saving technologies – such conditions ought therefore to be created. In the wealthiest countries, it is desirable that mechanization and chemicalization be limited in favour of energy-saving technologies.
European agriculture made the greatest progress in reducing quantities of GHG emitted per kilocalorie of agricultural production. In 1991–2009, pollution from agricultural production in Europe per capita fell by 22%, and in spite of the high density of population, this is relatively low on a global scale. This results from changes made to agriculture in the EU and may be linked both to the implementation of pro-environmental policy in agriculture, particularly since the early 1990s when the MacSharry reform realigned the principles of the Common Agricultural Policy, and to technological progress expressed by energy-saving technologies and renewable energy sources (Poczta and Baer-Nawrocka, 2016). However, our models show that in the countries of clusters B and A, the market conditions do not encourage the investment of capital in such solutions of environmental-friendly technologies . Thus, it shall be recalled that support for income will not suffice to achieve sustainable development.
Footnotes
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
