Articles | Volume 14, issue 1
Research article
01 Feb 2021
Research article |  | 01 Feb 2021

Estimating real driving emissions from multi-axis differential optical absorption spectroscopy (MAX-DOAS) measurements at the A60 motorway near Mainz, Germany

Bianca Lauster, Steffen Dörner, Steffen Beirle, Sebastian Donner, Sergey Gromov, Katharina Uhlmannsiek, and Thomas Wagner

In urban areas, road traffic is a dominant source of nitrogen oxides (NOx=NO+NO2). Although the emissions from individual vehicles are regulated by the European emission standards, real driving emissions often exceed these limits. In this study, two multi-axis differential optical absorption spectroscopy (MAX-DOAS) instruments on opposite sides of the motorway were used to measure the NO2 absorption caused by road traffic at the A60 motorway close to Mainz, Germany. In combination with wind data, the total NOx emissions for the occurring traffic volume can be estimated. Hereto, the ozone-dependent photochemical equilibrium between NO and NO2 is considered. We show that for 10 May 2019 the measured emissions exceed the maximum expected emissions calculated from the European emission standards for standardised test cycles by a factor of 11±7. One major advantage of the method used here is that MAX-DOAS measurements are very sensitive to the integrated NO2 concentration close to the surface. Thus, all emitted NO2 molecules are detected independently from their altitude, and therefore the whole emission plume originating from the nearby motorway is captured, which is a key advantage compared to other approaches such as in situ measurements.

1 Introduction

Nitrogen oxides (NOx) is a collective term for nitrogen dioxide (NO2) and nitric oxide (NO). In the troposphere, a photochemical reaction with ozone leads to an equilibrium state between NO2 and NO (Pandis and Seinfeld2006). About three-quarters of the global emissions of NOx originate from anthropogenic sources (IPCC2013). Moreover, nitrogen oxides not only play a major role in atmospheric chemistry but are also important in terms of air quality. The World Health Organization reports negative short-term as well as long-term exposure effects in pulmonary function and in other organs (World Health Organization2000). For this reason, the limitation of the concentration of nitrogen oxides is part of the European programme regarding ambient air quality and cleaner air (European Parliament and Council of the European Union2008).

Fossil fuel combustion from road traffic is a major contributor to NOx emissions. Hence, the European emission standards were introduced to regulate the exhaust emissions of new vehicles in the EU since 1998 (European Parliament and Council of the European Union1998) and tightened in 2007 by a new regulation bringing into force the so-called Euro 5 and Euro 6 norms (European Parliament and Council of the European Union2007) for passenger cars. New vehicles sold in the EU need to undergo a type-approval procedure, which verifies the compliance with these regulations. This procedure is standardised depending on the emission class, e.g. by the New European Driving Cycle (NEDC; European Parliament and Council of the European Union1970) and since 2017 by the Worldwide harmonised Light-duty vehicles Test Procedure (WLTP; Council of the European Union2017). These include the measurement of exhaust emissions on a chassis dynamometer. Similarly, there are standards applying to heavy-duty vehicles.

However, various studies (Carslaw et al.2011; Chen and Borken-Kleefeld2014) have shown that the real driving conditions are more dynamic than the tested driving cycles. In addition, it is known that several manufacturers have installed software that manipulates the test results by reducing emissions specifically during the test procedure (Borgeest2017). This results in increased exhaust emissions during normal driving operation.

In situ measurements such as those used in vehicle chasing experiments, e.g. performed by Pöhler and Engel (2019), directly measure the exhaust plume of individual vehicles. Others use remote sensing techniques (Carslaw et al.2011; Chen and Borken-Kleefeld2014) to measure exhaust gases across-road. Both approaches are able to resolve the emission of individual vehicles, but it is difficult to derive representative fleet-average emission factors, e.g. to compare these with expected emissions from models, as large data sets would be required.

Nevertheless, in the atmosphere NO and NO2 form an equilibrium state which is mainly influenced by the ozone concentration and solar irradiance but not the primary composition and amount of the exhaust gases. Thus, the multi-axis differential optical absorption spectroscopy (MAX-DOAS) yields a key advantage when operated at some distance from the emission source. This method is described in more detail in the next section. The presented results are based on one measurement day (10 May 2019) for proof of concepts. Further measurements could then be used to analyse, for example, different driving conditions in more detail.

2 Method

The MAX-DOAS method (Platt and Stutz2008) allows the measurement of the differential slant column density (DSCD) of different trace gases (Hönninger et al.2004). Hereto, spectra of scattered sunlight are recorded at different elevation angles using ground-based instruments. To convert the slant column density (SCD), which represents the integrated concentration along the slant light path, into the vertical column density (VCD), the so-called air mass factor (AMF) is needed. For trace gas layers close to the ground the geometric approximation for the AMF can be used (Hönninger et al.2004). The integrated trace gas concentration along the vertical path is then given by

(1) VCD = SCD AMF sin ( α ) SCD ,

where α is the elevation angle.

In order to remove the Fraunhofer lines, the logarithm of a so-called Fraunhofer reference spectrum with preferably minimal trace gas absorption is subtracted from the logarithm of the measured spectra. To fulfil this criterion, the reference spectrum is usually recorded with an elevation angle α=90, i.e. in zenith direction. It can be assumed that for a given solar zenith angle, the stratospheric absorption is constant for measurements at different elevation angles. Then, the differential SCD yields the integrated tropospheric concentration of a specific trace gas along the light path for an altitude range from the surface up to about 2 to 3 km (Frieß et al.2019, and references therein), i.e. the column density relative to the reference spectrum.

In this study, the MAX-DOAS method is used to quantify the NOx emissions of vehicles on a motorway. Using two MAX-DOAS instruments on the two sides of the motorway allows us to measure the background NO2 DSCDs on the upwind side and additionally the traffic-induced NO2 on the downwind side. The background NO2 DSCD is then subtracted from the NO2 DSCD on the downwind side and thus yields the NO2 SCD caused by the traffic emissions. In a final step, the derived NO2 SCD is converted into NOx emissions by combining it with wind data and assuming a steady-state NOx-to-NO2 ratio. These steps are described in detail below.

2.1 Experimental set-up

To retrieve the amount of NOx emitted by road traffic, two Tube MAX-DOAS instruments (Donner2016) were set up on each side of a motorway. With these instruments, it is possible to measure the NO2 DSCD of the ambient air along the viewing direction. The chosen measurement site is located along the heavily used A60 motorway close to Mainz, Germany, and has a long straight section which provides an advantageous geometry for the measurement set-up. The exact alignment of the instruments for the presented measurement day is depicted in Fig. 1 and shows that the viewing direction is northward and parallel to the lane of traffic.

Figure 1Alignment of the two Tube MAX-DOAS instruments on the measurement day, 10 May 2019. The instruments are located on both sides of the A60 motorway in Mainz, Germany, with a viewing direction parallel to the lane of traffic. The area between both instruments encloses the motorway and the railway track. Our measurements have highest sensitivity within the shaded green area. On the measurement day, continuous wind from westerly directions was present. Created with © Google Earth Pro (2018).

The chosen motorway section has a speed limit of 100 km h−1. The next access and exit are about 1 km in one direction and 1.5 km in the other direction. Acceleration and deceleration should, therefore, only have a minor effect at the measurement site.

On the measurement day, continuous westerly wind was present so that the air mass transport was almost perpendicular to the motorway as well as to the viewing direction of the instruments. From the difference between the upwind (Instrument A, west side) and downwind (Instrument B, east side) signals, the emissions of the motorway are estimated. The locations of the instruments were about 160 and 220 m to the west and east side of the motorway, respectively. Therefore, the area enclosed by the two Tube MAX-DOAS instruments contains the motorway section and a railway track. Possible sources of NOx are thus traffic emissions from cars, trucks and trains since no other sources (e.g. fires) were detected in the area.

Measurements were taken at an elevation angle of 20 and with a total integration time of 2 s. The short integration time favours a high temporal resolution even if the quality of the spectral fit (Sect. 2.2) decreases slightly at the same time. Assuming a plume height of up to 200 m, the area of highest sensitivity can be estimated. It is also depicted in Fig. 1 by the green shading. The choice of a rather high elevation angle not only constrains the sensitivity region but also decreases the influence of variations in the background signal by reducing the light path length in the lowermost atmosphere. It should be noted that there were broken clouds on the measurement day which possibly induce differences between the two instruments. This effect is further analysed in Sect. 2.3.

In addition, a camera and a weather station were positioned on the upwind side to obtain further information. Taking videos with this set-up makes it possible to observe the traffic density on the motorway. The weather station records the wind direction and wind velocity as well as several other meteorological parameters such as pressure and temperature every second.

2.2 Spectral analysis

The spectral analysis of the obtained spectra is performed using the QDOAS software (version 2.112.2; Danckaert et al.2012). As a reference, a series of 90 measurements were taken simultaneously with both Tube MAX-DOAS instruments at the upwind measurement site. In order to categorise differences between the two instruments (Sect. 2.3), reference measurements were taken at the same location after the measurement series was completed on both sides of the motorway. The wavelength calibration is accomplished using a high-resolution solar spectrum (Chance and Kurucz2010). For the analysis, a wavelength range of 400 to 460 nm was selected. The DOAS fit settings are summarised in Table 1. The spectral analysis is run separately for each instrument, yielding the NO2 DSCD time series for both measurement sites. Exemplary, the QDOAS fit for one spectrum (Instrument A, upwind) is depicted in Fig. A1.

Vandaele et al. (1998)Thalman and Volkamer (2013)Serdyuchenko et al. (2014)Rothman et al. (2010)(Kraus2003)

Table 1Settings of the spectral analysis.

Download Print Version | Download XLSX

2.3 Instrumental differences

To estimate the influence of instrumental differences between the two Tube MAX-DOAS instruments on the NO2 results, the reference spectra are investigated in more detail. These measurements were taken simultaneously with both instruments on the upwind side with an elevation angle of 90 (zenith view). Figure 2 shows the time series of the NO2 results for these spectra. The first 90 measurement of each instrument is taken as a reference.

Figure 2Time series of the NO2 results for the 90 measurements of both instruments on the upwind side. The spectra are analysed using the first 90 spectrum as a reference. The shaded grey area denotes the range where both measured similar NO2 DSCDs. In the lower panel the difference between the two Tube MAX-DOAS instruments is depicted, zoomed into the shaded grey area.


As can be seen, in the shaded grey area the standard deviation of the difference in the two instruments only amounts to

(2) Δ ( NO 2  DSCD ) = 0.4 × 10 14 molec cm - 2 .

Hereto, the data points of Instrument B were interpolated to the time axis of Instrument A. For the spectra after 15:05 UTC, the signal differs widely where the difference between both instruments is characterised by a standard deviation of 7.9×1014 molec cm−2. This increased deviation is due to clouds passing by (see Sect. B1). Thus, the measurements in the shaded grey area show that both instruments measure similar NO2 DSCDs for the same measurement conditions, i.e. the same set-up, viewing direction and cloud conditions. Therefore, these spectra are being accumulated to minimise noise and are used as fixed references, which assures that both instruments are analysed under the same conditions. The offset between the two instruments, which is visible in Fig. 2, originates from instrument properties and is constant in time. Thus, choosing the same reference period also compensates for this difference in the further processing of the data.

2.4 Integration time

In order to investigate the influence of the integration time on the spectral analysis, the fitting procedure is performed for spectra with different integration times but the same fit settings. Therefore, two or more spectra are added before performing the DOAS fit. The result of the NO2 retrieval as well as the fit error and the average root mean square (RMS) over each measurement series is depicted in Fig. 3. As expected, the NO2 fit error, which is given by the QDOAS analysis, shows the same trend as the RMS. For the short integration time of our measurements, the spectral residual of the fit is dominated by photon shot noise. This is also clearly demonstrated by the observed dependence of the RMS (and the fit error) on integration time. The RMS decreases for longer integration times as the ratio of the measured signal to the photon shot noise increases. In contrast to the fit error decreasing with integration time, the NO2 retrieval yields the same average NO2 DSCDs for different integration times. Nonetheless, for longer integration times with decreased noise levels (as expected), also systematic structures appear in the residuals for one of the instruments. Detailed analyses showed that they are almost identical for all spectra, and there is no systematic evolution with time. This is also consistent with the fact that there is no change in the result of the NO2 retrieval regardless of the integration time. Therefore, it can be concluded that these systematic errors are unproblematic for this set-up and analysis procedure. The standard deviation between the results for different integration times amounts to less than 0.2×1014 molec cm−2 for the east side and 0.7×1014 molec cm−2 for the west side measurements, which is 2 orders of magnitude smaller than the NO2 signal. Consequently, the measurements taken with an integration time of 2 s give sufficient results above the detection limit. This is preferable as high temporal resolution is necessary to resolve specific traffic events. In the following we focus on time-averaged emissions as for the presented measurement day no individual emission plumes could be identified. Nevertheless, it is conceivable that a detection is possible for lower traffic volume (e.g. on Sundays) or higher workload (e.g. motorway sections with higher slopes).

Figure 3Average NO2 DSCD, fit error and RMS for both measurement sites (Instrument A: west side, upwind; Instrument B: east side, downwind) for different integration times.


3 Results

3.1 Measurement results

The measurement results for 10 May 2019 are shown in Fig. 4. Panel (a) depicts the time series of the measured NO2 DSCDs for both the upwind and downwind side, analysed as described in Sect. 2.2. Panel (b) shows the difference between both signals

(3) SCD traffic = DSCD downwind - DSCD upwind .

A persistent offset is found with a mean value of

(4) SCD traffic = ( 0.18 ± 0.04 ) × 10 16 molec cm - 2

as represented by the orange line. The error is calculated using error propagation of the standard errors in the mean for both instruments and additionally includes the deviation Δ(NO2 DSCD) between both instruments as derived in Sect. 2.3. Moreover, an uncertainty of 16 % is added to account for the impact of the broken clouds on the measurement time series. In order to investigate the effect, a cloud filter is applied as discussed in Sect. B2.

Figure 4Time series of the measurement results of 10 May 2019. Panel (a) depicts the measured NO2 DSCD for both measurement sites (Instrument A: west side, upwind; Instrument B: east side, downwind). In panel (b) SCDtraffic, i.e. the difference between both signals, is shown. The orange line symbolises the average value. Panel (c) presents the traffic volume during the measuring duration. The number of vehicles was retrieved by counting from the videos over 1 min intervals on a sample basis. The dashed grey lines represent the times of passing trains. Panel (d) depicts the wind direction as measured by the weather station. Panel (e) shows the wind velocity vwind, perpendicular to the viewing direction of the Tube MAX-DOAS instruments. The orange line denotes the mean value over the whole measurement period. Here, the light-grey values depict the error Δvwind, of the calculated wind velocity.


As there are no large sources of NO2 other than the motorway close to the measurement site, the background NO2 DSCDs in both measurements can be assumed to be the same. Therefore, the difference between both sides is most likely due to traffic emissions. There seems to be no significant additional emission due to the passing trains (marked by the dashed grey lines in Fig. 4), although the railway next to the measurement site is only used by diesel trains. Temporal variations can be found in the derived difference in addition to the constant offset. Panel (c) depicts the amount of traffic observed during the measurement period for which the number of vehicles was counted over 1 min intervals on a sample basis using the recorded videos. Panel (d) presents the wind direction measured by the weather station on the upwind side. Panel (e) shows the wind velocity vwind, perpendicular to the viewing direction. The calculation is detailed in the next section.

3.2 Plume age

For a better understanding of the retrieved signal, the wind field needs further investigation. The quantity of interest is the wind velocity vwind, perpendicular to the viewing direction of the Tube MAX-DOAS instruments, whose viewing directions are assumed to be parallel to the motorway. Thereby, the age of the measured plume can be quantified, which is needed to retrieve the total emission (Sect. 3.3). The perpendicular wind velocity vwind, is shown in Fig. 4e and is calculated using the wind velocity and direction as measured by the weather station. From the alignment of the two Tube MAX-DOAS instruments as depicted in Fig. 1, it can be seen that the viewing direction corresponds to approx. 330. The perpendicular wind velocity is thus

(5) v wind , = v wind,meas cos ( ϕ wind ) ,


(6) ϕ wind = ϕ wind,meas - 330 + 90 ,

where ϕwind,meas is the measured wind direction at the weather station. The error can be calculated using the propagation-of-uncertainty principle. Hereto, the error in the wind velocity is estimated using the minimum and maximum values over 1 s (with a sampling rate of 4 Hz). An additional error for a possible misalignment of the weather station with regard to the viewing direction of the telescopes of 2 is considered. During the measurement period, the wind velocity perpendicular to the viewing direction is at ground level on average

(7) v wind , = ( 2.8 ± 1.0 ) m s - 1 .

Effects such as turbulence, especially in the vicinity of the motorway, and changing wind fields at plume height lead to uncertainties which can, however, not be readily quantified.

Taking into account the average distance between the motorway and the downwind instrument's viewing direction x=(195±25) m estimated from Fig. 1 within the main area of high sensitivity, an average age of an air parcel of

(8) t = ( 1.2 ± 0.4 ) min

can be obtained. However, variations in the wind velocity and wind direction on short timescales affect the transport of an air parcel. Therefore, the plume age cannot always be correctly represented by Eq. (8). The correlation between the wind field and the measured NO2 SCDs is further discussed in Sect. B3. Concluding, a constant wind velocity is favourable when applying this method.

3.3 Estimation of real driving emissions

To estimate the real driving emissions, first the mean NO2 SCD must be converted into a VCD using the geometric approximation as given in Eq. (1). Thus, for the elevation angle of (20±2), the AMF amounts to 2.9±0.3, and the measurement yields

(9) VCD traffic = ( 0.6 ± 0.1 ) × 10 19 molec m - 2 .

Multiplying this value by the average wind velocity perpendicular to the viewing direction, the measured emission of NO2 amounts to

(10) E meas , NO 2 = ( 1.8 ± 0.7 ) × 10 19 molec ( m s ) - 1 .

This value now describes the number of molecules emitted per metre and second along the motorway section. It is a direct quantity of the measurements and can be converted into emissions per vehicle per second by dividing by the number of vehicles per length of the motorway.

In combustion processes, N2 is mainly oxidised into NO, and in the atmosphere it is further oxidised into NO2 and other oxides of nitrogen (Pandis and Seinfeld2006), forming an equilibrium between NO and NO2. Especially diesel vehicles also directly emit NO2 (Carslaw et al.2011, and references therein). Therefore, to retrieve the total NOx emissions from the observed NO2 levels, the share of NO2 in total NOx at the measurement site has to be known.

In order to estimate the rate of NO-to-NO2 conversion, we used the CAABA (Chemistry As A Boxmodel Application) box model simulation with representative environment conditions and a road traffic source for the measurement period. CAABA uses the atmospheric chemistry model MECCA (Module Efficiently Calculating the Chemistry of the Atmosphere) that includes state-of-the-art chemical mechanisms (Sander et al.2019). A fraction of the traffic-emitted NO is photochemically equilibrated with air NO2 under daytime near-surface conditions. Hereto, the solar radiation is calculated for clear sky using the solar inclination at the measurement location.

One important factor regarding the conversion is the ambient ozone level as it regulates the photochemical NOx cycling and influences the resulting NO-to-NO2 repartitioning dynamics. Where the emission fluxes are very high, the titration of ozone stops further conversion of NO to NO2. However, turbulent mixing with ambient air increases with distance from the source, and ozone in the air parcel containing the plume is replenished. Thereby, the conversion of NO to NO2 continues. Our observations confirm that, for the presented measurement, ozone titration only prevails close to the emission source and thus has no significant influence on our measurements (see Sect. C1). Sufficiently high ambient ozone concentrations were measured at local environmental monitoring stations (42 to 44 ppb at Mainz-Mombach and Wiesbaden-Süd, distance to the measurement site approx. 5 km and 9 km, respectively; Umweltbundesamt2019b).

The corresponding NOx-to-NO2 conversion factor for the time t=(1.2±0.4) min an air parcel needs to get from the vehicle exhaust to the sensitivity region of the Tube MAX-DOAS instrument can be deduced to be f=2.4±1.0 (Sect. C1). The NOx emission is then derived using

(11) E meas , NO x = f E meas , NO 2 ,

which equals

(12) E meas , NO x = ( 4.3 ± 2.5 ) × 10 19 molec ( m s ) - 1 .

In case the equilibrium is already reached, a conversion factor of feq=1.5 needs to be applied instead. Then, the total NOx emission would amount to

(13) E meas , NO x , eq = ( 2.7 ± 1.1 ) × 10 19 molec ( m s ) - 1 .

The determination of the conversion factor f relies on the rather rough estimate of the age of an air parcel as well as the ozone concentration and chemical processes during the measurement period. Therefore, the equilibrium value gives an estimate which is independent of these factors. However, it is rather unlikely that the equilibrium state is reached so close to the emission source (as also found for airborne measurements of emission fluxes from power plants; Meier2018). Nonetheless, the emission value Emeas,NOx,eq is within the error in Emeas,NOx. In the following, the more realistic value of Emeas,NOx is taken for the comparison with the expected traffic emissions.

3.4 Expected traffic emissions and comparison to real driving emissions

To calculate the expected traffic emissions, the emission per vehicle needs to be computed. The limiting values for NOx emissions, as given by the European emission standards, are summarised in Table 2. The limiting values for passenger cars are given in NO2 equivalents per kilometre depending on the fuel type. For trucks, the values are reported in NO2 equivalents per kilowatt-hour. To undertake the following calculation, the emission standards of trucks need to be converted into limiting values per kilometre. Therefore, the values are multiplied by a conversion factor of 1.5±0.5 kWh km−1. This is composed of the fuel value 10.4 kWh L−1 of diesel fuel, the efficiency of a diesel engine of about 40 % and an average consumption for trucks of 36 L (100 km)−1 (Hilgers2016). The error accounts for varying fuel consumption of ±10 L (100 km)−1 and the uncertainty in the efficiency of the vehicle engine.

Table 2European emission standards for NOx emissions (Umweltbundesamt2019a).

Download Print Version | Download XLSX

The European emission standards are theoretical values for the allowed emissions of different pollutants. They are, however, not the expected emissions under real driving conditions. In order to bring the values in line, so-called real driving emission (RDE) conformity factors are used for new emission norms (Euro 6d-temp; Council of the European Union2016). To avoid inconsistencies, in the following only the European emission standards serve to estimate the theoretically expected emissions.

For the calculations, the statistical composition of the vehicle fleet is considered (see Table 3). The passenger car fleet is broken down by registration districts, fuel types and emission groups. To analyse the emission per vehicle, the statistical distribution of Rheinhessen-Pfalz is chosen. This also includes the city of Mainz and the Mainz-Bingen region. Note that in this area more cars with old emission standards (Euro 3 and 4) are registered compared to the average in Germany. The relative number of trucks is broken down by emission group only and relates to the distance travelled by German trucks. Attention should be paid to the fact that non-German trucks account for about 35 % of the total distance travelled in Germany (Kraftfahrt-Bundesamt2017).


Table 3Vehicle fleet composition by emission group in per cent.

Download Print Version | Download XLSX

From the emission standards and the statistical composition of the vehicle fleet, a theoretical emission per vehicle can be calculated. The weighted average of the emission limits amounts to

(14) E limit, cars = ( 116 ± 5 ) mg km - 1


(15) E limit, trucks = ( 1248 ± 277 ) mg km - 1

for passenger cars and trucks, respectively. The observed amount of traffic is deduced by counting the vehicles as shown in Fig. 4 and shows average values of

(16) N cars = ( 91 ± 4 ) min - 1


(17) N trucks = ( 6 ± 2 ) min - 1 .

The error estimation accounts for miscounting the number of vehicles on the video, e.g. when a truck shields the view of the other traffic lanes. Taking into account the average traffic volume, the theoretical total emission for the measuring period is given by

(18) E calc , NO x = N cars E limit, cars + N trucks E limit, trucks ,

which yields

(19) E calc , NO x = ( 0.4 ± 0.1 ) × 10 19 molec ( m s ) - 1 .

Here, NO2 equivalents are given for NOx emissions. Thus considering the molar mass of NO2 of 46 g mol−1 (Haynes2014), 1 mg of NOx emissions corresponds to 1.3×1019 molec.

The theoretical emissions calculated from the European emission standards for standardised test cycles can now be compared to the measured NOx emissions. Evidently, for 10 May 2019 the measured amount of NOx is larger than theoretically expected by a factor of 11±7. Even if an equilibrium state between NO and NO2 for the measured traffic emissions were assumed, the measured NOx emissions still show a higher value (by a factor of 7±3) compared to the calculated emissions. Moreover, in the very unlikely case that the exhaust gases primarily consist of NO2, and the measured NO2 difference directly equals the NOx emissions, this discrepancy remains unexplained. Possible error sources in the measurement cannot completely explain these differences.

As the traffic volume was relatively constant throughout the measurement period, it is more likely that the statistics do not reflect the vehicle fleet well enough, and/or a large part of the vehicles do not meet the emission standards. Here, it should be noted that the deviations of the actual vehicle composition from the assumed one cannot be the sole reason for this factor.

Assuming that only Euro 3 diesel cars and Euro III trucks, i.e. the technical status quo of the year 2000, were driving during the measurement period, the expected traffic emission would amount to

(20) E calc , NO x , Euro3/III = ( 2.0 ± 0.5 ) × 10 19 molec ( m s ) - 1 ,

which is still lower than the measured emission. As today only a minor fraction of all vehicles are registered as Euro 3 cars and Euro III trucks, this worst-case scenario is highly unlikely. Considering that especially non-German trucks more often drive with defective exhaust gas systems, these could lead to large emissions even exceeding the Euro III norm. Thereby, the discrepancy between the theoretical and measured emissions could be partly explained. However, trucks only account for parts of the total traffic volume. This again implies an excess of the European emission standards regarding NOx emissions also for a significant number of passenger cars.

3.5 Comparison to the Handbook Emission Factors for Road Transport database

The Handbook Emission Factors for Road Transport (HBEFA; version 4.1;  Notter et al.2019) provides emission factors for all current vehicle categories as weighted average values for Germany. To draw a comparison to the results deduced in the previous sections, the vehicle categories “passenger cars” and “heavy-duty vehicles” are used as these can be readily identified in the camera recordings of the motorway section. The aggregated emission factors for NOx especially show higher emissions of passenger cars as compared to the theoretical emission limits (see Sect. C2). This results in an average emission flux of

(21) E HBEFA , NO x = ( 1.1 ± 0.1 ) × 10 19 molec ( m s ) - 1 ,

which is roughly 3 times larger than expected from the European emission standards. Although the database provides modelled real driving emissions, there remains a discrepancy in the measurements of a factor of 4±2. In conclusion, our measurement method yields reasonable results and is able to quantify average emissions of the motorway section. Nonetheless, differences remain which cannot easily be attributed to a specific error source.

4 Conclusions

The measurement of NOx emissions at the A60 motorway close to Mainz, Germany, gives an estimate of real driving emissions. With two MAX-DOAS instruments set up on each side of a motorway, it is possible to retrieve the NO2 signal caused by road traffic and calculate the total NOx emissions for the occurring traffic volume.

The most uncertain aspect during the analysis of the data was the age of the measured plume on the downwind side. It directly affects the conversion factor f of the NOx-to-NO2 ratio and thus the final result of the measured emission (Eq. 11). To further investigate the effect of the plume age, it is favourable to set up several MAX-DOAS instruments downwind with different distances to the motorway. Thereby, the set-up of the instruments could be optimised, and the equilibrium state of NO2 for the given weather conditions can be measured. Hereto, a stable wind field is advantageous. This yields a more accurate conversion factor.

Other aspects such as the high ozone concentration and relatively constant wind are uncritical for the presented measurement day and allow us to apply a constant conversion factor f to the average emission. Although the changing cloud cover caused large fluctuations in the NO2 DSCDs, filtering the data leads to only slightly lower emissions. Consequently, this effect cannot explain the difference between the measured and expected emissions.

The main possible error source regarding the derivation of the expected NOx emissions is the difference from the assumed vehicle fleet to the measured vehicle fleet. Although the statistics are relevant to the Mainz region, the exact composition remains unknown. However, the worst-case calculation showed that the uncertainty in the vehicle fleet cannot explain the deviation from the measured emission. Presumably, a considerable number of vehicles did not meet the European emission standards. Moreover, it must be assumed that a substantial number of trucks are non-German vehicles. Recent studies showed that a large fraction of these vehicles had conspicuously high emissions, which indicates deactivated fuel cleaning units (Pöhler and Engel2019). These could also explain the temporal variations in the measured time series. Different driving conditions (e.g. the slope of the motorway section, the allowed speed limit, road works etc.) and the impact of the composition of the vehicle fleet could be investigated in more detail by applying this method at different measurement sites.

It can be concluded that the measured emissions on 10 May 2019 exceed the maximum expected emissions calculated from the European emission standards for standardised test cycles (Umweltbundesamt2019a) by a factor of 11±7. The comparison to the HBEFA database also indicates elevated emissions on that motorway section. This observation is in line with the work of other groups (Carslaw et al.2011; Chen and Borken-Kleefeld2014; Pöhler and Engel2019). In particular, the whole plume originating from the nearby motorway was measured rather than individual vehicle plumes, and hence the possibility that parts of the plume get overlooked can be neglected, which is a key advantage compared to other approaches such as in situ measurements.

Appendix A: QDOAS analysis

This section exemplarily includes a fit result (Fig. A1) of the QDOAS analysis for a spectrum of Instrument A (west side, upwind) at an elevation angle of 20 using 2 s integration time. The fit settings are specified in Table 1.

Figure A1Example fit result of the QDOAS analysis (Instrument A: west side, upwind; 11:59:11 UTC, 10 May 2019). The measured optical densities of different absorbers are depicted in red, whereas the fit results are depicted in black. The values in the titles refer to the resulting slant column densities in molec cm−2. The error in the NO2 fit amounts to 0.53×1015 molec cm−2.


Appendix B

B1 Effect of clouds on the reference spectra

Clouds can have a great impact on MAX-DOAS measurements. A change in light path is caused by the increased scattering probability in clouds as there are more particles compared to the ambient air. Furthermore, the wavelength dependency of the scattered light changes for particle-scattering processes compared to pure Rayleigh scattering. This effect already occurs for aerosols and is even more pronounced for clouds.

There are different methods to identify and classify clouds. Here, the temporal variation in the colour index (Wagner et al.2014) is used. The colour index (CI) is defined as the ratio of two radiance values at different wavelengths. In this case, the wavelengths 320 and 440 nm are chosen. Thereby, the wavelengths cover a large range to pronounce the effect of the wavelength dependency.

The CI is calculated for the 90 measurements (compare to Fig. 2), and the obtained temporal evolution is given in Fig. B1. An almost constant CI is expected for cloud-free conditions in this time period. It can be seen that measurements after 15:05 UTC were affected by clouds. This leads to larger deviations in the retrieved NO2 signal as shown in Sect. 2.3. The offset of the CI between the two instruments can be ascribed to the specific instrumental properties as the instruments are not absolutely radiometrically calibrated. Accordingly, the CI analysis also encourages the approach to use only 90 measurements in the shaded grey area as a reference.

Figure B1The temporal evolution of the colour index (CI; intensity ratio 320 nm  440 nm) for the 90 measurements, which were taken simultaneously on the upwind side, is depicted. The shaded grey area depicts the range where both instruments measured the same NO2 signal (compare to Fig. 2).


B2 Effect of clouds on the measurement result

Calculating the CI as described in Sect. B1 for all spectra, a characteristic behaviour can be seen (Fig. B2). As high temporal variation indicates cloud cover, all spectra where the CI is below the reference CIref are filtered. The reference was inferred by fitting a second-order polynomial to the data and is depicted as a dashed line. The filtered time series are displayed in Fig. B3. Recalculating the mean difference between the two measurement sites yields

(B1) SCD traffic, filtered = ( 0.156 ± 0.009 ) × 10 16 molec cm - 2 ,

which is about 16 % smaller compared to the unfiltered case.

Figure B2The colour index (CI; intensity ratio 320 nm  440 nm) for both measurement series. The dashed line (CIref) indicates the filter threshold.


Figure B3Analysis result of the NO2 DSCDs for both sides (blue: west side, upwind; red: east side, downwind) with applied cloud filter based on the colour index (CI; intensity ratio 320 nm  440 nm). The grey data points are filtered out. The resulting difference SCDtraffic is depicted in the lowermost panel, yielding slightly lower NO2 SCDs compared to the unfiltered case.


B3 Correlation to the wind field

Assuming a constant emission, the NO2 difference is expected to be reciprocal to the wind velocity. However, an air parcel is also affected by obstacles such as trees and follows the turbulent flow of air. Furthermore, the wind varies on timescales of less than 1 min, whereas the transport of the air parcel from the emission location to the sensitivity region of the MAX-DOAS instrument happens on larger timescales of 1 min or more. This means that the time of the wind measurement and the time of the NO2 measurement are shifted by a time difference in which the wind might change strongly. Hence, the age of the air parcel cannot always be correctly represented by the simple calculation in Eq. (8).

To test this hypothesis, both the wind measurements and the time series of the NO2 differences are averaged over a time period of 12 min. Figure B4 shows the correlation between both quantities (R2=0.365). The data points are fitted using the linear least squares (LLS) method (orange line) as well as using the orthogonal distance regression (ODR; green line). Here, ODR is able to take into account the standard errors in the mean values in the fitting procedure (Cantrell2008). In doing so, the slope of the fit increases and at the same time the intercept decreases. Comparing the fit results with the obtained emission Emeas,NO2 over the complete NO2 measurement series as described in Sect. 3.3, a slope of about 5000±2000 molec (m s)−1 is expected. The fits from Fig. B4 show slopes of 4230±208 molec (m s)−1 for the LLS and 7539±2013 molec (m s)−1 for the ODR method which are in agreement with the expected value.

Figure B4Correlation between the inverse of the wind velocity vwind,-1 perpendicular to the viewing direction of the Tube MAX-DOAS instruments and the NO2 signal (SCDtraffic) for a 12 min averaging time span. The data points were fitted using the linear least squares (LLS) method and orthogonal distance regression (ODR).


Nevertheless, the weak correlation is not completely surprising because of the low variability in the wind velocity. Moreover, a constant wind velocity is generally advantageous for the measurements.

Appendix C

C1 CAABA–MECCA simulation

Figure C1 presents the results of the plume emission simulation using the CAABA–MECCA box model (Sander et al.2019). Applying representative environmental conditions for the measurement period, road emission was approximated with a 10 s emission pulse of NO into the box. From the evolution of changes in NO and NO2 abundances in the air parcel with time, the NO2-to-NOx ratio of the plume is deduced.

To analyse the possible effect of ozone titration, a Gaussian dispersion model is applied. It uses Pasquill stability classes (Pandis and Seinfeld2006) based on the atmospheric stability of the measurement day. With this dispersion model, the extent of the emission plume is estimated, and the NO2 mixing ratio from our measurements is calculated. While turbulence induced by the local topography and obstacles like trees is neglected, it helps to estimate the evolution of the NO2 mixing ratio between emission source and measurement location. From the comparison of the dispersion model and the observations, it can be concluded that the ozone-poor chemical regime only prevails close to the emission source.

In order to consider this in the emission estimate calculation, the transport of the air parcel containing the plume is subdivided into two sections: (1) close to the emission source we assume that only negligible amounts of NO are converted into NO2, and no further conversion takes place as ozone is depleted; (2) turbulent mixing with ambient air refills the ozone reservoir, and NO-to-NO2 conversion can be described by the CAABA box model simulations. For simplicity, the distance which corresponds to the NO2 mixing ratio matching the one simulated in the box model is chosen as the transition between both sections. Thereby, the time for NO-to-NO2 conversion is shorter than without consideration of ozone limitations. The resulting NOx-to-NO2 ratio of the measured air parcel is estimated to be f=2.4±1.0.

C2 The HBEFA database

The HBEFA (version 4.1; Notter et al.2019) database provides emission factors for various common vehicle types (passenger cars, light- and heavy-duty vehicles, buses and coaches, and motorbikes). Here, they differentiate by emission standard (Euro 0 to Euro VI) and by different traffic situations. The database includes not only all regulated pollutants but also a number of non-regulated pollutants, including CO2 and fuel and energy consumption.

In Sect. 3.5 the aggregated NOx emission factors are considered for the vehicle categories “passenger cars” and “heavy-duty vehicles” regarding the year 2020. These yield values of

(C1) E HBEFA, cars = 499 mg km - 1


(C2) E HBEFA, trucks = 1426 mg km - 1

for the two vehicle categories, respectively. Applying the same calculations as in Sect. 3.4, i.e. taking into account the number of vehicles during the measuring period and using Eq. (18) with the emission values given by the HBEFA database, the emission flux can be estimated at

(C3) E HBEFA , NO x = ( 1.1 ± 0.1 ) × 10 19 molec ( m s ) - 1 .

Figure C1CAABA–MECCA box model simulation for the presented measurement day using representative environmental conditions. At 15:00 simulated local time, a 10 s emission of NO into the box is performed, representing the emission from road traffic. Panel (a) shows changes in NOx=NO+NO2 and ozone (O3) compared to the background values. Panel (b) depicts the NO2-to-NOx as well as NO-to-NOx ratio in the plume.


Data availability

Measurement data are provided in the Supplement.


The supplement related to this article is available online at:

Author contributions

TW, SDö and BL designed the experiment. Adaptation of the instruments to the measurement set-up was implemented by SDo, SDö and BL. SDö, BL and KU performed the measurements. SG developed and performed the simulations. BL prepared the manuscript with contributions from all co-authors. TW, SB, SDö and SDo contributed with valuable feedback and supervised the study.

Competing interests

Thomas Wagner is a member of the editorial board of Atmospheric Measurement Techniques.


We acknowledge the electronics workshop (i.e. Thomas Klimach, Christian Gurk, Mark Lamneck and Frank Helleis) and the mechanical workshop (i.e. Michael Dietrich and Ralf Wittkowski) of the Max Planck Institute for Chemistry, Mainz, for the continuous support in the development of the Tube MAX-DOAS instrument. We are also thankful to Denis Pöhler (Airyx GmbH) for sharing his expertise on traffic emission estimation.

Financial support

The article processing charges for this open-access publication were covered by the Max Planck Society.

Review statement

This paper was edited by Michel Van Roozendael and reviewed by two anonymous referees.


Borgeest, K.: Manipulation von Abgaswerten: Technische, gesundheitliche, rechtliche und politische Hintergründe des Abgasskandals, Springer-Verlag,, 2017. a

Cantrell, C. A.: Technical Note: Review of methods for linear least-squares fitting of data and application to atmospheric chemistry problems, Atmos. Chem. Phys., 8, 5477–5487,, 2008. a

Carslaw, D. C., Beevers, S. D., Tate, J. E., Westmoreland, E. J., and Williams, M. L.: Recent evidence concerning higher NOx emissions from passenger cars and light duty vehicles, Atmos. Environ., 45, 7053–7063,, 2011. a, b, c, d

Chance, K. and Kurucz, R. L.: An improved high-resolution solar reference spectrum for earth's atmosphere measurements in the ultraviolet, visible, and near infrared, J. Quant. Spectrosc. Ra., 111, 1289–1295,, 2010. a

Chen, Y. and Borken-Kleefeld, J.: Real-driving emissions from cars and light commercial vehicles–Results from 13 years remote sensing at Zurich/CH, Atmos. Environ., 88, 157–164,, 2014. a, b, c

Council of the European Union: Commission Regulation (EU) 2016/427 of 10 March 2016 amending Regulation (EC) No. 692/2008 as regards emissions from light passenger and commercial vehicles (Euro 6), Official Journal of the European Union, L82, 1–98, 2016. a

Council of the European Union: Commission Regulation (EU) 2017/1151 of 1 June 2017 supplementing Regulation (EC) No. 715/2007 of the European Parliament and of the Council on type-approval of motor vehicles with respect to emissions from light passenger cars and commercial vehicles (Euro 5 and Euro 6), Official Journal of the European Union, L175, 1–643, 2017. a

Danckaert, T., Fayt, C., Van Roozendael, M., De Smedt, I., Letocart, V., Merlaud, A., and Pinardi, G.: QDOAS Software user manual, BIRA-IASB, Brussels, Belgium, 2012. a

Donner, S.: Mobile MAX-DOAS measurements of the tropospheric formaldehyde column in the Rhein-Main region, Master's thesis, University of Mainz, available at: (last access: 28 January 2021), 2016. a

European Parliament and Council of the European Union: 70/220/EEC of 20 March 1970 on the approximation of the laws of the Member States relating to measures to be taken against air pollution by gases from positive-ignition engines of motor vehicles, Official Journal of the European Union, L76, 1–22, 1970. a

European Parliament and Council of the European Union: 98/69/EC of the European Parliament and of the Council of 13 October 1998 relating to measures to be taken against air pollution by emissions from motor vehicles and amending Council Directive 70/220/EEC, Official Journal of the European Union, L350, 12, 1998. a

European Parliament and Council of the European Union: Regulation (EC) No 715/2007 of the European Parliament and of the Council of 20 June 2007 on type approval of motor vehicles with respect to emissions from light passenger and commercial vehicles (Euro 5 and Euro 6) and on access to vehicle repair and maintenance information, Official Journal of the European Union, L171, 1–16, 2007. a

European Parliament and Council of the European Union: Council Directive 2008/50/EC on ambient air and cleaner air for Europe, Official Journal of the European Union, L151, 1–44, 2008. a

Frieß, U., Beirle, S., Alvarado Bonilla, L., Bösch, T., Friedrich, M. M., Hendrick, F., Piters, A., Richter, A., van Roozendael, M., Rozanov, V. V., Spinei, E., Tirpitz, J.-L., Vlemmix, T., Wagner, T., and Wang, Y.: Intercomparison of MAX-DOAS vertical profile retrieval algorithms: studies using synthetic data, Atmos. Meas. Tech., 12, 2155–2181,, 2019. a

Google Earth Pro: Mainz, 495855.48′′ N 81234.81′′ E, eye altitude 4.98 km, Google 2018, GeoBasis-DE/BKG 2009, 2018. a

Haynes, W. M.: CRC handbook of chemistry and physics, CRC press, Boca Raton, Florida, USA; London, UK; New York, USA, 2014. a

Hilgers, M.: Kraftstoffverbrauch und Verbrauchsoptimierung, Springer,, 2016. a

Hönninger, G., von Friedeburg, C., and Platt, U.: Multi axis differential optical absorption spectroscopy (MAX-DOAS), Atmos. Chem. Phys., 4, 231–254,, 2004. a, b

IPCC: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., Cambridge, UK and New York, USA, 2013. a

Kraftfahrt-Bundesamt: Lastfahrten im Inlandsverkehr nach Hauptverkehrsbeziehungen mit europäischen Lastkraftfahrzeugen im Jahr 2017, available at: (last access: 18 July 2019), 2017. a

Kraftfahrt-Bundesamt: FZ 1 Bestand an Kraftfahrzeugen und Kraftfahrzeuganhängern nach Zulassungsbezirken, 1. Januar 2019, available at: (last access: 29 January 2021), 2019a. a, b

Kraftfahrt-Bundesamt: VD 5 Verkehr deutscher Lastkraftfahrzeuge Gesamtverkehr Mr̈z 2019, available at: (last access: 29 January 2021), 2019b. a, b

Kraus, S.: DOAS Intelligent System, institute of Environmental Physics, University of Heidelberg, Cooperation with Hoffmann Messtechnik GmbH, Shaker, Aachen, Germany, 2003. a

Meier, A. C.: Measurements of horizontal trace gas distributions using airborne imaging differential optical absorption spectroscopy, PhD thesis, Universität Bremen, Bremen, Germany, 2018. a

Notter, B., Keller, M., Althaus, H.-J., Cox, B., Knörr, W., Heidt, C., Biemann, K., Räder, D., and Jamet, M.: HBEFA (Handbook Emission Factors for Road Transport) 4.1 Development Report, available at: (last access: 18 August 2020), 2019. a, b

Pandis, S. N. and Seinfeld, J. H.: Atmospheric chemistry and physics: From air pollution to climate change, John Wiley & Sons, Inc., New York, USA, 2006. a, b, c

Platt, U. and Stutz, J.: Differential Optical Absorption Spectroscopy: Principles and Applications, Springer Science & Business Media, Berlin, Heidelberg, Germany, 2008. a

Pöhler, D. and Engel, T.: Bestimmung von realen Lkw NOx-Emissionen (Real Driving Emissions) und hohen Emittern auf deutschen Autobahnen, Universität Heidelberg, Heidelberg, Germany, 2019.  a, b, c

Rothman, L., Gordon, I., Barber, R., Dothe, H., Gamache, R., Goldman, A., Perevalov, V., Tashkun, S., and Tennyson, J.: HITEMP, the high-temperature molecular spectroscopic database, J. Quant. Spectrosc. Ra., 111, 2139–2150,, 2010. a

Sander, R., Baumgaertner, A., Cabrera-Perez, D., Frank, F., Gromov, S., Grooß, J.-U., Harder, H., Huijnen, V., Jöckel, P., Karydis, V. A., Niemeyer, K. E., Pozzer, A., Riede, H., Schultz, M. G., Taraborrelli, D., and Tauer, S.: The community atmospheric chemistry box model CAABA/MECCA-4.0, Geosci. Model Dev., 12, 1365–1385,, 2019. a, b

Serdyuchenko, A., Gorshelev, V., Weber, M., Chehade, W., and Burrows, J. P.: High spectral resolution ozone absorption cross-sections – Part 2: Temperature dependence, Atmos. Meas. Tech., 7, 625–636,, 2014. a

Thalman, R. and Volkamer, R.: Temperature dependent absorption cross-sections of O 2–O 2 collision pairs between 340 and 630 nm and at atmospherically relevant pressure, Phys. Chem. Chem. Phys., 15, 15371–15381,, 2013. a

Umweltbundesamt: Emissionsstandards, available at:, last access: 6 July 2019a. a, b

Umweltbundesamt: Aktuelle Luftdaten, Fachgebiet II 4.2, Beurteilung der Luftqualität, ozone concentration of 10 May 2019 at Mainz-Mombach (DERP007) and Wiesbaden-Süd (DEHE022), available at:, last access: 22 May 2019b. a

Vandaele, A. C., Hermans, C., Simon, P. C., Carleer, M., Colin, R., Fally, S., Merienne, M.-F., Jenouvrier, A., and Coquart, B.: Measurements of the NO2 absorption cross-section from 42 000 cm−1 to 10 000 cm cm−1 (238–1000 nm) at 220 K and 294 K, J. Quant. Spectrosc. Ra., 59, 171–184, 1998. a

Wagner, T., Apituley, A., Beirle, S., Dörner, S., Friess, U., Remmers, J., and Shaiganfar, R.: Cloud detection and classification based on MAX-DOAS observations, Atmos. Meas. Tech., 7, 1289–1320,, 2014. a

World Health Organization: Air quality guidelines for Europe, Regional Office for Europe, edited by: Theakston, F., World Health Organization, Bonn, Germany, 2000. a

Short summary
In urban areas, road traffic is a dominant source of nitrogen oxides. In this study, two multi-axis differential optical absorption spectroscopy (MAX-DOAS) instruments on opposite sides of a motorway were used to measure the nitrogen dioxide absorption near Mainz, Germany. Total nitrogen oxide emissions are estimated for the occurring traffic flux. We show that the measured emissions systematically exceed the maximum expected emissions calculated from the European emission standards.