First on-line isotopic characterization of N 2 O emitted from intensively managed grassland

Introduction Conclusions References


Introduction
Atmospheric nitrous oxide (N 2 O) mole fraction is increasing since pre-industrial times predominately due to increased agricultural activity (Davidson, 2009;Mosier et al., 1998).Owing to the approximately 300 times higher global warming potential (GWP) compared to CO 2 , this greenhouse gas (GHG) currently accounts for 6 % of total anthropogenic radiative forcing (Myhre et al., 2013).Recent estimates showed that N 2 O is in addition the single most important ozone-depleting substance (Ravishankara et al., 2009).Because at least 60 % of total anthropogenic N 2 O emissions is attributed to food production (Syakila and Kroeze, 2011), growing human population and meat consumption per capita as well as biofuel production will accelerate the rate of increase in atmospheric N 2 O concentration.Hence, the development of adequate mitigation strategies is pertinent and requires a better understanding of the processes driving N 2 O fluxes.To date, nitrification, nitrifier denitrification and denitrification are considered to constitute the dominant N 2 O producing processes, especially in agricultural soils (Wrage et al., 2001).Other N 2 O source-processes such as abiotic N 2 O production, co-denitrification and heterotrophic nitrification have also been observed; a concise overview of observed processes is given elsewhere (Butterbach-Bahl et al., 2013).This complexity inherent in the N cycle and associated transformation processes is a major challenge in developing mitigation strategies, as attribution of N 2 O production to the respective processes is required to tailor target-oriented actions (Baggs, 2008).Approaches for apportioning of N 2 O emissions to nitrification, denitrification, and N 2 O reduction to N 2 (source partitioning) have mostly relied on acetylene (C 2 H 2 ) inhibition and isotope labeling (Groffman et al., 2006), but denitrification rates are underestimated by the C 2 H 2 method (Butterbach-Bahl et al., 2013;Groffman et al., 2006;Watts and Seitzinger, 2000).Isotope labeling approaches are vulnerable to Figures incomplete diffusion of the tracer and to stimulation of process rates by the addition of the labeled substrates themselves (Groffman et al., 2006).Changes in natural abundance of 15 N and 18 O in N 2 O have been explored to investigate N 2 O production processes, but the determined δ 15 N and δ 18 O depend on both fractionation factors and isotopic composition of precursors, which in turn exhibit strong variations (Baggs, 2008;Bedard-Haughn et al., 2003;Heil et al., 2014;Toyoda et al., 2011).N 2 O is a linear molecule and four main isotopic species can be discerned: 14 N 14 N 16 O, 14 N 15 N 16 O, 15 N 14 N 16 O and 14 N 14 N 18 O.The isotopic species 14 N 14 N 16 O, 14 N 14 N 18 O and 14 N 15 N 16 O (or 15 N 14 N 16 O) are isotopologues, while 14 N 15 N 16 O and 15 N 14 N 16 O are isotopomers and will be termed 15 N α -N 2 O and 15 N β -N 2 O (Toyoda and Yoshida, 1999).The umbrella term isotopocule is used for both isotopomers and isotopologues.The intra-molecular distribution of 15 N in N 2 O ("site preference"; SP = δ 15 N α − δ 15 N β ) has been reported to be independent of the substrate's isotopic composition, as SP remained constant even though δ 15 N and δ 18 O values of both N 2 O and substrates changed markedly during experiments with pure cultures (Heil et al., 2014;Sutka et al., 2003Sutka et al., , 2006Sutka et al., , 2008;;Toyoda et al., 2005).Therefore, SP can be considered as a tracer conserving the source process information (Ostrom and Ostrom, 2011).The SP of different processes has been characterized in a number of pure-culture, mixed culture (Ostrom et al., 2007;Sutka et al., 2003Sutka et al., , 2006;;Toyoda et al., 2005;Wunderlin et al., 2012Wunderlin et al., , 2013)), and soil-incubation studies (Köster et al., 2011(Köster et al., , 2013a;;Lewicka-Szczebak et al., 2014;Well et al., 2006Well et al., , 2008) ) with a compilation of data inToyoda et al. (2011).A recent review on source partitioning and SP (Decock and Six, 2013)  and Toyoda, 2000).Owing to the temporal and spatial variability of isotopomer ratios, it is indispensable to derive flux-weighted average values from different sources (such as ecosystems) for later use in budget analysis using box models (Kim and Craig, 1993;Perez et al., 2001;Yoshida and Toyoda, 2000).Intramolecular distrubution of 15 N in N 2 O can be measured by mass spectrometry, but it requires discrete flask sampling with subsequent laboratory analysis.Hence, this approach is limited in temporal and spatial resolution.Additionally it is indirect, as information on the site-specific isotopic composition is derived from the analysis of the NO + fragment and N 2 O + molecular ion.Recently, a quantum cascade laser absorption spectrometer (QCLAS) capable of selective analysis of the three most abundant N 2 O isotopocules has been presented (Waechter et al., 2008) and its potential for in situ measurements in conjunction with an automated pre-concentration unit has been shown (Mohn et al., 2010(Mohn et al., , 2012)).Here we present the results obtained from a, to our knowledge worldwide first, campaign in which the isotopic composition of N 2 O (SP, δ 15 N, δ 18 O) in the atmospheric surface layer was determined on-line by using an optimized state-of-the-art laser spectrometer.With the combination of N 2 O isotopic analysis by QCLAS, accompanying eddy-covariance based N 2 O flux measurements as well as monitoring of environmental conditions and inorganic nitrogen concentrations, our specific objectives for this study were: (i) to demonstrate the capability of QCLAS systems for high precision isotopic analysis of (soil emitted) N 2 O in ambient air, (ii) to investigate management and weather effects on isotopic composition and source processes, (iii) to test the capability of the N 2 O isotopic composition for source partitioning, and (iv) to characterize the flux-averaged isotopic composition of N 2 O emitted from an intensively managed grassland.Introduction

Conclusions References
Tables Figures

Back Close
Full 2 Material and methods

Study site
The agricultural research station Chamau (CHA) is located in Central Switzerland at an elevation of 400 m a.s.l.The experiment was conducted on an intensively managed grassland belonging to CHA which is primarily used for fodder production and occasional winter grazing by sheep (Zeeman et al., 2010).The soil type is a cambisol with a pH of 5-5.5 (Roth, 2006).Mean annual temperature and annual precipitation are 9.1 • C and 1151 mm, respectively (Merbold et al., 2014).Management practices aim at fodder production and consist of mowing followed by slurry application, with up to six mowing/slurry applications per year and occasional grazing of sheep and cattle in October and November.During the campaign in summer 2013, three management cycles were carried out.Harvest dates were 6 June, 11 July and 21 August and slurry was applied within 10 days after each mowing event.Nitrogen input was calculated from the applied amount of slurry brought to the field and the N concentration determined (Labor für Boden-und Umweltanalytik, Eric Schweizer AG, Thun, Switzerland) in a sample drawn from the supply to the trailing hose applicator.The applied N amounted to 30, 40 and 43.3 kg N ha −1 for the first, second and third application, respectively.The grassland is re-established via ploughing and resowing approximately every 10 years.The last re-establishment event took place in 2012.

Instrumental setup for analysis of N 2 O isotopocule ratios
The four most abundant N 2 O isotopic species were quantified using a modified QCLAS (Aerodyne Research Inc., Billerica MA, USA) equipped with a continuous wave quantum cascade laser (cw-QCL) with spectral emission at 2203 cm −1 , an astigmatic Herriott multi-pass absorption cell (204 m path length, AMAC-200), and reference path with a short (5 cm) N 2 O-filled cell to lock the laser emission frequency (Tuzson et al., 2013).During the campaign, the QCLAS was operated in an air-conditioned trailer Introduction

Conclusions References
Tables Figures

Back Close
Full located 60 m west of the eddy-covariance (EC) tower.This trailer position contributes < 20 % to the main flux and is at the far side of prevailing wind direction (Zeeman et al., 2010).The sample air inlet was installed next to the inlet of the EC tower (2 m height).Sample air was drawn through a PTFE tube (4 mm ID) by a membrane pump (PM 25032-022, KNF Neuberger, Switzerland).Upstream of the pump, the sample air was pre-dried with a permeation drier (MD-050-72S-1, PermaPure Inc., USA).Following the pump, the pressure was maintained at 4 bar overpressure using a pressure relieve valve.Humidity, as well as CO 2 , were quantitatively removed from the gas flow by applying a chemical trap filled with Ascarite (7 g, 10-35 mesh, Fluka, Switzerland) bracketed by Mg(ClO 4 ) 2 (2 × 1.5 g, Fluka, Switzerland).Finally, the sample gas was passed through a sintered metal filter (SS-6F-MM-2, Swagelok, USA) and directed to a preconcentration unit described in detail previously (Mohn et al., 2010(Mohn et al., , 2012)).For an increase of N 2 O mixing ratios from ambient level to around 50 ppm N 2 O, approx.8 L of ambient air were preconcentrated.Afterwards, the preconcentrated N 2 O was introduced into the evacuated multi-pass cell of the QCLAS.Isotopic fractionation during preconcentration (0.31 ± 0.10, 0.34 ± 0.16 and 0.29 ± 0.07 ‰ for δ 15 N α , δ 15 N β and δ 18 O, respectively) was quantified by preconcentration of N 2 O with a known isotopic composition and subsequently corrected.Compatibility of N 2 O isotopomer analysis by QCLAS with isotope ratio mass spectrometry (IRMS) laboratories was recently demonstrated in an inter-laboratory comparison campaign (Mohn et al., 2014).

Measurement and calibration strategy
To ensure high accuracy and repeatability of the analytical system, a measurement and calibration strategy similar to the one presented by Mohn et al. (2012)  Full Schweiz AG) resulted in pressurized gas mixtures with 90 ppm N 2 O (parts per million, 10 −6 mol of trace gas per mole of dry air).Both standards were calibrated against primary standards which were previously measured by the Tokyo Institute of Technology (TIT, Toyoda and Yoshida) to anchor δ values to the international isotopic standard scales.The first standard (S1, Table 1) was used as an anchor point to the international δ scale and used as input data for data analysis algorithms (see data processing).Therefore, the N 2 O isotopic composition of S1 was targeted to closely resemble background air.As the N 2 O isotopic composition of surface layer air is mainly a mixture of soil-derived and background composition, the second standard (S2, Table 1) used for span correction was depleted in δ 15 N α , δ 15 N β and δ 18 O compared to background air in accordance with the expected terrestrial source signatures.
The measurement protocol started with the injection of S1, dynamically diluted to 50 ppm, the mole fraction of ambient N 2 O after preconcentration.After flushing the absorption cell with synthetic air, S2 was injected, also diluted to 50 ppm.For determination of the slight concentration dependence already reported (Mohn et al., 2012), S1 was injected again but at a higher mole fraction of 67 ppm (later referred to as S1 h ).This mole fraction represents the mole fraction expected after preconcentration of high concentration surface layer air.Subsequently, S1 was injected again, diluted to 50 ppm, before the cell was filled with preconcentrated ambient N 2 O (A).This subroutine (S1 + A) of injection of S1 and preconcentrated ambient N 2 O took 35 min and was repeated three times.For an independent determination of repeatability, the fourth sample was preconcentrated compressed air (target gas).During the campaign, two compressed air cylinders (C1 and C2, referred to as target gas) were used.Isotopic composition and N 2 O mixing ratio of both cylinders were determined in the laboratory prior to campaign start (Table 1).N 2 O mole fractions and isotopic composition analysed in the laboratory and at the field site agreed within their analytical uncertainty.Following target gas analysis, S1 and S1 h were analyzed again.Another set of three subroutines S1 + A completed one run.One complete cycle of 6 ambient air samples and one compressed air sample took 340 min, leading to approx.25 ambient air samples Introduction

Conclusions References
Tables Figures

Back Close
Full being analysed during 24 h.N 2 O mole fractions were determined according to Mohn et al. (2012).

Data processing
Data processing is based on individual mixing ratios of the four main N 2 O isotopic species and spectrometer characteristics as recorded by the instruments's software (TDLWintel, Aerodyne Research Inc., Billerica, MA, USA).In the first step, variations in the isotope ratios induced by drifts in the instrument working parameters during the field operation were corrected.A linear additive model explaining the deviation of isotope ratios R α , R β and R 18 O for repeated measurements of standard S1 from their mean value by absorption cell temperature (T1), laser temperature (T2), line position (LP) and pressure (p) was calibrated based on S1 injections.For isotope ratios of S1, S1 h , S2, sample air and compressed air, these systematic deviations were corrected based on the respective values of T1, T2, LP and p.In a second step, concentration dependence of isotope ratios, determined using the measurements of S1 and S1 h , was addressed with corrections (0.013, 0.028 and 0.004 ‰ ppb −1 for δ 15 N α , δ 15 N β and δ 18 O) being in the same range as described earlier (Mohn et al., 2012).Subsequently, remaining drifts were corrected based on analysis of S1.Finally, isotope ratios were converted to δ values using a 2-point calibration derived from corrected values of S1 and S2.

Determination of soil-emitted N 2 O isotopic composition
Isotopic composition of the source process "soil N 2 O emission" was derived using the Keeling plot approach (Keeling, 1958) air to the atmospheric background.At night, the surface layer becomes more stable and the N 2 O mole fraction increases, shifting isotopic composition towards its source composition.As a consequence, Keeling plots were based on noon-to-noon periods.This approach is discussed in Sect.4.6.

N 2 O Flux measurement
At CHA, greenhouse gas mole fractions, including N 2 O, are measured continuously since 2012 by means of the eddy covariance (EC) method (Baldocchi and Meyers, 1998).The system consists of a three-dimensional sonic anemometer to measure wind speed and direction (2 m height, Solent R3, Gill Instruments, Lymington, UK) and a QCLAS (mini-QCLAS, Aerodyne Research Inc., Billerica, MA, USA) to determine N 2 O mole fractions at a temporal resolution of 10 Hz.Both data streams are merged near-real time within a data acquisition system (MOXA embedded Linux computer; Moxa, Brea, CA, USA) via an RS-232 serial data link (Eugster and Plüss, 2010).The setup has been described in detail previously (Merbold et al., 2014).Post-processing of N 2 O fluxes included screening for obvious out-of-range values (±100 nmol m −2 s −1 ).
N 2 O fluxes were further aggregated to noon-to-noon daily averages to smoothen the large variability in the 30 min flux averages.Daily averages were calculated for days where more than 30 half-hour values were available, with this filter excluding three days from analysis.

Soil inorganic N and environmental conditions
Ammonium (NH + 4 ) and nitrate (NO − 3 ) concentrations were determined from soil (0-20 cm depth) sampled at 10 positions along a transect within the footprint of the EC measurements following the predominant wind direction.Samples were taken weekly throughout the campaign or daily during mowing and slurry application events.Per sample, ∼ 15 g of fresh soil were added to specimen vessels containing 50 mL 1 M Introduction

Conclusions References
Tables Figures

Back Close
Full KCl.After 1 h on a shaker, the supernatant was filtered (Whatman no.42 ashless filter paper, 150 mm diameter) and analysed colorimetrically for NH + 4 and NO − 3 .Soil temperatures and volumetric soil moisture contents at 10 cm depth were measured at the same 10 locations along the transect (5TM-sensors, Decagon Devices Ltd., Pullman, USA).Data were stored as 10 min averages on a data logger (EM50, Decagon Devices Ltd., Pullman, USA).The volumetric water content was converted to water filled pore space (wfps) using a bulk density of 1.09 g cm −3 .Precipitation was measured with a tipping bucket rain gauge (Type 10116, Toss GmbH, Potsdam, Germany) and stored as 10 min averages on a data logger (CR10X-2M, Campbell Scientific Inc., Logan, USA).

Long term precision for target gas analysis
System performance for N 2 O mole fractions and isotopic composition was determined based on repeated analysis of compressed air from target gas tanks (C1, C2).There was no significant drift in the δ values and N 2 O mole fractions, indicating stability of the applied measurement technique.Repeatability, calculated as the SD (σ) of 331 target gas measurements, amounted to 0.20, 0.12, 0.10, 0.12 and 0.22 ‰ for δ 15 N α , δ 15 N β , δ 18 O, δ 15 N bulk and SP, respectively (Fig. 1).SD for the N 2 O mole fraction of the target gas was 0.25 ppb.

Auxiliary measurements
Half hourly N 2 O fluxes were averaged from noon-to-noon (f N 2 O ), and ranged from −1 mole fractions (∆N 2 O) as analysed by QCLAS and could not be attributed to slurry application events alone (Fig. 3).Among the correlations of f N 2 O and auxiliary variables, only the one with nitrate concentration (r 2 = 0.18) was significant (p < 0.01).Soil water content (wfps) was modulated by precipitation and two clear states could be identified.
During the "wet" part of the campaign lasting until 7 July, average wfps was with 62±4 % significantly (t test, p < 0.001) higher than the average of 37 ± 4 % calculated for the remainder of the campaign (referred to as the "dry" part).Soil temperature did not show such a clear two-phase pattern, however temperatures during the first, "wet" part were with 16.7 ± 4 • C significantly (p < 0.001) lower than during the "dry" phase with Background NH + 4 and NO − 3 concentrations were smaller than 3 µg g −1 soil and clearly responded to mowing and slurry application in the second and third management events.The NO soil for these two management events.In contrast, dissolved organic carbon concentrations (DOC) did not respond to management events, but were higher during the "dry" phase of the campaign (p < 0.001).Introduction

Conclusions References
Tables Figures

Back Close
Full The uncertainty of the determined source isotopic composition was estimated based on the standard error of the Keeling plot intercept and depends on the degree to which soil air accumulated in the surface layer (∆N 2 O, Fig. 4).For instance, the intercept (source) standard error ranged from 0.3 to 82 ‰ for SP.To apply the Keeling plot approach only to situations in which soil air accumulated in the surface layer, only source isotopic compositions for overnight increases in N 2 O mole fractions of more than 12 ppb were considered in this study.This filter lead to a maximum and average (µ) standard error of 6.8 (µ = 2.2), 4.5 (µ = 1.4) and 2.2 (µ = 1) ‰ for SP, δ 15 N bulk and δ 18 O isotopic source signatures, respectively.
During the field campaign, soil temperature and NO − 3 were significant.However, the adjusted r 2 for all regressions was below 0.4; in addition, multiple explanatory variables such as NH + 4 and NO − 3 or wfps and temperature (Fig. 5) did not increase the explained variance above this value.

Event-based data aggregation
As already described in Sect.3.3, there was a "wet" phase (n = 27) in the beginning of the campaign, which lasted about one month and a "dry" phase lasting about two months (n = 38).Therefore, the dataset was split in two corresponding parts with averages of 7.4 ± 3.6 vs. 11.1 ± 4.2 ‰ for SP, −19 ± 3.8 vs. −12.5 ± 5.9 ‰ for δ 15 N bulk and 28.7 ± 2.2 vs. 29.7 ± 3.4 ‰ for δ 18 O in the wet vs. the dry phase, respectively.
Averages of SP and δ 15 N bulk were significantly different (p < 0.001) but δ 18 O averages were not.Based on this simple classification, the dry phase contains rewetting events.

BGD Introduction Conclusions References
Tables Figures

Back Close
Full A rewetting event was defined as a two day period starting at the day for which wfps increased.Exclusion of these rewetting events during the dry phase increased average δ-values (n = 30) as well as decreased SDs for SP, δ 15 N bulk and δ 18 O to 12.5±3.4,−10.8±4.5 and 30.7±2.8 ‰.Moreover the difference in δ 18 O was significant (p < 0.001).
In addition to the dry/wet classification, we also defined three subsets representing the N 2 O emission associated with management events of mowing followed by fertilization ("Mana I"-"Mana III"), one subset representing a rewetting event between Mana II and III ("Rewetting") and one subset representing background ("BG", all remaining measurements).There were two distinct rewetting events between management events II and III, but N 2 O isotopic composition is only available for the first one (29-31 July 2014).Isotopic compositions of soil-emitted N 2 O were assigned to subsets of management or rewetting if the associated flux or nutrient concentration was elevated.This classification scheme led to 3-7 measurements for management and rewetting events (Fig. 3, underlaid in transparent blue) while 47 measurements were assigned to class BG.Boxplots for SP, δ 15 N bulk , δ 18 O, and wfps (Fig. 6) showed characteristic δ-values and wfps for management and rewetting, but not for subset BG.Measurements assigned to BG covered practically the whole range of values observed across all the other classes.Therefore, SDs for class BG were one order of magnitude larger than for the four other classes.
Statistical analysis is confounded by low and unequal sample size so that we compared exclusively the subsets management and rewetting using multiple nonparametric Wilcoxon tests after having checked homogeneity of variances using Bartlett test.For all investigated δ-values, only differences between groups Mana II and Mana III were significant.Introduction

Conclusions References
Tables Figures

Back Close
Full To avoid bias to the flux-weighted average of emitted N 2 O by either one of the above mentioned possible reasons, the weighted averages were calculated for positive flux events only.Flux weighted averages were 6.9 ± 4.3, −17.4 ± 6.2 and 27.4 ± 3.6 ‰ for SP, δ 15 N bulk and δ 18 O respectively (n = 48).

Analytical performance
To our knowledge, only two pilot studies exist demonstrating the potential of QCLAS based analytical techniques for on-line and high-precision analysis of N 2 O mole fractions and isotopic composition in surface layer air.Full analytical performance was determined, in accordance with the presented study, based on repeated analysis of compressed air samples.Thereby, the analytical precision reached in the presented study, was distinctly higher for δ 15 N β and δ 18 O and similar for δ 15 N α compared to these two previous studies, even though the measurements were done under field-conditions and over a much longer, three months, period.This confirms the high level of precision associated with QCLAS based determination of N 2 O isotopic composition.Standard errors for Keeling plot intercepts (Fig. 4) confirm that this precision is sufficient to resolve the variability of atmospheric N 2 O sampled close to the ground.As our instrument was located directly at the field site and measurements were conducted over a period of more than three months, our study indicates that this level of repeatability can be achieved both at long time scales and in the field.

N 2 O isotopic composition in the atmospheric surface layer (2 m height)
In our study, δ-values of single preconcentrated air samples were between atmospheric background and 14.3 ‰ (SP) and −4.7 ‰ (δ 15 N bulk ).Mohn et al. (2012) reported similar values between atmospheric background and 12 ‰ (SP) and −4 ‰ (δ 15 N bulk ).Therefore the variation observed in both studies is much higher compared to the measurements by Harris et al. (2014)  combustion processes indicate that some of these sources might be less depleted or even enriched in 15 N-N 2 O compared to ambient N 2 O (Harris et al., 2015;Ogawa and Yoshida, 2005).

Isotopic composition of soil-emitted N 2 O
SP of soil-emitted N 2 O observed in our study (1-17 ‰) is within the ranges expected for a mixture of the two process groups N 2 O N and N 2 O D , and does not necessarily indicate significant contribution of N 2 O reduction, an effect which is discussed further below.Isotopic composition of soil-emitted N 2 O has been predominately determined in laboratory incubation studies (Köster et al., 2013a, b;Perez et al., 2006;Well and Flessa, 2009b;Well et al., 2006Well et al., , 2008)).Additionally, results from field experiments using static chambers (Opdyke et al., 2009;Ostrom et al., 2010;Toyoda et al., 2011;Yamulki et al., 2001) and N 2 O accumulation below a snowpack have been published (Mohn et al., 2013).Based on pure culture studies SP values from 19.7 to 40 and −8.7 to 8.5 ‰, were observed for N 2 O N and N 2 O D , respectively (Decock and Six, 2013).In field experiments SP was found to range between −1 and 32 ‰ (Opdyke et al., 2009), −3 and 18 ‰ (Yamulki et al., 2001), −14 and 90 ‰ (Toyoda et al., 2011) and 0 and 13 ‰ (Ostrom et al., 2010).The very high SP values detected by Toyoda et al. (2011) may have resulted from extensive N 2 O reduction to N 2 , a process increasing SP, δ 15 N bulk and δ 18 O (Ostrom et al., 2007).For δ 15 N bulk and δ 18 O, a much wider variation as compared to SP is expected, because these variables depend both on fractionation factors, which vary among different microbial communities and depend on reaction conditions, as well as on the isotopic composition of the substrate (Baggs, 2008).Under field conditions, δ 15 N bulk was reported to range between −17 and 9 ‰ (Opdyke et al., 2009), −27 and 1 ‰ (Yamulki et al., 2001), −44 and 34 ‰ (Toyoda et al., 2011)  For identification of processes determining N 2 O isotopic composition, isotopocule maps were suggested in which site preference is plotted vs. the difference in substrate and product isotopic composition (Koba et al., 2009).Determination of isotopic composition in the substrates is time consuming and additionally confounded in our study by the large and varying footprint area.Therefore, we present a modified isotope map of SP vs. δ 15 N bulk (Fig. 7) instead of ∆δ 15 N, the δ 15 N differences between substrate and product (i. in the isotope values from denitrification may have been caused by variations in the extent to which N 2 O was reduced to N 2 .Additionally, δ 18 O was found to be positively correlated with δ 15 N bulk , which enforces the interpretation that varying shares of N 2 O reduction occurred because it acts on both N and O isotopic composition (Koehler et al., 2012).It is noteworthy that based on such modified isotope maps, systematic changes in δ 15 N bulk induced by systematic changes in N isotopic composition of one of the precursors NH + 4 or NO − 3 could be misinterpreted as reduction events.The ratios of fractionation factors for δ 18 O and δ 15 N bulk (r o-n ) and SP and δ 18 O (r sp-o ) during N 2 O reduction were suggested for estimation of the share of N 2 O reduction to N 2 since these ratios were found to be 2.5 and 0.2 to 0.5, respectively in laboratory incubation experiments (Jinuntuya-Nortman et al., 2008;Ostrom et al., 2007;Well and Flessa, 2009a).We calculated these ratios for a subset of data for which all δvalues (SP, δ 15 N bulk and δ 18 O) increased for two consecutive days, indicating that N 2 O reduction occurred.Such events were observed on 8 occasions.If source processes (N 2 O D , N 2 O N ) contributed constantly over two consecutive measuring days, changes in the isotopic composition of emitted N 2 O were solely attributed to changes in the fraction of N 2 O reduction.Under such conditions one would expect that the ratio of the changes in δ 18 O and δ 15 N bulk (r o-n ) is around 2.5 and that the ratio of the changes in SP and δ 18 O (r sp-o ) is between 0.2 and 0.5.The mean (median) ratios for r o-n and r sp-o for these selected events were 0.69 (0.44) and 2.1 (1.16), respectively.While the high values of r sp-o indicate that for instance changing physical conditions such as soil moisture may play a role in field measurements, the deviation of r o-n from the value of 2.5 could either indicate that the fractionation factor for 18 O might be smaller than the one for 15 N or that there is no correlation of fractionation factors in natural environments.This is in line with recent findings showing that apparent isotope effects associated with N 2 O reduction are sensitive to experimental conditions which influenced diffusive isotope effects (Lewicka-Szczebak et al., 2014).The same study also showed that fractionation factors during N 2 O reduction for 15 N and 18 O were variable (from −11 to +12 ‰ and from −18 to +4 ‰, respectively), and not predictable for field conditions yet.

BGD Introduction Conclusions References
Tables Figures

Back Close
Full Therefore, to date, the amount of N 2 O reduction prior to emission cannot be inferred with sufficient robustness from field measurements alone, without the knowledge of isotopic composition of the substrates.

Controls on isotopic composition and event based data aggregation
The high temporal resolution of N 2 O isotopic and auxiliary measurements allowed us to investigate controls on N 2 O isotopic composition over the 3 months campaign period.Correlations with isotopic composition were highest and positive for DOC and soil temperature (Table 2).The significant correlation with temperature for the whole campaign was due to a significant correlation during the "dry" part of the campaign.
If the increase in SP was due to increased contribution of nitrification, δ 15 N bulk should decrease due to the higher isotopic fractionation during this process.The simultaneous increase in SP, δ 15 N bulk and δ 18 O revealed in Fig. 7, however, indicates an increased share of N 2 O reduction to N 2 which might have been triggered by increased substrate availability (DOC) for heterotrophic denitrification.The reported effect of temperature on the N 2 O : N 2 ratio is not without any doubt, but a decrease has been observed with increasing temperature, supporting the hypothesis that N 2 O reduction increased as temperature rose throughout the measurement period (Saggar et al., 2013).
Though substrate availability has been identified as a major control on N 2 O source processes (see references in Saggar et al., 2013), correlations between N 2 O isotopic composition and NO

Short term variation of isotopic composition
The Keeling plot approach is based on conservation of mass and assumes that the atmospheric concentration of a gas in the surface layer is a mixture of background atmospheric concentration and a variable amount of gas added by a source, raising the atmospheric concentration above background.The source's isotope value can be determined given that its isotope value remains constant during the observation period.In this study, we used noon-to-noon data in the Keeling plots to determine isotope values of soil-derived N 2 O for the respective noon-to-noon period.Hence, the source processes underlying these N 2 O emissions have to be constant on this time scale.Currently, little is known about the rate of change of N 2 O source processes over time-steps of minutes to hours.However, changing relative contributions of source processes, which change the isotopic composition in soil-emitted N 2 O, would be reflected by deviations from a linear relation between inverse concentration and isotopic composition.As the Keeling plots showed no obvious deviations from a linear relation within our measurement precision, we conclude (1) that the use of the Keeling plot approach was valid in our study, and (2) that changes in N 2 O source processes in our study site occurred at a time step of one day or more.While our data suggests that there are little or no changes in source processes underlying N 2 O emissions within a noon-to-noon period, clear and distinct day-to-day variation in isotope values of soil derived N 2 O, especially in SP, were observed.Such changes were often strong and abrupt following management events (Mana I and III, Rewetting), indicating a significant response of microbial processes to the imposed disturbance.Larger than expected variability in isotope values was observed in-between management events (class BG), Introduction

Conclusions References
Tables Figures

Back Close
Full when no obvious variation in environmental drivers occurred.Since noon-to-noon concentration increases were very small during these periods, part of this variability may be attributed to increased uncertainty around the intercept of the Keeling plot.This is also reflected in the relatively large error bars around isotope values on days when N 2 O fluxes were low (Fig. 3).Alternatively some of the variation in isotope values associated with these small fluxes may result from air masses not representative of the grassland site as the concentration footprint influencing the N 2 O source signature is larger than the flux footprint (Griffis et al., 2007).

Flux weighted averages of source isotopic compositions
N 2 O isotopic composition can be used to calculate and further constrain the global N 2 O budget (Kim and Craig, 1993;Yoshida and Toyoda, 2000).The analysis of emissions from different sources such as agricultural soils or managed grasslands based on box models and isotopic composition is complicated by distinct temporal and spatial variability of isotopic composition (Kim and Craig, 1993;Toyoda et al., 2011;Yoshida and Toyoda, 2000); hence, flux weighted averages are required to obtain representative values for agricultural N 2 O (Perez et al., 2001).Our flux weighted averages of 6.9±4.3,−17.4 ± 6.2 and 27.4 ± 3.6 ‰ for SP, δ 15 N bulk and δ 18 O are well within the range of values 2.9 to 36.6, −41.5 to −1.9 and 23.2 to 51.7 ‰ for agricultural soils (Park et al., 2011;Toyoda et al., 2011), but the comparison with other grassland soils (Opdyke et al., 2009;Park et al., 2011) indicates that the variability of isotopic composition within a group, such as grassland, may be considerable (for SP: 2.2 to 11.1 ‰).Part of the variability might be also explained by a limited compatibility of laboratory results, as recently demonstrated in an inter-laboratory comparison campaign (Mohn et al., 2014).
The uncertainty in budgets derived by isotopic composition depends on the uncertainty of the representative isotopic composition for a single source, which can be reduced by a quasi-continuous measurement approach, as shown in this study.

BGD Introduction
Full    Full  Full Discussion Paper | Discussion Paper | Discussion Paper | concluded that SP is capable of distinguishing between the process groups N 2 O N (NH 2 OH-oxidation, fungal denitrification and abiotic N 2 O production; SP = 32.8 ± 4.0 ‰) and N 2 O D (nitrifier-denitrification and denitrification; SP = −1.6 ± 3.8 ‰).In addition, the intramolecular distribution of 15 N can be used as an independent validation of the global, measurement-based bottom-up N 2 O budget and has already confirmed that the isotopically light sources such as agriculture and industry contribute to the increase in atmospheric N 2 O (Toyoda et al., 2013; Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | , where δ values measured (here in 2 m height) are plotted vs. the inverse of N 2 O mole fractions.The intercept of the linear regression line can be interpreted as the isotopic composition of soil emitted N 2 O (Pataki et al., 2003).Therefore, determination of soil N 2 O isotopic composition requires an increase in N 2 O mole fraction.During the day, turbulence mixes surface layer Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

N 2 O
isotopic composition of the surface layer (lowest tens of meters above ground) air samples (n = 2130) ranged from 2.5 to 16.1, −11.9 to −2.4, 37.6 to 44.6, −4.6 to 6.6, and 14.3 to 19.3 ‰ for δ 15 N α , δ 15 N β , δ 18 O, δ 15 N bulk and SP, respectively (Fig. 2).Surface layer N 2 O mole fractions varied between 325 and 469 ppb and followed Introduction Discussion Paper | Discussion Paper | Discussion Paper | a diurnal cycle with highest values during the night when the boundary layer became more stable.Increasing N 2 O mole fractions were associated with decreasing δ values, indicating that soil emitted N 2 O that mixed into the surface layer was depleted in 15 N as compared to N 2 O in the atmospheric background.
Discussion Paper | Discussion Paper | Discussion Paper | of N 2 O isotopic signature for intensively managed grassland Simple averages of daily isotopic composition of soil-emitted N 2 O were 9.6 ± 4.4, −15.2 ± 6.0 and 29.3 ± 3 ‰ for SP, δ 15 N bulk and δ 18 O, respectively (n = 62).Representative isotopic signature for agricultural land can be estimated based on fluxweighted averages of daily signatures.For some noon-to-noon periods included in the above average, thus with an overnight increase in N 2 O mole fractions of at least 12 ppb, negative N 2 O fluxes were detected by the EC system (−0.17±2.1 nmol m −2 s −1 ; n = 14).This might be due to the uncertainty of N 2 O flux measurements, temporal averaging over positive and negative fluxes in a noon-to-noon period or different footprint regions for N 2 O flux and isotopic analysis (flux vs. concentration footprint).
Discussion Paper | Discussion Paper | Discussion Paper | where the N 2 O isotopic composition deviated only slightly from atmospheric background.A consistent decrease in δ 15 N bulk in parallel with increasing N 2 O mole fractions (accumulation of soil-derived N 2 O) confirms that the soil N 2 O source is depleted in 15 N-N 2 O relative to ambient N 2 O (Toyoda et al., 2013).A similar pattern was found for δ 18 O; an increase in N 2 O mole fraction was associated with a decrease in 18 O-N 2 O, again indicating that soil emissions were depleted in 18 O-N 2 O with respect to the atmospheric background.In contrast, Harris et al. (2014) reported a decoupling of δ 18 O and δ 15 N bulk .This may have been due to only marginal influence of soil-emitted N 2 O since the measurements were carried out in urban area and approx.95 m above the ground.Studies on N 2 O derived from Discussion Paper | Discussion Paper | Discussion Paper | and −18 and −15 ‰(Ostrom et al., 2010), covering the range of −29 to −3 ‰ observed in this study.With respect to δ 18 O, the values of 22.6 to 34.8 ‰ detected for grassland in this study are at the lower end of measurements under field conditions (4-82 ‰).Quantitative source partitioning between process groups N 2 O N and N 2 O D based on SP is possible only when no other processes except those contained in the process groups have an influence on the site-specific N 2 O isotopic composition.However, in the terminal step of denitrification, namely the reduction of N 2 O to N 2 , where N 2 O is the substrate, the lighter isotopic species is consumed, leading to an increase in SP,δ15 N bulk and δ 18 O.Consequently, part of the N 2 O originating from a combination of the two process groups, i.e.N 2 O N and N 2 O D , may have been consumed by N 2 O to N 2 reduction prior to emission.
e. N 2 O gas).Rectangles for process groups N 2 O N and N 2 O D are defined by SP values given by Decock and Six (2013) and by δ 15 N bulk values calculated based on process fractionation factors and substrate isotopic composition.For nitrification and denitrification minimum and maximum fractionation factors of −90 to −40 and −40 to −15 ‰ were assumed (Baggs, 2008), for the isotopic compositions of the N 2 O precursors (i.e., NH + 4 and NO − 3 ) a range of −20 to +10 and −25 to 15 ‰ were assumed.Koba et al. (2009) attributed a concurrent decrease in δ 15 N bulk with increasing SP values as indicative for an increasing contribution of N 2 O N .In contrast, an increase in δ 15 N bulk in parallel to increasing SP values (enrichment of 15 N in the α position relative to the β position), as observed in the present study, was allocated to a substantial increase in N 2 O reduction to N 2 .For ε 15 N bulk /εSP of N 2 O reduction, Koba et al. (2009) assumed a factor of 1.2 based on previous publications.Our results (Fig. 7) indicate that N 2 O is predominately formed by denitrification, and that deviations Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | low, except for the correlation with δ 15 N bulk .The reason might be both the number of measurment points for substrate concentrations being lower compared to other explanatory variables and substrate concentrations not necessarily reflecting process or turnover rates (Wu et al., 2012).The low explanatory power of all linear regressions underlines that drivers for N 2 O emissions are highly variable and may even change from event to event.In absence of management or rewetting events (group BG), isotopic composition covered the whole range of measured values, while management or rewetting events were characterized Discussion Paper | Discussion Paper | Discussion Paper | by lower variability in isotopic composition.Values for SP, δ 15 N bulk and δ 18 O were low for Mana I, rewetting and Mana III, whereas event Mana II showed increased SP, δ 15 N bulk and δ 18 O.This indicates that processes must have been different for Mana II, although management was almost identical.
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

Acknowledgements.
We are grateful to Hans-Ruedi Wettstein and his team for the collaboration with the ETH research station Chamau and Christoph Zellweger for support with determination of GHG mixing ratios in our target gases.Antoine Roth is acknowledged for his support during the field campaign.This project was funded by the State Secretariat for Education and Research (SER) within COST Action ES0806.The QCLAS used for EC measureents was funded by the R'Equip Project (206021 133763) by the Swiss National Science Foundation.Funding form GHG-Europe (FP7, EU contract No. 244122) and COST-ES0804 ABBA is gratefully acknowledged.Instrumental developments at Empa were supported by the Swiss National Science Foundation (SNSF).Technical support on the eddy covariance station has been provided by Thomas Baur and Peter Plüss.Preparation of N 2 O isotope standards and inter-laboratory comparison measurements were supported by the EMRP ENV52 project "Metrology for high-impact greenhouse gases".The EMRP is jointly funded by the EMRP participating countries within EURAMET and the European Union.Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |
Keeling plot derived isotopic composition of soil-emitted N 2 O ranged from 1.4 to 17.3, −29 to −3 and 22.6 to 34.8 ‰ for SP, δ Table 2).For δ 15 N bulk , correlations with ∆N 2 O, wfps, soil temperature, DOC and NO 18O correlations of f N 2 O , ∆N 2 O, precipitation, Our field observations indicate that nitrifier-denitrification and denitrification (process group N 2 O D ) dominated throughout the measurement period and that variation in isotopic composition was more likely due to variation in the extent of N 2 O reduction rather than contributions of NH 2 OH oxidation.High temporal resolution of isotopic composition in soil-emitted N 2 O showed that at the beginning of the growing season, medium wfps and low temperature induced low isotope values (representative for process group N 2 O D ), whereas in the second part of the measurement period, higher temperature and DOC stimulated N 2 O reduction to N 2 , although wfps was lower.Management or rewetting events were mostly characterized by low SP, δ 15 N bulk and δ 18 O, but the event Mana II indicated that processes underlying N 2 O emissions can vary even under similar management conditions.With this study, a new method is available that can provide real-time datasets for various single N 2 O emitting (eco)systems, such as as grasslands or agricultural soils, which will help in further constraining the global N 2 O budget based on box model calculations.

Table 2 .
Adjusted r 2 and p values for regression analysis of Keeling-plot derived isotopic compositions in soil-emitted N 2 O vs. auxiliary variables N 2 O flux (f N 2 O ), difference of maximum and minimum concentration over a noon-to-noon period (∆N 2 O), precipitation (prcp), soil moisture (wfps), soil temperature (T ) and nutrient concentrations (NO − 3 , NH + 4 and DOC).