Background

Nitrogen (N) is a key input for crop production, but has adverse effects on the environment through the emissions of reactive N. Since the development of the Haber–Bosch process for technical ammonia synthesis around a century ago, humans have intervened in the N cycle more than in any other geochemical elemental cycle [1]. The production of N fertilisers on an industrial scale strongly increased global agricultural productivity in most regions of the world [2]. However, a large proportion of this N is lost to the environment: on the global scale less than 20% of the fertiliser N applied to agricultural land is actually consumed by humans in meat, but over 80% is lost to the environment, much of it into the air and into rivers, groundwater and coastal seas [3, 4]. Thus, it is a big challenge to produce enough food for the growing global population while at the same time minimising the emission of reactive N to the environment.

Given this situation, agriculture N budgets have been introduced to clarify the nutrient flows, to quantify the potential N losses, to indicate the pressure on environmental media, to raise awareness in nutrient management, and to enforce and monitor nutrient mitigation regulations in practice [5]. For example, “Gross nitrogen balance” and “Mineral fertiliser consumption” are two (of 28) agri-environmental indicators used by EUROSTAT and the European Environmental Agency to evaluate agricultural impact on the environment and to monitor the effects of agri-policy measures in the EU [6]. The European Environment Agency [7] regarded the surplus of N as the best overall environmental indicator for the changes in the agricultural impact on the environment. De Vries et al. [8] also stated that N budgets of agricultural systems give important information for assessing the impact of N inputs on the environment, and for identifying levers for action. Since the 1990s, agriculture in Western Europe has managed to reduce its N surpluses, due to national and EU policies. The reduction in N budget surplus and in groundwater nitrate concentration in Denmark illustrates the success that can be achieved if effective measures are taken [9].

The literature offers numerous methods, definitions of system boundaries and calculation schemes to estimate an N budget surplus. According to the Eurostat/OECD Nutrient Budget Handbook [10], an N budget quantifies all major N flows across all sectors and media within given boundaries, and flows across these boundaries, within a given time frame (typically one year), as well as the changes of N stocks in the relevant sectors and media. An N balance is ideally a closed N budget, i.e. all N flows can be explained and the balance equation is then [10]:

$$N\left( {\text{input}} \right) \, - \, N\left( {\text{output}} \right) \pm N\left( {\text{stock change}} \right) \, = \, 0.$$
(1)

In practice, however, essential flows of an N balance such as gaseous losses, nitrate leaching and soil stock changes cannot be determined. Therefore, the N budget surplus, as given by total inputs minus total outputs, represents the N quantity which is volatilised to the air, leached to groundwater, discharged to surface water, or stored in the soil. Therefore, an N budget surplus value includes all N components emitted into the environment or accumulated in the organic soil substance.

According to the objective and the system boundary there are basically three approaches to estimate nutrient budgets [8, 10, 11]: (i) farm budget (called farm-gate budget by Oenema et al. [5]), which records the amounts of N in all kinds of products that enter and leave the farm via the farm-gate; the equivalent on the national scale is the agricultural system budget, also termed sector balance [12]; (ii) Land budget (gross N balances by the OECD [13]), which records all N that enters a farmland and leaves the farmland by crop products; and (iii) Soil budget (soil surface budget by Oenema et al. [5]), where in contrast to the land N budget, the N inputs via animal manure are adjusted for gaseous losses of N emissions in housing and manure management systems. The latter budget is estimated in this work. Due to the lack of an exact definition of the terms N balance and N budget as it is described here, in older studies these terms are not always used consistently.

While the N budget surplus indicates the potential for N losses from soils into the atmosphere and the hydrosphere quantitatively, the Nitrogen Use Efficiency (NUE), or N recovery efficiency as termed by Conant et al. [14], evaluates the efficiency of the use of the resource N in crop production [11].

Over the last three decades, a huge number of agricultural N budgets and N balances have been presented. The approaches differ by many issues such as system boundaries, geographical entities, budget period, input and output budget terms, data sources, as well as factors and coefficients used for calculation. With regard to the regional focus of the N budget presented in our study, the following examples highlight the manifold approaches of spatialised N budgets and differences in employed data basis. Eriksson et al. [15] drew up N budgets of Polish agriculture 1960–2000 based on the administrative units of voivodeships (provinces). The regional N surpluses for the ten Danish geo-regions were derived by Hansen et al. [9]. Grizzetti et al. [16] produced EU-15 maps of spatialised inputs from agriculture and gross N balance by combining the EU Farm Structure Survey 2000 data with the CORINE Land Cover map. Wang et al. [17] estimated N balances for 2426 counties in China from statistical survey data.

Leip et al. [11] calculated soil N budgets for agriculture in the EU-27 using the agro-economic model CAPRI. Özbek et al. [18] assessed N budgets of agriculture in Turkey at the level of 27 NUTS 2 regions by an improved Eurostat/OECD method. Lassaletta et al. [19] calculated N budgets for the river Ebro catchment detailed for territorial units derived by overlaying maps of main agricultural use, livestock density, and hydrographic characteristics. Bach et al. [20] assessed agricultural nutrient balances from agricultural census data exemplary for 1109 municipalities of Baden-Wuerttemberg (Germany). A comparison of N budgets for the year 2000 of agro-ecosystems is made for the EU-27 countries by four models with different complexity and data requirements [8, 21]: The IDEAg model uses 206,000 HSMUs (Homogeneous Spatial Mapping Units), while the INTEGRATOR model operates for 39,052 NCUs (NitroEurope Calculation Units, clusters of 1-km2 grid cells that are characterised by similar environmental and agronomic conditions). While the total N surplus on the EU level is similar between the four models, the difference in N output fluxes between models is large on a regional scale.

The N soil surface budget surplus is applied in particular as a core indicator for assessing the nitrate input into groundwater and surface waters from agricultural land use. The Nitrate Atlas for Germany [22] presented one of the first such model approaches, using the N surplus of the utilised agricultural area (UAA) for 3 km × 3 km grids. Campling et al. [23] implemented spatialised soil surface N balances as pressure indicator for river basins in the EU-15 using the CORINE Land Cover map. The MITERRA-EUROPE model assesses the agricultural N balances and the loss of reactive N species to the environment for the EU Member States for the 234 NUTS 2 regions [24]. Bouraoui and Grizzetti [25] reviewed some of the models applied in Europe to assess the contribution of diffuse N water pollution and the effectiveness of mitigation measures. The risk of water contamination was modelled for the 3500 polygons of the 1:1,000,000 Map of Soil Landscapes of Canada by calculating the “residual soil nitrogen” (RSN) [26]. Howarth et al. [27] introduced the quantity “net anthropogenic nitrogen input” (NANI) as a well-correlated measure of N transported in rivers, which was also used by Lassaletta et al. [19] for the river Ebro catchment (Spain).

The listed publications illustrate the huge diversity of the spatial resolution and underlying data basis of approaches introduced to calculate an N budget. Given the heterogeneity of definitions, approaches and methods, it is not surprising that the results of N budgets from various studies may differ substantially and can only be compared with great limitations or not at all. N surplus data should therefore be interpreted less in terms of its absolute value than in terms of its change over time and/or its spatial dispersion within a region. As an example, the German Federal Government construes the trend in national N surplus as evidence of the effectiveness of the measures taken to reduce groundwater nitrate pollution, and reports these every 4 years to the EU Commission on the implementation of the EU Nitrates Directive. The N surplus of Germany’s agricultural N budget also forms a core indicator for evaluating the sustainability of agriculture in Germany [28]. The EU Water Framework Directive calls on all Member States to submit an inventory of nutrient emissions and river loads for their river basins. In Germany, the river basin management models MONERIS [29] and MoRE [30] are frequently applied for Germany to estimate N emissions into water bodies from various sources and to prioritise measures to reduce N losses into the hydrosphere. In both models, the N soil surface surplus of agricultural land is a core quantity to evaluate the spatial risk of nitrate losses into groundwater. For this purpose, the model MoRE calculates the N soil surface surplus for districts [31], while the model MONERIS involves the N input and output on the community level for Germany by the RAUMIS approach [32].

We present the methodology and the results of N soil surface budget calculation for the district administrative units (Nomenclature des unités territoriales statistiques, NUTS 3 level, “Kreise”) in Germany for the time series 1995 to 2017. To our knowledge, it is the first time that the N surplus for a territory is determined for a time series of more than 20 years with such high spatial resolution. We use the data to focus on five aspects: the status and trend of nitrogen input and budget surplus in German agriculture, the contribution of biogas production to the N budget surplus, the regional spread of N budget surplus and NUE, the change in regional nitrogen budget surplus, and finally the potential effect of the German Fertiliser Ordinance [33] on the regional distribution of the N surpluses.

Materials and methods

In general, nutrient budgets can be established for each agricultural production system: crop production, livestock production, and biogas (energy) generation. The latter has gained increasingly relevance in German agriculture since the Renewable Energy Sources Act was revised in 2004. The N budget surpluses are linked as follows:

$${\text{Agricultural system budget surplus }} = {\text{soil surface budget surplus }} + {\text{ livestock budget surplus }} + {\text{ biogas budget surplus}} .$$
(2)

Table 1 presents the N budget terms for the individual budgets and their allocation in the overall N budget sheet. For N budgets below the national level, i.e. for federal states, districts or municipalities, only N soil surface budgets can be determined due to limited data availability in Germany. The methodology and the input and output terms differ in some details from the calculation scheme for the N soil budget according to the OECD [12, 13] especially with regard to the biogas component. The NUE is finally calculated as the ratio of total output with harvest withdrawal to total input.

Table 1 Scheme of the N agricultural system (sectoral) budget and the linking between its components soil surface budget, livestock budget and biogas budget

Data base

Most of the data required to calculate the input, internal flow and output of N in agricultural production (e.g. cropping acreages, livestock numbers, and mineral fertiliser sales) are provided by the Federal Statistical Office (www.destatis.de). Over the past two decades, complete Farm Structure Surveys (FSS) covering all or a representative sample of agricultural holdings are available for 1995, 1996, 1999, 2001, 2003, 2007, 2010, 2013, and 2016. Full data sets on UAA, arable land by crop types and livestock numbers for districts are limited to the years with a FSS. For budget calculation of interim years, the figures for the districts were interpolated based on annual updates for the federal states. Data on harvest yields of most field crops are reported annually for districts. Coefficients of N contents in harvested crops and leguminous N fixation rates are taken from the German Fertiliser Ordinance [33]. Livestock N excretion rates, N input with bedding and gaseous N losses from manure management (housing and storage) (partial N emission factors) for the calculation of the N supply with manure are provided by the German Greenhouse Gas Emissions Inventory [34] as annual and district-specific values for livestock categories. Sources of statistical data and N coefficients as well as the budget calculation methodology are documented in detail in Additional file 1.

Mineral fertiliser input

Reliable statistical data on mineral fertiliser sales are published only for Germany as a whole and therefore estimates have to be made for regional N budgets entities. Assessing the appropriate mineral fertilising quantity is a major problem in N balancing of districts and the most sensitive factor of regional N surplus estimations. In our approach the amount of mineral fertiliser for districts are derived from the national total. The approach attempts to stimulate the farmer’s N management planning in a simplified way and starts with the N demand of the cultivated crops. Crop N demand is determined from the N harvest withdrawal multiplied by a fertiliser input factor (IF) to allow for the fact that regularly more N is applied with inorganic and organic fertilisers and legume N fixation than the N harvest withdrawal in Germany corresponds to. The crop N demand is met by mineral fertilisers (min_fert), on-farm organic N fertilisers (org_fert, i.e. liquid manure, solid manure, digestate) and biological N fixation by legumes (fix). However, only a proportion of N in organic fertilisers and from biological N fixation is available for plant nutrition and thus effective for fertilising, which respect the utilisation factors (UF):

$${\text{N}}\left( {{\text{tot}}\_{\text{fert}}} \right) \, = {\text{ N}}\left( {\text{harvest}} \right) \, \times {\text{ IF }} = {\text{ N}}\left( {{ \hbox{min} }\_{\text{fert}}} \right) \, + {\text{ N}}\left( {{\text{org}}\_{\text{fert}}} \right) \, \times {\text{ UF1 }} + {\text{ N}}\left( {\text{fix}} \right) \, \times {\text{ UF2}} .$$
(3)

Other organic N fertilisers (biowaste compost, sewage sludge, meat-and-bone meal) are not taken into account in Eq. (3), because the N supply with these substances is generally very low (on the national level) and the N is only slowly released, so that the contribution of these N sources to the current plant N supply is negligible. Due to the broad variability of real data, in the following UF1 is set as conservative assumption to 0.6 and UF2 to 0.8. This means, 40% and 20%, respectively, of the N from animal manure and biological N fixation are not taken into account into the N demand analysis according to the German Fertiliser Ordinance [33]. On the national level the annual values of the variables N(min_fert), N(harvest), N(org_fert) and N(fix) are given and then IF was calculated by Eq. (3) for each year of the time series 1995 to 2017, ranging from 0.98 to 1.35 with a mean value of 1.13. For details of the derivation of UF1, UF2 and IF, see Section 2.4.1 in Additional file 1.

Subsequently, for district N budgets the individual district mineral N fertilising as target value was calculated as remaining fertilising quantity after a share of the district total N needs are met with organic N fertilisers and biological N fixation, calculated from the livestock numbers and the legume cropping area of the given district:

$${\text{N}}\left( {{ \hbox{min} }\_{\text{fert}}} \right) \, = {\text{ N}}\left( {\text{harvest}} \right) \, \times {\text{ IF }}{-}{\text{ N}}\left( {{\text{org}}\_{\text{fert}}} \right) \, \times {\text{ UF1 }}{-}{\text{ N}}\left( {\text{fix}} \right) \, \times {\text{ UF2}} .$$
(4)

The variables on the right-hand side of Eq. (4) are known for each district (see below), the parameters IF, UF1 and UF2 were kept identical for all districts. This approach ensures that the amount of N(min_fert) summed over all districts is equal to the national total of mineral fertiliser sales.

Manure input (on-farm production and transfer)

The N amount in liquid and solid manure from livestock farming is calculated by multiplying the number of livestock units by the N excretion rate per animal place and year, for dairy cows the N excretion is a function of the annual milk yield. Thereof the gaseous N losses during livestock housing and manure storage are subtracted (separately for ammonia and for the other gaseous N species), whereas gaseous N losses from manure spreading on the field are not subtracted, to get the N input with manure to the N soil surface budget.

Since 2013, the federal states with highest regional concentration of livestock production (North Rhine-Westphalia, Lower Saxony, and Schleswig-Holstein) report the nutrient transfer within the districts of the federal state, between federal states and between these states and foreign countries. The import of manure from the Netherlands to Germany was assessed based on data of the Statistical Office of the Netherlands in combination with the nutrient reports of North Rhine-Westphalia and Lower Saxony.

Biogas digestate input

The German Federal Network Agency, the four German power transmission system operators and the German Energy Agency (Dena) provide inventory data on the location of biogas plants and their electric energy (GWhel) annually fed into the grid (for plants with on-site power generation) or their methane volume (Mm3) annually fed into the natural gas grid (for biomethane feed-in plants). Mean composition of the used biogas substrates (separately for six biomass plants, manure, and co-substrates) are taken from a national review of biogas activity data [35]. Substrate quantities are derived from the specific biogas yield of the substrates resulting in coefficients of the substrate N input per GWhel of electric energy fed-in (equivalent per Mm3 biomethane fed-in) for energy crops, manure and biowaste. The N input to crop production with the application of digestates is then taken as N input with substrates into the biogas plants minus the overall gaseous N emissions in the biogas process chain according to [34]. The contribution of biogas production to the N soil surface budget surplus (Fig. 2) was calculated based on the assumption that each N input source contributes to the N soil surface budget surplus according to their respective share from the total N input into the soil.

Other organic fertiliser input

Other organic fertilisers include compost from urban green waste, sewage sludge and meat-and-bone meal for use as fertiliser or soil amelioration. Only national data are available on compost and meat-and-bone meal use in agriculture; sewage sludge statistics are published for the federal states. The N quantity of these fertilisers is generally very low, thus the authors consider the error resulting from the even distribution over the UAA of all districts as acceptable.

Biological N fixation

The amount of N fixed by symbiotic bacteria in the soil is calculated as the acreage of legumes (pulses, clover, alfalfa and grass–legume mixtures) multiplied by the crop-specific N fixation rates [33].

Seeds and planting material input

Input with seeds of coarse grains (cereals, maize, and legumes) and potato planting material is accounted [33].

Atmospheric N deposition

N input with atmospheric deposition is derived from the PINETI-3 project (Pollutant INput and EcosysTem Impact) modelling the annual total deposition of NOx and NHy in Germany with 1 km x 1 km resolution for ten vegetation classes, from 2000 to 2015 [36]. PINETI-3 grid maps were superimposed with a land cover model of Germany and for each grid cell the fractions for the vegetation classes “arable land” and “grassland” were determined. Then this map was superimposed on the administrative boundary map and the annual NOx-N and NHy-N deposition rates per hectare UAA were averaged for the districts.

Gaseous NH3 emissions are almost entirely caused by agriculture. Thus, the NHy deposition on agricultural land originating from NH3 emissions from organic and mineral fertiliser application to the soil represents an internal N cycle within the soil surface N budget. In Germany approx. 42% of total N volatilisation from agriculture occurs during livestock housing and manure storage, 57% results from the application of all kinds of fertilisers, and 1% is released from biogas plants [37]. The NH3 emission attributed to fertiliser application must be subtracted from the gross NHy deposition as otherwise it is double counted on the input side and the assessed soil surface surplus would be overestimated. For the spatialised N soil surface budget it is assumed as a simplification that the NHy deposition on the UAA of each district originates exclusively from NH3 emissions within the same district. Thus, for each district the fraction of NH3 emission from fertiliser application from the total agricultural NH3 emissions was calculated. Doing this, non-agricultural NH3 emissions were neglected. To calculate the net NHy deposition, this fraction was multiplied with the gross NHy deposition. The net deposition corresponds to the NHy input to the soil surface budget from livestock production and biogas production via the atmosphere. In Table 2 this internal flow is taken into account as ‘NH3 emissions from fertiliser application returned to UAA’ as part of the withdrawal.

Table 2 N soil surface budget in Germany (averaged over all districts, mean 2015–2017)

Harvest withdrawal

To calculate the N harvest withdrawal with marketed crops, fodder crops and energy plants, the acreage of the particular crops is multiplied either by the yield per hectare and the N content in the harvested crop or for field crops without yield data by a fixed withdrawal quantity per hectare [33]. The only by-product considered is straw for bedding in livestock production. No data are available on the thermal use of straw, animal feeding with straw and the harvest and utilisation of other by-products. For soil surface N balancing it is not necessary to split the harvest withdrawal for utilisation as fodder crops (on-farm use) and energy plants. However, to compare the soil surface N budget for districts with the German sectoral budget the withdrawal was attributed as follows: all cereals and legumes (harvested as grains), potatoes, sugar beet and rape seed are marketed crops, grassland and green harvested plants are considered as fodder crops. The withdrawal used as energy crops in biogas plants was calculated based on electric work specific substrate mixes for each federal state according to KTBL [35] and the fed-in electric power from biogas plants. The amounts were subtracted from the harvest withdrawal in marketed and fodder crops.

Territorial status of districts

Territorial status of the districts refers to 1 January 2016 for the entire time series 1995 to 2017. It comprises 402 regional entities with city federal states (Berlin, Hamburg, and Bremen) treated as districts. For N budget calculation, the 402 entities were aggregated to 299 so-called “district regions” (map see Additional file 1). The aggregation should avoid methodological errors when calculating the N budget for urban districts for which the FSS often specifies only little UAA. The FSS applies the “holding location principle”, i.e. all cropping acreage and livestock figures are attributed to the district where the holding is registered, wherever the cropland or the livestock housing is actually located. As a result, the true values of acreages and livestock numbers within a district can deviate from the FSS data, distorting the spatial allocation of the N surpluses to the district’s UAA. Additionally, due to the large-scale farming typical for eastern Germany, for a number of districts no UAA or livestock numbers are published in the FSS for reasons of data privacy. These distortions can be at least partially avoided by merging small districts with the larger adjoining district.

Results

Status and trend of nitrogen input and surplus in Germany

At present (mean 2015–2017), the largest input of 104 kg N ha−1 a−1 is made with mineral fertiliser followed by organic fertilisers from livestock and biogas production with nearly 89 kg N ha−1 a−1 (Table 2). Note that the N budget terms always refer to the UAA, totalled over the districts. N fixation by legumes contributes about 13 kg N ha−1 a−1 to the N input, and the atmospheric gross deposition amounts to 14 kg N ha−1 a−1. The other inputs are of minor importance with a total of about 8 kg N ha−1 a−1. Just under half of the N withdrawal leaves the agricultural land with marketed field crops (67 kg N ha−1 a−1), while the remaining 77 kg N ha−1 a−1 is used on-farm as fodder in livestock production or as energy crops. 5 kg N ha−1 a−1 NH3 emissions from fertiliser application, which are returned to the soil with NHy deposition are also taken into account into the withdrawal. An input total of around 227 kg N ha−1 a−1 is opposed by a withdrawal from the UAA of around 149 kg N ha−1 a−1, which results in an N soil surface budget surplus of 77 kg N ha−1 a−1. The NUE in German plant production, thus amounts to 0.658.

Both, N input and output are subject to considerable fluctuations over time (Fig. 1). From 1995 to 2017, the N supply in the N soil surface budget varies between 201 and 232 kg N ha−1 a−1. The fluctuation mainly reflects the annual variation in mineral fertiliser sales volumes, but cannot necessarily be interpreted as actual changes in annual mineral fertiliser application to the UAA. Rather, it can be assumed that farmers and trading companies react to price alterations for mineral fertilisers by building up or reducing stocks. However, the N budget cannot capture this issue. N withdrawal with harvested crops (figures not shown) varies significantly more due to weather variability over the period 1995–2017. The span is given by the extremes of the drought year 2003 with an N withdrawal of only 115 kg N ha−1 a−1 and the overall record-yield year 2014 with an N withdrawal of 162 kg N ha−1 a−1.

Fig. 1
figure 1

Input terms of the annual nitrogen soil surface budget of districts aggregated for Germany 1995 to 2017

The time series of the N soil surface budget surplus depicts the annual fluctuation of both, the N input and output, with a maximum of 91 kg N ha−1 a−1 in 2003 and a minimum of 53 kg N ha−1 a−1 in 2009 (Fig. 2). The entire period 1995 to 2017 can be separated into two parts. From 1995 up to 2009 a slight downward tendency in the N budget surplus can be observed, while for the recent period from 2009 onwards no tendency is visible; but none of the regression slopes are statistically significant.

Fig. 2
figure 2

Annual nitrogen soil surface budget surplus for districts and thereof from digestates aggregated for Germany 1995 to 2017

Contribution of biogas production to N budget surplus

The share of N surplus originating from biogas production is proportional to the fraction of energy cropping acreage of the total UAA. Figure 2 shows the increasing share of the N turnover in biogas production for the N soil surface budget from the year 2000 onwards. Not accounting for the return of N with digestate from energy crops in biogas plants to the crop production would lead to a systematic underestimation of the N soil surface budget surplus.

In 2017, the input to the UAA with digestates contributes 562,900 t N equivalent to 15% of the total N input in crop production, thereof 345,900 t N are added by energy crops and biowaste and 217,000 t N by manure used in biogas plants. The share of manure from livestock farming passing a biogas plant before spreading reached 16% in 2017.

The N turnover in biogas plants is virtually loss-free, this means no N withdrawal takes place with biogas production. Theoretically, the amount of 308,000 t N from energy crops (without biowaste) in digestates, assuming an utilisation factor of 0.6, could replace 184,800 t of mineral N fertiliser. However, such a decline cannot be observed in the time series of the annual mineral fertiliser input (Fig. 1).

Regional spread of N budget surplus and NUE

The N budget surpluses (mean 2015–2017) of districts spreads over a wide range from 26 kg N ha−1 a−1 to 162 kg N ha−1 a−1 (Fig. 3a), which is more than sixfold between the minimum and the maximum. The NUE of crop production in the districts spans from 0.53 to 0.79 (Fig. 3b). The districts with a high N surplus show a low NUE of crop production and vice versa (Fig. 3a, b). The highest values of N surplus (and lowest NUE values) appear in regions with intensive livestock farming such as in north-western Germany and parts of Bavaria. In contrast, the eastern and central districts in Germany, dominated by cash crop farming, are characterised by medium to low N surplus values.

Fig. 3
figure 3

Nitrogen budget surplus of nitrogen soil surface budgets of districts in Germany, mean 2015–2017 (a), and nitrogen use efficiency, mean 2015–2017 (b), change in N surplus from 1995–1997 to 2015–2017 per utilised agricultural area (UAA) (c), and nitrogen surplus change due to manure transfer and import from the Netherlands, mean 2015–2017 (d)

Change in regional nitrogen budget surplus

The change in N soil surface budget surplus from mean 1995–1997 to mean 2015–2017 shows a further widening of the regional spread. The districts with an N surplus below average show nearly constant (change ± 5 kg N ha−1 a−1) or decreasing N surpluses (decline more than − 5 kg N ha−1 a−1) over the 23-year period (Fig. 3c). However, in the districts in north-west Germany with high livestock density and therefore a high to very high N surplus in most cases the N surplus even increased. In addition, in the districts of Lower Saxony and North Rhine-Westphalia adjacent to the border to the Netherlands, manure imports from the Netherlands lead to an even higher N surplus. However, increases of N surplus are partially transferred to the arable farming regions in the eastern part of Lower Saxony, Mecklenburg-Western Pomerania and Saxony-Anhalt via manure export. This means, the huge N surplus from intensive pig and poultry production in north-west Germany is distributed over a larger area.

Effect of German Fertiliser Ordinance on regional N surplus distribution

The German Fertiliser Ordinance limits the organic N application to 170 kg N ha−1 a−1 for the individual farms, which is approximately equivalent to a livestock density of 1.8 livestock units (LU) per hectare on district average. Thus, for districts with a livestock density of > 1.8 LU ha−1 this regulation forces the farmers to export their excess manure > 170 kg N ha−1 a−1 to other regions with capacity to uptake additional organic fertiliser N under the restrictions as given by the German Fertiliser Ordinance. Figure 4 illustrates that the N soil surface budget surplus (mean 2015–2017) is closely correlated to the livestock density for all districts up to 1.8 LU ha−1 with an increase of N surplus by nearly 60 kg N ha−1 a−1 for an increase in livestock density by one LU ha−1 UAA (grey symbols in Fig. 4). According to the nutrient reports of the states North Rhine-Westphalia and Lower Saxony, in 11 districts (with > 1.8 LU ha−1) the net manure transfer reduces the N surplus by more than − 10 kg N ha−1 a−1 (mean 2015–2017), the maximum reported reduction is -99 kg N ha−1 a−1 for the district Vechta. For these 11 districts the N surplus would hypothetically grow by nearly 60 kg N ha−1 a−1 per LU ha−1 and reach up to nearly 240 kg N ha−1 a−1 (red symbols) without the limitation of organic N fertilising. Thus, the threshold of 170 kg N ha−1 a−1 organic N application set by the Fertiliser Ordinance forces the farms to discount a part of their organic N by manure transfer and therefore confines the N surplus increase to a moderate rise in these districts (blue symbols). Due to the small number of districts with > 1.8 LU ha−1 these results should be interpreted cautious. Changes in N soil surface surplus might be driven too by other political measures and local water protection actions in the districts and might also be partly caused by market-driven changes of the agricultural structure.

Fig. 4
figure 4

Nitrogen soil surface budget surplus of districts in Germany as a function of livestock unit (LU) density (mean 2015–2017) and the effect of N surplus reduction by manure transfer (withdrawal) for the districts with more than 1.8 livestock units per hectare utilised agricultural area (UAA)

Discussion and conclusions

Several published results on the N soil budget surplus in Germany agree well with the N surplus presented here. Leip et al. [11] calculated an N surplus of 76.1 kg N ha−1 UAA for Germany based on the agri-economic model CAPRI (years 2001–2003), while our figure for these years is 79.8 kg N ha−1 UAA. With a coupling of several high resolution models of N fluxes in Europe, Leip et al. [21] further estimated an N surplus for agricultural soils of 89.8 kg N ha−1 UAA (year 2002) for Germany. Among other things, the difference of 10 kg N ha−1 UAA to our result can be attributed to the fact that essential quantities of the N soil budget such as N mineral fertilisation and crop yield have also been model-generated by Leip et al. [21]. In contrast, we used available statistical data, as did Leip et al. [11] in the CAPRI-based approach, where a country is treated as an aggregated (single) farm. The agricultural soil surface N surplus of 71.9 kg N ha−1 UAA for Germany (year 1997) by Campling et al. [23] is almost identical to the corresponding value of 72.9 kg N ha−1 UAA (mean 1996–1998) from our study. Our values of the N soil surface budgets summed up for the districts match also very well with the figures of the national N budget which is published annually by the Federal Ministry of Agriculture [37] and which serves as reference values for spatialised N budgets in Germany. Averaged over the time period 1995–2017 the difference is less than 1 kg N ha−1 a−1 UAA, due to a harmonisation of the underlying calculation methods between the two responsible working groups.

Overall, the long-time series of the N district soil surface surplus in Germany presented here forms a reliable basis for further interpretation, mainly with regard to the following questions: (a) How has the N surplus developed over time, both in Germany as a whole and regionally? Does the N surplus show a reaction to agri-environmental policies and actions? (b) In which districts are the highest N surpluses to be found (“hot spots”)? How pronounced are the regional differences?

Obviously the time series of Germany’s N soil surface budget surplus show no decline in the past ten years (Fig. 1). Apparently, the agri-environmental action programmes over the last twenty years were not sufficient to reduce the altogether high and regionally excessive N surpluses in German agriculture. That was judged by the EU Court of Justice on 21 June 2018 (case C-543/16) declaring that the Federal Republic of Germany failed to fulfil its obligations under Council Directive 91/676/EEC concerning the protection of ground water against pollution caused by nitrates and high N loads into North Sea and Baltic Sea from agricultural sources as documented in the first Nitrate Report published in 1996 and every Nitrate Report published since then thereafter.

This development is partly caused by the expansion of biogas production, because an increasing amount of N is retained within the agricultural cycle, as the N in harvested energy crops returns (almost) completely to the UAA with the digestates and less N leaves the agricultural sector in the form of marketed crops. Additionally, since the fermentation of manure in (gas-tight) biogas plants causes lower gaseous N losses compared to storage in uncovered slurry tanks, more N is applied with the digestate compared to the equivalent amount of manure and thus the N surplus tends to increase. These effects can be determined only when biogas production is included in the N budget scheme.

The development of the N soil surface surplus in the districts over 20 years shows a distinct spatial disproportionation. In the majority of the districts which had a low initial N surplus in 1995–1997, the values have decreased, while those districts with a high initial N surplus show a clear further increase (Fig. 3c). We assume that the spatial disproportionation is a consequence of the specialisation of German agriculture: animal husbandry, in particular pork and poultry production, has been concentrated in a few regions. Contrasting, livestock production is abandoned in large parts of Germany and farmers there mainly practise arable farming. Nevertheless, comparing the changes 1995–1997 to 2015–2017 of the N soil surface surplus and NUE, respectively, with the changes of the livestock density in the districts no significant correlation can be observed (R2 = 0.06 and 0.07, respectively). The exclusion of manure transfers, however, lead to a higher correlation between the change of livestock densities and the change of N surpluses 1995–1997 to 2015–2017 (R2 = 0.47). This lead us to the conclusion, that the impact of livestock density changes on the N soil surface surplus on the district level is largely offset by manure transfers.

Given this situation, the most noticeable effect of the Fertiliser Ordinance, which implemented Council Directive 91/676/EEC in Germany (first in 1996), appears to be that the application of manure was limited to 170 kg N ha−1 UAA. As a result, from the hot spot districts significant quantities of manure are transferred to adjacent regions with low livestock density which still have capacity for additional N input in their N budget, given the restrictions by the German Fertiliser Ordinance, as shown by Fig. 4.

Compared to the N soil surface budget surplus of major European countries (Fig. 5), our results show the same temporal development as most the N soil surface budget surplus except for the United Kingdom. Since these data sets are based on different, country-specific methodologies [38] a comparison between the results is only viable with limitations, it is recognisable, that the results for Germany do not deviate strongly from those of other European countries. However, it strongly deviates from the Eurostat data set for Germany (data not shown), which shows a decline since 1995, that partly can be explained by not including biogas production into this data set.

Fig. 5
figure 5

Comparison of the N soil surface budget surplus (1995 to 2017) of selected European countries (Spain, France, Italy, Poland and United Kingdom) according to Eurostat [38] to the N soil surface budget surplus Germany according to present methodology

Irrespective of the evaluation that the N budget provides robust and meaningful results, some weaknesses and options for improvement of the N budget approach presented here are also addressed in the following.

  1. 1.

    The fluctuations in the N input and output and also in the budget surplus impair the informative value of the figures for individual years. Thus, moving averages (three-yearly or longer) of budget values should always be used as the basis for evaluation and trend analysis of the N surplus.

  2. 2.

    The surplus of the N soil surface budget is a proxy value to characterise the nitrate groundwater pollution potential averaged over the entire UAA of a district, i.e. over all soil types and climate zones, types of farming (cash crop farming, livestock farming), crop types, and cultivation methods (conventional, organic) [39, 40]. The results are therefore only suitable for large-scale comparisons and do not permit any statements with a higher spatial resolution (within the districts). This leads to an average value that probably underestimates the actual N surplus for intensively managed UAA within a district, but overestimates it for extensively used areas. As an example, the districts along the Upper Rhine Plain in Baden-Wurttemberg can be cited, each comprising in almost equal proportions cash crop farms in the Upper Rhine Plain and cattle farms in the Black Forest. These problems can partly be solved by calculating future N surpluses separately for arable land, permanent crops and grassland.

  3. 3.

    The district N surplus refers to spatial units of very different size. In 2016, the UAA in the districts in Germany (district-free cities not included) ranged from 7000 ha UAA (district Sonneberg) to about 341,000 ha UAA (district Börde). Thus, the spatial representativeness of the individual N surplus figures span over a large range.

  4. 4.

    All federal states offer agri-environmental programmes which stipulate a limitation of N fertilisation and therefore reduce the N surplus of the included UAA. However, the N fertilisation cap varies substantially between the federal states and furthermore the programmes are subscribed to varying degrees in the districts. These factors modify the district N surpluses, but are not reflected in the approach for calculating the N mineral fertilising for the districts (Eq. 4).

  5. 5.

    Organic farms do not use mineral fertilisers and therefore have typically a lower N budget surplus [41, 42]. The proportion of organically managed UAA varies considerably between districts. This factor has not yet been taken into account in the N budget approach up to date; to improve the surplus results, the N soil budget for conventionally and organically managed UAA should be calculated separately.

In Germany, mineral fertiliser sales statistics are available at the state level. The approach to derive the regional mineral fertilising amounts remains the biggest methodological challenge and at the same time the most sensitive factor for the district N budget surpluses. As long as no detailed representative statistical data on mineral fertilising at farm level are available, we consider that the method applied here is a viable approach. Other publications presenting spatialised N soil surface budgets used different approaches, e.g. Özbek and Leip [18] employed regional mineral fertiliser sales statistics for 26 NUTS 2-regions in Turkey as well as Eriksson et al. [15] for two river basins in Poland. Grizzetti et al. [16] intersected culture-, year- and country-specific mineral fertiliser application rates as published by the International Fertiliser Association (IFA) with CORINE Land Cover maps to derive specific mineral fertiliser application rates for NUTS 2-regions of 15 EU member states. Lassaletta et al. [19] applied a similar approach for the calculation of spatialised N soil surface budgets in the Ebro river basin in Spain. Gamer et al. [43] used mineral fertiliser sales data from the German Farm Accountancy Data Network (FADN) to calculate spatialised N, P, and K balances in Baden-Württemberg.

Since 2014 district-level data on the transfer of farm fertilisers are available in the nutrient reports of four federal states in Germany [44,45,46]. However, there is a lack of data on the distribution of the manure released among the districts in the federal states without nutrient reports. To improve the transparency of the manure transfer, appropriate data should be collected in the federal states concerned. The next complete FSS will be conducted in 2020. Compared to earlier surveys, it will provide more detailed information on the storage and application of farm manure, bedding materials, organic fertiliser transfer, and other items. We expect that the calculation of spatialised N soil surface budgets can be further developed and improved on the basis of these data and considering the weak points listed above.

Uncertainty and sensitivity analysis as well as the evaluation of measures to reduce N budget surpluses will be published in two separate articles.