Brought to you by:
Letter The following article is Open access

Aligning agri-environmental subsidies and environmental needs: a comparative analysis between the US and EU

, , , , , , , , , , , and

Published 11 May 2021 © 2021 The Author(s). Published by IOP Publishing Ltd
, , Citation Sofia Biffi et al 2021 Environ. Res. Lett. 16 054067 DOI 10.1088/1748-9326/abfa4e

Download Article PDF
DownloadArticle ePub

You need an eReader or compatible software to experience the benefits of the ePub3 file format.

1748-9326/16/5/054067

Abstract

The global recognition of modern agricultural practices' impact on the environment has fuelled policy responses to ameliorate environmental degradation in agricultural landscapes. In the US and the EU, agri-environmental subsidies (AES) promote widespread adoption of sustainable practices by compensating farmers who voluntarily implement them on working farmland. Previous studies, however, have suggested limitations of their spatial targeting, with funds not allocated towards areas of the greatest environmental need. We analysed AES in the US and EU—specifically through the Environmental Quality Incentives Program (EQIP) and selected measures of the European Agricultural Fund for Rural Development (EAFRD)—to identify if AES are going where they are most needed to achieve environmental goals, using a set of environmental need indicators, socio-economic variables moderating allocation patterns, and contextual variables describing agricultural systems. Using linear mixed models and linear models we explored the associations among AES allocation and these predictors at different scales. We found that higher AES spending was associated with areas of low soil organic carbon and high greenhouse gas emissions both in the US and EU, and nitrogen surplus in the EU. More so than successes, however, clear mismatches of funding and environmental need emerged—AES allocation did not successfully target areas of highest water stress, biodiversity loss, soil erosion, and nutrient runoff. Socio-economic and agricultural context variables may explain some of these mismatches; we show that AES were allocated to areas with higher proportions of female producers in the EU but not in the US, where funds were directed towards areas with less tenant farmers. Moreover, we suggest that the potential for AES to remediate environmental issues may be curtailed by limited participation in intensive agricultural landscapes. These findings can help inform refinements to EQIP and EAFRD allocation mechanisms and identify opportunities for improving future targeting of AES spending.

Export citation and abstract BibTeX RIS

1. Introduction

Global concerns surrounding the negative effects of modern agricultural practices are growing, due to associated biodiversity loss, worsening water and air quality, and increased nutrient loading, soil erosion, and climate change (Laurance et al 2014, Leip et al 2015, UNESCO 2015, Beckmann et al 2019). Maintaining agricultural productivity amidst rapid environmental change requires integrated policy responses, and the widespread adoption of sustainable agricultural practices is a pathway to increase agricultural outputs while reducing environmental impacts (UN 2015, Bennett 2017, OECD 2017, Kremen and Merenlender 2018, Pe'er et al 2020, Weber et al 2020).

In the United States (US) and the European Union (EU)—two of the largest agricultural producers and markets globally—a number of policies seek to ameliorate the negative effects of agricultural production and to support environmentally friendly management. Among these, agri-environmental subsidies (AES) are designed to improve the environmental quality of agricultural landscapes through monetary compensation to farmers. Here, we use AES as an umbrella term covering public schemes devoted to the voluntary uptake and implementation of environmentally sustainable agricultural practices. These are generally known as 'best management practices' in the US, and 'agri-environment schemes' in the EU. Specifically, in the US we refer to the Environmental Quality Incentives Program (EQIP), designed to 'help agricultural producers in a manner that promotes agricultural production and environmental quality as compatible goals' (NRCS 2019). In the EU, we consider targeted measures of the European Agricultural Fund for Rural Development (EAFRD) under Common Agricultural Policy (CAP) Pillar II aimed at 'fostering agricultural competitiveness and ensuring suitable management of natural resources and climate action' (European Parliament 2020).

Overall, the ecological efficacy of many AES has been established; however, previous literature suggests highly context-dependent results (Wallander and Hand 2011, Uthes and Matzdorf 2013, Moxey and White 2014). In both the US and the EU, lack of spatial targeting is thought to challenge AES effectiveness. In the US, Wardropper et al (2015) and Qiu et al (2017) showed spatial incongruence in the allocation of water-related AES and areas of poor water quality. In Germany, Früh-Müller et al (2019) found that AES showed varied effectiveness in matching environmental needs and were not targeting high risk areas of nitrogen imbalance or peatland protection. Similarly, Uthes et al (2010) identified a mismatch between AES targeting and environmental needs vis-à-vis erosion control and grassland extensification. While several studies have assessed the spatial targeting of AES, they typically investigate specific countries or regions, rather than the continental scale. Additionally, studies in the US have more frequently focused on farmer uptake of AES (e.g. Reimer et al 2013, Carlisle 2016).

In this study we thus take a broader perspective, analysing AES in the US (EQIP) and the EU (EAFRD), asking: are AES going where they are most needed to achieve environmental goals? As previous empirical work indicates a potential mismatch of payments and environmental needs, we explore to what extent this might differ in the US and the EU, and expand on previous research by analysing the extent to which key socio-economic factors may moderate the spatial allocation of AES.

The spatial distribution of AES is affected by two primary drivers—firstly, by the two-tier allocation mechanisms which determine the amount of funds directed towards US states and EU member states (MS) and, subsequently, towards finer-scale areas, such as counties in the US and NUTS2 regions in the EU (see appendix A for more details). Secondly, it is driven by farmers' participation, which is motivated by social and economic drivers (Lastra-Bravo et al 2015, Malek et al 2019, Brown et al 2020, Piñeiro et al 2020). Thus, while the overarching goals of AES are environmental, programs may provide additional subsidies to underserved groups to support rural development and social equity. For example, the 2014 US Farm Bill enables farmers from historically underserved groups to access additional resources through EQIP (NRCS 2014b). Similarly, EAFRD promotes social aspects through sub-program themes, providing avenues to address social disadvantage and foster inclusion (EU 2013). Farmer's decisions and abilities to apply for AES are also influenced by a variety of non-environmental factors (Malek et al 2019)—here we include age, gender, and tenancy status based on previous literature (e.g. Barbercheck et al 2014, Giannakis 2014, Adusumilli and Wang 2019).

We identified indicators of environmental need that are supported by the literature (table 1) and are at the core of EQIP and EAFRD goals (table 2). We define 'need' in the context of the intended impacts of AES to enhance environmental and socio-economic sustainability of rural landscapes, as stated in their programmatic goals. Then, we use linear mixed models and linear models to explore the associations between our indicator variables and AES funding. Our consideration of US and EU subsidies, as well as our interdisciplinary approach spanning socio-economic and environmental dimensions, present a comprehensive framework to assess AES targeting and its potential to achieve sustainability goals.

Table 1. Indicators selected for the analysis, including descriptions and justification. For each environmental indicator, the table outlines our definition of need for AES. For socio-economic and contextual factors, the table outlines the relationship to AES distribution.

 DefinitionEU data sourceUS data sourceJustification and interpretation of need
Environmental indicators
Nitrogen (N) and phosphorus (P) balanceMass balance in kg ha−1 between elemental input and output of 140 crops.Global map of nutrient balances (West et al 2014) through earthstat.org (accessed 2019).Excessive fertilisation is associated with soil degradation and water eutrophication (Phoenix et al 2012). We define N or P excesses as indicating a higher need for AES. Deficits also indicate need for sustainable management, but are not expected to occur in the EU/US
Local biodiversity intactnessProportion of terrestrial biodiversity that remains intact after human pressure.PREDICTS database of local terrestrial biodiversity (Newbold et al 2015) accounting for various human-induced pressures. Values were inverted for the purpose of the analysis.Agriculture is largely responsible for biodiversity declines. We contend that areas of biodiversity declines should be prioritized by AES support to minimize further biodiversity loss and advance restoration of degraded habitats.
Soil erosionMean estimated agricultural soil loss in tonnes ha−1 yr−1.RUSLE2015 model of soil erosion by water in 2010 (Panagos et al 2015, EC JRC 2020).Global Soil Erosion map values for 2012 (Borrelli et al 2017, EC JRC 2020).Erosion of fertile topsoil high in organic matter affects crop production, water pollution and flood risk, and degradation of terrestrial and aquatic fauna & flora (Panagos et al 2015, 2016). We associate areas with high erosion rates with higher need for AES.
Soil organic carbon (SOC)Mean estimated stock of SOC in agricultural topsoil (0–30 cm) in tonnes ha−1.LUCAS model of topsoil SOC in 2010 (EC JRC 2020). Values were inverted for the purpose of analysis.Global soil organic carbon map (FAO 2019a). Values were inverted for the purpose of analysis.Well-established indicator of soil health and productivity, and an important component of carbon sink (Franzluebbers 2010). We interpret a low level of SOC as a high need for AES support.
Greenhouse gases (GHG) emissionsMean annual CO2 emissions in tonnes yr−1 ha−1 from agricultural sectors.Emissions Database for Global Atmospheric Research for 2014 (Crippa et al 2019). CO2, CH4 and N2O emissions were converted to total CO2 based on their 100-year Global Warming Potential (Pachauri and Reisinger 2007).Agriculture is a significant source of GHG emissions, which contribute to air pollution and climate change (Poore and Nemecek 2018). We define areas with high GHG emissions as having a high need for AES support.
Agricultural water stress (AWS)Strain placed on renewable water sources by crop production.Global Rivers Threat dataset (Vörösmarty et al 2010).Effects on water security, due to shortages of water for individual, industrial and agricultural use, and on biodiversity, due to the distortions of natural flow patterns, and habitat degradation (Vörösmarty et al 2010). We define high AWS as constituting a high need for AES allocation.
Manure nutrients runoffRunoff Vulnerability Index for manure nitrogen (N) by county and NUTS2 calculated following Kellogg (2000), see appendix C.Derived from N excretion estimates (DEFRA 2019) for cattle, pigs, sheep, goats, and poultry in 2016 (EC 2019) and runoff estimates (Fekete et al 2002).Derived from N excretion estimates (DEFRA 2019) by cattle, pigs, sheep, goats, and poultry (USDA National Agricultural Statistics Service, NASS) and 2012 runoff (USGS 2020).Together with excessive chemical inputs, N surface runoff from animal waste is a threat to water bodies (Smith et al 2007, Mallin et al 2015, Motew et al 2018). We define a high N runoff vulnerability index as indicating a high need for AES allocation.
Socio-economic indicators
Proportion of female farmersProportion of female producers per county and NUTS2 region.Derived from Eurostat (EC 2019) number of farms with female managers in 2013 and the number of farms.2012 Census of Agriculture from USDA NASS.Women may be less likely to receive financial support given historical inequality in access to land, capital, credit, and training (e.g. Luhrs 2016), however, they may also be more likely to adopt sustainable practices than their male counterparts (Barbercheck et al 2014).
Proportion of young farmersProportion of under 44 producers per county and NUTS2 region.Eurostat data for 2013 (EC 2019).2012 Census of Agriculture from USDA NASS.An aspect of farming sustainability is the prospect of a new generation taking over from those retiring, and younger farmers may be more likely to participate in AES programs due to longer planning horizons (Giannakis 2014).
Proportion of rented landProportion of rented agricultural land per county and NUTS2 region.Farm Accountancy Data Network data averaged for 2014–2016. Where only lower spatial resolution (NUTS1) was available, the same value was attributed to corresponding NUTS2 regions.2012 Census of Agriculture from USDA NASS.Farmers operating on their property may be more likely to invest in sustainable practices than tenants (Nyaupane et al 2012, Adusumilli and Wang 2019).
Contextual variables
Agricultural land coverProportion of agricultural land cover per county and NUTS2 region.Corine land cover map from 2012, produced by the European Environment Agency. All agricultural classes (211-244), excluding forestry.Cropland and pastureland layers of the National Land Cover Database 2013 (see Yang et al 2018).Higher proportions of agricultural land cover indicate rural regions, which we expect to be associated with higher AES allocation.
Farm sizeAverage size of a farm (ha) per county and NUTS2 region.Derived from agricultural land cover and number of farm operations (EC 2019) for 2013.Derived from number of farm operations sourced from 2012 Census of Agriculture provided by USDA NASS.As farm size increases, land homogenisation and mechanisation and associated externalities increase, with negative repercussions on the environment. Farm size is highly related to participation in AES, although the relationship is heterogeneous (Ahnström et al 2009, Baumgart-Getz et al 2012, Schroeder et al 2013).
Farm incomeAverage farm income by county (USD) and NUTS2 region (EUR)Derived from number of farm operations in 2013 and average income from agriculture, fishery, and forestry in 2014-2016 from Eurostat (EC 2019).2012 Census of Agriculture by the USDA NASS net cash farm income per operation, which includes cash receipts from farming, including government payments, minus expenses.Proxy of socio-economic stability of farmers, which can influence likelihood of AES adoption due to the related monetary investment (Reimer et al 2013, Daloğlu et al 2014).
Gross domestic product (GDP)GDP per capita by county and NUTS2 region.Eurostat (EC 2019) GDP per capita expressed in Purchasing Power Standard (PPS).Derived from the State Science & Technology Institute (SSTI 2020) Real GDP in USD in 2014 and population per county in 2014 from the U.S. Census Bureau.Reflects patterns of economic activity and urbanisation, with high GDP values associated to high industrialization and low ones associated to rural areas.

Table 2. Relevance of the environmental and socio-economic indicators selected for the analysis within the policy frameworks of the US Environmental Quality Incentives Program (EQIP, NRCS 2019) and the EU European Agricultural Fund for Rural Development (EAFRD, EC 2013). The relevance of each indicator is categorized as D = direct; I = indirect; N = not specified.

 RelevancePolicy details
 EQIPEAFRDEQIPEAFRD
Environmental indicators
Nitrogen (N) and phosphorus (P) balanceIIImplied within the priority for improving soil constituents, such as organic matter, contaminants, and nutrients.Implied within the priority for improving soil management, Article 5.4(c) and Article 4(b).
Local biodiversity intactnessDDSupport for at-risk species habitat conservation including development and improvement of wildlife habitat. 10% of funding is dedicated to wildlife.Restoring, preserving, and enhancing biodiversity is one of the Union set priorities, Article 5.4(a).
Soil erosionDDReduction in soil erosion and sedimentation from unacceptable levels are prioritized on eligible land.Article 5.4(c) emphasizes the importance of soil erosion prevention.
Soil organic carbon (SOC)DISupport for producers implementing conservation practices to improve soil health and increase carbon levels in the soil.Implied as part of the objective related to climate action (Article 4(b)) and the priority for promoting a 'low carbon and climate resilient economy'.
GHG emissions from agricultureDDEmissions reduction is prioritized, as well as energy conservation to save fuel and improve efficiency. GHG capture and storage are referenced in the Background section.Climate action is a Union priority and MS 'should be required to spend a minimum of 30% of the total contribution [...] on climate change mitigation and adaptation as well as environmental issues'.
Agricultural water stressDDReductions of nonpoint source pollution, such as nutrients, sediment, pesticides, or excess salinity is a priority, as well as reduction of surface and ground water contamination.Referred to with reference to the Water Framework Directive. To reduce water stress, 'half of the gain in terms of water efficiency should be translated into a real reduction in water use'.
Manure nutrients runoffDI50% of funds are mandated towards livestock producers to target nonpoint pollution sources. During the years analysed in this study, this target was set to 60%.Manure management is addressed for GHGs reduction and in the Nitrates Directive, with an emphasis on extensification.
Socio-economic indicators
Proportion of female farmersID5% of funds are dedicated to socially disadvantaged producers, however, women are not considered socially disadvantaged albeit being prioritized by other USDA subsidy programs.Gender is not included in the Union priorities, but can be prioritized at the MS level. Support for women is listed under 'Thematic sub-programs'.
Proportion of young farmersID5% of funds are dedicated to beginning producers and ranchers, however, young farmers are not targeted explicitly.Support for the entrance of 'young farmers' (under-40 s at time of application and new entrants to the sector).
Proportion of renting farmersNNNot mentionedNot mentioned

2. Comparability of AES programs

The environmental need indicators in table 2 reflect similarities between the overarching goals of EQIP and EAFRD, which prioritize soil health, water quality, biodiversity, and greenhouse gases (GHG) emissions reduction. Aside from their programmatic objectives, these policies exhibit comparable efforts to balance top-down decision-making processes with local autonomy. EQIP's goals are outlined nationally by the USDA, a federal agency. States, however, can incorporate local priorities, a system strengthened in the 2018 Farm Bill (e.g. State Conservationists can designate 'high-priority' practices eligible for increased payments, NRCS 2019). Similarly, EAFRD priorities are set by the European Commission, but parts of the EAFRD are implemented nationally and sub-regionally based on Community-Led Local Development models (EC 2014, Palmisano et al 2016).

While EQIP and EAFRD have similarities, they also differ in their targets and design. For example, while EQIP focuses on reducing negative environmental externalities (Baylis et al 2008), the re-distributive goals of EAFRD place more emphasis on rural development and provision of public goods (e.g. AES for women, young farmers, and cultural landscapes, table 2). Moreover, they increase contributions to poorer regions by reducing direct payments to larger farms (Johnson et al 2010) and grant greater autonomy to local entities than EQIP (Hanrahan and Zinn 2005, EC 2013).

It is important to note that EQIP and EAFRD are part of broader systems of subsidy programs that may influence trends in AES payments by affecting farmer opportunity costs. These include direct payments (Brady et al 2009, Matthews 2013), crop insurance, and land retirement schemes, such as the conservation reserve program (CRP) in the US (McLean-Meyinsse et al 1994, Hellerstein 2017, appendix B). This study focuses on working farmland AES to enhance comparability between programs.

3. Methods

3.1. Agri-environmental subsidy data

We analyzed average US EQIP spending between 2012 and 2014, as these are the most recent years with complete data available (EWG 2020). For the EU (including the UK), we used European Commission data on EAFRD spending for the years 2014 and 2015 reported at NUTS2 level (EC 2018). These years were chosen to reflect the new funding mechanisms in effect since the 2013 CAP reform, but are close enough to the EQIP spending years to warrant a comparison with the US.

EQIP focuses on a set of over 170 practices, ranging from nutrient and crop residue management to irrigation and grazing. To approximate the EAFRD spending allocated for comparable measures, we corrected total EAFRD spending by the proportion of spending for measures M4 ('Investments in physical assets') and M10 ('Agri-environment-climate') using an additional dataset (EC 2020a). These measures were deemed analogous to those of EQIP as they promote sustainable working land and livestock agricultural practices while largely excluding forestry investments. If spending was reported at lower resolution (often NUTS1 level), we used the same proportion for all underlying NUTS2 regions. Since our analysis focuses on spatial targeting rather than temporal analysis of AES impacts, we calculated AES as the average spending per hectare of agricultural land (table 1) by county in the US ($ ha−1) and by NUTS2 region in the EU ($e$ ha−1), modelling 2835 US counties and 211 EU NUTS2 regions.

3.2. Indicators of environmental need and contextual variables

We selected 14 indicators based on their environmental, socio-economic and policy relevance (table 1). These included: (1) six environmental indicators targeted by AES and capturing environmental need; (2) three socio-economic variables related to how subsidies are allocated and applied for in practice; and (3) four contextual variables describing regional agricultural systems. Table 2 describes their policy relevance within AES frameworks. See appendix D for indicators' maps.

Spatial data were masked to agricultural areas and averaged per region (US county or EU NUTS2). High indicator scores reflect high environmental need across all variables. The direction of soil organic carbon (SOC) and local biodiversity intactness were inverted to achieve this. Socio-economic and contextual indicators were included as drivers of allocation and farmer participation that shape environmental outcomes. Predictors were scaled to zero mean and unit standard deviation.

3.3. Statistical analysis

We conducted two sets of analyses for the continental US (48 states) and most regions of the EU (23 MS) to reflect the two-tier allocation system of the subsidy programs. First, we evaluated AES spending patterns among states ('multi-State models'). Then, we assessed spending patterns within states ('individual-State models') to examine how local spending related to the general patterns observed in the multi-State models. All analyses were conducted in R (R Core Team 2020).

This study covers a limited timespan and historical AES allocation may have influenced the predictors due to changes in agricultural management (Taylor and Morecroft 2009, Meals et al 2010, MacDonald et al 2012). Thus, our findings regarding AES spending are relevant primarily for the years analysed.

3.3.1. Multi-state analysis

We analysed the relationship between AES spending and the 14 predictors with a generalized linear mixed model (GLMM) that included state factors (US state, or EU MS) as a random effect using the nlme package in R (Pinheiro et al 2020). We first ran full GLMMs including all indicators as described above. We used backward variable selection based on AIC to identify the most important variables using the MASS package (Venables and Ripley 2002). We then compared the most important indicators between the US and EU. RGLMM were calculated as summary statistics using the MuMIn package (Bartoń 2020) to summarize the variance explained by fixed effects (marginal R2) and fixed and random effects together (conditional R2).

We assessed the stability of the selected variables using 1000 bootstrap replications (Heinze et al 2018) and tested for spatial autocorrelation in the residuals by running a Monte-Carlo simulation of Moran's I statistic using the spdep package (Bivand et al 2020).

3.3.2. Individual-state analysis

We ran linear regressions for a selected group of US states and EU MS, to exclude potential biases caused by allocation preferences towards certain states at the federal or EU level and compare allocation within different agri-economic systems. Only states with at least 15 counties (US) or NUTS2-regions (EU) were included in this analysis. In the EU, this restricted the analysis to the UK, France, Germany and Italy. In the US, out of all states with more than fifteen counties, we picked the five states with the greatest and with the least absolute distance from the average Pearson's correlation coefficient, as summed across all variables. Thus, we selected states with the most extreme and most representative predictor effect across all variables. For each of these states, we ran linear models with the same predictors as the reduced model from the multi-State analysis.

4. Results

4.1. Spatial patterns of environmental subsidies

In the US, high EQIP spending per hectare was concentrated in the coastal regions, while allocation in the Great Plains region was strikingly lower (figure 1). In the EU, most regions with the highest AES spending per hectare are found in Austria, Netherlands, Germany, Italy, Portugal, the Czech Republic and the UK (figure 1). Total average yearly EAFRD spending was 3.2 times higher than EQIP spending for the years of this study (table 3).

Figure 1.

Figure 1. Average spending of selected agri-environment subsidies in (a) counties in the US and (b) NUTS2 regions in the EU among the years 2012–2014 and 2014–2015, respectively. Uncoloured counties and NUTS2 regions were excluded due to missing data.

Standard image High-resolution image

Table 3. Spending comparison of selected agri-environment subsidy programs between the US and the EU: total average spending across the years of this study and average spending by US county and EU NUTS2 region. US spending is in USD ($), EU spending is in EUR ($e$).

SpendingUSEU
Total 0.74 M2.59 M
By US county and EU NUTS2   
  per capita 19.5 ±0.95.8 ± 0.3
  per agricultural area (ha) 5.9 ± 0.322.9 ± 2.1
  per full-time farm worker 501 ± 13.3564 ± 61.6

4.2. AES allocation patterns in the multi-state models

For the US, N and P balance, soil erosion, manure runoff, proportion of young farmers and farm income did not remain in the final model (figure 2, appendix E). The total variance explained was relatively high (R2 conditional = 0.65), but only partially related to the fixed effects (R2 marginal = 0.11). SOC and GHG emissions were positively associated with receipt of subsidies, meaning areas with low SOC and high GHG emissions received more subsidies. However, agricultural water stress (AWS) and local biodiversity were negatively associated with spending, suggesting a mismatch with these environmental needs. All socio-economic variables and contextual factors were negatively related to receipt of subsidies. Agricultural land cover ratio showed the largest effect size.

Figure 2.

Figure 2. Direction of estimates and effect size of the selected multi-state models and inclusion of each variable in all 1000 bootstrap replications. Crossed variables were not included in the final models. For full results see appendix E. *$P \lt 0.05$;**$P \lt 0.01$; ***$P \lt 0.001$.

Standard image High-resolution image

In the EU, average AWS, biodiversity, manure runoff, farm size, and proportion of rented land did not remain in the final model (figure 2, appendix E). Together, fixed and random effects explained half of the variance (R2 marginal = 0.25, R2 conditional = 0.54). Of the environmental need indicators, GHG emissions, N balance and SOC were positively associated with receipt of subsidies, indicating a match between subsidies and these environmental needs, while soil erosion and P balance showed a negative association, suggesting a mismatch in allocation. Similarly to the US, the proportion of agricultural land cover had the largest effect size, with a negative relationship to AES payments. Of the socio-economic variables and contextual factors, agricultural income, GDP, and proportion of young farmers were negatively associated with receipt of subsidies. The opposite was true for the proportion of female farmers.

The US and EU models both had some degree of spatial autocorrelation in the residuals, with Moran's I of 0.16 and 0.15, respectively. However, the indicators remained relatively constant when different spatial modelling methods were used (see appendix H for results).

See appendix F for the distribution of correlation coefficients between spending and predictors.

4.3. AES allocation patterns in individual-State models

Average model fit in the US was higher for states with highest (mean R2 = 0.57 ± 0.27, see appendix G) than for states with lowest distance from multi-state correlation (mean R2 = 0.19 ± 0.07). GHG emissions were significant and positively associated with subsidies allocation in Arkansas and North Carolina, but negatively associated in Montana. Water stress and biodiversity showed negative associations with most selected states, and were significantly negatively associated with allocation in Oregon and Wyoming, and New Mexico and Oklahoma, respectively. Subsidies allocation in Nebraska, Idaho, and Tennessee showed no significant association with the indicator variables.

For the EU MS models, few indicators showed significant associations to AES spending per area (see appendix G). In the UK, only the proportion of agricultural land cover was significantly negatively associated with subsidies, and the positive association with GHG emissions bordered on significance (P = 0.051). In France, average farm income was significantly positively associated with subsidies allocation. Models for Italy and Germany showed no significant associations.

5. Discussion

Our results suggest some level of inefficacy of AES in targeting areas with the highest environmental needs under current program designs in the US and the EU. While we found agreement between spending and low SOC and high GHG emissions in both the US and the EU, and spending and N surplus in the EU, our findings suggest a mismatch of payments with other environmental variables, including AWS and local biodiversity loss in the US and soil erosion and P balance in the EU. This partially confirms our hypothesis of a mismatch between AES and environmental need when accounting for socio-economic factors that may moderate allocation. These findings support the growing body of literature advocating for more robust spatial targeting of AES, as well as increased consideration of landscape multi-functionality, synergies and trade-offs in AES allocation (Shortle et al 2012, Galler et al 2015, Zasada et al 2018, Pe'er et al 2020, Seppelt et al 2020).

Selected models for the US and the EU retained many of the original environmental and socio-economic indicator variables, confirming that both AES programs rely on similar underlying principles and utilise comparable funding mechanisms (Baylis et al 2008). However, differences were found in the amount of variance explained by the models, with differences among states accounting for more variation in the US than in the EU. This suggests greater variability in allocation formulas among US states than EU MS. This may be explained by the EU's mandate that 25% of the MS budget must be allocated to AES (EU 2005). The uniformity of AES spending in the EU was also confirmed by the lack of within-MS variation in allocation emerging from the individual-State models. These models may have benefitted from finer resolution AES data; however, data for NUTS3 regions is only available from MS-specific databases and was thus not accessible here.

5.1. Successful AES allocation for environmental goals

Our analyses showed success of both EQIP and EAFRD in targeting areas of high GHGs, reflecting the programs' goals regarding agricultural emissions reduction and carbon sequestration. This alignment was also reflected in the models of North Carolina and Arkansas, large producers of livestock, rice, soybeans, feed grains, and cotton—all strong drivers of agricultural GHGs (EPA 2015, FTM 2016, NASS 2017). Cropland and livestock emissions reduction are fundamental for meeting Paris Climate Agreement of mitigation targets (Reisinger and Clark 2018, Rogelj et al 2018), and our results suggest that AES may be successfully targeting areas of greatest need. Although both programmes showed a positive association with SOC, this was only significant in the US. As improved SOC sequestration and agricultural GHG emissions mitigation are linked (Frank et al 2017), this may be a potential win-win scenario from the same AES measures. These results are promising and warrant further investigation to understand the drivers of this successful targeting, which could be applied to other areas of environmental need.

The EU model illustrated EAFRD's success in targeting areas of greatest N surplus. The EU has been addressing issues of N imbalance since the 1990s with the nitrates directive (EU 1991), and our results indicate that the spatial distribution of AES funding may help further reduce N surplus. In contrast, this match was not mirrored by P surplus, which was negatively associated with AES spending. This result warrants further exploration, as the EU lacks a common P management strategy (Ronchi et al 2019).

5.2. Mismatches of AES allocation and environmental need

More so than successes, clear mismatches of funding and environmental needs emerged from our results, highlighting potential misdirection of AES and gaps in environmental targeting of EQIP and EAFRD. Although indicators of soil erosion, biodiversity loss, water stress, and nutrient management are explicitly featured in AES programs' goals (table 2), both the US and the EU showed mismatches between policy priorities and funding allocation. Soil erosion was negatively associated with AES in the EU, and did not feature in the US model, while the opposite was true for local biodiversity.

In the EU, the discrepancy between funding allocation and soil erosion comports with previous studies highlighting how the lack of a common binding strategy has inhibited soil conservation efforts (Turpin et al 2017, Helming et al 2018), a shortcoming that the new Healthy Soil Initiative strives to mitigate (EC 2020b). Moreover, the absence of an association with biodiversity conforms with existing literature highlighting insufficient spatial targeting of conservation measures (Pe'er et al 2014, Batáry et al 2015). In the EU, 0.57% of EAFRD funding is reserved for Natura 2000 sites (Measure 12, Dwyer et al 2016); however, the agricultural matrix surrounding protected areas is crucial in supporting biodiversity (Gonthier et al 2014) and biodiversity conservation is among the supporting ecosystem services financed within Measure M10 ('Agri-environment climate').

We found that subsidy allocation in the US concentrated in areas of reduced water stress need, while the variable was excluded from the EU model. Agriculture is a major cause of freshwater ecosystems degradation (Allan 2004, Poole et al 2013) and water efficiency is becoming increasingly important in light of climate change and groundwater depletion (Marshall et al 2015, Cotterman et al 2018). Thus, this mismatch warrants further policy attention. While it is possible that high water stress and at-risk biodiversity areas are more likely to be enrolled in the CRP program in the US due to its focus on environmentally sensitive lands, or in complementary initiatives like the Conservation Reserve Enhancement Program, this result also highlights the lack of geographical targeting by EQIP since its 2002 reform (Shortle et al 2012, Drevno 2016, Hellerstein 2017). At the state level, water stress emerged as a strong negative predictor in Oregon, indicating a mismatch between EQIP spending and environmental need. Oregon's water supply is threatened by droughts and irrigation demands, with repercussions for agriculture and public health (OEC 2012, Schimpf and Cude 2020). In the EU, the post-2020 CAP reform revised the water exploitation index to relate water stress to availability of renewable water resources, possibly enabling more efficient AES allocation in the future (EEA, 2020).

Nutrient surplus and runoff, especially when tied with water stress and flooding, have serious socio-environmental repercussions that can be ameliorated with sustainable practices (Jones et al 2017, Blanco-Canqui 2018, Roy et al 2021); however, they did not remain in the US model. As the US accounts for large proportions of global livestock production and cropland N surplus (West et al 2014, FAO 2019b), nutrients are known threats to the US water supply (Grant et al 2002, Howarth et al 2002). EQIP subsidizes nutrient management plans and infrastructure development to improve surface water quality, and mandates 50% of total spending to livestock producers (NRCS 2017). Notably, it provides waste management assistance to concentrated animal feeding operations, leading sources of nonpoint water pollution due to chemical inputs and manure runoff from animal feed crops (Burkholder et al 2007, Martin et al 2018). While conservation practices have been shown to reduce nonpoint source loadings locally (Poudel 2016, Liu et al 2018, Sneeringer et al 2018), significant amounts of agricultural land are over-fertilised and need improved nutrient management (Jackson et al 2000, NRCS 2011, 2013, 2014a, Long et al 2018). Moreover, these 'pay-the-polluter' initiatives have been criticised due to insufficient resources, as well as their full reliance on voluntary compliance (Collins 2012, Shortle et al 2012, Shortle and Uetake 2015, Drevno 2016). It should be noted that demand for EQIP funding consistently exceeds allocation (Stubbs 2010), thus, funding limitations may constrain AES' ability to optimize spatial targeting; our results highlight a potential gap in AES targeting which should be further investigated, particularly due to the severe impacts of nutrient surplus and manure runoff associated with US high livestock densities.

5.3. AES allocation and socio-economic variables

In some cases, social context may explain the mismatch between funding and environmental need. In the EU, AES were allocated to areas with higher proportions of female producers, suggesting that EAFRD program goals regarding gender equality and rural development successfully influenced funding allocation, despite, or perhaps in competition with, environmental need. In the US, instead, the proportion of female producers was negatively associated with county spending. This could be due to gender inequality not being a main focus of EQIP per se, but only of overarching USDA policies (NRCS 2019). There may be opportunities for EQIP to learn from EAFRD's focus on social inclusion, especially considering the propensity of women farmers to farm in sustainable, conservation-oriented ways in line with AES program goals (Paul and Fremstad 2016).

Neither program was successful at targeting younger farmers, although only the EU includes additional benefits for young farmers in policy targets (table 2). This supports the findings from an EU audit that despite policy priorities, young farmer participation is declining due to poorly defined interventions and unsatisfactory monitoring systems (ECA 2017). Rented land did not remain in the EU model, and our findings suggest a negative association between tenancy and AES in the US. While Reimer et al (2013) argued that EQIP funding may be preferred over CRP in states with larger proportions of rented land, other studies indicated land ownership as a positive predictor of EQIP participation at local scales (Parker et al 2007, Nyaupane et al 2012, Zhong et al 2016). Rented land is not an explicit target of either AES program; however, previous studies have evidenced lower uptake of sustainable practices by tenant farmers, particularly in conventional systems (Sklenicka et al 2015, Walmsley and Sklenička 2017, Ranjan et al 2019). Local EQIP initiatives and selected EAFRD measures could thus be tailored towards tenants to increase widespread adoption of sustainable practices in conventional landscapes.

Given EAFRD's re-distributive goals, it is somewhat surprising that farm income was not significantly associated with AES allocation in the EU, although a negative association did exist. This could be due to our focus on measures M4 ('Investments in physical assets') and M10 ('Agri-environment-climate'), while excluding other EAFRD funds. Furthermore, there are myriads of drivers which could motivate farmer AES applications that we do not capture here, including personal beliefs, previous experiences and crop prices (McCracken et al 2015, Pavlis et al 2016, Holland et al 2020). Without detailed information on the applications received and approved, we cannot isolate drivers of farmer adoption from those of EQIP and EAFRD's selection criteria when interpreting socio-economic patterns of AES allocation.

5.4. Influence of farming context on AES allocation

Contextual factors related to agricultural production systems may also explain mismatches between AES spending and environmental need, as previous studies have suggested that production system type and structure are a larger driver of subsidy allocation than environmental conditions (Reimer et al 2013, Reimer and Prokopy 2014, Zasada et al 2018). Our findings suggest that the potential for AES to remediate environmental issues may be curtailed within policy frameworks by limited participation from farmers engaged in highly intensive and expansive operations. This may be due to higher utilisation costs of certain conservation practices in intensive farming, or a greater predisposition for these practices in lower intensity areas (see Früh-Müller et al 2019). This was evidenced by significant negative associations with the proportion of agricultural land cover in the US and EU, and with farm size in the US. In the US, AES spending per area was strikingly lower in central regions typically dominated by large farms, while coastal areas had higher spending, despite being less dominated by agriculture. Similarly, in the EU, predominantly agricultural regions such as central Germany and France received lower payments per area.

Our results align with Zasada et al (2018), who also found that smaller farms were more likely to receive higher AES in the EU. However, previous research on the influence of farm size on US conservation practice adoption reports contrasting results, with some finding larger farms more (Baradi 2009) and other less (McLean-Meyinsse et al 1994) likely to adopt conservation practices. Others assert that the influence of farm size varies depending on the management practice and conservation program in question (Soule et al 2000, Lambert et al 2007, Reimer 2015). Future research should further explore whether these farms are less likely to engage with AES—particularly since they are drivers of sustainability issues and experience environmental risks that AES may help mitigate.

6. Conclusions

Our analysis is novel in its approach to compare subsidy systems across continents using an interdisciplinary, human-environment systems perspective. The differences between US and EU AES systems are well documented; however, by developing a consistent framework to assess environmental need, we have identified successes and mismatches in subsidies allocation, and common challenges and opportunities for future policy development and research in both the US and the EU. These findings can help inform refinements to EQIP and EAFRD allocation mechanisms and identify opportunities for improved spatial targeting of AES spending. Furthermore, we identify several socio-economic factors associated with AES allocation that bear further investigation, including the relationship between production system and likelihood of applying for and receiving payments. Finer-scale analyses could assess further indicators of particular interest for the US or the EU, which we did not include for comparability—for example, demographic indicators to account for underserved producers, including racial minorities and farmers with income at or below the national poverty level in the US, and traditional farmers in the EU. Moreover, we did not account for environmental and socio-economic climate change risks in this analysis, a critical area for the sector, and a strong opportunity for future research. Finally, although we investigate spatial targeting of subsidies, we did not analyse the temporal impact of AES. Long-term maintenance of conservation practices is critical to ameliorate the negative externalities of agriculture. For example, changes in soil carbon can take decades to manifest. Additionally, historical AES payments may show varying associations to environmental needs, due to changes in policy frameworks and allocation formulas that have occurred in recent decades.

This research contributes to the growing evidence-base surrounding spatial targeting of AES programs, with implications for farmer engagement and environmental quality. Identifying mismatches in allocation is particularly relevant in light of recent and upcoming reforms of both subsidy programs. While EQIP is increasing its focus on addressing soil health and climate resilience, biodiversity loss, nutrient management, and water stress warrant further attention. Similarly, as the CAP 2020 reform recognises the need for fundamental approaches to sustainable agricultural management, we identified loss of biodiversity, P surplus, and water stress as EAFRD sustainability goals that may need additional targeting.

Acknowledgments

This paper is the outcome of the ARAGOG ('An area of conflict: the effects of agricultural subsidization on conservation goals worldwide') programme. This work was supported by the National Socio-Environmental Synthesis Center (SESYNC) under funding received from the National Science Foundation DBI-1639145 and by the Helmholtz Research School for Ecosystem Services under Changing Land Use and Climate (ESCALATE, VH-KO-613). GZ acknowledges the support of the BESTMAP project funded by the European Union's Horizon 2020 research and innovation programme under Grant Agreement No. 817501. DES was supported by NRT-INFEWS: UMD Global STEWARDS (STEM Training at the Nexus of Energy, WAter Reuse and FooD Systems), NSF Grant Number 1828910. The authors thank the Environmental Working Group for generously sharing the EQIP data for the US. In addition, the authors thank Bridget Kerner, Zora van Leeuwen, and Tim Winter for discussion during the early stages of the research, Jonathan Kramer for his leadership, and two anonymous reviewers for helpful comments on previous versions of this paper.

Data availability statement

The data that support the findings of this study are available upon reasonable request from the authors.

Contribution

All authors conceptualized the study and collaboratively defined the methodology. S B led analysis and data curation with support from L E, R T, B C, K T. S B, L E, K T, B C led code development process. S B, R T, B C, E S led manuscript writing, with support from L E, N M, M O, D E S, H W and review from M B, R S, G Z. M B, N M, R S, G Z, S B, R T provided project administration support including co-organizing workshops to advance the study. M B, N M, R S, G Z acquired funding to support the work.

Appendix A.: Brief overview of funding allocation mechanisms

A.1. EQIP

Allocation is administered by the USDA's Natural Resource Conservation Service (NRCS), distributed from the federal government to each state first, and then distributed within each state to its counties. For the national allocation to states, NRCS determines each state's EQIP funds using a formula that reflects national priorities and available natural resources, including: the significance of environmental and natural resource concerns and the opportunity for environmental improvement, the ways the program can best assist producers in complying with Federal, State, local, and Tribal environmental laws, and the amount of agricultural land in different land use categories. For state-level fund distribution, the State Conservationist develops an allocation formula, considering State and local level resource concerns, science-based information on environmental status, and relevant local programs and specialized farming operations (e.g. specialty crops, livestock, organic, small-scale), among other factors (NRCS 2019). The incentives provided include technical assistance and cost-shares of up to 75% of implementation costs.

A.2. EAFRD

Funds are provided in the form of less favoured area payments, agri-environment schemes, and investment support towards rural development in MS. The framework for fund allocation is based on a two-tier process: at the central stage (EU level), the overall framework is established, financial modalities are outlined and eligibility criteria are defined. Additionally, the EU outlines a set of six priority areas, including fostering knowledge transfer and innovation; enhancing viability and competitiveness of agriculture; promoting food chain organization, animal welfare and risk management; promoting resource efficiency, and low-carbon and climate resilient agriculture; preserving and enhancing ecosystems; and promoting social inclusion, poverty reduction, and economic development in rural areas. The distribution of the overall amount for rural development between MS is based on objective criteria and past performance. At the second level, MS develop national strategic plans and set quantitative objectives for priority areas. At least four of these national level priorities must address the EU level priorities. A minimum value of 30% of rural development funds must be set aside for environmental management measures, falling primarily within the measures considered in this study: M4 'Investments in physical assets', which receives the highest share of EAFRD budget (24%), and M10 ('Agri-Environment Climate'), which receives the third highest share (20%, ENRD CP 2015). There are several AESs across different EU MS, and the allocation mechanisms for specific subsidies tend to vary depending on the focus of the scheme (EC 2013).

Appendix B.: Brief overview of the conservation reserve program

Administered by the USDA's Farm Service Agency (FSA), the CRP is a land retirement program. In exchange for a yearly rental payment for 10–15 years, farmers enrolled in the program agree to remove environmentally sensitive land from agricultural production and plant species that will improve environmental health and quality (Hellerstein 2017). The long-term goal of the program is to re-establish valuable land cover to help improve water quality, prevent soil erosion, and reduce loss of wildlife habitat. The main mechanism of farmer enrolment in this program is through a competitive process known as CRP General Sign-up. During a bidding period, any farmer with highly erodible or environmentally sensitive cropland can apply for the program by indicating the parcels they wish to enrol and the annual payments they require together with the contract length. FSA then determines allocation based on national rankings of an Environmental Benefits Index (EBI) score and based on an overall budget that varies year by year. All parcels with an EBI score above the critical national cutoff are accepted while all parcels with an EBI score below the cutoff are rejected (Hellerstein 2017). Because CRP is largely a land retirement program, and because of its distinctly bottom-up funding allocation mechanism and decision-making process, this program is fundamentally different from the focus of our study—conservation practices used on working farmland supported by EQIP and EAFRD with largely top-down allocation strategies.

Appendix C.: Manure nitrogen runoff vulnerability index

For each county, the index was calculated following Kellogg (2000):

Surface runoff is a measure of the potential for water soluble nutrients to run off fields during precipitations and flooding. Thus, the index estimates the maximum potential for livestock manure nitrogen to move from farms to the water supply. We estimated the amount of nitrogen excreted by each livestock type following the nitrate vulnerable zones guidance (DEFRA 2013, 2019).

Livestock typeTotal N produced by 1 livestock unit (kg N year−1)
Adult bovine101
Adult swine88
Sheep7.6
Goat15
Poultry231

Appendix D.: Maps of predictors used in the analysis

Figure D1.

Figure D1. Maps of the predictors used in this study for 48 US states.

Standard image High-resolution image
Figure D2.

Figure D2. Maps of the predictors used in this study for 23 EU member states and the UK.

Standard image High-resolution image

Appendix E.: Multi-state models results

Table E1. Regression coefficients of the selected multi-State models in the US and EU and inclusion of each variable in all 1000 bootstrap replications and the frequency of the same direction of the effect in the selected model and the bootstrap replications. *$P \lt 0.05$; **$P \lt 0.01$; ***$P \lt 0.001$.

 USEU
IndicatorEstimateSESig.InclusionDirectionEstimateSESig.InclusionDirection
Intercept1.270.15***100.0100.02.570.12***100.0100.0
Environmental          
Nitrogen balance     0.180.08*66.699.8
Phosphorous balance     −0.200.09*65.098.8
Local biodiversity−0.130.04**92.4100.0     
Erosion     −0.140.07*67.798.4
Soil organic carbon0.170.03***99.9100.00.130.08ns49.998.8
GHGs emissions0.230.02***100.0100.00.340.08**99.7100.0
Water stress−0.140.04**98.0100.0     
Manure N runoff          
Socio-economic          
Prop. female farmers−0.100.02***99.7100.00.210.07**94.6100.0
Prop. young farmers     −0.130.09ns41.196.8
Prop. renting farmers−0.070.02**86.6100.0     
Contextual           
Farm size−0.110.02***100.0100.0     
Farm income     −0.080.05ns88.157.9
GDP−0.060.02***99.3100.0−0.100.06ns75.6100.0
Prop. agricultural land−0.330.04***100.0100.0−0.450.09***99.4100.0

Appendix F.: Distribution of correlation coefficients with spending

Figure F3.

Figure F3. Violin plots of Pearson correlation coefficients between (A) EQIP spending per hectare of agricultural area in the US and (B) EAFRD spending for selected measures per hectare of agricultural area in the EU and the indicator variables. Only correlation coefficients for states and MS with ≥4 regions are shown.

Standard image High-resolution image

Appendix G.: Individual-state models results

Results of the individual-state linear regressions for the selected US states and EU MS. Figure 4 compares direction and estimate size across the US and EU individual-state models. Table 2 shows the estimate model results the ten selected US states, while table 3 shows the estimate model results for the four selected EU MS.

Figure G4.

Figure G4. Direction of estimates and effect size in the linear models of individual-state AES allocation in (a) the US states with the greatest and shortest distance from the federal average and (b) the EU for four member states with at least 15 NUTS2-regions. Shaded bars indicate statistical significance with α = 0.05.

Standard image High-resolution image

Table G2. Estimate model results for individual-state analysis in the US for the five states with the greatest and least absolute distance from the average Pearson's correlation coefficient summed across all model variables. *$P \lt 0.05$: **$P \lt 0.01$: ***$P \lt 0.001$.

StateVariableEstimateSE t-valueSig.
Greatest distance
New Mexico (n = 29) R2 = 0.85Intercept5.712.202.59*
 Local biodiversity−6.681.34−5.00***
 Soil organic carbon9.632.433.96***
 GHG emissions0.900.980.92ns
 Water stress0.970.651.49ns
 Prop. female farmers−0.030.13−0.21ns
 Prop. renting farmers0.430.162.69*
 Farm size0.100.100.94ns
 GDP per capita−0.190.12−1.58ns
 Prop. agricultural land−0.550.19−2.91**
Montana (n = 52) R2 = 0.56Intercept−1.090.99−1.10ns
 Local biodiversity−0.580.89−0.65ns
 Soil organic carbon1.020.971.05ns
 GHG emissions−2.850.98−2.90**
 Water stress−0.490.49−1.00ns
 Prop. female farmers0.160.220.74ns
 Prop. renting farmers0.620.222.78**
 Farm size−0.150.09−1.77ns
 GDP per capita−0.120.14−0.85ns
 Prop. agricultural land0.470.301.60ns
Nebraska (n = 92) R2 = 0.17Intercept2.120.703.01**
 Local biodiversity−0.890.80−1.11ns
 Soil organic carbon0.200.670.30ns
 GHG emissions0.220.161.42ns
 Water stress−0.600.36−1.65ns
 Prop. female farmers−0.150.28−0.54ns
 Prop. renting farmers0.050.390.13ns
 Farm size−0.020.24−0.10ns
 GDP per capita0.300.221.36ns
 Prop. agricultural land−0.350.53−0.65ns
Oregon (n = 33) R2 = 0.49Intercept1.853.860.48ns
 Local biodiversity−7.094.39−1.61ns
 Soil organic carbon−1.652.78−0.59ns
 GHG emissions7.264.371.66ns
 Water stress−14.685.71−2.57*
 Prop. female farmers−0.901.28−0.71ns
 Prop. renting farmers0.610.940.65ns
 Farm size−1.361.01−1.34ns
 GDP per capita0.500.361.38ns
 Prop. agricultural land−0.141.34−0.10ns
Wyoming (n = 21) R2 = 0.61Intercept−4.993.17−1.57ns
 Local biodiversity−2.443.94−0.62ns
 Soil organic carbon6.595.231.26ns
 GHG emissions3.693.691.00ns
 Water stress−3.311.42−2.33*
 Prop. female farmers−0.610.50−1.20ns
 Prop. renting farmers−1.040.41−2.53*
 Farm size−0.280.16−1.73ns
 GDP per capita0.030.200.14ns
 Prop. agricultural land−0.950.73−1.31ns
Least distance
Arkansas (n = 75) R2 = 0.14Intercept0.854.680.18ns
 Local biodiversity−1.063.80−0.28ns
 Soil organic carbon2.227.500.30ns
 GHG emissions2.641.242.12*
 Water stress−3.174.27−0.74ns
 Prop. female farmers−1.611.93−0.83ns
 Prop. renting farmers−1.351.74−0.78ns
 Farm size−4.135.25−0.79ns
 GDP per capita−1.021.28−0.80ns
 Prop. Agricultural land−2.224.62−0.48ns
Idaho (n = 40) R2 = 0.14Intercept4.541.722.65*
 Local biodiversity−5.684.44−1.28ns
 Soil organic carbon0.361.940.19ns
 GHG emissions−0.320.76−0.41ns
 Water stress−4.434.50−0.98ns
 Prop. female farmers−0.761.02−0.74ns
 Prop. renting farmers−0.130.78−0.17ns
 Farm size−0.021.82−0.01ns
 GDP per capita−1.411.03−1.38ns
 Prop. Agricultural land−1.851.68−1.10ns
North Carolina (n = 92) R2 = 0.13Intercept−0.103.94−0.03ns
 Local biodiversity1.513.270.46ns
 Soil organic carbon0.551.060.52ns
 GHG emissions1.880.623.03**
 Water stress−1.353.09−0.44ns
 Prop. female farmers0.430.890.48ns
 Prop. renting farmers0.911.010.90ns
 Farm size−11.206.38−1.76ns
 GDP per capita−2.210.97−2.27*
 Prop. Agricultural land−2.673.35−0.80ns
Oklahoma (n = 77) R2 = 0.20Intercept1.430.512.80**
 Local biodiversity−1.920.66−2.90**
 Soil organic carbon2.110.802.64*
 GHG emissions0.660.391.67ns
 Water stress−0.650.40−1.63ns
 Prop. female farmers−0.240.34−0.70ns
 Prop. renting farmers0.500.321.57ns
 Farm size0.180.710.26ns
 GDP per capita−0.230.19−1.21ns
 Prop. Agricultural land−0.140.40−0.34ns
Tennessee (n = 93) R2 = 0.15Intercept2.131.731.23ns
 Local biodiversity−1.611.54−1.05ns
 Soil organic carbon2.672.291.17ns
 GHG emissions1.981.201.65ns
 Water stress−0.941.95−0.48ns
 Prop. female farmers−1.700.97−1.75ns
 Prop. renting farmers−1.280.96−1.33ns
 Farm size3.034.000.76ns
 GDP per capita−0.140.64−0.22ns
 Prop. Agricultural land−1.421.70−0.84ns

Table G3. Estimate model results for individual-State analysis in the EU for the member states with more than 15 NUTS2. *$P \lt 0.05$: **$P \lt 0.01$:***$P \lt 0.001$.

Member stateVariableEstimateSE t-valueSig.
United Kingdom (n = 31) R2 = 0.81Intercept74.2431.632.35*
 Phosphorous balance−28.1820.96−1.35ns
 Nitrogen balance−16.8716.11−1.05ns
 Soil erosion8.3412.290.68ns
 Soil organic carbon14.379.811.47ns
 GHG emissions40.3919.422.08ns
 Prop. of female farmers−8.6417.63−0.49ns
 Prop. of young farmers−20.8519.21−1.09ns
 Farm income−5.2413.75−0.38ns
 GDP per capita−0.099.46−0.01ns
 Prop. agricultural land−53.209.58−5.55***
France (n = 21) R2 = 0.74Intercept6.942.672.60*
 Phosphorous balance4.302.241.92ns
 Nitrogen balance−1.301.40−0.93ns
 Soil erosion−0.231.86−0.13ns
 Soil organic carbon0.192.080.09ns
 GHG emissions−0.642.60−0.25ns
 Prop. of female farmers−1.901.82−1.05ns
 Prop. of young farmers0.072.000.04ns
 Farm income13.554.293.16*
 GDP per capita−1.661.32−1.26ns
 Prop. agricultural land3.602.471.46ns
Italy (n = 19) R2 = 0.64Intercept87.53101.060.87ns
 Phosphorous balance−107.2187.54−1.23ns
 Nitrogen balance84.4683.711.01ns
 Soil erosion−0.0310.130.00ns
 Soil organic carbon87.56105.450.83ns
 GHG emissions46.5036.991.26ns
 Prop. of female farmers46.1236.301.27ns
 Prop. of young farmers52.1828.671.82ns
 Farm income28.36185.640.15ns
 GDP per capita−6.0630.49−0.20ns
 Prop. agricultural land23.3336.640.64ns
Germany (n = 34) R2 = 0.42Intercept23.9715.051.59ns
 Phosphorous balance6.335.571.14ns
 Nitrogen balance−5.215.60−0.93ns
 Soil erosion14.9617.920.84ns
 Soil organic carbon4.416.450.68ns
 GHG emissions−1.048.21−0.13ns
 Prop. of female farmers−4.1612.00−0.35ns
 Prop. of young farmers−2.605.91−0.44ns
 Farm income7.606.461.18ns
 GDP per capita−5.984.64−1.29ns
 Prop. agricultural land6.799.540.71ns

Appendix H.: Spatial models

We used Moran's I statistic to determine the relationship between the US and EU model residuals and their surrounding values. We tested this with a Monte-Carlo simulation with 1000 permutations using moran.mc in the spdep R package. This compares the observed value of Moran's I with a simulated distribution to assess the likelihood that the observed values could be observed at random. We tested several spatial models on our data, and present the spatial error models obtained with errorsarlm from the spatialreg package in R.

The US model had a Moran's I of 0.16 indicating that there was some spatial autocorrelation in the residuals. There was a correlation of each county and the adjacent counties of 0.31, and the random intercepts for each state revealed a pattern of coastal states having positive random intercepts. Figure H5 shows the residuals of the model, and figure H6 the random intercepts. Table H4 and figure H7 present the spatial error model.

Figure H5.

Figure H5. Residuals of the US model by county.

Standard image High-resolution image
Figure H6.

Figure H6. Random intercepts of the US model by state.

Standard image High-resolution image
Figure H7.

Figure H7. Residuals of the spatial error model by county.

Standard image High-resolution image

Table H4. Spatial error model for the US. *$P \lt 0.05$;**$P \lt 0.01$; ***$P \lt 0.001$.

VariableEstimateSE z-valueSig.
Intercept1.1230.147.993***
Local biodiversity−0.0970.05−2.004*
Soil organic carbon0.1640.044.423***
GHG emissions0.2530.039.636***
Water stress−0.1070.05−2.344*
Prop. female farmers−0.0690.02−2.928**
Renting farmers−0.0680.02−2.878**
Farm size−0.1110.03−4.367***
GDP per capita−0.060.02−3.678***
Prop. agricultural land−0.3180.04−7.906***

The EU model also had a Moran's I of 0.15. The random intercepts did not show a clear geographic pattern, however, there was a correlation of each NUTS2 region and the adjacent regions of 0.26. Figure H8 shows the residuals at NUTS2 level, figure H9 the random intercept for each country, and table H5 the spatial error model output shown in figure H10.

Figure H8.

Figure H8. Residuals of the EU model by NUTS2.

Standard image High-resolution image
Figure H9.

Figure H9. Random intercepts of the EU model by member State.

Standard image High-resolution image
Figure H10.

Figure H10. Residuals of the EU spatial error model by NUTS2.

Standard image High-resolution image

Table H5. Spatial error model for the EU. $*P \lt 0.05$;**$P \lt 0.01$; ***$P \lt 0.001$.

VariableEstimateSE z-value Sig.
Intercept2.7510.3398.113***
Nitrogen balance0.1660.0812.062*
Phosphorous balance−0.1730.084−2.051*
Soil erosion−0.1420.062−2.302*
Soil organic carbon0.1970.0822.401*
GHG emissions0.3130.0774.046***
Prop. of female farmers0.2580.0843.056**
Prop. of young farmers0.0010.1040.007ns
Farm income−0.1270.048−2.631**
GDP per capita−0.0970.061−1.592ns
Prop. agricultural land−0.4870.089−5.449***
Please wait… references are loading.