Development of a passive doas system to retrieve atmospheric pollution columns in the 200 to 355 nm region
© Galicia et al.; licensee BioMed Central Ltd. 2013
Received: 13 December 2012
Accepted: 15 December 2012
Published: 8 January 2013
In recent years several techniques have been developed to measure and monitor the pollution of the air. Among these techniques, remote sensing using optical methods stands out due to several advantages for air quality control applications. A Passive Differential Optical Absorption Spectroscopy system that uses the ultraviolet region from 200 to 355 nm of the solar radiation is presented. The developed system is portable; therefore it is practical for real time and in situ measurements. The enhanced wavelength range of the system is intended to detect the ultraviolet light penetration in the Mexican Valley considering the solar zenith angle and the altitude. The system was applied to retrieve atmospheric SO2 columns emitted either by anthropogenic (power plant) or natural sources (volcano), reaching a detection limit of about 1 ppm. The measurement of the penetrating solar radiation on the earth surface at the UVC range is presented and the possibility to measure pollution traces of some contaminants as O3, NO2 and aromatic compounds in real time and in situ in the ultraviolet region is discussed.
The production of thousands of chemicals has contributed to industrial and economic development in many parts of the world. This trend however has been associated with the release of new chemicals and possibly toxic substances into the environment and food chain and adversely affecting human health in many instances. The pollutants associated with the anthropogenic activities could be inorganic as well organic compounds . The air pollution is very complex dynamic phenomena. It could be associated with the industrial activities as well as with the concentration of the population in large towns and cities .
The power plants are using resources like fuel and water to provide electricity that is one of the essential needs for sustainable development and life. This activity produces and discharges all different kinds of pollutants such as, gaseous, liquid, electromagnetic fields, and noise which endanger our lives and environments. The amounts of pollutants discharged by power plants, especially air pollution are more than assimilation capacity of nature because now it is clear that sustaining and assimilative capacity of the biosphere though tremendous, is after all finite. The first step is to understand the magnitude of emissions from each source .
The Differential Optical Absorption Spectroscopy (DOAS) has been successfully used for air pollution detection in several research works . Important advantages over other techniques like gas chromatography, chemiluminescence, and electrochemical, gravimetric, and matrix isolation, are the lower cost and the time saving, especially when the sunlight is used as the source (passive DOAS). This technique is particularly interesting in the monitoring of pollution emissions by refineries or thermoelectric power stations.
Due to the absorption by the O3 layer above the earth in the UVC interval (200 to 290 nm), the DOAS technique has been based on the UVB-VIS range absorption measurements. It has been used to detect chemicals such as CO, SO2, O3, NO2 HNO2, HCHO, CH3, CHO, benzene, toluene, p- and m-Xylene .
The absorption lines of these pollutants interfere with eachother making more laborious the processing of the measured spectra. Usually, the measurement of SO2 columns by the passive DOAS technique is done using the 290 to 330 nm wavelength range, which is well-known penetrates at sea level. In this region the SO2 absorption lines interfere with the absorption lines of less common species as HCHO, BrO, O3 and ClO .
In the wavelength range from 230 to 260 nm the ozone is the principal absorber. Benzene and toluene have cross sections one order of magnitude lower than O3. Ozone atmospheric measurements have been previously reported using active DOAS in the range from 265 to 345 nm [6, 7]. Also, atmospheric O3 measurements using passive DOAS have been done in the range from 438.5 to 540 nm. However due to the overlapping of its absorption spectrum with the NO2 absorption spectrum the concentration of both elements is measured .
It is very well documented that the solar radiation in the 190–280 nm range is strongly absorbed by the stratospheric ozone. Hitherto every report states that almost all incoming solar UVC and 90% of UVB are absorbed by stratospheric ozone . As a result, there are no published experiments measuring the absorption levels in the UVC range. However, if there is any irradiance left, once the radiation passes to the troposphere the O3 concentration decreases and it would be possible to detect the remaining radiation in UVC with a sufficiently sensitive detector.
The solar irradiance at the earth’s surface varies greatly depending on factors such as latitude, time of the day, month of the year, cloud cover, and haze (aerosols). The spectral irradiance in the UVC range on the earth surface also depends on the ozone column, on the altitude above sea level and the solar zenith angle. Very sensitive measurements from 280 to 300 nm at sea level, 33° North latitude and 25.14° Solar Zenith Angle (SZA), using PMT detectors, reported the detection of solar spectral irradiance of 3x10-4 Wm-2 nm-1. Similar measurements have reported a value of two orders of magnitude higher at 2500 m altitude, 220 North latitude and 7.8° SZA .
The Mexican Valley is a high plateau in central Mexico with a minimum altitude of 2200 m and latitudes between 19o and 20o North Latitude. Under such conditions the solar zenith angle is almost 0o at noon. By what, these values propitiate a higher probability to detect some irradiance at UVC wavelength range with very sensitive detectors as photomultiplier tubes used in spectroradiometers or charged coupled devices (CCD) used in minispectrometers . In Mexico City the O3 column amounts between 240 and 320 DU , if there is a possibility to detect some spectral irradiance in the UVC region it could be possible to detect traces absorbing in this wavelength interval. So, in order to prove the usefulness of the passive DOAS technique in the range from 200 to 350 nm in the Mexican Valley, we have developed a portable system to measure pollution slant columns in this spectral range of several pollutants as SO2, O3, NO2, aromatic compounds, etc.
In this issue the SO2 emissions around the Tula's Industrial Complex in the Mexican Valley and around the Popocatepetl Volcano are reported. These calculations were performed using the developed system. Some discussion is presented regarding the observed well-structured absorption spectra in the range of 200 to 290 nm.
Materials and methods
The system consists of a fused silica lenses telescope 84-UV-25 (Ocean Optics) with a 25.4 mm diameter, focal length of 100 mm. This telescope is especially suitable for light collimation from a distant source. The collimated light is guided to the spectrometer using an optical fiber (QP600-2-SR, Ocean Optics) with a 600 m diameter core and 2 m long. The spectra measurement was performed using a miniature spectrometer (USB 2000, Ocean Optics) with a diffraction grating of 2400 lines/mm. A cylindrical lens is incorporated to the detector to increase light collection efficiency focusing a 1 mm high slit onto a considerably shorter (200 μm) detector elements. The CCD detector (ILX511A, Sony) is coated with an UV anti-reflecting film to reduce the reflection losses and therefore increasing their sensitivity. The spectral range of the spectrometer is in the 190 to 355 nm range with a resolution of 0.5 nm. This system shows a high signal to dark current random noise ratio in the 190 to 290 nm range (over 3 times for measurements under pollution measurements and about 12 times under clean sky measurements). The system was checked and no overlapping of different diffraction orders and the scattering of the incident light which could affect the readings in the UVC region was detected. To reduce the background noise of the measurement an averaging of 3 to 5 spectra was performed.
The coordinates of the measuring spots were determined with a GPS locator (GT 37231). All the information inputs a laptop were the developed software processes the data. The coordinates given by the GPS locator were superimposed on the Google satellite images using the Google-Earth-KH-80.llb libraries to show the actual measuring points on a satellite photograph.
The attenuation of the irradiance I0(λ) is described by the Beer-Lambert’s Law I(λ) = I0(λ)e− σ(λ) ∫ c(s)ds. The attenuation depends of the concentration c(s) of the different compounds of the light path media . The absorption cross section σ(λ) defines the light absorption to a specific wavelength. The dependence of the cross sections σ(λ)on the temperature and pressure has not been considered in this work. The gases columns are evaluated by using the mixing ratio method . The light absorption cross section of a specific pollutant can be obtained from reported data on literature or from laboratory measurements.
Usually, it is very difficult to obtain separately the pollutant concentrations c(s) and the optical length l=∫ds. Therefore, these parameters are combined to form what is called slant column density S=∫c(s)ds expressed in molecules/cm2[16, 17].
Actually, the measured absorption involves all the present molecules that absorb in that specific wavelength. The absorption of several substances i will be then described by the following expression: , where σ i (λ), c i (S) and S i are their absorption cross section, its specific concentration and its slant column density, respectively. The sunlight spectrum that reaches the Earth surface shows slow and fast variations. The slow variations of the spectrum are determined by the light source and by the scattering. The fast variations are determined by the absorption of molecules. The spectrum of the cross sections can be divided in a slow variations σ i S (λ) and fast variations parts σ i F (λ). Therefore, the absorption cross section for each component will be σ i (λ) = σ i S (λ) + σ i F (λ).
The slow variations part of the measured light intensity, which depends on I0(λ), can be approximated by a polynomial. Using some polynomial adjustment method or by using high-pass digital filters, for instance the Savitzky-Golay  a solution for fast variations part σ i F (λ) could be obtained [4, 19].
Where D(λ) is the Optical Density of the medium [4, 15]. The cross sections of the gas trace at the atmosphere could be considered independent from their placement; this means that they are independent from their altitude .
High resolution absorption spectrum for the SO2 has been reported . The reference spectrum for the evaluation of the measurements is needed and it must consider only the fast variations σ i F (λ) and the spectral deformations produced by the spectrometer and determined by the spectrometer transfer function (STF). In practice, to obtain the STF the spectral lines emitted by a Hg lamp could be used because, these are narrower (~10 pm)  than the typical resolution of the used spectrometer (0.1–1 nm). The measured spectrum is a good approximation of the spectrometer transference function . Mathematically the measured spectrum I*(λ) can be expressed as the convolution of the Hg lamp emission spectrum I(λ') with the spectrometer transference function I*(λ) = I(λ′) H sp (λ′), where H sp (λ') is the spectrometer transfer function. The STF was measured at 334.1 nm because at this wavelength the “cleanest” signal is obtained.
A 2.5 mbar, 3 cm long, SO2 cell was used to validate our DOAS system. We used a Xe-lamp as radiation source. The measured column agrees within a +− 4% error. To determine the density of a measured column a non-linear adjustment is done by using the Levenberg-Marquardt method  using the measured I0(λ) and the reference spectra I0(λ). A dark spectrum, previously captured, is subtracted from I(λ) and I0(λ). A residual Δψ is obtained from the fitting process, in our case in the order of 2.5x1017 particles/cm2, which can be used to calculate the detection limit of our system trough the relation ; hence the detection limit is the order of 1 ppm.
The spectrometer is connected to a laptop. A program was developed in Labview to capture and filter the measured spectra to perform the non-linear adjustment and to identify the amount of polluting gas. This software consists of several programming stages. First of all, a driver for the Windows operative system is installed from the spectrometer that comes from the OOIDDrv32.dll (Ocean Optics, 2011) library. Once the hardware is recognized a Sarvitzky-Golay filter program filters the reference I0(λ) and measured I(λ) spectra. A dark spectrum is subtracted from and I(λ) to reduce the thermal noise of the instrument itself. Immediately, the relation of these spectra is calculated and is applied the logarithm function. Spectra are geographically referenced throughout the journey to integrate all the SO2 measured columns. A mouse-type GPS was used with the NMEA protocol to link it with the program, to do so, an API for Windows is used which is provided by Labview under the VISA name, which allows to read the data from the GPS (; Visa, ).
When the DOAS technique is used, the path length of light depends on the number of scattering events that take place until the photon reaches the spectrometer. Some methods to obtain at least the statistics of photon propagation are based on the simulation Monte Carlo method  to be able to validate statistically what happens with a large amount of photons it is necessary to simulate such events unfortunately this process requires a computing robust system and the simulation process is not performed in real time. That is not practical for insitu measurements. Taking this into consideration, in this work the column densities are indicated in ppm*m, using mixing ratios [27, 28]. The SO2 flow is calculated in tones by day , using column densities measured by the system, which must be multiplied by the distance of the traverse route using the following formula: Here X is the traverse distance in meters from the point x 1 to x n , V is the wind speed in m/s, θ is the angle between vectors X and V. For the real time in situ measurements around the industrial zone and the volcano, the telescope was fixed on a car for the traverse. Five spectra were averaged in each data entry to minimize the noise effects. The spectra are taken every 3 seconds. The system takes firstly the dark spectrum to reduce the noise produced by the spectrometer and the skylight spectrum I skl (λ) is measured. Finally, the sample spectra are taken during the traverse around the industrial zone.
The reference and measured spectra used for the determination of SO2 content in the 280–320 nm wavelength range are shown in Figure 1. The reference spectrum was obtained according to the method described earlier. The absorption bands of SO2 are very well defined in this interval.
The developed system was used to measure both, anthropogenic and natural SO2 emissions. Two sites with significant SO2 concentration were chosen near Mexico City.
Measurements at Tula hidalgo refinery
Measurements at Popocatepetl volcano
At some point of the traverse, indicated with an arrow on Figure 6A, the direction changes and it becomes nearly perpendicular to the volcano’s direction. In this second part of the traverse the measured values of SO2 columns remain practically constant.
Usually, the measurement of SO2 columns by the passive DOAS technique is done using the 290 to 330 nm wavelength range. In the wavelength range from 230 to 260 nm the ozone is the principal absorber. Albeit the arrangement of our spectrometer shows less resolution than a similar one with a range from 290 to 355 nm, from Figures 5 and 6 it can be clearly observed that the developed system has still enough resolution to determine the SO2 pollutant distribution and, after Figure 3 (for altitude and latitude conditions as in the Mexican Valley), also could be used to simultaneously measure O3 columns in the 230 to 260 nm range.
Still is necessary the measurement of the wind direction and speed in the plume at the same time with the column measurement to estimate the flow of the SO2; this is an open problem. If the wind conditions remain steady as the columns measurements along an appropriate traverse are performed, a good approximation of the emitted concentration and the pollutant flow could be calculated.
In this work we describe a passive DOAS system developed to determine SO2 polluting columns. Portability of the developed system makes it practical for real time measurements and it allows the SO2 columns measurements in situ. It could be applied to develop some models of pollutants distribution considering the source, the diffusion and the wind conditions. The direct measurement of the pollutant columns provides useful information to evaluate the pollution level in real time. Also, DOAS technique avoids contact with the samples providing a better interpretation of the SO2 columns measurements.
Attenuation in the plume of the reference spectrum I0(λ) in the wavelength interval from 200 to 290 nm is presented in Figure 2B. It should be emphasized that a certain level of irradiance in the UVC region is detected on the earth surface. It could be related to the altitude and latitude of Mexican Valley, as well as with the continuous improvement of optical detectors. The absorption in the UVC range could be related to the O3 concentration. It leads to a possibility to perform DOAS in this wavelength band applied to ozone detection and monitoring at Mexico City region, due to its altitude and latitude. However, much more research work should be done to corroborate this hypothesis.
This study was supported and founded by IPN-México. Rubén Galicia was supported by COTEPABE-IPN and COFFA department.
- Ebtekar M: Effects of persistent organic pollutants on the immune system: The case of dioxins. Iranian J Env Health Sci Eng. 2004, 1 (2): 1-7.Google Scholar
- Makhelouf A: The effect of green spaces on urban climate and pollution. Iranian J Env Health Sci Eng. 2009, 6 (1): 35-40.Google Scholar
- Ghiasseddin M: Comparison of Iran power plants air pollutants before and after shifting to natural gas. Iranian J Env Health Sci Eng. 2004, 1 (1): 39-41.Google Scholar
- Platt U, Stutz J: Differential Optical Absorption Spectroscopy, Principles and Applications. 2008, Springer, Berlin Heidelberg, 135-285.View ArticleGoogle Scholar
- Mellqvist J, Samuelsson J, Rivera C, Lefer B, Patel M: Measurements of industrial emissions of VOC`s, NH3, NO2 and SO2 in Texas using the solar occultation flux method and mobile DOAS. 2007, Final Report HARC proyect H-53: Available at http://files.harc.edu/Projects/AirQuality/Projects/H053.2005/H053FinalReport.pdfGoogle Scholar
- Axelsson H, Edner H, Galle B, Ragnarson P, Rudin M: Differential optical absorption spectroscopy (DOAS) measurements of ozone in the 280–290 nm wavelength region. Appl Spectrosc. 1990, 44: 1654-1658. 10.1366/0003702904417625.View ArticleGoogle Scholar
- Lee JS, Kuk BJ, Kim YJ: Development of a differential optical absorption spectroscopy (DOAS) system for the detection of atmospheric trace gas species; NO2, SO2, and O3. J Korean Physical Soc. 2002, 41: 693-698.Google Scholar
- Raponi MM, Jímenez R, Wolfram E, Tocho JO, Quel EJ: Remote sensing of stratospheric O3 and NO2 using a portable and compact DOAS spectrometer. J Physics: Conf Series. 2011, 274: 012094-Google Scholar
- Dalback A: Global monitoring of atmospheric ozone and solar UV radiation in: Solar Radiation and Human Health. Edited by: Bjertness E. 2008, The Norwegian Academy of Science and Letters, Oslo, 23-34.Google Scholar
- Gueymard CA, Meyers D, Emery K: Proposed reference irradiance spectra for solar energy systems testing. Solar Energy. 2002, 73: 443-467. 10.1016/S0038-092X(03)00005-7.View ArticleGoogle Scholar
- Pinedo JL, Mireles F, Rios C, Quirino LL, Davila JI: Spectral signature of ultraviolet solar irradiance in Zacatecas. Geofisical Int. 2006, 45: 263-269.Google Scholar
- Abbott TMC, Kleinman SJ: CCD'S and PMT's: A Simultaneous comparision. Interacting Binary Stars: a symposium held in conjution with the 105th Meeting of the Astronomical Society of the Pacific, San Diego State University, 13-15 July 1993. Edited by: Shafter AW. 1994, Astronomical Society of the Pacific, San Francisco, Calif, 407-4012. 56Google Scholar
- Juárez A, Bravo JL: Influence of urban zone in the measurements of the total ozone column in Mexico City. Atmosfera. 1995, 8: 35-43.Google Scholar
- Barrancos J, Rosello JI, Calvo D, Padrón E, Melián G, Hernández PA, Pérez NM, Millán MM, Galle B: SO2 emission from active volcanoes measured simultaneously by COSPEC and mini-DOAS. Pure Appl Geophys. 2008, 165: 115-133. 10.1007/s00024-007-0290-8.View ArticleGoogle Scholar
- Wang HS, Zhang YG, Wu SH, Lou XT, Zhang ZG, Qin YK: Using broadband absorption spectroscopy to measure concentration of sulfur dioxide. Appl Phys B. 2010, 100: 637-641. 10.1007/s00340-010-4151-2.View ArticleGoogle Scholar
- Johansson M, Rivera C, de Foy B, Lei W, Song J, Zhang Y, Galle B, Molina L: Mobile mini-DOAS measurement of the outflow of NO2 and HCHO from Mexico City. Atmos Chem Phys. 2009, 9: 5647-5653. 10.5194/acp-9-5647-2009.View ArticleGoogle Scholar
- Kraus SG: DOASIS, A Framework Design for DOAS.: Dissertation. 2006, Editorial Shaker-Verlag, Manheim Germany, 1-184. ISBN 3-8322-5452-8Google Scholar
- Barak P: Smoothing and differentiation by an adaptive-degree polynomial filter. Anal Chem. 1995, 67: 2758-2762. 10.1021/ac00113a006.View ArticleGoogle Scholar
- Jähne B: Digital Imagen Processing. 2002, Springer Verlag, Berlin, Heidelberg, 618, 5View ArticleGoogle Scholar
- Rozanov VV, Rozanov AV: Differential optical absorption spectroscopy (DOAS) and air mass factor concept for a multiply scattering vertically inhomogeneous medium: theoretical consideration. Atmos Meas Tech Discuss. 2010, 3: 697-784. 10.5194/amtd-3-697-2010.View ArticleGoogle Scholar
- Vandaele AC, Simon PC, Guilmot JM, Carleer M, Colin R: SO2 absorption cross section measurement in the UV using a Fourier transform spectrometer. J Geophys Res. 1994, 99 (D12): 25599-25605. 10.1029/94JD02187.View ArticleGoogle Scholar
- Sansonetti CJ, Salit ML, Reader J: Wavelengths of spectral lines in mercury lamps. Appl Opt. 1996, 35: 74-77. 10.1364/AO.35.000074.View ArticleGoogle Scholar
- Kraus S: DOASIS: DOAS Intelligent System. 2010, Available at http://www.iup.uni-heidelberg.de/bugtracker/projects/doasis/Google Scholar
- Manolis IA: A brief description of the Levenberg-Marquardt algorithm Implemened by levmar. Institute of Computer Science, Foundation for Research and Technology - Hellas (FORTH) Vassilika Vouton, P.O. Box 1385, GR 711 10. 2005, Available at http://www.ics.forth.gr/~laurakis/levmar/levmar.pdfGoogle Scholar
- van Diggelen F: Assited GPS, GNSS, and SBAS. 2009, Artech House Boston, London, 25-90.Google Scholar
- Visa. 2011, http://www.ni.com.visa.joule.ni.com/nidu/cds/view/p/id/3342/lang/es,
- McGonigle AJS, Aiuppa A, Giudice G, Tamburello G, Hodson AJ, Gurrieri S: Unmanned aerial vehicle measurements of volcanic carbon dioxide fluxes. Geophys Res Lett. 2008, 35: LO6303-View ArticleGoogle Scholar
- Rivera C, García JA, Galle B, Alonso L, Zhang Y, Johansson M, Matabuena M, Gangoiti G: Validation of optical remote sensing measurement strategies applied to industrial gas emissions. Int J Remote Sensing. 2010, 30: 3191-3204.View ArticleGoogle Scholar
- Sosa IG, Rivera C, Ortega L: Medición del flujo de emisiones totales de SO2 y NO2 generadas por fuentes urbanas mediante detección remota. Tecnol ciencia Ed (IMIQ). 2006, 21 (2): 103-109.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.