Intercomparison of Ozone Vertical Profile Measurements by Differential Absorption Lidar and IASI / MetOp Satellite in the Upper Troposphere – Lower Stratosphere

This paper introduces the technique of retrieving the profiles of vertical distribution of ozone considering temperature and aerosol correction in DIAL sounding of the atmosphere. The authors determine wavelengths, which are promising for measurements of ozone profiles in the upper troposphere–lower stratosphere. An ozone differential absorption lidar is designed for the measurements. The results of applying the developed technique to the retrieval of the vertical profiles of ozone considering temperature and aerosol correction in the altitude range 6–15 km in DIAL sounding of the atmosphere confirm the prospects of ozone sounding at selected wavelengths of 341 and 299 nm with the proposed lidar. The 2015 ozone profiles retrieved were compared with satellite IASI data and the Kruger model.


Introduction
Laser remote sounding techniques with the use of the lidar (Light Detection and Ranging) technology are widely used for the study and monitoring of the atmosphere.Among the techniques used for measuring the spatial distribution of the concentration of an atmospheric gas, the most sensitive is the differential absorption method (DIAL).The essence of the method is that radiation is transmitted in the atmosphere simultaneously at two wavelengths: one of them (λ on ) is on a strong absorption line or band of the gas to be measured and another (λ off ) is off the absorption line and is weakly absorbed in the atmosphere or not absorbed at all.At present, DIAL sounding of trace atmospheric gases including ozone, is implemented [1][2][3].
Laser sounding of the ozonosphere became routine at some observatories since the second half of the 1980s [4][5][6][7].It allows data to be received on vertical distribution of ozone (VDO), which successfully supplement similar data received by in situ methods with the use of ozonesondes, rockets, or satellites (TOMS, SAGE-II, TERRA, MetOp, etc.).
It should be noted that the profiles of IASI/MetOp-measured atmospheric ozone concentrations were earlier compared with the profiles measured in the troposphere with compact airborne lidars, ozone sondes, ground-based lidars, and Brewer-Dobson spectrophotometers [8][9][10].
Laser sounding of the vertical distribution of stratospheric ozone has been carried out at the Siberian Lidar Station (SLS) of Institute of Atmospheric Optics SB RAS (Tomsk, 56.5 • N, 85.0 • E) since 1989.A long period of lidar observations of stratospheric ozone has shown that the most important part of the ozonosphere for the study is located in the lower stratosphere, where ozone is affected by the dynamic factor.DIAL measurements of ozone profiles in different altitude ranges with different ozone content are carried out at different combinations of wavelengths [4][5][6][7].
The aims of this work are the development of an algorithm and program for VDO retrieval with temperature and aerosol correction; selection of the VDO sounding wavelengths; the design of an ozone lidar; and analysis of the measurement results, including the comparison between the ozone vertical profiles retrieved from the ground-based DIAL lidar system data and the profiles retrieved from IASI/MetOp data.

Selection of Wavelengths
Lidar measurements of VDO are carried out on the basis of the method of the differential absorption of the backscattered laser radiation in the ultraviolet spectral range 200-370 nm (Hartley-Huggins band) [11].In practice, several pairs of wavelengths can be implemented in ozone lidars with the help of different lasers.Table 1 represents the specifications of some lidars used for ozone measurements and operating at different combinations of wavelengths.More than 85% of atmospheric ozone is located in the stratosphere.To measure the much smaller concentrations of tropospheric ozone, sounding wavelengths should be selected from the short wavelength range, closer to the ozone absorption band center, to increase the concentration sensitivity of the method.In this spectral range, the absorption cross-section σ is several times larger than that for wavelengths used in stratospheric measurements (e.g., σ 299 = 4.4 × 10 −19 cm 2 for the wavelength λ on = 299 nm and σ 308 = 1.4 × 10 −19 cm 2 for λ on = 308 nm).
KrF laser radiation (248 nm) or the fourth harmonics of an Nd:YAG laser (266 nm) is usually used for measurements of tropospheric ozone, in combination with a technique based on the stimulated Raman scattering (SRS) in H 2 , D 2 , CO 2 , and other gases [4,5,7].The most common are hydrogen and deuterium.Table 2 represents possible sets of wavelengths that correspond to the first, second, and third Stokes (C) frequencies of SRS conversion in H 2 , D 2 , and CO 2 .Different wavelength combinations are used in practice in different altitude ranges in the troposphere and lower stratosphere.Thus, the wavelength pairs 289/316 and 287/299 nm allow measurements of ozone profiles up to altitudes of about 10 km [4,5]; the pair 292/319 nm, up to 14-16 km [4]; the pairs 277/313 and 292/313 nm, up to altitudes of 8-12 and 15 km, respectively [7].
We have estimated possibilities of the vertical ozone profile sounding in the upper troposphere-lower stratosphere at the wavelength pair 299/341 nm.During the calculations, actual lidar parameters were used: radiation energy of 20 mJ at both wavelengths, pulse repetition frequency of 15 Hz, receiving mirror diameter of 0.5 m, signal accumulation time of 1.5 h.To determine the transmitter-receiver efficiency, actual values of transmittance of optical elements of spectral selection and of efficiency of photomultiplier tubes (PMTs) were used; noises were taken from actual measurements.The calculations have shown that these wavelengths allow the sounding up to about 22 km (the ozone maximum is located in the altitude range 19-21 km in Tomsk) and ozone measurements in the troposphere.The ozone concentration error is within 4-10% limits up to altitudes of about 20 km.
The sounding altitude maximum is determined, first, by the range of signal detection at λ on , which is always shorter than the range of signal detection at λ off due to stronger ozone absorption.In view of this, λ on = 299 nm is preferable to 277 or 292 nm.In addition, wavelengths of 299 and 341 nm are implemented in one sounding beam (from one laser source in one SRS cell), in contrast to, e.g., the 292/313 nm pair (see Table 2).
A system on the basis of a SRS cell filled with hydrogen is cheaper than with deuterium.Thus, the wavelength pair 299/341 nm is of higher information content for VDO measurements in the upper troposphere-lowerstratosphere (5-22 km altitude range).

Theoretical Base of the VDO Retrieval
Initial equations for the calculation of the ozone concentration during DIAL lidar sounding of the atmosphere have the form where N(H) is the echo-signal recorded at corresponding wavelengths (on at an absorption line and off out of the absorption line), C is the instrumental constant, α a is the aerosol extinction coefficient, β a is the aerosol backscattering coefficient, k is the absorption cross-section of ozone, n(H) is the ozone concentration.
Let Equation (1) be divided by Equation (2): where k on (H) and k o f f (H) are the absorption coefficients and are off the ozone absorption line, dependent on the temperature.Let us transform Equation (3) to the form and derive Equation (4): The final equation for the ozone concentration is derived with the use of mathematical transformation or Actual variations in the air temperature can cause variations in the ozone absorption cross-section, which results in systematic errors in VDO retrieval.Therefore, it is reasonable to carry out a correction to the temperature dependence k on (H,T), k off (H,T) in the algorithm for VDO retrieval.
A model of the behavior of ozone absorption cross-sections, presented in Table 3 and based on data from [18,19], was used in the technique.
The absorption coefficients k on (H) and k o f f (H) are used in term A of Equation (6).Real temperature variations in the atmosphere can significantly change the ozone absorption coefficient.Therefore, it is reasonable to use the correction to the temperature dependence k on (H, T), k o f f (H, T) in the VDO retrieval algorithm such that A takes the form where T(H)-is the model Kelvin temperature distributed with altitude.Let us consider term C in Equation ( 6): and transform it to the form where, at the corresponding wavelengths λ (on at an absorption line and off out of the absorption line), R off (H) is the real distribution of the scattering ratio, x is the parameter that characterizes the particle size, β α off (H) is the aerosol backscattering coefficient, β m off (H) is the molecular backscattering coefficient.

From the equalities
Finally, Now consider term D in Equation ( 6).Since then Using Let us also transform term F in Equation ( 6) using the relations Aerosol scattering exceeds molecular scattering by several times in the case of high atmospheric aerosol content, which significantly distorts ozone profiles retrieved under unconsidered scattering and attenuating properties of the atmosphere at the sounding wavelengths.In the algorithm for VDO retrieval described, the aerosol correction is considered in the Equations ( 9) and ( 15) by means of introduction of a real distribution of the scattering ratio R off (H), while VDO in the usual, undisturbed atmosphere can be calculated at R off (H) = 1.
The vertical profile of the backscattering coefficient β a π (H) is derived from laser sounding data.The coefficient decreases with an increase in altitude.For clearer representation of aerosol stratification, the scattering ratio is used, where β M π (H) is the molecular backscattering coefficient.

Software for VDO Retrieval
Software for altitude ozone profile retrieval from laser sounding data developed on the basis of the above algorithm (see the block diagram in Figure 1) allows: 1.
Reading the lidar data; 2.
Recording the retrieval results in ASCII format; 3.
Moving average smoothing of lidar signals; 4.
Smoothing of the VDO retrieval results.To decrease the retrieval errors, the temperature correction of ozone absorption coefficients is used in the software.A high aerosol concentration in the 0-20 km altitude range should be considered during lidar signal retrieval when sounding at 272/289 nm and 299/341 nm wavelengths; therefore, aerosol correction is considered in the software.Model seasonal (winter and summer) midlatitude values of altitude distributions of the temperature and molecular backscattering coefficient have been introduced in the software for calculations.
Linear smoothing is used in the software for both input lidar data and retrieval results.The linear smoothing (moving average smoothing) is a well-known technique and is widely used in processing of experimental data in different fields of natural sciences.The linear smoothing is a special case of numerical filtering of signals with random errors using a rectangular window and unit weight coefficients.
The technique and software developed were used for VDO retrieval in the upper stratosphere-lower stratosphere at 299/341 nm wavelengths.

SLS Ozone Lidar
The sounding wavelengths selected (299/341 nm) were used in the SLS ozone lidar designed.Its block diagram is shown in Figure 2.
Remote Sens. 2017, 9, 447 7 of 15 midlatitude values of altitude distributions of the temperature and molecular backscattering coefficient have been introduced in the software for calculations.Linear smoothing is used in the software for both input lidar data and retrieval results.The linear smoothing (moving average smoothing) is a well-known technique and is widely used in processing of experimental data in different fields of natural sciences.The linear smoothing is a special case of numerical filtering of signals with random errors using a rectangular window and unit weight coefficients.
The technique and software developed were used for VDO retrieval in the upper stratosphere-lower stratosphere at 299/341 nm wavelengths.

SLS Ozone Lidar
The sounding wavelengths selected (299/341 nm) were used in the SLS ozone lidar designed.Its block diagram is shown in Figure 2. The fourth harmonics (266 nm) of the fundamental frequency of a Nd:YAG laser (LS-2134UT laser, LOTIS TII company, Minsk) is used as a laser radiation source, which is then SRS converted in hydrogen in the first (299 nm) and second (341 nm) Stokes components.The receiving telescope has been designed according to the Newton scheme on the basis of a primary mirror, 0.5 m in diameter, with a focal length of 1.5 m.The recording channel of the lidar is equipped with PMTs (R7207-01) and HAMAMATSU amplifier/discriminators (C3866).Lidar signals are recorded in the photocurrent pulse counting mode.To support PMT linear modes, a mechanical shutter is used, which cuts off a high-power optical signal from the nearest sounding zone.An automated unit for the output rotary mirror adjustment has been designed on the basis of computer-driven step motors.
The SRS cell is made from a tube (stainless steel) 3 cm in inner diameter and 1 m in length.Input and output windows are made of quartz Quartz Ultraviolet (KU).The pumping pulse energy at a wavelength of 266 nm is 60 mJ.The pumping power density required for SRS conversion is provided by lens L1 with a focal length of 1 m.It is mounted before the SRS cell and focuses radiation at its center.Confocal collimating lens L2 is mounted behind the cell.
Basic specifications of laser sources and receiving optical elements of the SLS ozone lidar: The fourth harmonics (266 nm) of the fundamental frequency of a Nd:YAG laser (LS-2134UT laser, LOTIS TII company, Minsk) is used as a laser radiation source, which is then SRS converted in hydrogen in the first (299 nm) and second (341 nm) Stokes components.The receiving telescope has been designed according to the Newton scheme on the basis of a primary mirror, 0.5 m in diameter, with a focal length of 1.5 m.The recording channel of the lidar is equipped with PMTs (R7207-01) and HAMAMATSU amplifier/discriminators (C3866).Lidar signals are recorded in the photocurrent pulse counting mode.To support PMT linear modes, a mechanical shutter is used, which cuts off a high-power optical signal from the nearest sounding zone.An automated unit for the output rotary mirror adjustment has been designed on the basis of computer-driven step motors.
The SRS cell is made from a tube (stainless steel) 3 cm in inner diameter and 1 m in length.Input and output windows are made of quartz Quartz Ultraviolet (KU).The pumping pulse energy at a wavelength of 266 nm is 60 mJ.The pumping power density required for SRS conversion is provided by lens L1 with a focal length of 1 m.It is mounted before the SRS cell and focuses radiation at its center.Confocal collimating lens L2 is mounted behind the cell.
Basic specifications of laser sources and receiving optical elements of the SLS ozone lidar: The efficiency of SRS conversion was measured versus the hydrogen pressure in the Raman cell, which was varied from 1 to 9 atm. Figure 3 shows the relative intensities of pumping radiation (266 nm), the first (299 nm) and the second (341 nm) Stokes components of SRS conversion as functions of hydrogen pressure at the SRS cell exit.Mirror diameter, m 0.5 8 Focal length, m 1.5 The efficiency of SRS conversion was measured versus the hydrogen pressure in the Raman cell, which was varied from 1 to 9 atm. Figure 3 shows the relative intensities of pumping radiation (266 nm), the first (299 nm) and the second (341 nm) Stokes components of SRS conversion as functions of hydrogen pressure at the SRS cell exit.The intensities of 299 and 341 nm lines become equal at a hydrogen pressure of 2 atm, which allows ozone sounding under equal radiation energies at these wavelengths.However, to increase the upper sounding limit, a pressure of 1 atm is more efficient, since the energy is redistributed toward the 299 nm line, which is absorbed by ozone stronger than the 341 nm line.

IASI/MetOp
The IASI is mounted onboard the European Space Agency meteorological satellite MetOp.The satellite monitors CO2, CH4, N2O, CO, O3, and HNO3 atmospheric gases; measures the temperature and humidity profiles in the troposphere and lower stratosphere in the near real time within the European Program "European Polar System".IASI provides for spectra of high radiometric quality with a resolution of 0.5 cm −1 in the range from 625 to 2760 cm −1 [20].The ozone profiles were retrieved from satellite sounding data in the range 1025-1075 cm −1 .
The satellite data were received by the 2.4 XLB satellite data reception station (Orbital Systems, USA) put into operation at IAO SB RAS in 2011 [21].The information from the station allows comparison between the satellite data and SLS lidar sounding data.

Results
The VDO profiles retrieved from lidar sounding data are compared with the IASI/MetOp measured profiles.
VDO was measured at SLS IAO SB RAS during 2015.Using the above described technique, the retrieved profiles of ozone were calculated for the upper troposphere-lower stratosphere and compared with the IASI/MetOpt measured profiles.The range of lidar signal detection was from 6 to 15 km.The ozone vertical profile error was from 6 to 11%.The retrieval error is estimated by the sum The intensities of 299 and 341 nm lines become equal at a hydrogen pressure of 2 atm, which allows ozone sounding under equal radiation energies at these wavelengths.However, to increase the upper sounding limit, a pressure of 1 atm is more efficient, since the energy is redistributed toward the 299 nm line, which is absorbed by ozone stronger than the 341 nm line.

IASI/MetOp
The IASI is mounted onboard the European Space Agency meteorological satellite MetOp.The satellite monitors CO 2 , CH 4 , N 2 O, CO, O 3 , and HNO 3 atmospheric gases; measures the temperature and humidity profiles in the troposphere and lower stratosphere in the near real time within the European Program "European Polar System".IASI provides for spectra of high radiometric quality with a resolution of 0.5 cm −1 in the range from 625 to 2760 cm −1 [20].The ozone profiles were retrieved from satellite sounding data in the range 1025-1075 cm −1 .
The satellite data were received by the 2.4 XLB satellite data reception station (Orbital Systems, USA) put into operation at IAO SB RAS in 2011 [21].The information from the station allows comparison between the satellite data and SLS lidar sounding data.

Results
The VDO profiles retrieved from lidar sounding data are compared with the IASI/MetOp measured profiles.
VDO was measured at SLS IAO SB RAS during 2015.Using the above described technique, the retrieved profiles of ozone were calculated for the upper troposphere-lower stratosphere and compared with the IASI/MetOpt measured profiles.The range of lidar signal detection was from 6 to 15 km.The ozone vertical profile error was from 6 to 11%.The retrieval error is estimated by the sum where e 1 is the absorption cross-section error, e 2 is the standard error of measurements in the photon counting mode, and e 3 is the scattering ratio error.The absorption cross-section error e 2 does not exceed 2% [18].The standard error of measurements in the photon counting mode e 2 is defined as The scattering ratio error e 3 is defined as where H calib -calibration height, N noise -noise signal, K-constant value that appears due to the assumptions of the processing technique and the estimated value 3 • (0.01) 2 .
The aerosol impact on the ozone profile in the troposphere and lower stratosphere is strong; therefore, the aerosol correction is to be used in the VDO retrieval algorithm.
Before the VDO retrieval, the scattering ratio is calculated by the lidar signal at the wavelength 341 nm.This allows the aerosol correction with the use of a real scattering ratio instead of its model values, and thus minimizes the aerosol impact on the ozone profile in the dynamic gaseous and aerosol medium at tropospheric and lower stratospheric altitudes.Figure 4 shows the mean error of VDO retrieval over all the measurement days.It should be noted that the lidar operates in the photon counting mode at a spatial resolution of 100 m, while the IASI profile has a step of 150 m and larger as the altitude increases, which is connected with the ozone profile retrieval algorithm of this interferometer [22].IASI O 3 retrieval errors, according to the estimates given in paper [23], reach values of 30%.The IASI curve is smoother and does not describe VDO in detail; the lidar measurements show the ozone behavior in more detail.
The coordinates and time of the lidar and satellite sounding are given in Table 4.We have analyzed the ozone profiles retrieved for April 2015 corrected to the temperature and aerosol and without corrections.One can see from Figure 5a that the efficiency of the corrections is significant.For the retrieval without temperature correction, the absorption cross-section constant σ299= 4.4 × 10 −19 cm 2 was used.High saturation of the aerosol component at tropospheric and lower stratospheric altitudes is corrected to both aerosol (terms C and D or Equations ( 12) and ( 15)) and temperature (term A or Equation ( 8)).Each correction used approaches the profiles to the IASI data.We have analyzed the ozone profiles retrieved for April 2015 corrected to the temperature and aerosol and without corrections.One can see from Figure 5a that the efficiency of the corrections is significant.For the retrieval without temperature correction, the absorption cross-section constant σ 299 = 4.4 × 10 −19 cm 2 was used.High saturation of the aerosol component at tropospheric and lower stratospheric altitudes is corrected to both aerosol (terms C and D or Equations ( 12) and ( 15)) and temperature (term A or Equation ( 8)).Each correction used approaches the profiles to the IASI data.Thus, the combination of corrections provides for a reliable ozone profile close to the IASI data (Figure 5b).The differences between profiles 3 and 4 are caused by the absorption cross-section values: the absorption cross-section in profile 4 is smaller than in profile 3; therefore, the values of profile 4 are higher than the values of profile 3 (Figure 5a).
Figure 6 shows all the considered cases of comparison of ozone profiles measured at the SLS and retrieved from IASI data (see Table 4).The measurements are reduced to the total altitude range 6-15 km for convenience.Increased ozone concentrations measured at the SLS as compared to IASI are seen in the range from 8.5 to 12.5 km throughout the whole observation period.
Thus, the combination of corrections provides for a reliable ozone profile close to the IASI data (Figure 5b).The differences between profiles 3 and 4 are caused by the absorption cross-section values: the absorption cross-section in profile 4 is smaller than in profile 3; therefore, the values of profile 4 are higher than the values of profile 3 (Figure 5a).
Figure 6 shows all the considered cases of comparison of ozone profiles measured at the SLS and retrieved from IASI data (see Table 4).The measurements are reduced to the total altitude range 6-15 km for convenience.Increased ozone concentrations measured at the SLS as compared to IASI are seen in the range from 8.5 to 12.5 km throughout the whole observation period.

Discussion
The following conclusions can be drawn from the analysis of the ozone profiles average over the period under study (Figures 7 and 8).
Figure 7a shows the mean lidar and satellite ozone profiles, and Figure 7b, the total difference between them (lidar-IASI) over all measurement days, as well as the standard deviation of this difference, minimum and maximum, and the mean.The error is calculated with the use of the standard deviation of the difference (lidar-IASI) and the ratio (lidar-IASI)/lidar over all days of measurements.Figure 7c shows the total difference between the lidar and IASI profiles normalized to the lidar profile for each measurement day, and the standard deviation with the minimum, maximum, and mean.
Figure 8a shows the mean lidar and Krueger model profiles, and Figure 8b, the total difference between them (lidar-Krueger model) over all measurement days, as well as the standard deviation of this difference, minimum and maximum, and the mean.Figure 8c shows the total difference between the lidar and Krueger model profiles normalized to the lidar profile for each measurement day, and the standard deviation with the minimum, maximum, and mean.
Thus, in Figure 7b, one can trace the variability of the difference between the lidar and satellite VDO in absolute units.The difference minima show how IASI data exceed the lidar data on VDO, and the minima, vice versa.The mean difference shows the difference between the lidar and IASI data over all measurement days.The normalized difference in Figure 7c shows more clearly the deviations between the values over all measurement days in percentage.
The difference between the lidar and IASI profiles of the ozone concentration grows with altitude.Hence, the mean difference in the ozone concentrations (lidar-IASI) varies from −1.56 × 10 12 mol.cm −3 at an altitude of 15 km to 0.53 × 10 12 mol.cm −3 at 13 km (Figure 6b).The maximal differences over all the profiles are from −0.01 × 10 12 mol.cm −3 at an altitude of 7.2 km to 1.08 × 10 12 mol.cm −3 at 14.3 km.The minimal differences over all the profiles vary from −0.77 × 10 12 mol.cm −3 at 15 km to 0.39 × 10 12 mol.cm −3 at 12.3 km.

Discussion
The following conclusions can be drawn from the analysis of the ozone profiles average over the period under study (Figures 7 and 8).
Figure 7a shows the mean lidar and satellite ozone profiles, and Figure 7b, the total difference between them (lidar-IASI) over all measurement days, as well as the standard deviation of this difference, minimum and maximum, and the mean.The error is calculated with the use of the standard deviation of the difference (lidar-IASI) and the ratio (lidar-IASI)/lidar over all days of measurements.Figure 7c shows the total difference between the lidar and IASI profiles normalized to the lidar profile for each measurement day, and the standard deviation with the minimum, maximum, and mean.
Figure 8a shows the mean lidar and Krueger model profiles, and Figure 8b, the total difference between them (lidar-Krueger model) over all measurement days, as well as the standard deviation of this difference, minimum and maximum, and the mean.Figure 8c shows the total difference between the lidar and Krueger model profiles normalized to the lidar profile for each measurement day, and the standard deviation with the minimum, maximum, and mean.
Thus, in Figure 7b, one can trace the variability of the difference between the lidar and satellite VDO in absolute units.The difference minima show how IASI data exceed the lidar data on VDO, and the minima, vice versa.The mean difference shows the difference between the lidar and IASI data over all measurement days.The normalized difference in Figure 7c shows more clearly the deviations between the values over all measurement days in percentage.
The difference between the lidar and IASI profiles of the ozone concentration grows with altitude.Hence, the mean difference in the ozone concentrations (lidar-IASI) varies from −1.56 × 10 12 mol.cm −3 at an altitude of 15 km to 0.53 × 10 12 mol.cm −3 at 13 km (Figure 6b).The maximal differences over all the profiles are from −0.01 × 10 12 mol.cm −3 at an altitude of 7.2 km to 1.08 × 10 12 mol.cm −3 at 14.3 km.The minimal differences over all the profiles vary from −0.77 × 10 12 mol.cm −3 at 15 km to 0.39 × 10 12 mol.cm −3 at 12.3 km.
Ozone shows pronounced annual variations; therefore, to find relative errors of its measurements with IASI, the difference in the concentrations was normalized to the lidar data: (Lidar-IASI)/Lidar.These data are shown in Figure 7c; it is seen that the mean relative difference is positive in the altitude ranges 6-6.5 and 8.4-14.6 km and attains 23.2% at an altitude of 12.4 km.The mean relative difference is negative in the altitude range from 6.5 to 8.4 km and 14.6-15 km: it attains the negative maximum of −33.8% at 7.5 km.The maximal relative difference changes from 3.45 to 60% in the range 6-15 km.The minimal relative difference over all the profiles in this range is −120.6% at 7.7 km; it attains 15.4% at 12.3 km.The relative difference varies in the range from −19.3 to 60% at 6 km and from −70.5 to 28.8% at 15 km.
Thus, the intercomparison performed shows that the absolute differences in the lidar and IASI measured ozone concentrations can change from −0.77 to 1.08 × 10 12 mol.cm −3 ; therefore, the relative difference is in the range from −120.6 to 60%.
It should be noted that the retrieved profiles of altitude distribution of the ozone concentration tend to IASI satellite data profiles more than to the Krueger model.Ozone shows pronounced annual variations; therefore, to find relative errors of its measurements with IASI, the difference in the concentrations was normalized to the lidar data: (Lidar-IASI)/Lidar.These data are shown in Figure 7c; it is seen that the mean relative difference is positive in the altitude ranges 6-6.5 and 8.4-14.6 km and attains 23.2% at an altitude of 12.4 km.The mean relative difference is negative in the altitude range from 6.5 to 8.4 km and 14.6-15 km: it attains the negative maximum of −33.8% at 7.5 km.The maximal relative difference changes from 3.45 to 60% in the range 6-15 km.The minimal relative difference over all the profiles in this range is −120.6% at 7.7 km; it attains 15.4% at 12.3 km.The relative difference varies in the range from −19.3 to 60% at 6 km and from −70.5 to 28.8% at 15 km.
Thus, the intercomparison performed shows that the absolute differences in the lidar and IASI measured ozone concentrations can change from −0.77 to 1.08 × 10 12 mol.cm −3 ; therefore, the relative difference is in the range from −120.6 to 60%.
It should be noted that the retrieved profiles of altitude distribution of the ozone concentration tend to IASI satellite data profiles more than to the Krueger model.

Conclusions
The results of using the developed techniques for VDO retrieval with temperature and aerosol correction in the altitude range 6-15 km during the DIAL lidar sounding of the atmosphere confirm the prospects of the wavelengths chosen (299 and 341 nm) for ozone lidar sounding.
Results of lidar measurements at 299 and 341 nm agree with model estimates, which point towards acceptable accuracy of ozone sounding in altitudes near 6-15 km.
At present, works are being carried out on optimization of optical and photoelectronic elements of the lidar signal detection system with the aim of increasing the upper limit of the sounding and improving the measurement accuracy.A more effective and comparatively easy-to-use ВаВ2O4 (BBO) crystal has been mounted as a fourth harmonic converter.

Conclusions
The results of using the developed techniques for VDO retrieval with temperature and aerosol correction in the altitude range 6-15 km during the DIAL lidar sounding of the atmosphere confirm the prospects of the wavelengths chosen (299 and 341 nm) for ozone lidar sounding.
Results of lidar measurements at 299 and 341 nm agree with model estimates, which point towards acceptable accuracy of ozone sounding in altitudes near 6-15 km.
At present, works are being carried out on optimization of optical and photoelectronic elements of the lidar signal detection system with the aim of increasing the upper limit of the sounding and improving the measurement accuracy.A more effective and comparatively easy-to-use BaB 2 O 4 (BBO) crystal has been mounted as a fourth harmonic converter.

Figure 1 .
Figure 1.Structure of the software for vertical distribution of ozone (VDO) retrieval from laser sounding data.

Figure 3 .
Figure 3. Relative intensities of pumping (266 nm), the first (299 nm) and the second (341 nm) Stokes conversion components of SRS conversion as functions of hydrogen pressure.

Figure 3 .
Figure 3. Relative intensities of pumping (266 nm), the first (299 nm) and the second (341 nm) Stokes conversion components of SRS conversion as functions of hydrogen pressure.

Figure 4 .
Figure 4. Mean error of ozone profile retrieval over 2015.

Figure 5
Figure 5 exemplifies the comparison of lidar and satellite vertical profiles of ozone measured on 26 April 2015; the VDO profiles are compared with the Krueger model [24].15 H, km 15 H, km

Figure 4 .
Figure 4. Mean error of ozone profile retrieval over 2015.

Figure 4 .
Figure 4. Mean error of ozone profile retrieval over 2015.

Figure 5
Figure 5 exemplifies the comparison of lidar and satellite vertical profiles of ozone measured on 26 April 2015; the VDO profiles are compared with the Krueger model [24].

Figure 5 .
Figure 5. VDO profiles retrieved in comparison with the Kruger model: ozone profile without temperature and aerosol corrections (1), ozone profile with temperature and without aerosol corrections (2), ozone profile without temperature and with aerosol corrections (3), combination of corrections (aerosol and temperature) (4) (a); VDO profiles retrieved in comparison with the Kruger model and IASI satellite data (b).

Figure 5 .
Figure 5. VDO profiles retrieved in comparison with the Kruger model: ozone profile without temperature and aerosol corrections (1), ozone profile with temperature and without aerosol corrections (2), ozone profile without temperature and with aerosol corrections (3), combination of corrections (aerosol and temperature) (4) (a); VDO profiles retrieved in comparison with the Kruger model and IASI satellite data (b).

Figure 6 .
Figure 6.Intercomparison of vertical profiles of ozone for 2015.

Figure 6 .
Figure 6.Intercomparison of vertical profiles of ozone for 2015.

Acknowledgments:
The work was supported by the Russian Science Foundation (Agreement No. 15-17-10001 in carrying out lidar measurements of ozone) and President of the Russian Federation within the Program for Support of Leading Scientific Schools (grant No. NSh-8199.2016.5).uthors.Alexey V. Nevzorov and Oleg Romanovskii defined the research theme.All authors performed the

Table 1 .
Specifications of lidars used for measurements of stratospheric ozone.

Table 2 .
Sets of wavelengths that correspond to Stokes (C) frequencies of SRS conversion in H 2 , D 2 , and CO 2 .

Table 4 .
Coordinates and time of satellite and lidar sounding used for comparison of VDO profiles for 2015.Figure 5 exemplifies the comparison of lidar and satellite vertical profiles of ozone measured on 26 April 2015; the VDO profiles are compared with the Krueger model [24].