Research Article

 

Nitrogen migration in crop rotations differing in fertilisation

 

Saulius Guzys

Water Resources Engineering Institute, Aleksandras Stulginskis University, Universiteto 10, LT-53361 Akademija, Kaunas Distr., Lithuania

Stefanija Miseviciene

Water Resources Engineering Institute, Aleksandras Stulginskis University, Universiteto 10, LT-53361 Akademija, Kaunas Distr., Lithuania

 

Abstract

Inappropriate use of nitrogen fertilisers is becoming a global problem; however, continuous fertilisation with N fertiliser ensures large and constant harvests. To evaluate the relationships of differently fertilised cultivated plant rotation with N metabolism in the agroecosystem the research was conducted between 2006 and 2013 at Lipliūnai, Lithuania, in fields with calcareous gley brown soil, i.e. Endocalcari Endohypogleyic Cambisol (CMg-n-w-can). The research area covered three drained plots where crop rotation of differently fertilised cereals and perennial grasses were applied. The greatest productivity was found in a higher fertilisation (TII, 843 kg N/ha) cereals crop rotation. With less fertilisation (TI, 540 kg N/ha) crop rotation productivity of cereals and perennial grasses (TIII, 218 kg N/ha) was 11-35% lower. The highest amount of mineral soil N (average 76 kg/ha) was found in TI. It was influenced by fertilisation (r=0.71) and crop productivity (r=0.39). TIII tended to reduce Nmin (12.1 mg/L) and Ntotal (12.8 mg/L) concentrations in drainage water and leaching of these elements (7 and 8 kg/ha). Nmin and Ntotal concentrations in the water depended on crop productivity respectively (r=0.48; r=0.36), quantity of mineral soil N (r=0.65; r=0.59), fertilisation (r=0.59; r=0.52), and N balance (r=0.26; r=0.35). Cereal crop rotation increased N leaching by 12-42%. The use of all crop rotations resulted in a negative N balance. Nitrogen balance depended on fertilisation with N fertiliser (r=0.55). The application of perennial grasses crop rotation in agricultural fields was the best environmental tool, reducing N migration to drainage.

Additional key words: drainage; leaching; nitrogen; balance; yield.

Abbreviations used: a.m. (active matter); r (pair correlation); LSD05 (limit of reliable (95%) difference); tfact. (estimated Student’s test statistic); ttheor.95% (Student’s test value at 95% confidence level); Nmin (mineral nitrogen); Ntotal (total nitrogen).

Citation: Guzys, S.; Miseviciene, S. (2015). Nitrogen migration in crop rotations differing in fertilisation. Spanish Journal of Agricultural Research, Volume 13, Issue 2, e0303, 13 pages. .http://dx.doi.org/10.5424/sjar/2015132-6672.

Received: 12 Aug 2014. Accepted: 18 May 2015

Copyright © 2015 INIA. This is an open access article distributed under the Creative Commons Attribution License (CC by 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Funding: The authors received no specific funding for this work.

Competing interests: The authors have declared that no competing interests exist.

Correspondence Correspondence should be addressed to Stefanija Miseviciene: kestmis@rygveda.lt


 

CONTENTS

Abstract

Introduction

Material and methods

Results

Discussion

References

IntroductionTop

Nitrogen is the main crop nutrition element that has some impact on its productivity (Oenema et al., 2009). Continuous fertilisation with N fertiliser ensures large and constant harvests. Nitrogen excess in agricultural production has a negative impact on the environment (leached nitrates, evaporating ammonia and N oxide) (Ahlgren et al., 2008). The abuse of N fertiliser is becoming a global problem and affects the environment in a number of ways: (1) human health problems caused by nitrates (Powlson et al., 2008); (2) livestock health problems (Lundberg et al., 2008); (3) surface water eutrophication (Smith & Schindler, 2009); (4) formation of nitrous acid, which is linked to acid rain (Menz & Seip, 2004); (5) depletion of the ozone layer in the atmosphere because of atmospheric nitrate oxides, which strengthen ultraviolet radiation (Savci, 2012; Rosenstock et al., 2013); and (6) causes the global warming effect (Martinez-Blanco et al., 2014; Payen et al., 2015).

In terms of non-point source pollution, N constitutes 52 to 61% (Klaushal et al., 2011). Agricultural pollution constitutes on average 30 to 35% of N leaching (Klaushal et al., 2011). Due to particular agricultural operations in Lithuania, 36,000 tonnes of N per year reach the Baltic Sea waterborne load (Larsson & Granstedt, 2010). Previous research has shown that leaching of elements and compounds from the soil is essentially determined by the hydrothermal regime, the amount of these element in the soil and their fertilisation (Beeston et al., 2010; Komonweeraket et al., 2011). A higher rate of Nmin leaching occurs in humid years (Sogbedji et al., 2001). Views on fertilisation systems constituting human influence on water quality are rather mixed in the scientific literature. While some authors claim that mineral fertilisation systems are more environmentally dangerous (Kokkora & Hann, 2007; Ulen et al., 2012), others claim the opposite, that organic systems are more dangerous (Hossain et al., 2007). Nonetheless, it is known that unproductive fertiliser losses and water pollution increase considerably through the irrational and unbalanced use of both mineral and organic fertiliser (Hongxu et al., 2011; Hyytiäinen et al., 2011). It was determined in Finland that the yield of crops increases substantially, while N leaching falls with balanced fertilisation with N carried out several times during the vegetation of crops (Salvioni et al., 2011). As tests have shown, NO3- N leaching from the soil can be reduced significantly by changing scattered fertilisation into local (Nakamura et al., 2004).

Current knowledge suggests that the scale of leaching of chemical elements and compounds should be inversely proportional to the agroecosystem’s total biomass synthesis. Cabrera et al. (2007) determined that a bigger biomass synthesis and dryer year limit N leaching by drainage and ground water. Securing a rich harvest (crop rotation with cover crops and continuously occupied soil, plants of high biological potential, balance and rational fertilisation and application of organic production) reduces the potential water pollution risk (Thorup-Kristensen et al., 2012; Buckley & Carney, 2013).

The quality of water which filters though the soil and in the surface runoff depends on the changes that occur in the soil when a variety of crops is grown. The fact that around 67% of the leached N is from the soil has been known for a very long time (Ulen & Johansson, 2009). Research data on N migration and its compounds’ dynamics in the soil in the agroecosystems affected by human activity is not substantial.

Based on statistical data, cereals crop areas have increased since 1991 in Lithuania and constitute 34.3 to 41.8% of the total crop area. Perennial and annual grass areas reduced from 46 to 41%. For many years, large areas of perennial grasses, which are used for a longer time without inclusion in crop rotations, are used to grow cereals successively. The farmer’s decision to grow more of one or another kind of crop, as well as the general farming culture, is affected by economic interest. Such a one-sided approach can worsen the indicators of the agro landscape and become a potential non-point pollution source, even when other agrarian measures are used correctly and rationally.

On the basis of the data obtained, we have proposed the hypothesis that non-point source pollution of the agro-landscape is inversely proportional to the synthesis of total biomass of the ecosystem, and that all of the available means ensuring high crop productivity actually reduce the potential risk of water pollution. The novelty of this paper is the evaluation of the effect of crop rotations under different fertilisation levels on the environment. The aim of the research was to evaluate the relationship between crop rotations of different intensity and biological potential with the migration of N.

Material and methodsTop

Site description

The experiment was performed in the Middle Lithuania lowland in the basin of Graisupis (area 16.6 km2) and in the village of Lipliūnai (N 55°18′, E 23°51′). A subsurface drainage site was used for the investigations. Investigations were carried out from 2006 to 2013. The experimental site included three plots (Fig. 1).

Figure 1. Geographical location and experimental layout of the Lipliūnai site.

The soil in the site is non-acid Endocalcari - Endohypogleyic Cambisol (CMg-n-w-can) sandy light loam and sandy loam on sandy loam and sandy light loam (Liekis, 2001). According to the size distribution of soil particles, the research soil is rather homogenous. The arable layer of TI and TII is sandy light loam. Its composition in deeper layers lightens until sandy loam. The soil in TIII is slightly lighter (Table 1). In the upper (0-20 cm) soil layer, the soil was average and high in humus (1.79-4.32 %). Agrochemical data analysis showed that differences in chemical elements of the soil among the treatments at the beginning of the research were not significant, except TIII, where we determined that N (tfact.=2.82>ttheor.95%=2.3) and humus (tfact.=2.47>ttheor.95%=2.3) were significantly lower than in TI. However, as we projected to use the perennial grass crop rotation in TIII, where perennial grasses fix N uses legume bacteria, we decided to start the research.


Table 1. Soil properties in different layers before the start of the research


Meteorological conditions are described using the data of the Dotnuva meteorological station, which is 8 km away from the study area.

Experimental design

In each of these plots, a field crop rotation with different fertilisations using mineral fertilisers was used. This was a production experiment; therefore, in different years at the beginning of the experiment, the amount of fertiliser used depended on the financial state of the farm. The following crops were grown: sugar-beet (Beta vulgaris L. var. sacharifera) ‘Belmonte’, spring barley (Hordeum vulgare L.) ‘Ūla’, winter wheat (Triticum aestivum Host) ‘Portal’, perennial grass red clover (Trifolium pratense L.) ‘Liepsna’ + feeding timothy (Phleum pratense L.) ‘Gintaras II’, summer wheat (Triticum aestivum L.) ‘Nandu’ and winter rape (Brassica napus L.) ‘Valesca’. The long-term field experiment in Lipliūnai consisted of three crop rotations with 75, 87.5 and 50% proportions of cereals (Table 2).


Table 2. Field crop rotations (Rot) with proportions of 75, 87.5 and 50% of cereals and its N fertilisation (Fert, kg/ha a.m.) and crop yield (t/ha, dry matter)


Only mineral fertilisers were applied to crop fields. Nitrogen fertiliser (as well as phosphorus and potassium) was applied in the spring. The fertilisation schedule is presented in Table 2. The N fertiliser used was ammonium nitrate (NH4NO3). A dispersed fertilisation method was adopted, using a suspended fertiliser spreader “Amazone”. The scattered fertiliser amount in higher or lower fertilisation fields was determined by academic recommendations and the financial position of the production farm.

Soil, plant and water analyses

Before starting the research, there were two soil profiles dug out in each of the treatments, and their soil genetic horizons were determined as well as samples from each of the layer were taken to determine physical and agrochemical characteristics from 0-24, 24-55 and 55-91 cm depth. The size distribution of soil particles (or soil texture) was determined by the Fere triangle method (Liekis, 2001), and soil density by the weight method in intact structure samples (Vadiunina & Korcagina, 1986). In order to determine the agrochemical soil characteristic, samples were taken in the autumn, after the harvest. Each year, in order to determine the agrochemical properties of the soil, samples were taken in autumn after harvest. Soil samples were collected using a drill up to 60 cm depth at every 20 cm. In TI and TII there were 15 boreholes drilled, while in TIII there were 20 boreholes. In each experiment treatment there were 3 collective samples created, which were analysed in the laboratory.

The soil pH was determined potentiometrically in a soil suspension in 1 M KCl (ISO 10390:2005); the humus content was determined by the method of Tyurin (1937); the Ntotal was determined by Kjeldahl’s procedure (ISO 11261:1995); available P2O5 and K2O were determined by the AL method (soil:solution ratio – 1:20, 0.1 M ammonium lactate + 0.4 N acetic acid, pH 3.74, 2 h shaking) (Egner et al., 1960); and Nmin amount (NH4+ -N + NO3-N) contained in the soil was measured from the extraction of 0.2 M KCl using a flow injection analyser (FIA star 5012, AN 551, ASN 65-31-84, and ASN 50-01-84) (ISO/TS 14256-1:2003).

The yield was harvested mechanically and crop samples were taken from all replications. The cereals and grasses crop yield was harvested mechanically. The yield of row crops was determined manually. The stacionary area of the test fields of crops and grass was 30 m2, row crops covered an area of 45 m2. Crop harvest was determined by 6 replications in TI and TII, and by 8 in TIII. The main and secondary production of crops were weighed separately. The samples of the main and secondary production of yield were taken at the same time. Twenty different plants of sugar beet were taken in each test field (in a diagonal way). Leaves were sampled taking 2 kg of leaf mass. During the crop harvest time 3 kg of grain and 2 kg of straw were taken from each test field. Grass samples were taken during the grass haymaking time when a certain amount of grass was removed from 10 (in TI and TII) to 15 (in TIII) places in each test field. The amount of dry matter was determined in all samples used for yield production by drying at 105°C and reweighing to constant weight. Total crop yield energy was calculated using equivalents, which were multiplied with dry matter yield (Jankauskas et al., 2000). Ntotal in crop yield was determined after burning with sulphuric acid (ISO 5983-1:2005/Cor. 1:2008).

Drainage water discharge was measured by a volumetric method at every three days. Daily drainage was calculated by linear interpolation. Samples of drainage water were taken every 10 days. Monthly Nmin and Ntotal leaching was calculated having multiplied the average monthly nutrient concentration by the monthly drainage value. N concentration in water was determined colorimetrically, NH4 + by the gas diffusion method (ISO 11732:2005), Ntotal was determined after organic matter mineralisation with potassium persulphate and NO3 by the cadmium reduction method (ISO 13395:1996).

Data analysis

The data were processed by statistical methods. Disperse analysis and correlation and regression methods were applied (Shibusawa & Hache, 2009). The statistical significance of the data was computed by a multiple-way analysis of variance (ANOVA). In addition, the LSD multiple range test was used, the means of which differed significantly using a significance level of 0.01 and 0.05%. Student’s test statistic was used to determine the reliability of differences between treatments.

ResultsTop

Data show that perennial grasses and sugar beets produced the biggest yield (3.7-10.9 and 10.3-10.4 t/ha of dry material respectively). In all years, the higher fertilisation (843 N kg/ha) cereal crop rotation yield was determined to be larger than in the lower fertilisation (540 N kg/ha) cereal crop rotation (Table 2).

As different types of crops were grown in the research field, in order to determine which crop rotation field was the most productive, each of the yields of dry material amount was converted to GJ/ha and cumulative curves were drawn. After eight years of research, the biggest productivity was found in crop rotation of cereal with higher fertilisation, which accumulated 1483 GJ/ha of total energy harvest (Fig. 2).

Figure 2. The influence of field crop rotations differing in intensity on cumulative field crops yield.

Less fertilised cereal crop rotation productivity was 35% lower, and it was the lowest in the research range (966 GJ/ha). Perennial grass crop rotation productivity was 11% lower than the higher fertilised cereal crop rotation and constituted 1318 GJ/ha. With respect to humus and N, the soil of the site was not very homogeneous, as the variation coefficients of these indicators were 17 and 14, respectively. At the end of the research the highest humus reserves were found in the field, where perennial grasses crop rotation was applied (4.11%) (Fig. 3). A slightly lower soil humus is obtained by applying higher fertilisation cereal crop rotation (3.7%), while the lowest is obtained using a lower fertilisation cereal crop rotation (2.94%). The data analysis revealed that the effect of different crop rotations on humus accumulation was significant: between TI and TII tfact.=|–5.07|>ttheor.95%=2.11, between TI and TIII tfact.=|–8.02|>ttheor.95%=2.11, and between TII and TIII tfact.=|–2.28|>ttheor.95%=2.11. An increasing tendency in Ntotal accumulation in the soil was determined when using different fertilisation rates in the crop rotations (–0.25%, in TIII, –0.23% in TII, and –0.19% in TI). However, the differences between the treatments were not significant, between TI and TII tfact.=|–0.96|<ttheor.95%=2.11; between TI and TIII – tfact.=|–1.53|<ttheor.95%=2.11; and between TII and TIII - tfact.=|–0.53|<ttheor.95%=2.11.

Figure 3. Variation in humus content (a) and total nitrogen (b) in the top layer (0-20 cm). x (arithmetical mean); Rv (amplitude of variation); Sx (error of average); V% (variation’s coefficient).

The most stable Nmin concentration and accumulated amount were in TII, where variation coefficients of this element were 22% and 18%, respectively. With TI and TIII, Nmin concentration and accumulated amount in the soil varied more (from –41 to 35% and from 51 to 57%, respectively; Table 3). The average perennial data for Nmin concentration and accumulated amount in the soil are also rather uneven. The data show that in TI Nmin concentration and its accumulation in the soil was lower than in TII, however, this difference was not statistically significant (tfact.=|–0.94|<ttheor.95%=2.14 and tfact.=|–0.79|<ttheor.95%=2.20, respectively). Statistically significant differences between Nmin concentration and its accumulation in the soil were found in TI and TIII and in TII and TIII (tfact.=2.98>ttheor.95%=2.14 and tfact.=2.23>ttheor.95%=2.14; and tfact.=5.27>ttheor.95%=2.14 and tfact.=3.52>ttheor.95%=2.14, respectively).


Table 3. The impact of different crop rotations on Nmin concentration (mg/kg)and Nmin amount (kg/ha) in the soil (0-60 cm)


Nmin concentration and its accumulated amount in the soil in higher and lower fertilisation cereal crop rotations were not statistically significant, even though the amount of this element in a lower fertilisation cereal crop rotation was determined to be lower, 13 and 9% up to 6.6 mg/kg and 69 kg/ha, respectively. The lowest Nmin concentration and its accumulated amount in the soil had the most productive perennial grasses crop rotation, and these indicators reached 3.3 mg/kg and 42 kg/ha,respectively. Correlation analysis results showed that Nmin amount was closely related to fertilisation. Also, its amount in the soil had a tendency to fall when the yield was high (Fig. 4).

Figure 4. Dependence of mineral nitrogen amount in soil kg/ha (y) (0–60 cm) on crop yieldGJ/ha (x1) and nitrogen fertilisation kg/ha a.m. (x2).

A weak reverse linear was determined between Nmin amount in the soil and crop productivity. Increasing field crop productivity reduces the amount of Nmin in the soil, whereas increasing field crop fertilisation by N fertiliser increases the amount of Nmin in the soil. Field crop productivity determined the fluctuation of Nmin quantity by around 15%, and fertilisation by around 55%.

Precipitation dispersion was uneven during the research period, and it had effect on the formation of drainage (Table 4). It can be seen from the results that four of the years were dry (2006, 2008, 2011, and 2013), when precipitation was 80, 97, 97, and 91%, respectively, of the climate normals; 2007, 2009, 2010, and 2012 were damp years, because precipitation was 113, 125, 122, and 123%, respectively, of the climate normals. During almost the entire research period, the average air temperature was higher than the climate normals, except in 2010, when it corresponded to it.


Table 4. Meteorological conditions and drainage of the study periods


Drainage quantity depended on the precipitation level and changed according to the direct linear pattern, as a higher level of precipitation increased the drainage.

According to the research data, on average during 8 years in a more productive, higher fertilisation cereal (TII: 55 ± 23 mm) and grass crop rotation there was a falling tendency in drainage determined, while in the least productive and lower fertilisation cereal (TI: 71 ± 29 mm) crop rotation there was an increasing tendency observed.

Even though drainage was 29% higher in TI and 25% higher in TIII than in TII, the differences between these treatments were not significant: between TI and TII, tfact.=1.20<ttheor.95%=2.14; between TI and TIII, tfact.=0.12<ttheor.95%=2.14; and between TII and TIII, tfact.=–0.98<ttheor.95%=2.14.

Drainage essentially depended on precipitation (r=0.68; tfact.=4.31>ttheor.95%=2.15). More intensively grown, larger biological potential crops, which use greater amounts of water, reduce not only soil humidity, but also drainage (r=0.43); tfact.=2.23>ttheor.95%=2.15) (Fig. 5). In TIII, Nmin and Ntotal concentrations were found to be 12.1 and 12.8 mg/L, respectively, in drainage water (Table 5).

Figure 5. Drainage dependency on precipitation and crop yield.


Table 5. Different crop rotation impact on average yearly N concentration in drainage water and leaching (arithmetical mean ± average error)


Using a lower fertilised cereal crop rotation (TI) meant that Nmin and Ntotal concentrations had a tendency to increase by 13 and 24% to 13.7 and 15.9 mg/L, respectively. An increasing tendency was also determined for Nmin and Ntotal concentrations in drainage water in TII (15.6 and 18.5 mg/L). According to correlation and regression data analysis (Table 6), N concentration in water was determined by a number of factors.


Table 6. The dependence of Nmin (y1) and Ntotal (y2) concentrations (mg/L) in drainage water on environmental factors (x)


Thus, as the field crops‘ productivity increased, N concentration in the drainage water changes according to the parabolic pattern and it was at its lowest when the field crops‘ productivity was 150-170 GJ/ha. When a higher productivity is achieved by fertilising field crops, both Nmin and Ntotal concentrations in the drainage water increase. The relationship between soil Nmin and fertilisation by Nmin fertiliser with N concentration in drainage water can be expressed by this function: y = a – b / x (Fig. 6); as soil Nmin as well as fertilisation with N fertiliser increases, the concentration of both N forms in drainage water increases.

Figure 6. Mineral (Nmin) and total (Ntotal) leaching by drainage dependent on fertilisation with N fertiliser and the obtainable harvest.

The relationship of annual and cumulative N balance with N concentration in drainage water is of a direct linear nature and, as N balance increases, the concentration of both N forms in drainage water increases. Direct linear consistency means that Nmin concentration in drainage water increases as the soil N and humus increase. Due to the impact of field occupation coefficient N concentration in drainage water reduces.

Research data of leaching of N compounds by drainage showed (Table 5), that both Nmin and Ntotal leaching had a downward trend when using grass crop rotation (7.0 and 8.0 kg/ha respectively).

An increasing tendency was found for N leaching in lower fertilised cereal crop rotation conditions (12 and 13%) and in higher fertilised cereals crop rotation conditions, 10.0 kg Nmin/ha and 11.5 kg Ntotal/ha, respectively. Correlation and regression data analysis has shown that both Nmin and Ntotal leaching is related to the productivity of agroecosystem and factors affecting it (Fig. 6). Crop fertilisation with mineral fertiliser had a very big effect on Nmin leaching (r=0.78), while it had only average effect on Ntotal leaching (r=0.65). The yield affected Nmin leaching on average (r=0.57), while Ntotal was weakly affected (r=0.46); however, it was statistically significant.

Nitrogen leaching through the drainage relationship between the fertilisation of an agroecosystem with mineral fertiliser and its productivity is parabolic. Fertilisation with N fertiliser until 66-71 kg/ha a.m. does not increase N leaching losses. Increasing the harvest by fertilisation leads to the increase in N leaching. The lowest Nmin leaching losses were found when productivity reached 213 GJ/ha, whereas the lowest Ntotal leaching losses were when productivity reached 212 GJ/ha. The relationship between Nmin and Ntotal leaching by drainage with the balance of agroecosystem can be expressed by the formula y = a – b / x. As N balance changes from –300 to –20 kg N/ha, Nmin and Ntotal leaching barely increases. It increases more significantly when the N balance becomes zero and positive. The relationship of Nmin leaching (kg/ha) with soil N % (x1) and humus % (x2) is as follows:

The usage of all crop rotations did not ensure a positive N balance in the agroecosystem. The lowest deficit of N was gained with the higher fertilised cereal crop rotation (–14 ± 48.1 kg N/ha/yr). Under less fertilised cereal crop rotation conditions, the N balance deficit increased by 2.2 times and reached -31±56.7 kg N/ha/yr. It was found to be the highest with perennial grass crop rotation (–100 ± 37.4 kg N/ha/yr). Significant N balance differences were found between TI and TIII as well as between TII and TIII (tfact.=2.60>ttheor.95%=2.11 and tfact.=3.02>ttheor.95%=2.07, respectively). Nitrogen balance was not statistically significant between TI and TII (tfact.=–1.27<ttheor.95%=2.07).

An average annual N balance over the eight research years was estimated, excluding N fixation with legume crops and denitrification losses in the soil, as we did not have the tools to record these indicators. An average correlation was determined between fertilisation with N fertiliser (x) and N balance (y):

According to average data, in order to create 1 GJ energy, TIII outstood due to high costs (23.5). Slightly lower costs were reported using TII (18.1), while the lowest costs were using TI (13.9). This is related to the productivity of field crops (dry matter), expressed as GJ/ha. Costs to create vegetative produce include fertiliser, fertilisation, field crops’ agricultural engineering, pesticides, etc. estimated as GJ/ha. LSD05 of named indicators reached 7.6 GJ/ha.

DiscussionTop

During the 2006-2013 period, we were analysing several research questions, such as how Nmin affects the yield, how fertilisation affects Nmin in the soil, the effect of different crop rotations on the annual N balance, and the differences in mineral and Ntotal concentrations in drainage water with fertilisation, as well as how they depend on the soil.

The study showed that a higher fertilised cereal crop rotation had the highest cumulative productivity (1483 GJ/ha), while cereals with less fertilisation and grass crop rotations had a productivity which tended to be lower by 11 to 35%, respectively, suggesting that the crop yield depends on fertilisation rate. Similar results were found by Křen et al. (2014), who analysed spring barley yields applying different N fertilisation rates. The yield increase was determined to be higher using N rate from 0 kg N/ha to 45 kg N/ha, rather than the rate from 45 kg N/ha to 90 kg N/ha. Grover et al. (2009) claimed that crop rotation is a prerequisite for a large crop yield, especially when legume crops are included in the rotation. However, the cultivation of corn monoculture produced the same yield as the crop rotation (Grover et al., 2009). Babuliková (2014) contributed to the findings using wheat in the crop rotation. It was found that in order to obtain the highest wheat yield the best N rate was 96-168 kg N/ha. There was a small difference by fertilising with a rate of 150-200 kg N/ha; however, by increasing the N rate to 240 kg N/ha, the grain yield and N effectiveness dropped. In addition to this, a long-term research revealed that crop rotation and different fertilisation had only a slight effect on yield and soil fertility parameters (Cuvardic et al., 2004). On the other hand, according to Coulter et al. (2011), a high yield in the “extended crop rotations for corn and soybean” can be obtained, even by reducing the fertilisation rates. Additionally, Øgaard (2014) found that excessive fertilisation with N leads to a fall in crop yields, an increase in N losses, and the formation of high N surplus.

In general, the results suggest that higher fertilisation increases mineral N concentration and accumulated amount in the soil. This is evident from our research, which showed that the highest mineral N concentration and accumulated amount in the soil was found in the higher fertilisation cereal crop rotation (respectively 7.6 mg/kg and 76 kg/ha), while the lowestamounts and concentrations were found in the grass crop rotation (3.3 mg/kg and 42 kg/ha), with values that were 57 and 45%, respectively, lowerthan for the higher fertilised cereal crop rotation. Similar results were found by Cuvardic et al. (2004), who observed an increasing trend in soil N (from 0.32 to 0.36%). Additionally, our research showed an increase in humus (4.11%) in perennial grasses crop rotation (Fig. 3). Other researchers found that using N-fixing legumes and red clover increase mineral N concentrations in the soil and might even accumulate the amount of 200 kg/ha of biological N for further crop rotation (Lapinskas, 1998; Riedell et al., 2009). In addition, Röing et al. (2005) claimed that 50% more N mineralises in crop rotations with ley rather than in crop rotations where only cereals are grown. According to Sonneveld & Bouma (2003), on old pastures, reducing N application levels can lower the probability of exceeding the environmental threshold for nitrate by up 20%, whilst hardly affecting N uptake. Köhler et al. (2006), in sand soils, showed that nitrate leaching depended on N fertilisation and different crop rotations as well as mineralisation from soil organic matter. Catch crops were the most efficient way to reduce the NO3 concentrations in the groundwater recharge of sandy soils. These results suggest that the Nmin concentrations in the soil depend not only on the fertilisation rate, but also on the type of crop rotation.

We found that the application of all crop rotations generally resulted in a negative N balance. However, the use of different cereals crop rotations showed a lower N balance deficit (–31 and –14 kg N/ha), while the grass crop rotation resulted in a negative N balance, which increased by 3.2 and 7.1 times and reached -100 kg N/ha. According to Kutra & Aksomaitiene (2003), the most negative N balance (-426 kg N/ha) was found in the perennial grass fields; the most positive balance (+ 174 kg N/ha) was in fields with sugar beets rotation.

The results from a number of countries show that the main environmental problem related to agriculture is N leaching. According to Arregui & Quemada (2006), the amount of leached N highly depends on the amount of it found in drainage and soil before planting crops. When N balance is positive, there is a high possibility that N concentrations in drainage water will be above the maximum allowable values. When N balance is negative, its concentrations are determined to be low. Crop rotation where perennial grasses are grown is regarded to be the best in terms of the impact to the environment, as Ntotal concentrations in drainage water were determined to be the lowest(12.8 mg/L) and not to exceed the required environmental regulations (VZ, 2011). This was also observed in our research, as the Nmin and Ntotal leaching had a tendency to fall in grass crop rotation (7.0 and 8.0 kg/ha, respectively), while in cereal crop rotations with lower and higher fertilisation there was an increasing tendency in N leaching (by 13 to 42%) observed. Additionally, Kutra & Aksomaitiene (2003) showed that in fields where cereals were dominating, N surplus was + 74 kg N/ha. Weber & Kubiniok (2013) also found that in the soils sown with perennial grasses, chemical elements’ leaching losses were lower than in cereals. Thus, due to the fact that grass crop rotations reduce N concentrations in the soil, they also reduce the risk of it leaching into the drainage water. However, there are more factors that affect N leaching, as fertilization rate not higher than plant need and minimized soil tillage systems (Kutra et al., 2006).

Lysimetric research revealed that nitrate concentrations in lysimeter water were highly dependent on fertilisation rates. Moreover, according to a 30-year research work (Adomaitis et al., 2008), by fertilising crop rotation plants (winter wheat, sugar beet, spring barley, annual grasses and perennial grasses), using the fertilisation rates 112 and 224 kg N/ha, the nitrate concentrations observed in lysimeter water at a depth of 0-80 cm were 112.1 and 187.2 mg/L, respectively. Tripolskaja & Verbylienė (2014) determined that by fertilising with N fertiliser (on average 120 kg/ha), compared to unfertilised soil, the average nitrate concentration in lysimeter water was increased by 37.5 to 54.3%. However drainage, due to more dense and more exuberant crops, which reduce precipitation filtration, fell during the year on average by 13.7-16.2%. This is also supported by our research results, which showed that larger biological potential crops reduced drainage by 23% compared to lower biological potential crops. Tiemeyer et al. (2006) determined that N concentrations in drainage water were increasing by increasing flow rates, high loss rates always occurred at high flow rates, and higher loss rates are to be expected in wetter years. This was also confirmed by our research: when the precipitation level was higher, the drainage was also higher. Drainage essentiallydepended on meteorological conditions. Previous research showed that agroecosystem´s hidrological regime is determined not only by natural, but also by human farming activity factors (the crops grown, the farming level). More intensively grown, higher biological potential crops use up more water and thus reduce not only the soil moisture but also drainage (Aksomaitienė et al., 2002). The problematic nature of the research is that the majority of results were obtained from small drainage sites or even by lysimetric research methods. To determine similar patterns by performing tests in large drainage systems is also difficult, not only because of their large amortisation effect (Cornelese et al., 2001), but also because of the research background promiscuity (soil cover variegation, ground water stratification depth, etc.); however, under some circumstances, patterns remain (Sileika & Guzys, 2003).

Both Nmin and Ntotal leaching by drainage, due to the crop productivity and their fertilisation with N fertiliser, change according to y = a – bx / cx2. The smallest N leaching losses were obtained by fertilising 66-71 kg N/ha and crop productivity of 212 GJ/ha. Furthermore, Nmin and Ntotal concentration in drainage water depended on crop rotation, its productivity and the factors determining productivity. Both Nmin and Ntotal concentration in drainage water had a tendency to fall in grass crop rotation (12.1 and 12.8 mg/L respectively). The use of cereal crop rotations with higher and lower fertilisation increased Nmin concentration by 13 and 29%, while total N by 24 and 44%. According to Adomaitis et al. (2008), a larger amount of nitrates is leached when annual crops are grown, compared with perennial grasses. Regular grass sown in arable land and not fertilised with N reduced nitrate concentration to the level of perennial pasture in just one year. The risk of exceeding the maximum allowable limit of nitrates in crop rotation in arable land was significant and did not depend on the level of fertilisation with N (Mašauskas et al., 2006).

Moreover, Beaudoin et al. (2005), determined that nitrate concentrations change subject to the year, crops and soil type; in a deep loamy soil, concentrations fluctuated between 31 mg/L and 92 mg/L in shallow sand. Correa et al. (2005) and Tonitto et al. (2006) observed that the biggest leaching losses occurred in the bare, fallow soil. The lowest concentrations were determined in crop rotation with sugar beet/wheat (38 mg/L), whereas the highest concentrations were found in the pea/wheat crop rotation (66 mg/L) (Beaudoin et al., 2005). Kutra et al. (2006), found that the majority of N was leached from row crops (22.4 kg/ha), slightly less from cereal fields (16.6-18.9 kg/ha) and least of all from grasslands (10.5 kg/ha). On the contrary, Bohm et al. (2009) showed that legume crop cultivation can considerably increase Nmin concentration in the soil solution and leaching as well.

The data imply that the amount of Nmin in the soil is a very labile indicator. Loosening during the mineralisation period and spreading during fertilisation ensure that Nmin is quickly used up by crops. The longer the crop vegetation period and the greater their biological potential, the stronger this preventative effect (Kutra et al., 2006; Askegaard et al., 2011; Doltra et al., 2011).

Finally, as the data collected by many researchers show, the advanced agricultural technologies and efficient application (not the rates) of mineral fertilisers are among the essential preconditions for reduced non-point source water pollution with N (Goulding et al., 2000; Sileika & Guzys, 2003).

The obtained results showed that the application of perennial grasses crop rotation in agricultural fields was the best environmental tool, reducing N migration to drainage. By fertilising with other treatments using higher N rates, mineral N accummulation in the soil increased, stimulating a higher yield and reducing N leaching by drainage. By fertilising with higher N rates than it is neccesary for plants, both mineral and total N leaching to drainage had a tendency to increase.


ReferencesTop

Adomaitis T, Vaisvila Z, Mazvila J, Staugaitis G, Fullen MA, 2008. Influence of mineral fertilizer on nitrogen leaching. Acta Agric Scand B - Soil Plant Sci 58(3): 199-207.doi: DOI: 10.1080/09064710701593012.
Ahlgren S, Baky A, Bernesson S, Nordberg Å, Norén, Hansson PA, 2008. Ammonium nitrate fertilizer production based on biomass. Environmental effects from a life cycle perspective. Bioresour Technol 99(17): 8034-8041http://dx.doi.org/10.1016/j.biortech.2008.03.041.
Aksomaitienė R, Kutra G, Petrokienė Z, 2002. Productivity of crop rotations and their impact on NPK migration in agrosystems. Water Manage Eng Trans 19(41): 3-13. [In Lithuanian].
Alosevičienė G (ed), 1992. Lithuanian climate reference: weather temperature. Lithuanian Hydrometeorological Board, Lithuanian Hydrometeorological Centre, Vilnius, 140 pp. [In Lithuanian].
Arregui LM, Quemada M, 2006. Drainage and nitrate leaching in crop rotation under different N-fertilizer strategies: application of capacitance probes. Plant Soil 288(1-2): 57-69. http://dx.doi.org/10.1007/s11104-006-9064-9.
Askegaard M, Olesen JE, Rasmussen IA, Kristensen K, 2011. Nitrate leaching from organic arable crop rotations is mostly determined by autumn field management. Agr Ecosyst Environ 142(3): 149-160. http://dx.doi.org/10.1016/j.agee.2011.04.014.
Babuliková M, 2014. The influence of fertilization and crop rotation on the winter wheat production. Plant Soil Environ 60(7): 297-302.
Beaudoin N, Saad JK, Van Laethem C, Machet JM, Maucorps J, Mary B, 2005. Nitrate leaching in intensive agriculture in Northern France: effect of farming practices, soils and crop rotations. Agr Ecosyst Environ 111(1-4): 292-310. http://dx.doi.org/10.1016/j.agee.2005.06.006.
Beeston MP, Pohar A, Elteren JT van, Plazl I, Štejkovec Z, Veber M, Glass HJ, 2010. Assessment of physical leaching processes of some elements in soil upon ingestion by continuous leaching and modelling. Environ Sci Technol 44(16): 6242-6248. http://dx.doi.org/10.1021/es1006725.
Bohm Ch, Landgraf D, Makeschin F, 2009. Effects of two contrasting agricultural land–use practice on nitrogen leaching in a sandy soil of Middle Germany. J Plant Nutr Soil Sc 172(3): 408-417. http://dx.doi.org/10.1002/jpln.200700104.
Buckley C, Carney P, 2013. The potential to reduce the risk of diffuse pollution from agriculture while improving economic performance at farm level. Environ Sci Policy 25: 118-126. http://dx.doi.org/10.1016/j.envsci.2012.10.002.
Cabrera VE, Jagtab SS, Hilderbrand PE, 2007. Strategies to limit (minimize) nitrogen leaching on dairy farms driven by seasonal climate forecasts. Agr Ecosyst Environ 122(4): 479-489. http://dx.doi.org/10.1016/j.agee.2007.03.005.
Cornelese AA, Linden AMA van der, Verschoor AJ, 2001. Pesticide behaviour in soils and water. British Crop Protection Council, Brighton, UK, Nov 13-15. pp: 95-100.
Correa RS, White RE, Weatherley AJ, 2005. Modelling the risk of nitrate leaching from two soils amended with five different biosolids. Rev Bras Cienc Solo 29(4): 619-626. http://dx.doi.org/10.1590/S0100-06832005000400014.
Coulter JA, Sheaffer CC, Wyse DL, Haar MJ, Porter PM, Quiring SR, Klossner LD, 2011. Agronomic performance of cropping systems with contrasting crop rotations and external inputs. Agron J 103(1): 182-192. http://dx.doi.org/10.2134/agronj2010.0211.
Cuvardic M, Tveitnes S, Krogstad T, 2004. Long-term effects of crop rotations and different fertilization systems on soil fertility and productivity. Acta Agric Scand B - Soil Plant Sci 54(4): 193-201. http://dx.doi.org/10.1080/09064700410030258.
Doltra J, Lægdsmand M, Olesen JE, 2011. Cereal yield and quality as affected by nitrogen availability in organic and conventional arable crop rotations: A combined modeling and experimental approach. Eur J Agron 34: 83-95. http://dx.doi.org/10.1016/j.eja.2010.11.002.
Egner H, Riehm H, Domingo WR, 1960. Untersuchungen uber die chemische Bodenanalyse als Grundlage fur die Burtreilung des Nahrstoffzus-tandes der Boden. II. Chemische Extractionsmethoden zur Phosphor uns Kalumbestimmung. Kunliga Landbouksshogskolans Annaler 26: 199-215.
Goulding KWT, Poulton PR, Webster CP, Howe MT, 2000. Nitrate leaching from the Broadbalk Wheat Experiment, Rothamsted, UK, as influenced by fertilizer and manure inputs and the weather. Soil Use Manage 16(4): 244-250. http://dx.doi.org/10.1111/j.1475-2743.2000.tb00203.x.
Grover KK, Karsten HD, Roth GW, 2009. Corn grain yields and yield stability in four long-term cropping systems. Agron J 101(4): 940-946. http://dx.doi.org/10.2134/agronj2008.0221x.
Hongxu W, Chengyang X, Lvyi M, Lini J, Yuan CH, 2011. Leaching of mineral nitrogen and available phosphorus during culture of Bareroot Larix olgensis seedlings under native fertilizer management in Northeastern China. Proc Int Conf on computer distributed control and intelligent environment monitoring, Changsha (China), Feb 19-20. pp: 1247-1252.
Hossain TS, Sugimoto H, Ueno H, Rafiul Huque SM, 2007. Adoption of organic rice for sustainable development in Bangladesh. J Org Syst 2: 27-37.
Hyytiäinen K, Niemni JK, Koikkalainen K, Palosuo T, Salo T, 2011. Adaptive optimization of crop production and nitrogen leaching abatement under yield uncertainty. Agr Syst 104(8): 634-644. http://dx.doi.org/10.1016/j.agsy.2011.06.006.
ISO 11261:1995. Soil quality. Determination of nitrogen. Modified Kjeldahl method. Available in http://www.iso.org/iso/home/store/catalogue_detail.htm?csnumber=19239. [2 December 2014].
ISO 13395:1996. Water quality. Determination of nitrate nitrogen and nitrate nitrogen and the sum of both by flow analysis (CFA and FIA) and spectrometric detection. Available in http://www.iso.org/iso/catalogue_detail.htm?csnumber=21870. [2 December 2014].
ISO/TS 14256-1:2003. Soil quality. Determination of nitrate, nitrite and ammonium in field moist soils by extraction with potassium chloride solution. Part 1: Manual method. Available in http://www.iso.org/iso/catalogue_detail.htm?csnumber=36706. [2 December 2014].
ISO 5983-1:2005/Cor 1:2008. Animal feeding stuffs. Determination of introgen content and calculation of crude protein content. Part 1: Kjeldahl method. Available in http://www.iso.org/iso/catalogue_detail.htm?csnumber=52374. [2 December 2014].
ISO 10390:2005. Soil quality. Determination of pH. Available in http://www.iso.org/iso/catalogue_detail.htm?csnumber=40879. [2 December 2014].
ISO 11732:2005. Water quality. Determination of ammonium nitrogen. Method by flow analysis (CFA and FIA) and spectrometric detection. Available in http://www.iso.org/iso/catalogue_detail.htm?csnumber=38924. [2 December 2014].
Jankauskas B, Jankauskienė G, Švedas A, 2000. The comparison of yield energetic estimation calculation methods. Zemdirbyste 72: 239-251. [In Lithuanian].
Juozapaitis A, Berankienė L, Zelionkienė V, 1995. Hydrometric Annual. Drainage 1962-1990. Parts I and II. Kėdainiai, Vilainiai. [In Lithuanian].
Klaushal SS, Groffman PM, Band LE, Elliot EM, Shields CA, Kendall C, 2011. Tracking nonpoint source nitrogen pollution in human-impacted watersheds. Environ Sci Technol 45 (19): 8225-8232. http://dx.doi.org/10.1021/es200779e.
Köhler K, Duynisveld WHM, Böttcher J, 2006. Nitrogen fertilization and nitrate leaching into groundwater on arable sandy soils. J Plant Nutr Soil Sci 169(2): 185-195. http://dx.doi.org/10.1002/jpln.200521765.
Kokkora MI, Hann MJ, 2007. Crop production and nitrogen leaching resulting from biowaste and onion compost amended sand. Proc XI Int Symp on waste management and land till, S. Margherita di Pula – Cagliari, Sardinia (Italy), Oct 1-5. pp: 258-263.
Komonweeraket K, Benson CH, Edil TB, Blem WF, 2011. Leaching behaviour and mechanisms controlling the release of elements from soil stabilized with fly ash. Proc Conf Geo-Frontiers 2011: Advances in Geotechnical engineering. Dallas, TX, USA, Mar 13-16. pp: 1101-1110.
Korkutis P (ed), 1991. Lithuanian climate reference: precipitation. Lithuanian Hydrometeorological Board, Lithuanian Hydrometeorological Centre, Vilnius, 122 pp. [In Lithuanian].
Křen J, Klem K, Svobodová I, Míša P, Neudert L. 2014. Yield and grain quality of spring barley as affected by biomass formation at early growth stages. Plant Soil Environ 60(5): 221-227.
Kutra G, Aksomaitiene R, 2003. Use of nutrient balances for environmental impact calculations on experimental field scale. Eur J Agron 20(1-2): 127-135. http://dx.doi.org/10.1016/S1161-0301(03)00082-0.
Kutra G, Gaigalis K, Šmitienė A, 2006. Lands use influence on nitrogen leaching and options for pollution mitigation. Zemdirbyste-Agr 93 (4): 119-129.
Lapinskas E, 1998. Biological nitrogen fixation and nitragine. Lithuanian Institute of Agriculture, Akademija. 218 pp. [In Lithuanian].
Larsson M, Granstedt A, 2010. Sustainable governance of the agriculture and the Baltic Sea. Agricultural reforms, food production and curbed eutrophication. Ecol Econ 69(10): 1943-1951. http://dx.doi.org/10.1016/j.ecolecon.2010.05.003.
Liekis A (ed), 2001. Soils of Lithuania. Lithuanian Academy of Sciences, Spauda, Vilnius. 1244 pp.
Lundberg JO, Weitzberg E, Gladwin MT, 2008. The nitrate-nitrite-nitric oxide pathway in physiology and therapeutics. Nature Rev Drug Discov 7: 156-167. http://dx.doi.org/10.1038/nrd2466.
Martínez-Blanco J, Rieradevall J, Antón A, Mu-oz P, 2014. Multifunctionality-solving approaches of compost application in crop rotations. J Clean Prod 64: 384-396. http://dx.doi.org/10.1016/j.jclepro.2013.08.014.
Mašauskas V, Antanaitis Š, Lazauskas S, Mašauskienė A, 2006. Content of nitrates in drainage ang groundwater from permanent pasture, grassland and arable crop rotation. Ekologija 4: 83-88.
Menz FC, Seip HM, 2004. Acid rain in Europe and the United States: an update. Environ Sci Policy 7(4): 253-265. http://dx.doi.org/10.1016/j.envsci.2004.05.005.
Nakamura K, Harter T, Hirono Y, Horino H, Mitsuno T, 2004. Assessment of root zone nitrogen leaching as affected by irrigation and nutrient management practices. Vadose Zone J 3: 1353-1366. http://dx.doi.org/10.2136/vzj2004.1353.
Oenema O, Witzke HP, Klimont Z, Lesschen JP, Velthof GL, 2009. Integrated assessment of promising measures to decreace nitrogen losses in agriculture in EU-27. Agr Ecosyst Environ 133: 280-288. http://dx.doi.org/10.1016/j.agee.2009.04.025.
Øgaard AF, 2014. Nitrogen balance and nitrogen use efficiency in cereal production in Norway. Acta Agric Scand B - Soil Plant Sci 63(2): 146-155. http://dx.doi.org/10.1080/09064710.2013.843718.
Payen S, Mens CB, Perret S, 2015. LCA of local and imported tomato: an energy and water trade-off. J Clean Prod 87: 139-148. http://dx.doi.org/10.1016/j.jclepro.2014.10.007.
Powlson DS, Addiscott TM, Benjamin N, Cassman KG, de Kok TM, van Grinsven H, L’hirondel JL, Avery AA, van Kessel C, 2008. When does nitrate become a risk for humans? J Environ Qual 37: 291-295. http://dx.doi.org/10.2134/jeq2007.0177.
Riedell WE, Pikul JL, Jr, Jaradat AA, Scumacher TE, 2009. Crop rotation and nitrogen input effects on soil fertility, maize mineral nutrition, yield, and seed composition. Agron J 101(4): 870-879. http://dx.doi.org/10.2134/agronj2008.0186x.
Röing K, Andren O, Mattsson L, 2005. Long-term management effect on plant N uptake and topsoil carbon levels in Swedish long-term field experiments: cereals and ley, crop residue treatment and fertilizer N application. Acta Agric Scand B - Soil Plant Sci 55(1): 16-22. http://dx.doi.org/10.1080/09064710510008739.
Rosenstock TS, Liptzin D, Six J, Tomich TP, 2013. Nitrogen fertilizers in California: Assessing the data, trends and a way forward. Calif Agr 67(1): 68-79. http://dx.doi.org/10.3733/ca.E.v067n01p68.
Salvioni C, Jaskho N, Kauko K, 2011. Adaptive optimisation of crop production and nitrogen leaching abatement under yield uncertainty. Proc XVIII Annu Conf of European association of environmental and resource economists, Rome (Italy), June 29-July 2. pp: 283-291.
Savci S, 2012. An agricultural pollutant: chemical fertilizer. Int J Environ Sci Dev 3: 77-80. http://dx.doi.org/10.7763/ijesd.2012.v3.191.
Shibusawa S, Hache C, 2009. Data collection and analysis methods for data from field experiments. Systems analysis and modelling in food and agriculture. In: UNESCO’s Digital Encyclopaedia of Life Support Systems (Ting KC, Fleisher DH, Rodriquez LF, eds). EOLSS Publishers, Oxford (UK), pp: 312-338.
Sileika AS, Guzys S, 2003. Drainage runoff and migration of mineral elements in organic and conventional cropping systems. Agronomie 23: 633-641. http://dx.doi.org/10.1051/agro:2003428.
Smith VH, Schindler DW, 2009. Eutrophication science: where do we go from here? Trends Ecol Evol 24(4): 201-207. http://dx.doi.org/10.1016/j.tree.2008.11.009.
Sogbedji JM, van Es HM, Klausner SD, Bouldin DR, Cox WJ, 2001. Spatial and temporal processes affecting nitrogen availability at the landscape scale. Soil Till Res 58: 233-244. http://dx.doi.org/10.1016/S0167-1987(00)00171-9.
Sonneveld MPW, Bouma J, 2003. Effects of combinations of land use history and nitrogen application on nitrate concentration in the groundwater. NJAS-Wagen J Life Sc 51(1-2): 135-146. http://dx.doi.org/10.1016/S1573-5214(03)80030-1.
Thorup-Kristensen K, Dresbøll DB, Kristensen HL, 2012. Crop yield, root growth, and nutrient dynamics in a conventional and three organic cropping systems with different levels of external inputs and N re-cycling through fertility building crops. Eur J Agron 37(1): 66-82. http://dx.doi.org/10.1016/j.eja.2011.11.004.
Tiemeyer B, Kahle P, Lennartz B, 2006. Nutrient losses from artificially drained catchments in North-Eastern Germany at different scales. Agr Water Manage 85: 47-57. http://dx.doi.org/10.1016/j.agwat.2006.03.016.
Tonitto C, David MB, Drinkwater LA, 2006. Replacing bare follows with cover crops in fertilizer-intensive cropping systems: a meta-analysis of crop yield and N dynamics. Agr Ecosyst Environ 112 (1): 58-72. http://dx.doi.org/10.1016/j.agee.2005.07.003.
Tripolskaja L, Verbylienė I, 2014. The effect of different forms of nitrogen fertilizers on nitrogen leaching. Zemdirbyste-Agr 101(3): 243–248. http://dx.doi.org/10.13080/z-a.2014.101.031.
Tyurin IV, 1937. Soil organic matter and its role in pedogenesis and soil productivity. Study of soil humus. Sel’skozgiz, Moskva. 491 pp. [In Russian].
Ulen B, Johansson G, 2009. Long-term nutrient leaching from a Swedish arable field with intensified crop production against a background of climate change. Acta Agric Scand B - Soil Plant Sci 59: 157-169. Doi: 10.1080/09064710802040541.
Ulen B, Brommsen C, Johansson G, Torstensson G, 2012. Trends in nutrient concentrations in drainage water from single fields under ordinary cultivation. Agr Ecosyst Environ 151: 61-69. http://dx.doi.org/10.1016/j.agee.2012.02.005.
Vadiunina AF, Korcagina ZA, 1986. Research methods of physical soil characteristics. Moskva, Agropromizdat, 416 pp. [In Russian].
VZ, 2011. On the ratification of the farming entities’ environment monitoring regulations. Lithuanian Minister of the Environment order No. Dl-911 as of 25th November 2011 On the change of Minister of the Environment order No. D1-546 as of 16th September 2009. Valstybes Zinios No. 148-6962. [In Lithuanian].
Weber G, Kubiniok J, 2013. Investigation of spring waters to assess the geo-ecological state of landscape and soil. J Agr Sci Technol A 3: 592-602.