Laboratory measurements of stomatal NO2 deposition to native California trees and the role of forests in the NOx cycle

Both canopy-level field measurements and laboratory studies suggest that uptake of NO2 through the leaf stomata of vegetation is a significant sink of atmospheric NOx. However, the mechanisms of this foliar NO2 uptake and their impact on NOx lifetimes remain incompletely understood. To understand the leaf-level processes affecting ecosystem-scale atmosphere–biosphere NOx exchange, we have conducted laboratory experiments of branch-level NO2 deposition fluxes to six coniferous and four broadleaf native California trees using a branch enclosure system with direct laser-induced fluorescence (LIF) detection of NO2. We report NO2 foliar deposition that demonstrates a large degree of inter-species variability, with maximum observed deposition velocities ranging from 0.15 to 0.51 cms−1 during the daytime, as well as significant stomatal opening during the night. We also find that the contribution of mesophyllic processing to the overall deposition rate of NO2 varies by tree species but has an ultimately inconsequential impact on NOx budgets and lifetimes. Additionally, we find no evidence of any emission of NO2 from leaves, suggesting an effective unidirectional exchange of NOx between the atmosphere and vegetation.


Introduction
Nitrogen oxides (NO x ≡ NO + NO 2 ) are a form of reactive nitrogen that play a major role in the chemistry of the atmosphere. NO x catalyzes tropospheric ozone formation, contributes to the production of photochemical smog, and influences the oxidative capacity of the atmosphere (Crutzen, 1979). NO x is primarily emitted as NO through fossil fuel burning, lighting, and soil microbial activity (Seinfeld and Pandis, 2006). The latter source is of particular importance in remote forested and agricultural regions, where emission from soils is the primary source of NO x (e.g., Jacob and Wofsy, 1990;Lerdau et al., 2000;Seinfeld and Pandis, 2006;Romer et al., 2018;Almaraz et al., 2018).
Understanding the fate of atmospheric NO x , in addition to its emission sources, is essential for interpreting the impact of NO x on atmospheric chemistry. Prior studies have demonstrated that NO 2 can directly deposit to foliage after diffusion through stomata (e.g., Teklemmariam and Sparks, 2006;Chaparro-Suarez et al., 2011;Breuninger et al., 2013;Delaria et al., 2018). The currently understood mechanism of this uptake process is as follows: NO 2 enters through the stomatal cavity and dissolves into the apoplastic fluid, forming nitrate, which then is reduced to ammonium by the enzyme nitrate reductase (Park and Lee, 1988;Ammann et al., 1995;Tischner, 2000;Lillo, 2008;Heidari et al., 2011). There is evidence that NO 2 may also be directly scavenged by antioxidants, most notably ascorbate (Ramge et al., 1993;Teklemmariam and Sparks, 2006). These processes may be impacted by the leaf pH, which is known to change under conditions of limited water availability (Bahrun et al., 2002). Experiments using 15 N as an isotopic tracer have demonstrated that absorbed NO 2 is eventually assimilated into amino acids (Rogers et al., 1979;Okano and Totsuka, 1986). Although the role of stomatal conductance (g s ) in controlling the deposition of NO 2 is well-documented, the impact of mesophyllic processes remains poorly resolved. These mesophyllic mechanisms are complex and include any process taking place between the intercellular air space and the ultimate nitrogen assimilation site. The question of whether and how much mes-ophyllic processes affect NO x budgets at the canopy scale thus persists.
The most divisive example of the mesophyll quandary is the sometimes-reported emission of NO x from plants, mostly in the form of NO, at low NO x mixing ratios that would be relevant to remote forested regions (Johansson, 1987;Rondón and Granat, 1994;Hereid and Monson, 2001;Sparks et al., 2001;Teklemmariam and Sparks, 2006). This would, under many conditions, indicate that trees instead serve as a constant source, rather than sink, of NO x . However, this idea has been called into question by a number of recent studies including Chaparro-Suarez et al. (2011), Breuninger et al. (2013, and Delaria et al. (2018). It is possible that the magnitude and direction of the NO x flux to leaves may vary depending on the species and conditions. One such factor that has been suggested to impact foliar emission and deposition of NO x is elevated soil nitrogen. Soil nitrate fertilization has been documented to lead to an increase in nitrate reductase activity in the needles of Scots pine seedlings (Andrews, 1986;Pietilainen and Lahdesmaki, 1988;Sarjala, 1991). It is possible that as a result of abundant nitrate fertilization, nitrate accumulates in leaves, leading to emission or a reduction in uptake. For example, Chen et al. (2012) observed an increase in NO emission and Teklemmariam and  detected an increase in NO 2 emission under conditions of elevated soil nitrate. By contrast, Joensuu et al. (2014) found no evidence of fertilization-induced NO x emissions. No influence of soil nitrogen on either NO 2 or NO uptake has been documented (Okano and Totsuka, 1986;Teklemmariam and Sparks, 2006;Joensuu et al., 2014).
In this study we present results from laboratory measurements of NO 2 fluxes of 10 native California tree speciessix conifers and four broadleaf trees -using a branch enclosure system and laser-induced fluorescence (LIF) detection of NO 2 . Here we investigate the relative influence of stomatal and mesophyllic processes on the total uptake rate of NO 2 under atmospherically relevant conditions. Our aim is to assess the factors controlling NO 2 foliar deposition and their ultimate impact on the NO x cycle. To test this, we measured the NO 2 deposition velocity over a range of stomatal conductances and considered evidence for additional limits on the uptake rate. We also conducted experiments under drought and elevated soil nitrogen and tested for indications of NO 2 emission or changes in the apparent mesophyllic uptake limit.

Tree specimens
Foliar deposition of NO 2 was investigated in the laboratory using 10 native California tree species -Pinus sabiniana, Pinus ponderosa, Pinus contorta, Pseudotsuga menziesii, Calocedrus decurrens, Sequoia sempervirens, Arbutus men-ziesii, Acer macrophyllum, Quercus agrifolia, and Quercus douglasii. Three to six individuals of each species were purchased from a local native California plant nursery (Native Here Nursery) or Forestfarm, where the plants were grown from seeds and cuttings. The tree specimens were grown in a nutrient-rich commercial soil mixture of Sun Gro Sunshine no. 4 and Supersoil potting soil in 20-40 L pots in an outdoor section of the Oxford facility greenhouse at the University of California, Berkeley. The trees were 2-3 yr old when measurements were taken. No additional fertilizers or pesticides were used on the plants. Trees were transported into the lab for experimentation, where they were exposed to a 12 h light-dark cycle. Trees were illuminated with an LED diode array of 430-475 and 620-670 nm lights (Apollo Horticulture). For the deciduous trees (Q. douglasii and A. macrophyllum) experiments were run between May and September 2019. For all other species experiments were conducted yearround, between October 2018 and November 2019.

LIF measurement of NO 2 deposition fluxes
Measurements were made with a dynamic chamber and LIF detection of NO 2 . A full description of our apparatus can be found in Delaria et al. (2018). Briefly, an NO 2 standard was mixed with humidified zero air (air filtered to remove NO x and reactive species) and delivered to a ∼ 10 L chamber enclosing the branch of a tree at a total flow rate of ∼ 6000 cm 3 min −1 (Fig. 1). The lifetime of air within the chamber was ∼ 2 min. Humidity was adjusted by controlling the fraction of zero air that passed through a bubbler filled with distilled water. The mixing ratios of NO 2 entering the chamber were typically between 0 and 10 ppb. Some of the air entering the chamber was diverted to cell no. 1 of the NO 2 LIF analyzer and to two LI-COR instruments (6262 and 7000), such that the flow rate of air directly into the chamber was ∼ 5000 cm 3 min −1 . Cell no. 1 of the LIF analyzer measured NO 2 and the two LI-COR instruments measured CO 2 and H 2 O in the in-flowing air stream. Air from the chamber was simultaneously pumped out to cell no. 2 of the NO 2 LIF analyzer and the LI-COR-7000 instrument for measuring the mixing ratio of NO 2 within the chamber and the change in CO 2 and water vapor between the incoming and outgoing air streams, respectively (Fig. 1). A slight positive pressure was maintained within the chamber to ensure lab air did not leak into the chamber.
Fluxes of NO 2 to leaves were calculated according to (Eqs. 1 and 2): where [NO 2 ] in and [NO 2 ] out are concentrations of NO 2 entering and exiting the chamber, respectively, at chamber equilibrium. Chamber equilibrium is achieved when the flow rates in and out of the chamber are equal and can be iden- tified by a constant concentration of [NO 2 ] out . [NO 2 ] comp is the compensation point concentration -the ambient concentration of NO 2 at which the tree acts as a source rather than a sink of NO 2 . Q is the flow rate (cm 3 s −1 ), A is the enclosed one-sided leaf area, and V d is the deposition velocity. The leaf area was determined using the ImageJ software package (Schneider et al., 2012), and the flow rate was measured at the beginning of each experimental run (Mesa Laboratories 510-M Bios Defender). Peroxyacetyl nitrate (PAN) and acetone were also delivered to the chamber for simultaneous measurements of PAN stomatal deposition. Negligible thermal production of NO 2 was observed. The results of PAN deposition experiments are discussed in Place et al. (2020). The NO 2 mixing ratio was also corrected for the differences in collisional quenching of the excited state NO 2 by water vapor in cells no. 1 and no. 2, caused by transpiration of the tree within the chamber (Thornton et al., 2000).
where X H 2 O is the difference in the water vapor mole fraction between the chamber and the incoming air stream. Experiments to an empty chamber were conducted approximately every 2 months during this study to calculate the deposition of NO 2 to the chamber walls. The wall loss was at maximum ∼ 2 % of the [NO 2 ] in concentration and was background-subtracted from our flux calculations. Deposition velocities were determined using the method described in Delaria et al. (2018): a weighted orthogonal distance linear regression was performed on NO 2 fluxes (determined using Eq. 1) against [NO 2 ] out to obtain a slope equal to V d . A positive x intercept was interpreted as evidence for a possible compensation point. During each day of experimentation we stepped through at least eight different NO 2 concentrations, with each concentration step lasting for 40 min. Uncertainty in V d was obtained through propagating uncertainty in measured NO 2 concentrations, Q, and A. The uncertainty in NO 2 concentrations was estimated as 1 standard deviation of variation in measurements during the last 10 minutes of each concentration step. The uncertainty in Q was estimated as < 1 % and a 10 % uncertainty was estimated for the enclosed one-sided leaf area.
The deposition velocities measured can be related to the resistance-model framework for deposition of trace gases developed by Baldocchi et al. (1987) (Eqs. 4-6).
R is the total resistance to deposition, R a is the aerodynamic resistance, R b is the boundary layer resistance, and R leaf is resistance to uptake by the leaf. R a was assumed to be negligible under our chamber conditions (Pape et al., 2009;Breuninger et al., 2012;Delaria et al., 2018). R leaf is made up of R cut , R s , and R m . Respectively, these refer to the cuticular resistance (resistance to deposition to the surface of the leaf), stomatal resistance (1/g s ), and mesophyllic resistance (resistance associated with all processes taking place within the leaf that limit uptake).

Measurement of stomatal conductance
CO 2 and water vapor exchanges were measured using the LI-COR-6262 and LI-COR-7000 instruments. Measurements of water vapor exchange were used to calculate the transpiration rate (E) and total conductance to water vapor (g w t ) using Eqs. (7) and (8), according to von Caemmerer and Farquhar (1981).
where w a and w e are the mole fractions of water vapor to the outgoing and incoming airstreams, respectively, and ω i is the internal leaf water vapor mole fraction. ω e was measured with the LI-COR-6262 with dry air as a reference, and ω (ω a − ω e ) was measured with the LI-COR-7000 with incoming air as the reference. ω e was kept constant throughout a day of measurements and was varied between days. Measurements of an empty chamber were also used to calculate and correct for the water vapor deposition to the chamber at varying relative humidity. The difference between ω a and ω e for an empty chamber was not statistically significant and at E. R. Delaria et al.: Stomatal NO 2 deposition to native California trees all relative humidity levels was within instrumental uncertainty of the LI-COR-6262. ω i was assumed to be the saturation vapor pressure at the leaf temperature, which was measured with a thermocouple at the surface of an enclosed leaf. The chamber temperature was measured with a second thermocouple and was typically 20 ± 3 • C. Over the course of a day the temperature and humidity varied by a maximum of 2 • C and 5 %, respectively. These deviations were not found to be significantly correlated with stomatal opening. The photosynthetic photon flux density (PPFD) was monitored outside the chamber with a LI-COR quantum sensor (LI-COR LI-190SA) and was 1190 µmol m −2 s −1 , approximately the PPFD for Berkeley, California, at noon during the month of October. We performed calculations based on von Caemmerer and Farquhar (1981) to confirm this is above the photon flux required to achieve maximal stomatal aperture for tree types relevant to this study. Total conductance was calculated as the average over the light or dark period of an experiment. The uncertainty in our calculation of total conductance to water vapor was primarily influenced by uncertainty in the leaf temperature and the assumption of leaf water vapor saturation. We observed fluctuations in the temperature of enclosed leaves of ±2 • C. Total uncertainty in g w t was determined by propagating this uncertainty in leaf temperature, which resulted in larger estimated uncertainties at larger chamber humidities, usually coinciding with higher stomatal conductances. Chamber relative humidity was maintained at less than 90 % to minimize this effect. Variations in stomatal conductance were achieved by varying the mole fraction of water vapor in the air delivered to the chamber. The LI-COR-6262 instrument was calibrated weekly using standard CO 2 cylinders and a LI-COR-610 dew point generator. The LI-COR-7000 instrument was calibrated daily.
The stomatal conductance (g w s ) could then be calculated from Eq. (9): where 1/g w b is the boundary layer resistance to water vapor. The boundary layer resistance to water vapor was estimated to be negligible under our experimental conditions, with an upper bound of 0.6 s cm −1 . This was calculated by measuring the deposition of NO 2 to a 30 cm 2 tray of activated charcoal and confirmed by measuring the evaporation from a watersoaked Whatman no. 1 filter paper. (Delaria et al., 2018). A detailed description of our assumption of negligible R b can be found in Sect. 3.1. Stomatal (g s ) and total (g t ) conductances to NO 2 were calculated by scaling the values for water vapor by the ratio of diffusivities in air (D NO 2 /D H 2 O ) according to Massman (1998).

Nitrogen measurements
To test the influence of excess soil nitrogen on the ability of trees to take up nitrogen through their stomata in the form of NO 2 , we fertilized three individuals of both Quercus agrifolia and Pseudotsuga menziesii with a 20 mM ammonium nitrate solution. The trees were watered with 250 mL of this ammonium nitrate solution 3 d per week. Three individuals of each species were watered with deionized (DI) water as the control group. The trees underwent this fertilization treatment for 120 d before beginning dynamic chamber measurements on NO 2 foliar deposition. NO 2 deposition experiments were conducted for 70 d, during which time the soil fertilization treatments were continued.

Soil nitrogen
Approximately 5 mg of a soil core sample was taken each day from the individual on which we conducted an NO 2 deposition experiment. The soil was sifted through a mesh 2 mm sieve. Soil nitrate and ammonium were extracted by shaking ≈ 2.5 mg of the soil sample in 30 mL of ≈ 2 M KCl for 1 h, followed by filtering the samples through a Whatman no. 1 filter paper. The other ≈ 2.5 mg was dried in a drying oven at 60 • C for at least 48 h. The mass of the soil after drying was measured to determine the percentage dry mass of the extracted soil sample. Twelve quality control (QC) samples were made and carried through the extraction process. These samples consisted of six 30 mL KCl blanks and six KCl samples spiked with either 5 mL (low QC) or 10 mL (high QC) of 1 ppm ammonium nitrate solution. NH + 4 and NO − 3 were measured using a colorimetric synthesis following the method of Sims et al. (1995) and Decina et al. (2017). Briefly, a standard 1 ppm stock solution of ammonium nitrate was made from ammonium nitrate solid dissolved in Milli-Q water and was diluted to 0, 0.1, 0.2, 0.3, 0.4, and 0.5 mg L −1 in 1 cm, 2.5 mL cuvettes. These standard solutions served as the calibration standards; we made three sets of calibration standards for both ammonium and nitrate analysis. All glassware was acid washed in a 1 M solution of HCl prior to all measurements and extractions to prevent contamination.
For ammonium analysis, 160 µL of each soil extraction sample from the control group, 10 µL from the fertilizertreated group, and 1.6 mL of the QC samples were pipetted into individual cuvettes. In total, 100 µL of 0.2 M citrate, 200 µL of 5 mM nitroprusside, 100 µL of 0.3 M hypochlorite reagents, and 500 µL of Milli-Q water were then added sequentially into each cuvette. The cuvettes were filled to a final volume of 2.5 mL with KCl, and the samples were allowed to sit for 30 min. For nitrate measurements, 320 and 10 µL of soil samples from the control and fertilized groups, respectively, and 1550 µL of the QC samples, were pipetted into separate cuvettes. In total, 950 µL of a regent containing 1 g L −1 vanadium chloride and 25 mg L −1 N-(1-naphthyl)ethylenediamine (NEDD) were subsequently added to each cuvette; these were then filled to a final volume of 2.5 mL with KCl and allowed to sit for 24 h. To test the effects of the soil matrix on the calibration, 160 and 320 µL of a control Q. agrifolia soil extraction sample were added to one set of calibration standards for ammonium and nitrate analysis, respectively.
Concentrations of ammonium and nitrate in each sample were determined with colorimetric measurements using a custom-built spectrophotometer. The spectrophotometer light source was a broad spectrum quartz tungsten halogen lamp (QTH10 Thorlabs Inc.). The absorption of each sample and standard was measured with the light source passing through a 540 ± 2 nm bandpass filter (FB570-10 Thorlabs Inc.) for nitrate analysis or a 670 ± 2 nm bandpass filter (FB540-10 Thorlabs Inc.) for ammonium analysis.

Uncertainty analysis
Concentrations of ammonium and nitrate in the soil extraction samples were determined from the slope in their respective calibration curves. The calibrations for ammonium and nitrate analysis had respective uncertainties of 7 % and 5 %. The slopes of the calibration curves with added sample from a Q. agrifolia soil extraction were not statistically different from those containing only standards, allowing us to exclude the possibility of interference from the soil matrix.
The accuracy uncertainty in the high and low QC samples were 3 % and 11 %, respectively, for ammonium measurements, and 3 % and 12 % for nitrate measurements. We estimated the resulting uncertainty for cuvette samples with less than 0.15 mg L −1 NH + 4 or NO − 3 (≈ 1.8 µg mg −1 soil NH + 4 or NO − 3 ) to be 15 %. Samples with larger concentrations were estimated to have 5 % uncertainty. The blank quality control standards contained 0.04 mg L −1 ammonium and nitrate. This was blank-subtracted from each sample.

Leaf nitrogen
After deposition experiments were completed the leaves were removed from the trees and dried for 48 h in a drying oven. The leaves were then ground to a fine powder and the percent nitrogen, hydrogen, and carbon content were measured with an inductively coupled plasma-optical emission spectrometry instrument (Optima 7000 DV ICP-OES)

Drought stress
Calocedrus decurrens and Pinus ponderosa were drought stressed to study the impact of drought on NO 2 deposition. Three individuals of each species were watered daily (control group) and three individuals of each species were watered with 250 mL once every 4 weeks (drought group). Limitedwater treatment of the drought group was carried out for 60 d before conducting dynamic chamber experiments for NO 2 foliar deposition. NO 2 deposition experiments were run for 30-40 d. During the experiments, the control group was watered 50 mL daily and the experimental group was watered 50 mL once every 2 weeks. The P. ponderosa drought stress experiments took place between March and June 2019. The C. decurrens drought stress spanned from August to December 2019. The xylem water potential ( p ) of the trees was monitored to measure the drought stress level of the trees, using a Scholander pressure chamber (Model 670 PMS Instr. Comp.). Leaves were cut, wrapped in aluminum foil, and then inserted into the pressure bomb. The p of cuttings were measured around 11:00 LT each day. A p measurement lower than −1.0 MPa indicated signs of drought stress in P. ponderosa. C. decurrens did not show evidence of drought stress in p measurements, while in the greenhouse, however, early signs of embolism were observed.

Results
V d was calculated for each day of measurements with a weighted linear regression of measured fluxes and chamber NO 2 concentrations (Delaria et al., 2018). No statistically significant compensation point was observed under any experimental condition for the majority of the species studied, in agreement with previous work (Chaparro-Suarez et al., 2011;Breuninger et al., 2013;Delaria et al., 2018). Only P. menziesii was found to have a compensation point, estimated to be 20 ppt, but this concentration is below the limit of quantification for our instrument, so we consider this measurement to be consistent with a compensation point of zero. V d and g s measurements allowed for consideration of whether the deposition of NO 2 is exclusively stomatally controlled or is also affected by the internal processing in the mesophyll. We rarely observed total closing of the stomata when the chamber lights were turned off at night. All of the deposition observed at night could be explained by deposition to these partially open stomata. This is consistent with previous studies observing only partial closing of stomata at night in a variety of plant species Drake et al., 2013). The results of experiments are shown in (Table 2).

Measurements of mesophyllic resistance
We utilized two methods for analyzing the importance of the mesophyllic resistance on the deposition of NO 2 . Figure 2 shows the predicted stomatally limited NO 2 deposition fluxes, assuming negligible R b , R c , and R m (Flux = b Statistically significant (α = 0.01) compensation point. Compensation point listed is at limit of detection for the instrument. All other compensation points are not statistically significant (α = 0.05). c Statistically significant (α = 0.05) correlation. Correlations not indicated are not statistically significant (α = 0.05). d Maximum deposition velocity that was observed during our experiments and the error associated with that measurement. e Listed maximum g w t the maximum stomatal conductance to water vapor that was observed during our experiments and the error associated with that measurement. Units in mmol m −2 s −1 for ease of comparison with other stomatal conductance studies. f Total slope of measured vs. predicted fluxes (Fig. 2). g Individual slopes of predicted vs. measured fluxes from each day an experiment was run. h Calculated including data in helium. g t [NO 2 ] out ) plotted vs. the measured NO 2 fluxes. Our upper bound measurement of R b for NO 2 was 1 s cm −1 (0.6 s cm −1 for water vapor). Assuming g s = g t would lead to a maximum of a 60 % or 10 % error in the calculated g s with a g t = 0.6 cm s −1 or g t = 0.1 cm s −1 , respectively. However, R b decreases with the enclosed leaf area according to Pape et al. (2009), which at a minimum was 200 cm 2 . The maximum R b in the chamber should have thus been ≈ 0.1 s cm −1 . Assuming g s = g t would lead to a maximum of a 6 % error at g t = 0.6 cm s −1 in this case. Any deviation from unity in the observed slope of predicted vs. measured fluxes can thus be attributed to R m . Any error in our assumption of negligible R b may partially mask the effect of R m . We do not expect that variation in R b due to changes in leaf morphology, micrometeorology, and leaf movement would substantially change the effect of R b , although we cannot rule out the possibility that this was partially responsible for day-to-day fluctuations in NO 2 fluxes. We confirmed the validity of our assumption of negligible R b by comparing measurements of total conductance to water vapor, g w t , in the chamber to measurements of stomatal conductance for the enclosed branch with a LI-COR-6800 instrument under identical environmental conditions of light irradiation, humidity, and temperature. This test was performed on one individual of three different tree species, and in all cases the chamber g w t measurements were found to be approximately equal to the LI-COR-6800 measurements of g w s within the range of uncertainty in g w t . Significant deviations from unity in the slope of g t [NO 2 ] out vs. measured fluxes could be seen in several species, most notably S. sempervirens (Table 2 and Fig. 2). Figure 2 shows each flux measurement as a single data point. For each day of experiments a slope of predicted vs. measured fluxes was obtained from a least squares cu-bic weighted fit on the 8-12 fluxes measured at varying NO 2 concentrations. The reported slope for a given species (Table 2, shown in blue in Fig. 2) was calculated using a weighted average of the slopes from all experiment days. This was done to minimize the contribution of systematic errors potentially introduced by the LI-COR-7000 instrument, which was calibrated daily. All data points for a given day were excluded (shown in red in Fig. 2) if the calculated slope on that day was determined to be an outlier by a generalized extreme studentized deviate test for outliers. Identified outliers were excluded both to account for potentially erroneous deviations in the V d /g t ratio (most likely due to systematic error in calibration of the LI-COR-7000 instrument) and to avoid overweighting of days with abnormally large stomatal conductances. These latter instances normally coincided with low V d /g t ratios, and if these data were also subject to some systematic error, they would bias our analysis of R m .
R m was also explicitly calculated using the relationship of V d and g t . Figure 3 shows V d from each day of experiments plotted against the measured g t . Positive y intercepts are indications of cuticular deposition, and curvatures in the fit away from the 1 : 1 line are implications of mesophyllic resistance. R m was calculated with a weighted fit of the resistance model: No evidence of cuticular deposition was observed, so only results of R m are recorded ( Table 2). The deposition observed with the chamber lights turned off could be explained completely by the measured stomatal conductance. Fits of the resistance model (Eq. 10) typically resulted in cuticular resistances on the order of 1000 s cm −1 . R m was calculated both Figure 2. Measured fluxes (mflux) plotted against stomatally limited predicted fluxes (pflux = g t [NO 2 ] out ). Drought data and nitrogen fertilization data are included. Blue solid lines are the linear fit to data. Red lines are the 1 : 1 line. Error bars for the measured fluxes are calculated by propagating uncertainty in the measured NO 2 mixing ratios, the flow rate, and the leaf area (Eq. 1). Error bars for the predicted fluxes are calculated by propagating uncertainties in the measured NO 2 mixing ratios and the total conductance (Eq. 8). Red markers indicate data determined to be outliers by a generalized extreme studentized deviate test for outliers. Figure 3. Deposition velocities (V d ) plotted against measured total conductances to NO 2 (g t ). Black markers represent measurements in zero air, and red-yellow markers are measurements in helium. Measurements in helium are subject to less uncertainty introduced by potential systematic error in the leaf temperature. Solid blue lines are the 1 : 1 line, and dashed blue lines are error-weighted fits to the resistance model using only measurements in zero air, assuming the boundary layer resistance is negligible (Eq. 4). Fits to the resistance model including data from helium measurements are shown as dashed red lines.
assuming negligible R b (g s = g t ) and R b = 1 s cm −1 . There were no significant differences between these two calculations (Table 2).

Effects of excess soil nitrogen
The impact of soil fertilization on the foliar uptake of NO 2 by two tree species, Q. agrifolia and P. menziesii, was examined by watering a control group of both species with deionized water and a fertilized group with 20 ppm ammonium nitrate. On average, the soil nitrogen concentrations of data. The V d /g t ratio is plotted against the leaf nitrogen : carbon ratio for (b) Q. agrifolia and (d) P. menziesii. V d /g t ratios less that 1 imply contributions from the mesophyll to the NO 2 uptake rate. On each panel the Pearson's correlation coefficient and the p value for the slope are shown. The amount of soil and leaf nitrogen has no significant impact on the V d /g t ratio.
NH + 4 and NO − 3 were 100× larger for the fertilized groups than the control groups (Table 1). The percentage of leaf nitrogen content approximately doubled between the control groups and the fertilized groups (Table 1).
The effect of soil nitrogen fertilization and leaf nitrogen content on the ratio of V d /g t is shown in Fig. 4. No significant relationship (α = 0.01) was observed for either Q. agrifolia or P. menziesii, suggesting the mesophyllic processing of NO 2 is unaffected by soil or leaf nitrogen content . We also observe no increase in the compensation point of NO 2 as a result of higher leaf nitrogen content or elevated soil nitrogen (Fig. 5).

Drought stress measurements
The impact of drought stress on NO 2 foliar uptake for C. decurrens and P. ponderosa was observed by regularly watering a control group and watering an experimental drought group at much lower frequency (once every 4 weeks in the greenhouse and once every 2 weeks in lab). The median p measured was lower for the drought groups than the control groups (Table 3). C. decurrens drought median p was −0.80 MPa compared to control median of −0.30 MPa, and P. ponderosa drought median was −1.05 MPa compared to control median of −0.60 MPa. The first quartiles of the con- trol groups and third quartiles of the drought groups did not overlap, reflecting a significant difference between the p measurements of the two groups. We also observed a strong correlation between measured p and stomatal conductance. We found a more substantial impact of drought on the water potentials and of the water potentials on the stomatal conductance, in P. ponderosa trees than C. decurrens. Both these California conifer species are quite drought resistant (Pharis, 1966;Kolb and Robberecht, 1996;Maherali and DeLucia, 2000), but these results may indicate C. decurrens is particularly protected against water loss.
The mesophyllic resistance (R m ) calculated showed a statistically significant difference for both C. decurrens and P. ponderosa between drought-stressed and control groups. R m in drought-stressed C. decurrens increased from 0.37 to 1.17 s cm −1 , while in P. ponderosa R m decreased from 0.86 to 0 s cm −1 (Fig. S5 in the Supplement).

Effects of mesophyll resistance on the lifetime of NO x
The mesophyllic resistances (R m ) for each of the 10 tree species measured are calculated from Fig. 3 and Eq. (10) and are tabulated in Table 2, assuming either g s = g t or the upper bound for R b . The slopes of predicted fluxes vs. measured fluxes, calculated in Fig. 2, are also tabulated in Table 2. The importance of the mesophyllic resistance and internal processing of NO 2 can be evaluated by examining both R m and the slope of measured vs. predicted fluxes. We also examined the potential impact of the mesophyllic processing of NO 2 by considering the Pearson's correlation coefficient between g t and the slope for an individual experiment (1 d of light or dark data) of measured vs. predicted fluxes (Fig. S3 in the Supplement). These correlation coefficients can be found in Table 2. The more negative this correlation, the greater the deviation in the slope from unity for higher values of g t , consistent with a larger impact of the mesophyll on the NO 2 uptake rate. All tree species except for C. decurrens, Q. agrifolia, and Q. douglasii show statistically significant correlations (α = 0.05) ( Table 2). R m becomes more important at larger stomatal conductances (lower stomatal resistances), as can be seen with the increasing deviations from 1 : 1 in some species at higher values of g t in Fig. 3. Thus, even for trees with higher calculated R m , the impact of mesophyllic processing is unlikely to be large if the maximum stomatal conductance observed is relatively small, resulting in a slope in the measured vs. predicted flux that does not deviate greatly from unity. This is the case for Q. agrifolia and P. ponderosa. Alternatively, P. sabiniana demonstrates a case of a relatively small R m but also a smaller slope in measured vs. predicted fluxes, driven by consistently larger stomatal conductances (lower R s ) (Fig. 3). The most sizable impacts of mesophyllic NO 2 processing ares seen in S. sempervirens, P. sabiniana, and A. macrophyllum. These species have the largest maximum observed g t (Fig. 3, Table 2) and slopes of measured vs. predicted fluxes of 0.79±0.04, 0.84±0.03, and 0.84±0.03, respectively. However, the greater uncertainty in measurements of stomatal conductance at a larger chamber humidity calls in to question the accuracy of many g t measurements larger than approximately 0.4 cm s −1 .
To evaluate with greater certainty the relationship of V d and g t , we conducted a set of experiments in helium to raise the stomatal conductance by increasing the gas diffusivities while maintaining relatively lower chamber humidity. These experiments were conducted on four of the tree species: P. sabiniana, S. sempervirens, Q. agrifolia, and A. menziesii. In these experiments the V d /g t ratio for P. sabiniana remained close to 1 : 1 up to 1.3 cm s −1 stomatal conductance (Fig. 3). Experiments in helium for this species thus suggest a smaller contribution of the mesophyll (red dashed line in Fig. 3). R m calculated including helium experiments was not statistically different for S. sempervirens, Q. agrifolia, or A. menziesii.
Our laboratory measurements of mesophyllic resistance address the uncertainty in the literature for whether reactions in the mesophyll may be consequential for NO 2 deposition velocities. To our knowledge, no previous studies have explicitly calculated the mesophyllic resistance. Differences between leaf-level deposition velocities and stomatal conductances measured by Breuninger et al. (2013) and ob- servations of leaf ascorbate impacts on uptake rates by Teklemmariam and  have indicated mesophyllic reactions may be important. Additional studies (Gut et al., 2002;Eller and Sparks, 2006;Chaparro-Suarez et al., 2011) have also shown some evidence that between 20 % and 40 % of NO 2 deposition is under mesophyllic control. Our findings, however, suggest nearly 90 % of uptake is controlled by the stomata. Currently, atmospheric models incorporate a mesophyllic resistance to NO 2 of 0.1 s cm −1 (Zhang et al., 2002). This would result in a slope of measured vs. predicted fluxes of 0.94, even with a relatively large average g t of 0.6 cm s −1 . The median slope measured in our study was 0.89. Using the multibox canopy model presented in Delaria and Cohen (2020), we investigated whether our results could possibly imply a more important impact of the mesophyllic resistance on the atmospheric fate of NO x at the canopy level. This model takes into account in-canopy processes (e.g., vertical transport, chemistry) to scale leaf-level processes to the canopy level. The model was run using meteorological conditions for June measured during the BEARPEX-2009 campaign, located at a ponderosa pine forest in the western foothills of the Sierra Nevada (38 • 58 42.9 N, 120 • 57 57.9 W; elevation 1315 m). The model was initialized over 2 d and data from the third day were analyzed. We conducted two model runs at a stomatal conductance (g s ) to NO 2 deposition of 0.3 cm s −1 : the median measured maximum stomatal conductance excluding P. sabiniana -with an R m of either 0.1 or 0.6 s cm −1 ; the median measured R m excluding P. sabiniana. For a stomatal conductance to NO 2 of 0.3 cm s −1 (≈ 0.5 cm s −1 to water vapor) the model predicts only a 2.5 % decrease in NO x lost to deposition with an R m of 0.6 compared with an R m of 0.1 s cm −1 . The lifetime to deposition with an R m of 0.1 and 0.6 s cm −1 was 30.5 and 32.2 h, respectively, representing only a 6 % difference. The total atmospheric lifetime of NO x in the boundary layer with an R m of 0.1 and 0.6 s cm −1 was 4.86 and 4.89 h, respectively, representing only a 0.6 % difference. Even the observed seemingly substantial mesophyllic resistance of S. sempervirens is therefore likely to be irrelevant at the canopy scale. Con-tributions from mesophyllic processing, though mechanistically important at a cellular level, are likely to not matter at the canopy scale in California forests. We therefore suggest that on canopy and regional scales, mesophyllic processes within leaves of trees represent a negligible contribution to NO x budgets and lifetimes in California. More studies on crops, grasses, and North American tree species from outside of California are needed.

Effects of excess soil nitrogen
We observed no effects of soil nitrogen, in the form of NH + 4 and NO − 3 , or the leaf nitrogen content on the ratio of V d /g t (Fig. 4) for either Q. agrifolia or P. menziesii. Changes in this ratio would indicate an effect on the mesophyllic resistance. We did observe declines in g t in the fertilized group relative to the control group during the later stages of experimentation, which coincided with observable evidence of plant stress (e.g., browning, wilting, and beginning signs of embolism). All variation in the uptake rates (V d ) could be explained exclusively with deviations in g t . These results are supported by previous studies which have also found a negligible impact of nitrogen fertilization on NO 2 uptake (Teklemmariam and Joensuu et al., 2014). If the fertilizer results in increased NO − 3 and NO − 2 in the leaf, this suggests that the mechanism of NO 2 uptake via dissolution and subsequent reduction in NO − 3 and NO − 2 is likely not reversible and not influenced by accumulation of NO − 3 and NO − 2 within the mesophyll. Alternatively, if the increase in soil nitrogen leads only to an accumulation of organic nitrogen in the leaf, this increase has no effect on the uptake rates. Numerous studies indicate nitrate reductase activity is affected by the presence of ammonium, nitrate, and organic nitrogen in the form of amino acids in a variety of plant species (e.g., Datta et al., 1981;McCarty and Bremner, 1992;Woodin et al., 2006). Based on our current understanding of the mechanism of NO 2 mesophyllic processing, if reactions in the mesophyll indeed affect the rate of stomatal uptake, our fertilization experiments should have succeeded in changing NO 2 uptake rates, given that they succeeded in changing leaf nitrogen content. Because we observed no ef-fect of nitrogen fertilization on NO 2 uptake, we believe that this finding further supports that reactions within the mesophyll may be atmospherically unimportant. It is also possible that the disproportionation of NO 2 to form nitrate and nitrite and scavenging by antioxidants (e.g., ascorbate) are the ratelimiting steps in the mesophyllic processing of NO 2 rather than enzyme activity. More leaf and cellular-level studies are needed to elucidate the uptake mechanism.
We also did not observe any evidence for a relationship between the NO 2 compensation point and the soil nitrogen content nor the leaf nitrogen content (Fig. 5) for either Q. agrifolia or P. menziesii. In general, we only observed uptake and no emission of NO 2 . We also conducted measurements of NO uptake and emission, but the fluxes measured were so small they were below the limit of quantification for our instrument. Chen et al. (2012) observed a strong relationship between NO emissions from stomata and soil nitrate fertilization. However, the maximum NO emissions they measured were a factor of 50 lower than the deposition of NO 2 measured here. NO emission from leaves is therefore not likely to be an important source of atmospheric NO x . P. menziesii was the only tree examined in our experiments that demonstrated any evidence for emission of NO 2 at low mixing ratios, with a compensation point of ≈ 20 ppt. This concentration is much lower than has been observed in previous studies that have detected an NO 2 compensation point (Hereid and Monson, 2001;Teklemmariam and Sparks, 2006). However, this concentration is near the limit of detection for our instrument (Delaria et al., 2018) so should be taken with a grain of salt. A possible cause for discrepancy between our study and those that have measured significant NO 2 compensation points is that our experiments are conducted only using photosynthetically active radiation. Some past work has demonstrated that UV light may cause photolysis of nitrate at the leaf surface and subsequent emission of NO x (Hari et al., 2003;Raivonen et al., 2006). The lack of a relationship between NO x emission and soil N fertilization contrasts with the results of Teklemmariam and Sparks (2006) but is consistent with the nitrogen fertilization experiments conducted by Joensuu et al. (2014).

Effects of drought stress
Although there was a statistically significant impact of drought stress on R m , this is unlikely to be important to the overall uptake rates of NO 2 at the canopy scale for reasons discussed in Sect. 4.1. The differing effects of drought on R m between P. ponderosa and C. decurrens is surprising, with the drought group having a smaller R m in P. ponderosa and larger R m in C. decurrens. However, in the case of P. ponderosa, the lack of measurements at larger g t is likely to mask any existing mesophyllic effects, leading to minimal deviation in the total slope of predicted vs. measured fluxes from unity (Fig. S5). Despite a calculation of significant mesophyllic resistance in both drought and control C. decurrens individuals, the lack of a statistically significant (α = 0.05) correlation between g t and the slopes of predicted vs. measured fluxes casts doubt on this relationship. The control group of P. ponderosa is the only one for which this correlation is significant. The impact of drought on NO 2 uptake at the leaf level is thus primarily its effect on the stomatal conductance. At the canopy level, documented effects of drought on leaf area also require consideration (Pharis, 1966;Kolb and Robberecht, 1996;Maherali and DeLucia, 2000).

Effects of nighttime stomatal deposition
Most atmospheric chemical transport models, such as the abundantly utilized WRF-Chem and GEOS-Chem, use the Wesley model for parameterizing dry deposition of gaseous species (e.g., Skamarock et al., 2008;Fast et al., 2014;Amnuaylojaroen et al., 2014;Ng et al., 2017). The Wesley model implicitly assumes the stomata are fully closed at night, despite more recent studies demonstrating many species of vegetation maintain partially open stomata at night (Musselman and Minnick, 2000;Dawson et al., 2007;Fisher et al., 2007;Drake et al., 2013). We find minimal cuticular deposition of NO 2 , in agreement with several other studies (Sparks et al., 2001;Chaparro-Suarez et al., 2011). However, field observations have shown that substantial leaf-level nighttime deposition of NO 2 is necessary to explain nighttime levels of NO x (Jacob and Wofsy, 1990). The same phenomenon has been seen with other gaseous molecules, most notably PAN, which has also been suggested by a number of field observations to have significant non-stomatal deposition at night (Turnipseed et al., 2006;Wolfe et al., 2009;Crowley et al., 2018). Sparks et al. (2003) did not observe any evidence of non-stomatal deposition in the laboratory, but more recently Sun et al. (2016) implicated non-stomatal deposition in accounting for over 20 % of PAN leaf-level deposition. Our PAN deposition experiments, however, discussed in Place et al. (2020), also did not identify any significant non-stomatal deposition. Despite the existing differences regarding the importance of non-stomatal PAN deposition, we suggest that a significant portion of the "missing" deposition sink of NO 2 and peroxyacyl nitrates at night may be due to non-total closure of the stomata.
To assess the impact of nighttime stomatal opening on the atmospheric fates and lifetimes of NO x at night, we ran our 1-D multibox canopy model, under the conditions described above, at the minimum, maximum, 25th percentile, and 75th percentile of the median nighttime deposition velocities measured in this study (0.004, 0.087, 0.009, and 0.038 cm s −1 , respectively). At such low stomatal conductances, we found these deposition velocities to be not significantly different (α = 0.05) from the stomatal conductance to NO 2 . The fractions of NO x loss to deposition and chemistry to these levels of stomatal opening at night are shown in Fig. 6. Here chemistry represents loss to HNO 3 , RONO 2 , and PAN, and nighttime is defined from 20:00 to 05:00. The range between the Figure 6. Fraction of NO x loss to deposition and chemistry (nitric acid, alkyl nitrate, and peroxyacyl nitrate) at night (20:00-05:00). The four dashed lines between the deposition and chemistry fractions show NO x loss with a nighttime NO 2 deposition velocity of 0.004, 0.009, 0.038, and 0.087 cm s −1 . These deposition velocities, respectively, represent the minimum, first quartile, third quartile, and maximum of the median nighttime deposition velocities measured for the native California trees examined in this study.
first and third quartile of the nighttime deposition observed results in a range in the fraction of NO x loss to deposition from 13 % to 25 % (Fig. 6) and a range in total NO x lifetime from ≈ 7.5 to 5 h.
The relatively large impact of the nighttime stomatal conductance on the fate of NO x , coupled with the large degree of inter-species variation in nighttime stomatal opening, indicates a need for more extensive studies of the nighttime deposition of NO 2 . Deposition is a permanent sink of atmospheric NO x , contrasting with the chemical nighttime sink of NO x to peroxyacyl nitrates (Russell et al., 1986;Cantrell et al., 1986;Perring et al., 2009). Heterogenous reactions at aerosol surfaces involving the NO x reservoir N 2 O 5 and alkyl nitrate formation are among the other major nighttime chemical NO x sinks (Perring et al., 2009;Stavrakou et al., 2013;Kenagy et al., 2018). The relative fractions of nighttime NO x loss to deposition and chemistry is likely to have a substantial impact on the fate of atmospheric NO x and the cycling of reactive nitrogen.

Impacts on the nitrogen cycle in California
To our knowledge, this is the first study conducted on NO 2 stomatal deposition to native California tree species, except for Q. agrifolia (Delaria et al., 2018). However, there are many measurements of the stomatal conductance of California trees (Table 4) with which to compare our maximum total conductance to water vapor measurements (max g w t ). Murray et al. (2019) examined patterns in maximum stomatal conductance to water vapor (max g w s ) across bioclimatic zones. Among the species they looked at were A. menziesii, A. macrophyllum, and Q. agrifolia, for which they measured an average max g w s of 550, 420, and 390 mmol m −2 s −1 , respectively. In comparison, our measurements of max g w t for these species were, respectively, 210 ± 10, 400 ± 100, and 90 ± 20 mmol m −2 s −1 . Our estimates of max g w t for A. menziesii and Q. agrifolia are substantially lower. Matzner et al. (2003) report larger conductances than we do for Q. douglasii as well (Table 4). Maire et al. (2015) determined a maximum stomatal conductance for A. menziesii of 150 mmol m −2 s −1 , in better agreement with our measurements. Henry et al. (2019) measured a similar maximum stomatal conductance of Q. agrifolia to our study of 95 mmol m −2 s −1 , also in better agreement with our results than Murray et al. (2019). Maire et al. (2015) measured a maximum stomatal conductance to water vapor for P. ponderosa and S. sempervirens of 124 and ∼ 91 mmol m −2 s −1 , respectively -considerably smaller than the values measured in this study. Ambrose et al. (2010) measured a max g w s for S. sempervirens of 240 m −2 s −1 , in better agreement with our measurements. C. decurrens max g w t reported here is in good agreement with previous measurements of max g w s (Grantz et al., 2019). For Quercus and Acer species in similar climate regions to California, Maire et al. (2015) calculated max g w s ranging from 103 to 890 and 112 to 320 mmol m −2 s −1 , respectively. The median of max g w t for all four angiosperms we measured was 200 mmol m −2 s −1 , in good agreement with the 250 mmol m −2 s −1 median of all angiosperms in Mediterranean climate regions found by Murray et al. (2019) and the 215 mmol m −2 s −1 median found by Maire et al. (2015). Our median for the six gymnosperms measured was 230 mmol m −2 s −1 , considerably larger than the median 100 mmol m −2 s −1 max g w s found by Maire et al. (2015) in Mediterranean climate regions (defined as warm temperature steppe regions as classified by Kottek et al., 2006).
Overall, the total conductances to water vapor measured in our laboratory experiments fall within the ranges of maximum stomatal conductances measured in previous studiesalthough inconsistencies exist in the current literature. (We also consider this to further support our conclusion that the boundary layer resistance in our chamber is negligible.) Possible discrepancies may have resulted from the location at which each species was measured, growing conditions, ages of the trees, etc. Nevertheless, our NO 2 deposition resultsand their applicability to California forests -are bolstered by the fact that our max g w t measurements fall within the ranges of max g w s measured for mature trees in the field. To assess the impact of the lab-measured deposition velocities on the NO x cycle in California, we used our measurements of maximum V d during the day and median V d at night (V max d and V med d (night), respectively) to estimate the flux and lifetime of NO x to deposition in forests throughout the state during the day and night, respectively (Figs. 7 and 8).
The average deposition flux to trees in California was calculated via Eq. (11): Leaf area index (LAI) data for June 2018 were obtained from MCD15A2H Version 6 Moderate Resolution Imaging Spectroradiometer (MODIS) Level 4 product (Myneni et al., 2015) (Fig. S6 in the Supplement). The NO 2 surface concentrations and planetary boundary layer heights over California were obtained from a WRF-CHEM simulation for June 2014 (Fig. S6) . The month of June was chosen because in California this is when forests have a large LAI, large gross primary productivity (GPP), and the greatest sunlight availability, and ecosystems often experience water limitations in the later summer (Turner et al., 2020). Land cover data were obtained from NLCD Land Cover (CONUS) for 2016 (Yang et al., 2018) (Fig. S1 in the Supplement). Only forested sites were considered. Although the use of products from different years may introduce some error into our calculations, this will not qualitatively change our conclusion. Tree counts were obtained from the USDA Forest Service Forestry Inventory Analysis Database (Forest Inventory and Analysis, 2014) ( Fig. S2 in the Supplement). For each approximately 24 km 2 hexagonal plot (Bechtold, 2005) in the Forest Service Inventory that contained more than 50 % of the trees measured in our study, an effective deposition velocity to NO 2 (V eff d ) was calculated as a weighted (by tree species abundance) average from the V max d values listed in Table 2 (Fig. S6). Plots that contained less than 50 % of the trees measured were not considered. Data were interpolated to a 500 m grid. The resulting midday fluxes throughout California are shown in Fig. 7 and midnight fluxes are shown in Fig. 8. The greatest fluxes predicted are south of the San Francisco Bay Area, where there are high NO x concentrations and also a relatively high forest LAI for an urban region (Fig. S6). Similar hotspots can be seen near Los Angeles in the inland chaparral regions. Large fluxes are also predicted in the foothill forest region of the Sierra Nevada, where there is a large LAI, and frequent occurrences of P. sabiniana, the tree having the largest V d (Figs. S2 and S6). Relatively large fluxes occur in this region particularly during the nighttime.
The resulting lifetime of NO 2 to deposition was calculated via Eq. (12): where PBL is the planetary boundary layer height. The lifetimes to deposition during the day are shown in Fig. 7. In forested regions the lifetime to deposition is approximately 10 h. This relatively short lifetime may be especially consequential south of San Francisco Bay, where deposition could be competitive with the chemical sinks of HNO 3 and RONO 2 formation, which typically represent a lifetime to NO x loss of 2-11 h (e.g., Nunnermacker et al., 2000;Dillon et al., 2002;Alvarado et al., 2010;Valin et al., 2013;Romer et al., 2016;. The deposition fluxes and lifetimes to deposition during the night are shown in Fig. 8. With reduced deposition velocities at night, the nighttime deposition flux and the resulting total loss of NO 2 to deposition is small. However, with a reduced boundary layer during the night, the lifetime of NO x to deposition at night is on the same order as the deposition lifetime during the day (10-100 h) and the overall NO x lifetime at night. This indicates this loss pathway may be an important nighttime sink of NO x from the atmosphere and may affect the nighttime chemical NO x sinks of alkyl nitrate formation and N 2 O 5 chemistry (Brown et al., 2004(Brown et al., , 2006Crowley et al., 2010).
The estimations provided here are intended only to suggest qualitative indications of where NO x deposition may be important. Because we are ignoring effects of vertical transport and light attenuation through the canopy and because we are using maximum measured deposition velocities, the   a References refer to values in the reported max g w s column. b Study did not report value as a maximum stomatal conductance. The conductances shown are the maximum of the stomatal conductances reported in the cited study. c Theoretical calculation. deposition reported here is likely to be an upper-bound estimate. We recommend areas where this estimated deposition is highest as regions that should be the subject of future field and large-scale modeling studies.

Conclusions
We present measurements assessing the relative effects of stomatal diffusion and mesophyllic processing of NO 2 on the uptake rate of NO 2 . We find that the deposition velocity of NO 2 is essentially equal to the stomatal conductance to NO 2 under conditions of drought, excess soil nitrogen, variations in relative humidity, and during both the day and night. We find no evidence of any emission of NO 2 from leaves. NO 2 foliar exchange is thus unidirectional and variations are driven -from an atmospheric perspective -nearly entirely by the rate of diffusion through open stomata. This opens the possibility of using direct measurements of stomatal con-ductance -coupled with models and measurements of chemical transport, known relationships of the effects of environmental conditions on stomatal opening, measurements of canopy conductance -as well as indirect measurements, such as satellite solar-induced fluorescence data, to infer NO x foliar exchange. Additionally, we find significant differences in deposition velocities between species, reflecting differences in maximum stomatal conductance measurements that have been found by a number of previous studies (e.g., Ambrose et al., 2010;Maire et al., 2015;Henry et al., 2019;Murray et al., 2019). This diversity is not reflected in current atmospheric models and may have a meaningful impact on estimates of regional NO x fluxes and lifetimes. Our observations of stomatal opening in the absence of light also suggest foliar deposition may represent as much as 25 % of the total NO x loss at night, with stomatal deposition velocities as high as 0.038 cm s −1 . These findings not only have important implications for NO x chemistry but are also relevant for the atmosphere-biosphere exchange of other gases, such as CO 2 and biogenic volatile organic compounds.
Data availability. The data collected in this study can be obtained from the authors upon request.
Author contributions. ERD and BKP designed the experimental setup, and ERD, BKP, and AXL collected all NO 2 exchange data. BKP and ERD designed methods and collected data for nitrogen fertilization experiments. ERD and AXL designed methods and collected data for drought stress experiments. ERD performed data analysis, with assistance from AXL. ERD prepared the paper in consultation with RCC. RCC supervised the project.