Nitrous Oxide Emission from Grazing Is Low across a Gradient of Plant Functional Diversity and Soil Conditions

: Nitrous oxide (N 2 O) emissions from pastures can vary signiﬁcantly depending on soil and environmental conditions, nitrogen (N) input, as well as the plant’s ability to take up the N. We tested the hypothesis that legume-based N sources are characterized by signiﬁcantly lower emission factors than mineral N based dairy systems. Therefore, this study monitored N 2 O emissions for a minimum of 100 days and up to two growing seasons across a gradient of plant species diversity. Emissions were measured from both, grazed pastures and a controlled application of urine and dung using the static chamber method. About 90% of the accumulated N 2 O emissions occurred during the ﬁrst 60–75 days. The average accumulated N 2 O emissions were 0.11, 0.87, 0.99, and 0.21 kg ha − 1 for control, dung, urine patches, and grazed pastures, respectively. The N uptake efﬁciency at the excreta patch scale was about 70% for both dung and urine. The highest N 2 O-N emission factor was less than half compared with the IPCC default (0.3 vs. 0.77), suggesting an overestimation of N 2 O-N emissions from organically managed pastures in temperate climates. Plant diversity showed no signiﬁcant effect on N 2 O emission. However, functional groups were signiﬁcant ( p < 0.05). We concluded that legume-containing pasture systems without a fertilizer addition generally appear capable of utilizing nitrogen inputs from excreta patches efﬁciently, resulting in low N 2 O emissions.


Introduction
Agriculture is susceptible to climate change, as global food systems are threatened by an increased risk of intra-and inter-annual yield variability due to the less predictable growing conditions [1]. In addition, it is a contributor to climate change by greenhouse gas (GHG) emissions. Among the GHGs, nitrous oxide (N 2 O) is about 265 times stronger in terms of the global warming potential compared with CO 2 [2], as it destroys the ozone layer in the troposphere [3]. The N 2 O in the atmosphere is increasing linearly at a rate of 0.26% per year [4] and was estimated for the period 2007-2016 to be about 17.0 Tg N year −1 [5]. Agricultural soils are the foremost important contributor to the global anthropogenic N 2 O emissions (about 3.8 Tg N year −1 ) due to the mineral or organic fertilizer application [5]. Similarly, pastures emit N 2 O, as grazing animals return high loads of N (200 to 2000 kg ha −1 year −1 ) in urine and dung patches [6], which exceed the use capacity of grasslands. The excess N results in losses via leaching, volatilization, and N 2 O emissions [7]. organic carbon of 2.2% (Appendix A, Table A1). The study was located in a region of high grazing intensity due to the well-distributed precipitation across the year. The long-term (1981-2019) mean temperature and annual precipitation of the site are 9.7 • C and 859 mm, respectively.
Previously, the site was used for arable cropping in a 5-year crop rotation until 1993 when it was converted to "Bio-land" (prohibiting the use of chemical fertilizers and pesticides) and subsequently managed as grassland (grass-clover mixtures). The grassland was managed in a mixed system (1-2 silage cuts followed by 3-4 grazing cycles by cattle). Two years before the start of the experiment, the site was ploughed to establish the pastures. They were managed organically and without any external nutrient input. Therefore, the legumes were the primary source of nitrogen supply.

Environmental Conditions
Weather parameters, including daily precipitation and soil temperature, were obtained from a meteorological station located within a 70 m radius of the experimental area. The daily maximum and minimum air temperature, daily precipitation, and filled pore space (WFPS) in the soils during the studies are shown in Figure 1. The mean annual temperature at the site was similar for 2018 and 2019 (10.3 • C) but 6% higher than the long-term average (LTA, 1981(LTA, -2019. The mean soil temperature was relatively higher in summer compared with spring or autumn in both 2018 and 2019. In both experimental years, soil temperature following excreta deposition was a little below 10 • C in spring but rose quickly and was between 18-20 • C during summer excreta application, and decreased to about 15 • C during autumn excreta applications ( Figure 1). Precipitation was 25% lower (647.4 vs. 859.2 mm) and 7% higher (922.4 vs. 859.2 mm) for 2018 and 2019, respectively, compared with the LTA. Whereas, the spring of both 2018 and 2019 had 1.6-2 times higher precipitation relative to the LTA, their summers and autumn of 2018 had reduced (20-40%) precipitation compared with the LTA. WFPS ranged between 21.3-71.4% (mean ± SD; 43.5 ± 16.0) in the simulated grazing plots following the excreta deposition in 2018. Accordingly, the urine and dung patches experienced different environmental conditions after each excreta application event ( Figure 1). Therefore, WFPS ranged between 58.2-86.8% (mean ± SD = 77.8 ± 7.7%) across seasons following the excreta application ( Figure 1). Consequently, the mean WFPS was about 1.6 times higher for the irrigated pastures compared with the non-irrigated.

Experiments
Three experiments were conducted to test the hypotheses set out for the study. Experiment one was set up to quantify the N 2 O emission from the three diverse pastures (prwc, prwrc, and prwrch) under grazing stress. The measurements spanned the period of April 2018 to January 2020. The second experiment was simulated grazing trial which also started in 2018. Therefore, cow urine and dung with a known amount of N were applied to the three diverse pastures to establish their response to the N 2 O emission and to calculate the N 2 O emission factors (EF). The third experiment was also a simulated grazing trial established in April 2019 and was meant to establish the N 2 O EF under good soil moisture conditions (which the study site is known for), as 2018 was unusually dry. For this, only the binary mixture (prwc) was used. The pasture was irrigated using the sprinkler method Atmosphere 2021, 12, 223 4 of 26 to simulate rain at least a day before gas sampling each week, depending on the WFPS (i.e., when WFPS fell below 60%), to maintain adequate soil moisture for pasture growth and denitrification. The summary of the experiments, including the specific hypotheses, study design, and measurements are shown in Table 1.

Experiments
Three experiments were conducted to test the hypotheses set out for the study. Experiment one was set up to quantify the N2O emission from the three diverse pastures (prwc, prwrc, and prwrch) under grazing stress. The measurements spanned the period of April 2018 to January 2020. The second experiment was simulated grazing trial which also started in 2018. Therefore, cow urine and dung with a known amount of N were applied to the three diverse pastures to establish their response to the N2O emission and to calculate the N2O emission factors (EF). The third experiment was also a simulated grazing trial established in April 2019 and was meant to establish the N2O EF under good soil moisture conditions (which the study site is known for), as 2018 was unusually dry. For this, only the binary mixture (prwc) was used. The pasture was irrigated using the sprinkler method to simulate rain at least a day before gas sampling each week, depending on the WFPS (i.e., when WFPS fell below 60%), to maintain adequate soil moisture for pasture growth and denitrification. The summary of the experiments, including the specific hypotheses, study design, and measurements are shown in Table 1. Table 1. Summary of the experimental designs and measurements.

Experimental Design
Measurements Experiment 1: Grassland diversity affects N2O emission from pastures under grazing stress Static chambers established on the three diverse pastures (prwc, prwrc, and prwrch); in a randomized complete block design, laid out in three replicate blocks (18 experimental units in all); each replicate plot measured 12 × 9 m; Jersey dairy cows rotationally grazed on a 3-5-week cycle depending on Gas sampling done at least once a week to measure N2O; gas sampling started from April 2018 to January 2020, weather parameters were monitored.  Gas sampling was done at least once a week to measure N 2 O; gas sampling started from April 2018 to January 2020, weather parameters were monitored.

Experiment 2: Botanical and environmental factors determine N 2 O emission from soils
Static chambers arranged in a split-split-split plot design; the diverse grasslands as main plots, excreta treatments as split-split plots, and season as split-split-split plots, laid out in three replicated blocks (27 experimental units per each set-up); urine and dung applied once in spring, summer and autumn; each application was measured for 1 year; duplicate plots for soil sampling.
Gas sampling started immediately after excreta application; was done at least once a week for 1 year, except in December where only three samples were taken; weather parameters, soil mineral N concentration, biomass dry matter and N yields, species composition, water-filled pore space (WFPS) were monitored. Static chambers were arranged in a split-split plot design, using the prwc swards; the season of excreta application as a split-plot and excreta treatment as a split-split plot, all laid out randomly in one block with four replicates (12 experimental units per each set-up); duplicate plots were used for soil sampling.
Gas sampling was done at least once a week for 100 days; weather parameters, mineral N concentration, biomass dry matter and WFPS were monitored; sampling duration was 100 days. We collected dung and urine from dairy cows kept under confinement systems (fed conserved fodder and concentrates) to ensure an almost constant N content of dung and urine throughout the year. The urine and dung specimens were collected immediately after morning milking (~05:00 h). Urine samples were taken midstream after manual stimulation of the vulva, while faeces samples were collected as the animal defecated. Faeces and urine were homogenized and stored immediately under refrigeration (<4 • C) and were applied within 7 days. The urine and dung, when removed from the refrigerator, were allowed to attain ambient temperature before application. Homogenized and representative quantities of fresh urine and dung were applied once per measurement period (1 year or 100 days) in spring, summer, and autumn to avoid introducing new effects due to differences in excreta volume or mass [35]. Urine was applied at a constant rate of 8.8 L m −2 from a height of 40 cm and was expected to cover the entire collar by percolation [36], whereas each dung patch of 2.2 kg covered an approximate area of 0.041 m 2 . The timing of excreta application, N concentrations, and the total amount of N applied per ha are shown in Table 2.

N 2 O Sampling
We sampled gas fluxes using the closed static chamber method (Hutchinson and Mosier, 1981). This consisted of polyvinyl chloride (PVC) collar (h = 15 cm, d = 60 cm) installed to a depth of 10 cm into the soil on each replicate plot. The chamber was made of PVC, 30 and 40 cm high for grazed pastures and simulated grazing, respectively, with an internal diameter of 60 cm. In each sampling session, the chamber was deployed onto the plastic collar between 10:00 and 12:00 h [9], with the chamber-collar sealed with a rubber belt. One air sample was taken with a 30-mL polypropylene syringe from the chamber-top at 0, 20, 40, and 60 min after closure. The air sample in the syringe was transferred to 12 mL pre-evacuated glass vials (Labco, High Wycombe, UK), and stored for a maximum of 2 weeks before analyses.

Soil and Biomass Sampling
Each micro plot for gas measurement had a duplication adjacently placed to evaluate WFPS, ammonium (NH 4 ), and nitrate (NO 3 ) concentrations in the top 25 cm of soil. At least six soil samples were taken in this duplicated plot within the first 2 months following the excreta application and at least once in the subsequent months, using a core sampler. The soil samples were stored immediately at −17 • C until further processing. The biomass was harvested from the rings and micro-plots and was in synchrony with the grazing cycle in experiment 1. Harvesting was done by cutting all forage within each subplot at a 5 cm height from the soil surface. The total fresh matter (FM) weight of the biomass from each plot was recorded immediately after harvesting. After determining the fresh weight, two sub-samples (80-100 g) were taken, one for dry matter determination and nutrient analyses, while the second sub-sample was sorted into three functional groups (grass, legumes, and non-legume herbs).

Chemical Analyses
Urine and dung samples were analyzed in duplicates by dry combustion using the C/N analyzer (Vario Max CN, Elementar Analysensysteme, Hanau, Germany) to determine total carbon and nitrogen concentrations, using helium as a purge and carrier gas. Two sub-samples of dung were dried at 60 • C to determine the dry mass content. Biomass sub-samples were oven-dried to a constant weight at 55 • C to determine their dry-matter. The dried samples were then milled through a 1-mm sieve (Cyclotech mill, Foss analytical, Hilleröd, Denmark) and later used to determine the N concentration by near-infrared reflectance spectroscopy (NIRS) using a NIR-System 5000 monochromator (Foss Analytical, Hillerød, Denmark). Soil samples were extracted with 0.01 M CaCl 2 ; NO 3 and NH 4 concentrations were determined photometrically, using a dual-channel continuous flow analyzer (Skalar Analytical Instrument, Breda, the Netherlands). We used 50-80 g of each soil sample to determine the gravimetric moisture by oven-drying the fresh samples at 105 • C to a constant weight. This index was used to correct for soil moisture and to calculate WFPS. The N 2 O concentration in the sampled gas was analyzed by gas chromatography (model 7890a, Agilent Technology Inc., Santa Clara, CA, USA) equipped with a 63 Ni electron-capture detector. Helium was used as the carrier gas and argon-methane as the make-up gas. Three certified gas standards (300, 620, 1510 ppb) were measured regularly to determine the analytical precision (standard deviation of 10 repeated measurements <3 ppb). The gas chromatograph procedure used a detector temperature of 320 • C, a column temperature of 40 • C, and an injector temperature of 200 • C.

Calculations
We estimated biologically fixed N (BFN) using a model (Equation (1)) by Høgh-Jensen et al. [37]: where DM legume is the accumulated amount of legume shoot dry matter above the standard defoliation height; N% is the concentration of N in the dry matter of the legume (kg kg −1 ); P fix is the fixed N 2 as the proportion of total N in the shoot dry matter of the legume; P root+stubble is the fixed N 2 in the root and stubble as the proportion of totally fixed shoot N at the end of the growing period; P transsoil is the belowground transfer of fixed legume N 2 located in the grass mixtures as the proportion of total fixed shoot N at the end of the growing period; P immobile is the fixed N 2 immobilized in an organic soil pool at the end of the growing period as a proportion of fixed shoot N; P transanimal is the above-ground transfer (by grazing animals) of fixed legume N 2 located in the grass in mixtures as the proportion of total fixed shoot N at the end of the growing period. To parameterize the model for N%, P fix , P root+stubble , P transoil , and P immobile , we used 0.0375, 0.845, 0.25, 0.15, and 0.315, respectively for swards with white and red clover, and 0.04, 0.80, 0.25, 0.25, and 0.38, for white clover only swards [37], P transanimal was set to zero. We considered the N applied, BFN, and annual deposition by precipitation (12.5 kg N ha −1 ; [38]) as the N input and N yield in harvested biomass as the N output. Patch residual N (PRN) was calculated as the difference between the N input and output [39]. This parameter indicated the available N liable to denitrification and nitrification losses, leaching, ammonia emissions, and immobilization. Volumetric water content (VWC) was calculated as the product of the gravimetric water and soil bulk density (BD) ( Table A1); and WFPS was calculated (Equation (1)) after considering VWC, the BD, and a particle density (PD) of 2.65 g cm −3 : We calculated NO 3 and NH 4 intensities in the 0-25 cm soil depth for 100 days sampling duration by submitting the daily NO 3 and NH 4 data to linear interpolation; the NO 3 or NH 4 intensity is an integrated measure of the accumulation of mineral N in the soil over time [27]. N 2 O-N fluxes were calculated for each treatment and replicate by linear regression between measured N 2 O-N concentrations and time [39]. Hourly (Equation (3)), daily, cumulative N 2 O emissions, as well as emission factors, were estimated according to methods described by Krol et al. [30]. The hourly N 2 O-N flux (µg m −2 h −1 ) was estimated by taking account of the slope (∆c/∆t) of the line that describes the concentration increase of N 2 O inside the chamber during the 60 min deployment; the temperature, pressure (assumed as 1 atm), and internal volume of the chamber (V); and the area of the micro plot delimited by the metal base (A); the molar weight (M) of N in N 2 O; and the ideal gas constant (R): Daily N 2 O-N emission (g ha −1 d −1 ) was calculated by multiplying the hourly fluxes by 24 h. The cumulative emission of N 2 O-N in each season (kg ha −1 ) was estimated by integrating the fluxes over each season's monitoring period (area under the curve). The N 2 O-N emission factor (EF) was calculated using the IPCC [31] recommended method (Equation (4)), which only considered the amount of N applied. However, as all swards contained legumes, we also estimated N 2 O-N EF using (Equation (5)) based on the total N input (modelled BFN + N applied):

Statistics
We performed all the statistical procedures and produced all the graphs using R [40]. For Experiment 1, the effect of grassland diversity on N 2 O emission from grazed pastures was tested using a generalized linear model with the diverse grasslands as a fixed effect and the block as a random factor. For Experiments 2 and 3, N 2 O emission, soil, and biomass variables were analyzed using a mixed-effect model with grassland, the season of excreta application and excreta type or control as fixed factors and blocks and replicates as random factors, except for 2019 data where the season of excreta application and excreta type were considered as fixed factors and replicates as random factors. The treatment factors, as described, were used to analyze all the variables in a full factorial arrangement. In each case, the adequacy of the model was assessed by examining the appropriateness of residual plots. Variables were log-transformed where necessary before the analysis to ensure a normal distribution of residuals. We separated the effective means at p < 0.05 with Tukey's post hoc tests, using the "lsmeans" function of the "multcomp" package [41]. Regression tests (linear and quadratic) were performed to establish relationships between N 2 O emission/EFs and the measured/estimated variables. The best-fitted models were selected based on their Akaike information criterion (AIC) with a small sample second-order bias correction [42].

Plant Yields, Botanical Composition, and Soil Residual N
All fixed factors (grassland diversity, the season of excreta application, and excreta type) had significant effects on the measured biomass yields and species compositions (p < 0.05), except for the herbs and unsown species which were not significantly affected by the differences in seasons (Table 3). Whereas, grass proportions in the swards decreased with the increasing diversity and legume proportions. In addition, BFN were highest, intermediate, and lowest for prwrc, prwrch, and prwc, respectively. The dry matter yield and herbage N were lower for the binary swards compared with the tertiary or the multispecies swards. The urine and dung application significantly affected most of the measured parameters but not the dry matter yield or herb proportions (Table 3).  Grass proportions and herbage N were highest, intermediate, and lowest for urine, dung, and the control plots, respectively. However, legume proportions and the resulting BFN were lowest, intermediate, and highest for urine, dung, and the control plots, respectively. The seasonal differences or time of excreta application significantly affected legume proportions, dry matter yield, herbage N and BFN (Table 3). Whereas, the legume proportion only differed between spring and summer excreta treatments, dry matter yield and herbage N differed significantly between the three seasonal treatments, being highest, intermediate, and lowest for spring, summer, and autumn, respectively. BFN was higher for spring treatments compared with summer or autumn.
Grassland diversity, the excreta application season and excreta type significantly affected (p < 0.05) N input, N uptake, and PRN (Table 3). N input was highest, intermediate, and lowest for prwrc, prwrch, and prwc, respectively. However, N uptake was lower for the binary swards compared with the tertiary or the multi-species swards, while PRN was only higher for prwrc than prwc or prwrch. Similar to the legume proportion and resulting BFN, N input was highest, intermediate, and lowest for urine, dung, and the control plots, respectively. N uptake was only higher for urine-treated plots compared with the control plots, while the resulting PRN was lower for the controls compared with urine or dung treated swards. N input and PRN differed significantly (p < 0.05) between the three seasonal treatments, being highest, intermediate, and lowest for spring, summer, and autumn, respectively, but N uptake was higher in autumn compared with summer or spring.

N 2 O Emissions from Grazed Pastures with a Gradient in Diversity
Daily N 2 O-N fluxes from the grazed grasslands ( Figure 2) ranged from 0.04-3.24 mg m −2 d −1 (mean ± SD = 0.31 ± 0.33). The 8-month (May-December) mean(SD) flux was slightly higher in 2019 (0.35 ± 0.46) compared with 2018 (0.29 ± 0.23). Grassland diversity had no apparent impact on the daily N 2 O-N fluxes, except for 2019 where the highest peak was observed from the prwrc swards ( Figure 2). The peaks usually followed precipitation events and appeared to occur a few days after a grazing activity. Despite no grazing events, some appreciable emissions (up to about 1.0 mg N m −2 d −1 ) also occurred during the winter months. The mean (±se) annual N 2 O-N emissions were 0.11 ± 0.01 kg N ha −1 and 1.0 ± 0.08 kg N ha −1 for the ungrazed and grazed pastures, respectively, with the differences between the pastures being insignificant (Figure 3; p > 0.05). The 8-month (May-December) average was slightly higher for 2018 compared with 2019 (0.71 ± 0.06 vs. 0.56 ± 0.11 kg N ha −1 ) with the differences between pastures and years or their interaction being insignificant (Figure 3; p > 0.05).

N 2 O Emission from Dung and Urine Patches (Simulated Grazing) Temporal N Fluxes
The excreta deposition increased the soil mineral N concentration in soils of both nonirrigated and irrigated pastures across seasons, being higher for urine compared to dung ( Figure 4A,B). Each excreta application resulted within the next days in a large increase in soil NH 4 and subsequently, in soil NO 3 , but with varying magnitudes depending on the season and excreta type ( Figure 4). As expected, the incidence of NH 4 peaks occurred earlier than NO 3 peaks as urine is first transformed into NH 4 and then into NO 3 . NO 3 flux from spring-applied urine patches was higher for irrigated pastures compared with the non-irrigated ones. For the non-irrigated pastures, the summer application resulted in the highest and most extended peaks of NH 4 concentrations when compared with spring and autumn ( Figure 4A(a)). NO 3 concentrations returned to background levels within the same timeframe as NH 4 in spring but appeared to persist over a more extended period in summer and autumn ( Figure 4A(b)). Daily N 2 O-N fluxes ranged between −0.2-24.2 mg m −2 d −1 with the highest peaks in spring from urine patches ( Figure 4A(c)). The fluxes were back to background levels between 1-3 months, depending on the season ( Figure 4A(c)). Subsequent fluxes that occurred following the excreta application in 2018, particularly in summer and autumn, coincided with precipitation events after the initial peak fluxes (Figure 1).
The highest NH 4 peak from the irrigated pastures occurred in spring but with a shorter lifespan compared with the other seasons ( Figure 4B(a)). NO 3 concentrations in summer levelled to the background within the same time frame as NH 4 ( Figure 4B(a,b)). Whereas, NO 3 concentrations in soils were affected by the excreta type and excreta deposition timing, NH 4 concentrations were only affected by the excreta type (Table 3). Both NO 3 and NH 4 concentrations were highest, intermediate, and lowest for urine, dung, and control plots. NO 3 concentrations were three-four times lower when excreta was applied in spring or summer than the autumn applications (p < 0.05). Daily N 2 O-N fluxes ranged from −5-293 g ha −1 d −1 with the highest peak in spring from urine patches ( Figure 4B(c)). Generally, peak N 2 O-N fluxes' incidences were earlier from urine patches than from dung patches following the excreta application in all seasons. The fluxes were back to background levels between 1-3 months, depending on the season ( Figure 4B    background levels within the same timeframe as NH4 in spring but appeared to persist over a more extended period in summer and autumn ( Figure 4A(b)). Daily N2O-N fluxes ranged between −0.2-24.2 mg m −2 d −1 with the highest peaks in spring from urine patches ( Figure 4A(c)). The fluxes were back to background levels between 1-3 months, depending on the season ( Figure 4A(c)). Subsequent fluxes that occurred following the excreta application in 2018, particularly in summer and autumn, coincided with precipitation events after the initial peak fluxes (Figure 1). The highest NH4 peak from the irrigated pastures occurred in spring but with a shorter lifespan compared with the other seasons ( Figure 4B(a)). NO3 concentrations in summer levelled to the background within the same time frame as NH4 ( Figure 4B(a,b)).

Cumulative N 2 O Emissions
The annual N 2 O-N emission from Experiment 2 (non-irrigated grasslands) following 2018 excreta applications ranged from 0.071-0.130, 0.288-1.69, and 0.515-1.70 g ha −1 d −1 for control, dung, and urine patches, respectively. All the N 2 O emissions were similar across grasslands, but N 2 O emissions were affected by the type of excreta or time of application (p < 0.05; Table 3; Appendix A, Table A2). Accordingly, cumulative N 2 O-N emissions were 8-10 times higher for urine or dung treatments than the control and were highest in spring, and lowest in autumn, with 420 and 230 g N ha −1 fewer emissions in summer and autumn, respectively. However, grassland and season or treatment interaction effects were not significant (p > 0.05; Appendix A, Table A2). In addition, N 2 O-N emissions per unit N uptake was significantly affected (p < 0.05) by the excreta type and season of excreta application (Table 3; Appendix A, Table A2). Accordingly, this parameter was similar for urine and dung but higher than the control, and was highest, intermediate, and lowest in spring, summer, and autumn, respectively.
Cumulative N 2 O-N emissions from Experiment 3 (irrigated pastures, sampled for 100-days) ranged from 0.02-0.09, 0.256-1.118, and 0.958-2.615 kg ha −1 from control, dung, and urine treatments (Table 3) with mean ± SD of 0.81 ± 0.91 kg ha −1 . The cumulative N 2 O-N emissions from non-irrigated similar swards and a similar sampling duration were 0.46 ± 0.55 kg ha −1 (mean ± SD). The cumulative N 2 O emissions were significantly (p < 0.01) affected by the treatment type and treatment-season interaction effects (Appendix A, Table A2). Therefore, whereas cumulative N 2 O emissions were higher for urine patches than dung in spring, the emissions were similar in summer and autumn but higher than the control plots (Appendix A, Figure A1). Whereas, cumulative urine N 2 O emissions were higher in spring relative to summer or autumn, the dung N 2 O emission only differed between spring and summer (Appendix A, Figure A1).

Drivers of N 2 O Emission
We fitted linear and quadratic models around the measured cumulative N 2 O emission from excreta patches using the botanical or soil N variables as predictors. For the nonirrigated pastures, an exponential model involving N input as the predictor variable explained 49% of the variations in annual N 2 O emission (Figure 5a). However, predicting annual N 2 O-N emissions using PRN (Figure 5b) yielded a relatively higher coefficient of determination (R 2 = 0.77). The cumulative N 2 O-N emissions from the irrigated pastures were best described by the mineral N intensity across seasons, with a higher coefficient of determination for the NH 4 -N intensity than the NO 3 -N intensity (R 2 = 0.84 vs. 0.65; Figure 5c,d).

N2O Emission Factors
N2O-N EFs were calculated based on 100 days of sampling, higher than the IPPC recommended 30 days threshold. Calculating N2O-N EF based on (Equation (5)) reduced the N2O-N EFs by 16% across grasslands, excreta type, and seasons of its application (Table 2). N2O-N EFs for the non-irrigated and irrigated pastures ranged from 0.01-0.35% (mean ± SD = 0.16 ± 0.11%) and 0.04-0.45% (0.25 ± 0.14%), respectively. For the nonirrigated swards, N2O-N EF was significant for the pasture and N source interaction (p < 0.05; Appendix A, Table A2, Figure 6a), with prwc swards having lower EF for dung compared with prwrc. In the case of the irrigated swards, there was a treatment and season interaction effect (p < 0.01; Appendix A, Table A2). Therefore, whereas the treatment effect was only evident in spring, N2O-N EF was only higher for dung in summer than spring  (5)) reduced the N 2 O-N EFs by 16% across grasslands, excreta type, and seasons of its application (Table 2). N 2 O-N EFs for the non-irrigated and irrigated pastures ranged from 0.01-0.35% (mean ± SD = 0.16 ± 0.11%) and 0.04-0.45% (0.25 ± 0.14%), respectively. For the nonirrigated swards, N 2 O-N EF was significant for the pasture and N source interaction (p < 0.05; Appendix A, Table A2, Figure 6a), with prwc swards having lower EF for dung compared with prwrc. In the case of the irrigated swards, there was a treatment and season interaction effect (p < 0.01; Appendix A, Table A2). Therefore, whereas the treatment effect was only evident in spring, N 2 O-N EF was only higher for dung in summer than spring (Figure 6b).

Drivers of N2O-N EF
We tested the effect of the measured/estimated parameters on the calculated N2O-N EF using data from Experiments 2 and 3 (the results are shown in Appendix A, Table A3). We found that the soil temperature, precipitation, excreta N and C content, N loading, NO3-N intensity, legume and grass proportions in the swards, sward N uptake, and PRN significantly (p < 0.05) affected N2O EFs for the non-irrigated swards. For the irrigated swards, excreta C and N, N applied, and soil mineral N intensity, as well as the dry matter yield, significantly affected N2O-N EF (p < 0.05; Appendix A, Table A3). The increasing grass proportions resulted in lower N2O-N EF (Figure 7a), while legume proportions >40% resulted in increased N2O-N EF (Figure 7b).
As expected, the N uptake (Figure 7c) was negatively correlated with N2O-N EF, while PRN associated positively with a higher N2O-N EF (p < 0.05; Table 3; Figure 7d).

Drivers of N 2 O-N EF
We tested the effect of the measured/estimated parameters on the calculated N 2 O-N EF using data from Experiments 2 and 3 (the results are shown in Appendix A, Table A3). We found that the soil temperature, precipitation, excreta N and C content, N loading, NO 3 -N intensity, legume and grass proportions in the swards, sward N uptake, and PRN significantly (p < 0.05) affected N 2 O EFs for the non-irrigated swards. For the irrigated swards, excreta C and N, N applied, and soil mineral N intensity, as well as the dry matter yield, significantly affected N 2 O-N EF (p < 0.05; Appendix A, Table A3). The increasing grass proportions resulted in lower N 2 O-N EF (Figure 7a), while legume proportions >40% resulted in increased N 2 O-N EF (Figure 7b).  Table A4; Figure 7e,f). It appears N2O-N EF only increased when WFPS increased above 60% (Figure 7f) and when soil temperatures increased above 15 °C (Figure 7f).

Plant Yields, Botanical Composition, and Soil Residual N
We did not consider yields and forage uptake from the grazed plots (Experiment 1). However, a previous report by Lorenz [43]showed high forage production from these fields (9.86-11.64 Mg ha −1 year −1 ). As expected, dry matter forage production and N uptake were higher in the diverse prwrch pastures than the binary prwc pastures in Experiment 2. This difference might be partly explained by the higher BFN inputs in the prwrch pastures due to their higher legume (N) yield ( Table 3). The similar dry matter production of prwrch and prwrc despite the higher N input in the mixture without herbs might be due to the higher ability of multispecies swards to efficiently utilize resources for growth. Previous studies reported higher below and above-ground biomass yield by more diverse compared with less diverse swards [20]. The particularly low dry matter yield of prwc pastures can additionally be at least partly explained by the inability of white clovercontaining swards to withstand drought [44], the summer of 2018 was especially dry but had an extremely productive spring and an extremely productive long autumn, with excellent growing conditions for white clover. White clover suffered only 6 weeks of As expected, the N uptake (Figure 7c) was negatively correlated with N 2 O-N EF, while PRN associated positively with a higher N 2 O-N EF (p < 0.05; Table 3; Figure 7d). Whereas, soil temperature was significantly and linearly associated with non-irrigated N 2 O-N EFs, the environmental factors did not associate significantly with irrigated N 2 O-N EFs. However, using pooled data from both experiments for conducted regression approaches showed that WFPS and soil temperature were most related to N 2 O-N EF within the first 15 and 60 days following excreta application, respectively, and were quadratic functions (Appendix A, Table A4; Figure 7e,f). It appears N 2 O-N EF only increased when WFPS increased above 60% (Figure 7f) and when soil temperatures increased above 15 • C (Figure 7f).

Plant Yields, Botanical Composition, and Soil Residual N
We did not consider yields and forage uptake from the grazed plots (Experiment 1). However, a previous report by Lorenz [43] showed high forage production from these fields (9.86-11.64 Mg ha −1 year −1 ). As expected, dry matter forage production and N uptake were higher in the diverse prwrch pastures than the binary prwc pastures in Experiment 2. This difference might be partly explained by the higher BFN inputs in the prwrch pastures due to their higher legume (N) yield ( Table 3). The similar dry matter production of prwrch and prwrc despite the higher N input in the mixture without herbs might be due to the higher ability of multispecies swards to efficiently utilize resources for growth. Previous studies reported higher below and above-ground biomass yield by more diverse compared with less diverse swards [20]. The particularly low dry matter yield of prwc pastures can additionally be at least partly explained by the inability of white clover-containing swards to withstand drought [44], the summer of 2018 was especially dry but had an extremely productive spring and an extremely productive long autumn, with excellent growing conditions for white clover. White clover suffered only 6 weeks of drought. The range of dry matter yield observed in this study agrees with the 9-12 Mg ha −1 reported for similar grasslands under similar climatic and soil conditions, depending on the age or soil N status [43,45].
As the N concentration of both urine and dung could not be controlled, N loading rates varied across seasons and treatments with means (±SD) of 493 ± 43 and 397 ± 47 kg ha −1 for urine and dung patches, respectively, but were similar across grasslands. Still, these differences appear negligible compared to the range of N loading rates of 200-2000 kg N ha −1 year −1 from urine and dung patches reported for dairy cows [6]. The excreta application increased the grass proportion in the swards but reduced the legume yield due to the availability of N [46]. With this, the partner grass outcompetes the legume easily [47], especially the shallow-rooted white clover. This same argument could be offered to explain the low legume yield by prwc (Table 3).
BFN in grass/legume swards worldwide ranges 13-682 kg N ha −1 year −1 , depending on the legume persistence, yield, associated grass(es), or measurement method. Therefore, the BFN range (51-286 kg N ha −1 year −1 ) observed in this study could be considered normal. Whereas, the N applied and N input correlated positively (but weakly; r = 0.27-0.33), while the dry matter yield, legume proportion, and BFN correlated more strongly (r = 0.63-0.76). This observation confirmed the critical role of BFN in the dry matter production and supported the assertion that legumes deliver a more synchronized N to meet plant demand compared with the N supply from external sources and stimulate the over-yielding of biomass [19,48]. The relatively high dry matter yield from the control plots compared with dung and urine treatments (Table 3) further corroborates this assertion. Indeed, Ledgard [49] contended that pastoral systems could rely solely on BFN for moderate to high dry matter production with moderate N losses.

N 2 O Emissions and Emission Factors
N 2 O emissions from grazed pastures are associated with substantial heterogeneity due to the random nature of emission hotspots from urine and dung depositions, making the eddy covariance method preferred to the chamber method [29] for non-simulated grazing experimental designs. However, the chamber method can capture the magnitude and temporal patterns from emission hotspots (i.e., urine and dung patches), and it makes the comparison of different treatments in parallel possible. Moreover, the accuracy of measurements can be improved if the interval of measurements is increased [10]. Nevertheless, the high coefficient of variation (about 23%) associated with the current measurements might have undermined the detection of significant differences. Due to this, the additional experiment with controlled urine and dung depositions was designed, which allowed creating more uniform N loadings and tracking the emissions from exactly one urine or dung patch.
The range of daily N 2 O-N fluxes observed from our grazed pastures was similar to the 1.1-2.9 mg N m −2 d −1 reported for grazed Irish pastures on sandy soils [50] and 0.03-2.5 mg N 2 O m −2 d −1 reported in New Zealand [51]. However, the range of annual N 2 O-N emission observed for the grazed pastures (mean ± SD = 1.0 ± 0.2 kg N ha −1 ) appears lower than the values reported by these earlier authors due to the infrequent incidence of high peaks (Figure 2). The pastures in these earlier studies received additional chemical N fertilizers and were conducted in locations with slightly higher temperatures compared with our site. These factors might explain the lower fluxes in our case since farm management contributes to the effect of site characteristics and climate variability on the N 2 O emission [52]. Nevertheless, the result agrees with the 0.04-21.21 kg N 2 O-N ha −1 year −1 reported for European grasslands (Jena Experiment), depending on the site and management [52].
The range of cumulative N 2 O-N fluxes (0.05-1.59 kg N ha −1 ) observed in this study from experiment 2 was slightly lower than the 0.08-3.2 kg N ha −1 reported for similar simulated grazing experiments [7,22], probably since these previous authors used crop swards. The N 2 O EFs observed for urine and dung in the current study might be related to the soil, climate, and management practices. Previous studies conducted in similar climates [25,[27][28][29][30] reported low EFs for urine and dung in line with our findings. De Klein et al. [53] and Marsden et al. [26] reported low (0.03-0.3%) N 2 O-N EF for animal excreta deposited on free-draining soils. The soils at the site are free-draining sandy loam (Appendix A, Table A1), known to have increased infiltration and percolation of urine, leading to increased N losses via leaching [27] and reduced N 2 O emission. Recently, low N 2 O-N EFs (0.12-0.41%), depending on the age of sward and fertilization (slurry) rate, was observed for grass-legume grasslands at the study site [45]. As the current study measured N 2 O emissions for, at least, 100 days, the N 2 O-N EFs could be considered site-specific. However, the summer of 2018 was unusually dry, though spring was typical, and autumn very productive.
When N 2 O-N EF was estimated based on the total N input (i.e., N 2 O-N EF 1 = BFN + N applied), to account for the contribution from legumes in the swards, the mean EF reduced by 16% but was not significantly (p > 0.05) affected by the sward diversity nor excreta type (Table 3). However, this method of estimating EF contrasts with the IPCC 2006 guidelines, which requires that background emissions are considered. Nevertheless, this method provides an opportunity to illustrate the impact of legumes on N 2 O emission and might guide the maximal N input for an appreciable level of N 2 O emission from grass-legume swards [45]. The N 2 O EFs from our study were relatively low compared with the revised IPCC default EF of 0.77%.
Holtan-Hartwig et al. [54] observed that some soils generally have the intrinsic potential for low N 2 O emission irrespective of environmental conditions, and this has been attributed to differences in denitrifying communities [55]. Moreover, previous and current soil management practices might affect microbial dynamics in soils. The initial N status of soils is likely to have reverberating effects on the microbial soil community's functions, including N transformations, abundance, and activity of nitrifying and denitrifying microorganisms [56], and by extension the production of N 2 O. The microbial community's composition is an essential factor that determines the N 2 O:N 2 ratio emitted from temperate soils [57]. Therefore, the soils were not fertilized and must have had a low residual N, as shown by the control plots' negative PRN (Table 3).

The Effect of Grassland Diversity
We expected lower N 2 O emissions from the most diverse pasture (prwrch) under both grazing and simulated grazing experiments, based on the assumption that diverse pastures increase plant productivity through increased niche complementarity, resource uptake efficiency, and N transfer from legumes to grasses/herbs, and particularly, as these also contained plants reported to reduce N 2 O emission via nitrification inhibitors [7,22]. A decreasing trend in N 2 O emission as grassland diversity increased was observed from grazed plots in 2019 ( Figure 3) and when excreta N 2 O emitted was related to the N uptake (Table 3). However, the annual N 2 O, N 2 O EF or N 2 O emission intensity were statistically similar for the grasslands. Moreover, the timing of peak N 2 O production appeared similar across grasslands, seasons, and treatments. This observation contradicts reports by Hoeft et al. [58] who observed alterations in the timing of peak N 2 O production after urine addition due to species differences. Although the N input differed among the grasslands, obviously due to the differences in legume proportions, N 2 O emission was similar across the grasslands. The mineral N concentration was similar for the grassland soils ( Table 3), suggesting that grasslands might not be too different functionally to influence N cycling. Therefore, we reject the hypothesis that a more diverse grassland would necessarily reduce the N 2 O emission.
The effect of sward diversity on N 2 O emission is not conclusive. In previous studies, Hoeft et al. [58] reported no difference in the N 2 O emission between grass pastures and more diverse pastures. In contrast, Niklaus et al. [59] observed reduced N 2 O emission from more diverse pastures (species richness from 1-16) but not in the presence of fertilizer or large proportions of legumes. In this study, herbs constituted about 18% of the most diverse swards and forage plantain proportion was only around 5% (results not shown). This is likely too low for any effects on N 2 O emissions. Additionally, the soil-environment interaction might affect plantain BNI effectiveness [7,23]. According to Experiment 2, prwc yields less dry matter and less nitrogen (Table 3), and had possibly a higher uptake by the animals, which would have resulted in higher droppings of urine and dung. Therefore, similar N 2 O emissions from higher N-yields in the more complex grasslands could be considered positive.

Effect of N Source
The effect of excreta type on the N 2 O emission was not consistent as the interaction effects were evident for both non-irrigated and irrigated pastures. However, PRN and all the N 2 O emission indices were generally similar for both dung and urine in non-irrigated pastures, but with irrigation N 2 O emissions/EFs were 1-2 times higher for urine patches compared with dung. Luo et al. [7] also reported a higher N 2 O emission from urine patches relative to dung, and this could be ascribed to the slow mineralization of dung N. Whereas, urine contains high amounts of available N, dung is rich in C and organic N, and undergoes slower mineralization rates and higher fixation of N by adsorption to soil particles, depending on soil properties and environmental conditions [60]. However, the N mineralization from dung patches and subsequent emission of N 2 O might be a function of biological, environmental conditions, and time, as over time immobilized N becomes available for mineralization resulting in a later peak, as occurred in experiment 2. It is not clear if all N from dung patches mineralized during the period as dung N were bound in very slowly decaying organic bindings, or more complex molecules might need more time than urine, in which most N is in mineral form. The relatively dry soil conditions in 2018 might have limited N 2 O emission from the urine patches. The drought must have slowed down mineralisation even further in dung patches as they remained visibly dry in the collars.

Effect of Timing of Excreta Deposition
The peak N 2 O emissions were visually associated with a decline in NH 4 and an increase in NO 3 concentration in the soils. The highest N 2 O-N peak was obtained in spring, followed by summer, and then autumn in both years, similar to previous reports [7,30,61]. The seasonal effect on PRN, N 2 O-N emission, and N 2 O-N EF might be due to the prevailing soil moisture and temperature interaction effects at or following excreta deposition [61]. These environmental differences influence the N uptake and residual N in the soils, and hence N-leaching to the groundwater, which is a function of timing of excreta deposition, soil inorganic N status, and climatic conditions [62]. Although the leaching load was not measured in this study, Biernat et al. [62] reported higher leaching loads when soil mineral N was high in late summer, autumn or early winter. This observation might partly explain the low residual N or N 2 O emission associated with summer and autumn excreta applications relative to spring.

Effect of Irrigation
The irrigated pastures in Experiment 3 had relatively high WFPS compared with Experiment 2 (Figure 1). It is worth mentioning that irrigation was adopted to provide conducive soil conditions for denitrification losses. Such high soil moisture condition could promote plant growth and enhance plant N uptake but could also promote N leaching to the underground water pool. Similar to reports by Smit et al. [63], the irrigation increased (almost doubled) both cumulative N 2 O-N emitted (0.46 vs. 0.81) and N 2 O-N EF (0.25 vs. 0.16%) compared with non-irrigation emissions. The dry conditions in 2018 must have reduced N mineralization, particularly in dung patches as they were visibly dry, and would have limited denitrification. With the irrigated swards, both nitrification and denitrification were significant processes by which N 2 O fluxes were produced; nitrification seemed to be more dominant as shown by the more robust relationship between N 2 O and NH 4 compared with NO 3 (R 2 = 0.84 vs. 0.65; Figure 5c,d) [27].

Botanical Characteristics
The measured/estimated botanical variables were significantly associated with N 2 O-N EFs in this study. The impact of plants on N 2 O emissions and underlying mechanisms are not always conclusive [21] as several factors might play confounding roles. However, grasslands' ability to influence N 2 O emission might depend on factors including N uptake from soils, legume proportion, and the production of soluble sugars [21,59]. Regression analyses showed a significant effect of grass proportion, legume proportion, and N uptake on N 2 O-N EFs (Appendix A, Table A3, Figure 7a-c). The higher grass proportion and N uptake appeared to depress the N 2 O emission factor, but the higher legume proportion associated positively with N 2 O-N EF. Accordingly, Bracken et al. [64] reported a significant functional group identity effect on cumulative N 2 O emissions, with emissions increasing with the increasing white clover proportion and decreasing with the increasing perennial ryegrass proportion for both irrigated and non-irrigated pastures. Similarly, Niklaus et al. [59] reported a significant legume effect on N 2 O emission. On the contrary, Barneze et al. [65] found no legume effect on N 2 O emission, total mineral N or mineral N in fertilized or unfertilized soils.
Under monoculture, legumes appeared to emit more N 2 O than other species [59]. However, this might not be necessarily additive in mixed pastures [66], but increases in the proportion of legumes might influence C and N cycling processes [65]. This legume effect on N 2 O emission might be related to the increased N supply (Table 3), leading to nitrification and denitrification losses. Therefore, the increased NO 3 availability due to the release of fixed N from legume roots via decomposition and legumes' inefficiency in acquiring mineral N [65] may explain the high association of legumes with the N 2 O emission. The inclusion of legumes in mixed-swards at proportions of 30-50% is described as most effective, providing multiple outputs including reduced greenhouse gas emission, higher productivity and increased protein self-sufficiency, and reduced N leaching [48,67]. In this study, legume proportions ranging from 30-50% resulted in N 2 O-N EFs ranging from 0.17-0.25% (Figure 7b). At the same time, legume proportions up to 30% have been widely promoted in New Zealand and Switzerland as a norm [19].

Soil Mineral N Availability and Plant N Uptake
In our study, soil mineral N intensity, PRN, and N uptake were significantly associated with N 2 O-N EF. The N input was 116-289 kg ha −1 more than the 300 kg N ha −1 required by highly productive grasslands to meet plant demand [50]. At the same time, the N uptake (320-540 kg ha −1 ) was lower than the N input, resulting in patch residual N ranging −87-206 (mean = 103 kg ha −1 ), being positive for excreta patches (Table 3) and negative for the control plots. This estimate might not be precise as lateral flows from urine and dung patches outside the rings were possible [36]. Also, N losses via nitrification, denitrification, and leaching were not considered. Moreover, these pastures can sequester about 1 Mg C per year, equivalent to~100 kg N (assuming a CN ratio of 10), captured in soil organic matter. Nonetheless, the residual soil N on the entire pasture must be far lower considering that excreta deposition by grazing livestock affects about 60% of the paddock. In any case, N 2 O production is primarily driven by the amount of residual N in the soil [57], as corroborated by the high association between N 2 O emission and PRN (R 2 = 0.77; p < 0.001; Figure 5b), which is a function of N input and uptake of the swards. The higher residual N will potentially result in higher N 2 O emission as N becomes available for nitrifies and denitrifies to act on. It appears values up to 250 kg ha −1 might still keep N 2 O-N EF from cow excreta patches below 0.2% (Figure 5b) at this site, but the high surplus N might be leached to pollute underground water pool and/or become a secondary source of N 2 O emission.
PRN was about 67% higher for prwrc than prwc or prwrc, leading to a relatively low N uptake for prwrc. The similar N uptake for prwrc and prwrch might be due to the high grass proportion (Table 3) as grasses have a higher potential to utilize soil N for growth, hence, reducing soil N and the nitrification/denitrification potential of the soils [22]. The range of N uptake by the grasslands in this study appears to be on the upper side of the range of values (10-450 kg N ha −1 year −1 ) reported by Bessler et al. [68]. It appears grass-legume swards are associated with the high N uptake, as corroborated in this study and by Reinsch et al. [45]. This has been attributed to species complementarity evident by the development of high root length and root density by grass-legume species [69], leading to low PRN.
The high N uptake (>70%) in this study might partly explain the low N 2 O-N EFs observed in this study, as the higher N uptake is associated with low N 2 O-N EF (Figure 7c). The higher N uptake reduced N 2 O EF, and the higher N input might stimulate the higher N uptake, as observed in this study (results not shown; r = 0.44, p < 0.05), similar to reports by Lü et al. [70]. This observation supports the general concept that the nutrient uptake by roots increases in response to the increasing nutrient supply [71]. Besides, this highlights the critical role of plant N uptake and its potential usage in mediating the effects of N addition from grazing animals. It must, however, be noted that the plant N uptake/productivity depends also on the availability of other nutrients and has a physiological maximum. Although N 2 O emissions were similar across the grasslands, the most diverse grassland had a relatively high N uptake compared with the binary-sward grassland, suggesting a higher potential for reduced N 2 O emission as corroborated by the 9% lower N 2 O emission intensity by prwrch compared with prwc (Table 3).

Soil Moisture and Temperature
As shown by regression analyses, soil moisture and temperature were significant predictors of N 2 O-N EF for both non-irrigated and irrigated pastures but appeared to be most vital within the first two weeks and two months, respectively, following excreta application (Appendix A, Table A4). Soil moisture and ambient temperature, are known to influence N 2 O emission significantly by modulating the microbial activity and transport of gases in soil [10,11]. Higher N 2 O-N EFs associated positively with higher WFPS. However, it appears that increases in N 2 O-N EF only occurred when WFPS was above 60%, respectively (Figure 7e). Soil moisture affects nitrification and denitrification processes [11], and for that matter, ammonia-oxidizing species [12], reduce their population below 60% of field capacity. The soil water threshold value differs according to the soil type [57], and the critical WFPS for denitrification in many soils is equivalent to the field capacity or above [72]. At higher WFPS where denitrification becomes the primary process, the N 2 produced might exceed the N 2 O production, leading to reduced N 2 O emissions at high WFPS, especially when temperatures are high. Davidson [73] explained that both nitrification and denitrification processes co-occur within soils, but the low soil water content and coarse soil texture favour nitrification, while the high soil water content and fine soil texture with a high organic matter content promotes denitrification [74]. The seemingly high occurrence of nitrification might be due to the sandy nature of soils at the site. This might partly explain the relatively low N 2 O emission observed in this study.
The higher N 2 O-N EF associated positively with the higher soil temperature, but it appears that increases in N 2 O-N EF only occurred when the soil temperatures went above 15 • C (Figure 7f). Maag and Vinther [75] reported maximum nitrification rates at 20 • C and small responses to soil temperature changes at lower temperatures. Although denitrification depends on the soil type, with sandy loam soils responding significantly to both increases in soil temperature and soil moisture, temperatures above 5 • C are often considered to be required for effective denitrification rates [76]. However, we observed an appreciable N 2 O emission from grazed fields in winter and early spring (Figure 2) similar to reports by Flessa et al. [77]. In winter, the observed N 2 O emissions, despite the absence of grazing activity, might be due to substrates accumulations from the proceeding grazing season. Such winter emissions have been attributed to specific freeze/thaw phenomena [78]. However, Holtan-Hartwig et al. [54] explained that decreasing temperatures increase the N 2 O/N 2 product ratio under anaerobic soil conditions as a result of differences in the temperature response of N 2 O diffusion versus biological process rates.

Conclusions and Implications of the Study
In this study, the grasslands' species richness did not significantly affect the N 2 O emission, although there was a tendency for low N 2 O emission from the prwrch swards when the N 2 O emission was expressed as a fraction of the N uptake. Differences in the functional composition (i.e., grass and forage legumes) of the swards were significantly related to the N 2 O emission. We observed a high N uptake efficiency by the grasslands, partly explaining the low N 2 O emission and EFs at the site. In addition, soil temperature and WFPS, botanical factors, including legume or grass proportions in the swards, as well as N uptake were dominant factors that drove the N 2 O emission from dung and urine patches. The high N uptakes observed for grasslands show their high potential to optimize the N use. An elevated WFPS almost doubled the N 2 O-N EFs, however, the EF values were below the default IPCC [32] for wet temperate climates, suggesting a further disaggregation of the default value to cater for organic dairy systems in wet temperate regions.
The excreta CN ratio, N applied, total N input, PRN, grass and legume proportions, as well as N uptake were factors that drove the N 2 O emission from dung and urine patches. The high N uptake (>70%), legume proportions (16-46%), as well as the soil type might explain the relatively low N 2 O EFs observed in this study. Results suggest that simple regression approaches using grass or legume proportions and N uptake can be used to predict the N 2 O emission more precisely over a wide range of environments. We further concluded that an efficient proportion of legumes without additional mineral fertilizer use accelerates an efficient N-cycling and is more important than a complex species diversity to reduce undesired N 2 O emissions from pastures. However, the year 2018 was unusually dry and results might be or were probably affected by this. Hence, the experiments or at least the simulated grazing studies should be continued to see the difference to "normal" years. Acknowledgments: We wish to thank Rita Kopp, Petra Voß, and Katrin Helmich of the Grassland and Forage Science/Organic Agriculture Department at the CAU for their immense technical assistance during the analyses of soils, plants, and gas samples, respectively. We are indebted to Thomas Ehmsen and Keanu Heuck of the Organic Research Farm Lindhof for their unflinching support during the project implementation and data collection.

Conflicts of Interest:
The authors declare that they have no competing interests.       Level of significance, ns = p > 0.05, * < 0.05, ** < 0.01, *** < 0.001; WFPS: Water-filled pore space. Table A4. Effect of sampling days on the relationships (R 2 ) between N 2 O-N emission and environmental factors (soil temperature and water-filled pore space).