Tropospheric NO2 measurements using a three-wavelength optical parametric oscillator differential absorption lidar

The conventional two-wavelength differential absorption lidar (DIAL) has measured air pollutants such as nitrogen dioxide (NO2). However, high concentrations of aerosol within the planetary boundary layer (PBL) can cause significant retrieval errors using only a two-wavelength DIAL technique to measure NO2. We proposed a new technique to obtain more accurate measurements of NO2 using a three-wavelength DIAL technique based on an optical parametric oscillator (OPO) laser. This study derives the three-wavelength DIAL retrieval equations necessary to retrieve vertical profiles of NO2 in the troposphere. Additionally, two rules to obtain the optimum choice of the three wavelengths applied in the retrieval are designed to help increase the differences in the NO2 absorption crosssections and reduce aerosol interference. NO2 retrieval relative uncertainties caused by aerosol extinction, molecular extinction, absorption of gases other than the gas of interest and backscattering are calculated using two-wavelength DIAL (438 and 439.5 nm) and three-wavelength DIAL (438, 439.5 and 441 nm) techniques. The retrieval uncertainties in aerosol extinction using the three-wavelength DIAL technique are reduced to less than 2 % of those when using the two-wavelength DIAL technique. Moreover, the retrieval uncertainty analysis indicates that the three-wavelength DIAL technique can reduce more fluctuation caused by aerosol backscattering than the two-wavelength DIAL technique. This study presents NO2 concentration profiles which were obtained using the HU (Hampton University) threewavelength OPO DIAL. As a first step to assess the accuracy of the HU lidar NO2 profiles, we compared the NO2 profiles to simulated data from the Weather Research and Forecasting Chemistry (WRF-Chem) model. This comparison suggests that the NO2 profiles retrieved with the threewavelength DIAL technique have similar vertical structure and magnitudes typically within±0.1 ppb compared to modeled profiles.


Introduction
Nitrogen dioxide (NO 2 ) plays a critical role in the tropospheric chemistry and is one of the reactive gases collectively referred to as "nitrogen oxides" (NO x = nitric oxide and nitrogen dioxide, NO + NO 2 ) (U.S. EPA, 2020). The sources of NO x emissions include transportation (on-road vehicles, airplanes, trains, ships), wood burning, industrial and chemical processes, activities for oil and gas development, soil emissions, lightning, and wildfires (see Nitrogen Oxides Emissions indicator) (U.S. EPA, 2020). Once emitted, NO reacts rapidly in the presence of ozone to form NO 2 . In US urban locations, most measured airborne NO 2 comes from the reaction of these two precursors rather than from Published by Copernicus Publications on behalf of the European Geosciences Union.
direct NO 2 emissions (Bertram, et al., 2005;Beirle, et al., 2011). Scientific evidence indicates that short-term NO 2 exposure, ranging from 30 min to 24 h, can cause the exacerbation of asthma symptoms, in some cases resulting in hospitalization (Berglund, et al., 1993). Long-term NO 2 exposure is likely to have a causal relationship with respiratory effects, based on evidence for the development of asthma (U.S. EPA, 2016). And NO 2 will be included in future cycles of the Global Burden of Disease as global exposure estimates and evidence of their role as independent risk factors accumulate (Larkin et al., 2017). Additionally, atmospheric processing of NO 2 leads to the formation of nitrogen-bearing particles that can eventually deposit to the surface, causing acidification, nitrogen enrichment and other ecological effects (Russell et al., 2012). Local or global NO 2 monitoring is essential for understanding atmospheric chemistry as well as for humanhealth and environmental management and control.
Measurements of the intensity of ultraviolet or visible absorption spectra from the ground or from satellites are commonly used to retrieve the column density of NO 2 Valks et al., 2011;Berg et al., 2012). Satellitebased instruments such as the Ozone Monitoring Instrument (OMI), Global Ozone Monitoring Experiment (GOME and GOME-2) and SCanning Imaging Absorption SpectroMeter for Atmospheric CHartographY (SCIAMACHY) can provide global-scale NO 2 column measurements during daytime Bucsela et al., 2008). Moreover, plumes of NO 2 by cities, power plants and even ships have been able to be tracked using the recent high-spatialresolution observations of NO 2 from the TROPOspheric Monitoring Instrument (TROPOMI) on Sentinel-5P since 2017 (Lorente, et al., 2019;Georgoulias et al., 2020). However, they are unable to obtain local high-temporal-resolution NO 2 emissions such as variations in hourly NO 2 concentrations due to their long repeat cycle since the lifetime of tropospheric NO 2 is only about 6 h in summer and 18-24 h in winter due to the photochemical effect (Beirle, et al., 2003;Cui et al., 2016). In addition, measurements of tropospheric NO 2 from satellites or aircraft are also influenced and limited by clouds (Bovensmann et al., 1999;Liang et al., 2017). Ground-based measurements of column NO 2 from instruments such as Pandora using differential optical absorption spectroscopy (DOAS) are often used for the validation of satellite instruments (Herman et al., 2009;Lamsal et al., 2014;Kollonige et al., 2018). In situ measurements of near-surface NO 2 can best monitor local emissions. However, at this point in time, they cannot provide vertically resolved measurements. Balloon measurements using a NO 2 sonde can produce vertical profiles, but these measurements are very limited in time and space, especially in the Southern Hemisphere. The primary source of data on the vertical distribution of NO 2 comes from operational sites around the world. However, their operation can be expensive and laborintensive (Scott et al., 1999;Herman et al., 2009;Sluis et al., 2010).
The differential absorption lidar (DIAL) technique offers the potential for autonomous, 24 × 7 operation, with improved temporal resolution. Absorption of light by molecules is the basis for DIAL and numerous atmospheric constituents absorbing light. Conventional DIAL operates at two absorbing wavelengths, one stronger than the other, indicated by the on (λ on ) and off (λ off ) wavelength of the gaseous absorption feature of interest. Because of different absorption at λ on and λ off , the difference between the backscattered laser signals at the two wavelengths can be used to derive the number density of the absorption gas. Taking the log-ratio of these returns at closely spaced wavelengths removes system parameters and attenuation to and from the target of interest (Rothe et al., 1974;Sullivan et al., 2014). Thus, this technology provides measurements of the concentration of gas, such as NO 2 , O 3 and SO 2 , at a particular location and time (Fredriksson and Hertz, 1984;Newchurch et al., 2003;Volten et al., 2009;Kuang et al., 2013;Sullivan et al., 2017). The DIAL technique provides the unique capability of remotely monitoring localized NO 2 concentrations and emissions in urban and rural areas and profiling their tropospheric vertical NO 2 concentration. However, aerosols are abundant within the planetary boundary layer (PBL) and can cause significant retrieval errors in a two-wavelength DIAL technique to measure NO 2 . To better understand this aerosol problem and produce a more accurate NO 2 profile measurement, we described a new technique using a three-wavelength DIAL technique based on the intrinsic capabilities of using a multi-wavelength optical parametric oscillator (OPO) laser system. Hampton University (HU) has incorporated an OPO laser into its lidar system. The OPO laser enables researchers to optimize (tune) wavelength choices for specific measurements (Weibring et al., 2003). The three-wavelength DIAL retrieval equations are derived in this study. Our optimum choices for the three wavelengths to be used for our NO 2 retrievals are designed to help increase the difference in NO 2 absorption cross-section and reduce aerosol influence. NO 2 retrieval relative uncertainties are calculated using the two-wavelength DIAL (438 and 439.5 nm) and the three-wavelength DIAL (438, 439.5 and 441 nm). Tropospheric NO 2 profiles were obtained by applying the proposed technique to HU OPO DIAL lidar. As a first-order assessment, the HU lidar results were compared with simulated data from the Weather Research and Forecasting Chemistry (WRF-Chem) air quality model.

Method
To minimize aerosols-interference in the retrievals of NO 2 , a three-wavelength DIAL technique was proposed with λ 1 < λ 2 < λ 3 . Table 1 shows expressions for the extinction and backscatter of molecules and aerosols for these three wavelengths. In Table 1, β m and β a are backscatter from molecules and aerosols for the wavelength of λ 2 ; α m and α a are the extinction of molecules and aerosols for the wavelength of λ 2 ; and e is the aerosol Ångström exponent, assumed to be equal for the three wavelengths because the three wavelengths are very close.
The three elastic lidar equations can be expressed as where X is the lidar signal; C 1 , C 2 and C 3 are lidar constants; the subscripts a and m represent aerosol and molecule, respectively; σ N is the absorption cross-section for the gas of interest; N N is the molecular density of the gas of interest; O abs is absorption of gases other than the gas of interest; and z is the altitude. The molecular density of the gas of interest can be obtained using Eqs. (1), (2) and (3). The NO 2 density retrieval equation can be expressed as where AED, MED, OAD and B are the correction terms of aerosol extinction, molecular extinction, absorption of gases other than the gas of interest and backscattering, respectively. Because the atmospheric molecular density is relatively stable, MED can be corrected using a numerical model or local real-time radiosonde data. OAD can be removed by choosing appropriate wavelengths. However, aerosol is variable, especially in the PBL. For correction of AED and B, we need accurate aerosol measurements. However, accurate aerosol measurements are not easily obtained. From the above NO 2 retrieval relative equation, all of the correction terms are related to the three wavelengths, so how to choose the three wavelengths is very critical to reduce correction terms and improve the accuracy of NO 2 retrievals. We designed two rules to obtain the optimum choice for the three wavelengths: a. The three chosen wavelengths increase differences in the NO 2 absorption cross-section ( σ N ) to improve NO 2 retrieval.
According to Eq. (4), the greater σ N is, the less all of the correction terms are. So the three chosen wavelengths should help to increase σ N . Generally, researchers only used an increasing absorption method Fig. 1) to choose the three wavelengths (Wang, et al., 1997;Liu, et al., 2017). Wang used three wavelengths corresponding to the strong, medium and weak absorption of O 3 to obtain an accurate stratospheric ozone profile in the presence of volcanic aerosols. Liu used three wavelengths of 448.10, 447.20 and 446.60 nm corresponding to the strong, medium and weak absorption of NO 2 to retrieve NO 2 . Equations (10) and (11) are calculated values of σ N for the increasing absorption method and the decreasing absorption method using Eq. (5). Using the increasing absorption method and the decreasing absorption method to choose the three wavelengths, the values of σ N are both decreased according to Eqs. (10) and (11) compared to the conventional two-wavelength DIAL technique (the conventional two-wavelength DIAL technique retrieval equation is in the Supplement). According to characteristics of the NO 2 absorption spectrum shown in Fig. 2, a bumping absorption method (σ N (λ 1 ) < σ N (λ 2 ) & σ N (λ 3 ) < σ N (λ 2 )) is designed to choose the three wavelengths which can increase the value of σ N compared to the two-wavelength DIAL technique according to Eq. (12). However, for DIAL systems to measure other atmospheric gases like ozone, it is only practical to use wavelength selection Method B because of the shape of the ozone absorption spectrum (lacking narrow peaks). Table 1. Extinction and backscatter of molecule and aerosol for wavelengths of λ 1 , λ 2 and λ 3 .
Wavelength Molecular Aerosol Molecular Aerosol backscattering backscattering extinction extinction Increasing absorption method: Decreasing absorption method: Bumping absorption method: b. The three chosen wavelengths can reduce or remove AED.
It means the value of AED is equal or close to 0. Choosing the appropriate three wavelengths to make the value of K in Eq. (12) equal or close to 0, the value of AED will be equal or close to 0. The value of K in Eq. (12) changes with different aerosol Ångström exponents. For example, to remove boundary layer aerosol influence, we can set aerosol Ångström exponents to 1 to calculate the value of K to choose the three wavelengths because the size of aerosol in the boundary layer is typically large (Schuster et al., 2006;Li et al., 2017).

HU three-wavelength OPO DIAL system
The HU lidar is located on the campus of HU (  laser have recently been incorporated into the HU lidar system. The OPO laser enables researchers to optimize (tune) the wavelength choices and provides more flexibility than fixed-frequency wavelength shifters such as Raman cells.
The wavelength tuning range of our OPO extends from 192 to 2750 nm. This range is fully automated with precision scanning for true hands-free operation. Figure 3a and b show the Continuum Horizon II output energy and its parameters. The OPO laser energy outputs between 400 and 500 nm, which overlaps with the NO 2 strong-absorption spectral zone in Fig. 2 produced near the maximum possible power in the spectrum. Combining the OPO laser energy outputs, NO 2 absorption spectrum and two chosen three-wavelength rules, 438, 439.5 and 441 nm shown in Fig. 2 result in the wavelengths of the HU three-wavelength DIAL system because σ N of the three-wavelength pair is more than other threewavelength pairs in NO 2 strong-absorption spectral zone, and the K value of the three-wavelength DIAL system is 0.000023 (close to 0). The HU lidar system currently consists of a Continuum OPO laser system as the light source, a 48 in.
(1.22 m) non-coaxial Cassegrainian-configured telescope receiver, a light separation system that uses beam splitters and interference filters, a detecting system including photomultiplier tubes (PMTs) and avalanche photodiodes (APDs), and a Licel optical transient recorder. A schematic of the lidar system is shown in Fig. 4. The system can be configured to measure multi-wavelength aerosols and NO 2 density. High- resolution backscatter measurements extend from the boundary layer (1.2 km) to the free troposphere. The pump laser operates at three fixed wavelengths (1064, 532 and 354.7 nm). The 354.7 nm laser is mostly reflected into the OPO laser to produce three wavelengths (338, 339.5 and 441 nm). Steering mirrors whose axes are aligned with a receiving telescope axis direct these laser outputs into the atmosphere. The laser backscatter is collected by a 48 in. (1.22 m) diameter telescope and split into specific wavelength bands by a beam separation unit, which combines filters and beam splitters for dispersion of the return backscatter to various detection channels. Using filters and beam splitters makes the beamsplitting system simple, compact, and easy to change or add other spectral channels for other measurements. Currently, wavelengths of 438, 439.5, 441, 354.7, 532 and 1064 nm are focused on PMTs and APDs and recorded by a Licel datacollecting system for measurements of aerosol and NO 2 .
To demonstrate that the HU three-wavelength OPO DIAL system can effectively reduce aerosol influence and accurately retrieve NO 2 , retrieval correction terms of AED, MED, OAD and B in Eq. (4) are simulated using the twowavelength DIAL technique (438 and 439.5 nm) and the three-wavelength DIAL technique (438, 439.5 and 441 nm). Ozone was used for the simulation of OAD because only ozone absorption can produce a little influence on NO 2 retrieval based on the HITRAN 1.1.2.0 database. Atmospheric data of aerosol, molecules, O 3 and NO 2 for these simulations are from the HU local lidar aerosol measurements, radiosonde, NASA Tropospheric Ozone Lidar Network (TOL-Net), and NASA Deriving Information on Surface Conditions from COlumn and VERtically Resolved Observations Relevant to Air Quality (DISCOVER-AQ) measurements shown in Fig. 5. Extinction and backscatter of aerosol at 438, 439.5 and 441 nm can be calculated from the aerosol extinction profile at 532 nm in Fig. 5a with the setting of the lidar ratio = 50 and e = 1, 2 and 3. The lidar ratio is wavelengthdependent, and its value in the visible band is in general smaller than in the UV band for the same type of aerosols (Kuang et al., 2020;Reid et al., 2017). Absorption of NO 2 and O 3 at 438, 439.5 and 441 nm can be calculated using their mixing ratio profiles in Fig. 5b and their absorption crosssections from the HITRAN 1.1.2.0 database. MED, AED, OAD, B and absorption difference in NO 2 (NAD) are simulated using the two-wavelength DIAL technique with different aerosol Ångström exponents (e = 1, 2 and 3) shown in Fig. 6a, c and e and the three-wavelength DIAL technique shown in Fig. 6b, d and f. In Fig. 6, red lines are NAD, black lines are MED, deep-blue lines are AED, and light-blue lines are OAD. In Fig. 6, all OAD is far less than NAD. It is concluded that ozone absorption has a negligible influence on the retrieval of NO 2 . In Fig. 6a, c and d, MED and AED in the PBL are both more than NAD using the two-wavelength DIAL technique. Because atmospheric molecules are relatively stable, MED can be corrected using local model or real-time radiosonde data. However, aerosol is variable, so aerosols are a significant uncertainty for retrieving NO 2 with the conventional two-wavelength DIAL technique. In Fig. 6b, d and f, MED and AED in the boundary layer are both much smaller than NAD using the proposed three-wavelength DIAL technique. It is proven that the three-wavelength DIAL technique can effectively decrease retrieval errors caused by aerosol extinction. From Fig. 5, we can see that AED using the three-wavelength DIAL technique can be reduced to less than 2 % of AED using the twowavelength DIAL technique at least. Therefore, even if AED is not corrected, NO 2 can still be accurately retrieved. Moreover, simulated B using the two-wavelength DIAL technique and the three-wavelength DIAL technique is shown in Fig. 6 with green lines. The sharp change in vertical adjacent aerosol backscatter can cause drastic changes in the B term. In Fig. 6, the value of the B term using the three-wavelength DIAL technique is far less than using the two-wavelength DIAL technique. So the three-wavelength DIAL technique can reduce more fluctuation caused by aerosol backscattering than the two-wavelength DIAL technique.

Uncertainty analysis
According to NO 2 retrieval Eq. (4), the NO 2 measurement uncertainty is due to molecules, absorption of gases other than the gas of interest, aerosol and noise of lidar signals. The total relative uncertainty can be expressed as Eq. (13) (Leblanc et al., 2016).
where U NO 2 is NO 2 total retrieval relative uncertainty using the three-wavelength DIAL technique; U MED , U OAD , U AED , U B and U S are NO 2 retrieval relative uncertainty caused by molecules, absorption of gases other than the gas of interest, aerosol (extinction and backscattering) and noise of lidar signals expressed as Eqs. (14), (15), (16), (17) and (18); u is the uncertainty function; N m and N o are the number density   (ND) of air and ozone; σ m is the Rayleigh scattering crosssection; σ o is the absorption cross-section of ozone; and S is the lidar ratio. From Eqs. (14) and (15) In our measurements, profiles of temperature and pressure from local radiosondes are used to calculate N m . Usually, one radiosonde is launched for about 8 h of measurement. One profile of air number density from local radiosondes is used to correct 8 h NO 2 measurements. According to statistics of 8 h variation in temperature and pressure in four local seasons, the uncertainty in N m is between 1 % and 3 %. U MED using the two-wavelength DIAL technique and the three-wavelength DIAL technique is calculated according to Eq. (14), with the uncertainty in N a as 1 %, 2 % and 3 %, as shown in Fig. 7a. U MED using the three-wavelength DIAL technique is far less than using the two-wavelength DIAL technique. N o is obtained from local measurements. Because of very low values of ozone absorption cross-section differentials, with the uncertainty in N o as 50 % and 100 %, U OAD using the two-wavelength DIAL technique and U OAD using the three-wavelength DIAL technique are both less than 0.5 %, as seen in Fig. 7b. Ozone absorption correction is neglected in NO 2 retrieval. From Eqs. (16) and (17), U AED and U B are determined by uncertainties in a a , β a and e. For the HU lidar system, 532 nm elastic signals are used to calculate a a and β a with the Fernald method to correct NO 2 retrieval Figure 12. NO 2 number density total relative uncertainty with e = 1 (a), e = 2 (b) and e = 3 (c). (Fernald et al., 1972). A value of 50 sr is usually chosen as a lidar ratio to retrieve a a and β a . The lidar ratio is variable, so uncertainties in a a and β a are caused by the chosen lidar ratio. The range of the lidar ratio is about 30 to 70 sr for 532 nm. The uncertainty in the lidar ratio is 40 % for 50 sr. The uncertainties in a a and β a are calculated with uncertainty in the lidar ratio as 40 %. Finally, U AED and U B using the twowavelength DIAL technique and using the three-wavelength DIAL technique are calculated with the Ångström exponent as 1, 2 and 3, as shown in Figs. 8,9 and 10. From these figures, U AED and U B using the three-wavelength DIAL technique are both less than 4 %. However, U AED below 2 km using the two-wavelength DIAL technique is more than 90 % after correction of aerosol extinction. From Eq. (18), U S is determined by uncertainties in three-wavelength lidar signals. The uncertainties in lidar signals with an average integration time of 1 and 2 min are derived from Poisson statistics associated with the probability of detection of a repeated random event (Megie, 1985;Leblanc et al., 2016). NO 2 number density relative uncertainty owing to the noise of lidar signals with an average integration time of 1 and 2 min is obtained, as shown in Fig. 11. We can see that U S using the two-wavelength DIAL technique is smaller than using the three-wavelength DIAL technique. With an increase in average integration time from 1 to 2 min, U S can be effectively reduced. At last, U NO 2 (the total relative uncertainties in NO 2 ) with e as 1, 2 and 3 is calculated, as shown in Fig. 12a, b and c.

Results
The three-wavelength DIAL technique was implemented by the HU lidar measurements during two cases at night, and the resulting vertical profiles are presented in Fig. 13. All NO 2 lidar measurements presented here are obtained at times with less than 10 % cloud coverage below 8 km. HU lidar 438 nm (blue line), 439.5 nm (red line) and 441 nm (black line) elastic signals measured at 21:00 LT on 13 May 2020 and 22:00 LT on 27 July 2020 are shown in Fig. 13a and c, respectively. The average integration time for these signals is 2 min. Determined from the lidar elastic signals in Fig. 13a and c, there is an existing aerosol layer between 2.2 and 3.5 km on 13 May, while 27 July presented a clean atmosphere. Figure 13b and d show retrieved NO 2 profiles using the three-wavelength DIAL technique (red line). The black error bars in Fig. 13b and d indicate the uncertainty Figure 13. HU lidar 438, 439.5 and 441 nm elastic signals measured at 21:00 LT on 13 May 2020 (a) and 22:00 LT on 27 July 2020 (c); NO 2 profiles obtained using the three-wavelength DIAL technique, two-wavelength DIAL technique and WRF-Chem model at 21:00 LT on 13 May 2020 (b) and 22:00 LT on 27 July 2020 (d).
in NO 2 retrieval calculated using Eq. (13). In Fig. 13b, the retrieved NO 2 profile between 2.2 and 3.5 km on 13 May is smooth and not affected by the aerosol layer. The NO 2 profiles (sky-blue line and purple line) were also retrieved using the conventional two-wavelength DIAL technique without and with aerosol correction, as shown in Fig. 13b, resulting in a bump between 2.2 and 3.5 km in the NO 2 profile retrieved using the two-wavelength DIAL technique. This inconsistency suggests that the two-wavelength DIAL technique cannot remove AED of the aerosol layer between 2.2 and 3.5 km, and the retrieved NO 2 profile contains AED interference. Moreover, the NO 2 retrievals below 2 km using the two-wavelength DIAL technique shown in Fig. 13b and d are more than the three-wavelength DIAL technique, suggesting that the AED of boundary aerosol was not correctly removed. Aerosol correction is very important for NO 2 retrieval using the conventional two-wavelength DIAL technique (Sasano et al., 1985). These results suggest that the proposed three-wavelength DIAL technique can effectively remove the influence of aerosol on the retrieval of NO 2 . As a first-order assessment of the HU lidar NO 2 profiles, we compare the retrieval results to simulated data from the Weather Research and Forecasting Chemistry (WRF-Chem) model (Grell et al., 2005) at 12 km × 12 km spatial resolu-tion and 200 m vertical resolution. Past studies have demonstrated that WRF-Chem-simulated NO 2 results show good agreement between the OMI satellite measurements and aircraft measurements (Amnuaylojaroen et al., 2019;Barten et al., 2020), providing a data source to examine the accuracy of the HU retrievals using both the two-wavelength DIAL technique and the three-wavelength DIAL technique. The HU local NO 2 profiles for these two cases are simulated using the WRF-Chem model and shown in Fig. 13b and d. WRF-Chem-simulated NO 2 magnitudes tend to be lower compared to HU-retrieved NO 2 profiles using the three-wavelength DIAL technique (typically within ±0.1 ppb), except above 3.5 km on 13 May 2020; however, the comparison demonstrates a consistent vertical profile shape between observations and the model simulation. And retrieval results using the three-wavelength DIAL technique are much closer to simulated values compared to using the two-wavelength DIAL technique. These figures also demonstrate that the reduced fluctuations caused by aerosol backscatter when using the three-wavelength DIAL technique result in vertical profiles of NO 2 which are much more consistent with simulated data when compared to results of the two-wavelength DIAL retrievals. Both the WRF-Chem-simulated profiles and the HU retrievals of NO 2 using the three-wavelength DIAL technique are associated with uncertainties which could result in the differences in magnitude; however, the consistent nature in the vertical profile shapes from both data sources provides confidence that the HU lidar is retrieving NO 2 vertical profiles using the three-wavelength DIAL technique in the troposphere.

Conclusion
This study describes a lidar retrieval technique using three wavelengths simultaneously emitted from an OPO laser to measure tropospheric NO 2 profiles. The three-wavelength DIAL retrieval equations describe how the retrievals decrease errors caused by aerosol interference. Aerosol extinction differences using this proposed technique can be decreased to less than 2 % of the aerosol extinction differences resulting from a conventional two-wavelength DIAL technique. Comparing the HU lidar results to WRF-Chem model output demonstrates that the NO 2 magnitudes and vertical structure are in much better agreement with simulated data when applying the three-wavelength DIAL technique compared to using the two-wavelength technique. In the future, we will add new filters to obtain daytime NO 2 measurements. We also plan to purchase NO 2 balloon sondes for acquiring true validation data to evaluate HU lidar NO 2 results.
Code availability. The software code for this paper is available from the first author.
Data availability. The HU lidar data are archived at http://cas. hamptonu.edu/data-products/ (Su, 2020). DISCOVER-AQ data and WRF-Chem model output data are available from the first authors upon request.
Author contributions. JS designed the three-wavelength OPO DIAL system, did simulation work for the system, developed the NO 2 retrieval algorithm and prepared the original manuscript. PMM was responsible for funding acquisition. MSJ provided simulated data from the WRF-Chem model to verify NO 2 retrieval in this study. JTS, MJN, TAB, SK and GPG contributed to the analysis of NO 2 retrieval uncertainty. All listed authors contributed to the review and editing of this paper.