Farm structure and environmental context drive farmers’ decisions on the spatial distribution of ecological focus areas in Germany

Ecological Focus Areas (EFAs) were designed as part of the greening strategy of the common agricultural policy to conserve biodiversity in European farmland, prevent soil erosion and improve soil quality. Farmers receive economic support if they dedicate at least 5% of their arable farmland to any type of EFA, which can be selected from a list of options drawn up at the European Union level. However, EFAs have been criticized for failing to achieve their environmental goals and being ineffective in conserving farmland biodiversity, mainly because they are not spatially targeted and because they promote economic rather than ecological considerations in farm management decisions. We used a spatially explicit approach to assess the influence of farm and field context as well as field terrain and soil conditions on the likelihood of whether or not a particular EFA type was implemented in a field. We used a multinomial model approach using field-level land use and management data from 879 farms that complied with the EFA policy in 2019 in the Mulde River Basin in Saxony, Germany. Geospatial environmental information was used to assess which predictor variables (related to farm context, field context or field terrain and soil conditions) increased the probability of a field being assigned to a particular EFA. We tested the hypothesis that productive EFAs are more often implemented on fields that are more suitable for agricultural production and that EFA options that are considered more valuable for biodiversity (e.g. non-productive EFAs) are allocated on fields that are less suitable for agricultural production. We found that farms embedded in landscapes with a low proportion of small woody features or nature conservation areas mainly fulfilled the EFA policy with productive EFAs (e.g. nitrogen fixing crops). Conversely, farms with a higher proportion of small woody features or nature conservation areas were more likely to adopt non-productive EFAs. As predicted, large and compact fields with higher soil fertility and lower erosion risk were assigned to productive EFAs. Non-productive EFAs were placed on small fields in naturally disadvantaged areas. EFA options considered particularly beneficial for biodiversity, such as fallow land, were allocated far away from other semi-natural or nature protection areas. Our results highlight that the lack of spatial targeting of EFAs may result in EFA options being assigned to areas where their relative contribution to conservation goals is lower (e.g. farms with higher shares of protected areas) and absent in areas where they are most needed (e.g. high intensity farms). To ensure that greening policies actually promote biodiversity in European agriculture, incentives are needed to encourage greater uptake of ecologically effective measures on intensively used farms. These should be coupled with additional measures to conserve threatened species with specific habitat requirements.

they are not spatially targeted and because they promote economic rather than ecological considerations in farm management decisions. Objectives We used a spatially explicit approach to assess the influence of farm and field context as well as field terrain and soil conditions on the likelihood of whether or not a particular EFA type was implemented in a field. Methods We used a multinomial model approach using field-level land use and management data from 879 farms that complied with the EFA policy in 2019 in the Mulde River Basin in Saxony, Germany. Geospatial environmental information was used to assess which predictor variables (related to farm context, field context or field terrain and soil conditions) increased the probability of a field being assigned to a particular EFA. We tested the hypothesis that productive EFAs are more often implemented on fields that are more suitable for agricultural production and that EFA options that are considered more valuable for biodiversity (e.g. non-productive EFAs) are allocated on fields that are less suitable for agricultural production. Results We found that farms embedded in landscapes with a low proportion of small woody features or nature conservation areas mainly fulfilled the EFA policy with productive EFAs (e.g. nitrogen fixing crops). Conversely, farms with a higher proportion of small woody features or nature conservation areas were more likely to adopt non-productive EFAs. As predicted, large and compact fields with higher

Introduction
To achieve greater environmental sustainability, the Common Agricultural Policy (CAP) of the European Union (EU) introduced a direct payment scheme (greening) in 2013, aimed at protecting and enhancing biodiversity, water, and soils in European farmland (European Union 2013;European Commission 2022). Mandatory greening practices include: (i) the maintenance of permanent grassland, (ii) the obligation of agricultural diversification of farmland (cultivation of at least two crops if the arable farmland exceeds 10 hectares), and (iii) the allocation of at least 5% of the arable area to ecological focus areas (EFAs) if the arable area of conventional farms exceeds 15 hectares (German Federal Environmental Agency 2014;European Commission 2017. Farmers must comply with these requirements to qualify for the greening payments, except farms already implementing sustainable farming practices such as organic farming (European Commission 2022).The ecological value of EFAs varies across EFA types (Cole et al. 2012;Cormont et al. 2016). While non-productive options such as buffer strips, fallow land, and landscape features are considered most beneficial for conserving biodiversity Sutter et al. 2018;Tarjuelo et al. 2020aTarjuelo et al. , 2020b, less conservation-effective (namely productive) options, such as planting nitrogen-fixing or catch crops, represent more than 70% of the registered EFAs Nilsson et al. 2019;Pe'er et al. 2020). Results from interviews with stakeholders, farmers, and agricultural representatives have shown that farmers consider the EFA regulation complicated (Oppermann 2015;Bonke et al. 2021). Overall, economic and administrative motivations appear to be more important than ecological considerations in defining farm management practices concerning EFA implementation Nitsch et al. 2018;Brown et al. 2021;Bonke et al. 2021). Nevertheless, farmers' motivations for managing their land include the environmental context in which their farm is embedded (Brown et al. 2021;Santos et al. 2021). Farm management decisions may also be heavily constrained by the farm context (e.g. farm size), the field context (e.g. field size), or the field terrain and soil conditions (e.g. soil fertility) (Nitsch et al. 2018). Also, the success of particular land uses in enhancing biodiversity or ecosystem services may highly depend on their proximity to and connectivity with other semi-natural areas (Tscharntke et al. 2005(Tscharntke et al. , 2012. Hence, the consideration of a broader set of factors influencing farm management practices, including the ecological context of the farm, is urgently needed to broaden our understanding of how to manage agricultural landscapes in more sustainable ways (Tzilivakis et al. 2016;Singh and Leppanen 2020;Brown et al. 2021).
Farms across the EU differ in their environmental context according to their geographic location. Additionally, various policy tools exist to support the CAP greening policy and reach the EU environmental targets across member states. For instance, farmers can receive further economic support for areas considered challenging to cultivate (European Commission 2022). In our study region, locations near water bodies (SMEKUL 2021a) and within water protection (SMEKUL 2021b) or Natura 2000 (European Environment Agency 2022) areas have special restrictions regarding management practices or nutrient/pest control inputs. To account for regional differences, member states can select a subset of EFA options considered more suitable for their specific context. Furthermore, there is a weighting factor for each type of EFA to reflect their ecological benefits and implementation costs. Weighting factors can range from 0.3 (e.g. for catch crops) to two (e.g. for hedges) and are used to calculate the exact area per farm relevant for the policy subsidy payments. In Germany, farmers can select among 20 different land use options eligible as EFA Lakes et al. 2020, see Table S1).
In the present study, we used field-level information on land use and management from the Integrated Administration and Control System (IACS) (SMEKUL 2019), as well as geospatial information on small woody features and nature protection areas, to assess which structural, terrain and soil conditions influence the probability of a field being assigned to a specific EFA. Since EFA policies focus on the farm level, our predictors included variables reflecting the farm context in terms of the diversity of additional habitats within the farm boundaries, the relative usability of a specific field for production, and the terrain and soil conditions of the field compared to the other fields of the same farm. Overall, we expected that fields whose structural (size, compactness, proximity to other landscape elements) or terrain and soil conditions (slope, soil fertility) are below the farm's average are more likely to be assigned to an EFA. Additionally, we expected that large farms with a smaller proportion of protected areas or small woody features, whose fields are predominantly suitable for production, would most likely meet the policy requirements with productive EFA options. Since landscape elements are not easily movable and can be used to receive EFA policy payments, we anticipated that farms with landscape features within their boundaries are likely to register them under the policy. Conversely, farms without landscape features within their boundaries are more likely to allocate productive EFA options than other EFA types.

Study region
The study area is the Mulde River Basin, located in the western part of the federal state of Saxony, Germany ( Fig. 1). It covers an area of 5814 km²; the terrain varies from flat to hilly, and the elevation ranges between 24 and 1214 m.a.s.l. (Sachsen Staatsbetrieb Geobasisinformation und Vermessung 2016). The climate is predominantly continental, with total annual precipitation between 570 and 1260 mm and mean annual temperatures between 7.4 and 14.1 °C (DWD 2020). 38% of the study area is covered by arable land, with winter wheat, oilseed rape, winter barley, and maize as predominant crops, followed by permanent grassland (mowing pastures and meadows).

Data preparation
Our data included field-level information on land use and management of 879 farms (41,936 fields) that complied with EFA policy in 2019 (SMEKUL 2019). The median number of fields per farm was 25. The median farm size was 100 hectares, and the median field size was 6 hectares. Farmers could select among 20 land use types to be registered as EFA (Table S1). In our study, a farm consisted of a group of fields registered in the IACS database (SMEKUL 2019) under the same anonymized farm identification number. A field consisted of a parcel with uniform land use (a certain type of crop or EFA). We grouped our predictor variables into three categories; predictors related to the farm context, the field context, or the terrain and soil conditions (Table 1).

Farm context predictors
This category included: the farm area (FARM AREA) calculated as the area of all fields registered under the same anonymized farm identification number, the percentage of the farm covered by Natura 2000 areas (NATURA), the percentage of the farm covered by landscape features registered as EFA (LAND FEAT), the percentage of the farm area covered by fields near small woody features (i.e. in a buffer of 50 m, SWF), the percentage of the farm area designated as naturally constrained (ANC), and the percentage of the farm area designated as areas of importance for drinking water and medicinal spring protection areas (WATERPROT) ( Table 1).

Field context predictors
We used z transformation to standardize the values of each field concerning the farm context. For instance, a positive value of the field area means that the respective field is larger than the average field of the farm to which it belongs. The field context predictors category included: the standard z-score of the area of the field (FIELDAREA), the standard z-score of the compactness index of the field (COMPACTNESS), the standard z-score of the distance of the field to the farm centroid (DISTANCE), a binary predictor of the presence of a small woody feature in the proximity of 50-meter buffer (SWF), information on whether the field belonged to an area designated as naturally constrained (ANC), and the presence of a water body in the proximity of 20-meter buffer (WATERBODY) Fig. 1 Depiction of the geographical location of the study area and the distribution of productive, non-productive, and landscape features as EFAs used in our analysis. The square in the top right represents an inset of the data  (Büchi 2002).

Soil and terrain conditions
This category included values of soil fertility (SOIL FERTILITY), soil moisture (SOIL MOIS-TURE), wind (WINDEROSION), and water erosion (EROSION(KLSR) extracted from the digital soil map at a scale of 1:50,000 from the Saxon State Ministry for Energy, Climate Protection, Environment, and Agriculture (SMEKUL 2021c) ( Table 1). The soil fertility index considered the capillarity, risk of waterlogging, cation exchange capacity, stone content, and plantavailable water in the root space of the soil (Table S2). The soil moisture index was calculated based on parameters such as the soil pore system, climate, and terrain (Table S3). The wind erosion indicator described the erosion hazard depending on the soil type and longterm average wind speed (Table S4). The water erosion risk considered the soil type, slope length, slope inclination, and regenerative erosivity of the soil (Table S5). All predictors were tested for collinearity using pairwise correlation analysis. Since no correlations of Pearson's |r| > 0.7 occurred (Dormann et al. 2013), no variables were excluded (Fig. S1-3).

Statistical analysis
We used multinomial regression models to evaluate the influence of variables related to the farm and field context and the field terrain and soil conditions on EFA allocation in farmland. Multinomial regression is a logistic regression employed to predict the probability that a case is in a particular category when the dependent outcome includes more than two categories. Multinomial models break the outcome variable into a series of comparisons between pairs of categories, on which the probability (odds) of an observation being a member of a specific group rather than the reference group is calculated. In multinomial models, a logistic transformation of the odds (Log (OR)) is the dependent variable. The odds value can range from 0 to infinity and express the likelihood that an observation is a member of the reference group rather than a member of the group used for comparison. It holds that odds = p/ (1 − p). The odds ratio (OR) estimates the change in the odds of membership in the target group for a one-unit increase in the predictor and it is calculated by using the regression coefficient of the predictor as the exponent: log(odds) = logit(p) = log(p∕(1 − p)) = 0 + 1 x 1 + … Soil and Terrain conditions SOIL FERTILITY Standardized (z-score) value of the field soil fertility concerning the mean field soil fertility of the respective farm.

SMEKUL (2021c)
SLOPE Standardized (z-score) value of the field slope concerning the mean field slope of the respective farm.

SMEKUL (2021c)
SOILMOISTURE Standardized (z-score) value of the field soil moisture concerning the mean field soil moisture of the respective farm.

SMEKUL (2021c)
WIND EROSION Standardized (z-score) value of the field wind erosion concerning the mean field wind erosion of the respective farm.

SMEKUL (2021c)
EROSION (KLSR) Standardized (z-score) value of the field water erosion risk concerning the mean field erosion risk of the respective farm.

SMEKUL (2021c)
Where p = the probability that a case is in a particular category, β 0 is the intercept, and β 1 , … are the regression coefficients of x 1 . In order to identify specific responses according to the EFA group (e.g. productive, non-productive, landscape features) or EFA type (e.g. catch crops, nitrogen-fixing crops, etc.), we ran four multinomial models. All models included the same predictor variables (described in Table 1) and had fields dedicated to agricultural production (defined as NOEFA) as a reference value. We used Model 1 to identify differences in drivers of EFA allocation according to their productive value. Besides the reference value (NOEFA), Model 1 included three additional response categories: (i) productive EFA, (ii) non-productive EFA, and (iii) landscape features registered as EFA. We used Model 2, Model 3, and Model 4 to identify important drivers for specific EFAs (See Figs. S4, S5, S6). Model 2 incorporated only productive EFAs (c = 6) as the response categories. Model 3 included response non-productive EFAs (c = 6), while Model 4 had EFA landscape features (c = 8) as response categories. Given that we only considered farms with at least one field dedicated to one of our focal EFAs (defined by the model), all models differed in the number of farms and fields used as input (Table 2). All data analysis and data processing was conducted in R version 4.0.5 (R Core Team 2021) using the packages nnet (Fox and Weisberg 2019), car (Fox and Weisberg 2019), caret (Kuhn 2008), dplyr (Wickham et al. 2019), raster (Hijmans and van Etten 2012), sf (Pebesma 2018), and terra (Hijmans 2023). Details of the code can be found in the GitHub repository related to this manuscript (https:// github. com/ vivialarc on/ EFAal locat ionin farml and. git).

Results
7% of all fields were allocated with a productive EFA, 6% with a non-productive EFA, and 5% to landscape features registered as EFA. Among productive EFAs, catch crops, and nitrogen-fixing were the most commonly adopted options. Fallow land and buffer strips were the most frequently allocated non-productive EFAs, and hedges and rows of trees were often registered as landscape features. Farms covered with a low percentage of Natura 2000 (NATURA), landscape features registered as EFA (LAND FEAT), small woody features (SWF), or water protection areas (WATER PROT) mostly allocated productive EFAs (Fig. 2). Conversely, farms with higher values on the same variables were more prone to allocate non-productive EFAs or landscape features. Small and linear fields (FIELD AREA, COM-PACTNESS) had higher probabilities of being allocated with non-productive EFAs or landscape features (Fig. 2). Non-productive EFAs and landscape features were allocated further away from small woody features (SWF) or water bodies (WATER BODY), but they were frequently present in fields designated as naturally constrained for agricultural production (ANC). Conversely, large and rounded fields were most prone to be assigned with productive EFAs. Contrary to our expectations, EFA fields were not allocated in fields with low soil fertility (SOIL FERTILITY). Nevertheless, both productive and nonproductive EFAs were allocated in fields with lower risks of water erosion (EROSION (KLSR)).
Among productive options, catch crops and nitrogen-fixing crops responded differently to the same variables (Table 3, Fig. S4). While nitrogen crops were frequently allocated on small farms (FARM AREA) with a lower percentage of landscape features (LAND FEAT), small woody features (SWF) and water protection areas (WATER PROT), catch crops were allocated in farms with higher values of Natura 2000 (NATURA), landscape features (LAND FEAT (%)), and water protection areas (WATER PROT) ( Table 3, Fig. S4). Catch crops were frequently assigned to large and rounded fields (FIELD AREA, COMPACTNESS) in naturally constrained (ANC) areas. Fallow land was more likely to be allocated in farms with lower values of NATURA, LAND FEAT, and WATER PROT), and in small and linearly shaped fields (Table 3, Fig. S5). Conversely, buffer strips were most often established on farms with higher values of NATURA, LAND FEAT, SWF or WATER PROT. Hedges and rows of trees were most abundant in farms with low percentages of SWF and ANC, but within these farms, to fields located in ANC areas (Table 3, Fig. S6). Most variables related to terrain and soil conditions had non-significant effects for the majority of EFA types, with the exception of slope and erosion KLSR.

Discussion
In this study, we used a spatially-explicit approach to analyze how the farm context, the field context, and the terrain and soil conditions influenced the allocation of EFAs in farmland. We demonstrated that the farm and the field contexts, rather than the terrain and soil conditions, are decisive for the assignment of a particular EFA type to a particular field. In our study region, EFA types known to benefit biodiversity, such as fallow land, buffer strips, and landscape features (Martin et al. 2019;Pe'er et al. 2022), were used on farms already embedded in contexts with a higher cover of nature conservation and non-crop habitats (Fig. 2). In contrast, productive EFAs such as catch crops or nitrogen-fixing crops were often found on farms where cropland dominated (Fig. 2). These results have important implications for conserving farmland biodiversity. Considering that the effectiveness of conservation measures in agricultural areas is highest in landscapes with intermediate levels of land use heterogeneity ("intermediate landscape-complexity hypothesis", Tscharntke et al. 2012), the potential of EFAs to contribute to biodiversity conservation in farmland can be strongly moderated by their spatial distribution (Sutter et al. 2018;Concepción and Díaz 2019;Concepción et al. 2020).
EFAs on farms embedded in complex landscapes may contribute little to restoring existing biodiversity on this land. On the other hand, EFAs on farms embedded in simplified landscapes may also contribute little to restoring the biodiversity that has already been lost in these areas. The allocation of EFAs in an intensively used area, however, has the potential to increase landscape complexity, for Fig. 2 Model coefficients (logarithm of odd ratio) of the multinomial regression used to evaluate the drivers of EFA allocation according to the EFA group (model 1). Explanatory variables are grouped into three categories (farm context, field context, soil and terrain conditions). Colored circles indicate p-values lower than or equal to 0.05. Farm context indicators (except farm area) reflected the percentage of farm area covered by: Natura 2000 designated area (NATURA (%)), landscape features registered as EFA (LAND FEAT (%)), small woody features (SWF (%)), naturally constrained areas (ANC (%)) or water protection areas (WATER PROT (%)). SWF indicates the presence (1) or absence (0) of a small woody feature in a buffer of 50 m around the field. ANC indicates whether the field is in a naturally constrained area. WATER BODY indicates the presence (1) or absence (0) of a water body in a buffer of 50 m around the field. EROSION (KLSR) refers to the water erosion risk instance by increasing crop diversity or the number of non-crop elements within farm boundaries (Concepción et al. 2012). Nevertheless, EFAs should not be used as the only tool to conserve biodiversity in farmlands, as their benefits for threatened or specialist species depend highly on their ecological quality (Herzog et al. 2017;Pfiffner et al. 2018). Further, many EFAs only benefit farmland generalist species (Kleijn et al. 2006;Aviron et al. 2007Aviron et al. , 2009) and need to be complemented with additional nature protection efforts or improved and tailored agri-environmental schemes in order to protect and enhance a wider range of biodiversity (Aviron et al. 2009).
Farm area in our study region was a determining factor only for a subset of the EFA options (Table 3). Fallow land and nitrogen-fixing crops were most often adopted by small farms, while hedges were most frequently used by larger farms (Table 3). Large farms with a higher cover of landscape features in their boundaries used these elements to fulfill the policy requirements. Interestingly, farms that did not fulfill their EFA quota with landscape elements tended to allocate productive instead of non-productive EFA (Fig. 2). These results support previous findings on the prevalence of economic rather than ecological considerations on farmers' decisions on greening policy allocation . Productive EFAs (e.g. catch crops and nitrogen-fixing crops) are mainly targeted to reduce soil erosion and facilitate nutrient uptake (Cerdà et al. 2022;Quintarelli et al. 2022) and may help to facilitate crop production in the same field on the following year (Brown et al. 2021;Wittstock et al. 2022). On the other hand, the placement of non-productive EFAs and landscape features is not spatially flexible and may appear disadvantageous from the farmer's point of view, given that working with large machinery may be more difficult if individual trees or hedges are present at the field boundaries. While productive EFA options were commonly assigned to large, rounded fields with lower soil erosion hazard, non-productive EFAs were allocated in small, linear fields in areas considered disadvantaged for crop production (Fig. 2). A similar pattern was found by Paulus et al. (2022) for Agri-environmental schemes allocation in our study region. Paulus et al. (2022) demonstrated that Agrienvironmental schemes are more often allocated in areas with low potential for agricultural intensification, supporting the notion that the lack of careful spatial targeting of agricultural policies may diminish their ecological contribution.

Conclusions
Farmers' economic considerations play an essential role in deciding which type of EFA to choose (e.g. productive versus non-productive) and where to place it (e.g. large and rounded fields versus small, linear, and less fertile fields). Our results demonstrate that the farm and the field context have a significant impact on the spatial distribution of greening measures within the farm boundaries. These findings are especially important for possible future spatial targeting of conservation and nature protection measures (including eco-schemes) on farmland. This relates to several measures that are spatially flexible, such as fallow land and buffer stripsoptions that are widely recognized in expert-based assessments as highly valuable for biodiversity and ecosystem services (Pe'er et al. 2017; Traba and Morales 2019).