Nitrogen oxides in the global upper troposphere : interpreting cloud-sliced NO 2 observations from the OMI satellite instrument

Nitrogen oxides (NOx o NO + NO2) in the upper troposphere (UT) have a large impact on global tropospheric ozone 20 and OH (the main atmospheric oxidant). New cloud-sliced observations of UT NO2 at 450-280 hPa (~6-9 km) from the OMI satellite instrument produced by NASA and KNMI provide global coverage to test our understanding of the factors controlling UT NOx. We find that these products offer useful information when averaged over coarse scales (20° × 32°, seasonal), and that the NASA product is more consistent with aircraft observations of UT NO2. Correlation with LIS/OTD satellite observations of lightning flash frequencies suggests that lightning is the dominant source of NOx to the upper troposphere except for 25 extratropical latitudes in winter. We infer a global mean NOx yield of 280 moles per lightning flash, with no significant difference between the tropics and mid-latitudes, and a global lightning NOx source of 5.5 Tg N a-1. There is indication that the NOx yield per flash increases with lightning flash footprint and with flash energy.


Introduction
Nitrogen oxides (NO x ≡ NO + NO 2 ) in the upper troposphere (UT) have profound effects on the oxidizing capacity of the atmosphere and on climate, but the factors controlling their concentrations are poorly understood.NO x in the UT impacts climate by efficiently producing ozone where it is a potent greenhouse gas (Dahlmann et al., 2011;Worden et al., 2011;Rap et al., 2015) and by increasing the concentration of OH (the main tropospheric oxidant) (Murray et al., 2012(Murray et al., , 2014)).Primary NO x sources in the UT include lightning, aircraft, convective injection, and downwelling from the stratosphere (Ehhalt et al., 1992;Jaeglé et al., 1998;Bertram et al., 2007).NO x cycles chemically with reservoir species including nitric acid (HNO 3 ), pernitric acid (HNO 4 ), dini-Published by Copernicus Publications on behalf of the European Geosciences Union.E. A. Marais et al.: Nitrogen oxides in the global upper troposphere trogen pentoxide (N 2 O 5 ), peroxyacetyl nitrate (PAN), and other organic nitrates, thus defining the NO y chemical family (NO y ≡ NO x + reservoirs).Effective loss of NO x from the UT is through subsidence of NO y to lower altitudes where deposition of HNO 3 provides the ultimate sink.The residence time of NO y in the UT is 10-20 days (Prather and Jacob, 1997).The lifetime of NO x against conversion to shortlived reservoirs varies from ∼ 3 h in the convective outflow of thunderstorms to 0.5-1.5 days in background air (Nault et al., 2016).Chemical recycling from these reservoirs maintains relatively high UT NO x background concentrations (Bradshaw et al., 2000;Baehr et al., 2003;Nault et al., 2016).
Representation of lightning NO x in chemical transport models (CTMs) is particularly uncertain.Physically based parameterizations relating lightning frequency to deep convective cloud tops, convective mass flux, convective precipitation, or high-cloud ice content have poor predictive capability (Tost et al., 2007;Allen et al., 2010;Murray et al., 2012;Finney et al., 2014), limiting our ability to estimate the response of lightning NO x to future climate (Finney et al., 2016(Finney et al., , 2018).An alternative is to prescribe flash densities from space-based observations and static NO x production rates per flash (Sauvage et al., 2007;Allen et al., 2010;Murray et al., 2012).NO x production efficiencies per flash in the literature vary from < 10 to 5000 moles nitrogen per flash (mol N fl −1 ) (Schumann and Huntrieser, 2007;Murray, 2016).Global CTMs typically use 100-500 mol N fl −1 , sometimes assuming higher production rates at midlatitudes than in the tropics (Hudman et al., 2007;Ott et al., 2010), and a global lightning NO x source of 3-7 Tg N a −1 to match observations of tropospheric ozone and NO y species (Sauvage et al., 2007).
Our understanding of UT NO x has so far been evaluated with observations from aircraft campaigns (Drummond et al., 1988;Jacob et al., 1996;Crawford et al., 1997;Jaeglé et al., 1998;Bradshaw et al., 2000;Hudman et al., 2007;Stratmann et al., 2016).There are also long-term NO x measurements from instruments on board commercial aircraft dating back to the 1990s, but these are mostly over the north Atlantic and the NO 2 measurements have low precision and interference from thermally unstable NO x reservoir compounds (Brunner et al., 2001).A number of studies have used satellite observations of tropospheric NO 2 columns from solar backscatter to infer lightning NO x emissions (Beirle et al., 2010;Pickering et al., 2016;Bucsela et al., 2010), including in combination with global models (Boersma et al., 2005;Martin et al., 2007;Miyazaki et al., 2014).These studies estimate global lightning NO x emission of 1 to 8 Tg N a −1 .
New cloud-sliced satellite products of tropospheric NO 2 mixing ratios at 280-450 hPa (∼ 6-9 km) offer additional vertical resolution by retrieving partial NO 2 columns above clouds and exploiting differences in heights of neighboring clouds to calculate NO 2 mixing ratios (Choi et al., 2014;Belmonte-Rivas et al., 2015).There are two new products of seasonal mean UT NO 2 mixing ratios retrieved from Ozone Monitoring Instrument (OMI) partial NO 2 columns by research groups at the Royal Netherlands Meteorological Institute (KNMI) and NASA.The KNMI product has been evaluated against UT NO 2 from the Tracer Model version 4 (TM4) CTM.Large regional differences between OMI and TM4 are attributed to model deficiencies in lightning NO x and uplift of anthropogenic pollution (Belmonte- Rivas et al., 2015).The NASA UT product is new to this work and follows a similar retrieval approach to the midtropospheric (900-650 hPa) product of Choi et al. (2014).That product was evaluated with aircraft observations of NO 2 and interpreted with the Global Modeling Initiative (GMI) CTM (Choi et al., 2014).Choi et al. (2014) identified large discrepancies between modeled and observed NO 2 seasonality over regions influenced by pollution and lightning.
Here, we compare the two UT NO 2 products, obtained with distinct retrieval methods, and use aircraft observations of NO 2 from multiple NASA DC8 aircraft campaigns to arbitrate and evaluate the information that can be derived from the satellite datasets.We go on to test current understanding of UT NO x and the implications for lightning emissions using the GEOS-Chem CTM.

OMI observations of upper troposphere NO 2
OMI is on board the NASA Aura satellite launched into Sunsynchronous orbit in July 2004.It has an overpass time of about 13:30 LT, a swath width of 2600 km, and a horizontal resolution of 13 km × 24 km at nadir (Levelt et al., 2006).Columns of NO 2 along the instrument viewing path (slant columns) are obtained by spectral fitting of solar backscattered radiation in the 405-465 nm window (Boersma et al., 2011).Standard products of total and tropospheric column NO 2 are screened for cloudy scenes using a cloud radiance fraction threshold of 0.5.Partial columns of NO 2 above cloudy scenes can be used to estimate vertically resolved NO 2 mixing ratios, as was first demonstrated with satellite observations of ozone (Ziemke et al., 2001).This approach, so-called cloud slicing, assumes a uniform trace gas concentration between two horizontally nearby clouds at different altitudes, so that the UT NO 2 mixing ratio is proportional to the slope of the partial columns versus the corresponding cloud pressures at the optical center of the cloud.Two products of seasonal mean UT NO 2 have been retrieved from OMI following distinct retrieval steps detailed below: a product from KNMI at 330-450 hPa for 2006 (Belmonte- Rivas et al., 2015) and from NASA at 280-450 hPa for 2005-2007 following an approach similar to that used to retrieve midtropospheric NO 2 (Choi et al., 2014).In what follows, we distinguish between the two OMI NO 2 products as KNMI and NASA.
The KNMI product uses DOMINO v2.0 slant columns (Boersma et al., 2011) and cloud fractions and altitudes from the O 2 -O 2 absorption cloud product, OMCLDO2 (Acarreta et al., 2004) over partially to very cloudy scenes (cloud radiance fraction > 0.5).Contamination due to NO 2 from below (up to 66 % over polluted land masses) is estimated using the TM4 model and removed.Stratospheric NO 2 from an assimilated product (Belmonte- Rivas et al., 2014) is also removed.An air mass factor (AMF) (detailed in Boersma et al., 2004) that accounts for viewing geometry, surface albedo, light attenuation by gases along the viewing path, and sensitivity to NO 2 vertical distribution is applied to the resultant partial slant columns to convert to vertical columns.Additional data filtering removes scenes with solar zenith angle (SZA) ≥ 70 • and surface albedo ≥ 30 %. Resultant daily vertical partial columns are aggregated on consistent pressure and horizontal (1 • ×1 • ) grids and used to determine seasonal mean UT NO 2 mixing ratios for grid squares with at least 30 measurements.UT NO 2 centered at 380 hPa (range 330-450 hPa) is estimated as the difference between partial tropospheric columns retrieved above two neighboring clouds with cloud pressures in the ranges 330-450 and 380-500 hPa, respectively.Biases from sampling cloudy scenes, such as the effect of clouds on photochemistry, are corrected using TM4.These are small (typically < 20 %) in the UT (Belmonte- Rivas et al., 2015).
The NASA UT NO 2 product for 2005-2007, centered at 350 hPa (∼ 280-450 hPa), uses updated version 3 slant columns (OMNO2 v3.0) (Krotkov et al., 2017) that correct for a positive bias in the DOMINO v2.0 product with improved spectral fitting (Marchenko et al., 2015;van Geffen et al., 2015).Partial columns from the cloud height to the top of the atmosphere are retrieved for individual OMI pixels above very cloudy scenes (cloud radiance fraction > 0.7) to minimize contamination from below.Cloud fraction and height are from the OMCLDO2 product (Acarreta et al., 2004).The AMF accounts for viewing path geometry and light scattering by clouds with uniform scatter that is optically thick and geometrically thin (near-Lambertian clouds).Data filtering is applied to remove scenes with SZA > 80 • , snow/ice cover, and severe aerosol pollution that could be misclassified as clouds.Daily UT NO 2 is estimated for neighboring partial columns with sufficient cloud variability (cloud pressure distance > 160 hPa) and well-mixed NO 2 (NO 2 vertical gradient < 0.33 pptv hPa −1 diagnosed with the GMI CTM).The stratospheric column is assumed uniform above neighboring clouds and thus is removed when differencing two nearby partial columns.Daily values of UT NO 2 are gridded to obtain seasonal means at 5 • × 8 • (latitude × longitude) for scenes with at least 50 measurements.Gaussian weighting is applied to assign higher weighting to UT NO 2 closest to 350 hPa.Choi et al. (2014) used a similar approach to retrieve midtropospheric NO 2 except that cloud fraction and height were from the rotational Raman scattering product (OM-CLDRR), and successful retrieval required a stricter cloud radiance fraction of 0.9, a minimum of 30 measurements, and a wider minimum cloud pressure distance of 200 hPa.A shift in cloud radiance fraction threshold from 0.9 (Choi et al., 2014) to 0.7 (this work) only introduces a small (< 5 %) difference in the retrieved partial columns due to contamination from below, as estimated by Pickering et al. (2016) for OMI scenes over the Gulf of Mexico.
Figure 1 compares seasonal mean UT NO 2 from the two satellite products in December-February and June-August.KNMI NO 2 is gridded to the NASA coarse grid.Data for March-May and September-November are in the Supplement (Fig. S1).KNMI NO 2 has greater coverage than the NASA product due to a lower cloud fraction threshold in the retrieval.The two products exhibit very different spatial features.Spatial correlation between the two products (Pearson's correlation coefficient between coincident grid squares) is R = 0.41 in December-February and R = 0.38 in June-August.There is marginal improvement in the correlation with further spatial averaging.At 20 • × 32 • , we find R = 0.50 in December-February and R = 0.45 in June-August.The correlation only increases substantially in September-November from R = 0.49 at 5 • × 8 • (Fig. S1) to R = 0.66 at 20 • × 32 • .KNMI is systematically lower than NASA in all seasons for coincident grid squares, varying from 16 % lower in June-August to 48 % lower in December-February at 20 • × 32 • .Contamination of UT NO 2 from below the cloud may still be present in the datasets despite attempts to correct for this using the TM4 model in the case of KNMI and by only considering very cloudy scenes in the case of NASA.These include a large enhancement in KNMI NO 2 (> 90 pptv) over southern Africa in June-August when there is intense biomass burning, and the NO 2 hotspot over northeast China  as the OMI products but we consider interannual variability to be a small source of error.All NO 2 measurements are from thermal-dissociation laser-induced fluorescence (TD-LIF) instruments (Day et al., 2002).These are susceptible to interference from decomposition of thermally unstable reservoir compounds including methyl peroxy nitrate (CH 3 O 2 NO 2 ) and HNO 4 (Browne et al., 2011).Publicly available DC3 and SEAC 4 RS TD-LIF NO 2 are already corrected for this interference.We apply a correction for the other campaigns using the relationship between temperature and percentage interference from Browne et al. (2011).Observed mean ambient air temperature in the UT during INTEX-A was 246 K, corresponding to 20 % interference.That for INTEX-B was 241 K (30 % interference) and 236 K for ARCTAS (38 % interference).
There are also NO 2 observations from the recent NASA ATom campaign and from the In-service Aircraft for a Global Observing System (IAGOS) commercial aircraft campaign (Berkes et al., 2018).These use chemiluminescence instruments that are also susceptible to interference.Chemiluminescence and TD-LIF NO 2 are consistent during the SEAC 4 RS campaign for the altitude range considered in this work (6-9 km) (Travis et al., 2016), but the interference from chemiluminescence is challenging to quantify due to dependence on the operator and instrument design that varies across campaigns (Reed et al., 2016).
Figure 2 shows the sampling extent of TD-LIF UT NO 2 over North America and Greenland in spring-summer at 450-280 hPa around the satellite overpass (11:00-16:00 LT) for scenes not influenced by the stratosphere (diagnosed with collocated ozone / CO > 1.25 mol mol −1 ; Hudman et al., 2007).Concentrations of UT NO 2 exceed 80 pptv over the eastern US due to lightning NO x emissions and convective transport of boundary layer pollution, and are < 30 pptv over the rest of the domain.
Figure 3 shows the spatial correlation between the March-August mean gridded aircraft data and the OMI UT NO 2 from the KNMI and NASA products as a function of horizontal resolution.There is no significant spatial correlation between the OMI products and aircraft NO 2 at 5 • × 8 • (R < 0.1) and 10 • × 16 • (R < 0.2).The correlation improves with further spatial averaging, peaking at 20 • × 32 • (R = 0.56 for KNMI, R = 0.64 for NASA).The satellite products are also spatially consistent at this resolution (R = 0.89), but KNMI is 43 % lower than NASA.
Figure 4 compares the spatial distribution of OMI and aircraft UT NO 2 at 20 • × 32 • over North America.Domain mean KNMI UT NO 2 is 38 % lower than the aircraft observations, compared to 2.2 % higher for NASA UT NO 2 .Both products exhibit less variability (reduced major axis, RMA, regression slopes < 1) and high bias in background NO 2 compared to the aircraft observations (positive RMA intercepts of 5.9 ± 1.4 pptv for KNMI and 9.2 ± 2.7 pptv for NASA).We proceed with the NASA UT NO 2 product at 20 • × 32 • , as correlation peaks at this resolution and the NASA product is more consistent with domain mean aircraft UT NO 2 than the KNMI product.

Constraints on upper tropospheric NO x
The NASA product provides near-global coverage of UT NO 2 to assess current understanding of regional UT NO x sources and dynamics by comparing to UT NO 2 from the GEOS-Chem CTM (version 10-01; http://wiki.seas.harvard.edu/geos-chem/index.php/GEOS-Chem_v10-01, last access: 28 November 2017) driven by NASA Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) reanalysis meteorology.The model horizontal resolution is 2 • × 2.5 • and the output is regridded to 20 • × 32 • for comparison with OMI.GEOS-Chem is sampled under all-sky conditions in the satellite overpass window (12:00-15:00 LT).We find that the effect on NO 2 of sampling the model under cloudy conditions is small.Isolating NO 2 under very cloudy conditions using MERRA-2 cloud fractions decreases modeled UT NO 2 by no more than 5 pptv in the tropics/subtropics and less at higher latitudes.We use output from the model for 2006 following a 1-year spin-up for chemical initialization.Interannual variability in UT NO 2 , determined as the difference between modeled 2006 and multi-year mean (2005)(2006)(2007) UT NO 2 , is small (< 4 pptv) everywhere except central Africa year round (4-12 pptv), the Arctic north of 60 • N (up to 25 pptv), and www.atmos-chem-phys.net/18/17017/2018/Atmos.Chem.Phys., 18, 17017-17027, 2018 the Middle East in June-August and northern India in March-May (both 10-20 pptv).Recent evaluation of model NO 2 with observed vertical profiles from the SEAC 4 RS aircraft campaign shows no significant bias in the 6-9 km range of the OMI product (Travis et al., 2016;Silvern et al., 2018).
Local GEOS-Chem emissions of NO x in the UT include aircraft and lightning.Aircraft emissions from the Aviation Emissions Inventory Code (AEIC) inventory (Stettler et al., 2011) total 0.82 Tg N in 2006, which is much less than lightning in the same year (6.5 Tg N).Transport from the stratosphere is simulated using a climatology of NO y species concentrations from the GMI model above the tropopause (Murray et al., 2012) and is very small (0.4 Tg N a −1 as total NO y ).Lightning in the model is estimated using the parameterization implemented by Murray et al. (2012).This includes an initial estimate of lightning flashes using the Price and Rind (1992, 1993, 1994) relationship between cloud-top height and lightning flashes.These are then scaled to the same annual global flash frequency (46 fl s −1 ) and regional distribution as the climatology from the combined Lightning Imaging Sensor (LIS) and Optical Transient Detector (OTD) high-resolution monthly climatology (LIS/OTD HRMC) (Cecil et al., 2014).The standard GEOS-Chem model has higher NO x yields per flash at northern midlatitudes (north of 35 • N) than in the tropics (500 mol N fl −1 versus 260 mol N fl −1 ), but we find that this overestimates observed OMI UT NO 2 by 10-20 pptv (20 %-40 %) at northern midlatitudes in summer when the lightning source is dominant.Here, we address this overestimate by assuming a NO x yield of 260 mol N fl −1 everywhere.This decreases global lightning NO x emissions by 15 % from 6.5 to 5.5 Tg N a −1 .The lightning parameterization in GEOS-Chem does not distinguish lightning NO x production from flashes within or between clouds (intra-or inter-cloud) or from the cloud to the Earth's surface (cloud to ground).
Figure 3 shows the spatial correlation between the model and aircraft observations.The model is more consistent with the aircraft observations than OMI at fine spatial resolution.Like OMI, GEOS-Chem correlation with the aircraft observations improves with spatial averaging, peaking at 20 • × 32 • (R = 0.75).Figure 4 also shows comparison of March-August GEOS-Chem UT NO 2 with the aircraft observations at 20 • × 32 • .The model is sampled over the same pressure range as NASA (280-450 hPa) around the OMI overpass (12:00-15:00 LT) and is filtered for stratospheric influence using model ozone / CO > 1.25 mol mol −1 .Domain average UT NO 2 from the model is 19% lower than the aircraft measurements, and the model also overestimates background UT NO 2 (intercept of 7.5 ± 1.0 pptv) and underestimates the variability (slope of 0.45 ± 0.09).
Figure 5 compares seasonal mean OMI and GEOS-Chem UT NO 2 in December-February and June-August.The other seasons are shown in the Supplement (Fig. S2).Formation of PAN, HNO 4 , and CH 3 O 2 NO 2 accounts for over 75 % of NO x loss in the model in all seasons.Lower concentrations of UT NO 2 in the Northern Hemisphere winter compared to summer in the model is mostly because lightning activity is at a minimum.The model underestimates UT NO 2 in the northern midlatitudes in winter by 20-40 pptv, suggesting misrepresentation of another process in the model, such as excessive NO x loss by N 2 O 5 hydrolysis in aerosols (Kenagy et al., 2018).The particularly large bias over polluted regions in winter could also be due to contamination of the UT NO 2 retrievals by enhanced boundary layer NO 2 .
Figure 6 shows the log-log relationship between seasonal mean LIS/OTD lightning flash climatology and seasonal mean UT NO 2 from OMI and GEOS-Chem, and the corresponding reduced major axis linear regression fits.Data are divided into northern midlatitudes and tropics.We exclude the contaminated observations over northeast China and the wintertime northern midlatitude grid squares that show no correlation with lightning flashes (R < 0.1).Results from multi-model sensitivity studies indicate that UT NO x in winter is predominantly from surface sources, with a smaller contribution from extratropical lightning (Grewe et al., 2001).Background concentration of UT NO 2 (intercepts in Fig. 6) from non-lightning sources is 10-20 pptv and is 3-5 pptv higher in the northern midlatitudes than in the tropics.The slopes for the linear fits to lightning flash frequency are consistent between the OMI observations and GEOS-Chem, and show similar slopes for northern midlatitudes and the tropics.Fitting the ratio between OMI obser-  vations and GEOS-Chem on the 20 • × 32 • grid implies a NO x yield per flash of 280 ± 80 mol N fl −1 with no significant difference between midlatitudes and the tropics, and no significant difference with the GEOS-Chem prior estimate of 260 mol N fl −1 .Our prior estimate of global lightning source was 5.5 Tg N a −1 , and the improved estimate is 5.9 ± 1.7 Tg N a −1 .
Properties of lightning flashes, including energy, duration, and footprint area, have been retrieved from the OTD and LIS sensors (Beirle et al., 2014).The flash footprint area is the spatial extent of lightning detection events contributing to the flash (collection of local events) diagnosed by the satellite data.Figure 7 shows the relationship between OMI and GEOS-Chem derived lightning NO x production rates and LIS lightning properties from Beirle et al. (2014).The strongest correlation is with lightning extent (R = 0.50), followed by energy (R = 0.40).The correlation with flash duration is weak (R = 0.25).The relationships in Fig. 7 suggest a dependence of lightning NO x production rates on lightning flash energy of 510 ± 80 mol N (J m −2 sr −1 µm −1 ) −1 and on flash footprint area of 2.0 ± 0.3 mol N km −2 , possibly offering guidance for relating NO x yields to physical properties in global models rather than the current approach of assigning static values.

Conclusions
Measurements of NO x in the UT have mainly been from aircraft campaigns that are limited in space and time.Two new cloud-slicing UT NO 2 products from OMI produced by KNMI and NASA offer the potential to address uncertainties in our understanding of UT NO x sources.Here, we intercompared these products, evaluated them with aircraft observations, and used them together with GEOS-Chem model simulations to demonstrate a dominance of lightning as a source of UT NO x .
The KNMI and NASA UT NO 2 products use very different retrieval methods.Seasonal mean concentrations from the two products show weak global correlation at the 5 × 32 • (R = 0.5-0.7).At that resolution, they show correlation with in situ aircraft observations of UT NO 2 over North America for different years (R = 0.56-0.64).The KNMI product is biased low by 38 % relative to the aircraft observations, while the NASA product has no significant bias.
We find from the relationship of OMI UT NO 2 with LIS/OTD flash rates that most NO x in the upper troposphere is from lightning, except in the midlatitudes in winter.The background NO 2 concentration in the absence of lightning is 10-20 pptv.The relationship suggests no difference in NO yields per flash between the midlatitudes and the tropics, in contrast to the higher yields at midlatitudes often assumed in models.We derive a global mean lightning NO x production rate per flash of 280 ± 80 mol N fl −1 , from which we infer a best estimate for the global lightning NO x emission of 5.9 ± 1.7 Tg N a −1 .

Figure 1 .
Figure 1.Upper troposphere (UT) NO 2 from the OMI satellite instrument.Seasonal mean UT NO 2 from KNMI in 2006 at 330-450 hPa (a, b) is compared to NASA in 2005-2007 at 280-450 hPa (c, d).Data are at 5 • × 8 • horizontal resolution for December-February (a, c) and June-August (b, d).Grey areas indicate no data and, for NASA, scenes with fewer than 50 measurements.

Figure 4 .
Figure 4. March-August upper troposphere NO 2 over North America.All data are at 20 • × 32 • .Background colors in the different panels show concentrations from KNMI, NASA, and GEOS-Chem (GC).Circles show the aircraft observations (same in all panels).Aircraft observations are for 11:00-16:00 LT.The model is sampled in the satellite overpass time window (12:00-15:00 LT).Model and aircraft data are at 280-450 hPa and screened for stratospheric influence using ozone / CO > 1.25 mol mol −1 .Inset boxes show reduced major axis (RMA) regression statistics and mean NO 2 for coincident grid squares.Grey grid squares indicate no coincident observations.

Figure 5 .
Figure 5. Observed and modeled upper troposphere NO 2 .The figure shows NASA OMI seasonal mean UT NO 2 for 2005-2007 (a, b) and corresponding GEOS-Chem model values (c, d).The model is sampled at 280-450 hPa during the satellite overpass (12:00-15:00 LT) and filtered for stratospheric influence.Data are at 20 • ×32 • horizontal resolution for December-February (a, c) and June-August (b, d).Grey grid squares in the top panels indicate no OMI data.

Figure 6 .
Figure 6.Log-log relationship between upper troposphere NO 2 and lightning flash frequencies, and linear regression fits between the two.Individual points are coincident seasonal mean UT NO 2 from OMI (a) and GEOS-Chem (b) versus seasonal mean LIS/OTD lightning flash climatologies for coincident 20 •× 32 • grid squares in the northern midlatitudes (> 30 • N; blue) and tropics (< 30 • N; orange).Northern midlatitude points exclude December-February, which show poor correlation with lightning flashes (see text for details).Lines and legends show reduced major axis linear regression fits to the data with corresponding Pearson's correlation coefficients.The regression lines plot as curves on the log-log scale, highlighting the NO 2 background at low lightning flash rates and the correlation of NO 2 with lightning at high flash rates.

Figure 7 .
Figure 7. Relationship between OMI and GEOS-Chem derived lightning NO x production rates and LIS lightning properties: energy (radiance), duration, and footprint area.Individual points are seasonal mean 20 • × 32 • grid squares at 40 • N-40 • S.
Belmonte-Rivas et al. (2015)5)caution that the contamination correction in the KNMI product relies on accurate simulation of NO 2 vertical distribution.