Full Text

Turn on search term navigation

Atmos. Meas. Tech., 9, 11811205, 2016 www.atmos-meas-tech.net/9/1181/2016/ doi:10.5194/amt-9-1181-2016 Author(s) 2016. CC Attribution 3.0 License.

Anatoli Chaikovsky1, Oleg Dubovik2, Brent Holben3, Andrey Bril1, Philippe Goloub2, Didier Tanr2,Gelsomina Pappalardo4, Ulla Wandinger5, Ludmila Chaikovskaya1, Sergey Denisov1, Jan Grudo1, Anton Lopatin1, Yana Karol1, Tatsiana Lapyonok2, Vassilis Amiridis6, Albert Ansmann5, Arnoud Apituley7,Lucas Allados-Arboledas8, Ioannis Binietoglou9, Antonella Boselli10,4, Giuseppe DAmico4, Volker Freudenthaler11, David Giles3, Mara Jos Granados-Muoz8, Panayotis Kokkalis6,12, Doina Nicolae9, Sergey Oshchepkov1,Alex Papayannis12, Maria Rita Perrone13, Alexander Pietruczuk14, Francesc Rocadenbosch15, Michal Sicard15, Ilya Slutsker3, Camelia Talianu9, Ferdinando De Tomasi13, Alexandra Tsekeri6, Janet Wagner5, and Xuan Wang10

1Institute of Physics, NAS of Belarus, Minsk, 220072, Belarus

2LOA, Universite de Lille, Lille, 59650, France

3NASA Goddard Spaceight Center, Greenbelt, MA 20771, USA

4Consiglio Nazionale delle Ricerche Istituto di Metodologie per lAnalisi Ambientale (CNR-IMAA), Potenza, 85050, Italy

5Leibniz Institute for Tropospheric Research, Leipzig, 04318, Germany

6Institute for Astronomy, Astrophysics, Space Applications and Remote Sensing,

National Observatory of Athens, Athens, 15236, Greece

7KNMI Royal Netherlands Meteorological Institute, De Bilt, 3731, the Netherlands

8Andalusian Institute for Earth System Research (IISTA-CEAMA), University of Granada, Autonomous Government of Andalusia, Granada, 18071, Spain

9National Institute of R&D for Optoelectronics, Magurele, 77125, Romania

10Consorzio Nazionale Interuniversitario per le Scienze Fisiche della Materia, Naples, 80138, Italy

11Ludwig-Maximilians Universitt, Meteorological Institute, Mnchen, 80539, Germany

12National Technical University of Athens, Department of Physics, Athens, 15780, Greece

13Consorzio Nazionale Interuniversitario per le Scienze Fisiche della Materia (CNISM) and Universita del Salento, Lecce, 73100, Italy

14Institute of Geophysics, Polish Academy of Sciences, Warsaw, 01-452, Poland

15Remote Sensing Laboratory (RSLAB), Department of Signal Theory and Communications, Universitat Politcnica de Catalunya (UPC)/Institute for Space Studies of Catalonia (IEEC), Barcelona, 08034, Spain

Correspondence to: Anatoli Chaikovsky ([email protected])

Received: 19 October 2015 Published in Atmos. Meas. Tech. Discuss.: 7 December 2015 Revised: 17 February 2016 Accepted: 29 February 2016 Published: 21 March 2016

Abstract. This paper presents a detailed description of

LIRIC (LIdar-Radiometer Inversion Code) algorithm for simultaneous processing of coincident lidar and radiometric (sun photometric) observations for the retrieval of the aerosol concentration vertical proles. As the lidar/radiometric input data we use measurements from European Aerosol Research Lidar Network (EARLINET) lidars and collocated

sun-photometers of Aerosol Robotic Network (AERONET).

The LIRIC data processing provides sequential inversion of the combined lidar and radiometric data. The algorithm starts with the estimations of column-integrated aerosol parameters from radiometric measurements followed by the retrieval of height dependent concentrations of ne and coarse aerosols from lidar signals using integrated column characteristics of

Published by Copernicus Publications on behalf of the European Geosciences Union.

Lidar-Radiometer Inversion Code (LIRIC) for the retrieval of vertical aerosol properties from combined lidar/radiometer data: development and distribution in EARLINET

1182 A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC)

aerosol layer as a priori constraints. The use of polarized lidar observations allows us to discriminate between spherical and non-spherical particles of the coarse aerosol mode.

The LIRIC software package was implemented and tested at a number of EARLINET stations. Intercomparison of the LIRIC-based aerosol retrievals was performed for the observations by seven EARLINET lidars in Leipzig, Germany on25 May 2009. We found close agreement between the aerosol parameters derived from different lidars that supports high robustness of the LIRIC algorithm. The sensitivity of the retrieval results to the possible reduction of the available observation data is also discussed.

1 Introduction

The aerosol impact on the radiation balance of the atmosphere is an important climate forcing factor. In addition, aerosol particles are among the unhealthiest air pollutants.This is made more severe by rapid propagation of pollutants in the atmosphere that expands local ecocatastrophes to a global scale. Therefore, the monitoring of the aerosol evolution and transport in the atmosphere is an obligatory prerequisite for predicting climatic and ecological changes.

Sun-radiometer and lidar networks contribute to aerosol remote sensing. The global Aerosol Robotic Network (AERONET) of ground-based sunsky-scanning radiometers (e.g. Holben et al., 1998) provides reliable data on columnar aerosol properties from more than 200 globally distributed sites. The results of AERONET observations are the aerosol optical thickness (AOT) obtained from direct sun observations and additional microphysical and optical properties of aerosol particles (single scattering albedo, volume distribution of aerosol particles, complex refractive index, fraction of spherical particles, etc.) derived by the inversion of direct and scattered radiation measurements (Dubovik and King, 2000; Dubovik et al., 2002, 2004). The regional radiometer network SKYNET was established in the southeastern Asian regions (Takamura et al., 2004), and it employs its own equipment and processing procedure (Hashimoto et al., 2012).

The lidar measurements are used to provide information on the vertical variability of the aerosol characteristics. Currently, lidar networks, such as the European Aerosol Research Lidar Network (EARLINET) (Bsenberg et al., 2000;Pappalardo et al., 2014), the micro-pulse lidars network (MPL-Net) (Welton et al., 2002), the Asian dust network (AD-Net) (Murayama et al., 2001), the lidar network in former Soviet Union countries CIS-LiNet (Chaikovsky et al., 2005), the northeast American CREST Lidar Network (CLN) (Hoff et al., 2009), and the Latin America Lidar Network LALINET (Antua et al., 2012), monitor aerosol vertical distributions in the atmosphere over vast regions of the Earth.The Global Atmosphere Watch (GAW) Aerosol Lidar Ob-

servation Network (GALION), also known as the network of networks (e.g. Bsenberg and Hoff, 2007), was established under the aegis of GAW to coordinate lidar activity all over the world. The outcome of the lidar observations are presented in the lidar network databases as vertical proles of aerosol backscatter and extinction coefcients.

Aerosol columnar properties from AERONET and aforementioned vertical proles of aerosol parameters from lidar networks are complementary pieces of information characterizing aerosol properties. Nowadays, lidars and sunsky-scanning radiometers are among the basic tools in comprehensive experiments aimed at studying the transformation and transport of smoke (e.g. Lund Myhre et al., 2007; McK-endry et al., 2011; Colarco et al., 2004), dust (e.g. Ansmann et al., 2009; McKendry et al., 2007; Mller et al., 2003; Papayannis et al., 2008), and volcanic ash (e.g. Ansmann et al., 2010, 2011, 2012; Papayannis et al., 2012; Gasteiger et al., 2011). A number of SKYNET sites (Takamura et al., 2004) and most of the EARLINET stations are equipped with lidar and radiometer instruments. Further enhancement of the aerosol characterization is expected from the synergy of co-located radiometer and lidar observations. Namely, the coordination of measurement procedures of the two systems and the derivation of aerosol parameters from combined measurements results in advanced characterization of the aerosol layer with a superior performance compared to the aerosol information that would have been obtained from independent processing of lidar and radiometer data.

The idea of combined lidar and radiometer sounding (LRS) for retrieving vertical distributions of aerosol characteristics was rst proposed by Chaikovsky et al. (2002), and it gave rise to the development of the lidarradiometer synergetic algorithms (e.g. Chaikovsky et al., 2004a, b).Later, in 2012 under the ACTRIS Research Infrastructure project within the European Union Seventh Framework Programme, the algorithm and software package, named LIRIC (LIdar-Radiometer Inversion Code), was developed for processing data of EARLINET measurements. LIRIC is based on processing co-located lidar and radiometer measurements by using a two-step sequential inversion. First, the radiometer data was processed according to the standard AERONET inversion algorithm. Then, rst-step results are used as a priori constrains on aerosol properties for lidar data processing.

First application of LIRIC technique to the actual data processing was presented by Chaikovsky et al. (2004a). In that study, the technique was adapted to the EARLINETAERONET stations in Minsk (Belarus) and Belsk (Poland) (e.g. Chaikovsky et al., 2004c, 2010a; Pietruczuk and Chaikovsky, 2007). Results of the LRS observations were of interest for the study of long-range aerosol transport in the eastern European region (Kabashnikov et al., 2010;Chaikovsky et al., 2010b; Papayannis et al., 2014).

Another algorithm for data processing in combined lidarand-radiometer experiments exploits the decomposition of the AERONET column-integrated aerosol size distribution

Atmos. Meas. Tech., 9, 11811205, 2016 www.atmos-meas-tech.net/9/1181/2016/

A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC) 1183

into log-normal modes and selection of some of these modes for the characterization of aerosol layers using measured lidar data (Cuesta et al., 2008).

The LRS technique for retrieving the aerosol concentration proles from single-wavelength lidar measurements at the MPLNET (Micro-Pulse Lidar Network) stations collocated with the sunsky radiometer sites of AERONET was developed by Ganguly et al., 2009a. Then this method was applied to processing of the combined AERONET and space CALIOP lidar data (Ganguly et al., 2009b).

Besides, the single-wavelength POLIPHON technique was developed as an alternative (e.g. Tesche et al., 2009;Ansmann et al., 2012). This technique retrieves particle volume concentration proles of aerosol separately for ne and coarse fractions. The algorithm relies on the measured proles of the particle linear depolarization ratio and lidar ratio, and it does not require the assumption of a specic particle shape. Columnar concentrations of aerosol modes retrieved by AERONET are used in POLIPHON as additional input data. The algorithm POLIPHON is designed for the data processing in lidar sounding of the aerosol layers with coarse non-spherical particles (dust, volcano ash).

In recent years, the LRS technique has been implemented within the advanced research network ACTRIS in the frame of EU 7th Framework Programme project. To date, a number of joint EARLINET/AERONET stations have implemented regular atmospheric observations using LIRIC for processing combined sun-radiometer and lidar-measured data (e.g.Chaikovsky et al., 2012; Papayannis et al., 2014; Tsekeri et al., 2013). The aerosol model and mathematical basis of the LIRIC algorithm became the prerequisite for further development of algorithms for simultaneous inversion of combined lidarradiometer measurements, e.g. GARRLIC (Generalized Aerosol Retrieval from Radiometer and Lidar Combined data) (Lopatin et al., 2013) and the results of ground/satellite closure experiments (Dubovik et al., 2014). Note that LIRIC technique should not be regarded only as a basis for new algorithms (e.g. POLIPHON or GARRLIC). LIRIC might be superior to them for many aerosol scenarios: it allows one, for example, to distinguish between ne and coarse spherical fractions (unlike POLIPHON) or distinguish between spherical and non-spherical coarse particles (unlike GARRLIC).At the same time, a comprehensive description of the LIRIC algorithm has not been yet documented in detail.

This paper describes the basic physical and mathematical aspects of LIRIC algorithm with all necessary equations, thus lling up this gap. The appendices contain the details of the inversion scheme and can be useful for advanced users to modify and improve this code.

2 The algorithm concept and structure

The aerosol retrievals from combined lidar and radiometer measurements belong to a class of ill-posed inverse problems that, in particular, is characterized by non-unique and

highly unstable solutions arising even under small measurement or simulation errors. In practice, the solution of the ill-posed problems requires the introduction a priori information (e.g. Turchin et al., 1971; Tikhonov and Arsenin, 1977;Twomey, 1977; Tarantola, 1987; Rodgers, 2000). LIRIC algorithm was designed on the basis of multi-term LSM (least square method) (Dubovik, 2004). This method was implemented in AERONET data processing (Dubovik and King, 2000) and then it was rened in the retrieval algorithms for the data processing of the combined optical measurements (e.g. Dubovik, 2004; Dubovik et al., 2011, 2014; Lopatin et al., 2013).

The inversion algorithm LIRIC can be divided into three key procedures (e.g. Tarantola, 1987): (i) parameterization of the object under study (i.e. development of the aerosol layer model); (ii) forward modelling, i.e. derivation of the equations that relate observed signals with specied parameters of the aerosol model; and (iii) inverse modelling or retrieval of the target parameters of the aerosol model that minimize discrepancies between the measured and the calculated input signals.

2.1 Combined lidar/radiometer experiment and aerosol model

The lidar/radiometer input data assumed to come from measurements of EARLINET lidars (e.g. Matthias et al., 2004; Freudenthaler et al., 2010) and spectral-scanning sunradiometers of AERONET (Holben et al., 1998). The majority of EARLINET lidar stations provides daytime measurements of elastic backscatter signals at three wavelengths (355, 532, and 1064 nm) and two cross/parallel-polarization components of the signal at a single wavelength. Additional information on aerosol parameters is expected to come from daytime Raman lidar measurements. (The current version of LIRIC algorithm is not designed for using Raman lidar data.)

Radiometric data includes results of direct-sun and almucantar (scanning) measurements (Holben et al., 1998;Dubovik and King, 2000). Direct-sun measurements are carried out in 15 min intervals. Almost-clear-sky measurements are required to obtain almucantar data, and approximately two to six successful measurements are made during the daytime under favourable meteorological conditions at EARLINET/AERONET stations. Under these circumstances, time synchronization of lidar and radiometric observations usually means nearly simultaneous measurements within the same 1 h interval.

These radiometric measurements enable the retrieval of the aerosol properties over the entire atmospheric column.Thus, except for volcanic events, the maximum lidar sounding height, hmax, can be limited to the tropopause level because the stratospheric aerosol layer does not signicantly contribute to columnar aerosol optical parameters. In contrast, aerosols in the lower troposphere are key contributors to the observed columnar characteristics. Consequently, it is

www.atmos-meas-tech.net/9/1181/2016/ Atmos. Meas. Tech., 9, 11811205, 2016

1184 A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC)

desirable to perform the lidar sounding from the lowest possible altitude. Likewise, the contribution of the bottom layer (which is not observed by lidar) to the columnar optical parameters must be small enough to be modelled by a homogeneous layer with the same aerosol parameters as at the lowest level of lidar sounding. In practice, the lower sounding limit for most of the lidar measurements in EARLINET is about 200 m or more that can be too high especially for low boundary layers in winter seasons. It should be decreased in winter to compensate reduction of the boundary layer height. Therefore, lidar measurements in the lower layer have to be carried out by a second, complementary receiving system with smaller objective and larger eld of view or by sounding the atmosphere along a slant trajectory.

The choice of the optical aerosol model is a key step of the retrieval algorithm. The optical model should be constructed following the principle of parsimony or Occams razor: the number of aerosol parameters has to be minimal but complete in order to provide unbiased retrieval from available measurements.

In this work, we use the AERONET model approach to characterize the aerosol layer of the atmosphere (Dubovik and King, 2000): aerosols are modelled by several modes with a certain aerosol particle size distribution, wherein each mode is a mixture of homogeneous spherical particles and randomly oriented spheroids (Dubovik et al., 2002, 2006).The distribution of the spheroid aspect ratio is xed. The number of aerosol modes, K, depends on specication of the lidar data. If we use only total (scalar) backscatter lidar measurements, the aerosol model includes ne and coarse modes (K = 2). There is boundary size between ne and coarse

fractions in the algorithm, which is determined as the value in 0.1940.576 m range that corresponds to a minimum of the column particle volume size distribution, dV (r)/dlnr.If measurements of cross and parallel co-polarized components are available, spherical and non-spherical particles of the coarse mode are considered as two different fractions (K = 3).

Thus, two sets of parameters characterize the aerosol layer.

1. The rst is a number of columnar aerosol parameters retrieved from radiometer measurements (Dubovik and King, 2000; Dubovik et al., 2002, 2006). This set of parameters is formed by the following: (1) the total content of each aerosol mode, (i.e. columnar volume concentrations),

CVk =

rmax,k

(3) complex refractive indices at the wavelength , m( ) = n( ) + i ( ); (4) the sphericity, sph (the

ratio of spherical particles volume to the total volume); (5) aerosol optical thickness (AOT) of the kth aerosol mode,[notdef]k j

; (6) the single scattering albedo for the kth aerosol mode, [pi1]k( ); (7) the elements of the backscattering matrix, P kx,x( ,180 ); and (8) coefcients ak and bk, which determine optical extinction and backscatter characteristics of aerosol particles for the k-aerosol mode (see Sect. 3.1). Parameters (1)(4) are the independent state variables, whilst parameters (5)(8) are derived from the state variables. Parameters m( ) and sph are assumed the same for particles of all sizes.

Denitions and a detailed description of the columnar aerosol parameters are available at the AERONET information system; cloud screening and quality control algorithms were described by Holben et al. (2006).

2. The second set of parameters is represented by the height, h, distributions of particle volume concentrations (PVC) for each of aerosol mode, ck(h), which dene the vertical variability of the aerosol features.

A lack of lidar data to resolve height-variation of aerosol microstructure motivates the assumption of altitude-independent microphysical parameters of the aerosol modes.

2.2 Algorithms structure

Two options of the retrieval procedure for the processing LRS data have been developed:

1. First one deals with sequential inversion of lidar and radiometer data. It is carried out by preliminary calculation of the column parameters dened in Sect. 2.1 from radiometric measurements by using the AERONET inversion algorithm (Dubovik and King, 2000), followed by subsequent inversion of the height distribution parameters by using lidar data with columnar characteristics of aerosol layer passed as a priori data (Chaikovsky et al., 2012);

2. Second option suggests simultaneous inversion approach for retrieving optimal parameters of the aerosol model by using a joint inversion procedure from combined lidar and radiometer data.

While the sequential algorithm could be considered as an unsophisticated inversion procedure to combine lidar and AERONET data, the parallel inversion method leads, in principle, to more effective estimation of aerosol parameters because it allows simultaneously retrieved columnar aerosol parameters to be specied in accordance with the additional lidar data. Currently, the simultaneous inversion algorithm for a two-component aerosol model is implemented in GARRLIC (Lopatin et al., 2013). Similar aerosol mode concentration proles and residual discrepancies between measured

Atmos. Meas. Tech., 9, 11811205, 2016 www.atmos-meas-tech.net/9/1181/2016/

[integraldisplay]

dVk (r)dlnr dlnr, (1)

where rmin,k and rmax,k is the minimum and the maximum radius of the kth aerosol mode (k =

1,...,K) , respectively; (2) the particle volume size distributiondVk(r)/dlnr for each aerosol mode;

rmin,k

A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC) 1185

and calculated input signals are obtained from both retrieval procedures when processing experimental data (Lopatin et al., 2013).

Advantages of the simultaneous inversion approach are expected for more involved measurements, such as in the unied algorithm GRASP (Generalized Retrieval of Aerosol and Surface Properties), which aimed at characterizing atmospheric properties from remote ground and satellite observations (Dubovik et al., 2014).

LIRIC algorithm described below was created on the base of the sequential inversion approach. Figure 1 shows the structure of the algorithm.

The algorithm is divided into several rather independent modules to provide exibility of the software package. Module 1 (preprocessing of lidar data) creates a set of smoothed and normalized lidar signals, L , covariance matrix, [Omega1]L, and setting parameters (type of lidar measurement, sounding wavelength, geographical coordinates of lidar station and date of measurement, etc.) for modelling aerosol and molecular layers. Module 2 (recalculation of radiometer data) estimates columnar parameters of the aerosol model for lidar sounding wavelengths. Level 1.5 or Level 2.0 AERONET data are acceptable as input data in LIRIC. (These data are inputs to Module 2). Initial proles of the aerosol-mode concentrations, c0k(h), as well as molecular (Rayleigh) extinction, r( ,h), and molecular backscatter coefcients, r( ,h), are generated by Module 3 (atmospheric model).Module 4 (forward model) calculates arrays of lidar signals, Lj

cm1k(h)

S j j ,h

, (2)

S j ( ,href)

exp 2r( j ,h,href)

where S j j ,h

[parenrightbig]

= P j( j ,h)h2; j j ,href

[parenrightbig]

is the value of

S j j ,h

[parenrightbig]

, and columnar volume concentrations,

CV,m1k, given aerosol concentration proles, cm1k(h), in iterative inversion procedure, where m stands for the mth retrieval iteration and j is the number of the receiving channel. Inversion parameters, constraints on the smoothness characteristics, and error signals for the sensitivity test are passed to the algorithm by Module 5 (inversion settings & error modelling). The sensitivity test (see Sect. 6) was designed to estimate the response of the retrieval results to measurement errors and/or uncertainties of input data. Module 6 (numerical inversion) is responsible for tting aerosol-mode concentration proles for the retrieved aerosol model, cm1k(h), given measured data and a priori information.

3 Forward modelling of LRS experiment

Range-corrected normalized lidar signals and columnar-aerosol parameters retrieved from radiometer measurements are the input data to the LRS processing procedure (see Fig. 1). Below, we dene a set of basic equations that are needed for the forward modelling of the measured quantities as well as to estimate the error-covariance matrix.

3.1 Basic lidar equations

The multichannel lidar carries out J different lidar measurements (j 2 1,...J ) that yields a set of lidar signal

records, P j, j 2 1,...J . The term different means that dif

ferent kinds of lidar measurements are performed, such as total intensity as well as cross- and parallel-polarized signal components at different wavelengths. Here we consider that each different lidar measurement is provided by a specic j-th channel. Parameter J stands for the number of lidar channels irrespective of the actual implementation of the lidar system.

Range-corrected normalized lidar signals are calculated at the preprocessing stage of the inversion procedure (Module 1 in Fig. 1):

L j (h) =

at the reference point, href is usually dened in the end of the sensing range, r( j ,h,href) is the molecular optical thickness related to the range of (h,href), j is the wavelength, and h is the height. The set of lidar signals, L j (h), constitutes the input lidar vector, L .

The lidar system provides measurements from the lowest to the highest altitude levels specied by hmin and hmax, respectively. Currently, it is assumed that the radiometer is co-located at a height of h0 < hmin, so columnar aerosol optical properties of the layer h0 < h < hmin are to be taken into consideration. If there is no information on the aerosol parameters in the surface layer, this layer is assumed to be homogeneous. Under this assumption, scattering parameters for the altitude range h0 < h < hmin of the lidar vector L are set equal to the values at hmin .

The relationship between the measured lidar signals L ( ) and the aerosol mode concentration, ck(h), can be written as follows:

L = L( ,ck(h),ak,bk) + [Delta1]L, (3) where [Delta1]L is the vector of measurement uncertainties. Here, an asterisk (*) denotes measured and no-asterisk denotes model estimated.

Since function L(...) in Eq. (3) depends on the type of lidar measurement, it is expedient to introduce special parameter, pj 2 1,2,...,U, that indicates the type of measurement

associated to the j-channel of the lidar, and U is a number of the types. In our case, pj 2 1,2,3, indicates total intensity,

cross-polarized, and parallel-polarized measurements, correspondingly.

The lidar functions, Lj,pj (...), for the pj -type measurements are dened by the following equations:

www.atmos-meas-tech.net/9/1181/2016/ Atmos. Meas. Tech., 9, 11811205, 2016

1186 A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC)

Figure 1. Flowchart of LIRIC algorithm. Details are in Sect. 2.2.

by the equation for the total backscatter signal,

Lj,1 j ,h

[parenrightbig]

a,1( j ,h) + r( j ,h) Rj,1( j ,href) r( j ,href)

exp

[integraldisplay]

a,2( j ,h) + [notdef] a,3( j ,h) + ~+[notdef]~+1 r( j ,h)[parenrightBig] ~+[notdef]

~+1 r( j ,h)Reff( j ,href)

exp

a( j ,h)dh

, (4)

href

where

Rj,1( j ,h) =

[integraldisplay]

a( j ,h)dh

, (8)

href

a,1( j ,h) + r( j ,h) r( j ,h) ; (5)

by the equation for the parallel-polarized signal component,

Lj,3 j ,h

[parenrightbig]

a,3( j ,h) +

where

Reff( j ,h) =[parenleftBigg]

(~ + [notdef])

a,2( j ,h) +

~+1 ,r( j ,h)

~+1 r( j ,h)

[notdef]

(~ + [notdef])

a,3( j ,h) +

11+~ r( j ,h)

11+~ r( j ,h) [parenrightBigg]

. (9)

In Eqs. (4)(9), a,1, a,3, and a,2 denote the aerosol backscatter coefcient and its parallel- and cross-polarized components, respectively; a( j ,h) is the aerosol extinction coefcient; ~ j

[parenrightbig]

11+~ r( j ,h)

11+~ r( j ,href)Rj,3( j ,href)

exp

[integraldisplay]

a( j ,h)dh

, (6)

href

= r,2( j) r,3( j) is the ratio of cross- and parallel-polarized components of the molecular backscatter coefcient.

Different cross-talk factors contribute to the spurious signal in the cross-polarized receiving channel. These factors in-

Atmos. Meas. Tech., 9, 11811205, 2016 www.atmos-meas-tech.net/9/1181/2016/

where

Rj,3( j ,h) =

a,3( j ,h) +

11+~ r( j ,h)

; (7)

and by the equation for the cross-polarized signal component,

Lj,2 j ,h

[parenrightbig]

A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC) 1187

clude the residual of cross-polarized component of the laser beam, non-ideal adjustment of the polarization planes between transmitter/receiver channels and depolarization by optical elements. Equations (6) and (8) allow for these crosstalk effects in a similar manner to Chaikovskii (1990) and Biele et al. (2000). Thus, parameter [notdef] in Eqs. (8)(9) represents the leakage of the parallel component of the sounding beam into the cross-polarized lidar receiving channel. Parameter [notdef] is an instrument characteristic that is assumed to be a known quantity; i.e. it is not updated by the retrieval procedure.

The aerosol extinction and backscatter coefcients in the Eqs. (3)(9) are expressed as a function of the parameters of the following aerosol modes:

a( j ,h) =Xk

ck(h)ak( j ), (10)

a,1( j ,h) =Xk

ck(h)bk,1( j ), (11)

a,2( j ,h) =Xk

ck(h)bk,2( j ), (12)

and

a,3( j ,h) =Xk

ck(h)bk,3( j ). (13)

The coefcients ak( j ) and bk,x( j ), pointed out in Sect. 2.1, are determined by columnar optical parameters of aerosol modes:

ak( j ) =

Xi=1ck(hi)[Delta1]hi + [Delta1]V,k. (19)

The structure of the vectors V , c, and matrix H is consid

ered in Appendix C.

4 Numerical inversion

Statistical regularization technique (e.g. Turchin et al., 1971; Tarantola, 1987; Rodgers, 2000) considers errors, [Delta1]L and [Delta1]V , in Eqs. (3) and (18) as random variables. Under the additional assumption that errors have independent normal distributions, the multidimensional conditional probability density function (PDF) (or likelihood function) is dened by Chaikovsky et al. (2004a)

L , V [vextendsingle][vextendsingle][vextendsingle]

Ek( j )

CVk

[bracketleftbigg][parenleftbigg]

, (14)

[parenrightBig]

exp

[parenleftbigg][parenleftBigg]

L L(c)

T [Omega1]1L

bk,1( j ) =

14 [pi1]k( j )ak( j )P k1,1( j ,180 ), (15)

L L(c)

[parenrightbig]

[parenleftBig]

+ C V Hc[parenrightBig]T

bk,3( j ) =

14 [pi1]k( j )ak( j )

P k1,1( j ,180 ) + P k2,2( j ,180 )2 , (16)

[Omega1]1V

[parenleftBig]

C V Hc[parenrightBig][parenrightbigg][parenrightbigg][bracketrightbigg]

. (20)

[vextendsingle][vextendsingle][vextendsingle]

is the PDF of measurement vectors

L and V , L(c) is the vector function in Eq. (3), H is the

matrix in Eq. (18), c is the target retrieval vector of aerosol modes concentration, and [Omega1]L and [Omega1]V are the covariance matrices of error vectors [Delta1]L and [Delta1]V , respectively.

An extensively used tool for the regularization of an ill-posed problem is the application of a priori constraint on the smoothness of retrieved characteristics. LIRIC restricts the norms of the second differences of functions ck(hi). Following the statistical regularization approach (Turchin et al., 1971) we included a priori probability function,

Fapr (c) exp

1 2

Here, F

L , V

[parenrightBig]

14 [pi1]k( j )ak( j )

P k1,1( j ,180 ) P k2,2( j ,180 )2 , (17)

where[notdef]k is aerosol optical thickness for the kth aerosol

mode, [pi1]k( ) is the single scattering albedo for the kth aerosol mode, and P kx,x( ,180 ) are the elements of the backscattering matrix.

3.2 Forward model of radiometer data

In accordance with the multi-term LSM approach (Dubovik, 2004), the columnar concentrations of aerosol modes,Vk ,

obtained from radiometer measurements are formally considered in LIRIC as a result of additional independent measurements.

The equation for the vector, V , which is dened as the

measured columnar volume concentrations of the aerosol modes given vector of aerosol modes concentration, c(hi), i 2 1,...I, can be written in the following form:

C V = Hc + [Delta1]V , (18) where H is convolution matrix for summing the height-resolved concentration over the column; [Delta1]V is the vector of

C V uncertainties.The k-th component of the vector V is dened by the

following equation:

C Vk (ck (hi)) =

bk,2( j ) =

cT [Omega1]Sc

[parenrightBig][parenrightbigg]

(21)

into the retrieval procedure as the additional constraint. Here, [Omega1]S = ST2 Q12S2 is the smoothing matrix, S2 is the matrix of

www.atmos-meas-tech.net/9/1181/2016/ Atmos. Meas. Tech., 9, 11811205, 2016

1188 A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC)

the second-order differences, and Q2 is diagonal weighting matrix (Twomey, 1977; Dubovik et al., 2011).

The Bayes strategy (Turchin et al., 1971; Tarantola, 1987;

Rodgers, 2000) for solving an ill-posed problem combined with multi-term LSM technique (Dubovik, 2004; Dubovik et al., 2011) denes the solution in accordance with the

maximum a posteriori rule

c = argmin c

{[Psi1](c)[notdef],

where the objective or cost function, [Psi1](c), has the following multi-term representation (Dubovik, 2004; Dubovik et al., 2011)

[Psi1](c) = L L(c)

T [Omega1]1L

The coefcients V and S are so-called Lagrange multipliers that determine the weight of different contributors from each source of information (i.e. measurements and a-priori contribution) to the retrieval solution relative to the contribution of the rst data source (since L = 1). Equa

tions (22) and (24) are equivalent; however, Eq. (24) is more convenient for the analysis of the relative contribution from different data source.

If V , S ! 0, we return to a non-regularized solution

for vector c that is based solely on measured lidar data with the minimum discrepancy between measured and calculated input signals. This solution, however, could be non-physical, multivalued, and unstable. The possible solution space should be restricted by increasing the Lagrange multipliers despite the fact that it results in increasing discrepancy between measured and model signal. The algorithms to determine the Lagrange multipliers by nding a reasonable compromise between the solution quality and the closeness of the measured and model signals are described in Hansen (2001), Vogel (2002), and Doicu et al. (2010). The set of Lagrange multipliers is provided to LIRICs users along with software package. However, we do not consider this set as the ultimate one, and we allow it to be modied to meet users specications.

The nal step of the retrieval procedure is calculation of the concentration proles ck(hi) for each aerosol mode. Initial approximations c0k(hi) are set and stepwise improved to provide the minimum of the objective function (Eq. 25).

Increments are calculated by means of the Levenberg-Marquardt method (Levenberg, 1944; Marquardt, 1963).

The analytical expressions of the terms of Eq. (25), the covariance matrices, and the details of the inversion procedure are described in Appendices A, B, and C.

5 Program package for processing combined lidar and radiometer data

Figure 2 shows the structure of the software package that implements the LIRIC algorithm. A set of specic programs are joined in three sub-packages.

The sub-package LiOpt implements module (2) of the LIRIC algorithm (Fig. 1), which provides preprocessing of the AERONET retrieval products. Program AERLID recalculates the columnar optical characteristics for the lidar sounding wavelengths, including the elements of the scattering matrices for the spherical and non-spherical particles as well as for ne and coarse aerosol modes. Then, this code transfers data to the radiometer database.

The preprocessing of lidar data is carried out by the SignalSuite sub-package. It contains several programs. Among them are the following:

ULIS an operational program that provides measurement procedures and record of raw lidar data to Microsoft ACCESS database;

Atmos. Meas. Tech., 9, 11811205, 2016 www.atmos-meas-tech.net/9/1181/2016/

L L(c)

[parenleftBig]

+ C V Hc[parenrightBig]T [Omega1]1V

[parenleftBig]

C V Hc[parenrightBig]

. (22)

We assume that the errors [Delta1]L in Eq. (3) and [Delta1]V in Eq. (18) are uncorrelated. In this case, the non-zero diagonal elements of the covariance matrices [Omega1]L and [Omega1]V are the variances of the elements of the vectors [Delta1]L and[Delta1]V , respectively.

Since the minimization procedure does not prescribe a residual value for[Psi1](c), it is convenient to reformulate weight matrices as follows (Dubovik, 2004):

[Omega1]L =

1 "2L

cT ST2 Q12S2c

[Omega1]L;

[Omega1]V =

1 "2V

[Omega1]V ;

[Omega1]S =

1 "2S

[Omega1]S , (23)

where "2L, "2W , and "2S are the rst elements of the corresponding covariance matrices.

After substitution of the covariance matrices expressed through the weight matrices into Eq. (22) and multiplication it by "2L, the [Psi1](c) takes the form of the sum of three components:

[Psi1](L ,V ,c) =

[Psi1]L(L ,c) + V

[Psi1]V (

V ,c) + S

[Psi1]S(c), (24)

where

[Psi1]L(L ,c) = L L(c)

[Omega1]1L L L(c)

, (25)

is related to lidar-measured data, Eq. (3),

[Psi1]V ( V ,c) = [parenleftBig]

C V Hc[parenrightBig]T

[Omega1]1V

C V Hc[parenrightBig]

, (26)

is related to radiometer-measured data, Eq. (18),

[Psi1]S(c) =

cTST

Q12Sc

, (27)

is related to a priori information, Eq. (21),

V =

"2L

"2V ;

[parenleftBig]

S =

"2L

"22

. (28)

A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC) 1189

Figure 2. Flowchart of the program package.

nc2mdb a program to convert EARLINET standard raw-lidar nc (network common data form) les into mdb (database le used by Microsoft Access) les to process by LIRIC;

Synthesizer a program to average the series of lidar signals, converts the proles to the optimal altitude scale and, then, glues signals (i.e. synthesizes single signal) for the upper and lower troposphere, which were measured with different receiving systems, as well as provides the dead-time correction, i.e. the correction for the nite time resolution of the photo-counting system;

Tropoexport a program to calculate a normalized smoothed lidar signal and its variance, and generates molecular and aerosol atmospheric models; this program aims at implementing modules 2 and 3 of the algorithm.

Finally, the main sub-package ProleRetriever implements the LIRIC inversion procedure. The program ConcentRetriever retrieves proles ck,mV(h) of the aerosol mode concentrations and writes data down to ACCESS database, DB-processed. The module inversion setting & errors modelling generates a set of noise-corrupted input data les by adding white noise and amplitude distortions to the initial lidar signals and perturbing aerosol model parameters retrieved from radiometer measurements in order to provide the error sensitivity analysis. The user can upgrade default

instrumental noise parameters to meet real measurement conditions and technical features of the lidar system; the accuracy of columnar aerosol parameters retrieved from the radiometer measurements (Dubovik and King, 2000; Dubovik et al., 2000) is also taken into account in setting parameters of the module. The program OutputViewer allows viewing the output data and their conversion from mdb-les into other formats.

6 Verication of operability and sensitivity tests

The LRS technique uses the aerosol model that was initially developed in AERONET to describe column-averaged aerosol properties and generalized it to the case of the height-resolved aerosol concentrations. This model assumes that aerosol consists of ne and coarse modes and that both are mixtures of spherical particles and randomly oriented homogeneous spheroids. The advanced T-matrix code (Mishchenko al., 2000, 2002) provides computation of scattering matrices of the aerosol particles. Thus, any optical characteristic of the aerosol layer can be calculated using data of the LRS experiment.

The applicability analysis of the AERONET spheroid model to aerosol particles is beyond the scope of this paper. We only note that this model was validated by the comparison of calculated optical parameters and laboratory measurements of light scattering matrices for mineral dust particles (Volten et al., 2001). Incorporation of the spheroid model into AERONET operational retrieval code has signicantly improved AERONET products when evaluating parameters of coarse non-spherical particles (Cattrall et al., 2005; Dubovik et al., 2006). This model has also been incorporated when processing data from ground-based polarimetric measurements (e.g. Li et al., 2009), lidar sounding data (e.g. Veselovskii et al., 2010; David et al., 2013; Mller et al., 2013), and satellite-base observations (e.g. Levy et al., 2007a, b; Dubovik et al., 2011; Schuster et al., 2012).

6.1 Verication of LIRIC program package: EARLI09 intercomparison experiment

EARLI09 intercomparison experiment was held in May 2009 at Leibniz Institute for Tropospheric Research in Leipzig, Germany (Wandinger et al., 2015). This campaign provided an excellent opportunity to validate the LRS technique for network measurements. The results of the LIRIC data processing for simultaneous measurements by seven lidars of different scientic teams on 25 May 2009 in Leipzig were compared.

Total optical depth distribution (Fig. 3a) and back-trajectory analysis (Fig. 3b) indicates that LRS measurements were carried out during the Saharan dust event in the Leipzig region and the dust was transported in the layer above 2 km.

www.atmos-meas-tech.net/9/1181/2016/ Atmos. Meas. Tech., 9, 11811205, 2016

1190 A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC)

Figure 3. EARLI09 intercomparison experiment: (a) NAAPS Total Optical Depth forecast, 25 May 2009 at 12:00 UTC; (b) 7-day back trajectories ending over Leipzig, Germany at 12:00 UTC on 25 May 2009.

(a) (b)

Altitude, km

0 0.01 0.02 0.03

0 0 0.01 0.02 0.03 PVC, ppb

0 1 2 3

PVC, ppb

mi ms po bh hh bu le

0 1 2 3

rel_dev

Altitude, km

0 0 0.01 0.02 0.03

0 0.01 0.02 0.03

Figure 4. Particle volume concentrations (PVC) proles, ck(h), and estimated deviations retrieved from data of EARLI09 intercomparison campaign, 10:2011:40 UTC, 25 May 2009, Leipzig, Germany, measured in Leipzig by six EARLINET lidars: mi Minsk, ms Mnchen, po Potenza, bh Bilthoven, hh Hamburg, bu Bucharest; (a, d) ne, (b, e) coarse spherical; (c, f) coarse non-spherical; 1 average PVC prole, 2 rms-deviation (rms_dev), 3 relative deviation (rel_dev). Measured data from four lidar channels (355, 532-parallel, 532-cross, 1064 nm) and three-mode aerosol model were used.

Figures 4 and 5 show PVC proles, ck(h), retrieved from lidar data of the different EARLINET teams combined with

PVC, r.m.s.dev; ppb

Figure 5. Identical to Fig. 4 except for data from three lidar channels (355, 532 intensity/parallel polarized component, and 1064 nm) and two mode aerosol model were used. Label le stands for lidar PollyXT of TROPOS, Leipzig: (a, c) ne, (b, d) coarse spherical aerosol mode.

PVC, r.m.s.dev; ppb

the same AERONET information, as well as their root mean square deviations and relative deviations for the two types of input data set, namely, with and without depolarization measurements.

Atmos. Meas. Tech., 9, 11811205, 2016 www.atmos-meas-tech.net/9/1181/2016/

A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC) 1191

It is evident from Figs. 4ac and 5ab that ck(h) proles have similar structure over the troposphere except for the lower layer. The relative deviations increase mainly when values of the aerosol concentration become negligible. The discrepancies are also possible in the near-surface atmospheric layer due to overlap effect (e.g. for the Hamburg lidar system, Fig. 4a).

We explain the discrepancy between ck(h) proles in the near-surface atmosphere by the uncertainty in geometrical overlap factors and the differences in lower-boundary heights of the considered lidar systems. Also some differences in the retrieved concentration proles ck(h) are due to measurement errors and uncertainties in aerosol modelling.

The potential errors in the PVC proles for the specic combined lidar/radiometer experiment were estimated by using the Errors modelling module of the LIRIC package (Fig. 2). Figures 6 and 7 illustrate the sensitivity of the retrieved aerosol concentration proles to the errors of the lidar measurements. The original lidar signals were taken as they measured by Mnchen lidar (curves 4 in Fig. 6) and have been perturbed by adding white noise with different root-mean-square deviations (rms-deviations), j , and have been distorted by multiplying them by the coefcient,

kj (hi) = 1 +

(a) (b)

1 [1/15]

2 [1/-20]

3 [0.5/-40]

4 [0/0]

Lidar signal

0 2 4 6

Altitude, km

1 [1/5]

2 [1/-30]

3 [0.5/-20]

4 [0/0]

Lidar signal

0 2 4 6

0 0 2 4 6

Altitude, km

1 [1/15]

2 [1/-20]

3 [0.5/-40]

4 [0/0]

Altitude, km

1 [1/-5]

2 [1/30]

3 [0.5/20]

4 [0/0]

href hi

href , (29)

where percentage parameter [Delta1]j determines the amount of non-linearity.

In response, the program module generated 12 disturbed lidar signal sets that allowed us to estimate the impact of measurement errors. As an illustration, Figs. 6 and 7 simulate higher errors than typical sets in most EARLINET lidars. Four realizations of the disturbed signals are shown in Fig. 6. Coefcient kj (hi) increases/decreases from referent to start point that results in divergence of the lidar signals in Fig. 6.

PVC proles, ck(h), corresponding to the lidar signals in Fig. 6 and their rms-deviations calculated for full ensembles of input data are shown in Fig. 7. Changes in the PVC proles of the dominant coarse non-spherical mode are shown by the Fig. 7 to be minor (Fig. 7c). Although proles ck(h)

of ne and coarse spherical particles (Fig. 7a and b) are not very stable, they qualitatively retain similarity with the initial distributions.

Figure 8 illustrates the effect of uncertainties in columnar aerosol parameters retrieved from radiometer data. Variations of the columnar aerosol characteristics lead to changes in coefcients a and b of lidar-related Eqs. (14)(17) (Sect. 3.1). Statistical characteristics of aerosol concentration proles retrieved with relative deviation of the parameter #jk,p (effective lidar ratio of the aerosol fraction, see Appendix B) in the range [notdef]20% (the full range) are presented in Fig. 8. Relative

deviation of aerosol concentration prole becomes signi-cant only for small values of the concentration.

Figure 6. Range-corrected normalized lidar signals, L , corrupted with noise and amplitude distortions. Original data are provided by the Mnchen lidar team in the frame of EARLI09 intercomparison campaign, 14:3015:30 UTC, 25 May 2009, Leipzig, Germany:(a) 355 nm, (b) 1064 nm, (c) 532 nm, parallel polarized, (d) 532 nm, cross polarized; 4 original signal, 13 corrupted signals. In square brackets distortion parameters j /[Delta1]j are given.

[Delta1]j 100

(a) (b) (c)

Figure 7. PVC proles, ck(h), and their rms-deviations retrieved in response to disturbed data from of the Mnchen lidar, EARLI09 intercomparison campaign, 14:3015:30 UTC, 25 May 2009, Leipzig, Germany: (a) ne, (b) coarse spherical, (c) coarse non-spherical modes; 4 for the original signal, 13 for disturbed signals; 5 rms-deviation.

6.2 Dependence of retrieved aerosol concentration proles on the content of the input data set

Three types of data set related to different sources of information compose the LIRIC input data-le: three or four

www.atmos-meas-tech.net/9/1181/2016/ Atmos. Meas. Tech., 9, 11811205, 2016

1192 A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC)

(a) (b) (c)

0 1 2

rel_dev

c(h), rms_dev, ppb

0 1 2

rel_dev

0 1 2

Altitude, km

0 0.005

0 0 0.005

0 0 0.015 0.03

c(h), rms_dev, ppb

Figure 8. Variations of PVC proles, ck(h), retrieved with 20 % uncertainties in the aerosol lidar ratios; data of Mnchen lidar,

EARLI09 intercomparison campaign, 14:3015:30 UTC, 25 May 2009, Leipzig, Germany are used; (a) ne, (b) coarse spherical, (c) coarse non-spherical modes; 1 average value, 2 rms-deviation, 3 relative deviation.

measured lidar signals, column-aerosol parameters from radiometer measurements, and a priori smoothness constraints. Two- or three-mode aerosol models are used according to the type of the measured lidar signals. Formally, we deal with redundant input information and, hence, for the accepted aerosol model, the number of input data set can be decreased. Consequently, the signicance of the different information components in retrieval procedure is of interest as well as variations of the retrieved proles, ck(h), in the absence of some input data

As pointed out in Sect. 4, the objective functions of LIRIC regularization algorithm (Eq. 22), consists of a set of terms that implement contribution of different types of input data into the retrieval process. Setting the variance of the specic kind of measurement to a large value implies neglecting the correspondent term in the objective function (Eq. 22) and the elimination of this part of the input data in estimation of the nal aerosol parameters. Program package implements this option and makes allowing one to analyze the contribution of different measured data in the processing procedure of a specic experiment.

Below we shortly examine sensitivity of the retrieved proles, ck(h), to the input data selection for the case of combined lidar/radiometer sounding of the atmospheric aerosol during the last period of Eyjafjallajkull volcano ash transport to the European area in Lille, France, on the 19 May 2010. Air mass back trajectories (Fig. 9) forecasted the possibility of appearance of volcanic ash in the layer between 1300 and 2500 m. The structure of the retrieved proles, ck(h), shown in Fig. 10a agrees well with the forecast. Deviations (by deviations hereinafter we mean standard deviation) (ck(hi)) associated to the proles ck(h) have been calculated by an error modelling procedure similar to the one described in Sect. 6.1.

c(h), rms_dev, ppb

Figure 9. Air-mass back trajectories for Lille at 08:00 UTC, 19 May 2010, (NOAA HYSPLIT model).

(a) (b)

finecoarse/sph coarse/nsph rms_dev(fine) rms_dev(coarse/sph) rms_dev(coarse/nsph)

D(1)

D(2) rms_dev(1) rms_dev(2)

Altitude, km

0 0.01 0.02

c(h), ppb

0 0 0.05 0.1 0.15

Par. Dep. Ratio

Figure 10. (a), PVC proles, ck(h), of the ne, course-spherical (coarse/sph) and coarse-nonspherical (coarse/nsph) aerosol modes, and their rms-deviations (rms_dev(ne), rms_dev(coarse/sph), and rms_dev(coarse/nsph)); (b), particle depolarization ratio, D(1) and(2), and their rms-deviations, rms_dev(1) and rms_dev(2). Proles were retrieved from the data measured in Lille, 19 May 2010, 09:1709:58 UTC. Proles D(1) and rms_dev(1) are the results of the direct calculation of depolarization ratio and their rms-deviations from lidar measurements, as well as D(2) and rms_ev(2) were calculated from retrieved aerosol mode concentrations, ck(h).

A mixture of spherical and non-spherical particles constitutes the aerosol layer at the height of about 2000 m. The prole of particle depolarization ratio at 532 nm and its deviation have been calculated from the retrieved aerosol mode

Atmos. Meas. Tech., 9, 11811205, 2016 www.atmos-meas-tech.net/9/1181/2016/

A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC) 1193

(a) (b) (c)

355

532

1064

532-cross

Cv Original

fine(3) coarse/sph coarse/nsph coarse(3) fine(2) coarse(2)

355

532

1064

Cv Original

Altitude, km

0 0 0.01 0.02 0.03 c(h), ppb

0 0 0.01 0.02 0.03

0 0 0.01 0.02 0.03 c(h), ppb

0 0 0.01 0.02 0.03

0 0 0.02 0.04 c(h), ppb

Figure 12. Comparison of PVC proles, ck(h), for the two- and three-mode aerosol models (a), and variations of concentration proles, ck(h), for ne (b) and coarse (c) aerosol modes of the two-mode aerosol model in response to elimination of different parts of input information. In Fig. 12a tags ne(2) and coarse(2) denote ne and coarse modes of two-mode aerosol model. Tags ne(2), coarse/sph, coarse/nsph and coarse(3) denote ne, coarse spherical, course non-spherical and total course mode of three-mode aerosol model, correspondingly. In Fig. 12b and c tag Original means complete set of input data; tag 355 (or 532, 1064) denotes that the lidar signal at 355 nm (or 532, 1064) wavelength is excluded; tag CV denotes that columnar volume concentrations of aerosol modes are excluded. Lille, 08:00 UTC, 19 May 2010.

for two-mode aerosol model, coarse (2), and the sum of two coarse components for three-mode aerosol model, coarse (3), are similar in shape but quantitatively are a bit different. The column concentrations of the course (2) and (3) modes are equal.

The curves in Figs. 12b and c show the deviations of the concentration proles, ck(h), for the two-mode aerosol model after reduction of the input data set. Deviations of ck(h) proles are rather similar to those for the three-mode aerosol model in Fig. 11. Deviations of ne-mode concentration prole are small, even if any single sub-set of input data is eliminated. Coarse-mode concentration proles preserve original forms when one of the lidar signals at the 355 or 532 nm wavelength is excluded from the processing procedure.

Generally, for measurement conditions that characterize the experiment under discussion, two-wavelength lidar sounding (at 355 and 1064 or at 532 and 1064 nm) combined with radiometer measurement provides retrieving concentration proles of ne and coarse aerosol modes for two-mode aerosol model.

7 Discussion and conclusions

The active process of dissemination of the LIRIC in EARLINET started in 2012. Nowadays, 11 EARLINET teams participate in implementation of LRS technique (see Fig. 13).

www.atmos-meas-tech.net/9/1181/2016/ Atmos. Meas. Tech., 9, 11811205, 2016

0 0.01 0.02 0.03

c(h), ppb

Figure 11. Variation of aerosol concentration proles, ck(h), for ne (a), coarse spherical (b) and coarse non-spherical (c) aerosol modes in response to elimination of different parts of input information. Tag Original denotes complete set of input data; tag 355 (or 532, 1064, 532-cross) denotes that lidar signal at 355 nm (or 532, 1064, 532-cross) wavelength is excluded; tag CV denotes that columnar volume concentrations of aerosol modes are excluded. Lille, 08:00 UTC, 19 May 2010.

concentrations, ck(h). The proles are shown in Fig. 10b, curves D(2) and rms_dev(2). The results of the direct calculation of depolarization ratio and their deviations from lidar measurements are presented by curves D(1) and rms_dev(1). It should be noted that the lidar measurements included additional calibration measurement that was not used by the retrieval procedure. Proles D(1) and D(2) show rather close agreement in magnitude and vertical structure that could conrm the efciency of the aerosol modelling used in this study.

The curves in Fig. 11 show the deviations in the retrieved concentration proles, ck(h), after elimination one of the lidar signals or columnar volume concentrations of aerosol modes, CV , from the input data set. As can be seen from Fig. 11, the concentration prole of the ne-particle mode undergoes minor changes upon elimination of a single lidar signal or the elimination of columnar volume concentrations. This could imply that our experiment initially included redundant input information, with respect to the ne-mode concentration. On the other hand, concentrations of coarse modes are sensitive to input information. Thus, lidar data at 1064 nm wavelength plays a crucial role in the retrieval of the coarse spherical mode. In the same manner, lidar depolarization measurement is the key factor in the retrieval of the coarse spheroid particle mode. Evaluations of columnar volume concentrations from radiometer measurement are necessary for all cases.

Figure 12a shows concentration proles, ck(h), which were retrieved for two- and three-mode aerosol models and characterized the aerosol layer in the same LRS experiment. The ne-mode concentration proles for two aerosol models are practically coincident. Proles ck(h) of coarse modes

1194 A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC)

Figure 13. Map of the EARLINET stations (red dots). Green dots indicate the stations where LIRIC program package has been implemented.

New scientic teams beyond EARLINET join the LIRIC user group. The detailed description of LIRIC algorithm and software in this paper should contribute to the effective implementation of the LRS technique by advanced users.

Retrieval of the aerosol parameters from the LRS measurements is an ill-posed inverse problem, and its solution should be tested on stability to the measurement errors and variations of the regularization parameters, which are set by the module Inversion setting and errors moduling of the software package (Fig. 2). Results of the EARLI09 intercom-parison experiment presented in Sect. 6.1 demonstrate rather small scatter in ck(h) proles that were retrieved from the data of different lidar systems with signicantly corrupted input lidar signals and big uncertainties of the aerosol lidar ratio. This scatter is characterized by standard deviations of 5 20 % of the maximum aerosol layer concentration. Increase in ck(h) deviation in the bottom layer results from uncertainties of the overlap function of the lidar systems.

The uncertainties in the retrieved aerosol parameters for different aerosol types, aerosol loads, overlap characteristics of the lidar systems and regularization parameters that are dened by the LIRIC operator were evaluated by Granados-Muoz et al. (2014). The analysis covered combined lidar and radiometer measurements that were carried out during dust, smoke, and anthropogenic pollution events. This analysis mostly supports our conclusions on the stability of LIRIC solutions that retrieve basic aerosol features even under signicant measurement errors. In particular, variations of the regularization parameters within one order interval from the original set lead to minor deviations of the retrieved ck(h)

proles. Usually, it is unnecessary to change recommended utility regularization parameters while homogeneous input data sets are processed. The requirements for pre-processing

lidar signals along with the set of recommended regularization parameters are provided in the LIRIC user guide. However, the utility parameters for error modelling menu should be dened by the LIRIC user with regard to the specic lidar system.

The requirement of having possibly minimal full overlap height of lidar sensing is an important technical problem for LRS measurements, because the near-surface aerosol layer contributes strongly to the radiometric data. In the absence of lidar data, the surface aerosol layer is assumed to be homogeneous in the LIRIC aerosol modelling. Obviously, aerosol parameters can vary within the near-surface layer resulting in signicant uncertainties in the LIRIC product, especially when the lidar dead zone becomes comparable to the boundary-layer thickness. The effective solution of this problem is the set-up of a double lidar receiving block with special near-range channels for the detection of near-ground aerosol.

The analysis of the aerosol parameters that are retrieved from the incomplete sets of lidar data in Sect. 6.2 supports the possibility to use LIRIC for processing data of two-wavelength lidar systems. Aerosol sounding by two-wavelength lidars, usually at 532 and 1064 nm wavelengths, is a widespread practice in atmospheric investigations. Simulation results in Sect. 6.2 show the possibility to retrieve ck(h) for two-mode aerosol model. The uncertainties of such evaluated ck(h) are expected to surpass ones of three-wavelength lidar sounding.

LIRIC implementation for the special lidar data set (532-cross, 532-parallel and 1064 nm) for retrieving parameters of the three-mode aerosol model is of interest for the satellite lidar CALIOP that provides similar lidar data (Winker et al, 2006).

Since the beginning of LIRIC dissemination in EARLINET community, experimental works on the validation of the LIRIC product for different aerosol types have being carried out. Comparisons of aerosol backscatter coefcients and depolarization ratios directly derived from lidar data against similar characteristics calculated from the aerosol optical and microphysical parameters retrieved by LIRIC (e.g. Tsekeri et al., 2012, 2013; Wagner et al., 2013; Kokkalis et al., 2013;Granados-Muoz et al., 2014) as well as LIRIC against modelled or airborne in situ measured proles of aerosol mode concentrations (e.g. Kokkalis et al., 2012, 2013; Nemuc et al., 2013) have shown reasonable agreement.

The LIRIC concentration proles of aerosol fractions during dust and volcano ash events have been compared with those for spherical and non-spherical particles derived from polarization measurements using the POLIPHON technique (e.g. Wagner et al., 2013; Nemuc et al., 2013, Papayannis et al., 2014). In spite of the noticeable difference between the aerosol models and independent processing algorithms, the retrieved aerosol concentration proles have proved to be similar. This is quite natural because both approaches use the

Atmos. Meas. Tech., 9, 11811205, 2016 www.atmos-meas-tech.net/9/1181/2016/

A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC) 1195

depolarization of backscatter signal to distinguish between spherical and non-spherical particles.

The number of aerosol studies using LIRIC algorithm increases. They focus on the investigation of the dynamics of aerosol microstructure during transport of air masses polluted by dust (e.g. Chaikovsky et al., 2010b; Tsekeri et al., 2013; Binietoglou et al., 2015; Granados-Muoz et al., 2015a), re smoke (e.g. Chaikovsky et al., 2010b; Granados-Muoz et al., 2015b), volcano ash (Kokkalis et al., 2013) and anthropogenic pollution (Granados-Muoz et al., 2014).LIRIC has become a tool for validation of the modelling

of aerosol transport in atmosphere (Binietoglou et al., 2015; Granados-Muoz et al., 2015b). EARLINET teams form the data-base of the results of combined lidar and radiometer sounding.

The list of lidar teams that take advantage of LIRIC is still expanding. The LIRIC software package is open and distributed both within the EARLINET community and beyond it. The EARLINET teams provide continuous improvement of the software and cooperate on the implementation of the LRS measurements at new sites.

www.atmos-meas-tech.net/9/1181/2016/ Atmos. Meas. Tech., 9, 11811205, 2016

1196 A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC)

Appendix A: General equation for received lidar signal

Using general formula for received lidar signal instead of Eqs. (4), (6), and (8) allows us to derive compact and explicit expression for the covariance matrices, [Omega1]L, and regularizing term,

[Psi1]L(L ,c) (Sect. 4).

We will use the utility function

jpj,u =

Therefore, the related to the lidar objective function,

[Psi1]L(L ,c), (Eq. 25), is given by the equation:

[Psi1]L(L ,c) =Xj Xi

[Delta1]hi

[Omega1]Lj (i,i)

0[parenleftBigg]L j,i [summationtext]

ck(hi)bk,pj ( j ) +

j p,2[notdef]

ck(hi)bk,pj ( j )

[parenrightBigg]

r ( j ,pj ,href) R

1...if .....pj = u

0...f ...pj [negationslash]=u

efj ( j ,pj ,href)

(A8)

, (A1)

along with the following denitions of combinations of aerosol and molecular optical parameters in Eqs. (4)(9):

efa( j ,pj ,h) = a,pj ( j ,h)

+ jpj,2[notdef] a,3( j ,h)

[parenrightbig]

=Xk

ck(h)bk,pj ( j )

+ jpj,2[notdef]

ck(h)bk,pj ( j )

exp 2

ck(hi)ak( j )[Delta1]hi

[parenrightBigg]

i 2 1,...I.

[Psi1]V ( V ,c), which brings radiometer data

into the processing procedure can be expressed as follows:

[Psi1]V ( V ,c) =

Xk1

Equation (26),

[Omega1]V (k,k)

[parenrightBigg]

(A2)

C V

Xi ck(hi)[notdef][Delta1]hi[notdef]

!2. (A9)

Calculation of the smoothness part of the objective function is described in details in Dubovik, 2004; and Dubovik et al., 2011.

Appendix B: Evaluation of covariance matrix [Omega1]L

The covariance matrixes, [Omega1]L, [Omega1]V , and [Omega1]2, dened in Sect. 4 characterize uncertainties of the complex input vector,

L , V , 0 , where 0 is zero vector that is dened to formalize a priori smoothness restrictions on concentration proles (e.g. Dubovik, 2004). These matrices determine the weights of different parts of input information through the minimization procedure of the objective function (Eq. 22).

In our case the measure of the smoothness for concentration proles, ck(hi), should be chosen as a priori evaluated parameters. Aerosol columnar volume concentrations, V , and variances, [Omega1]V (k,k), are the parts of input radiometer data. Thus, only evaluation of covariance matrix, [Omega1]L, is to be done.

The assumption of independent normal distribution for variations of lidar vector, L , at different heights implies the diagonal covariance matrix. The non-zero diagonal elements, [Omega1]Lj (hi,hi), of the covariance matrix are the variances of differences between the components, L j,i, of the lidar vector and the appropriate modelled function,

Lj ck,pj , j ,hi

, in Eq. (A7).

efr( j ,pj ,h) =[parenleftbigg]

jpj,2(pj )

[notdef] 1 ~ + 1

[parenrightbigg]

+ 1 + jpj,3~ [parenrightbigg]

r( j ,h) (A3)

ef( j ,pj ,h) = efa( j ,pj ,h)

+ efr( j ,pj ,h) (A4)

Refj( j ,pj ,h) =

efa( j ,pj ,h) + efr( j ,pj ,h) efr( j ,pj ,h)

(A5)

a( j ,h)dh. (A6)

This permits Eqs. (4), (5) and (8) to be written in general form:

Lj pj , j ,h

a( j ,h,href) =

href

[integraldisplay]h

exp 2a( j ,h,href)

efr ( j ,pj ,href) Refj( j ,pj ,href)

. (A7)

ef j ,pj ,h

[parenrightbig]

Atmos. Meas. Tech., 9, 11811205, 2016 www.atmos-meas-tech.net/9/1181/2016/

A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC) 1197

Given Eqs. (2), (3) and (A1)(A-8), the elements of vector, [Delta1]Lj , are dened by the following:

[Delta1]Lj (hi) = L j,i Lj pj , j ,hi [parenrightbig]

S j ,hi

[parenrightbig]

S j ,href

exp(2r( j ,hi,href))

[parenrightbig]

of a priori dened optical characteristics. We aim at approximate estimation of [Omega1]Lj (hi,hi) at the preprocessing stage without involving of retrieved parameters. This feedback-free approach greatly simplies the structure of the inversion algorithm.

Uncertainties of the optical parameters

The term 2 ef

ef j ,pj ,hi

exp 2a( j ,hi,href)

efr ( j ,pj ,href) Refj( j ,pj ,href)

. (B1)

Using the nite differences technique (e.g. Russell et al., 1979) one can expand [Delta1]Lj (hi) in Taylor series, and then neglect all the terms of the second or higher order. As a result, variation ([Delta1]Lj (hi)) can be expressed as a function of variations related with the input parameters, S j ,hi

[parenrightbig][parenrightbig],

[parenrightbig]

/ ef

2 in Eq. (B3) is the relative variance of the total backscatter coefcient. It can be transformed into the sum of relative variances of aerosol and molecular backscatter coefcients:

2 ef( j ,pj ,hi)

2 efr( j ,pj ,hi)

ef( j ,pj ,hi)

efr( j ,pj ,hi)

2 1

( ef( j ,pj ,hi)), (a( j ,hi,href), and (r( j ,hi,href)):

([Delta1]Lj (hi)) = 2L j,i (r( j ,hi,href))

+L j,i

[parenleftBig]

Refj( j ,pj ,hi)

2 efa( j ,pj ,hi)

S j ,hi

efa( j ,pj ,hi)

[parenrightbig]

[parenrightbig][parenrightbig]

ef j ,pj ,hi

[parenrightbig]

[parenleftBig]

exp 2a( j ,hi,href)

efr ( j ,pj ,href) Refj( j ,pj ,href)

ef j ,pj ,hi

[parenrightbig][parenrightbig]

Refj( j ,pj ,hi) 1[parenrightBig]2

[parenleftBig]

Refj( j ,pj ,hi)

ef j ,pj ,hi

[parenrightbig]

2 ef j ,pj ,hi

[parenrightbig]

exp 2a( j ,hi,href)

efr ( j ,pj ,href) Refj( j ,pj ,href)

(a( j ,hi,href)

L j,i[parenleftbigg]

2 . (B4)

The International Standard Atmosphere ISO 2533 and seasonal latitudinal changed model CIRA (Committee on Space Research (COSPAR), 2006; Fleming et al., 1988), as well as measurements by radiosondes are applied in LIRIC for the calculation of molecular optical parameters. The relative variance of calculated molecular backscatter coefcient

21 =

S j (hi)

ef( j ,pj ,hi) ef( j ,pj ,hi)

2 efr( j ,pj ,hi)

2 (B5) is assumed to be a constant and its value can be reduced to 1 = 0.01 (e.g. Russell et al., 1979) if data of coordinated

radiosonde measurements are available.

The aerosol backscatter coefcients, efa( j ,pj ,hi), are estimated by using Eqs. (10)(17). Uncertainties of efa( j ,pj ,hi) basically follow from estimation errors of the coefcient b( ,j,pj ) in Eqs. (15)(17) that can be written by the equation:

b(j,pj ,k) =

2 r( j ,hi,href)

[parenrightbig]

2 a( j ,hi,href)

efr( j ,pj ,hi)

[parenrightbig][parenrightbigg]

. (B2)

Under the assumption of independent variations of different parameters, the variance [Omega1]L(hi,hn) is expressed as follows

[Omega1]L(hi,hi) =

([Delta1]Lj (hi)) ([Delta1]Lj (hi))

[angbracketrightbig]

= L 2j,n[parenleftbigg]

P j,i

[parenrightBig]

P j,i

2 ef( j ,pj ,hi)

1 E k j

, (B6)

ef( j ,pj ,hi)

CVk

+4 2 r( j ,hi,href)

[parenrightbig]

+4 2 a( j ,hi,href)

[parenrightbig][parenrightbigg]

where

#jk,p =

14 [pi1]k( j )Ajk,p, (B7)

, (B3)

where [angbracketleft]...[angbracketright] denotes ensemble averaging over measurement

realizations, and P j,i = P j(hi).

The terms in the large round parentheses in Eq. (B3) determine contributions of measurement errors and uncertainties

1,1( j , = 180 ) if pj = 1

P k

1,1( j , = 180 ) P

Ajk,p =

8 >

k2,2( j , = 180 )

2 pj = 2

P k

1,1( j , = 180 ) + P

. (B8)

www.atmos-meas-tech.net/9/1181/2016/ Atmos. Meas. Tech., 9, 11811205, 2016

k2,2( j , = 180 )

2 pj = 3

1198 A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC)

Parameter #jk,p = ka( j ,hi)/ ef,ka( j ,hi) in Eq. (B8) is

the extinction-to-backscatter ratio or lidar ratio of the k-aerosol mode.

Parameters #jk,p are retrieved from the data of radiometric direct sun and almucantar measurements that are usually performed with the maximum scattering angle less than 150 . The range of the scattering angles decreases as the sun zenith angle decreases. Retrieval of optical parameters in the backscatter direction, in a certain sense, is an extrapolation procedure out of the measured range with possible increasing of estimation uncertainties. One assumes that the errors of the estimation of #jk,p are the main reason of the incorrect calculation of backscatter coefcients efa( j ,pj ,hi) and introduces parameter 2 for characterization of the standard deviation of coefcients 1/#jk,p in LRS measurements.

Thus, Eq. (B4) is transformed to

2 ef( j ,pj ,hi)

[parenrightbig]

ef( j ,pj ,hi)

2 Refj( j ,pj ,hi)

+ 22

Refj( j ,pj ,hi) 1[parenrightBig]2

Refj( j ,pj ,hi)

2 . (B9)

The backscatter ratio Refj( j ,pj ,hi) in Eq. (B9) under assumption [notdef] = 0 is approximately calculated at the pre

processing stage using the Klett algorithm (Klett, 1981).

Basically, the variance of aerosol optical thickness, 2 a( j ,hi,href)

, arises from altitude variations of aerosol modes that are not assumed by the aerosol model. Relative error of a( j ,hi,href) is zero at the reference point and is equal to 23 at the start pointh1, where 3 is close to the error of AOT calculation from radiometer measurements. Thus, the following approximation is used in the LIRIC algorithm:

2 a( j ,hi,href)

[parenrightbig]

!2j

P j,n + B j[parenrightBig]2

(P j,n)2

= 232a( j ,hi,href). (B10)

Term 2 r( j ,hi,href)

[parenrightbig]

G j

2 +q2j

in Eq. (B3) denotes the variance of molecular optical thickness of the atmospheric layer (hn,href). Only long-scale or systematic deviations of molecular density contribute to the variance 2 r( j ,hi,hN)

P j,n + B j[parenrightBig]

+ A(2M + 1)(P j,n)2

U j

, (B12)

where ! is the coefcient of nonlinearity, G j is the amplitude of electrical noise, q2j is the coefcient characterizing the power of the Schottky noise, U j is the amplitude of synchronized noise;

Similar to Eq. (B10),

(P j,n)2

2 r( j ,hi,href)

= 242r( j ,hi,href). (B11)

Measurement errors

Optical signals, detected by the lidar data acquisition system consists of backscatter P j,i and background B j components.

A suitable algorithm for estimating the measurement errors is described by Slesar et al. (2013, 2015). Regardless of the type of the photo-receiving sensor, three factors determine the measurement errors:

non-linearity of the recording channel, which consist of nonlinearity of the photodetector and electronic units;

non-synchronous noise (non-correlated with the sounding pulse);

synchronized noise (correlated with the sounding pulse).

Non-linearity of a receiving channel basically originates from saturation of an output signal at high incident light because of photo-sensor or electronic unit limitations. Likewise, deviations of an amplier gain cause linear distortions of the detecting signal within the working range of photo-receiving module.

Basic difference between two types of noise is that nonsynchronous noise can be reduced by accumulation of input signals or by decreasing frequency bandwidth of the receiving channel, while this method is ineffective for synchronized noise. The main type of the non-synchronous noise is the Schottky noise. Synchronous noise is basically caused by the interference of the electrical impulses from the laser power supply, synchronous with the sounding optical pulse. It is predominantly a low-frequency noise, and acceptable limitation of the frequency band of the photo-receiving channel does not lead to its decline.

We assume that the accumulation of the receiving lidar signal with A sounding pulses and the averaging of the lidar signal over 2M +1 bins are carried out at the measurement and

pre-processing stages.

Summing up the contributions of the noise components, one can write the following expression for the variances of the receiving analog and photon-counting signals (Slesar et al., 2013, 2015):

for the analog channel:

2(P j,i)

(P j,i)2 =

[parenrightbig]

Atmos. Meas. Tech., 9, 11811205, 2016 www.atmos-meas-tech.net/9/1181/2016/

A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC) 1199

for the counting channel:

2(N j,i)

(N j,i)2 =

!2j

N j,i + Nj,B[parenrightBig]2

(N j,i)2

N j,G

2 +N j,i + Nj,B A(2M + 1)(N j,i)2

N j,U

Cond { }

, (B13)

where N j,i is the detected lidar signal, Nj,B is the background signal, Nj,G is the external non-synchronous noise, and N j,U is synchronized noise.

Parameters !j , G j, N j,G,qj ,U j, N j,G, N j,U for specic photo-receiving module can be evaluated on a dedicated test bench by means of special calibration procedures (Slesar et al., 2013, 2015).

Appendix C: Details of inversion procedure

One can understand intuitively that optical parameters of aerosol modes, which constitute the aerosol model (see Sect. 1.1), should be different to allow retrieving aerosol mode concentrations by means of algorithm described in Sect. 4. More correct denition of this is that there should not be a linear relationship between the sets of coefcients

{ak,bk[notdef] which dene optical characteristics of kth aerosol

mode. This conclusion results from the linear approximation of the Eq. (2). It means that we seek the solution, c(hi), only from data of multi-wavelength lidar sounding. The linear least squares solution of Eq. (3) can be written as

c =

(N j,i)2

Figure C1. Cumulative distribution functions (CDFs) of parameter Cond

Uk,k

[bracerightbig](condition number) calculated from radiometer data of the AERONET station in Minsk for two- and three-fraction aerosol models, Model 2 and Model 3, respectively.

UL = KTL[Omega1]1LKL with dimensions KI [notdef]KI takes the block-

diagonal form

UL =

1 "21

Uk,k ...

...

... 1

, (C3)

"2i

Uk,k ...

...

... 1

"2I

Uk,k

where matrix Uk,k , (k 2 1,...,K), does not depend on

the superscript i. For 3-mode aerosol model (K = 3) and

4-channel lidar measurements (J = 4) matrix Uk,k can be

written

Uk,k =

Pj b 2

[j,pj ,1]

Pj b[j,pj ,1]b[j,pj ,2]

Pj b[j,pj ,1]b[j,pj ,3]

KTL[Omega1]1LKL

Pj b[j,pj ,1]b[j,pj ,2]

Pj b 2

[j,pj ,2]

Pj b[j,pj ,2]b[j,pj ,3]

1KTL[Omega1]1LL , (C1)

where KL is the Jacobi matrix of the rst partial derivatives

{KL[notdef]x,y = Lx/ cy[notdef]c. The following denitions are used in

Eq. (C1) for measured vector, L and state vector, c, with

dimensions J I [notdef] 1 and KI [notdef] 1 , correspondingly:

L j,i =

(C4)

Pj b[j,pj ,1]b[j,pj ,3]

Pj b[j,pj ,2]b[j,pj ,3]

Pj b 2

[j,pj ,3]

Thus, results of the retrieval depend on the specics of matrix Uk,k . The well-conditioned matrix Uk,k provides suitable solution of Eq. (3). On the analogy with Veselovskii et al. (2005), the eigenvalue decomposition technique has been used to evaluate the condition number of matrix Uk,k

Cond

Uk,k

[bracerightbig]

= [notdef] max[notdef]/[notdef] min[notdef], (C5)

where max and min are the maximum and minimum eigenvalues of matrix,Uk,k respectively. Parameter

qCond

Uk,k

[bracerightbig]is a coefcient of increasing relative error of ck,i as compared to the relative error of L n,j estimation (Trefethen and Bau, 1997).

The data of radiometric measurements in Minsk during 20022010 were used to calculate the parameters Cond

Uk,k

[bracerightbig]for the aerosol models with two and three

www.atmos-meas-tech.net/9/1181/2016/ Atmos. Meas. Tech., 9, 11811205, 2016

8 >

L 1(h1) L 2(h1)

...

L j(hi)

...

L J (hI )

, ck,i =

8 >

c1(h1) c2(h1)

... ck(hi)

... cK(hI )

. (C2)

The formula (C1) is valid if det

UL = KTL[Omega1]1LKL[parenrightBig] [negationslash]=0.We use additional requirements that optical thickness of the aerosol layer is small, and the variances of the measured errors, ("2L,1,"2L,i,"2L,I ), do not depend on hi. So the matrix

1200 A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC)

aerosol fractions (three and four measuring channels, correspondingly). The cumulative distribution functions (CDFs) of parameter Cond

Uk,k

[bracerightbig]is shown in Fig. (C1).

Matrix Uk,k is sufciently well conditioned for the two-fraction aerosol model, and solution (C1) is applicable for the calculation of aerosol mode concentrations. In the case of the three-fractional aerosol model, parameters Cond

Uk,k

increase approximately by 10, and the matrix Uk,k becomes ill- conditioned. In a case such as this, we have to involve the Eq. (18) in retrieving procedure, i.e. to use information on parameter V from radiometric measurements. With our

denitions, matrix H in Eq. (18) is written as

HK[notdef]IK = (C6)

[Delta1]h1 . . . 0 [Delta1]h2 . . . 0 . . . [Delta1]hI . . . 0

0 [Delta1]h1 0 0 . . . 0 . . . 0 . . . 0. . . . . . . . . . . . . . . . . . . . . . . . . . . . . .0 . . . [Delta1]h1 0 . . . [Delta1]h2 . . . 0 . . . [Delta1]hI

Finally, a priori smoothness restrictions are used as the additional factor for regularizing the ill-posed problem solution.

Atmos. Meas. Tech., 9, 11811205, 2016 www.atmos-meas-tech.net/9/1181/2016/

A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC) 1201

Acknowledgements. The nancial support by the European Unions Horizon 2020 research and innovation programme (ACTRIS-2, grant agreement no. 654109) is gratefully acknowledged. The background of LIRIC algorithm and software was developed under the ACTRIS Research Infrastructure project, grant agreement no. 262254, within the European Union Seventh Framework Programme, which nancial support is gratefully acknowledged.

The authors are very thankful to AERONET community for establishing and maintaining the radiometer network and its information system that provides the excellent tools for aerosol studying and became the base for our work. The authors also thank TomL. Kucsera (GESTAR/USRA) for back trajectories available at the http://aeronet.gsfc.nasa.gov

Web End =http://aeronet.gsfc.nasa.gov website.

The authors gratefully acknowledge the Naval Research Laboratory (NRL) for the provision of the NAAPS results and the NOAA Air Resources Laboratory (ARL) for the provision of the HYSPLIT transport and dispersion model used in this publication.I. Binietoglou received funding from the European Unions Seventh Framework Programme for research, technological development and demonstration under the grant agreement no. 289923 ITARS.

Edited by: N. Sugimoto

References

Ansmann, A., Baars, H., Tesche, M., Mller, D., Althausen, D., Engelmann, R., Pauliquevis, T., and Artaxo, P.: Dust and smoke transport from Africa to South America: Lidar proling over Cape Verde and the Amazon rainforest. Geophys. Res. Lett., 36, L11802, doi:http://dx.doi.org/10.1029/2009GL037923

Web End =10.1029/2009GL037923 http://dx.doi.org/10.1029/2009GL037923

Web End = , 2009.

Ansmann, A., Tesche, M., Gro, S., Freudenthaler, V., Seifert, P.,

Hiebsch, A., Schmidt, J., Wandinger, U., Mattis, I., Mller, D., and Wiegner, M.: The 16 April 2010 major volcanic ash plume over central Europe: EARLINET lidar and AERONET photometer observations at Leipzig and Munich, Germany, Geophys. Res. Lett., 37, L13810, doi:http://dx.doi.org/10.1029/2010GL043809

Web End =10.1029/2010GL043809 http://dx.doi.org/10.1029/2010GL043809

Web End = , 2010. Ansmann, A., Tesche, M., Seifert, P., Gro, S., Freudenthaler, V.,

Apituley, A., Wilson, K. M., Serikov, I., Linn, H., Heinold, B., Hiebsch, A., Schnell, F., Schmidt, J., Mattis, I., Wandinger, U., and Wiegner, M.: Ash and ne mode particle mass proles from EARLINET/AERONET observations over central Europe after the eruptions of the Eyjafjallajokull volcano in 2010, J. Geophys. Res., 116, D00U02, doi:http://dx.doi.org/10.1029/2010JD015567

Web End =10.1029/2010JD015567 http://dx.doi.org/10.1029/2010JD015567

Web End = , 2011. Ansmann, A., Seifert, P., Tesche, M., and Wandinger, U.: Proling of ne and coarse particle mass: case studies of Saharan dust and Eyjafjallajkull/Grimsvtn volcanic plumes, Atmos. Chem. Phys., 12, 93999415, doi:http://dx.doi.org/10.5194/acp-12-9399-2012

Web End =10.5194/acp-12-9399-2012 http://dx.doi.org/10.5194/acp-12-9399-2012

Web End = , 2012. Antua, J. C., Landulfo, E., Clemesha, B., Zaratti, F., Quel, E.,

Bastidas, A., Estevan, R., and Barja, B.: Lidar community in Latin America: a decade of challenges and successes: in Reviewed and Revised Papers Presented at the 26th International Laser Radar Conference (ILRC 2012), 2529 June 2012, Porto Heli, Greece, 323326, 2012.

Biele, J., Beyerle, G. and Baumgarten G.: Polarization lidar: corrections of instrumental effects, Opt. Express 7, 427435, 2000. Binietoglou, I., Basart, S., Alados-Arboledas, L., Amiridis, V., Argyrouli, A., Baars, H., Baldasano, J. M., Balis, D., Belegante,

L., Bravo-Aranda, J. A., Burlizzi, P., Carrasco, V., Chaikovsky,A., Comern, A., DAmico, G., Filioglou, M., Granados-Muoz,M. J., Guerrero-Rascado, J. L., Ilic, L., Kokkalis, P., Maurizi,A., Mona, L., Monti, F., Muoz-Porcar, C., Nicolae, D., Papayannis, A., Pappalardo, G., Pejanovic, G., Pereira, S. N., Perrone, M. R., Pietruczuk, A., Posyniak, M., Rocadenbosch, F., Rodrguez-Gmez, A., Sicard, M., Siomos, N., Szkop, A., Terradellas, E., Tsekeri, A., Vukovic, A., Wandinger, U., and Wagner, J.: A methodology for investigating dust model performance using synergistic EARLINET/AERONET dust concentration retrievals, Atmos. Meas. Tech., 8, 35773600, doi:http://dx.doi.org/10.5194/amt-8-3577-2015

Web End =10.5194/amt-8- http://dx.doi.org/10.5194/amt-8-3577-2015

Web End =3577-2015 , 2015.

Bsenberg, J. and Hoff, R. M.: Plan for the implementation of the

GAW Aerosol Lidar Observation Network GALION, World Meteorological Report #178, WMO Geneva, Switzerland, 2007. Bsenberg, J., Ansmann, A., Baldasano, J. M., Balis, D., Bck-mann, Ch., Calpini, B., Chaikovsky, A., Flamant, P., Hgrd,A., Mitev, V., Papayannis, A., Pelon, J., Resendes, D., Schneider J., Spinelli, N., Trickl, T., Vaughan, G., Visconti, G., and Wiegner, M.: EARLINET-A European Aerosol Research Lidar Network, Advances in Laser Remote sensing, in: Selected papers 20th Int. Laser Radar Conference (ILRC), 1014 July 2000, Vichy, France, 155158, 2000.

Cattrall, C., Reagan, J., Thome, K., and Dubovik, O.: Variability of aerosol and spectral lidar and extinction ratios of key aerosol types derived from selected aerosol robotic network locations, J. Geophys. Res, 110, D10S11, doi:http://dx.doi.org/10.1029/2004JD005124

Web End =10.1029/2004JD005124 http://dx.doi.org/10.1029/2004JD005124

Web End = , 2005. Chaikovsky, A., Bril, A., Barun, V., Dubovik, O., Holben, B.,

Goloub, P., and Sobolewski, P.: Methodology and sample results of retrieving aerosol parameters by combined multi-wavelength lidar and Sun-sky scanning measurements, Proc. SPIE 5397, Tenth Joint International Symposium on Atmospheric and Ocean Optics/Atmospheric Physics. Part II: Laser Sensing and Atmospheric Physics, 23 February 2004, 146, doi:http://dx.doi.org/10.1117/12.548588

Web End =10.1117/12.548588 http://dx.doi.org/10.1117/12.548588

Web End = , 2004a.

Chaikovsky, A., Bril, A., Barun, V., Dubovik, O., Holben, B., Thompson, A., Goloub, P., and Sobolewski P.: Studying altitude proles of atmospheric aerosol parameters by combined multi-wavelength lidar and sun sky radiance measurements, Reviewed and Revised papers presented at the 22nd International Laser Radar Conference (ILRC 2004), 1216 July 2004, Matera, Italy, 345348, 2004b.

Chaikovsky, A., Bril, A., Dubovik, O., Holben, B., Thompson, A., Goloub, P., ONeill, N., Sobolewski, P., Bsenberg, J., Ansmann,A., Wandinger, U., and Mattis, I.: CIMEL and multiwavelength lidar measurements for troposphere aerosol altitude distributions investigation, long-range transfer monitoring and regional ecological problems solution: eld validation of retrieval techniques, Optica Pura y Aplicada, 37, 32413246, 2004c.

Chaikovsky, A., Kabashnikov, V., Germenchuk, M., Goloub, P.,

Dubovik, O., Zhukova, O., Ivanov, A., Kozeruk, B., Korol, Y., Lopatsin, A., Asipenka, F., and Utochkina, S.: Monitoring of transboundary pollution transfer in the atmosphere with use of remote sounding systems and the global and regional measuring network data: I Operational monitoring of particular matter in Belarus region, Nature Resour., 1, 95108, 2010a (in Russian). Chaikovsky A., Dubovik, O., Goloub, P., Tanr, D., Lopatsin, A.,

Denisov, S., Lapyonok, T., and Karol, Y.: The retrieval of aerosol microphysical properties in the vertical column using combined

www.atmos-meas-tech.net/9/1181/2016/ Atmos. Meas. Tech., 9, 11811205, 2016

1202 A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC)

lidar/photometer data: a step to integrating photometer and lidar networks, Proceedings of the 25th International Laser Radar Conference, 59 July 2010, St.-Petersburg, Russia, 10871091. 2010b.

Chaikovsky, A., Dubovik, O., Goloub, P., Tanr, D., Pappalardo,G., Wandinger, U., Chaikovskaya, L., Denisov, S., Grudo, Y., Lopatsin, A., Karol, Y., Lapyonok, T., Korol, M., Osipenko, F., Savitski, D., Slesar, A., Apituley, A., Arboledas, L. A., Bini-etoglou, I., Kokkalis, P., Granados-Muoz, M. J., Papayannis,A., Perrone, M. R., Pietruczuk, A., Pisani, G., Rocadenbosch, F., Sicard, M., De Tomasi, F., Wagner, J., and Wang, X.: Algorithm and software for the retrieval of vertical aerosol properties using combined lidar/ radiometer data: Dissemination in EARLINET, in: Proceedings of the 26th International Laser and Radar Conference, 2529 June 2012, Porto Heli, Greece, 399402, 2012.Chaikovsky, A. P., Dubovik, O., Holben, B. N., and BrilA. I.: Methodology to retrieve atmospheric aerosol parameters by combining ground-based measurements of multi-wavelength lidar and sun sky-scanning radiometer, Eighth International Symposium on Atmospheric and Ocean Optics: Atmospheric Physics, 28 February 2002, Proc. SPIE 4678, 257268, doi:http://dx.doi.org/10.1117/12.458450

Web End =10.1117/12.458450 http://dx.doi.org/10.1117/12.458450

Web End = , 2002.

Chaikovsky, A. P., Ivanov, A. P., Balin, Y. S., Elnikov, A. V., Tulinov, G. F., Plusnin, I. I., Bukin, O. A., and Chen B. B.: CISLiNet lidar network for monitoring aerosol and ozone: methodology and instrumentation, Atmos. and Oceanic Opt., 18, 958963, 2005.

Chaykovskii, A. P.: Method for investigating the structure of the stratospheric aerosol layer based on laser echo depolarization measurements, Atmos. and Oceanic Opt., 3, 12211223, 1990.Colarco, P. R., Schoeberl, M. R., Doddridge B. G., Marufu, L. T.,

Torres O., and Welton, E. J.: Transport of smoke from Canadian forest res to the surface near Washington, D.C.: Injection height, entrainment, and optical properties, J. Geophys. Res., 109, D06203, doi:http://dx.doi.org/10.1029/2003JD004248

Web End =10.1029/2003JD004248 http://dx.doi.org/10.1029/2003JD004248

Web End = , 2004.

Committee on Space Research: NASA National Space Science Data Center: COSPAR International Reference Atmosphere (CIRA-86): Global Climatology of Atmospheric Parameters. NCAS British Atmospheric Data Centre, available at: http://browse.ceda.ac.uk/browse/badc/cira/data/

Web End =http://browse. http://browse.ceda.ac.uk/browse/badc/cira/data/

Web End =ceda.ac.uk/browse/badc/cira/data/ (last access: 13 March 2016), 2006.

Cuesta, J., Flamant, H. P., and Flamant, C.: Synergetic technique combining elastic backscatter lidar data and sunphotometer AERONET inversion for retrieval by layer of aerosol optical and microphysical properties, Appl. Optics, 47, 45984611, 2008.

David, G., Thomas, B., Nousiainen, T., Miffre, A., and Rairoux,P.: Retrieving simulated volcanic, desert dust and sea-salt particle properties from two/three-component particle mixtures using UV-VIS polarization lidar and T matrix, Atmos. Chem. Phys., 13, 67576776, doi:http://dx.doi.org/10.5194/acp-13-6757-2013

Web End =10.5194/acp-13-6757-2013 http://dx.doi.org/10.5194/acp-13-6757-2013

Web End = , 2013.

Doicu, A., Trautmann T., and Schreier, F. Numerical Regularization for Atmospheric Inverse Problem, Springer-Verlag, Berlin-Heidelberg, Germany, 426 pp., 2010.

Dubovik, O. and King, M.: A exible inversion algorithm for retrieval of aerosol optical properties from Sun and sky radiance measurements, J. Geophys. Res., 105, 2067320696, 2000.Dubovik, O., Smirnov, A., Holben, B. N., King, M. D., Kaufman, Y., Eck T. F., and Slutsker, I.: Accuracy assessments of

aerosol optical properties retrieved from Aerosol Robotic Network (AERONET) sun and sky radiance measurements, J. Geophys. Res., 105, 97919800, 2000.

Dubovik, O., Holben, B. N., Lapyonok, T., Sinyuk, A., Mishchenko,M. I., Yang, P., and Slutsker, I.: Non-spherical aerosol retrieval method employing light scattering by spheroids, Geophys. Res. Lett., 29, 54-154-4, doi:http://dx.doi.org/10.1029/2001GL014506

Web End =10.1029/2001GL014506 http://dx.doi.org/10.1029/2001GL014506

Web End = , 2002. Dubovik, O.: Optimization of Numerical Inversion in Photopolari-metric Remote Sensing, in: Photopolarimetry in Remote Sensing, edited by: Videen, G., Yatskiv, Y., and Mishchenko, M., Kluwer Academic Publishers, Dordrecht, the Netherlands, 65 106, 2004.

Dubovik, O., Sinyuk, A., Lapyonok, T., Holben, B. N., Mishchenko,M., Yang, P., Eck, T. F., Volten, H., Munoz, O., Veihelmann, B., van der Zande, W. J., Leon, J.-F., Sorokin, M., and Slutsker, I.: Application of spheroid models to account for aerosol particle nonsphericity in remote sensing of desert dust, J. Geophys. Res., 111, D11208, doi:http://dx.doi.org/10.1029/2005JD006619

Web End =10.1029/2005JD006619 http://dx.doi.org/10.1029/2005JD006619

Web End = , 2006.

Dubovik, O., Herman, M., Holdak, A., Lapyonok, T., Tanr, D., Deuz, J. L., Ducos, F., Sinyuk, A., and Lopatin, A.: Statistically optimized inversion algorithm for enhanced retrieval of aerosol properties from spectral multi-angle polarimetric satellite observations, Atmos. Meas. Tech., 4, 9751018, doi:http://dx.doi.org/10.5194/amt-4-975-2011

Web End =10.5194/amt-4- http://dx.doi.org/10.5194/amt-4-975-2011

Web End =975-2011 , 2011.

Dubovik, O., Lapyonok, T., Litvinov, P., Herman, M., Fuertes, D., Ducos, F., Lopatin, A., Chaikovsky, A., Torres, B., Derimian, Y., Xin Huang, Aspetsberger, M., and Federspiel, C.: GRASP: a versatile algorithm for characterizing the atmosphere, SPIE News-room, doi:http://dx.doi.org/10.1117/2.1201408.005558

Web End =10.1117/2.1201408.005558 http://dx.doi.org/10.1117/2.1201408.005558

Web End = , 2014.

Fleming, E., Chandra, S., Shoeberl, M., and Barnett, J.: Monthly Mean Global Climatology of Temperature, Wind, Geopotential Height, and Pressure for 0120 km, NASA Technical Memorandum 100697, 91 pp., 1988.

Freudenthaler, V., Gross, S., Engelmann, R., Mattis, I, Wandinger,U., Pappalardo, G., Amodeo, A., Giunta, A., DAmico, G., Chaikovsky A., Osipenko, F., Slesar, A., Nicolae, D., Belegante,L., Talianu, C., Serikov, I., Linne, H., Jansen, F., Wilson, K., Graaf, M., Apituley, A., Trickl, T., Giehl, H., and Adam, M.: EARLI09 Direct intercomparison of eleven EARLINET lidar systems, Proceedings of the 25th International Laser Radar Conference, 59 July 2010, St.-Petersburg, Russia, 891894, 2010. Ganguly, D., Ginoux, P., Ramaswamy, V., Dubovik, O., Welton,J., Reid, E. A., and Holben, B. N.: Inferring the composition and concentration of aerosols by combining AERONET and MPLNET data: Comparison with other measurements and utilization to evaluate GCM output, J. Geophys. Res., 114, D16203, doi:http://dx.doi.org/10.1029/2009JD011895

Web End =10.1029/2009JD011895 http://dx.doi.org/10.1029/2009JD011895

Web End = , 2009a.

Ganguly, D., Ginoux, P., Ramaswamy, V., Ganguly, D., Winker,D. M., Holben, B. N., and Tripathi, S. N.: Retrieving the composition and concentration of aerosols over the Indo-Gangetic basin using CALIOP and AERONET data, Geophys. Res. Lett., 36, L13806, doi:http://dx.doi.org/10.1029/2009GL038315

Web End =10.1029/2009GL038315 http://dx.doi.org/10.1029/2009GL038315

Web End = , 2009b.

Gasteiger, J., Gro, S., Freudenthaler, V., and Wiegner, M.: Volcanic ash from Iceland over Munich: mass concentration retrieved from ground-based remote sensing measurements, Atmos. Chem. Phys., 11, 22092223, doi:http://dx.doi.org/10.5194/acp-11-2209-2011

Web End =10.5194/acp-11-2209- http://dx.doi.org/10.5194/acp-11-2209-2011

Web End =2011 , 2011.

Granados-Muoz, M. J., Guerrero-Rascado, J. L., Bravo-Aranda, J.A., Navas-Guzmn, F., Valenzuela, A., Lyamani, H., Chaikovsky,

Atmos. Meas. Tech., 9, 11811205, 2016 www.atmos-meas-tech.net/9/1181/2016/

A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC) 1203

A., Wandinger, U., Ansmann, A., Dubovik, O., Grudo, J., and Alados-Arboledas, L.: Retrieving aerosol microphysical properties by Lidar-Radiometer Inversion Code (LIRIC) for different aerosol types, J. Geophys. Res., 119, 48364858, doi:http://dx.doi.org/10.1002/2013JD021116

Web End =10.1002/2013JD021116 http://dx.doi.org/10.1002/2013JD021116

Web End = , 2014.

Granados-Muoz, M. J., Bravo-Aranda, J. A., Baumgardner, D., Guerrero-Rascado, J. L., Prez-Ramrez, D., Navas-Guzmn, F., Veselovskii, I., Lyamani, H., Valenzuela, A., Olmo, F. J., Titos,G., Andrey, J., Chaikovsky, A., Dubovik, O., Gil-Ojeda, M., and Alados-Arboledas, L.: Study of aerosol microphysical properties proles retrieved from ground-based remote sensing and aircraft in-situ measurements during a Saharan dust event, Atmos. Meas.Tech. Discuss., 8, 92899338, doi:http://dx.doi.org/10.5194/amtd-8-9289-2015

Web End =10.5194/amtd-8-9289-2015 http://dx.doi.org/10.5194/amtd-8-9289-2015

Web End = , 2015a.

Granados-Muoz, M. J., Navas-Guzmn, F., Guerrero-Rascado, J.L., Bravo-Aranda, J. A., Binietoglou, I., Pereira, S. N., Basart,S., Baldasano, J. M., Belegante, L., Chaikovsky, A., Comern,A., DAmico, G., Dubovik, O., Ilic, L., Kokkalis, P., Muoz-Porcar, C., Nickovic, S., Nicolae, D., Olmo, F. J., Papayannis,A., Pappalardo, G., Rodrguez, A., Schepanski, K., Sicard, M., Vukovic, A., Wandinger, U., Dulac, F., and Alados-Arboledas,L.: Proling of aerosol microphysical properties at several EARLINET/AERONET sites during July 2012 ChArMEx/EMEP campaign, Atmos. Chem. Phys. Discuss., 15, 3283132887, doi:10.5194/acpd-15-32831-2015, 2015b.

Hansen, P. C.: The L-curve and its use in the numerical treatment of inverse problems. Computational Inverse Problems in Electro-cardiology, WIT Press, Southampton, UK, 119142, 2001.Hashimoto, M., Nakajima, T., Dubovik, O., Campanelli, M., Che,H., Khatri, P., Takamura, T., and Pandithurai, G.: Development of a new data-processing method for SKYNET sky radiometer observations, Atmos. Meas. Tech., 5, 27232737, doi:http://dx.doi.org/10.5194/amt-5-2723-2012

Web End =10.5194/amt- http://dx.doi.org/10.5194/amt-5-2723-2012

Web End =5-2723-2012 , 2012.

Hoff, R., Moshary, F., Ahmed, S., Gross, B., Mc-Cormick, M. P., and Parsiani, H.: Plan for the implementation of the CREST Lidar Network (CLN). CREST Publication Series vol 07, no. 01-2009, University of Maryland, Baltimore County/City College of New York/Hampton University/Mayaguez University, 2009.Holben, B. N., Eck, T. F., Slutsker, I., Tanr, D., Buis, J. P., Setzer, A., Vermote, E., Reagan, J. A., Kaufman, Y. J., Nakajima,T., Lavenu, F., I, J., and Smirnov, A.: AERONET A federated instrument network and data archive for aerosol characterization, Remote Sens. Environ., 66, 116, 1998.

Holben, B. N., Eck, T. F., Slutsker, I., Smirnov, A., Sinyuk, A., Schafer, J., Giles, D., and Dubovik, O.: AERONETs Version2.0 quality assurance criteria, in Remote Sensing of the Atmosphere and Clouds, edited by: Tsay, S.-C., Nakajima, T., Singh, R. P., and Sridharan R., Proc. SPIE, 6408, 64080Q, doi:http://dx.doi.org/10.1117/12.706524

Web End =10.1117/12.706524 http://dx.doi.org/10.1117/12.706524

Web End = , 2006.

Kabashnikov, V., Chaikovsky, A., Denisov, S., Dubovik, O., Goloub, P., Ivanov, A., V. Kusmin, V., Kazeruk, D., Korol,M., Karol, Y., Lopatsin, A., Miatselskaya, N., Osipenko, F., Pietruczuk, F., Slesar, A., Sobolewski, H., and Tanre, D.: Long range transport of air pollution in the east European regions: four years observations, in Proceedings of the 25th International Laser Radar Conference, 59 July 2010, St.-Petersburg, Russia, 1043 1046, 2010.

Klett, D.: Stable analytical inversion solution for processing lidar returns, Appl. Optics, 20, 211220, 1981.

Kokkalis, P., Papayannis, A., Amiridis, V., Mamouri, R. E., Chaikovsky, A., Dubovik, O., and Tsekeri, A.: Evaluation of ne mode lidar concentration retrievals using airborne in-situ measurements, Reviewed and revised papers of the 26th International Lidar Radar Conference, 2529 June 2012, Porto Heli, Peloponnesus, Greece, 617620, 2012.

Kokkalis, P., Papayannis, A., Amiridis, V., Mamouri, R. E., Veselovskii, I., Kolgotin, A., Tsaknakis, G., Kristiansen, N.I., Stohl, A., and Mona, L.: Optical, microphysical, mass and geometrical properties of aged volcanic particles observed over Athens, Greece, during the Eyjafjallajkull eruption in April 2010 through synergy of Raman lidar and sunphotometer measurements, Atmos. Chem. Phys., 13, 93039320, doi:http://dx.doi.org/10.5194/acp-13-9303-2013

Web End =10.5194/acp-13-9303-2013 http://dx.doi.org/10.5194/acp-13-9303-2013

Web End = , 2013.

Levenberg, K.: A method for the solution of certain problems in least squares, Q. Appl. Math., 2, 164168, 1944.

Levy, R. C, Remer, L. A., Mattoo, S., Vermote, E. F., and Kaufman, Y. J.: Second-generation operational algorithm: Retrieval of aerosol properties over land from inversion of moderate resolution imaging spectroradiometer spectral reectance, J. Geophys. Res, 112, D13211, doi:http://dx.doi.org/10.1029/2006JD007811

Web End =10.1029/2006JD007811 http://dx.doi.org/10.1029/2006JD007811

Web End = , 2007a.

Levy, R. C, Remer, L. A., and Dubovik, O.: Global aerosol optical models and application to Moderate Resolution Imaging Spec-troradiometer aerosol retrieval over land, J. Geophys. Res., 112, D13210, doi:http://dx.doi.org/10.1029/2006JD007815

Web End =10.1029/2006JD007815 http://dx.doi.org/10.1029/2006JD007815

Web End = , 2007b.

Li, Z., Goloub, P., Dubovik, O., Blarel, L., Zhang, W., Podvin, T., Sinyuk A., Sorokin, M., Chen, H., Holben, B. N., Tanre, D., Canini, M., and Buis, J.-P: Improvements for ground-based remote sensing of atmospheric aerosol properties by additional polarimetric measurements, J. Quant. Spectrosc. Ra., 110, 1954 1961, 2009.

Lopatin, A., Dubovik, O., Chaikovsky, A., Goloub, P., Lapyonok,T., Tanr, D., and Litvinov, P.: Enhancement of aerosol characterization using synergy of lidar and sun-photometer coincident observations: the GARRLiC algorithm, Atmos. Meas. Tech., 6, 20652088, doi:http://dx.doi.org/10.5194/amt-6-2065-2013

Web End =10.5194/amt-6-2065-2013 http://dx.doi.org/10.5194/amt-6-2065-2013

Web End = , 2013.

Lund Myhre, C., Toledano, C., Myhre, G., Stebel, K., Yttri, K. E., Aaltonen, V., Johnsrud, M., Frioud, M., Cachorro, V., de Frutos,A., Lihavainen, H., Campbell, J. R., Chaikovsky, A. P., Shiobara,M., Welton, E. J., and Trseth, K.: Regional aerosol optical properties and radiative impact of the extreme smoke event in the European Arctic in spring 2006, Atmos. Chem. Phys., 7, 5899 5915, doi:http://dx.doi.org/10.5194/acp-7-5899-2007

Web End =10.5194/acp-7-5899-2007 http://dx.doi.org/10.5194/acp-7-5899-2007

Web End = , 2007.

Marquardt, D.: An algorithm for least-squares estimation of nonlinear parameters, SIAM J. Appl. Math., 11, 431441, 1963. Matthias, V., Bsenberg, J., Freudenthaler, V., Amodeo, A., Balis,D., Chaikovsky, A., Chourdakis, G., Comeron, A., Delaval, A., De Tomasi, F., Eixmann, R., Hgrd, A., Komguem, L., Kreipl,S., Matthey, R., Mattis, I., Rizi, V., Rodriguez, J. A., Simeonov,V., and Wang, X.: Aerosol lidar intercomparison in the framework of the EARLINET project. 1. Instruments, Appl. Optics, 43, 961976, 2004.

McKendry, I., Strawbridge, K., ONeill, N., Macdonald, A. M., Liu, P., Leaitch, R., Anlauf, K., Jaegle, L., Fairlie, T. D., and Westphal, D.: Trans-Pacic transport of Saharan dust to western North America: A case study, J. Geophys. Res., 112, D01103, doi:http://dx.doi.org/10.1029/2006JD007129

Web End =10.1029/2006JD007129 http://dx.doi.org/10.1029/2006JD007129

Web End = , 2007.

McKendry, I., Strawbridge, K., Karumudi, M. L., ONeill, N., Macdonald, A. M., Leaitch, R., Jaffe, D., Cottle, P., Sharma, S., Sheri-

www.atmos-meas-tech.net/9/1181/2016/ Atmos. Meas. Tech., 9, 11811205, 2016

1204 A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC)

dan, P., and Ogren, J.: Californian forest re plumes over Southwestern British Columbia: lidar, sunphotometry, and mountain-top chemistry observations, Atmos. Chem. Phys., 11, 465477, doi:http://dx.doi.org/10.5194/acp-11-465-2011

Web End =10.5194/acp-11-465-2011 http://dx.doi.org/10.5194/acp-11-465-2011

Web End = , 2011.

Mishchenko, M. I., Hovenier, J. W., and Travis, L. D.: Light Scattering by Nonspherical Particles, Elsevier, New York, USA, 2000.Mishchenko, M. I., Travis, L. D., and Lacis, A. A.: Scattering, Absorption, and Emission of Light by Small Particles, Cambridge Univ. Press, New York, USA, 2002.

Mller, D., Mattis, I., Wandinger, U., Ansmann, A., Althausen, D., Dubovik, O., Eckhardt, S., and Stohl, A.: Saharan dust over a Central European EARLINET-AERONET site: Combined observations with Raman lidar and Sun photometer, J. Geophys.Res., 108, 4345, doi:http://dx.doi.org/10.1029/2002JD002918

Web End =10.1029/2002JD002918 http://dx.doi.org/10.1029/2002JD002918

Web End = , 2003.

Mller, D., Veselovskii, I., Kolgotin, A., Tesche, M., Ansmann,A., and Dubovik, O.: Vertical proles of pure dust and mixed smokedust plumes inferred from inversion of multiwavelength Raman/polarization lidar data and comparison to AERONET retrievals and in situ observations, Appl. Optics, 52, 31783202 doi:http://dx.doi.org/10.1364/AO.52.003178

Web End =10.1364/AO.52.003178 http://dx.doi.org/10.1364/AO.52.003178

Web End = , 2013.

Murayama, T., Sugimoto, N., Matsui, I., Liu, Z., Sakai, T., Shibata,T., Iwsasaka, Y., Won, J.-G., Yoon, S.-C., Li, T., Zhou, J., and Hu, H.: Lidar Network Observation of Asian Dust, in: Advances in Laser Remote sensing, Selected papers 20th Int. Laser Radar Conference (ILRC), 1014 July 2000, Vichy, France, edited by: Dabas, A., Loth, C., and Pelon, J., 169177, 2001.

Nemuc, A., Vasilescu, J., Talianu, C., Belegante, L., and Nicolae,D.: Assessment of aerosols mass concentrations from measured linear particle depolarization ratio (vertically resolved) and simulations, Atmos. Meas. Tech., 6, 32433255, doi:http://dx.doi.org/10.5194/amt-6-3243-2013

Web End =10.5194/amt-6- http://dx.doi.org/10.5194/amt-6-3243-2013

Web End =3243-2013 , 2013.

Papayannis, A., Amiridis, V., Mona, L., Tsaknakis, G., Balis, D., Bosenberg, J., Chaikovski, A., De Tomasi, F., Grigorov, I., Mattis, I., Mitev, V., Muller, D., Nickovic, S., Perez, C., Pietruczuk,A., Pisani, G., Ravetta, F., Rizi, V., Sicard, M., Trickl, T., Wiegner, M., Gerding, M., Mamouri, R. E., DAmico, G., and Pappalardo, G.: Systematic lidar observations of Saharan dust over Europe in the frame of EARLINET (20002002), J. Geophys.Res., 113, D10204, doi:http://dx.doi.org/10.1029/2007JD009028

Web End =10.1029/2007JD009028 http://dx.doi.org/10.1029/2007JD009028

Web End = , 2008.Papayannis, A., Mamouri, R. E., Amiridis, V., Giannakaki, E.,

Veselovskii, I., Kokkalis, P., Tsaknakis, G., Balis, D., Kristiansen, N. I., Stohl, A., Korenskiy, M., Allakhverdiev, K., Huseyinoglu, M. F., and Baykara, T.: Optical properties and vertical extension of aged ash layers over the Eastern Mediterranean as observed by Raman lidars during the Eyjafjallajokull eruption in May 2010, Atmos. Environ., 48, 5665, 2012.

Papayannis, A., Nicolae, D., Kokkalis, P., Binetoglou, I., Talianu,C., Belegante, L., Tsaknakis., G., Cazacu, M. M., Vetres, I., and Ilic, I.: Optical, size and mass properties of mixed type aerosols in Greece and Romania as observed by synergy of lidar and sunphotometers in combination with model simulations: A case study, Sci. Tot. Environ., 500501, 277294, 2014.Pappalardo, G., Amodeo, A., Apituley, A., Comeron, A., Freudenthaler, V., Linn, H., Ansmann, A., Bsenberg, J., DAmico,G., Mattis, I., Mona, L., Wandinger, U., Amiridis, V., Alados-Arboledas, L., Nicolae, D., and Wiegner, M.: EARLINET: towards an advanced sustainable European aerosol lidar network, Atmos. Meas. Tech., 7, 23892409, doi:http://dx.doi.org/10.5194/amt-7-2389-2014

Web End =10.5194/amt-7-2389- http://dx.doi.org/10.5194/amt-7-2389-2014

Web End =2014 , 2014.

Pietruczuk, A. and Chaikovsky, A. P.: Properties of re smoke in eastern Europe measured by remote sensing methods, Proceedings of SPIE, 6745, 67451T, doi:http://dx.doi.org/10.1117/12.740916

Web End =10.1117/12.740916 http://dx.doi.org/10.1117/12.740916

Web End = , 2007.

Rodgers C. D.: Inverse methods for atmospheric sounding. Theory and Practice. Word Scientic, Singapore, 238 pp., 2000.

Russell, P. B., Swissler, T. J., and McCormick, M. P.: Methodology for error analysis and simulation of lidar aerosol measurements, Appl. Optics, 18, 37833797, 1979.

Schuster, G. L., Vaughan, M., MacDonnell, D., Su, W., Winker,D., Dubovik, O., Lapyonok, T., and Trepte, C.: Comparison of CALIPSO aerosol optical depth retrievals to AERONET measurements, and a climatology for the lidar ratio of dust, Atmos. Chem. Phys., 12, 74317452, doi:http://dx.doi.org/10.5194/acp-12-7431-2012

Web End =10.5194/acp-12-7431-2012 http://dx.doi.org/10.5194/acp-12-7431-2012

Web End = , 2012.

Slesar, A. S., Chaykovskii, A. P., Ivanov, A. P., Denisov, S. V., Korol, M. M., Osipenko, F. P., Balin, Yu. S., Kokhanenko, G. P., and Penner I. E.: Fotodetector modules for CIS-LiNet lidar network stations, Atmos. and Oceanic Opt., 26, 10731081, 2013 (in Russian).

Slesar, A. S., Chaikovskii, A. P., Denisov, S. V., Korol, M. M., Osipenko, F. P., Balin, Yu. S., Kokhanenko, G. P., Penner, I.E., and Novoselov M. M.: Development of photodetectors for recording lidar signals in the photon counting and analog modes, Proc. SPIE 9680, 21st International Symposium Atmospheric and Ocean Optics: Atmospheric Physics, 19 November 2015, 96802Q, doi:http://dx.doi.org/10.1117/12.2203666

Web End =10.1117/12.2203666 http://dx.doi.org/10.1117/12.2203666

Web End = , 2015.

Takamura, T., Nakajima, T., and SKYNET community group: Overview of SKYNET and its activities, Optica Pura y Aplicada, 37, 33033308, 2004.

Tarantola, A.: Inverse Problem Theory: Methods for Data Fitting and Model Parameter Estimation, Elsevier, Amsterdam, the Netherlands, 614 pp., 1987.

Tesche, M., Ansmann, A., Mller, D., Althausen, D., Engelmann,R., Freudenthaler, V., and Gro, S.: Separation of dust and smoke proles over Cape Verde by using multiwavelength Raman and polarization lidars during SAMUM 2008, J. Geophys. Res., 114, D13202, doi:http://dx.doi.org/10.1029/2009JD011862

Web End =10.1029/2009JD011862 http://dx.doi.org/10.1029/2009JD011862

Web End = , 2009.

Tikhonov, A. N. and Arsenin, V. Y.: Solution of Ill-Posed Problems,Wiley, New York, USA, 300 pp., 1977.

Trefethen, L. N. and Bau, D.: Numerical Linear Algebra, SIAM,Philadelphia, USA, 361 pp., 1997.

Tsekeri, A., Amiridis, V., Kokkalis, P., Basart, S., Chaikovsky, A., Dubovik, O., Mamouri, R. E., Baldasano, J. M., and Papayannis,A.: Evaluation of dust modelling using a synergetic algorithm of lidar and sunphotometer data. Reviewed and Revised papers of the 26th International Lidar Radar Conference, 2529 June 2012, Porto Heli, Peloponnesus, Greece, 621624, 2012.

Tsekeri, A., Amiridis, V., Kokkalis, P., Basart, S., Chaikovsky, A., Dubovik, O., Mamouri, R. E., Papayannis, A., and Baldasano,J. M.: Application of a Synergetic Lidar and Sunphotometer Algorithm for the Characterization of a Dust Event Over Athens, Greece, British J. of Environment and Climate Change, 3, 531 546, doi:http://dx.doi.org/10.9734/BJECC/2013/2615

Web End =10.9734/BJECC/2013/2615 http://dx.doi.org/10.9734/BJECC/2013/2615

Web End = , 2013.

Turchin, V. F., Kozlov, V. P., and Malkevich, M. S.: The use of mathematical-statistics methods in the solution of incorrectly posed problems, Sov. Phys. Uspekhi, 13, 681703, 1971. Twomey, S.: Introduction to the mathematics of inversion in remote sensing and indirect measurements, Elsevier, Amsterdam, the Netherlands, 1977.

Atmos. Meas. Tech., 9, 11811205, 2016 www.atmos-meas-tech.net/9/1181/2016/

A. Chaikovsky et al.: Lidar-Radiometer Inversion Code (LIRIC) 1205

Veselovskii, I, Kolgotin, A., Mller, D., and Whiteman, D.: Information content of multiwavelength lidar data with respect to microphysical particle properties derived from eigenvalue analysis, Appl. Optics, 44, 52925303, 2005.

Veselovskii, I., Dubovik, O., Kolgotin, A., Lapyonok ,T., Di Girolamo, P., Summa, D., Whiteman, D. N., Mishchenko, M., and Tanre D.: Application of randomly oriented spheroids for retrieval of dust particle parameters from multiwave-length lidar measurements. J. Geophys. Res., 115, D21203, doi:http://dx.doi.org/10.1029/2010JD014139

Web End =10.1029/2010JD014139 http://dx.doi.org/10.1029/2010JD014139

Web End = , 2010.

Vogel, C. R.: Computational Methods for Inverse Problems, SIAM,Philadelphia, PA, USA, 183 pp., 2002.

Volten, H., Muoz, O., Rol, E., de Haan, J. F., Vassen, W., Hovenier,J. W., Muinonen, K., and Nousiainen, T.: Scattering matrices of mineral particles at 441.6 nm and 632.8 nm., J. Geophys. Res., 106, 1737517401, 2001.

Wagner, J., Ansmann, A., Wandinger, U., Seifert, P., Schwarz, A.,

Tesche, M., Chaikovsky, A., and Dubovik, O.: Evaluation of the Lidar/Radiometer Inversion Code (LIRIC) to determine micro-physical properties of volcanic and desert dust, Atmos. Meas.Tech., 6, 17071724, doi:http://dx.doi.org/10.5194/amt-6-1707-2013

Web End =10.5194/amt-6-1707-2013 http://dx.doi.org/10.5194/amt-6-1707-2013

Web End = , 2013.Wandinger, U., Freudenthaler, V., Baars, H., Amodeo, A., Engelmann, R., Mattis, I., Gro, S., Pappalardo, G., Giunta, A., DAmico, G., Chaikovsky, A., Osipenko, F., Slesar, A., Nicolae,D., Belegante, L., Talianu, C., Serikov, I., Linn, H., Jansen, F., Apituley, A., Wilson, K. M., de Graaf, M., Trickl, T., Giehl, H., Adam, M., Comern, A., Muoz, C., Rocadenbosch, F., Sicard,M., Toms, S., Lange, D., Kumar, D., Pujadas, M., Molero, F., Fernndez, A. J., Alados-Arboledas, L., Bravo-Aranda, J. A., Navas-Guzmn, F., Guerrero-Rascado, J. L., Granados-Muoz,M. J., Preiler, J., Wagner, F., Gausa, M., Grigorov, I., Stoyanov, D., Iarlori, M., Rizi, V., Spinelli, N., Boselli, A., Wang, X., Lo Feudo, T., Perrone, M. R., De Tomasi, F., and Burlizzi, P.: EARLINET instrument intercomparison campaigns: overview on strategy and results, Atmos. Meas. Tech. Discuss., 8, 10473 10522, doi:http://dx.doi.org/10.5194/amtd-8-10473-2015

Web End =10.5194/amtd-8-10473-2015 http://dx.doi.org/10.5194/amtd-8-10473-2015

Web End = , 2015.

Welton, E. J., Campbell, J. R., Berkoff, T. A., Spinhirne, J. D., Tsay,S.-C., Holben, B., and Shiobara, M.: The Micro-pulse Lidar Network (MPL-Net), in: Lidar Remote Sensing in Atmospheric and Earth Sciences, Reviewed and revised papers at the twenty-rst International Laser Radar Conference (ILRC21), 1705 Quebec, Canada, 812 July 2002, edited by: Bisonnette, L. R., Roy, G., and Valle, G., Defence R&D CanadaVacartier, 285288, 2002. Winker, D. M., Hostetler, C. A., Vaughan, M. A., and Omar,A. H.: CALIOP Algorithm Theoretical Basis Document. Part 1: CALIOP Instrument, and Algorithms Overview, PC-SCI-202, Part 1, Release 2.0, 9 September 2006, 29 pp., 2006.

www.atmos-meas-tech.net/9/1181/2016/ Atmos. Meas. Tech., 9, 11811205, 2016

Word count: 16781

Show less

Abstract

Translate

This paper presents a detailed description of LIRIC (LIdar-Radiometer Inversion Code) algorithm for simultaneous processing of coincident lidar and radiometric (sun photometric) observations for the retrieval of the aerosol concentration vertical profiles. As the lidar/radiometric input data we use measurements from European Aerosol Research Lidar Network (EARLINET) lidars and collocated sun-photometers of Aerosol Robotic Network (AERONET). The LIRIC data processing provides sequential inversion of the combined lidar and radiometric data. The algorithm starts with the estimations of column-integrated aerosol parameters from radiometric measurements followed by the retrieval of height dependent concentrations of fine and coarse aerosols from lidar signals using integrated column characteristics of aerosol layer as a priori constraints. The use of polarized lidar observations allows us to discriminate between spherical and non-spherical particles of the coarse aerosol mode.</br></br>The LIRIC software package was implemented and tested at a number of EARLINET stations. Intercomparison of the LIRIC-based aerosol retrievals was performed for the observations by seven EARLINET lidars in Leipzig, Germany on 25 May 2009. We found close agreement between the aerosol parameters derived from different lidars that supports high robustness of the LIRIC algorithm. The sensitivity of the retrieval results to the possible reduction of the available observation data is also discussed.

Details

Title

Lidar-Radiometer Inversion Code (LIRIC) for the retrieval of vertical aerosol properties from combined lidar/radiometer data: development and distribution in EARLINET

Pages

1181-1205

Publication year

2016

Publication date

2016

Publisher

Copernicus GmbH

ISSN

18671381

e-ISSN

18678548

Source type

Scholarly Journal

Language of publication

English

DOI

https://doi.org/10.5194/amt-9-1181-2016

ProQuest document ID

1774521381

Lidar-Radiometer Inversion Code (LIRIC) for the retrieval of vertical aerosol properties from combined lidar/radiometer data: development and distribution in EARLINET

Jump to:

Full Text

Abstract

Details

Suggested sources