Validation of urban NO2 concentrations and their diurnal and seasonal variations observed from space (SCIAMACHY and OMI sensors) using in situ measurements in Israeli cities

We compare a full-year (2006) record of surface air NO2 concentrations measured in Israeli cities to coinciding retrievals of tropospheric NO2 columns from satellite sensors (SCIAMACHY aboard ENVISAT and OMI aboard Aura). This provides a large statistical data set for validation of NO2 satellite measurements in urban air, where validation 5 is difficult yet crucial for using these measurements to infer NOx emissions by inverse modeling. Assuming that NO2 is well-mixed throughout the boundary layer (BL), and using observed average seasonal boundary layer heights, near-surface NO2 concentrations are converted into BL NO2 columns. The agreement between OMI and (13:45) BL NO2 columns (slope = 0.93, n=542), and the comparable results at 10:00 hrs for 10 SCIAMACHY, allow a validation of the seasonal, weekly, and diurnal cycles in satellitederived NO2. OMI and BL NO2 columns show consistent seasonal cycles (winter NO2 1.6–2.7× higher than summer). BL and coinciding OMI columns both show a strong weekly cycle with 45–50% smaller NO2 columns on Saturday relative to the weekday mean, reflecting the reduced weekend activity, and validating the weekly cycle 15 observed from space. The diurnal difference between SCIAMACHY (10:00) and OMI (13:45) NO2 is maximum in summer when SCIAMACHY is up to 40% higher than OMI, and minimum in winter when OMI slightly exceeds SCIAMACHY. A similar seasonal variation in the diurnal difference is found in the source region of Cairo. The surface measurements in Israel cities confirm this seasonal variation in the diurnal cycle. Using 20 simulations from a global 3-D chemical transport model (GEOS-Chem), we show that this seasonal cycle can be explained by a much stronger photochemical loss of NO2 in summer than in winter.

(OH) which is the main atmospheric oxidant.About 2/3 of the global source of NO x is from fossil fuel combustion (IPCC, 2007) and is concentrated in urban areas.Satellite observations of tropospheric NO 2 columns by solar backscatter from the GOME, SCIA-MACHY, and OMI sensors have provided increasing information over the past decade to quantify NO x emissions and their trends at various scales (e.g.Martin et al., 2003Martin et al., , 2006;;Richter et al., 2004;Konovalov et al., 2006;Wang et al., 2007;van der A et al., 2008).Validating these satellite observations is difficult, however.We show here that long-term surface NO 2 measurements from urban air quality monitoring networks provide valuable validation data and can test the seasonal, weekly, and diurnal variations seen from space by the SCIAMACHY and OMI sensors.
The atmospheric lifetime of NO x against oxidation is of the order of hours (Liang et al., 1998), resulting in large spatial and temporal variability around source areas.High concentrations are mainly restricted to the boundary layer up to 1-2 km, with rapid drop-off at higher altitudes.Drop-off is particularly sharp for NO 2 because of the decrease of the NO 2 /NO ratio with decreasing temperature (Martin et al., 2004).Validation of NO 2 column measurements from space has relied on aircraft vertical profiles coincident with satellite overpasses to capture the full extent of the column (Heland et al., 2002;Martin et al., 2004Martin et al., , 2006;;Bucsela et al., 2008).This is expensive and as a result validation data are few.Horizontal variability over the satellite footprint, aircraft traffic restrictions, and the need to extrapolate from the aircraft floor to the surface can further lead to significant errors in comparing aircraft-derived columns to the satellite column measurements.This is particularly problematic in urban areas which are of prime interest for validation (Boersma et al., 2008b).
Surface observations of NO 2 concentrations in urban areas offer an alternate vehicle for satellite validation.Most developed countries have dense urban monitoring networks for NO 2 because it is a regulated pollutant (Ellis, 1975;Hjelbrekke and Fjaeraa, 2008).The standard monitoring instrumentation is subject to interferences from NO x oxidation products (Gerboles et al., 2003;Steinbacher et al., 2007;Dunlea et al., 2007;Lamsal et al., 2008), and the lack of vertical information above the surface induces Introduction

Conclusions References
Tables Figures

Back Close
Full error.On the other hand, the high spatial and temporal density of the data greatly improve the validation statistics (Schaub et al., 2006;Ord ó ñez et al., 2006;Blond et al., 2007) over what is possible from aircraft profiles.The lack of vertical information in the surface data may in fact be less of a problem than the need to extrapolate aircraft observations down to the surface.We use the SCIAMACHY NO 2 retrievals (v1.1) available from the Tropospheric Emission Monitoring Internet Service (http://www.temis.nl)for the year 2006.The retrieval algorithm has been described in detail in Boersma et al. (2004), andBlond et al. (2007).In brief, slant column densities are determined with Differential Optical Absorption Spectroscopy (DOAS) in the 426.3-451.3nm spectral window.The stratospheric (background) slant column is estimated from data assimilation of slant columns in the global chemistry-transport model TM4 (Dentener et al., 2003).After subtraction of the stratospheric slant column, the residual tropospheric slant column is translated into a tropospheric vertical column by a tropospheric air mass factor that accounts for viewing geometry, the NO 2 and temperature profile shape, as well as the pressure and reflectivity of clouds and surface.The SCIAMACHY NO 2 retrieval (v1.1) is improved Introduction

Conclusions References
Tables Figures

Back Close
Full Screen / Esc

Printer-friendly Version
Interactive Discussion relative to earlier versions (v1.04) because it includes cloud parameters retrieved from FRESCO+.The FRESCO+ algorithm uses a more realistic physical model of the atmosphere by accounting for single Rayleigh scattering, and this leads to higher cloud fractions and cloud pressures.These in turn result in NO 2 columns that are lower for polluted situations with low clouds, leading to better agreement with independent NO 2 validation data (Wang et al., 2008).Errors in the slant column fitting, the stratospheric correction, and in the air mass factor lead to an overall error in the SCIAMACHY tropospheric NO 2 columns with an absolute component of approximately 0.65×10 15   Acarreta et al., 2004 for OMI).The a priori profile shapes are sampled from the same TM4 model, at 10:00 and 13:30 h for SCIAMACHY and OMI respectively (Boersma et al., 2008a).OMI retrieval errors have an absolute component of approximately 1.0×10 15 molecules cm −2 and a relative (AMF) component of 30% (Boersma et al., 2007).The spectral fitting error is smaller for SCIAMACHY than for OMI because of better SCIAMACHY signal-to-noise ratios, but OMI errors are generally much reduced by averaging over multiple days or over multiple pixels.

Surface NO 2 concentrations
We use here near-surface concentrations of NO 2 measured every 30 min at 8 sites in Israel for the whole year 2006.Figure 1 and Table 1 show the station locations.Measurements are with standard chemiluminescence instruments using molybdenum oxide converters.NO 2 is catalytically converted to NO on a heated ( 310• C) molybdenum surface and then measured as NO by chemiluminescence after reaction with ozone.The molybdenum converter also partly converts other oxidized nitrogen compounds such as nitric acid (HNO 3 ), alkyl nitrates, and peroxyacetyl nitrate (PAN) to NO (see e.g.Steinbacher et al., 2007;Dunlea et al., 2007 and references therein).This will cause NO 2 concentrations to be overestimated, especially in photochemically aged airmasses.Steinbacher et al. (2007) have shown biases up to +50% at rural locations downwind of pollution sources in Switzerland.For urban air in Mexico City, Dunlea et al. (2007) find no significant bias in the morning but increasing bias in the afternoon, presumably reflecting the diurnal formation of NO x oxidation products.
Recognizing that most NO y (with NO y defined as the sum of all reactive nitrogen oxides; NO+NO 2 +HNO 3 +organic nitrates) is NO x within urban areas at mid-morning with little prior photochemical activity, and given the lack of evidence for morninginterference in the study by Dunlea et al. (2007), we do not correct the surface concentrations measured at 10:00 (SCIAMACHY overpass).For measurements at 13:45 (OMI overpass) we show results from two approaches: Introduction

Conclusions References
Tables Figures

Back Close
Full Screen / Esc

Printer-friendly Version
Interactive Discussion 1. no correction, i.e. interpreting the NO 2 measurements as an upper limit for the true NO 2 concentrations; 2. correction based on the concurrently observed increase of O 3 from 10:00 to 13:45 as an indicator for the interference ∆[NO z ] (with NO z defined as NO y -NO x ), assuming∆[NO z ]=0.1•∆[O 3 ] on a molar basis based on Mexico City data (Dunlea et al., 2007), and with no correction if ∆[O 3 ] smaller or equal to 0.

Ensemble validation
Table 1 shows direct comparisons of OMI tropospheric NO 2 columns to surface NO 2 concentrations measured at 8 different stations in Israel at 13:30-14:00 h throughout 2006.From a reduced-major axis (RMA) regression analysis (Clarke, 1980), we find the lowest slopes (0.30, 0.59×10 15 molecules cm −2 /ppb) for the stations of Haifa and Ashdod that are situated immediately downwind of a large port (Haifa), and within 2 km of a large power plant and refinery complex (Ashdod).The other six stations are not influenced by such nearby NO x point sources.For Afula, no summertime data were available, which reduces the range of observed values, but other than that we did not see a good reason to exclude data from Afula from our analysis.Excluding data from Haifa and Ashdod (but maintaining Afula), OMI and surface NO 2 are well-correlated (r=0.63,n=542).Correcting the 13:45 h surface NO 2 observations for NO z interference as described above, leads to surface NO 2 concentrations that are on average smaller by 8%.We conclude that the effect of applying an interference correction for urban NO 2 in the early afternoon is modest, reducing surface-based NO 2 concentrations by approximately 0.7 ppb.
Table 2 shows comparisons of NO 2 concentrations measured at 10:00 to SCIA-MACHY tropospheric NO 2 .For the locations of Haifa and Ashdod, correlations and slopes are low, consistent with their situation downwind of strong point sources and with the results found for OMI (Table 1).The agreement between SCIAMACHY and Introduction

Conclusions References
Tables Figures

Back Close
Full surface NO 2 amounts to r=0.48 (n=173).Sharpening the spatial coincidence criterion to ±0.1 • , as for OMI, improves the correlation (r=0.55) between surface and SCIA-MACHY NO 2 , at the expense of the number of coincidences (n=29).We derive BL NO 2 columns from the near-surface NO 2 mixing ratios by assuming vertical uniformity in the BL and negligible concentrations above.We take the noontime mixing depths for Israel from a climatology of mixing depth observations (Table 2 in Dayan et al., 1988) with noontime mixing depths of [1.12, 0.87, 0.65, 0.95] km for winter, spring, summer, and fall respectively.Mixing depths in Israel are minimum in summer due to strong subsidence, and have an uncertainty of ≈25% (Dayan et al., 2002).Neglecting NO 2 above the BL is a reliable assumption (Martin et al., 2002;Boersma et al., 2008b).Figure 2 (left panel) compares tropospheric NO 2 columns retrieved from OMI and SCIAMACHY with interference-corrected BL NO 2 columns derived from near-surface measurements.BL NO 2 columns show a similar correlation with OMI (r=0.64,n=542) as the surface concentrations (r=0.63).The slope of the RMA regression is 0.93 ± 0.06 and the intercept is −0.37×10 15 molecules cm −2 , indicating that on average OMI does not deviate significantly from the BL NO 2 columns.Not correcting for interference gives similar agreement (r=0.67,slope 0.90±0.05)and a slightly larger negative smaller intercept (−0.72 ×10 15 molecules cm −2 ).The right panel of Fig. 2 shows the agreement between SCIAMACHY and BL NO 2 (r=0.46,n=173).Selecting pixels within ±0.1 We see that monthly mean NO 2 columns from the in situ observations and from OMI are similar within their error margins for all months.For the boundary-layer columns, error bars mainly reflect the uncertainty in the boundary-layer depths (25%, based on estimates from Dayan et al., 1988).If there is a bias in the BL-depths, it directly propagates into our BL column estimates, and this could explain the large differences between BL columns and OMI in winter.For OMI, error bars indicate the standard deviation of the mean.The seasonal cycle, expressed as the ratio of winter-to-summer mean NO 2 , is 2.7±0.6 in the boundary-layer NO 2 columns, and 1.6±0.5 in surface and OMI NO 2 .

Weekly cycle
The weekly cycle of coinciding ground-based and OMI NO 2 is shown in Fig. 4  and Jordan where sources are weaker and less likely to show a distinct weekend effect.In contrast, the small-pixel OMI observations used here are representative for the urban scale and more likely to show the local weekend emission reductions.

Diurnal cycle
Figure 5 shows average NO 2 columns measured at the locations of the Israeli stations on the same days at 10:00 and 13:45 in cloud-free conditions in the period March-November 2006.BL and satellite-retrieved NO 2 columns are consistent at 10:00 and at 13:45.The geometric mean of the ratio of SCIAMACHY and OMI NO 2 columns over Israeli cities is 1.22, and this tendency is quantitatively consistent with the geometric mean of 1.18 of 10:00 and 13:45 surface-based NO 2 columns.It is unlikely that these results are due to retrieval differences, as SCIAMACHY and OMI NO 2 columns have been retrieved with a consistent algorithm (Boersma et al., 2008a), and both datasets agree well with the independent data shown in Fig. 5.The results appear consistent with findings of Boersma et al. (2008a) over fossil fuel source regions at northern mid-latitudes, where SCIAMACHY observed 5-40% higher NO 2 than OMI in summer, explained by daytime photochemical loss of NO 2 , dampened by the diurnal cycle of anthropogenic NO x emissions that has a broad daytime maximum.Comparing 10:00 and 13:45 NO 2 observed on the same days during winter months in Israel did not show a clear diurnal cycle, and we will now investigate this further.lines in the inset.In summer, with shorter chemical lifetimes and smaller NO 2 columns than in winter, we see the reverse effect in Fig. 6c with 10:00-13:45 ratios of 1.29 over Israel and 1.16 over Cairo.Spring and autumn show diurnal variations in between the winter and summer extremes.Table 3 summarizes the mean 10:00-13:45 ratios for all seasons.
To evaluate this apparent seasonal cycle in the diurnal cycle, we compare ratios of NO 2 at 10:00 and 13:45 observed from satellites and from surface stations.We do not require spatial or temporal coincidence between satellite and surface measurements (as we did before), but only that both techniques measured under mostly-clear situations.The left panel of Fig. 7 shows significant 10:00-13:45 decreases in both satellite and surface-based observations over Israel in summer, and no decrease, or even a slight increase, in winter months.This is consistent with simulations from the 3-D chemistry-transport model GEOS-Chem (Bey et al., 2001) for the same location and time, as shown in the right panel of Fig. 7.In summer, GEOS-Chem simulates up to 21% higher NO 2 columns at 10:00, compared to smaller differences in winter.That the GEOS-Chem summertime ratios are smaller than observed by the two measurement techniques, is explained by the spatial extent of the GEOS-Chem grid cell ( 2• ×2.5 • ).
The GEOS-Chem grid cell covers not only the region of Israel where NO x sources are strong, but also a large part of the Mediterranean Sea where these sources are weak.
To bridge the gap between the simulation and the higher-resolution observations, we sampled the satellite measurements according to the boundaries of the GEOS-Chem grid cell.The seasonal cycle in the diurnal cycle observed at 2 • ×2.5 • (solid red curve) agrees much better with the GEOS-Chem simulations than the high-resolution curve (2 The diurnal variation of tropospheric NO 2 columns depends on the diurnal cycle of NO x emissions and chemical loss (see Eq. 3 in Boersma et al., 2008a).Over urban areas, NO x emissions show a daytime maximum that mainly reflects intense vehicle use between early morning and late afternoon.This pattern is very similar for all seasons.
In absence of strong seasonal changes in background NO x , the seasonal variation in the diurnal cycle of tropospheric NO 2 columns can therefore be understood mainly in terms of seasonally varying photochemistry.
Figure 8 shows a GEOS-Chem simulation over Israel of the diurnal variation in the NO 2 column and in the chemical loss rate constant k for a winter and a summer month in cloud-free conditions.k includes the chemical loss of NO x to HNO 3 through the gasphase reaction of NO 2 with OH and the hydrolysis of N 2 O 5 in aerosols.The GEOS-Chem model incorporating diurnally varying NO x emissions and photochemistry simulates up to 5× higher values and a much more pronounced diurnal cycle for k in August than in February 2006.In winter, the absence of a daytime maximum in NO x loss is expected to lead to an increase of NO x over the course of the day, because emissions are higher in the day than at night.The seasonal differences in k are consistent with higher OH concentrations in summer reflecting higher water vapour concentrations and more UV flux in that season (Spivakovsky et al., 2000).It is this strong daytime change in k (OH) in summer that is responsible for the distinct decrease in the NO 2 columns between 10:00 and 13:45 in summer.In winter, the chemical loss rate constant is generally too small throughout the day (low OH) to yield a net chemical loss, implying that NO x emissions may offset chemical loss during the daylight hours in that season.It is this mechanism that likely explains the observed wintertime 10:00-13:45 increases (expressed as ratios <1.0) in Fig. 7.

Conclusions
We showed that NO 2 concentration data from urban monitoring networks can provide useful validation for NO 2 column measurements from space.Previous validations of satellite retrievals relied primarily on profile measurements from aircraft, but these do not sample the lowest part of the atmosphere, where most of the NO 2 resides and strong gradients exist, and provide little statistics needed to overcome difficulties in vertical extrapolation and horizontal representativity.Ground-based measurements on the other hand sample NO 2 close to the source and offer the benefit of large, continu-Introduction

Conclusions References
Tables Figures

Back Close
Full Screen / Esc Printer-friendly Version Interactive Discussion ous datasets.
Our validation focused on urban locations in Israel in 2006.In situ measurements (10:00 and 13:45 h) of NO 2 concentrations show significant correlation with collocated SCIAMACHY and OMI tropospheric NO 2 column retrievals.We converted surface concentrations into columns representative for the boundary layer (BL) by assuming a well-mixed BL, and a climatology of observed Israeli BL depths.Afternoon surface NO 2 has been corrected for potential interference from oxidized nitrogen compounds (NO z ).This correction leads to 5-15% smaller NO 2 concentrations.The collocated 13:45 h BL and OMI NO 2 columns are in agreement (r=0.64,n=542, slope=0.93±0.07).A similar agreement is apparent in 10:00 h BL and SCIAMACHY NO 2 columns, provided they were collocated within 0.1 • .
The observed seasonal cycles in ground-based measurements and OMI NO 2 are consistent.BL and OMI NO 2 columns are 1.6-2.7×higher in winter than in summer, reflecting the strong seasonal variation in photochemistry and NO 2 lifetime.OMI observes a strong weekly cycle with approximately 50% smaller NO 2 columns on the Saturdays relative to the weekday mean.Ground-based techniques observe 46% lower NO 2 on Saturdays, similar as from space.
Outside the winter season, collocated satellite and ground-based measurements show higher NO 2 columns at 10:00 than at 13:45.In summer, the diurnal variation is strongest, with 10:00-13:45 ratios up to 1.4.In winter, the observed ratios are close to 1.0, and there is often a (small) increase in tropospheric NO 2 between mid-morning and early afternoon.Simulations with the global 3-D chemical transport model GEOS-Chem are consistent with the observed seasonal variation in the 10:00 vs. 13:45 difference.The simulations suggest that the high summertime and low wintertime ratios can be understood from a much stronger photochemical sink from oxidation by OH in summer than in winter.In winter, the photochemical sink is sufficiently small that continuous NO x emissions may offset chemical loss during the daylight hours, explaining the observed 10:00-13:45 increase in NO 2 .
The quality of the satellite retrievals at the urban scales shown here provides some Introduction

Conclusions References
Tables Figures

Back Close
Full confidence in top-down constraints on NO x emissions from the satellite data.But retrievals should be improved by using forward model parameters (surface albedo's, a priori profile shapes, and surface pressures) from data bases with horizontal resolutions similar to the satellite footprints.Aerosols contribute to retrieval uncertainty because they affect retrieved cloud fractions, and this complex interaction should be investigated further in order to reduce retrieval errors.

1
IntroductionNitrogen oxides (NO x ≡NO+NO 2 ) are of central importance for atmospheric chemistry.They are precursors of tropospheric ozone, aerosol nitrate, and the hydroxyl radical for five stations in the period between March and November 2006 (Afula was rejected because no summertime data were available).Wintertime observations have been excluded because in this season no coincidences were available for Saturdays.The weekly cycle is due to reduced NO x emissions on Fridays and Saturdays corresponding to the weekend pattern and days of rest.We observe reductions in NO 2 on Friday and Saturday in surface NO 2 concentrations (dashed line), corrected BL NO 2 columns (black diamonds), as well as in OMI NO 2 columns (red squares).The reduction on Saturday is significant, 50% relative to the weekday mean.Beirle et al. (2003) observed a similar Saturday NO 2 reduction over Israel on the regional scale from GOME.Their reduction (30%) is smaller because the relatively large GOME pixel (320×40 km 2 ) covers not just the strong source region of Israel, but also adjacent parts of the Mediterranean Sea Figures6a-dshow seasonally averaged differences between SCIAMACHY and OMI tropospheric NO 2 measured on the same day in mostly clear conditions over Israel and the surrounding regions in 2006 on a 0.5• ×0.5• grid.The inset shows the corresponding scatter plot with the Israeli grid cells marked in light blue.For comparison to a nearby industrialized region, we also plotted the Cairo grid cells, in red.In winter(DJF  2006), when NO 2 lifetimes are longest and NO 2 columns largest, we see in Fig.6athat SCIAMACHY tends to observe similar or lower NO 2 amounts at 10:00 than OMI at 13:45 h over the NO x source regions of Israel and Cairo.The geometrical mean ratios are <1.0 for both Israel and Cairo as indicated by the solid light blue and red

Fig. 7 .Fig. 8 .
Fig. 7. Seasonal cycle in the diurnal cycle of tropospheric NO 2 columns.Upper panel: cycle over the Israel source region from surface measurements (black line), and from SCIAMACHY and OMI measurements (red line).Lower panel: cycle over the 2 • ×2.5 • area of the GEOS-Chem grid cell from satellite measurements (solid red line) and from the GEOS-Chem model sampled at 10:00 and 13:30 h (solid blue line).
The black dashed line is for measurements coinciding with satellite data, i.e. measured under mostly-clear conditions only.Exclusion of cloudy scenes in the validation data subset decreases concentrations by 15% on average but does not affect the seasonal variation.Both OMI (red squares) and boundary-layer NO 2 columns (black symbols) are highest in the winter months, presumably reflecting the seasonal variation of NO 2 lifetime.
• , as for OMI, improves the correlation (r=0.53,slope=0.94)betweenBLandSCIAMACHYNO 2 columns.4Temporalvariability4.1 Annual cycleFigure3shows the seasonal variation of NO 2 over Israel in 2006 for the five stations with comparable NO 2 concentrations and seasonal variations and year-round coverage (for Afula, no summerdata were available).The grey dashed line is for all surface Introduction

Table 1 .
Summary of the comparison between surface NO 2 concentrations and OMI tropospheric NO 2 columns 1 at 8 Israeli cities in 2006.OMI tropospheric NO 2 columns (v1.02) had a pixel center within ±0.1• latitude/longitude of the station location, andwere measured under mostly-clear situations (cloud radiance fraction <50%).Only OMI observations with viewing zenith angles <35• were selected, ensuring that the satellite and surface observations coincide within 15 min.The OMI pixel centers were on average within 8 km of the station locations.Introduction 1

Table 2 .
Summary of the comparison between surface NO 2 concentrations and SCIAMACHY tropospheric NO 2 columns 1 at 8 Israeli cities in 2006.

Table 3 .
The 10:00-13:45 ratios 1 of tropospheric NO 2 columns for source regions in the Middle East.