1 Introduction
Nitrogen oxides () are pollutant trace gases that play a key role in the atmosphere by producing or destroying tropospheric ozone (), which can cause respiratory illness (Ciencewicki and Jaspers, 2007) and damage to plants (Emberson et al., 2018). Photolysis of nitrogen dioxide () (Reaction R1) is the primary source of tropospheric ozone (), and the nitric oxide (NO) product is oxidized back to by either (1) (Reaction R2) or (2) organic peroxy radicals (, under formation of alkoxy radicals (RO)) or hydroperoxyl radicals () (Reactions R3 and R4) (Lightfoot et al., 1992). The latter results in the formation of the hydroxyl radical (OH); thus, Reactions (R3) and (R4) represent routes to recycle the most important atmospheric radical initiator of oxidation (Hens et al., 2014). It is, therefore, essential to understand the sources and sinks of in the atmosphere. The dominant global sources of are anthropogenic in the form of the combustion of fossil fuels and, to a lesser degree, biomass burning and agricultural soils. The natural sources, which include lightning (Schumann and Huntrieser, 2007), wildfires (Val Martin et al., 2008), and unperturbed soil emissions from microbial activities (Davidson and Kingerlee, 1997), are important in regions remote from anthropogenic sources. and NO both react with peroxy radicals in the atmosphere to produce organic nitrates (Reactions R5 and R6), including peroxy nitrates () and alkyl nitrates (), which are important precursors for the formation of secondary organic aerosols (SOAs) (Hallquist et al., 2009; Kanakidou et al., 2005; Kiendler-Scharr et al., 2016). also reacts with OH radicals, , and nitrate radicals (), to form nitric acid () (Reaction R7); radicals (Reaction R8); and dinitrogen pentoxide () (Reaction R9), respectively. is in thermal equilibrium with , and and can interact with aqueous aerosol or moist surfaces to form (Reaction R10) (Kane et al., 2001) or nitryl chloride () (Phillips et al., 2013, 2012). Organic nitrates, SOAs, and are all removed from the boundary layer through dry and wet deposition, which also, in turn, removes from the atmosphere. In the planetary boundary layer, is also lost through dry deposition to surfaces, such as soil and leaves. Deposition takes places during both nighttime and daytime, but it is expected to be more efficient during daytime due to increased mixing through turbulence. When deposits onto humid surfaces, it can lead to the production of nitrous acid (HONO), which can be released to the atmosphere (Meusel et al., 2016; Elshorbany et al., 2012). uptake on leaves takes place through stomatal and non-stomatal processes, which have been reported to depend on multiple factors, such as the stomatal aperture and relative humidity. Stomatal uptake primarily occurs during daytime, when the stomata are open, leading to increased loss compared with nighttime, when the stomata are not fully open (Delaria et al., 2020, 2018; Chaparro-Suarez et al., 2011). Non-stomatal uptake occurs through the cuticles, although the importance of cuticular uptake has been reported to be small compared with the stomatal uptake (Delaria and Cohen, 2020; Delaria et al., 2020). uptake by leaves is reported to be enhanced in the presence of water films, which may exist when the relative humidity is % (Thoene et al., 1996; Weber and Rennenberg, 1996; Burkhardt and Eiden, 1994). There is, however, no consensus on this process, as other studies have not observed this effect (Gessler et al., 2000). Most recent work shows that the interactions with foliar surfaces are unidirectional, i.e. emissions are negligible (Delaria et al., 2020).
At nighttime, photolysis ceases; as a consequence, in the absence of combustion sources, the main sources of NO are emissions from soils (Jaeglé et al., 2005). As NO is oxidized efficiently by at night, its concentration will be highest at the surface and will decrease with altitude. The vertical profile of is the opposite, owing to its physical loss due to deposition near the surface and through chemical reaction with NO and/or alkenes combined with entrainment from the nocturnal residual layer. As is produced from the reaction between NO and , its vertical gradient is expected to be weaker than those of NO and (Geyer and Stutz, 2004; Stutz et al., 2004).
In this study, we use measurements from the ACROSS (Atmospheric ChemistRy Of the Suburban foreSt) campaign to investigate the nighttime sources and sinks of in a temperate forest. measurements are used to explain the observed NO features, and measurements of and total gas-phase nitrogen species () and particulate nitrate are used to investigate the lifetime and fate of in the forest environment.
2 The ACROSS campaign
The ACROSS campaign (13 June 2022 to 25 July 2022) was conducted in multiple locations in and around Paris, France (Cantrell and Michoud, 2022). Here, we present measurements from the Rambouillet Forest supersite located approximately 50 km southwest of Paris (48.687° N, 1.704° E). The forest consists of approximately 70 % oak, 20 % pine, and small fractions of beech and chestnut. The top of the forest canopy around the supersite is around 20–25 m. In this work, several instrumented containers were placed in a clearing ( ) along with a 41 m measurement tower. Most of the instruments used in this study were located in two different containers (MPIC and Orléans). The sampling inlets of the two containers were approximately 17 m apart, and the tower was approximately 9 m from the MPIC container and 16 m from the Orléans container. The soil measurements were carried out at the bottom of the tower, approximately 13 m from the MPIC container and approximately 17 m from the Orléans container. All of the instruments used in this study are described briefly below.
2.1 Measurements
2.1.1 Ground measurements
was measured using two different cavity ring-down spectroscopy (CRDS) instruments with co-located inlets sampling from a high-volume-flow stainless-steel tube (10 ; 15 cm diameter and 0.2 s residence time) that took air from a height of 5.4 m above the ground. One of the instruments (5CH-CRDS) consists of three cavities, operated at 408 nm, to measure and, via their thermal dissociation to , total peroxy nitrates (PNs, 448 K) and total alkyl nitrates (ANs, 648 K). Two additional cavities, operated at 662 nm, measured and (via thermal dissociation to at 373 K) (Sobanski et al., 2016). During this campaign, the cavity had a limit of detection (LOD) of 9.7 pptv (parts per trillion by volume) for a 1 min averaging time (3). The second instrument (k-NO3) primarily measured the reactivity, but it also had a cavity operated at 405 nm for the measurement of (Liebmann et al., 2018).
Another CRDS instrument was used to measure , , and particulate nitrate (p) from co-located inlets near the high-volume-flow stainless-steel tube. was measured by adding to the ambient sample, thereby oxidizing NO to , the latter of which was measured with CRDS at 405 nm (Friedrich et al., 2020). A judicious choice of and reaction time ensured that minimal ( %) was oxidized to . At times with low (or zero) NO, concentrations were in close agreement with both measurements. was measured by passing ambient air through a quartz inlet at K; this quantitatively converts reactive nitrogen trace gases to NO or . Exceptions are , HCN, and , which are not detected. At this location, is expected to mainly consist of .
Particulate nitrates (both organic and inorganic) were separately measured (as ) after denuding gas-phase reactive nitrogen species (Friedrich et al., 2020). To achieve this, problems involving the ineffective trapping of gas-phase by the denuder were eliminated, as will be described in a forthcoming technical paper.
was measured from the high-volume-flow stainless-steel tube with a commercial instrument (Model 205, 2B Technologies) using UV absorption at 254 nm. The LOD is 2 ppbv for a 10 s averaging time.
A spectral radiometer (METCON GmbH) was installed near the co-located inlets on top of the MPIC container to measure actinic fluxes, which were used to calculate photolysis frequencies, as described elsewhere (Meusel et al., 2016).
NO was measured from the Orléans container using a commercial chemiluminescence instrument (CLD 780 TR, Eco Physics; hereafter CLD) with a LOD of 10 pptv for a 1 min averaging time. The sampling height for NO measurements was about 0.6 and 3.2 m above the container top and the ground surface, respectively. The NO measurements required correction due to a change in the CLD sensitivity during the campaign caused by an interruption in the instrument's oxygen supply. The corrections and the corrective procedure are described in the Supplement.
HONO was measured by a commercial long-path absorption photometer (LOPAP-03, QUMA GmbH) with a sampling height of 2.0 m above ground level. Details about the LOPAP instrument can be found elsewhere (Heland et al., 2001; Kleffmann et al., 2006). During the campaign, the LOPAP was calibrated by diluted nitrite when any supporting solutions were changed. Zero calibration by measuring synthetic air was conducted two to three times per day. The detection limit was pptv.
The sum of peroxy radicals () was measured by their conversion to in the presence of NO and and detection of the generated using chemical ionization mass spectrometry (Kukui et al., 2008). The calibration coefficient is determined using actinometry and both OH and generation in a turbulent-flow reactor by photolysis of or at 184.9 nm. The calibration of , , and other species is performed by adding CO, , or other precursors into the calibration reactor and converting OH to . The overall estimated calibration accuracy (2) for is about 30 %, although the uncertainty of the measurements is typically higher due to uncertainty in the ambient-air composition. The lower limit of detection for radicals at a signal-to-noise ratio of 3 and a 4 min integration time is .
Time series of the most relevant measurements can be found in Figs. S1 and S2 in the Supplement. Due to missing total and measurements prior to 25 June and missing NO measurements after 18 July, the data analysis is focused on the time period in between these dates.
2.1.2 Tower measurements
Measurements at 41 m were conducted with instruments located on the tower as well as through a manifold with an inlet at the top of the tower. The manifold was built from BORODRAIN glass tubing (4.9 cm inner diameter) with a residence time in the manifold of 2.1 s. was measured using a cavity attenuated phase shift (CAPS) instrument on the tower with a LOD of 40 pptv; this instrument was zeroed every 1–2 h. NO and were both measured from the manifold using a chemiluminescence instrument with a LOD of 30 pptv and a HORIBA (APOA370) with a LOD of 2.5 ppbv, respectively. The NO measurements were corrected for losses due to the reaction of NO with in the manifold and the sampling line (total of 5.5 s), with corrections ranging from 1 % to 28 %. Time series of all three measurements are plotted in Fig. S3 in the Supplement.
2.1.3 Meteorology and soil measurements
Ambient temperature was measured at four different heights on the tower – 5, 13, 21, and 41 m – using temperature sensors from ATEXIS (PT1000) and Thermo Est (PT100). Relative humidity was measured at 5 m using a Vaisala humidity sensor (HMP45A). Soil temperature and moisture were measured at 5, 10, and 30 cm below the surface using probes from Thermo Est (PT100) and Delta-T Devices (ThetaProbe ML2X), respectively. Wind speed and direction were measured at 41 m using a wind monitor from Young Company. Time series of all of the meteorological and soil measurements are shown in Figs. S4 and S5 in the Supplement.
2.2 HYSPLIT
To identify different air masses, 48 h back trajectories were simulated every hour at a terminating height of 40 m using the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model (version 4, 2019; Draxler and Rolph, 2011). The back trajectories were modelled using meteorological data from the Global Data Assimilation System (GDAS) at a resolution of 1°. This led to the separation of the data into two periods, 25 June to 2 July and 3 July to 18 July, which are plotted in Fig. 1. The first phase was dominated by clean air from over the Atlantic Ocean (hereafter referred to as “Atlantic”). Back trajectories indicated that the vast majority of the air masses were transported within the boundary layer prior to reaching the site and, thus, may have a reasonably fresh “marine influence”. The second phase was dominated by air that had passed over urban locations, including Paris, Brussels, and the Ruhr area, within the prior 48 h (hereafter referred to as “Continental”).
Figure 1
The 48 h back trajectories from the Rambouillet Forest supersite using the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model (version 4, 2019).
[Figure omitted. See PDF]
3 Results and discussionTwo 24 h periods of temperature (at four different heights), NO, , relative humidity (RH), , and the photolysis rate constant () are plotted in Fig. 2. The left column shows the 24 h period with Atlantic air, whereas the right column shows the 24 h with Continental air. Immediately apparent in these datasets (and in Fig. S1) is the large diel cycle in the mixing ratios, with net daytime production resulting in mid-afternoon mixing ratios of between and 90 ppbv. In contrast, very low mixing ratios (often approaching zero) were observed at nighttime.
Figure 2
Measurements of temperature, NO, , RH, , and for two different nights during the campaign: one during the Atlantic phase (a, c, e, g, i) and one during the Continental phase (b, d, f, h, j). The four different colours in the legend symbolize four different heights: red represents 3–6 m, grey represents 13 m, orange represents 21 m, and black represents 41 m. Blue shows the RH at 5 m. The grey shaded areas in the plots show the time during which the MPIC container was in the shade in the afternoon and morning (light grey) and nighttime (dark grey). The vertical dashed black lines indicate the beginning and end of the observed temperature inversions in panels (a) and (b).
[Figure omitted. See PDF]
In Fig. 2i and j ( measurements), nighttime is shown in dark grey, while the two light-grey areas represent the time before sunset (about 5 h) and after sunrise (about 4.5 h) when very little direct sunlight reaches the ground at the site due to shading by the trees. This leaves about 6.5 h centred around midday when direct sunlight reaches the ground. The shading results in radiative cooling of the ground in the late afternoon, and associated temperature inversions begin to form prior to sunset, as can be observed in the right column of Fig. 2 and in more detail in Fig. S6 in the Supplement. The temperature inversions begin at approximately the same time as the ground temperature at 5 cm below the surface starts to decrease (see Fig. S6). These conditions of insolation were relatively consistent throughout the campaign.
Clear temperature inversions were observed for both of the nights shown in Fig. 2, the beginning and end of which are indicated by dashed lines. Vertical mixing can be significantly reduced during a temperature inversion, which is apparent from the and RH measurements in the right-hand column. In both examples, decreases at ground level (5.4 m) at the beginning of the temperature inversion and increases as the inversion breaks down in the morning. This behaviour is understood in terms of loss to soil surfaces and through stomatal and non-stomatal uptake on leaves (Zhou et al., 2017; Rannik et al., 2012; Altimir et al., 2006; Ganzeveld and Lelieveld, 1995) as well as through chemical reactions with, for example, NO, , and unsaturated (biogenic) organics (Kurpius and Goldstein, 2003). Reduced vertical mixing means that, during the inversion, is only slowly replenished by downward mixing of air masses above the canopy, where higher levels are observed. In contrast, the RH behaves in the opposite sense, as the air above the inversion is drier than that close to the ground, where evapotranspiration contributes to enhanced water vapour concentrations.
If the only source of NO was the photolysis of , NO mixing ratios would be expected to follow the photolysis rate during the day and tend toward zero at night, as NO is oxidized to by on a timescale of minutes (for ppb). This was not always the case during ACROSS. A pronounced NO peak (up to ppbv) was observed at ground level between 00:00 and 06:00 UTC (02:00 and 08:00 LT) during the phase dominated by Atlantic air (shown in Fig. 2); this peak is absent in the phase dominated by Continental air. The peak occurs prior to sunrise and is only observed by the ground-level measurements, suggesting a non-photolytic source of NO close to the ground; this is discussed further below. Very low (0–5 ppbv) mixing ratios coincide with the sustained nighttime NO peak observed, which is never reached in the Continental-phase example, although clear temperature inversions were seen in both cases. Additional examples of sustained NO peaks (i.e. lasting several hours at a level of between 1 and 2 ppbv) at night during the first phase are shown in Fig. S7 in the Supplement. Examples of additional nights with temperature inversions during phase 2, during which NO mixing ratios remained close to zero, are shown in Fig. S8 in the Supplement.
3.1 Nighttime ozone lossFor each night between 17 June and 22 July, the net loss-rate constant, (), was derived by fitting exponential expressions to the data for periods of 4.5 to 8 h. () was highly variable, with values between and , depending on the strength of the temperature inversion and the relative humidity (see discussion below). These values of () correspond to lifetimes of 1–15 h for at nighttime. Chemical losses of occur through reactions with NO, , and unsaturated biogenic volatile organic compounds (BVOCs) (Zhou et al., 2017). Rate coefficients of reactions of with NO ( at 298 K), ( at 298 K), limonene (a reactive terpene, at 298 K), and -caryophyllene (a sesquiterpene, at 298 K) are low; thus, mixing ratios in excess of 1 ppbv would be required for NO and -caryophyllene to explain the loss-rate constant (IUPAC, 2024). The required mixing ratios of terpenes or would be even larger (60–300 ppbv). As such high mixing ratios of NO and were not observed continuously and such BVOC levels are unlikely, we assume that chemical losses of are insignificant compared with deposition, as previously observed (Zhou et al., 2017). Ignoring entrainment from other heights, we can then equate () to (), where is the deposition velocity and is the boundary layer height; the factor 2 is used to account for a positive vertical gradient (Shepson et al., 1992). Using a boundary layer height of 20 m (arbitrarily set equal to the top of the canopy) gives net deposition velocities of between 0.018 and 0.3 . These values for are in broad agreement with other studies in temperate forests, where deposition velocities for at nighttime have been reported to be around 0.07–0.3 (Padro, 1996, 1993; Finkelstein et al., 2000; Wu et al., 2016).
In Fig. 3, the production rate (), RH, temperature at four different heights, and mixing ratio are plotted for two nights with high average RH to illustrate the impact of temperature inversions on the net loss-rate constants. The production rate of is used to identify periods in which production is negligible. In the left column, a night without a temperature inversion is plotted, where the average RH for the period used to fit the exponential decay is 93 % 3 %. These conditions resulted in a net loss-rate constant of . In contrast, the night depicted in the right column has the same average RH (92 % 3 %) and a very clear temperature inversion, which gives a net loss-rate constant of . This gives a factor of 5 between these two net loss-rate constants depending on whether a temperature inversion is observed or not. This can be understood in terms of the being replenished from above when there is no (or a weak) inversion, which is not the case when there is an inversion. Bearing this in mind, the use of () (a net loss constant) must result in a lower limit to unless strong temperature inversions (preventing entrainment from above) are present. The loss rate will also be enhanced under conditions of strong inversion if trace gases that are reactive towards are released into a very shallow boundary layer. However, as indicated above, chemical losses are not expected to compete with physical losses.
Figure 3
The production of (), temperature, RH, and plotted for two nights with high average RH: one without a temperature inversion (a, c, e) and one with a temperature inversion (b, d, f). The four different colours in the legend symbolize four different heights: red represents 5–5.5 m, grey represents 13 m, orange represents 21 m, and black represents 41 m. Blue shows the RH at 5 m. The net nighttime loss is fitted with an exponential decay curve (solid black line) in panels (e) and (f). The grey shaded areas represent nighttime.
[Figure omitted. See PDF]
To investigate the impact of RH on the net loss-rate constants, two nights with temperature inversions are plotted in Fig. S9 in the Supplement: one with a high RH (92 % 3 %) and one with a lower RH (63 % 6 %). Here, we see a large decrease in (), from to , when going from a high to a lower RH. The individually determined loss-rate constants are plotted as a function of RH in Fig. 4 and colour-coded depending on whether a temperature inversion is observed during the time period that was used for the exponential decay fit. A clear increase in loss-rate constants can be observed when the RH increased above 70 %–80 % when a temperature inversion was observed. A small increase at a RH higher than 70 %–80 % was also observed when temperature inversions were absent. The observed dependence of () on relative humidity is consistent with previous studies in forested regions, which have reported an increase in loss above 60 %–70 % RH (Altimir et al., 2006; Rannik et al., 2012; Zhou et al., 2017). Altimir et al. (2006) suggested an enhancement factor that is humidity dependent above 70 % RH: 1 at 70 % RH, 2 at 85 % RH, and a sharp increase to over 5 when moving towards 100 % RH. In a boreal forest, these observations have been explained by the formation of a “wet skin” on leaves that enhances surface losses by modifying (reducing) the surface-resistance to uptake (Zhou et al., 2017). This is in broad agreement with our observations during nights with a temperature inversion (see Fig. 4), and the discrepancies between the studies could be explained, for example, by different tree types, the height of the boundary layer, the strength of the inversion, and the temperature.
Figure 4
Net loss-rate constants at 5.4 m plotted against the average relative humidity measured during the time used to fit the exponential decay of . The error bars represent on the average RH. The dashed line symbolizes the observations made by Altimir et al. (2006).
[Figure omitted. See PDF]
The faster net rate of loss on nights with high RH and well-defined temperature inversions explains the differences observed in the mixing ratios at night during the Atlantic and Continental phases. The average nighttime (20:00–04:00 UTC) RH for the Atlantic phase was 87.4 % 7.6 % (1) compared with 68.4 % 12.7 % (1) for the Continental phase, indicating that higher loss-rate constants would be expected for the Atlantic phase on nights with temperature inversions. The high RH combined with the significantly lower average peak mixing ratio in the Atlantic phase (34.5 6.0 (1) ppbv between 14:00 and 15:00 UTC) compared with the Continental phase (52.7 13.6 (1) ppbv between 14:00 and 15:00 UTC) explains why the was essentially completely depleted on nights with temperature inversions during the Atlantic phase, as shown in Figs. 2 and S7.
3.2 Nitrogen oxide soil emissionsFigures 2c and S7 show nighttime periods in which NO was observed when was depleted during the Atlantic phase. The several-hour duration of the period during which NO was above the LOD excludes very local combustion as the source, leaving soil emissions resulting from microbial activity (Davidson and Kingerlee, 1997) as the most likely source of NO. At 293 K and 2 ppbv of , the lifetime of NO with respect to the reaction with is around 20 min. It is, therefore, reasonable to assume that NO is close to steady state when there is 2 ppbv or more of available. The NO emission rate () can, thus, be equated to the loss rate of NO, as described in Eq. (1), assuming that all peroxy radicals () react with the same rate coefficient as :
1 where and are the temperature-dependent rate constants for the reaction of NO with and , respectively (IUPAC, 2024), and [NO], [], and [] are the measured concentrations of NO, , and , respectively. In Fig. 5, NO and (when ppbv) at nighttime ( ) are separated by air masses and plotted against , where the outliers are defined as being outside 1.5 the interquartile range (IQR). While the nighttime NO mixing ratio increased rapidly when tended towards 0 ppbv during the Atlantic phase, was never depleted to less than 5 ppbv during the Continental phase and, therefore, no sustained periods of enhanced NO were observed at nighttime. In contrast, no significant trend is found when plotting against for either of the phases, which shows that the calculated soil emission of NO is not dependent on . This indicates that, while the soil is an important but variable source of NO, sustained nighttime NO peaks are only observed above the instrument LOD when is almost totally depleted so that the lifetime of NO is long enough to allow its concentration to build up sufficiently.
Figure 5
The NO (a) and NO emission (b) plotted against in a box and whiskers plot, where the outliers are defined as being outside . The colours represent the two different air masses.
[Figure omitted. See PDF]
Water content and temperature have previously been shown to impact the emission rate of NO from soil (Pilegaard, 2013; Rosenkranz et al., 2006). Rosenkranz et al. (2006) found a positive correlation between soil moisture and NO emission up to a 40 % water-filled pore space and an optimum between 12.5 and 15 soil temperature in a sessile oak forest in Hungary. In Fig. 6, NO and are plotted against the soil temperature and moisture at 5 cm below the surface. The measured NO mixing ratios peak towards the highest soil moisture and lowest soil temperature measured during this campaign; however, as with , there is no significant trend in the NO emission rates with soil moisture. At the low (11.5–12.5 ) and high (19.5–20.5 ) nighttime soil temperatures, very few measurements were made (around 2 h combined) compared with the rest of the temperature intervals. Across the remaining temperature intervals, no significant trend was observed in the estimated NO emission.
Figure 6
The NO (a, b) and NO emission (c, d) plotted against soil moisture (a, c) and temperature (b, d) at 5 cm below the surface in a box and whiskers plot, where the outliers are defined as being outside . The colours represent the two different air masses.
[Figure omitted. See PDF]
The average NO emission rate derived for the two phases is identical, with values of 1.45 1.61 (1, median of ) and 1.42 5.68 (1, median of ) for the Atlantic and Continental phases, respectively, when using data where ppbv. The Continental phase shows much higher variability, resulting from more spikes in the data during that period. When is completely depleted during the Atlantic phase, the increase in NO per hour results in NO emission rates of 0.3–1.8 , which is in reasonable agreement with the averages across each of the two phases when there is still present. By assuming a mixed nocturnal boundary layer (NBL) with a height of 20 m (top of the canopy), the average emission rates can be converted to NO emission fluxes of 16.6 18.5 (1) and 16.2 65.0 (1) , respectively. These values are within the range of previous measurements in different European forests with similar tree types to those found in the Rambouillet Forest (see Table 1). The measurements by Pilegaard et al. (2006) and Rosenkranz et al. (2006) were all performed using the chamber technique, whereas Schindlbacher et al. (2004) measured the emission from soil samples collected in the field and exposed to different temperature and humidity levels in the laboratory. The chamber-derived emission rates are all either lower than or, within combined uncertainties, equal to the values determined in this study, while emission rates from the soil samples were higher than or, within combined uncertainties, equal to the values derived in the present study. Davidson and Kingerlee (1997) modelled the global NO emission inventory from soil depending on the biome (e.g. temperate forest, agriculture, and savanna) and split the temperate forest category into regions affected or unaffected by nitrogen deposition. For temperate forests not affected by nitrogen deposition, those authors estimated a flux of 0.0–0.2 (0.0–2.3 ), which is in good agreement with the lower measurements by Pilegaard et al. (2006). In contrast, the temperate forests impacted by nitrogen deposition had estimated fluxes of 1.1–5.0 (12.6–57.1 ), which is in good agreement with our measurements at Rambouillet Forest, where nitrogen deposition is enhanced by pollution arriving from Paris and other surrounding urbanized/industrialized areas. While noting that our fluxes are broadly consistent with previous measurements, we recognize that the calculations are based on the assumptions of a well-mixed boundary layer of fixed height (arbitrarily set at 20 m) and should not be overinterpreted.
Table 1Measured NO soil emission in European forests with the same tree types as in the Rambouillet Forest.
Dominant tree type | Location | NO emission () | Reference |
---|---|---|---|
Oak | Mátra Mountains, Hungary | 2.1 | Pilegaard et al. (2006) |
(summer) | Rosenkranz et al. (2006) | ||
(autumn) | Rosenkranz et al. (2006) | ||
Pine | San Rossore, Italy | 5.4 | Pilegaard et al. (2006) |
Beech | Schottenwald, Austria | Schindlbacher et al. (2004) | |
4.2 | Pilegaard et al. (2006) | ||
Beech | Klausen-Leopoldsdorf, Austria | Schindlbacher et al. (2004) | |
0.7 | Pilegaard et al., 2006) | ||
Spruce–fir–beech | Achenkirch, Austria | Schindlbacher et al. (2004) | |
0.9 | Pilegaard et al. (2006) | ||
Mixed deciduous | Ticino Park, Italy | Schindlbacher et al. (2004) | |
Below LOD | Pilegaard et al. (2006) |
At nighttime, in the absence of its photolysis, may be expected to increase in concentration (via Reaction R2) when a constant NO source exists (e.g. from soil, as observed here) and when is present. For both the Atlantic and the Continental phases, an average diel profile, between 20:00 and 04:00 UTC, of (black) is plotted in Fig. 7. No obvious increase in can be observed in the Atlantic phase, whereas an average increase of around 1 ppbv can be observed in the Continental phase. The expected resulting from the reaction if there were no loss mechanisms of is plotted in red. This is determined by using the measured at 20:00 UTC and incrementally altering this value by the that would have been produced through NO oxidation by and peroxy radicals at each time step. In both phases, the simple assumptions of nighttime production through and and no loss result in the significant generation of with an overestimation of 10–12 ppbv of at the end of the night compared with the measured . Thus, a loss mechanism of around 1.4 of is necessary to explain the observed (lack of increase in) .
Figure 7
Average nighttime profiles of at 5.4 m for each of the two phases (black) plotted along with the expected with (grey) and without (red) loss.
[Figure omitted. See PDF]
3.3.1 Chemical lossesWhile is removed in a largely irreversible process through reaction with OH radicals to form during the daytime, this is unlikely to represent a significant sink at nighttime. In the absence of photochemical formation pathways, OH is generated at night via the ozonolysis of olefins and in the reaction of with and NO. In the forested environment in summer, emissions of BVOCs (e.g. olefinic terpenoids) will favour Reaction (R11) and, simultaneously, disfavour Reaction (R12), as will be reduced with respect to its concentration via its reactions with BVOCs. During the ACROSS campaign, ground levels were generally below the instrument detection limits of 2 pptv; hence, we can reasonably ignore Reaction (R12). Measurements of OH in forested environments are sparse, although they indicate that nocturnal OH levels are low, with concentrations generally lower than . Combining the rate coefficient for the reaction of OH with of (IUPAC, 2024) at ambient pressure and K with an upper limit (confirmed by measurements) to the OH concentration of results in an loss constant of , or (at the average nighttime concentration of pptv) a loss rate of , which is clearly insufficient to explain the observations.
is also lost via its reaction with to form the radical (Reaction R8). We have shown (Dewald et al., 2024) that the majority of NO3 formed in the forest will react with BVOCs rather than with NO (to re-form NO2) and, to a good approximation, that Reaction (R8) represents an irreversible loss of NO2, as the alkyl nitrates will not release nitrogen in the form of NO2 at nighttime. However, the rate coefficient for this process ( at 298 K; IUPAC, 2024) is very small; therefore, with average nighttime levels reduced by deposition (see above) to 23 ppbv, the lifetime of with respect to this reaction is 14 h and the loss rate (at the average nighttime concentration of pptv) is , which is again too slow to contribute significantly to the apparent loss rate of .
The chemical loss of via reaction with OH or via formation of and its further reactions with BVOCs to form alkyl nitrates is expected to result in the conversion of to . As described in Sect. 2.1.1, we operated an instrument to measure in both the gas- and particle-phases during the ACROSS campaign. Figure 8 displays the average diel profiles of (–) and p during the Atlantic and Continental phases. For both and p, the diel profiles show either a decrease or a stable mixing ratio across the period in which losses of 10–12 ppbv of are required to explain the observations. Clearly, the loss of at nighttime is not balanced by the formation of other forms of reactive nitrogen that were long-lived enough to be detected. Trace gases, such as or alkyl nitrates, may be lost via deposition to surfaces, especially at high relative humidity, and lifetimes for biogenic alkyl nitrates of a few hours have been reported (Liebmann et al., 2019; Farmer and Cohen, 2008; Browne et al., 2013; Romer Present et al., 2020). However, as shown above, the limiting step in the formation of organic nitrates is the slow reaction of with , which will not convert sufficient to to explain our observations. The formation of organic nitrates that do not require the intermediacy of (i.e. peroxy nitrates formed from ) would also have been detected by the instrument and can, thus, also be ruled out as major reservoirs of .
Figure 8
Average diel profiles of NO, , total gas-phase , and particulate nitrate (p) at 3–6 m above ground for the Atlantic (a, c, e, g) and Continental (b, d, f, h) phases. The grey shaded areas symbolize nighttime.
[Figure omitted. See PDF]
deposited to humid surfaces can be converted to HONO and released to the atmosphere (Elshorbany et al., 2012; Meusel et al., 2016). A time series of HONO can be found in Fig. S2 and reveals increases in HONO at nighttime. However, the HONO mixing ratios can account for only a small fraction of the loss described above. This may reflect the fact that, if formed at a moist surface, (soluble) HONO is unlikely to desorb quantitatively into the gas phase. The low HONO mixing ratios measured during the Atlantic phase, compared with the Continental phase, could potentially be explained by the difference in soil humidity; however, the factors influencing the formation and release of HONO are complex. The HONO observations will be analysed in detail in a separate publication from the ACROSS campaign.
In the absence of other known gas-phase mechanisms for the removal of at night and owing to the fact that very small amounts of other reactive nitrogen trace gases or nitrate particles are formed during nighttime, we conclude that the physical removal of (i.e. deposition) is responsible for its lack of build-up at night during ACROSS and that any transformation of at the surface does not lead to a quantitative release into the gas phase.
3.3.2 Physical lossesis known to be lost through dry deposition to surfaces such as soil and leaves, with the latter depending on whether the stomata are open (daytime) or not fully open (nighttime) (Delaria et al., 2020, 2018; Chaparro-Suarez et al., 2011). As for , the dry deposition of to surfaces can be described by an exponential decay with a first-order decay rate constant of , where is the deposition velocity and is the boundary layer height. This expression applies when gradients within the boundary layer are weak, as expected for (see above), even though vertical mixing is very slow at night. The net production (or loss) of is given by Eq. (2), where the first term on the right-hand side is the production rate from the reaction of NO with or (which is identical to the NO soil emission rate) and the second term is the loss rate assuming only depositional losses (see above) and ignoring the entrainment of from other heights. This will give an upper limit for the deposition rate, as a small fraction ( %) of is lost through chemical reactions with and OH (see above).
2 where []0 is the mixing ratio at 20:00 UTC. The concentration at any subsequent time can then be calculated as described in Eq. (3) via variation of the () term in order to match the observed mixing ratio. 3 In Fig. 7, the grey lines symbolize the calculated mixing ratios at nighttime using values of () between and . As expected, no single value of () can explain all of the measurements, as the height of the boundary layer will not be invariant during the whole night. However, for the Continental and Atlantic phases, the observed can be explained with and , respectively, resulting in respective lifetimes of –3 h and –110 min for at nighttime. As the deposition of in this environment represents a permanent loss of from the gas phase, this lifetime can be compared to, for example, the lifetime of with respect to its conversion to via the reaction of with OH, which is d (assuming an average of ). The low aerosol surface area during ACROSS combined with the low uptake coefficient for renders losses due to heterogeneous processes insignificant (IUPAC, 2024). Clearly, only depositional losses of in a forested environment contribute substantially to its lifetime at night and to the budget.
If we continue to assume that the NBL at the forest site is at the top of the canopy (20 m), the loss-rate constants that we determined can be converted to a deposition velocity of 0.4 0.2 and 0.55 0.25 for the Continental and Atlantic phase, respectively. These are comparable to previous measurements of deposition velocities of 0.15 (Dewald et al., 2022), 0.1–0.57 (Rondón et al., 1993), 0.098 (Breuninger et al., 2013), 0.2–0.5 (Horii et al., 2004), and 0.02–0.64 (Puxbaum and Gregori, 1998) for a mountain observatory surrounded by coniferous trees, boreal coniferous forests, a temperate coniferous forest, a temperate mixed deciduous forest, and a temperate oak forest, respectively, where a combination of soil and foliage deposition was measured. Horii et al. (2004) saw an increase in the deposition velocity with an increasing mixing ratio: from 0.2 at 1 ppbv to 0.5 at 30 ppbv. Puxbaum and Gregori (1998) reported monthly averages of 0.02–0.64 ; however, their average nighttime deposition velocities were below 0.05 . The deposition velocities determined here are a factor of 5–40 higher than what has been measured for nighttime foliage deposition velocities to the leaves of different trees native to California (Delaria et al., 2020, 2018), although they are in good agreement with measurements for daytime. It is, however, important to note that the deposition velocities estimated here are upper limits, as the estimation of the NO emission rate is an upper limit and the chemical loss of is not taken into account. Using an average nighttime mixing ratio of 1650 and 1450 pptv for the Continental and Atlantic phase, respectively, results in respective deposition rates of 13.6 6.8 and 18.7 8.5 , which are in reasonable agreement with that measured for soil deposition in a sessile oak forest (9.67 1.92 ) during the summer (Rosenkranz et al., 2006). The estimated NO soil emission rate and deposition rate are, within the uncertainties, identical, which means that the Rambouillet Forest is not a significant direct source or sink of .
4 ConclusionsMeasurements of NO, , , and during the ACROSS campaign (June–July 2022) in the Rambouillet Forest, southwest of Paris, France, have been used to gain insight into nighttime processes controlling in an anthropogenically impacted forest environment. Based on HYSPLIT back trajectories, two phases of the campaign were identified: one dominated by air originating over the Atlantic Ocean (“Atlantic”), which on average had high relative humidity and low mixing ratios, and one dominated by continental air masses from different urban/industrialized regions (“Continental”), which on average had a lower relative humidity than the Atlantic phase and higher mixing ratios. Strong diel profiles were observed in the measurements across the campaign, with daytime peak mixing ratios varying from to 90 ppbv and nighttime values tending towards 0–10 ppbv. The daily variation was driven by a variable but generally rapid deposition to soil and foliar surfaces, with a strong influence of relative humidity (influencing the surface resistance to uptake) and inversion (influencing the rate of entrainment of from above the canopy).
During the Atlantic phase, periods of sustained NO above the instrumental detection limit were observed at nighttime, when was sufficiently low (i.e. the NO lifetime was sufficiently long). This enabled the derivation of an average NO emission rate from the soil () of , which was confirmed by the approximately linear increase in NO observed in the absence of in the Atlantic phase. The estimated is in broad agreement with previous measurements in other European temperate forests with similar tree types to those found in the Rambouillet Forest. The uncertainty in the estimated NO emission rate is determined from the uncertainties in NO and at 3–5 m above the ground, which leads to higher relative uncertainties at low NO and mixing ratios. Measurements of either NO fluxes or highly resolved height profiles of NO and will improve the NO emission rate estimate during future field campaigns.
An increase in at night would be expected with a constant NO emission rate of in the presence of , as observed in this study; however, this was not the case. The lack of an increase in was used to estimate first-order decay constants of and , resulting in an effective lifetime of of –3 h. The loss of at nighttime is presumably driven by deposition to soil and foliar surfaces, as the lifetime of with respect to its reactions with OH and at night is and 14 h, respectively. By comparison, the daytime lifetime of with respect to loss by reaction with OH is about 1 d. We conclude that the nighttime deposition of is a major sink of boundary layer in this forested environment.
Data availability
All data from the ACROSS campaign, including (10.25326/705, Andersen and Crowley, 2023a; 10.25326/687, Crowley, 2024; 10.25326/708, Dusanter and Jamar, 2023; 10.25326/512, Xue et al., 2023b), (10.25326/707, Crowley, 2023), organic nitrates (10.25326/706, Andersen and Crowley, 2023b), HONO (10.25326/709, Xue et al., 2023a), meteorological quantities (10.25326/437, Denjean, 2023), and peroxy radicals (10.25326/509, Kukui, 2023), can be found at
The supplement related to this article is available online at:
Author contributions
All authors contributed measurements. Data analysis was conducted by STA, with contributions from JNC and PD. CC and VM organized the field campaign, with contributions from the individual group leads. STA and JNC developed the manuscript, with contributions from all co-authors.
Competing interests
At least one of the (co-)authors is a member of the editorial board of Atmospheric Chemistry and Physics. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.
Special issue statement
This article is part of the special issue “Atmospheric Chemistry of the Suburban Forest – multiplatform observational campaign of the chemistry and physics of mixed urban and biogenic emissions”. It is not associated with a conference.
Acknowledgements
Simone T. Andersen is thankful to the Alexander von Humboldt foundation for funding her stay at Max Planck Institute for Chemistry.
Patrick Dewald gratefully acknowledges the Deutsche Forschungsgemeinschaft (“MONOTONS” project; project no. 522970430).
The ACROSS project has received funding from the French National Research Agency (ANR), under the investment programme integrated into France 2030 (grant no. ANR-17-MPGA-0002), and was supported by the French National programme LEFE (Les Enveloppes Fluides et l'Environnement) of the CNRS/INSU (Centre National de la Recherche Scientifique/Institut National des Sciences de l'Univers). Data from the ACROSS campaign are hosted by the French national Atmosphere Data Centre (AERIS).
IMT Nord Europe acknowledges financial support from the CaPPA project, which is funded by the French National Research Agency (ANR) through the PIA (Programme d'Investissement d'Avenir; grant no. ANR-11-LABX-0005-01), the regional council of Hauts-de-France, and the European Regional Development Fund (ERDF).
Financial support
This research has been supported by the Agence Nationale de la Recherche (grant nos. ANR-17-MPGA-0002 and ANR-11-LABX-0005-01), the Alexander von Humboldt-Stiftung (postdoc fellowship), and the Deutsche Forschungsgemeinschaft (grant no. 522970430).The article processing charges for this open-access publication were covered by the Max Planck Society.
Review statement
This paper was edited by Amos Tai and reviewed by two anonymous referees.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2024. This work is published under https://creativecommons.org/licenses/by/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
The budget of reactive nitrogen species, which play a central role in atmospheric chemistry (e.g. in photochemical
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details










1 Atmospheric Chemistry Department, Max Planck Institute for Chemistry, 55128 Mainz, Germany
2 Institut de Combustion, Aérothermique, Réactivité Environnement (ICARE), CNRS, 45071 Orléans, France
3 CNRM, Universite de Toulouse, Meteo-France, CNRS, 31057 Toulouse, France
4 IMT Nord Europe, Institut Mines-Télécom, Université de Lille, Center for Energy and Environment, 59000 Lille, France
5 Univ Paris Est Creteil and Université de Paris Cité, CNRS, LISA, 94010 Créteil, France
6 Université Paris Cité and Univ Paris Est Creteil, CNRS, LISA, 75013 Paris, France
7 Laboratoire de Physique et Chimie de l'Environnement et de l'Espace (LPC2E), CNRS, Orléans, France
8 Institut de Combustion, Aérothermique, Réactivité Environnement (ICARE), CNRS, 45071 Orléans, France; University Mohammed VI Polytechnic (UM6P), Lot 660, Hay Moulay Rachid, Ben Guerir, 43150, Morocco
9 Wolfson Atmospheric Chemistry Laboratory, Department of Chemistry, University of York, York, YO10 5DD, UK