Next Article in Journal
Agro-Morphological and Biochemical Characterization of Wild Prunus spinosa L. Subsp. dasyphylla (Schur) Domin Genotypes Naturally Grown in Western Black Sea Region of Turkey
Next Article in Special Issue
Use of Biochar-Compost for Phosphorus Availability to Maize in a Concretionary Ferric Lixisol in Northern GHANA
Previous Article in Journal
Is Tillage a Suitable Option for Weed Management in Conservation Agriculture?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Productive and Environmental Consequences of Sixteen Years of Unbalanced Fertilization with Nitrogen and Phosphorus—Trials in Poland with Oilseed Rape, Wheat, Maize and Barley

Department of Plant Nutrition and Fertilization, Institute of Soil Science and Plant Cultivation—State Research Institute, Czartoryskich 8, 24-100 Puławy, Poland
*
Author to whom correspondence should be addressed.
Agronomy 2020, 10(11), 1747; https://doi.org/10.3390/agronomy10111747
Submission received: 22 September 2020 / Revised: 4 November 2020 / Accepted: 6 November 2020 / Published: 10 November 2020

Abstract

:
Two factorial field experiments were carried out between 2003 and 2018 in the Experimental Stations in Eastern and Western Poland using four crop rotations with winter oilseed rape, winter wheat, maize and spring barley. The initial value of phosphorus (P) in Grabów soil was 69.8 mg P·kg−1 soil and in Baborówko soil it was 111.3 mg P·kg−1 soil (Egner-Riehm Double-Lactate DL). P fertilizer was added annually at 39 kg P·ha−1 under winter oilseed rape, 35 kg P·ha−1 under maize and 31 kg P·ha−1 under wheat and barley using superphosphate and nitrogen (N), which was added at five levels (30–250 kg N·ha−1) per year as ammonium nitrate in addition to controls with no added fertilizer. Through the several years of the experiment, P fertilizer had no effect on crop N use efficiency (NUE) nor crop productivity. There was significant soil P mining particularly in the high-N fertilizer trials causing a reduction in the content of available soil P by up to 35%. This work recommends that, based on soil P analysis, P fertilizer should not be added to high-P soils. This practice may continue uninterrupted for several years (16 in this case) until the excess soil P has been reduced. This mechanism of removal of “legacy” P from soil has major implications in reducing runoff P into the Baltic Sea drainage area and other water bodies.

1. Introduction

Nitrogen (N) and phosphorus (P) are the two most essential nutrients ensuring food production. Both N and P are involved in vital plant functions such as photosynthesis protein formation, and symbiotic N fixation [1]. Simultaneously, N is the most important crop yield-limiting factor in the world [2]. As a main nutrient element, it is needed in a relatively large quantity for the production of proteins, nucleic acids and chlorophyll in plants [3]. The natural sources of N and P are very different: in principle, N availability is unlimited as an atmospheric gas, whereas P comes from rock phosphate, renewed with the uplift of continental rock [4]. In ecosystems, N and P are bound within waste organic products and dead organic matter and must be remineralized to release inorganic orthophosphates or dissolved and reduced to inorganic nitrate and ammonia before either element can be absorbed by plants.
However, since the 1960s, anthropogenic activities have dramatically modified the cycle of N and P in the biosphere, threatening natural ecosystems [4]. The invention of the Haber–Bosch process more than 100 years ago, allowed for the production of synthetic mineral fertilizers from ammonia on a large scale [5]. Thus, European N fertilizer use has grown from about 1.2 Tg in 1961 to 11.8 Tg in 2018 [6]. The increased use of N fertilizers has contributed greatly to the increased global food, feed and biofuel production [7]. Nevertheless, the great availability of N in agriculture increases losses to the wider environment, which can lead to problems related to human health and ecosystem degradation [8,9,10]. The volatilization of ammonia, leaching nitrate, and the emissions of di-nitrogen, nitrous oxide and nitrogen oxide are the main loss pathways of agricultural systems [11].
Numerous studies show that less than 50% N and 25% P introduced into agroecosystems in the form of mineral and natural fertilizers are effectively used while the rest is dispersed in the environment, and contributes to various negative ecological effects [12,13,14]. Nitrogen use efficiency (NUE) is considered quite low on average in conventional agricultural systems in the world, including developed countries. The NUE for world cereal production has been estimated for 33% of fertilizer N recovered by the crop [15].
The excessive eutrophication of waters, especially in marine coastal areas and estuaries, has been a growing problem worldwide for many decades. It is caused both on a global and local scale by the man-made intensification of biogeochemical cycles of N and P, which increases the loads of nutrients, introduced into the water from point and diffuse sources [16,17,18,19,20,21]. The main streams of nutrients, apart from natural sources, come from food consumption (point sources) and agriculture (diffuse nutrient outflow), which are broadly understood, where the significant increase in the N and P load is related to unbalanced mineral fertilization and the intensification of animal production [13,21].
To reduce eutrophication in the Baltic Sea, the HELCOM Contracting Parties agreed in 2007 and 2013 on the Baltic Sea Action Plan (BSAP), one of the main elements of which is the Nutrient Input Reduction Scheme, where the reduction targets are divided among the coastal countries [22,23]. The revised (2013) waterborne inputs for the entire Baltic Sea were determined as annual reductions of 89.26 t N and 14.37 t P [21,24]. The ecological state of the Baltic Sea is systematically improving as a result of international efforts, but according to HELCOM, excessive eutrophication still occurs in most parts of the basin [25].
The risk of P loss increases with its concentration in soil [26] and within 21 days of a fresh application of highly water-soluble P fertilizers [27]. The addition of P as reactive phosphate rocks has been shown to decrease P losses compared to superphosphate by 10%–80% at the plot, field and catchment scale [28], which is caused by a combination of low water solubility and a greater density of fertilizer particles, reducing the availability of P loss [29].
On the other hand, the mining of phosphate rocks for use in agricultural fertilizers has increased dramatically in the latter half on the past century. In the worst-case scenarios, about 40%–60% of the current resource base will be extracted by 2100 [30]; under the best estimates, the depletion would be around 20%–35% in the time. The environmental implications and the finite character of rock phosphate make it necessary to reconsider the views on this element in crop production.
Recent studies have suggested that the amount of reactive N compounds emitted into the environment is far too high and already exceeds the “safe operating space for humanity”. More appropriate management of N use is therefore of key importance, but it should be emphasized that the proper management of the element is linked to the balanced fertilization with other nutrients to achieve both economic and environmental goals [31,32]. The essential, irreplaceable role of P in plant nutrition requires that sufficient P is available in soil to meet the need of the growing crop in order to produce optimum, economic yields. The need varies through the life cycle of the plant and is usually greatest in the early stages of growth, when the root system is not fully developed, to explore the maximum volume of soil [33,34]. It has been found that, under insufficient P availability, crops do not respond to N fertilization [35,36]. There is evidence that high yielding crop varieties need more N and P to realize their potential and thus under insufficient fertilization with P the reserve of P in soil might be exhausted and lead to yield reduction. Nevertheless, an important issue is to recognize the possibility of reduced P fertilizer use on soils with a high concentration of available forms of the element.
This paper shows that the amount of P applied in agriculture could be considerably reduced by improving the fertilizer recommendation. In the data collected between 2003 and 2018, including four complete rotations of the most popular crops in Poland (winter oilseed rape, winter wheat, maize for corn and spring barley) fertilized with increasing N rates, the impact of long-term P cessation in the soil on crop productivity and changes of the available forms of P have been described.
The following hypothesis has been set: the cessation of P fertilization on soils with a very high and high available P content does not impact crop N use efficiency, and neither does it cause excessive soil P depletion in the long term.

2. Materials and Methods

The long-term field experiments were carried out between 2003 and 2018 at the Experimental Stations of the Institute of Soil Science and Plant Cultivation in Grabów (21°39′ E, 51°21′ N) and Baborówko (16°37′ E, 52°37′ N), Poland. Both Experimental Stations are located at the cold temperate dry climate zone in which potential evapotranspiration normally exceeds rainfall during the whole vegetation period. Average annual rainfall at Baborówko was 492 mm and at Grabów it was 587 mm during the period of investigation. In Experimental Station Grabów, soil under the experiment was heterogeneous—sandy loam (World Reference Base WRB: Stagnic Luvisols). In Baborówko the experiment was localized partly on the sandy loam (WRB: Albic Luvisol) and partly on the black earth (WRB: Gleyic Phaeozem) [37]. In two factorial, two replication experiments, four crops were grown each year in the rotation of: winter oilseed rape (Brassica napus L. var. oleifersa); winter wheat (Triticum aestivum L.); maize (Zea mays L.); spring barley (Hordeum vulgare L.). In Baborówko during 2015–2018 and in Grabów in 2015 and 2018, mustard (Sinapis alba L.) was grown instead of winter oilseed rape because of the extremely unfavorable conditions for the emergence or wintering of the crop.
During the experiment, no manure was applied, no leguminous plants were grown and straw was removed from the field. In the split-plot layout, the first variable was P fertilizer in the P-plus and P-minus (control) treatments and the second variable was six levels of N fertilizer, including a control with no N supply. In Grabów, the soil was characterized as slightly acid (pHKCl 6.2) (according to the classification of soil reaction used in Poland, based on the measurement of pH in 1M KCl 1:5, soil:solution), had a medium content of available potassium (88.8 mg·kg−1 soil) (Egner-Riehm Double-Lactate DL method) and a high content of available magnesium (30.4 mg Mg·kg−1 soil) (Schachthabel method). In Baborówko, the soil reaction was neutral (pHKCl 6.8), the content of potassium was high (116 mg K·kg−1 soil) and the concentration of magnesium was very high (54 mg Mg·kg−1 soil).
In Grabów, the initial content of available P in 2003 was recognized as high (69.8 mg P·kg−1 soil) and in Baborówko as very high (111.3 mg P·kg−1 soil), according to the classes of available P content in soil used in the Polish recommendation system. Since 2003, in the P-plus treatment, P fertilizer as a superphosphate was applied before sowing crops at the following rates: 39 kg P·ha−1·y−1 under winter oilseed rape, 35 kg P·ha−1·y−1 under maize and 31 kg P·ha−1·y−1·under wheat and barley. In the control, no P fertilizer was added past 2003. Potassium (66.0–116 kg K·ha−1·y−1) and magnesium (42 kg Mg·ha−1·y−1) were constant in the P-plus and in the control. Nitrogen fertilizers were applied as ammonium nitrate. Spring barley was fertilized with 0, 30, 60, 90, 120, 150 kg N·ha−1·y−1, winter wheat with 0, 40, 80, 120, 160, 200 kg N·ha−1·y−1, and winter oilseed rape and maize with 0, 50, 100, 150, 200 and 250 kg N·ha−1·y−1. For spring barley, the first dose of 30 kg N·ha−1·y−1 was applied at sowing, and the successive N doses were carried out over two-week intervals. For winter wheat, the first dose of 40 kg N·ha−1·y−1 was applied at the beginning of the spring vegetation, and the successive N doses were applied as for barley. For maize and winter oilseed rape, the total rates of 50 and 100 kg N·ha−1·y−1 were applied after emergence and in the early spring, respectively. In the N150-N250 treatments, the first rate of 100 kg N·ha−1·y−1 was applied after emergence and in the early spring, similarly, and the subsequent additions of 50 kg N·ha−1·y−1 were spread over 14-day intervals.
Soil samples from the plough layer were collected each year between 2007 and 2018 after the harvest of barley. The samples were dried, passed through a 2-mm sieve and analyzed for available P and potassium according to the Egner-Riehm DL method [38] with UV–Vis spectrophotometry (Nicolet Evolution 300), and for available magnesium according to the Schachtschabel method [39] with F-AAS spectroscopy (Varian Spectra AA-240 FS).
The grain and straw yields of winter oilseed rape, winter wheat and spring barley were harvested by a plot harvester; maize was harvested by hand. Wheat and barley were harvested over an area of 17.9 m2 and rape and maize over 19.1 and 20.3 m2, respectively. The grain yields of wheat, barley and corn were determined at a moisture level of 15%, and winter oilseed rape was determined at 12%. The plant samples of grain and straw for each crop were collected at full maturity from an area of 1 m2 and analyzed for total N and P as determined by the mineralization of the sample using the Kjeldahl method with Continuous Flow Analysis (CFA) and UV–Vis spectrophotometry (Skalar San+). The uptake of N and P in grain and straw was calculated from the percentage of the elements and the yields of the main and fore crop.
Following the method proposed by the EU Nitrogen Expert Panel (NEP), N use efficiency (NUE) was calculated for the crops to determine the effect of plants’ P supply on this parameter. According to the approach, NUE calculations based on N input and N output at different levels provide information about resource use efficiency, the economy of food production (N in harvested yield) and the pressure on the environment (N surplus). The desired value of NUE, N uptake (Yn) and N surplus (Ns) should be in the range NUE 50%–90%, Yn > 80 kg N·ha−1 and Ns < 80 kg N·ha−1 [40]. Nitrogen use efficiency was calculated according to the formula [40]:
NUE = Yn F × 100
where Yn is total N uptake by the crop (kg·ha−1); F is N fertilizer rate (kg·ha−1).
N surplus (Ns) was calculated according to the formula [40]:
Ns = F−Yn
P use efficiency (PUE) (%) by the crops in the P-plus treatment was calculated according to the formula [41]:
PUE ( % ) = Pn plus Pn control PF 100
where Pn plus is total P uptake by plants fertilized with P (kg P·ha−1), Pn control is total P uptake by plants without added P fertilizer (kg P·ha−1) and PF is the P fertilizer rate (kg·ha−1).
P surplus (Ps) was calculated according to the formula [41]:
Ps = PF−Pn
The data, derived over the sixteen years, were processed using an analysis of variance (ANOVA). Tukey’s test was applied to evaluate the significance of differences between the treatments. Statistical processing of the results were performed using the Statgraphics Centurion 18 Package (Statgraphics Centurion, Rockville, MA, USA).

3. Results

3.1. Effect of P Additions on Crop Productivity

The grain yields derived from the 2003–2018 P trials averaged over the six levels of N fertilizer are presented in Figure 1, Figure 2, Figure 3 and Figure 4. The general trend is that there were no significant increases in the growth of these crops due to added P. In general, higher yields were obtained in Grabów than in Baborówko. Within the analyzed period, comprising several years of experiments, there was found to be significant soil P mining by the crops in the controls where P fertilizer was not added.
The yield of winter oilseed rape oscillated between 0.85 and 4.68 t·ha−1 in Grabów, and in Baborówko it ranged from 1.34 to 3.15 t·ha−1 (Figure 1). In Grabów, the average yields obtained in the P-plus treatment were 2.77 and 2.63 t·ha−1 in the control. In 2017, when the crop was grown for the last time, a yield reduction of 0.85 t·ha−1 was observed in the treatment where P fertilization had ceased (Figure 1a). In Baborówko, the seed yields were comparable between the treatment and control achieving an average of 2.1 t·ha−1 (Figure 1b).
In Grabów, the average grain yields of winter wheat were 5.32 t·ha−1 in P-plus and 5.10 t·ha−1 in the control. In Baborówko, the mean grain yields were 4.40 t·ha−1 for both P-plus and the control (Figure 2).
In Grabów, the yields of maize showed a wide range from 1.9 to 12.4 t·ha−1. During 2015–2018, the crop yield was slightly lower in the control compared to the P-plus treatment. (Figure 3a). In the same years, maize in Baborówko yielded levels of 8.80 t·ha−1 and no differences between P-plus and the control were noted (Figure 3b).
The yields of spring barley in Grabów were more stable compared to the other crops and ranged from 2.5 to 5.1 t·ha−1 (Figure 4a). Similar to maize in this location during 2014–2018, barley growth was lower in the control. In Baborówko, the yields of barley were variable, with no P fertilizer effect (Figure 4b).

3.2. Crop Response to the Combined N and P Fertilizer Additions

The average grain yields for 2003–2016 of winter oilseed rape, winter wheat, maize and spring barley are presented in Figure 5, Figure 6, Figure 7 and Figure 8. The crop response to added N showed clear positive trends. In Grabów, for both P treatments, the increased yields of winter oilseed rape were statistically significant at 150 kg N·ha−1, winter wheat at 120 kg N·ha−1, maize at 100 kg N·ha−1 and barley at 120 kg N·ha−1. In Baborówko, crop growth increased with increasing N fertilizer additions, except for maize, which did not respond to N fertilization above 150 kg N·ha−1, regardless of the level of added P.
In Grabów, a tendency for grain yield reduction together with N rates was found, in particular for maize; however, the interaction between the two experimental factors has not been statistically proven. The average grain yields of maize in the P-plus treatment ranged from 5.90 to 8.36 t·ha−1 and in the P control between 5.40 and 8.07 t·ha−1, depending on the level of N fertilization (Figure 7). There was no similar relationship observed in Baborówko, where the average yields achieved in the analysis period were practically the same, regardless of the P supply.

3.3. Effects of N Additions on P Uptake and Use Efficiency

The following indices characterizing P use efficiency were calculated: P uptake (Pn), P surplus (Ps) and P use efficiency (PUE) by plants fertilized with P. The effect of N additions on these indices are presented in Table 1, Table 2, Table 3 and Table 4. In Grabów there were significant effects of P additions on P uptake, but no additional effect of N fertilization was seen. In Baborówko, P uptake by the crops was not affected by P additions.
The P uptake by winter oilseed rape, wheat, maize and barley increased together with N rates. The average P uptake by winter oilseed rape was 21 kg P·ha−1 in Grabów and 15 kg P·ha−1 in Baborówko (Table 1). The uptake of P in wheat reached 20 and 17 kg P·ha−1 (Table 2), in maize 29 and 24 kg P·ha−1 (Table 3) and in barley 20 and 15 kg P·ha−1 (Table 4), respectively.
In the treatment with P fertilizer, the surplus P calculated for the crops was positive over the whole range of N doses. The highest value of surplus P was found in the winter oilseed rape trial (18 kg P·ha−1 in Grabów and 24 kg P·ha−1 in Baborówko, on average) and the lowest was in the maize trial (7.0 kg P·ha−1 in Grabów and 10 kg P·ha−1 in Baborówko). Where no P fertilizer was added since 2003, the surplus P was negative with increasing negative values governed by the rate of N added. The most negative P surplus was obtained in the maize trials in both locations (−27 kg P·ha−1 in Grabów and −23 kg P·ha−1 in Baborówko on average) as a consequence of the enhanced P uptake by the crop. The P use efficiency (PUE) in the experiments was low, averaging 3%–7%.

3.4. Effects of P Additions on Nitrogen Uptake and Use Efficiency

The N use indices calculated for winter oilseed rape, winter wheat, maize and spring barley for each location are presented in Table 5, Table 6, Table 7 and Table 8.
Both N use efficiency and N uptake by all the crops exhibited higher values in Grabów compared to Baborówko. This left a greater N surplus in Baborówko (Table 5, Table 6, Table 7 and Table 8). Although NUE and Yn in the control were slightly lower than the P-plus treatment, the impact of P management on these indices was not statistically significant. Consistently, withholding P additions in the control did not significantly affect the N surplus.
Maize exhibited the highest values of NUE in the range of 71%−231%, with an average of 127%. The average N uptake by maize in both locations was greater than the minimal EU NEP [40] value of 80 kg N·ha−1 (Table 7). The crop removed the greatest quantity of soil N in the range of 50–150 kg N·ha−1 applied in the experiment (43 kg·ha−1 on average). Over the entire range of N fertilizer rates, winter oilseed rape showed the lowest indices of NUE (46%−145%), N uptake (50–155 kg N·ha−1) and there about 53 kg N·ha−1 remained in the soil (Table 5). The average value of NUE calculated for winter wheat was 110% on average, and the mean N uptake was 114 kg N·ha−1. In Grabów, the NUE for barley was higher than Baborówko and oscillated over the ranges 88–211 kg N·ha−1 and 66–135 kg N·ha−1, respectively. The uptake of N by barley in Baborówko was smaller, thus the average value of surplus N in this location in the control P treatment was positive.

3.5. Changes in Available P in the Soil

In Grabów, the average content of available P in soil in the P-plus treatment increased from the initial 2003 value of 69.8 to 90.3 mg P·kg−1 soil between 2007 and 2018. In Baborówko, a similar increase was observed—from 111.3 to 134.4 mg P·kg−1 soil. Through the process of P soil mining in the control treatment (no P added), the average concentration dropped to 63.3 mg P·kg−1 soil in Grabów and to 107.8 mg P·kg−1 soil in Baborówko. An analysis of variance showed a significant reduction in available P in soil at increasing rates of N fertilizer addition in Grabów as a result of the higher crop yields and P uptake. In Grabów, soil P levels decreased in both the P-plus and control as N doses increased reaching as low as 51.5 mg P·kg−1 soil (Figure 9a). In contrast, this trend was not seen in Baborówko (Figure 9b) which had much higher soil P levels than Grabów, even though the P content of the soil was also about 20% higher in the P-plus treatment than the control.

4. Discussion

The 16-year experiment, the results of which have been presented here, was conducted on soils with high (Grabów) and very high (Baborówko) contents of available P in the plough layer. Such soils should be considered as typical for Poland, as they account for almost 60% of the country. The achieved soil acidity in Grabów was slightly acid (pH KCl 6.2) and neutral in Baborówko (pH KCl 6.8). The Polish recommendation system for P fertilizer application is based on soil P supply and on crop P demand. Fertilizer recommendations using five classes of available P (determined by the Egner-Riehm DL method in topsoil) have been identified: very low, low, medium, high and very high. It is recommended to achieve at least a medium level of soil P (44.1–65.4 mg·kg−1 soil). In the experiments, the soil was characterized by a low content of soil organic carbon (5.0 g org C·kg−1 soil) and no exogenous organic matter has been applied for over 30 years.
The results of the experiments confirm the hypothesis that several years without P fertilizer addition did not exhaust the available soil P to levels that limit crop productivity. In Grabów, where the initial value of available soil P was 69.8 mg P·kg−1 soil, annual applications of P fertilizer led to an increase in available soil P over the whole range of N fertilizer rates, although the effect was suppressed at higher N levels. The P level in the soil remained within the high class of P content.
After P fertilizer addition was stopped in 2003, the topsoil P in the controls significantly dropped. It should be noted that this decrease in soil P was further enhanced by increases in N fertilizer rates. At no or low additions of N fertilizer, N0-N1 (0–50 kg N·ha−1), the content of available P in soil was close to the initial value of 69.8 mg P·kg−1 soil. In the range of medium rates, N2-N3 (60–120 kg N·kg−1) the P concentration was 63.2 on average. At the highest N4-N5 rates (120–200 kg N·kg−1), the content dropped further to 57.5 and 52.4 mg P·kg−1 soil, respectively. Going beyond the highest N5 rate could result in P limitation of crop productivity.
In Baborówko, where the experiment was located on soil with a very high content of available soil P (111.8 mg P·kg−1), the soil P content remained high and no interaction between the added N and P fertilizer was noted. Similarly, the P-plus treatment only increased available soil P to 133 mg P·kg−1. Crop yields were lower than in Grabów, and the lowest values of P uptake were noted. Hence, the P balance in Baborówko showed fewer negative values, and the N additions did not affect the soil P content.
Phosphorus in soil occurs in inorganic (Pi) and organic (Po) forms, which undergo various transformations, affecting their fate in soil and their bioavailability to plants [42,43]. The share of Pi in the total P (Pt) pool is estimated at 35% to 70% [44]. The soil solution contains only a small fraction of orthophosphate ions available to plants, which during plant growth are replenished from mobile reserves of Pi and organic matter. The remainder of Pt is the so-called “spare” P, the release of which is very slow [45,46]. Nevertheless, P is not irreversibly fixed in soil and is frequently mobilized during the growing period [47]. The mineralization of Po occurs under the influence of enzymes secreted by bacteria, fungi and plant roots. The rate of mineralization is influenced by soil environmental conditions, primarily pH and temperature [48,49]. In acid soils, P is fixed by free oxides and hydroxides of aluminum and iron. In alkaline soils it is greatly adsorbed onto lime surfaces and clay matter and precipitated in the form of different Ca-P minerals. All these fixation reactions contribute to reducing the solubility and plant availability of applied P [50,51]. In our experiments with the relatively high levels of soil P, plant availability was found to be optimal and had no constraints due to this soil parameter for P uptake being observed.
In Grabów, the P additions resulted in lower crop uptake, particularly at the highest N fertilizer rates. In Baborówko, where all the crops had lower yields, there were no significant differences between the P-plus and control treatments. The sixteen years of soil P mining in the zero-added P controls caused a negative P surplus. The average P surplus calculated for one rotation in Grabów was −82 kg P·ha−1 and it was −68 kg P·ha−1, in Baborówko. In effect then, over the sixteen years (four rotations), the capacity to reduce soil P was −328 and −272 kg P·ha−1, respectively. Nevertheless, such a negative balance has not affected the crop productivity in Baborówko, even at the highest rates of applied N fertilizer. This is a significant finding for managing P levels in agriculture soil. Soil mining by crops receiving only N fertilizer could have a major impact on reducing legacy P (from many decades of overfertilization) in the Baltic Sea drainage basin.
As reported in the literature, the P fertilizer inputs in the world are two to five times greater than the amount exported in the final products [52]. It was thought to be of great importance in the early years in the mid-19th century to apply more P than was removed in the harvested crop to achieve an acceptable yield. This concept led to the accumulation of large quantities of surplus P fixed in agriculture soils [53]. Johnson and Syers [54] found in experiments conducted between 1856 and 1901 in Rothamsted, that P fertilizers were applied year to year and the resulting P balance was positive. However, following 100 years after the last P application, residual P was still being recovered in winter wheat compared to treatments where additional P had been applied since 1856. Phosphorus accumulation in soils typically occurs when P fertilizer is added, whilst P depletion can occur under intensive cropping unless there is an external supply [55]. Sun et al. [56] reported that over a 24-year experiment with winter wheat and maize, despite P depletion, yields were maintained, albeit at low values.
In our experiments, the rate at which P has been released in soil each year has been sufficient for several years to provide enough P to produce the optimum yield for the modern high-demand for both N and P cultivars. The uptake of P by the crops was comparable between the P-plus and control treatments. The P use efficiency calculated as the difference between P uptake in plants fertilized with P and not fertilized with P was low and achieved at most 12%.
The data presented by Lasaletta et al. [57] show that currently only 47% of the reactive N added globally onto cropland is converted into harvested products, compared to 68% in the early 1960s. This means, that more than half of the N used in crop fertilizer is currently lost to the environment. The prime cause of N loss is through nitrate leaching or denitrification in connection with excessive rainfall. When the whole dose of N fertilizer is applied prior to planting, a significant risk of losses during the spring occurs. The time between N application and its active absorption by the crops, provides numerous opportunities for N leaching, clay fixation, immobilization, denitrification and volatilization [58]. Developing fertilizer strategies that increase NUE could reduce unnecessary input costs to farmers and environmental impact of N losses. Examples include adjusting N fertilizer additions based on knowledge about the soil mineral forms of N, phasing N fertilizer applications or using fertilizers with a nitrification inhibitor [59].
Nitrogen use efficiency is suggested to be a key indicator for agricultural systems. Several studies have estimated NUE in cropping systems, but these studies often use different definitions, input data and assumptions [60]. In this paper, NUE indices were calculated according to the method proposed by the EU NEP [40], using three related indicators: NUE, N surplus and N uptake. The impact of P additions on NUE indices for the most common crops in Poland was examined. Contrary to what might be expected for such a long period, the effect of increased soil P mining (by the crops) on the NUE indices was negligible. This means that P mining did not increase the risk of N losses due to insufficient N uptake by crops. Higher values of NUE were noted in Grabów compared to Baborówko due to a greater crop yield. Nitrogen uptake was enhanced in Grabów due to the more favorable rainfall conditions through the vegetation periods. Maize exhibited the highest values of NUE and N uptake and similarly the crop gave the lowest N surplus (average minus 7 kg N per ha). Over the entire range of N fertilizer additions, winter oilseed rape showed the lowest NUE and Yn values and left a soil surplus of about 50 kg N per ha per year. The cultivation of winter oilseed rape is problematic for the reason of considerable N surplus [61,62]; simultaneously the crop requires a high quantity of P for optimum growth [63,64] and P deficiency restricts both top and root growth [65,66]. In the presented experiments, surplus N has not been affected by long-term P fertilizer additions, although in Baborówko, under the rates of 200 and 250 kg N·ha−1, the surplus exceeded the desired value recommended by the EU NEP.
According to the results obtained by Quemada et al. [67], based on the data collected through surveys of 195 arable farms in Europe, with the prevalence of cereals grown in 2–3 crop rotations, the mean value of NUE was 60%. This value was obtained for a mean N input of 176 kg N·ha−1 and 105 kg N·ha−1 output, resulting in an N surplus of ca. 72 kg N·ha−1. In our experiments with inputs close to 170 kg N·ha−1—i.e., 150 kg N·ha−1 for winter oilseed rape, maize and barley—the mean NUE was 76%, N uptake was 120 kg N·ha−1 and N surplus was 30 kg N·ha−1. Similarly, at an input rate of 160 kg N·ha−1 for barley, the indices were NUE 84%, Yn 134 kg N·ha−1 and Ns 25 kg N·ha−1. The relatively high values of NUE and Yn in our experiments might be the result of the split application of N fertilizers.
In Poland, the average N fertilizer rates in the production of cereals over recent years do not exceed 90 kg N·ha−1 per year and for maize and winter oilseed rape do not exceed 145 kg N·ha−1·year−1 [68]. Considering that the annual consumption of mineral fertilizers between 2010 and 2018 oscillated around 130 kg NPK·ha−1 (of which 75 kg was N [69]) on soil with an average content of 73.6 mg P·ha−1 prevailing in the country, it might be recommended to forsake year-to-year P applications to avoid its further accumulation in the soil. Nevertheless, it is important to gauge the addition of P fertilizers to the considered N fertilizer rates since the high N uptake intensifies soil P mining/depletion. It should be emphasized, however, that these recommendations concern soils with optimal conditions and where acidity does not limit P uptake. Moreover, soil testing for available P in topsoil is critical for assessing the environmental consequences of unbalanced N and P applications. This highlights the importance of soil testing prior to making decisions on fertilizer use [70,71].

5. Conclusions

Long-term N and P fertilizer experiments over 16 years with winter oilseed rape, winter wheat, maize and spring barley revealed the presence of significant soil P mining causing a reduction in the content of available forms of P but without negative impacts on crop productivity. In the Grabów region, where the initial value of available P (Egner-Riehm) was classified as high, 69.8 mg P·kg−1 soil, the sixteen years of P soil mining led to a decrease in available topsoil P to an average of 63.2 mg P·kg−1 soil. This dropped the soil P classification level to medium. Additionally, there was a close relation between the intensity of soil P depletion and the range of applied N fertilizer rates. The smaller N doses of 30–50 kg N·ha−1 did not change the soil P levels over the long term but the highest rates, 150–250 kg N·ha−1, reduced available soil P by about 35%. This points to the need for monitoring P levels in soil, especially when intensified N fertilizer applications are considered. In contrast, in the Baborówko region, where soil P content was very high (111.2 mg P·kg−1 soil), the added N fertilizer did not affect the changes in soil P over the long term. Moreover, the year-to-year addition of P fertilizer resulted in an unnecessary accumulation of available P in the soil.
Increased grain yields, higher N use efficiency (NUE) and higher N uptake by all the crops were noted in Grabów but not in Baborówko where soil P levels were much higher to start with. At the relatively high levels of available soil P in these experiments, additional P fertilizer did not affect the NUE indices. Nitrogen use efficiency for each crop dropped systematically with increased N doses and the mean NUE was 127% for maize, 111% for winter wheat and 110% for winter oilseed rape and barley but P fertilizer additions had no impact. Therefore, due to economic and environmental concerns it can be recommended to refrain from year-to-year P fertilizer applications to soils that already have high or very high P soil content. This can be continued for several years (at least 16 in this case). At the same time, this deprivation of P fertilizer will have no significant negative impacts on crop yields. Still, over the long term, P fertilizer regimes should be closely linked to N fertilizer strategies to retain optimal conditions for productivity. It is also very important that the negative P balance surplus at the field level significantly increases the utilization of P from soil resources by plants and reduces the phosphorus pool potentially at risk of leaching to ground and surface waters. Limiting phosphorus fertilization is a measure recommended by HELCOM and aimed at controlling the P load introduced in the Baltic Sea. In the long-term perspective, it could contribute to the achievement of the reduction targets set out in the Baltic Sea Action Plan. The implications of this work on reducing soil P surpluses and runoff losses to the Baltic Sea drainage basin and other receiving water systems around the world are significant.

Author Contributions

Conceptualization, A.R.; Methodology, A.R.; Validation, A.R.; Formal analysis, A.R.; Investigation, A.R.; Resources, A.R.; Data curation, A.R.; Writing—Original Draft Preparation, A.R., P.S.; Writing—Review and Editing, A.R., P.S.; Visualization, A.R., P.S.; Supervision, A.R.; All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Marschner, H.; Kirkby, E.A.; Cakmak, I. Effects of mineral nutrition status on shoot-root partitioning of photo assimilates and cycling of mineral nutrients. J. Exp. Bot. 1996, 47, 1255–1263. [Google Scholar] [CrossRef] [PubMed]
  2. Muller, N.B.; Gerber, J.S.; Johnson, M.; Ray, D.K.; Ramankuty, N.; Foley, J.A. Closing yield gaps through nutrient and water management. Nature 2012, 490, 254–257. [Google Scholar] [CrossRef] [PubMed]
  3. Fernandez, M.S.; Rosiello, R.O.P. Mineral nitrogen in plant physiology and plant nutrition. Crit. Rev. Plant Sci. 1995, 14, 111–148. [Google Scholar] [CrossRef]
  4. Elser, J.J.; Bracken, M.E.S.; Cleland, E.E.; Gruner, D.S.; Harpole, W.S.; Hillebrand, H.J.; Ngai, T.; Seabloom, E.W.; Shurin, J.B.; Smith, J.E. Global analysis of nitrogen and phosphorus limitation of primary producers in freshwater marine and terrestrial ecosystems. Ecol. Lett. 2007, 10, 1135–1142. [Google Scholar] [CrossRef] [Green Version]
  5. Vojvodic, A.; Medford, A.J.; Studt, F.; Pedersen, F.A.; Khan, T.S.; Bligaard, T.; Nørskov, J.K. Exploring the limits: A low pressure, low-temperature Haber-Bosch process. Chem. Phys. Lett. 2014, 598, 108–112. [Google Scholar] [CrossRef]
  6. Fertilizers Europe. Industry Fact and Figures 2019. Brussels, Belgium. 2019. Available online: https://www.fertilizerseurope.com/wp-content/uploads/2019/07/Industry-Facts-and-Figures-2019-Digital-version.pdf (accessed on 14 October 2020).
  7. Smill, V. Enriching the Earth: Fritz Haber, Carl Bosch, and the Transformation of World Wood Production; MIT Press: Cambridge, MA, USA, 2004; p. 360. [Google Scholar]
  8. Galloway, J.N.; Townsend, A.R.; Erisman, J.W.; Bekunda, M.; Cai, Z.; Freney, J.R.; Martinelli, S.P.; Seitzinger, L.A.; Sutton, M.A. Transformation of the nitrogen cycle: Recent trends, questions, and potential solutions. Science 2008, 320, 889–892. [Google Scholar] [CrossRef] [Green Version]
  9. Sutton, M.; Howard, C.M.; Erisman, J.W.; Billen, G.; Bleeker, A.; Grennfelt, P.; van Grinsven, H.; Grizzetti, B. (Eds.) The European Nitrogen Assessment: Sources, Effects and Policy Perspectives; Cambridge University Press: Cambridge, UK, 2011; pp. 1–607. [Google Scholar]
  10. Jones, L.; Provins, A.; Holland, M.; Mills, G.; Hayes, F.; Emmetta, B.; Halla, J.; Sheppard, L.; Smith, R.; Sutton, M.; et al. A review and application of the evidence for nitrogen impacts on ecosystem services. Ecosyst. Serv. 2014, 7, 76–88. [Google Scholar] [CrossRef] [Green Version]
  11. Bouwman, L.; Goldewijk, K.K.; Van Der Hoek, K.W.; Beusen, A.H.W.; Van Vuuren, D.P.; Willems, J.; Rufino, M.C.; Stehfest, E. Exploring global changes in nitrogen and phosphorus cycles in agriculture induced by livestock production over the 1900–2050 period. Proc. Natl. Acad. Sci. USA 2013, 110, 20882–20887. [Google Scholar] [CrossRef] [Green Version]
  12. Erisman, J.W.; Galloway, J.N.; Seitzinger, S.; Bleeker, A.; Dise, N.B.; Petrescu, A.M.R.; Leach, A.M.; Vries, W. Consequences of human modification of the global nitrogen cycle. Phil. Trans. R. Soc. B 2013, 368, 20130165. [Google Scholar] [CrossRef] [Green Version]
  13. Howarth, R.W.; Billen, G.; Swaney, D.; Townsend, A.; Jaworski, N.; Lajtha, K.; Downing, J.A.; Elmgren, R.; Caraco, N.; Jordan, T.; et al. Regional nitrogen budgets and riverine N & P fluxes for the drainages to the North Atlantic Ocean: Natural and human influences. In Nitrogen Cycling in the North Atlantic Ocean and Its Watersheds; Howarth, R.W., Ed.; Springer: Dordrecht, Germany, 1996; pp. 75–139. [Google Scholar]
  14. Roberts, T.L.; Johnston, A.E. Phosphorus use efficiency and management in agriculture. Resour. Conserv. Recycl. 2015, 105, 275–281. [Google Scholar] [CrossRef]
  15. Raun, W.R.; Johnson, G.V. Improving nitrogen use efficiency for cereal production. Agron. J. 1999, 91, 357–363. [Google Scholar]
  16. Elser, J.; Bennett, E. A broken biogeochemical cycle. Nature 2011, 478, 29–31. [Google Scholar] [PubMed]
  17. Howarth, R.W. Coastal nitrogen pollution: A review of sources and trends globally and regionally. Harmful Algae 2008, 8, 14–20. [Google Scholar]
  18. Howarth, R.; Swaney, D.; Billen, G.; Garnier, J.; Hong, B.; Humborg, C.; Johnes, P.; Mörth, C.M.; Marino, R. Nitrogen fluxes from the landscape are controlled by net anthropogenic nitrogen inputs and by climate. Front. Ecol. Environ. 2012, 10, 37–43. [Google Scholar]
  19. Ilnicki, P. Emissions of nitrogen and phosphorus into rivers from agricultural land—Selected controversial issues. J. Water Land Dev. 2014, 23, 31–39. [Google Scholar]
  20. Kroeze, C.; Hofstra, N.; Ivens, W.; Lohr, A.; Strokal, M.; van Wijnen, J. The links between global carbon, water and nutrient cycles in an urbanizing world—The case of coastal eutrophication. Curr. Opin. Environ. Sustain. 2014, 5, 566–572. [Google Scholar]
  21. Pastuszak, M.; Kowalkowski, T.; Kopiński, J.; Doroszewski, A.; Jurga, B.; Buszewski, B. Long-term changes in nitrogen and phosphorus emission into the Vistula and Oder catchments (Poland)—Modelling (MONERIS) studies. Environ. Sci. Pollut. Res. 2018, 25, 29734–29751. [Google Scholar]
  22. HELCOM. Helcom Baltic Sea Action Plan Adopted on 15 November 2007 Krakow, Poland by the HELCOM Extraordinary Ministerial Meeting, Helcom, Helsinki. 2007. Available online: https://helcom.fi/media/documents/BSAP_Final.pdf (accessed on 21 September 2020).
  23. HELCOM. Copenhagen Ministerial Declaration. Taking Further Action to Implement the Baltic Sea Action Plan—Reaching Good Environmental Status for a Healthy Baltic Sea. 2013. Available online: http://www.helcom.fi/Documents/Ministerial2013 (accessed on 21 September 2020).
  24. HELCOM. Summary Report on Development of Revised Maximum Allowable Inputs (MAI) and Updated Country Allocated Reduction Targets (CART) of the Baltic Sea Action: Environmental Science; Springer: Berlin/Heidelberg, Germany, 2013; pp. 1–261. [Google Scholar]
  25. HELCOM. State of the Baltic Sea. Second HELCOM Holistic Assessment 2011–2016; Springer: Helsinki, Finland, 2018; pp. 1–155. [Google Scholar]
  26. Curran Cournane, F.; McDowell, R.W.; Condron, L.M. Effects of cattle, sheep and deer grazing on soil physical quality and losses of phosphorus and suspended sediment losses in surface runoff. Agric. Ecosyst. Environ. 2010, 140, 264–272. [Google Scholar]
  27. Nash, D.; Hannah, M.; Halliwell, D.; Murdoch, C. Phosphorus in runoff from a pasture based grazing system. In Proceedings of the National Soils Conference, Environmental Benefits of Soil Management, Brisbane, Australia, 27–29 April 1998; pp. 447–449. [Google Scholar]
  28. McDowell, R.W.; Littlejohn, R.P.; Blennerhassett, J.D. Phosphorus fertilizer form affects phosphorus loss to waterways: A paired catchment study. Soil Use Manag. 2010, 26, 365–373. [Google Scholar]
  29. Mc Dowell, R.W.; Catto, W. Alternative fertilisers and management to decrease incidental phosphorus loss. Environ. Chem. Lett. 2005, 2, 169–174. [Google Scholar]
  30. Van Vuuren, D.P.; Bouwman, A.F.; Beusen, A.H.W. Phosphorus demand for the 1970–2100 period: A scenario analysis of resource depletion. Glob. Environ. Chang. 2010, 20, 428–439. [Google Scholar]
  31. Brar, M.S.; Bijay-Singh; Bansal, S.K.; Srinivasarao, C. Role of Potassium Nutrition in Nitrogen Use Efficiency in Cereals. e-ifc 2011, 29, 20–27. [Google Scholar]
  32. Rossini, F.; Provenzano, M.E.; Sestili, F.; Ruggeri, R. Synergistic effect of sulphur and nitrogen in the organic and mineral fertilization of durum wheat: Grain yield and quality traits in the Mediterranean environment. Agronomy 2018, 8, 189–205. [Google Scholar]
  33. Rashid, A.; Awan, Z.I.; Ryan, J. Diagnosing phosphorus deficiency in spring wheat by plant analysis: Proposed critical concentration ranges. Commun. Soil Sci. Plant Anal. 2005, 36, 609–622. [Google Scholar] [CrossRef]
  34. Rutkowska, A. Sensitivity of plant and soil indices in evaluating the long-term consequences of soil mining from reserves of phosphorus, potassium, and magnesium. Commun. Soil Sci. Plant Anal. 2013, 44, 377–389. [Google Scholar] [CrossRef]
  35. Usherwood, N.R.; Segars, W.I. Nitrogen interactions with phosphorus and potassium for optimum crop yield, nitrogen use effectiveness and environmental stewardship. Sci. World 2001, 1, 57–60. [Google Scholar] [CrossRef]
  36. Reich, P.B.; Oleksyn, J.; Wright, I.J. Leaf phosphorus influences the photosynthesis-nitrogen relation: A cross-biome analysis of 314 species. Oecologia 2009, 160, 207–212. [Google Scholar] [CrossRef]
  37. IUSS Working Group WRB. World Reference Base for Soil Resources 2014, Update 2015 International Soil Classification System for Naming Soils and Creating Legends for Soil Maps; World Soil Resources Reports; FAO: Rome, Italy, 2015; p. 106. [Google Scholar]
  38. Egnér, H.; Riehm, H.; Domingo, W.R. Untersuchungen über die chemische Bodenanalyse als Grundlage für die Beurteilung des Nährstoffzustandes der Böden. II. Chemische Extraktionsmethoden zur Phosphor-und Kaliumbestimmung. K. Lantbr. Ann. 1960, 26, 199–215. [Google Scholar]
  39. Schachtschabel, P. Das pflanzenverfügbare Magnesiumdes Bodens und seine Bestimmung. J. Plant Nutr. Soil Sci. 1954, 67, 9–23. [Google Scholar]
  40. EU Nitrogen Expert Panel. Nitrogen Use Efficiency (NUE) an Indicator for the Utilization of Nitrogen in Food Systems; Wageningen University: Wageningen, The Netherlands, 2015. [Google Scholar]
  41. Fageria, N.K.; Baligar, V.C.; Moreira, A.; Portez, T. Dry bean genotypes evaluation for growth, yield components and phosphorus use efficiency. J. Plant Nutr. 2010, 33, 2167–2181. [Google Scholar]
  42. Hansen, J.C.; Cade-Menun, B.J.; Strawn, D.G. Phosphorus Speciation in Manure-Amended Alkaline Soils. J. Environ. Qual. 2004, 33, 1521–1527. [Google Scholar] [CrossRef] [PubMed]
  43. Shen, J.; Yuan, L.; Zhang, J.; Li, H.; Bai, Z.; Chen, X.; Zhang, W.; Zhang, F. Phosphorus Dynamics: From Soil to Plant. Plant Physiol. 2011, 156, 997–1005. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Harrison, A.F. Soil Organic Phosphorus—A Review of World Literature; CAB International: Oxon, UK, 1987; pp. 1–257. [Google Scholar]
  45. Oelkers, E.H.; Valsami-Jones, E. Phosphate mineral reactivity and global sustainability. Elements 2008, 4, 83–87. [Google Scholar] [CrossRef]
  46. Pierzynski, G.M.; McDowell, R.W.; Sims, J.T. Chemistry, Cycling, and Potential Movement of Inorganic Phosphorus in Soils. In Phosphorus, Agriculture and the Environment; American Society of Agronomy: Madison, WI, USA, 2005; pp. 53–86. [Google Scholar]
  47. Kruse, J.; Abraham, M.; Amelung, W.; Baum, C.; Bol, R.; Kühn, O.; Lewandowski, H.; Niederberger, J.; Oelmann, Y.; Rüger, C.; et al. Innovative methods in soil phosphorus research: A review. J. Plant Nutr. Soil Sci. 2015, 178, 43–88. [Google Scholar] [CrossRef] [Green Version]
  48. Turner, B.L.; Papházy, M.J.; Haygarth, P.M.; McKelvie, I.D. Inositol phosphates in the environment. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 2002, 357, 449–469. [Google Scholar] [CrossRef] [Green Version]
  49. Condron, L.M.; Turner, B.L.; Cade-Menun, B.J. Chemistry and dynamics of soil organic phosphorus. In Phosphorus: Agriculture and the Environment; Sims, J.T., Sharpley, A.N., Eds.; American Society of Agronomy, Crop Science Society of America, Soil Science Society of America, Inc.: Madison, WI, USA, 2005; pp. 87–121. [Google Scholar]
  50. Jones, D.A.; Smith, B.F.L.; Wilson, M.J.; Goodman, B.A. Solubilisator fungi of phosphate in rise soil. Mycol. Res. 2003, 95, 1009–1093. [Google Scholar]
  51. Lambers, H.; Shane, M.W.; Cramer, M.D.; Pearse, S.J.; Veneklaas, E.J. Root Structure and Functioning for Efficient Acquisition of Phosphorus: Matching Morphological and Physiological Traits. Ann. Bot. 2006, 98, 693–713. [Google Scholar] [CrossRef] [Green Version]
  52. Simpson, R.J.; Oberson, A.; Culvenor, R.A.; Ryan, M.A.; Veneklaas, E.J.; Lambers, H.; Lynch, J.P.; Ryan, P.R.; Delhaize, E.; Smith, F.A.; et al. Plant and microbial strategies to improve phosphorus efficiency of agriculture. Plant Soil 2011, 349, 89–120. [Google Scholar]
  53. Heffer, P.; Prud’homme, M.P.R.; Muirheid, B.; Isherwood, K.F. Phosphorus fertilisation: Issues and outlook. In Proceedings of International Fertilisers Society; IFS: York, UK, 2006; Volume 586, pp. 1–23. [Google Scholar]
  54. Johnson, A.E.; Syers, J.K. Changes in understanding the behaviour of soil and fertilizer phosphorus: Implications for their efficient use in agriculture. In Proceedings of International Fertilisers Society; IFS: York, UK, 2006; Volume 589, pp. 1–31. [Google Scholar]
  55. Benbi, D.K.; Biswas, C.R. Nutrient budgeting for phosphorus and potassium in a long-term fertilizer trial. Nutr. Cycl. Agroecosyst. 1999, 54, 125–132. [Google Scholar]
  56. Sun, B.H.; Cui, Q.H.; Guo, Y.; Yang, X.Y.; Zhang, S.L.; Gao, M.X.; Hopkins, D.W. Soil phosphorus and relationship to phosphorus balance under long-term fertilization. Plant Soil Environ. 2018, 64, 214–220. [Google Scholar]
  57. Lasaletta, L.; Billen, G.; Grizetti, B.; Garnier, J.; Leach, A.M.; Galloway, J.N. 50 years trends in nitrogen use efficiency of world cropping systems: The relationships between yield and nitrogen input to cropland. Environ. Res. Lett. 2014, 9, 10. [Google Scholar] [CrossRef]
  58. Scharf, P.C.; Lory, J.A. Calibration corn color from aerial photographs to predict side dress nitrogen nee. Agron. J. 2002, 94, 397–404. [Google Scholar] [CrossRef]
  59. Arregui, L.M.; Quemada, M. Strategies to improve nitrogen use efficiency in winter cereal production under rainfed conditions. Agron. J. 2008, 100, 277–284. [Google Scholar] [CrossRef]
  60. Sharma, L.K.; Bali, S.K. A review of methods to improve nitrogen use efficiency in agriculture. Sustainability 2017, 10, 51. [Google Scholar] [CrossRef] [Green Version]
  61. Bouchet, A.S.; Laperche, A.; Bissuel-Belaygue, C.; Snowdon, R.J.; Nesi, N.; Stahl, A. Nitrogen use efficiency in rapeseed. A review. Agron. Sustain. Dev. 2016, 36, 1–20. [Google Scholar] [CrossRef]
  62. Weiser, C.; Fuß, R.; Kage, H.; Flessa, H. Do farmers in Germany exploit the potential yield and nitrogen benefits from preceding oilseed rape in winter wheat cultivation? Arch. Agron. Soil Sci. 2017, 64, 1–13. [Google Scholar] [CrossRef]
  63. Bastani, S.; Hajiboland, R. Uptake and utilization of applied phosphorus in oilseed rape (Brassica napus L. cv. Hayola) plants at vegetative and reproductive stages: Comparison of root with foliar phosphorus application. Soil Sci. Plant Nutr. 2017, 63, 254–263. [Google Scholar] [CrossRef]
  64. Grzebisz, W.; Szczepaniak, W.; Barłóg, P.; Przygocka-Cyna, K.; Potarzycki, J. Phosphorus sources for winter oilseed rape (Brassica napus L.) during reproductive growth—Magnesium sulfate management impact on P use efficiency. Arch. Agron. Soil Sci. 2018, 64, 1646–1662. [Google Scholar] [CrossRef]
  65. Föhse, D.; Claassen, N.; Jungk, A. Phosphorus efficiency of plants: II. Significance of root radius, root hairs, and cation-anion balance for phosphorus influx in seven plant species. Plant Soil 1991, 132, 261–272. [Google Scholar] [CrossRef]
  66. Lickfett, T.; Mätthaus, B.; Velasco, L.; Möllers, C. Seed yield, oil and phytate concentration in the seeds of two oilseed rape cultivars as affected by different phosphorus supply. Eur. J. Agron. 1999, 11, 293–299. [Google Scholar] [CrossRef]
  67. Quemada, M.; Lassaletta, L.; Jensen, L.S.; Godinot, O.; Brentrup, F.; Buckley, C.; Foray, S.; Hvid, S.K.; Oenema, J.; Richards, K.G.; et al. Exploring nitrogen indicators of farm performance among farm types across several European case studies. Agric. Syst. 2020, 77, 1–14. [Google Scholar] [CrossRef]
  68. Fertilizers Europe. Forecast of Food Farming and Fertilizer Use in the European Union 2019–2029. Sustainable agriculture in Europe. 2019. Available online: https://www.fertilizerseurope.com/wp-content/uploads/2019/12/Forecast-of-food-farming-and-fertilizer-use-in-the-European-Union.pdf (accessed on 14 October 2020).
  69. GUS (Statistics Poland). Means of Production in Agriculture in the 2018/19 Farming Year. Warsaw, Poland; 2020. Available online: https://stat.gov.pl/download/gfx/portalinformacyjny/pl/defaultaktualnosci/5507/6/16/1/srodki_produkcji_w_rolnictwie_w_roku_gospodarczym_2018-2019.pdf (accessed on 1 November 2020).
  70. McLaughlin, M.J.; McBeath, T.M.; Smernik, R.; Stacey, S.P.; Ajiboye, B.; Guppy, C. The chemical nature of P accumulation in agricultural soils—Implications for fertiliser management and design: An Australian perspective. Plant Soil 2011, 349, 69–87. [Google Scholar] [CrossRef]
  71. Mühlbachová, G.; Čermák, P.; Vavera, R.; Káš, M.; Pechová, M.; Marková, K.; Kusá, H.; Růžek, P.; Hlušek, J.; Lošák, T. Boron availability and uptake under increasing phosphorus rates in a pot experiment. Plant Soil Environ. 2017, 63, 483–490. [Google Scholar]
Figure 1. Grain yields of winter oilseed rape in (a) Grabów and (b) Baborówko; 95% Tukey Honest Significant difference test (HSD) intervals; n = 336 for Grabów and n = 288 for Baborówko.
Figure 1. Grain yields of winter oilseed rape in (a) Grabów and (b) Baborówko; 95% Tukey Honest Significant difference test (HSD) intervals; n = 336 for Grabów and n = 288 for Baborówko.
Agronomy 10 01747 g001
Figure 2. Grain yields of winter wheat in (a) Grabów and (b) Baborówko; 95% Tukey HSD intervals; n = 384 for each location.
Figure 2. Grain yields of winter wheat in (a) Grabów and (b) Baborówko; 95% Tukey HSD intervals; n = 384 for each location.
Agronomy 10 01747 g002
Figure 3. Grain yields of maize in (a) Grabów and (b) Baborówko; 95% Tukey HSD intervals; n = 384 for each location.
Figure 3. Grain yields of maize in (a) Grabów and (b) Baborówko; 95% Tukey HSD intervals; n = 384 for each location.
Agronomy 10 01747 g003
Figure 4. Grain yields of spring barley in (a) Grabów and (b) Baborówko; 95% Tukey HSD intervals; n = 384 for each location.
Figure 4. Grain yields of spring barley in (a) Grabów and (b) Baborówko; 95% Tukey HSD intervals; n = 384 for each location.
Agronomy 10 01747 g004
Figure 5. Grain yields of winter oilseed rape at increasing additions of Nitrogen (N) fertilizer at the P-plus and control levels of phosphorus (P); 95% Tukey HSD intervals; n = 336 for (a) Grabów (P fertilization HSD 0.405) and n = 288 for (b) Baborówko (P fertilization HSD 0.321).
Figure 5. Grain yields of winter oilseed rape at increasing additions of Nitrogen (N) fertilizer at the P-plus and control levels of phosphorus (P); 95% Tukey HSD intervals; n = 336 for (a) Grabów (P fertilization HSD 0.405) and n = 288 for (b) Baborówko (P fertilization HSD 0.321).
Agronomy 10 01747 g005
Figure 6. Grain yields of winter wheat at increasing additions of N fertilizer at the P-plus and control levels of P; 95% Tukey HSD intervals; n = 384 for each location; (P fertilization (a) Grabów HSD 0.550, (b) Baborówko HSD 0.568).
Figure 6. Grain yields of winter wheat at increasing additions of N fertilizer at the P-plus and control levels of P; 95% Tukey HSD intervals; n = 384 for each location; (P fertilization (a) Grabów HSD 0.550, (b) Baborówko HSD 0.568).
Agronomy 10 01747 g006
Figure 7. Grain yields of maize at increasing additions of N fertilizer at the P-plus and control levels of P; 95% Tukey HSD intervals; n = 384 for each location (P fertilization (a) Grabów HSD 1.183, (b) Baborówko HSD 0.894).
Figure 7. Grain yields of maize at increasing additions of N fertilizer at the P-plus and control levels of P; 95% Tukey HSD intervals; n = 384 for each location (P fertilization (a) Grabów HSD 1.183, (b) Baborówko HSD 0.894).
Agronomy 10 01747 g007
Figure 8. Grain yields of spring barley at increasing additions of N fertilizer at the P-plus and control levels of P, 95% Tukey HSD intervals; n = 384 for each location (P fertilization (a) Grabów HSD 0.349, (b) Baborówko HSD 0.434).
Figure 8. Grain yields of spring barley at increasing additions of N fertilizer at the P-plus and control levels of P, 95% Tukey HSD intervals; n = 384 for each location (P fertilization (a) Grabów HSD 0.349, (b) Baborówko HSD 0.434).
Agronomy 10 01747 g008
Figure 9. Running averages of P content in the soil in Grabów (a) and Baborówko (b) (2007–2018); 95% Tukey HSD intervals; n = 288 for each location (P fertilizer Grabów HSD 4.026, Baborówko HSD 13.32).
Figure 9. Running averages of P content in the soil in Grabów (a) and Baborówko (b) (2007–2018); 95% Tukey HSD intervals; n = 288 for each location (P fertilizer Grabów HSD 4.026, Baborówko HSD 13.32).
Agronomy 10 01747 g009
Table 1. P use indices for winter oilseed rape; N fertilizer rate; P fertilizer rate; Pn (uptake) and Ps (surplus) (kg·ha−1); P use efficiency (PUE) (%).
Table 1. P use indices for winter oilseed rape; N fertilizer rate; P fertilizer rate; Pn (uptake) and Ps (surplus) (kg·ha−1); P use efficiency (PUE) (%).
N
Rate
GrabówBaborówko
P-PlusControlP-PlusControl
39 kg P·ha−10 kg P·ha−139 kg P·ha−10 kg P·ha−1
PUEPnPsPnPsPUEPnPsPnPs
03 a13.2 aA26 aA12.0 aA−12 a06.91 aA32 a6.87 aA−7 aB
503 a17.2 bA22 bA16.1 bA−16 b09.99 bA29 b10.6 bA−11 bB
1004 ab21.0 cA18 cA19.4 cA−19 c4 ab14.9 cA24 c13.2 cA−13 cB
1504 ab23.0 dA16 dA21.3 cdA−21 cd3 a17.8 dA21 d16.6 dA−16 dB
2007 b24.5 eA15 deA21.6 cdB−22 cd5 bc19.1 eA20 e17.2 eA−17d eB
25012 c26.5 fA13 fA22.5 dB−23 d5 ab21.2 fA18 f19.3 fA−19 fB
Mean6211819−186152414−14
Treatments with the same letter are not significantly different (p ≤ 0.05). The capital letters denote differences between the mean of Pn, Ps in P-plus and control treatments and the small letters denote difference between N rates; n = 280 for Grabów; n = 240 for Baborówko.
Table 2. P use indices for winter wheat; N fertilizer rate; P fertilizer rate; Pn (uptake) and Ps (surplus) (kg·ha−1); PUE (%).
Table 2. P use indices for winter wheat; N fertilizer rate; P fertilizer rate; Pn (uptake) and Ps (surplus) (kg·ha−1); PUE (%).
N
Rate
GrabówBaborówko
P-PlusControlP-PlusControl
31 kgP·ha−10 kgP·ha−131 kgP·ha−10 kgP·ha−1
PUEPnPsPnPsPUEPnPsPnPs
05 ab15.0 aA16 aA13.5 a−14 aB-8.85 aA22 aA9.00 aA−9 aB
406 bc17.3 bA14 bA15.5 bB−16 bB-13.0 bA18 bA13.9 bA−14 bB
804 a19.5 cA12 cA18.3 cA−18 cB4 a17.5 cA14 cA16.4 cA−16 cB
1204 a21.8 dA9 dA20.5 dA−2 dB-19.9 dA11 dA19.8 dA−20 dB
1606 bc22.6 deA8 deA20.6 dB−2 dB1 b21.4 eA10 eA21.2 eA−21 eB
2008 d24.1 eA7 eA21.5 dB−2 dB-21.5 eA10 eA21.5 eA−22 eB
Mean6201118−193171417−17
Treatments with the same letter are not significantly different (p ≤ 0.05). The capital letters denote differences between Pn, Ps in P-plus and control treatments and the small letters denote difference between N rates; n = 320 for each location.
Table 3. P use indices for maize; N fertilizer rate; P fertilizer rate; Pn (uptake) and Ps (surplus) (kg·ha−1); PUE (%).
Table 3. P use indices for maize; N fertilizer rate; P fertilizer rate; Pn (uptake) and Ps (surplus) (kg·ha−1); PUE (%).
N
Rate
GrabówBaborówko
P-PlusControlP-PlusControl
35 kgP·ha−10 kgP·ha−135 kgP·ha−10 kgP·ha−1
PUEPnPsPnPsPUEPnPsPnPs
011 e27.7 aA7 aA23.9 aB−24 aB5 c17.6 aA17 aA16.0 aA−16 aB
5012 e29.1 abA6 abA25.0 abB−25 abB2 a21.8 bA13 bA21.0 bA−22 bB
1004 b28.8 abcA6 abA27.3 bA−27 bB3 ab24.8 cA10 cA23.6 cA−23 cB
1501 a29.4 abcA6 abA28.9 bcA−29 bcB3 ab26.6 dA8 dA25.5 dA−26 dB
2005 bc30.2 bcA9 bA28.3 bcA−28 bcB5 c27.2 deA8 dA26.8 eA−27 eB
2505 bc29.9 bA6 abA28.2 bcA−28 bcB3 ab27.7 deA7 eA27.6 eA−28 fB
Mean729727−274241023−23
Treatments with the same letter are not significantly different (p ≤ 0.05). The capital letters denote differences between Pn, Ps in P-plus and control treatments and the small letters denote difference between N rates; n = 320 for each location.
Table 4. P use indices for spring barley; N fertilizer rate; P fertilizer rate; Pn (uptake) and Ps (surplus) (kg·ha−1); PUE (%).
Table 4. P use indices for spring barley; N fertilizer rate; P fertilizer rate; Pn (uptake) and Ps (surplus) (kg·ha−1); PUE (%).
N
Rate
GrabówBaborówko
P-PlusControlP-PlusControl
31 kgP·ha−10 kgP·ha−131 kgP·ha−10 kgP·ha−1
PUEPnPsPnPsPUEPnPsPnPs
03 a13.7 aA17 aA12.7 aA−13 aB3 b6.9 aA24 aA6.0 aA−6 aB
306 b17.6 bA14 bA15.8 bB−16 bB1 a10.5 bA21 bA10.2 bA−10 bB
604 ab20.4 cA11 cA19.2 cA−19 cB2 ab14.4 cA17 cA13.8 cA−14 cB
907 bc22.0 dA9 dA19.9 cdA−20 cdB3 b17.2 dA14 dA16.4 dA−16 dB
1207 bc22.9 dA8 dA20.6 deB−21 deB4 bc19.5 eA12 eA18.3 eA−18 eB
1509 c24.0 eA7 eA21.3 eB−2 eB3 b21.2 fA10 fA20.3 fA−20 fB
Mean6201118−183151614−14
Treatments with the same letter are not significantly different (p ≤ 0.05). The capital letters denote differences between Pn, Ps in P-plus and control treatments and the small letters denote difference between N rates; n = 320 for each location.
Table 5. N use indices for winter oilseed rape; N fertilizer rate, Yn (uptake), Ns (surplus) (kg·ha−1); N use efficiency (NUE) (%) for the P fertilizer addition (P-plus) and control (no P addition).
Table 5. N use indices for winter oilseed rape; N fertilizer rate, Yn (uptake), Ns (surplus) (kg·ha−1); N use efficiency (NUE) (%) for the P fertilizer addition (P-plus) and control (no P addition).
N RateGrabówBaborówko
P-PlusControlP-PlusControl
NUEYnNsNUEYnNsNUEYnNsNUEYnNs
50145 a72 a−22 a140 a70 a−20 a91 a46 a4 a100 a50 a0 a
100100 b102 b−2 a100 b101 b−1 a73 b74 b26 b67 b70 ab30 b
15080 ab120 bc30 b78 c117 bc33 b65 c98 c52 c59 c89 c61 c
20070 b140 cd60 c67 cd134 c66 c54 d107 c93 d52 cd104 d96 d
25062 b155 d95 d57 d141 c109 d47 d118 d132 e46 d114 e136 e
Mean91118328811237668961658681
The small letters denote difference between NUE, Yn, Ns in P-plus and control treatments in the same location; Treatments with the same letter are not significantly different (p ≤ 0.05); n = 280 for Grabów; n = 240 for Baborówko.
Table 6. N use indices for winter wheat; N fertilizer rate, Yn (uptake), Ns (surplus) (kg·ha−1); NUE (%) for the P fertilizer addition (P-plus) and control (no P addition).
Table 6. N use indices for winter wheat; N fertilizer rate, Yn (uptake), Ns (surplus) (kg·ha−1); NUE (%) for the P fertilizer addition (P-plus) and control (no P addition).
N RateGrabówBaborówko
P-PlusControlP-PlusControl
NUEYnNsNUEYnNsNUEYnNsNUEYnNs
40196 a78 a−38 a193 a77 a−37 a155 a62 a−22 a153 a61 a−21 a
80134 b107 b−27 b131 b105 b−25 b111 b89 b−9 b106 b85 b−5 b
120106 c128 c−8 c103 c123 c−3b c89 c107 c13 c89 c107 c13 c
16092 d147 d13 d91 cd146 d14 d78 c124 d36 d76 cd122 d38 d
20083 d166 e34 e81 d162 e38 e65 d130 d69 e64 d128 d72 e
Mean122125−5120123−3100102179810119
The small letters denote difference between NUE, Yn, Ns in P-plus and control treatments in the same location; Treatments with the same letter are not significantly different (p ≤ 0.05); n = 320 for each location.
Table 7. N use indices for maize; N fertilizer rate, Yn (uptake), Ns (surplus) (kg·ha−1); NUE (%) for the P fertilizer addition (P-plus) and control (no P addition).
Table 7. N use indices for maize; N fertilizer rate, Yn (uptake), Ns (surplus) (kg·ha−1); NUE (%) for the P fertilizer addition (P-plus) and control (no P addition).
N RateGrabówBaborówko
P-PlusControlP-PlusControl
NUEYnNsNUEYnNsNUEYnNsNUEYnNs
50231 a115 a−65 a227 a114−64 a195 a97 a−47 a199 a99 a−49 a
100154 b154 b−54 a145 b145−45 b145 b145 b−45 a138 b138 b−38 a
150107 c161 b−11 b99 c1482 c112 c168 c−18 b104 c156 c−6 b
20089 d179 c21 c85 cd17030 d87 d175 cd25 c87 d175 d25 c
25071 e175 c75 d69 d17377 e74 e185 d65 d76 d191 e59 d
Mean130156−71251500123154−4121151−2
The small letters denote difference between NUE, Yn, Ns in P-plus and control treatments in the same location; Treatments with the same letter are not significantly different (p ≤ 0.05); n = 320 for each location.
Table 8. N use indices for spring barley; N fertilizer rate, Yn (uptake), Ns (surplus) (kg·ha−1); NUE (%) for the P fertilizer addition (P-plus) and control (no P addition).
Table 8. N use indices for spring barley; N fertilizer rate, Yn (uptake), Ns (surplus) (kg·ha−1); NUE (%) for the P fertilizer addition (P-plus) and control (no P addition).
N RateGrabówBaborówko
P-PlusControlP-PlusControl
NUEYnNsNUEYnNsNUEYnNsNUEYnNs
30211 a63 a−33 a202 a61 a−31 a132 a39 a−9 a135 a41 a−9 a
60136 b81 b−21 b132 b78 b−19 b101 b61 b−1 b103 b62 b−2 b
90114 c102 c−12 c106 c95 c−5 c89 ab81 c9 c82 c74 c16 c
12099 d119 d1 d96 c114 d6 d77 bc92 d28 d76 cd91 d29 d
15088 d132 e18 e81 d121 d29 e70 cd105 e45 e66 d99 e51 e
Mean13099−1012394−4947614927317
The small letters denote difference between NUE, Yn, Ns in P-plus and control treatments in the same location; Treatments with the same letter are not significantly different (p ≤ 0.05); n = 320 for each location.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rutkowska, A.; Skowron, P. Productive and Environmental Consequences of Sixteen Years of Unbalanced Fertilization with Nitrogen and Phosphorus—Trials in Poland with Oilseed Rape, Wheat, Maize and Barley. Agronomy 2020, 10, 1747. https://doi.org/10.3390/agronomy10111747

AMA Style

Rutkowska A, Skowron P. Productive and Environmental Consequences of Sixteen Years of Unbalanced Fertilization with Nitrogen and Phosphorus—Trials in Poland with Oilseed Rape, Wheat, Maize and Barley. Agronomy. 2020; 10(11):1747. https://doi.org/10.3390/agronomy10111747

Chicago/Turabian Style

Rutkowska, Agnieszka, and Piotr Skowron. 2020. "Productive and Environmental Consequences of Sixteen Years of Unbalanced Fertilization with Nitrogen and Phosphorus—Trials in Poland with Oilseed Rape, Wheat, Maize and Barley" Agronomy 10, no. 11: 1747. https://doi.org/10.3390/agronomy10111747

APA Style

Rutkowska, A., & Skowron, P. (2020). Productive and Environmental Consequences of Sixteen Years of Unbalanced Fertilization with Nitrogen and Phosphorus—Trials in Poland with Oilseed Rape, Wheat, Maize and Barley. Agronomy, 10(11), 1747. https://doi.org/10.3390/agronomy10111747

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop