Grazing under Irrigation A ﬀ ects N 2 O-Emissions Substantially in South Africa

: Fertilized agricultural soils serve as a primary source of anthropogenic N 2 O emissions. In South Africa, there is a paucity of data on N 2 O emissions from fertilized, irrigated dairy-pasture systems and emission factors (EF) associated with the amount of N applied. A ﬁrst study aiming to quantify direct N 2 O emissions and associated EFs of intensive pasture-based dairy systems in sub-Sahara Africa was conducted in South Africa. Field trials were conducted to evaluate fertilizer rates (0, 220, 440, 660, and 880 kg N ha − 1 year − 1 ) on N 2 O emissions from irrigated kikuyu–perennial ryegrass ( Pennisetum clandestinum–Lolium perenne ) pastures. The static chamber method was used to collect weekly N 2 O samples for one year. The highest daily N 2 O ﬂuxes occurred in spring (0.99 kg ha − 1 day − 1 ) and summer (1.52 kg ha − 1 day − 1 ). Accumulated N 2 O emissions ranged between 2.45 and 15.5 kg N 2 O-N ha − 1 year − 1 and EFs for mineral fertilizers applied had an average of 0.9%. Nitrogen in yielded herbage varied between 582 and 900 kg N ha − 1 . There was no positive e ﬀ ect on growth of pasture herbage from adding N at high rates. The relationship between N balance and annual N 2 O emissions was exponential, which indicated that excessive fertilization of N will add directly to N 2 O emissions from the pastures. Results from this study could update South Africa’s greenhouse gas inventory more accurately to facilitate Tier 3 estimates. as function of N input and N balance . Statistical signiﬁcance of the tested factors, comparisons of means, and regression equations (intercepts and < 0.05.


Introduction
The main greenhouse gas (GHG) contributors towards net global warming potential from the agricultural sector are carbon dioxide (CO 2 ), methane (CH 4 ), and nitrous oxide (N 2 O) [1]. N 2 O is released into the atmosphere during the processes of denitrification (reduction of nitrate (NO 3 − ) to di-nitrogen gas (N 2 ) by soil microbes) under anaerobic conditions and nitrification (oxidation from NH 4 + to NO 3 − ) under aerobic conditions [2]. N 2 O has a 121-year atmospheric life span and a global warming potential of 265 times that of CO 2 compared over a 100-year period [3]. This makes mitigation strategies of N 2 O critical to control GHG emissions from the agricultural sector. Globally, agricultural soils serve as a primary source of anthropogenic N 2 O emissions [4,5]. The GHG national inventory report of South Africa [6] stated that the energy (10.7%) and Agriculture, Forestry and Other Land Use (AFOLU) (84.5%) sectors were the largest contributors to the total N 2 O emissions in 2015. It is predicted that South African agricultural soils are responsible for 28% of N 2 O emissions from the AFOLU sector as a result of fertilizer application, urea for animal feed, and animal excreta [7].
More consumers are becoming aware and concerned about the origin and impact that products have on the environment, which is expected to increase the demand for food products that generate low GHG emissions [30]. On the other hand, it should be important to consumers to also reduce food wastage and to encourage diets with a low N footprint which could help reduce the global demand and usage by using less fertilizers to grow food [15]. This would be in line with the European Nitrogen Assessment (ENA) which creates a better public awareness by identifying the challenges and threats associated with N pollution [29] and how to approach them as mitigation options.
In view of the arguments above, it is important to evaluate dairy-pasture systems in terms of associated greenhouse gas emissions. In the current study we aimed to obtain quantitative values for N 2 O emissions from pasture-based dairy production systems in South Africa. We evaluated the effect of different N fertilizer rates on N 2 O emissions. Subsequently, the emission factors (EF) associated with different N fertilizer rates were calculated and compared to the recommended IPCC values [1]. The EFs associated with different N fertilization rates will allow for development of strategies and policies which could assist in adaptation of farming practices that will reduce anthropogenic GHG emissions.
Accordingly, in this paper, we present results from field trials to address these important questions: (1) What would be the response of N 2 O emissions under managed soils in the southern Cape region of South Africa to N fertilization, under irrigation and intensive grazing practices, because of N levels exceeding plant requirements? (2) Can high stocking rates on intensively managed, highly fertilized and irrigated dairy pastures lead to high amounts of N returned through excreta to the soil and result in a high N surplus which could underestimate predicted N 2 O emissions? (3) Could the EFs, calculated from N 2 O emission values from intensively managed, grazed systems under irrigation, accurately predict N 2 O emissions compared to the suggested EF as set by the IPCC Tier 1 default value for N fertilizers?

Experimental Site Description
A field experiment was conducted at Outeniqua Research Farm (33 • 58 38" S; 22 • 25 16 E; 201 m.a.s.l.) of the Western Cape Department of Agriculture, South Africa. The research farm is located near the city of George in the southern Cape region of South Africa. The area has a temperate climate with mean monthly temperatures ranging between 7-18 • C and 15-25 • C in winter and summer, respectively. Rainfall is evenly distributed throughout the year and has a mean annual average precipitation of 728 mm. The soil type on the experimental site can be classified as a Podzol (IUSS Working Group 2015) which is locally known as a Witfontein soil form (Soil Classification Working Group 1991). Soil properties in different depths are shown in Table 1. Soil bulk densities for the 0-10, 10-20, and 20-30 cm layers were 1.3, 1.5, and 1.5 g cm −3 , respectively [31]. Winter was defined as months June-August, spring as September-November, summer as December-February, autumn as March-May.

Experimental Layout and Treatments
An experiment was laid out as a randomized block design with five N fertilizer rates as treatments, replicated in four blocks. Plots were 15 × 15 m. Different N fertilization treatments in the form of limestone ammonium nitrate (LAN) were applied in fixed rates of 0, 20, 40, 60, and 80 kg N ha −1 (respectively defined as N0, N20, N40, N60, and N80) approximately four days after each grazing. The study took place over a period of 372 productive days from June 2018 to June 2019. A total of 12 grazing cycles could be managed and therefore 11 fertilization events applied. Consequently, the total amounts of N fertilizer applied per productive year were 0, 220, 440, 660, and 880 kg N ha −1 for the treatments.

Pasture and Grazing Management
The experimental site consisted of kikuyu over-sown with ryegrass and managed under no-tillage practices. A previous trial, two years prior to the current study (2016-2017), was carried out using the same N fertilization rates as in the current study applied to the same experimental plots for the respective treatments. Thus, no artifact effects of previous N fertilizer management practices were expected. More information can be found in Viljoen et al. (2020) [27]. During over-sowing of ryegrass, pastures were first grazed to a height of ca. 50 mm above ground level and mulched afterwards. In autumn 2018 the kikuyu was over-sown with annual Italian ryegrass (Lolium multiflorum ssp. italicum) cv. Fox and in autumn 2019 the kikuyu pasture was over-sown with annual Italian ryegrass cv. Thabu Plus. The botanical composition of the pasture is mainly driven by seasonal growth with kikuyu dominating in summer and ryegrass in winter and it was comparable between treatments. However, the control plots (N0) had a higher component of unsown legumes (Trifolium spp.) compared to the other treatments and they contributed 14% to the sward. It is known that the contribution of legumes in a sward increases as the added amount of fertilizer decreases [27].
Pastures were under permanent sprinkler irrigation with 15 m spacing between sprinklers. Irrometer tensiometers (Calafrica SA, Nelspruit, South Africa) were installed at a depth of 15 cm and irrigation was scheduled to maintain a soil water matrix potential between −25 and −10 kPa. Pasture was intensively strip-grazed by Jersey cows with grazing cycles between ca. 28 days in summer and ca. 35 days in winter. Cows were allowed to voluntarily graze the experimental plots and were not only restricted to specific treatments. The cows grazed pastures to ground level so that residual effects in determining herbage yield were considered minimal. Cows were allocated between 8 and 10 kg DM pasture with an additional concentrate feeding in the milking parlor of 5-6 kg concentrate. The concentrate fed consisted of 12-15% crude protein (CP), a metabolizable energy content of 11 MJ kg −1 , 0.4% P, and 1% Ca, all on a DM basis.

N 2 O Measurements
N 2 O emissions were captured using the static chamber method [33]. Basal rings made from polyvinyl chloride (PVC) (55 cm diameter and 15 cm height) were installed in the soil two weeks prior to the onset of the experiment to a soil depth of 10 cm in each plot. This was done to avoid high measured gas fluxes which might be captured shortly after establishment. This period would also allow for the recovery of roots which might have been damaged during the insertion process [34]. The PVC chambers (60 cm diameter and 35 cm height) were deployed onto the basal rings during measurement and secured with a resistance rubber band to make an air-tight seal. The first measurement was taken immediately after deployment of chambers onto basal rings, and thereafter at 20, 40, and 60 min. Samples were obtained using a 30 mL syringe through a septum cap located at the top of the chamber and were immediately transferred into 12 mL pre-evacuated exetainers (Labco, High Wycombe, UK). Gas samples were analyzed for N 2 O through a gas chromatograph (SCION 456-GC, Bruker, Leiderdorp, The Netherlands) equipped with a 63 Ni electron-capture-detector using He as carrier gas. Samples were injected using an autosampler (model 271 LH, Gilson Inc., Middleton, WI, USA). Data were processed using the software Compass CDS (Version 3.0.1). N 2 O emissions were measured from the beginning of June 2018 and continued until the beginning of June 2019. Gas measurements were performed on a weekly basis unless in the case of a fertilization event, where it was done in three consecutive days after fertilization. Samples were collected according to the described interval in the morning at 11:00 a.m.

Flux Calculations
The N 2 O flux was calculated based on the increased chamber headspace concentration within 60 min of chamber deployment. Flux calculations were based on the calculation done by Beetz et al. [35] and can be found in Equation (1).
where F is the calculated flux (µg N 2 O-N m −2 h −1 ), k is a unit conversion factor (1.25 kg N m −3 for N 2 O), T is the mean temperature inside the chamber headspace ( • C), V is the chamber volume (m 3 ), A is the collar area (m 2 ), and ∆c/∆t is the change in concentration within the chamber headspace over time (N 2 O: ppb h −1 ). Fluxes were calculated for each treatment and replicated by linear regression between measured N 2 O concentrations and time. The accumulated N 2 O emissions were calculated by linear interpolation between the measured daily fluxes for each replication. The accumulated period included 12 months for all treatments.

Nitrogen Balance
For the calculation of N balance, the N excreted by the dairy animals was considered. To do so, the N excreted by grazing animals was calculated using Equation (2) as described by Nennich et al. (2005) [36].
where N exc is the N excreted by grazing animals (g N day −1 ), DMI is the dry matter intake of the cows (15 kg day −1 ), CP diet is the crude protein content in the diet (0.2 g g −1 of DM), and BW is the body weight of the animal (400 kg). During the study period, 12 grazing days with ≈114 animals were carried out leading to a total amount of excreted N of 450 kg N ha −1 year −1 . N input was calculated using a simple equation, adding the amount of N excreted by grazing animals and N applied as mineral fertilizer from different treatment rates, using Equation (3).
where N input (kg N ha −1 ) is the amount of N applied to different treatments, N exc is the nitrogen excreted by the grazing animals (kg N ha −1 ), and N fert is the different N fertilization rates applied (kg N ha −1 ). The N balance for all treatment rates was calculated using the different N inputs and N outputs as shown in Equation (4). N yield in the harvested biomass was considered as the nitrogen output. (4) where N balance is the difference of N (kg N ha −1 ) utilized/not utilized by subtracting N yield (kg N ha −1 ) from the N input (kg N ha −1 ).

Regression Analysis and Emission Factors
Both the linear and nonlinear regressions were used to predict annual N 2 O emissions (E1-E4) as a function of N input and N balance .
According to the Tier 1 approach of the IPCC guidelines the N 2 O-N emitted from soil can be calculated as a percentage (0.01) of N applied from mineral fertilizers (Equation (5)) (IPCCmin). However, as grazing was a part of our experiments a different EF had to be adopted for animal excreta N (0.02). The sum of N 2 O-N derived from mineral fertilizer application and excreta were defined as IPCC min+exc and calculated using Equation (6).
Atmosphere 2020, 11, 925 6 of 19 Moreover, we compared IPCC default values and regression-based predictions with an EF derived from our experimental data based on Equation (1) (EF data ); whereby EF data was calculated under the assumption that emissions from excreta of grazing animals are captured by the control treatment and after deduction EF data is referring to the amount of mineral fertilizer applied: where N 2 O fert is the accumulated N 2 O-N emissions from the different fertilization rates used in the current study (kg N 2 O-N ha −1 ), N 2 O control is the accumulated N 2 O-N emissions from the treatment which received no mineral fertilizer (kg N 2 O-N ha −1 ) and the mineral fertilizer applied is the different fertilization rates applied (kg N ha −1 ). The EF data was calculated for each treatment and block. An average value was used to predict annual N 2 O emissions from the EF data and the annual N-rate of mineral fertilizer applied in order to compare the different N 2 O prediction approaches. Based on the different approaches the measured values were compared with the predicted values. The performance of the different approaches was evaluated on the basis of the coefficient of determination (R 2 ), root mean squared error (RMSE), and Nash-Sutcliffe model efficiency coefficient (NSE) [37].

Herbage Yield and Forage Quality
Herbage yield was sampled by cutting the pasture prior to grazing. Metal rings with an area of 0.0985 m 2 were used to make 30 mm above ground cuttings. Ten ring samples per plot were cut by hand shears and collected in a bag. The rings were placed randomly within the different plots. Samples were placed in an oven (Scientific Manufacturing cc. (SMC) Oven, ODS, 1400 L, Killarney Gardens, South Africa) to dry at 60 • C for 72 h. The DM content and herbage yield (kg DM ha −1 ) could then be determined. Dried pasture herbage samples from each plot were dried and milled (SWC Hammer mill, 1 mm sieve), then photometrically analyzed using a Gallery™ Discrete Analyzer (Thermo Fisher Scientific, Waltham, MA, USA) to determine the N content of the pasture herbage from different plots. The N yield (kg N ha −1 year −1 ) could then be determined for the different fertilizer treatments.

Soil Sampling
Prior to grazing, soil samples from each plot were taken on the same day as herbage samples, which leads to a sample interval of approximately once per month. Three soil cores ( 7.5 cm) per plot were collected one day prior to every grazing to a 10 cm depth and seasonally sampled at three depths; 0-10, 10-20, and 20-30 cm. Subsamples were pooled and composited. Soil wet weight was recorded and then soil was left to dry in an oven (ODS, 1400 L, SMC Oven, Killarney Gardens, South Africa) at 30 • C for ca. 7 days. The weight was recorded after the soil had dried and gravimetric soil water (SWC) content determined.
The volumetric water content as a fraction of the total pore space was determined and used to calculate the water filled pore space (WFPS). A particle density of 2.65 g cm −3 was assumed.
Salicylic acid [38] and indophenol-blue [39] methods were used to analyze for NO 3 − -N and NH 4 + -N, respectively. N min content was defined as the sum of the NO 3 − -N and NH 4 + -N contents. N min (kg ha −1 ) was calculated using the appropriate bulk densities [31] for the corresponding soil depths.

Weather Data, Soil Temperature, and Soil Water
Air temperature was logged using a Decagon Em50 Weather Logger (Decagon Devices Inc., Pullman, Washington, DC, USA) which recorded daily temperatures in 15-min intervals. In addition, soil temperature was logged using a DFM capacitance probe (DFM Software Solutions cc, 2012) located in the experimental site. Rain gauges were used to manually record rainfall during the experimental period. Monthly rainfall and mean daily temperatures for the trial period are shown in Figure 1.

Weather Results and Environmental Variables
The mean monthly rainfall for the trial period was 49 ± 27.9 mm, compared to the long-term average of 60 ± 13.7 mm (30 years) (Figure 1). The winter of 2018 had below average rainfall and the highest rainfall occurred during spring in September, 101 mm ( Figure 1). January had a low average rainfall (24 mm) compared to the long-term average. Daily temperatures recorded during the trial period ranged from 7.9 to 29.2 °C with an average of 16.1 °C. As a result of temperatures being higher, irrigation applied to pasture was higher in spring and summer and had a mean value of 11 mm over the trial period. The mean values and standard error of the mean for Nmin, NH4 + -N, and NO3 − -N during the entire trial period are presented in Table 2 Supplementary Table S1). The mean annual Nmin ranged between 23.06 ± 1.19 (N20) and 60.02 ± 7.53 kg N ha −1 (N80) and the N60 and N80 treatments differed (p < 0.05) from the N0, N20, and N40 treatments (0-10 cm soil depth). The mean NH4 + -N ranged between 16.42 ± 0.87 (N80) and 23.99 ± 1.93 kg ha −1 (N80). Differences (p < 0.05) in NH4 + -N were found

Statistical Analyses
The statistical software R 3.6.1 (2019) [40] was used to analyze the data using the packages "nlme" [41] and "multcomp" [42]. The data evaluation started by defining an appropriate statistical mixed model [43,44]. Data distribution was visually assumed to be normal and heteroscedastic [45] with regard to the fertilization treatments. These assumptions were based on a graphical residual analysis. The statistical model included treatment as fixed factor. Block was regarded as the random factor. Based on this model a one-way analysis of variance (ANOVA) was conducted to test the hypothesis of the experiment. Furthermore, multiple contrast tests (e.g., see Bretz et al. 2011 [46]) were implemented in order to compare the several levels of the tested fertilization treatments. In addition, simple linear and nonlinear regression models were developed to predict N 2 O emissions as a function of N input and N balance . Statistical significance of the tested factors, comparisons of means, and regression equations (intercepts and slopes) were considered when p < 0.05.

Weather Results and Environmental Variables
The mean monthly rainfall for the trial period was 49 ± 27.9 mm, compared to the long-term average of 60 ± 13.7 mm (30 years) (Figure 1). The winter of 2018 had below average rainfall and the highest rainfall occurred during spring in September, 101 mm ( Figure 1). January had a low average rainfall (24 mm) compared to the long-term average. Daily temperatures recorded during the trial period ranged from 7.9 to 29.2 • C with an average of 16.1 • C. As a result of temperatures being higher, irrigation applied to pasture was higher in spring and summer and had a mean value of 11 mm over the trial period.
The mean values and standard error of the mean for N min , NH 4 + -N, and NO 3 − -N during the entire trial period are presented in Table 2. Soil inorganic N, NH 4 + -N, and NO 3 − -N were more concentrated at a shallow soil depth (0-10 cm). The minimum N min values for all treatments occurred in May 2019

Herbage and Nitrogen Yield
The herbage and N yield are shown in Table 3. Different fertilizer rates had an effect on herbage yield (p < 0.05). Herbage yield was within the range of 18.7 ± 1.1 (N0) to 21.5 ± 0.7 t DM ha −1 (N80). The only significant difference was found between the N80 treatment compared with the N40 and N20 treatments (p < 0.05). The N0, N20, N40, and N60 treatments did not differ (p > 0.05) from each other. Table 3. Means (± standard error) of herbage yield, N-yield, N input , N balance , and measured annual N 2 O-N emissions. The N0, N20, N40, N60, and N80 refer to the fertilizer rates used as treatments and were 0, 220, 440, 660, and 880 kg N ha −1 year −1 , respectively. N input and N balance were calculated using Equations (3) and (4) and were not independent from treatments, and therefore not subjected to statistical analysis. Means in the same row with no common superscript differed significantly (p < 0.05). The N yield differed significantly (p < 0.05) between treatments N20, N40, N60, and N80. No differences (p > 0.05) were found between the N0 and N20 treatments. The N yield from the different treatment plots increased as the amount of fertilizer applied increased and was in the range between 582 ± 38.9 (N0) and 900 ± 35.3 kg N ha −1 (N80).

Daily N 2 O Fluxes
The daily N 2 O fluxes over the trial period (June 2018 to June 2019) for the different N fertilizer treatments are shown in Figure 2. Variation in N 2 O emissions within different weeks were observed, especially in the weeks where N was added as mineral fertilizer. We observed N 2 O emissions to be episodic with small fluxes throughout the measurement period, except in warmer months where higher fluxes were observed. The minimum daily N 2 O-N fluxes were close to zero (0.01 kg N 2 O-N ha −1 day −1 ) in N0, whereas it ranged between 0.01 and 1.47 kg N 2 O-N ha −1 day −1 for the different N fertilizer rates during single measurements on sampling days. The highest N 2 O fluxes were observed for treatments N60 and N80, which occurred during spring and summer and reached peak fluxes of 0.99 and 1.47 kg N 2 O-N ha −1 day −1 , respectively. The lowest N 2 O fluxes over all treatments were observed during winter and oscillated between 0.01 and 0.49 kg N 2 O-N ha −1 day −1 . N 2 O fluxes peaked after fertilizer application in all treatments over the measured trial period (Figure 2). Higher N 2 O fluxes were observed in summer when there were fewer days between rotational grazing events (ca. 28 days) compared with winter grazing (ca. 35 days).
Atmosphere 2020, 11, x FOR PEER REVIEW 9 of 19 especially in the weeks where N was added as mineral fertilizer. We observed N2O emissions to be episodic with small fluxes throughout the measurement period, except in warmer months where higher fluxes were observed. The minimum daily N2O-N fluxes were close to zero (0.01 kg N2O-N ha −1 day −1 ) in N0, whereas it ranged between 0.01 and 1.47 kg N2O-N ha −1 day −1 for the different N fertilizer rates during single measurements on sampling days. The highest N2O fluxes were observed for treatments N60 and N80, which occurred during spring and summer and reached peak fluxes of 0.99 and 1.47 kg N2O-N ha −1 day −1 , respectively. The lowest N2O fluxes over all treatments were observed during winter and oscillated between 0.01 and 0.49 kg N2O-N ha −1 day −1 . N2O fluxes peaked after fertilizer application in all treatments over the measured trial period (Figure 2). Higher N2O fluxes were observed in summer when there were fewer days between rotational grazing events (ca. 28 days) compared with winter grazing (ca. 35 days).

Accumulated N 2 O-Losses
The mean accumulated N 2 O emissions averaged 6.8 ± 5.1 kg N 2 O-N ha −1 year −1 over the trial period ( Figure 3). Accumulated N 2 O emissions from fertilized treatments were found to be significantly higher (p < 0.05) than the control. Estimated accumulated N 2 O emissions ranged from 0.7 to 19.3 kg N 2 O-N ha −1 year −1 during the experimental period ( Figure 3). The N80 treatment resulted in the highest accumulated N 2 O emissions (15.5 kg N 2 O-N ha −1 year −1 ) whereas the control plot resulted in the lowest N 2 O emissions (2.45 kg N 2 O-N ha −1 year −1 ). The N20, N40, and N60 treatments resulted in higher emissions than N0 but lower emissions than N80 (3.85, 5.79, and 6.50 kg N 2 O-N ha −1 year −1 , respectively).

Regressions and Emission Factors
The accumulated N2O emissions showed a significant linear and nonlinear relationship with Ninput and Nbalance (p < 0.05) (Figure 4). However, the best fit was obtained using the nonlinear approach and the amount of total Ninput (Eq. 3). Considering the slopes of the linear regressions each kg of additional N will provoke a share of 1.3% and 1.8% for Ninput and Nbalance, respectively. different N fertilization treatments (N0, N20, N40, N60, and N80) over the trial period. Error bars denote standard error of the mean. Different letters indicate significant differences between the treatments (p < 0.05). The N0, N20, N40, N60, and N80 refer to the fertilizer rates used as treatments and were 0, 220, 440, 660, and 880 kg N ha −1 year −1 , respectively.

Regressions and Emission Factors
The accumulated N 2 O emissions showed a significant linear and nonlinear relationship with N input and N balance (p < 0.05) (Figure 4). However, the best fit was obtained using the nonlinear approach and the amount of total N input (Equation (3)). Considering the slopes of the linear regressions each kg of additional N will provoke a share of 1.3% and 1.8% for N input and N balance , respectively.
The EF data differed significantly (p < 0.05) between treatments where the N80 treatment differed from the N20 and N60 but not the N40 treatment (data not shown). On average between the different treatments, the EF data accounted for 0.9% ± 0.1%. Comparing the different approaches with regards to their ability to predict the annual N 2 O emission sufficiently, the nonlinear regression E3 performed the best (Table 4). In comparison, the IPCC approaches showed a good fit but slightly underestimating (IPCC min ) or heavily overestimating (IPCC min+exc ) the measured values. The regressions E1, E2, and E3 showed comparable R 2 and RMSE values but higher NSE compared to IPCC approaches. The EFdata differed significantly (p < 0.05) between treatments where the N80 treatment differed from the N20 and N60 but not the N40 treatment (data not shown). On average between the different treatments, the EFdata accounted for 0.9% ± 0.1%. Comparing the different approaches with regards to their ability to predict the annual N2O emission sufficiently, the nonlinear regression E3 performed the best (Table 4). In comparison, the IPCC approaches showed a good fit but slightly underestimating (IPCCmin) or heavily overestimating (IPCCmin+exc) the measured values. The regressions E1, E2, and E3 showed comparable R 2 and RMSE values but higher NSE compared to IPCC approaches.  IPCCmin+exc). The EFs are calculated based on the N2O measurements (EFdata) and linear (E1-E2) and nonlinear (E3-E4) regression approaches for the different N fertilization treatments (N0, N20, N40, N60, and N80 expressed as kg N2O-N ha −1 year −1 ). The R square (R 2 ), root mean square error (RMSE), and the Nash-Sutcliffe model efficiency coefficient (NSE) are shown as a result of the comparison between measured and predicted annual N2O emissions indicating the best fit of measured vs. predicted values. The N0, N20, N40, N60, and N80 refer to the different mineral fertilizer application rates used and were 0, 220, 440, 660, and 880 kg N ha −1 year −1 , respectively.   year −1 ). The R square (R 2 ), root mean square error (RMSE), and the Nash-Sutcliffe model efficiency coefficient (NSE) are shown as a result of the comparison between measured and predicted annual N 2 O emissions indicating the best fit of measured vs. predicted values. The N0, N20, N40, N60, and N80 refer to the different mineral fertilizer application rates used and were 0, 220, 440, 660, and 880 kg N ha −1 year −1 , respectively.

Nitrogen Yield
The herbage production from pastures in the current study was in the range of 18.7-21.5 t DM ha −1 year −1 . This is comparable to previously reported herbage production values for the southern Cape region of South Africa of between 13.5 and 22.9 t DM ha −1 year −1 [18,27,47,48]. Surprisingly, fertilization with mineral N-fertilizers had only a small effect on the herbage production during the current trial, achieving herbage yields and N yields of 19 t DM and 582 kg N ha −1 year −1 even in the non-N fertilized treatments. Viljoen et al. (2020) [27] demonstrated that the soil at the trial site has the ability to release a significant amount of plant-available N from soil organic matter following soil microbial mediated processes. Additionally, because of animal excreta in intensively grazed systems, the plant nutritional value is relevant. The predicted amount of N applied through excreta accounted for more than 400 kg N ha −1 year −1 , providing the majority of N in the controlled system.
The typical CP value for grazed kikuyu and ryegrass is affected by season and found to be in the range of 10.7-23.7% and 19-28%, respectively [14]. It is well documented that N content of pasture herbage increases when fertilizer is added [27]. This was also observed in the current study where CP and N yield were the highest in the N80 treatment (880 kg N ha −1 year −1 ). However, at the same time the high rate of N fertilization led to a N surplus of >400 kg N ha −1 year −1 , which means a potential economic loss and environmental threats such as increased GHG emissions and N-leaching to groundwater bodies [49]. High N fertilization rates led to CP contents of >25% with the consequence of high feeding costs and also associated with marginal benefits regarding CP in milk yield [50]. Moreover, high N fertilizer rates in pasture are associated with elevated ammonia volatilization through animal excreta as a result of high protein content found in the diet [50]. It is therefore important for common farming practices to consider the current N-status of the soil as well as the quantity of animal excreta N when fertilization is planned. From our experiment on a farm with a long history of dairy management and with high amounts of additional N inputs from concentrates fed in the milking parlor, a mineral fertilizer rate significantly below 200 kg N ha −1 year −1 is still sufficient to provide adequate herbage yields ≈19 t DM ha −1 with a high CP content (≈18%) for dairy, whilst still maintaining the sward N level. It has to be considered that the overall N input from livestock excreta is maybe slightly overestimated in the N0 and N20 treatments, but underestimated in the N60 and N80 treatments due to the experimental constraint that all plots were grazed in the same manner and thus excreta was potentially not evenly distributed at every grazing interval.

Daily N 2 O Fluxes
The use of mineral fertilizer and the effect thereof on N 2 O emissions is well documented in the literature [24,[51][52][53][54], and is also supported by our results. The highest fluxes were observed during the first two weeks after fertilization, which is also in agreement with other observations [55][56][57][58]. Fluxes of N 2 O are very closely related to environmental conditions as well as soil chemical and physical characteristics [59]. Regarding environmental variables, beyond the availability of N in the soil, the key drivers are soil moisture and temperature as they influence N 2 O diffusion to the atmosphere [17] and determine the relative rates of nitrification and denitrification [60,61] by influencing the microbial activity [17,62]. The WFPS in the measured period was in the range of 37-64% and averaged around 51% (Figure 2). Plants therefore did not experience water stress because irrigation was applied. N 2 O production is highest around 60% WFPS and lowest below 30% [63], thus irrigated pastures provided favorable conditions for elevated N 2 O emissions with regards to soil water. Moreover, a higher frequency of N 2 O peaks can be assumed as a result from the sum of naturally and artificial rainfall events. Accordingly, large fluxes of N 2 O coincide with irrigation events as reported by Liu et al. (2006) [57]. In the current study, rainfall and irrigation, together with higher temperatures, were highest during spring and summer (Figure 1). This allowed perfect soil conditions for N-mineralization to take place, accelerating nitrification and denitrification processes and, as a result, increased N 2 O fluxes [64][65][66]. Accordingly, we found a positive correlation for soil surface (0-1 cm) temperature (r = 0.22), WFPS (r = 0.18), and N min (NH 4 -N) (r = 0.11) and NO 3 -N (r = 0.33) on measured N 2 O-fluxes. However, even though we found increasing soil N concentration with increasing fertilization (Table 2 and Figure 2), correlations with N 2 O are flawed by the poor soil sampling intervals. However, this still indicates the importance of soil chemical properties and environmental variables, which needs to be taken into consideration when locally adapted dynamic modeling approaches are developed in the future.
In addition, sward management is also an explanatory reason for the variability of N 2 O-fluxes from grazed pastures. Kikuyu pastures in the southern Cape region of South Africa are usually mulched to ground level and over-sown during autumn [48]. The high N 2 O fluxes observed in April 2019 could be as a result of over-sowing practices and as a consequence, some soil disturbance caused by the furrows created by the seed-drill. A study done by Hao et al. (2001) [21] reported significantly higher N 2 O emissions during autumn under tillage practices, which also implies that soil disturbance is the cause of the elevated flux observed in the current study during autumn. Higher N mineralization occurs as a result of soil disturbance [67] and mulching [68] and leads to a higher amount of plant available N in soil. However, a small leaf surface area of the plants during autumn would limit the uptake of N as N requirements by plants during autumn are low, resulting in more available N in soil which causes higher N 2 O-N losses.

Accumulated N 2 O-losses
The relationship between N input and direct N 2 O emissions from managed agricultural soils is generally assumed to be linear [69,70]. However, others stated that a linear approach may be too conservative [52] and nonlinear approaches are more suitable to predict annual N 2 O emissions [26,49,52,54]. The nonlinear relationships are heavily dependent on the plant N demand leading to different slopes with regards to the crop and its specific N requirements. The plant N demand is also influenced by seasons [27]. As described by Kim et al. (2013) [52], direct N 2 O-N emissions are controlled by competition for available N between plants and microbes when N is provided. Accordingly, the N 2 O-N emissions would increase linearly as N input is increased. However, when N supply is high and exceeds the requirements of plants and microbes, the response would be better explained by an exponential function. As a result, the N surplus from excess N added would lead to lower plant N uptake efficiency [71] leaving residual N as a substrate for additional N 2 O-N emissions [72].
In our study, cows grazed pastures year-round which makes high N returns to pastures unavoidable. High levels of N surplus on pastures under intensive grazing are a result of N returned via urine and feces during grazing combined with the added N fertilizer (Table 2). This N surplus is the principal driver of N 2 O production from these managed soils [73]. Rafique et al. (2011) [49] found on highly productive grassland that an input of 300 kg N ha −1 year −1 (which equaled the N-plant demand) the annual N 2 O-N emissions were less than 5 kg N 2 O-N ha −1 year −1 , but that larger N 2 O-N emissions were observed when N input exceeded 300 kg N ha −1 year −1 , ranging from 5.0 to 12.6 kg N 2 O-N ha −1 year −1 . Accordingly, we found the lowest emissions in the control, which exclusively received only excreted N, resulting in a negative N-balance as no legumes were seeded in the grassland sward.

Predictions of Annual N 2 O Emissions from Pastures
The IPCC suggested a methodology to estimate N 2 O emissions from managed soils according to mineral and organic fertilizers as well as animal excreta. It is considered that 1% of added fertilizer and 2% of animal excreta N is converted to N 2 O-N [1]. Regarding this proposed linear approach, we also found a linear relationship between the annual accumulated N 2 O emissions and N input .
However, in this study we distinguished between the linear approaches by evaluating the effect of mineral fertilizers applied (after emissions from the grazed-only control treatment were deducted) and the N input (from mineral fertilizers and excreta). Measured emissions were compared with a calculated EF data value as well as with the IPCC default values. After deducting the background emissions from grazing, our calculated EF was on average 0.9% and thus in agreement with the IPCC default value of 1% for mineral fertilizers. Previous research on N 2 O emissions focused on one level of N fertilization but as a growing number of field experiments with multiple N fertilization rates have indicated, the relationship is better described as nonlinear [26]. The EFs from the current study follow the same relationship as N 2 O emissions increase nonlinearly when high levels of N input , taking fertilizer and animal excreta into account, exceed plant needs [74]. It was also suggested by Rafique et al. (2011) [49] that the IPCC tier 1 default value should consider a nonlinear relationship between N input and N 2 O emission in grazed grasslands. In accordance with that, our study indicated that the inclusion of the IPCC default value for animal excreta heavily overpredicted the measured N 2 O emissions. It is important to note that grazing animals excrete much of the consumed N as urine and manure back onto pastures, and when deposited unevenly this results in highly localized concentrations or so called "hot spots" of N in soils [75]. In the present study our use of the chamber method might, theoretically, have carried a possible risk that we could not consider all the N 2 O-N emissions emitted by these soils as excreta did not always fall within the chamber bases. However, this is not very likely as it can be assumed that an even distribution of 'hot spots' occurred during 11 grazing cycles and chambers were randomly spread. Recent studies also indicated that the IPCC EF for animal excreta is overpredicted in the grazing period, supporting the validity of our results regarding grazing impact [75][76][77].

Mitigation Strategies
It is evident that excreta are a major contributor to the N 2 O emissions on pasture-based systems [78]. Thus, mitigation can be achieved by taking animal excreta carefully into account when planning N fertilization. In our study, animal excreta-derived N already provided adequate herbage yield and CP contents without adding additional mineral fertilizer. However, when additional mineral fertilizer is used, it should be carefully planned and managed. It is suggested in a recent article by Viljoen et al. (2020) [27] that a variable N fertilization rate, when plant demands for N are high, could be adapted seasonally and be increased or decreased as the season dictates. This would further aid in mitigating N 2 O emissions. Such an approach would be more sustainable and more economically viable in pasture-based dairy systems in the southern Cape of South Africa depending on the practices on the specific farms.
Improved crop varieties and boosting soil health, which would lead to improved nitrogen-use efficiency, is one way of mitigation [15]. Using alternative forage species incorporated into pastures could be another mitigation strategy in reducing N 2 O emissions derived from animal excreta due to biological nitrification inhibition and/or changes in the soil water content [75]. The effective management of forage legumes in a pasture could lead to less N fertilizer used and as a result less N 2 O emissions [79]. Moreover, legume-based forage production shows an efficient N cycling which results in less N surpluses and lower environmental loads of reactive N [80]. However, with an increase in N fertilizer, the legume component of pasture is likely to decrease. A recent review highlighted the potential of plant species to reduce N 2 O emissions by looking at plant-effects on urinary-N 2 O emissions [81].
One of the most effective strategies considered for the reduction of N 2 O emissions under irrigation may be to combine the use of nitrification inhibitors [82] with reduced water volumes, spaced over more applications [83]. Different irrigation methods such as drip and sub-surface drip irrigation can also reduce N 2 O emissions, when compared to surface irrigation methods [84].
The final aim should be to maintain current yields or improve them to sustain a high milk production, all whilst reducing the negative environmental effects associated with dairy production.

Conclusions
Excessive use of mineral fertilizers on intensive rotational grazed pastures in South Africa resulted in high N surpluses (> 400 kg N ha −1 year −1 ) and provided the majority of N in the system. Our results indicated that the relationship between N 2 O-N losses and N input, when plant N demand is considered, can best be described by a nonlinear function rather than a linear function. The suggested EFs of the IPCC tier 1 default value for grazing systems led to an overestimation of N 2 O emissions. A better approach would be to replace EFs of the IPCC tier 1 default value with regional EF values which are dependent on the N balance. This could lead to more accurate N 2 O emissions from managed soils on a regional scale, where other environmental threats such as groundwater pollution and eutrophication are also addressed at the same time.