Kinetics of N2O production and reduction in a nitrate-contaminated aquifer inferred from laboratory incubation experiments

Knowledge of the kinetics of N2O production and reduction in groundwater is essential for the assessment of potential indirect emissions of the greenhouse gas. In the present study, we investigated this kinetics using a laboratory approach. The results were compared to field measurements in order to examine their transferability to the in situ conditions. The study site was the unconfined, predominantly sandy Fuhrberger Feld aquifer in northern Germany. A special characteristic of the aquifer is the occurrence of the vertically separated process zones of heterotrophic denitrification in the near-surface groundwater and of autotrophic denitrification in depths beyond 2-3 m below the groundwater table, respectively. The kinetics of N2O production and reduction in both process zones was studied during long-term anaerobic laboratory incubations of aquifer slurries using the 15N tracer technique. We measured N2O, N2, NO3-, NO2-, and SO42- concentrations as well as parameters of the aquifer material that were related to the relevant electron donors, i.e. organic carbon and pyrite. The laboratory incubations showed a low denitrification activity of heterotrophic denitrification with initial rates between 0.2 and 13 μg N kg-1 d-1. The process was carbon limited due to the poor availability of its electron donor. In the autotrophic denitrification zone, initial denitrification rates were considerably higher, ranging between 30 and 148 μg N kg-1 d-1, and NO3- as well as N2O were completely removed within 60 to 198 days. N2O accumulated during heterotrophic and autotrophic denitrification, but maximum concentrations were substantially higher during the autotrophic process. The results revealed a satisfactory transferability of the laboratory incubations to the field scale for autotrophic denitrification, whereas the heterotrophic process less reflected the field conditions due to considerably lower N2O accumulation during laboratory incubation. Finally, we applied a conventional model using first-order-kinetics to determine the reaction rate constants k1 for N2O production and k2 for N2O reduction, respectively. The goodness of fit to the experimental data was partly limited, indicating that a more sophisticated approach is essential to describe the investigated reaction kinetics satisfactorily.DF

Areal-averaged trace gas emission rates from long-range open-path measurements in stable boundary layer conditions

Measurements of land-surface emission rates of greenhouse and other gases at large spatial scales (10 000 m<sup>2</sup>) are needed to assess the spatial distribution of emissions. This can be readily done using spatial-integrating micro-meteorological methods like flux-gradient methods which were evaluated for determining land-surface emission rates of trace gases under stable boundary layers. Non-intrusive path-integrating measurements are utilized. Successful application of a flux-gradient method requires confidence in the gradients of trace gas concentration and wind, and in the applicability of boundary-layer turbulence theory; consequently the procedures to qualify measurements that can be used to determine the flux is critical. While there is relatively high confidence in flux measurements made under unstable atmospheres with mean winds greater than 1 m s<sup>−1</sup>, there is greater uncertainty in flux measurements made under free convective or stable conditions. The study of N<sub>2</sub>O emissions of flat grassland and NH<sub>3</sub> emissions from a cattle lagoon involves quality-assured determinations of fluxes under low wind, stable or night-time atmospheric conditions when the continuous "steady-state" turbulence of the surface boundary layer breaks down and the layer has intermittent turbulence. Results indicate that following the Monin-Obukhov similarity theory (MOST) flux-gradient methods that assume a log-linear profile of the wind speed and concentration gradient incorrectly determine vertical profiles and thus flux in the stable boundary layer. An alternative approach is considered on the basis of turbulent diffusivity, i.e. the measured friction velocity as well as height gradients of horizontal wind speeds and concentrations without MOST correction for stability. It is shown that this is the most accurate of the flux-gradient methods under stable conditions

Groundwater N2O emission factors of nitrate-contaminated aquifers as derived from denitrification progress and N2O accumulation

We investigated the dynamics of denitrification and nitrous oxide (N2O) accumulation in 4 nitrate (NO3-) contaminated denitrifying sand and gravel aquifers of northern Germany (Fuhrberg, Sulingen, Thulsfelde and Gottingen) to quantify their potential N2O emission and to evaluate existing concepts of N2O emission factors. Excess N-2 - N-2 produced by denitrification - was determined by using the argon (Ar) concentration in groundwater as a natural inert tracer, assuming that this noble gas functions as a stable component and does not change during denitrification. Furthermore, initial NO3- concentrations (NO3- that enters the groundwater) were derived from excess N-2 and actual NO3- concentrations in groundwater in order to determine potential indirect N2O emissions as a function of the N input. Median concentrations of N2O and excess N-2 ranged from 3 to 89 mu g N L-1 and from 3 to 10 mg N L-1, respectively. Reaction progress (RP) of denitrification was determined as the ratio between products (N2O-N + excess N-2) and starting material (initial NO3- concentration) of the process, characterizing the different stages of denitrification. N2O concentrations were lowest at RP close to 0 and RP close to 1 but relatively high at a RP between 0.2 and 0.6. For the first time, we report groundwater N2O emission factors consisting of the ratio between N2O-N and initial NO3--N concentrations (EF1). In addition, we determined a groundwater emission factor (EF2) using a previous concept consisting of the ratio between N2O-N and actual NO3--N concentrations. Depending on RP, EF(1) resulted in smaller values compared to EF(2), demonstrating (i) the relevance of NO3- consumption and consequently (ii) the need to take initial NO3--N concentrations into account. In general, both evaluated emission factors were highly variable within and among the aquifers. The site medians ranged between 0.00043-0.00438 for EF(1) and 0.00092-0.01801 for EF(2), respectively. For the aquifers of Fuhrberg and Sulingen, we found EF(1) median values which are close to the 2006 IPCC default value of 0.0025. In contrast, we determined significant lower EF values for the aquifers of Thulsfelde and Gottingen. Summing the results up, our study supports the substantial downward revision of the IPCC default EF5-g from 0.015 (1997) to 0.0025 (2006).DF

Coordinated optimization of visual cortical maps (II) Numerical studies

It is an attractive hypothesis that the spatial structure of visual cortical architecture can be explained by the coordinated optimization of multiple visual cortical maps representing orientation preference (OP), ocular dominance (OD), spatial frequency, or direction preference. In part (I) of this study we defined a class of analytically tractable coordinated optimization models and solved representative examples in which a spatially complex organization of the orientation preference map is induced by inter-map interactions. We found that attractor solutions near symmetry breaking threshold predict a highly ordered map layout and require a substantial OD bias for OP pinwheel stabilization. Here we examine in numerical simulations whether such models exhibit biologically more realistic spatially irregular solutions at a finite distance from threshold and when transients towards attractor states are considered. We also examine whether model behavior qualitatively changes when the spatial periodicities of the two maps are detuned and when considering more than 2 feature dimensions. Our numerical results support the view that neither minimal energy states nor intermediate transient states of our coordinated optimization models successfully explain the spatially irregular architecture of the visual cortex. We discuss several alternative scenarios and additional factors that may improve the agreement between model solutions and biological observations.Comment: 55 pages, 11 figures. arXiv admin note: substantial text overlap with arXiv:1102.335

Charge separation relative to the reaction plane in Pb-Pb collisions at sNN=2.76\sqrt{s_{\rm NN}}= 2.76 TeV

Measurements of charge dependent azimuthal correlations with the ALICE detector at the LHC are reported for Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV. Two- and three-particle charge-dependent azimuthal correlations in the pseudo-rapidity range $|\eta| < 0.8$ are presented as a function of the collision centrality, particle separation in pseudo-rapidity, and transverse momentum. A clear signal compatible with a charge-dependent separation relative to the reaction plane is observed, which shows little or no collision energy dependence when compared to measurements at RHIC energies. This provides a new insight for understanding the nature of the charge dependent azimuthal correlations observed at RHIC and LHC energies.Comment: 12 pages, 3 captioned figures, authors from page 2 to 6, published version, figures at http://aliceinfo.cern.ch/ArtSubmission/node/286

A note on comonotonicity and positivity of the control components of decoupled quadratic FBSDE

In this small note we are concerned with the solution of Forward-Backward Stochastic Differential Equations (FBSDE) with drivers that grow quadratically in the control component (quadratic growth FBSDE or qgFBSDE). The main theorem is a comparison result that allows comparing componentwise the signs of the control processes of two different qgFBSDE. As a byproduct one obtains conditions that allow establishing the positivity of the control process.Comment: accepted for publicatio

Multi-particle azimuthal correlations in p-Pb and Pb-Pb collisions at the CERN Large Hadron Collider

Measurements of multi-particle azimuthal correlations (cumulants) for charged particles in p-Pb and Pb-Pb collisions are presented. They help address the question of whether there is evidence for global, flow-like, azimuthal correlations in the p-Pb system. Comparisons are made to measurements from the larger Pb-Pb system, where such evidence is established. In particular, the second harmonic two-particle cumulants are found to decrease with multiplicity, characteristic of a dominance of few-particle correlations in p-Pb collisions. However, when a $|\Delta \eta|$ gap is placed to suppress such correlations, the two-particle cumulants begin to rise at high-multiplicity, indicating the presence of global azimuthal correlations. The Pb-Pb values are higher than the p-Pb values at similar multiplicities. In both systems, the second harmonic four-particle cumulants exhibit a transition from positive to negative values when the multiplicity increases. The negative values allow for a measurement of $v_{2}\{4\}$ to be made, which is found to be higher in Pb-Pb collisions at similar multiplicities. The second harmonic six-particle cumulants are also found to be higher in Pb-Pb collisions. In Pb-Pb collisions, we generally find $v_{2}\{4\} \simeq v_{2}\{6\}\neq 0$ which is indicative of a Bessel-Gaussian function for the $v_{2}$ distribution. For very high-multiplicity Pb-Pb collisions, we observe that the four- and six-particle cumulants become consistent with 0. Finally, third harmonic two-particle cumulants in p-Pb and Pb-Pb are measured. These are found to be similar for overlapping multiplicities, when a $|\Delta\eta| > 1.4$ gap is placed.Comment: 25 pages, 11 captioned figures, 3 tables, authors from page 20, published version, figures at http://aliceinfo.cern.ch/ArtSubmission/node/87

Multiplicity dependence of jet-like two-particle correlations in p-Pb collisions at sNN\sqrt{s_{NN}} = 5.02 TeV

Two-particle angular correlations between unidentified charged trigger and associated particles are measured by the ALICE detector in p-Pb collisions at a nucleon-nucleon centre-of-mass energy of 5.02 TeV. The transverse-momentum range 0.7 $< p_{\rm{T}, assoc} < p_{\rm{T}, trig} <$ 5.0 GeV/$c$ is examined, to include correlations induced by jets originating from low momen\-tum-transfer scatterings (minijets). The correlations expressed as associated yield per trigger particle are obtained in the pseudorapidity range $|\eta|<0.9$. The near-side long-range pseudorapidity correlations observed in high-multiplicity p-Pb collisions are subtracted from both near-side short-range and away-side correlations in order to remove the non-jet-like components. The yields in the jet-like peaks are found to be invariant with event multiplicity with the exception of events with low multiplicity. This invariance is consistent with the particles being produced via the incoherent fragmentation of multiple parton--parton scatterings, while the yield related to the previously observed ridge structures is not jet-related. The number of uncorrelated sources of particle production is found to increase linearly with multiplicity, suggesting no saturation of the number of multi-parton interactions even in the highest multiplicity p-Pb collisions. Further, the number scales in the intermediate multiplicity region with the number of binary nucleon-nucleon collisions estimated with a Glauber Monte-Carlo simulation.Comment: 23 pages, 6 captioned figures, 1 table, authors from page 17, published version, figures at http://aliceinfo.cern.ch/ArtSubmission/node/161

Transverse sphericity of primary charged particles in minimum bias proton-proton collisions at s=0.9\sqrt{s}=0.9, 2.76 and 7 TeV

Measurements of the sphericity of primary charged particles in minimum bias proton--proton collisions at $\sqrt{s}=0.9$, 2.76 and 7 TeV with the ALICE detector at the LHC are presented. The observable is linearized to be collinear safe and is measured in the plane perpendicular to the beam direction using primary charged tracks with $p_{\rm T}\geq0.5$ GeV/c in $|\eta|\leq0.8$. The mean sphericity as a function of the charged particle multiplicity at mid-rapidity ($N_{\rm ch}$) is reported for events with different $p_{\rm T}$ scales ("soft" and "hard") defined by the transverse momentum of the leading particle. In addition, the mean charged particle transverse momentum versus multiplicity is presented for the different event classes, and the sphericity distributions in bins of multiplicity are presented. The data are compared with calculations of standard Monte Carlo event generators. The transverse sphericity is found to grow with multiplicity at all collision energies, with a steeper rise at low $N_{\rm ch}$, whereas the event generators show the opposite tendency. The combined study of the sphericity and the mean $p_{\rm T}$ with multiplicity indicates that most of the tested event generators produce events with higher multiplicity by generating more back-to-back jets resulting in decreased sphericity (and isotropy). The PYTHIA6 generator with tune PERUGIA-2011 exhibits a noticeable improvement in describing the data, compared to the other tested generators.Comment: 21 pages, 9 captioned figures, 3 tables, authors from page 16, published version, figures from http://aliceinfo.cern.ch/ArtSubmission/node/308

Anisotropic flow of charged hadrons, pions and (anti-)protons measured at high transverse momentum in Pb-Pb collisions at sNN=2.76\sqrt{s_{\rm NN}}=2.76 TeV

The elliptic, $v_2$, triangular, $v_3$, and quadrangular, $v_4$, azimuthal anisotropic flow coefficients are measured for unidentified charged particles, pions and (anti-)protons in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV with the ALICE detector at the Large Hadron Collider. Results obtained with the event plane and four-particle cumulant methods are reported for the pseudo-rapidity range $|\eta|<0.8$ at different collision centralities and as a function of transverse momentum, $p_{\rm T}$, out to $p_{\rm T}=20$ GeV/$c$. The observed non-zero elliptic and triangular flow depends only weakly on transverse momentum for $p_{\rm T}>8$ GeV/$c$. The small $p_{\rm T}$ dependence of the difference between elliptic flow results obtained from the event plane and four-particle cumulant methods suggests a common origin of flow fluctuations up to $p_{\rm T}=8$ GeV/$c$. The magnitude of the (anti-)proton elliptic and triangular flow is larger than that of pions out to at least $p_{\rm T}=8$ GeV/$c$ indicating that the particle type dependence persists out to high $p_{\rm T}$.Comment: 16 pages, 5 captioned figures, authors from page 11, published version, figures at http://aliceinfo.cern.ch/ArtSubmission/node/186