Simultaneous estimation of aerosol cloud concentration and spectral backscatter from multiple-wavelength lidar data
Applied Optics, Vol. 47, Issue 24, pp. 4309-4320 (2008)
http://dx.doi.org/10.1364/AO.47.004309
Enhanced HTML Acrobat PDF (1231 KB)
Abstract
We present a sequential algorithm for estimating both concentration dependence on range and time and backscatter coefficient spectral dependence of optically thin localized atmospheric aerosols using data from rapidly tuned lidar. The range dependence of the aerosol is modeled as an expansion of the concentration in an orthonormal basis set whose coefficients carry the time dependence. Two estimators are run in parallel: a Kalman filter for the concentration range and time dependence and a maximum-likelihood estimator for the aerosol backscatter wavelength and time dependence. These two estimators exchange information continuously over the data-processing stream. The state model parameters of the Kalman filter are also estimated sequentially together with the concentration and backscatter. Lidar data collected prior to the aerosol release are used to estimate the ambient lidar return. The approach is illustrated on atmospheric backscatter long-wave infrared (
© 2008 Optical Society of America
OCIS Codes
(010.1100) Atmospheric and oceanic optics : Aerosol detection
(010.1120) Atmospheric and oceanic optics : Air pollution monitoring
(010.3640) Atmospheric and oceanic optics : Lidar
(280.1415) Remote sensing and sensors : Biological sensing and sensors
(010.1350) Atmospheric and oceanic optics : Backscattering
ToC Category:
Atmospheric and Oceanic Optics
History
Original Manuscript: May 2, 2008
Revised Manuscript: July 2, 2008
Manuscript Accepted: July 9, 2008
Published: August 12, 2008
Citation
Russell E. Warren, Richard G. Vanderbeek, Avishai Ben-David, and Jeffrey L. Ahl, "Simultaneous estimation of aerosol cloud concentration and spectral backscatter from multiple-wavelength lidar data," Appl. Opt. 47, 4309-4320 (2008)
http://www.opticsinfobase.org/ao/abstract.cfm?URI=ao-47-24-4309
Sort: Year | Journal | Reset
References
- The text Elastic Lidar, by V. A. Kovalev and W. E. Eichinger (Wiley, 2004) has an extensive treatment and bibliography of lidar analysis methods for aerosols distributed continuously over extended paths. [CrossRef]
- We are assuming a single aerosol material here. Generalization to more than one material present at a given time and range is currently being investigated.
- A note on units. From Eqs. we see that the physical units of αC must be those of β, 1/(m sr), for example. Since concentration is typically expressed in units of mg/m^{3}, this requires the units of α to be m^{2}/(mg sr). Setting g=1 of course means that these units no longer apply and, consequently, only relative backscatter and concentration can be estimated by this choice. This is described in more detail in Section .
- B. R. Frieden, Probability, Statistical Optics, and Data Testing (Springer-Verlag, 1983).
- A Gaussian model for Sg is not the only choice, of course. Using a step function with width r_{g} would give S_{g}(f)=σ_{g}2r_{g}sinc(fr_{g}). This model has many more ripples than the Gaussian and provides less noise damping for a given rg.
- Even with the choice of maximum likelihood, there is the decision of whether the estimates should be computed independently for each time step or be processed by a recursive algorithm, such as recursive least squares. We have adopted the former choice here since the estimates are to be used to train and test classifiers for which independent training samples are required. Actual deployment may be better served by a recursive algorithm.
- A related concern is the huge amount of background data needed to construct suitable sample estimates of the range-covariance matrices. This is not compatible with the limited amount of background data typically available.
- See S. Mallat, Wavelet Tour of Signal Processing, Section 9.1.3 and Appendix A6, for an excellent discussion of the KL expansion (Academic, 1998).
- If our covariance function, given by Eq. , had circular periodicity, Mallat shows that the Fourier basis would be the KL basis. We cannot claim this, but for target aerosols not close to either end of the sampling range, there is little error in making this assumption.
- See, for example, F. L. Lewis, Optimal Estimation (Wiley, 1986).
- A. P. Dempster, N. M. Laird, and D. B. Rubin, “Maximum likelihood from incomplete data via the EM algorithm (with discussion),” J. R. Stat. Soc. Ser. B. Methodol. 39, 1-38 (1977).
- R. H. Shumway and D. S. Stoffer, Time Series Analysis and Its Applications (Springer-Verlag, 2000), Chap. 4.
- S. M. Kay, Fundamentals of Statistical Signal Processing: Estimation Theory (Prentice Hall, 1993) has very readable treatments of ML and related topics.
Cited By |
Alert me when this paper is cited |
OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.
« Previous Article | Next Article »
OSA is a member of CrossRef.