255 research outputs found

    Concept design, analysis, and Integration of the new U.P.C. multispectral lidar system

    Get PDF
    The increasing need for range-resolved aerosol and water-vapour atmospheric observation networks worldwide has given rise to multi-spectral LIDARs (Light Detection and Ranging, a synonym of laser radar) as advanced remote sensing sensors. This Ph.D. presents the design, integration and analysis of the new 6-channel multispectral elastic/Raman LIDAR for aerosol and water-vapour content monitoring developed at the Remote Sensing Lab. (RSLAB) of the Universitat Polit ecnica de Catalunya (UPC). It is well known that the combination of at least three elastic and two Raman nitrogen channels are su cient to enable retrieval of the optical and microphysical properties of aerosols with a key impact on climate change variables. The UPC lidar is part of the EARLINET (European Aerosol Research Lidar Network) -GALION (Global Atmospheric Watch Atmospheric Lidar Observation Network), a ground-based continental network including more than 28 stations. Currently, only 8 of the 28 EARLINET stations are of such advanced type. This Ph.D. speci cally focuses on: (1) Concept link-budget instrument design and overlap factor assessment. The former includes opto-atmospheric parameter modelling and assessment of backscattered power and SNR levels, and maximum system range for the di erent reception channels (3 elastic, and 2 aerosol and 1 water-vapour Raman channels, ultraviolet to near-infrared bands). The latter studies the laser-telescope crossover function (or overlap function) by means of a novel raytracing Gaussian model. The problem of overlap function computation and its near-range sensitivity for medium size aperture (f=10, f=11) bi-axial tropospheric lidar systems using both detector and ber-optics coupling alternatives at the telescope focal-plane is analysed using this new ray-tracing approach, which provides a much simpler solution than analyticalbased methods. Sensitivity to laser divergence, eld-lens and detector/ ber positions, and ber's numerical aperture is considered. (2) Design and opto-mechanical implementation of the 6-channel polychromator (i.e., the spectrally selective unit in reception). Design trade-o s concerning light collimation, end-to-end transmissivity, net channel responsivity, and homogeneous spatial light distribution onto the detectors' active area discussed. (3) System integration and validation. This third part is two fold: On one hand, fi rst-order backscatter-coe cient error bounds (a level-1 data product) for the two-component elastic lidar inversion algorithm are estimated for both random (observation noise) and systematic error sources (user's uncertainty in the backscatter-coe cient calibration, and user's uncertainty in the aerosol extinction-to-backscatter lidar ratio). On the other hand, the multispectral lidar so far integrated is described at both hardware and control software level. Statistical validation results for the new UPC lidar (today in routine operation) in the framework of SPALI-2010 intercomparison campaign are presented as part of EARLINET quality assurance / optimisation of instruments' program. The methodology developed in the rst part of this Ph.D. has successfully been applied to the speci cation case study of the IFAE/UAB lidar system, which will be installed and operated at the Cherenkov Telescope Array (CTA) observatory. Finally, specs for automated unmanned unattended lidar operation with service times close to 365/24 are presented at the end of this Ph.D. in response to the increasing demand for larger observation times and availability periods of lidar stations

    Error analysis for the lidar backward inversion algorithm

    Get PDF
    Here we depart from the inhomogeneous solution of a lidar equation using the backward inversion algorithm that is nowadays generally referred to as the Klett method. In particular, we develop an error sensitivity study that relates errors in the user-input parameters boundary extinction and exponential term in the extinction-to-backscatter relationship to errors in the inverted extinction profile. The validity of the analysis presented is limited only by the validity of application of the inversion algorithm itself, its numerical performance having been tested for optical depths in the 0.01–10 range. Toward this end, we focus on an introductory background about how uncertainties in these two parameters can apply to a family of inverted extinction profiles rather than a single profile and on its range-dependent behavior as a function of the optical thickness of the lidar inversion range. Next, we performed a mathematical study to derive the error span of the inverted extinction profile that is due to uncertainties in the above-mentioned user calibration parameters. This takes the form of upper and lower range-dependent error bounds. Finally, appropriate inversion plots are presented as application examples of this study to a parameterized set of atmospheric scenes inverted from both synthesized elastic-backscatter lidar signals and a live signal.Peer Reviewe

    Current research in lidar technology used for the remote sensing of atmospheric aerosols

    Get PDF
    Lidars are active optical remote sensing instruments with unique capabilities for atmospheric sounding. A manifold of atmospheric variables can be profiled using different types of lidar: concentration of species, wind speed, temperature, etc. Among them, measurement of the properties of aerosol particles, whose influence in many atmospheric processes is important but is still poorly stated, stands as one of the main fields of application of current lidar systems. This paper presents a review on fundamentals, technology, methodologies and state-of-the art of the lidar systems used to obtain aerosol information. Retrieval of structural (aerosol layers profiling), optical (backscatter and extinction coefficients) and microphysical (size, shape and type) properties requires however different levels of instrumental complexity; this general outlook is structured following a classification that attends these criteria. Thus, elastic systems (detection only of emitted frequencies), Raman systems (detection also of Raman frequency-shifted spectral lines), high spectral resolution lidars, systems with depolarization measurement capabilities and multi-wavelength instruments are described, and the fundamentals in which the retrieval of aerosol parameters is based is in each case detailed.Peer ReviewedPostprint (published version

    Lidar and S-band radar profiling of the atmosphere : adaptive processing for boundary-layer monitoring, optical-parameter error estimation, and application cases

    Get PDF
    This Ph.D. thesis addresses remote sensing of the atmosphere by means of lidar and S-band clear-air weather radar, and related data signal processing. Active remote sensing by means of these instruments offers unprecedented capabilities of spatial and temporal resolutions for vertical atmospheric profiling and the retrieval of key optical and physical atmospheric products in an increasing environmental regulatory framework. The first goal is this Ph.D. concerns the estimation of error bounds in the inversion of the profile of the atmospheric backscatter coefficient from elastic lidar signals (i.e., without wavelength shift in reception when interacting with atmospheric scatterers) by means of the two-component inversion algorithm (the so-called Klett-Fernald-Sasano¿s algorithm). This objective departs from previous works at the Remote Sensing Lab. (RSLab) of the Universitat Politècnica de Catalunya (UPC) and derives first-order error-propagated bounds (approximate) and total-increment bounds (exact). As distinctive feature in the state of the art, the error bounds merge into a single body both systematic (i.e., user-calibration inputs) and random error sources (finite signal-to-noise ratio, SNR) yielding an explicit mathematical form. The second goal, central to this Ph.D., tackles retrieval of the Atmospheric Boundary Layer Height (ABLH) from elastic lidar and S-band Frequency-Modulated Continuous-Wave (FMCW) radar observations by using adaptive techniques based on the Extended Kalman Filter (EKF). The filter is based on morphological modelling of the Mixing-Layer-to-Free-Troposphere transition and continuous estimation of the noise covariance information. In the lidar-EKF realization the proposed technique is shown to outperform classic ABLH estimators such as those based on derivative techniques, thresholded decision, or the variance centroid method. The EKF formulation is applied to both ceilometer and UPC lidar records in high- and low-SNR scenes. The lidar-EKF approach is re-formulated and successfully extended to S-band radar scenes (Bragg¿s scattering) in presence of interferent noise sources (Rayleigh scattering from e.g., insects and birds). In this context, the FMCW feature enables the range-resolved capability. EKF-lidar and EKF-radar ABLH estimates are cross-examined from field campaign results. Finally, the third goal deals with exploitation of the existing UPC lidar station: In a first introductory part, a modified algorithm for enhancing the dynamic range of elastic lidar channels by ¿gluing¿ analog and photon-counting data records is formulated. In a second part, two case examples (including application of the gluing algorithm) are presented to illustrate the capabilities of the UPC lidar in networked atmospheric observation of two recent volcano eruption events as part of the EARLINET (European Aerosol Research Lidar Network). The latter is part of GALION (Global Atmospheric Watch Atmospheric Lidar Observation Network)-GEOSS (Global Earth Observation System of Systems) framework.La tesis doctoral aborda la teledetecció atmosfèrica amb tècniques lidar i radar (banda S) i llur tractament del senyal. La teledetecció activa amb aquests instruments ofereix resolucions espacials i temporals sense precedents en la perfilometria vertical de l'atmosfera i recuperació de productes de dades òptics i físics atmosfèrics en un marc de creixent regulació mediambiental. El primer objectiu d'aquesta tesi concerneix l'estimació de cotes d'error en la inversió del perfil del coeficient de retrodispersió atmosfèrica a partir de senyals lidar de tipus elàstic (és a dir, sense desplaçament de la longitud d'ona en recepció al interactuar amb els dispersors atmosfèrics) mitjançant l'algorisme d'inversió de dues components de Klett-Fernald-Sasano. Aquest objectiu parteix de treballs previs en el Remote Sensing Lab. (RSLab) de la Universitat Politècnica de Catalunya (UPC) i permet obtenir cotes de primer ordre (aproximades) basades en propagació d'errors i cotes (exactes) basades en el increment total de l'error. Característica diferencial en front l'estat de l'art és l'assimilació d'errors sistemàtics (per exemple, entrades de cal.libració d'usuari) i aleatoris (relació senyal-soroll, SNR, finita) en forma matemàtica explícita. El segon objectiu, central de la tesis, aborda l'estimació de l'altura de la capa límit atmosfèrica (ABLH) a partir de senyal lidar elàstics i d'observacions radar en banda S (ona continua amb modulació en freqüència, FMCW) utilitzant tècniques adaptatives basades en filtrat estès de Kalman (EKF). El filtre es basa en modelat morfològic de la transició atmosfèrica entre la capa de mescla i la troposfera lliure i en l'estimació continua de la informació de covariança del soroll. En el prototipus lidar-EKF la tècnica proposada millora clarament les tècniques clàssiques d'estimació de la ABLH como són les basades en mètodes derivatius, decisió de llindar, o el mètode de la variança-centroide. La formulació EKF s'aplica tant a mesures procedents de ceilòmetres lidar como de la pròpia estació lidar UPC en escenes d'alta i baixa SNR. Addicionalment, l'enfoc lidar-EKF es reformula i s'estén amb èxit a escenes radar en banda S (dispersió Bragg) en presència de fonts de soroll interferent (dispersió Rayleigh de, per exemple, insectes i ocells). En aquest context, la característica FMCW permet la capacitat de resolució en distància. L'estimació de la ABLH amb els prototipus lidar-EKF i radar-EKF s'intercompara en campanyes de mesura. Finalment, el tercer objectiu atén a l'explotació de l'estació lidar UPC existent: En una primera part introductòria, es formula un algorisme modificat de "gluing" per a la millora del marge dinàmic de canals lidar elàstics mitjançant combinació (o "enganxat") de senyals lidar adquirits analògicament i amb foto-comptatge. En una segona part, es presenten dos exemples (incloent l'aplicació de l'algorisme de "gluing") que il.lustren les capacitats del lidar de la UPC en l'observació atmosfèrica de dos recents erupcions volcàniques des de la xarxa d'observació EARLINET (European Aerosol Research Lidar Network). Aquesta última és part de GALION (Global Atmospheric Watch Atmospheric Lidar Observation Network)-GEOSS (Global Earth Observation System of Systems)

    Optical-microphysical Properties of Saharan Dust Aerosols and Composition Relationship Using a Multi-wavelength Raman Lidar, in Situ Sensors and Modelling: a Case Study Analysis

    Get PDF
    A strong Saharan dust event that occurred over the city of Athens, Greece (37.9° N, 23.6° E) between 27 March and 3 April 2009 was followed by a synergy of three instruments: a 6-wavelength Raman lidar, a CIMEL sun-sky radiometer and the MODIS sensor. The BSC-DREAM model was used to forecast the dust event and to simulate the vertical profiles of the aerosol concentration. Due to mixture of dust particles with low clouds during most of the reported period, the dust event could be followed by the lidar only during the cloud-free day of 2 April 2009. The lidar data obtained were used to retrieve the vertical profile of the optical (extinction and backscatter coefficients) properties of aerosols in the troposphere. The aerosol optical depth (AOD) values derived from the CIMEL ranged from 0.33-0.91 (355 nm) to 0.18-0.60 (532 nm), while the lidar ratio (LR) values retrieved from the Raman lidar ranged within 75-100 sr (355 nm) and 45-75 sr (532 nm). Inside a selected dust layer region, between 1.8 and 3.5 km height, mean LR values were 83 ± 7 and 54 ± 7 sr, at 355 and 532 nm, respectively, while the Ångström-backscatter-related (ABR 355/532) and Ångström-extinction-related (AER 355/532) were found larger than 1 (1.17 ± 0.08 and 1.11 ± 0.02, respectively), indicating mixing of dust with other particles. Additionally, a retrieval technique representing dust as a mixture of spheres and spheroids was used to derive the mean aerosol microphysical properties (mean and effective radius, number, surface and volume density, and mean refractive index) inside the selected atmospheric layers. Thus, the mean value of the retrieved refractive index was found to be 1.49( ± 0.10) + 0.007( ± 0.007)i, and that of the effective radiuses was 0.30 ± 0.18 μm. The final data set of the aerosol optical and microphysical properties along with the water vapor profiles obtained by Raman lidar were incorporated into the ISORROPIA II model to provide a possible aerosol composition consistent with the retrieved refractive index values. Thus, the inferred chemical properties showed 12-40% of dust content, sulfate composition of 16-60%, and organic carbon content of 15-64%, indicating a possible mixing of dust with haze and smoke. PM10 concentrations levels, PM10 composition results and SEM-EDX (Scanning Electron Microscope-Energy Dispersive X-ray) analysis results on sizes and mineralogy of particles from samples during the Saharan dust transport event were used to evaluate the retrieval

    Separation of the optical and mass features of particle components in different aerosol mixtures by using POLIPHON retrievals in synergy with continuous polarized Micro-Pulse Lidar (P-MPL) measurements

    Get PDF
    The application of the POLIPHON (POlarization-LIdar PHOtometer Networking) method is presented for the first time in synergy with continuous 24/7 polarized Micro-Pulse Lidar (P-MPL) measurements to derive the vertical separation of two or three particle components in different aerosol mixtures, and the retrieval of their particular optical properties. The procedure of extinction-to-mass conversion, together with an analysis of the mass extinction efficiency (MEE) parameter, is described, and the relative mass contribution of each aerosol component is also derived in a further step. The general POLIPHON algorithm is based on the specific particle linear depolarization ratio given for different types of aerosols and can be run in either 1-step (POL-1) or 2 steps (POL-2) versions with dependence on either the 2- or 3-component separation. In order to illustrate this procedure, aerosol mixing cases observed over Barcelona (NE Spain) are selected: a dust event on 5 July 2016, smoke plumes detected on 23 May 2016 and a pollination episode observed on 23 March 2016. In particular, the 3-component separation is just applied for the dust case: a combined POL-1 with POL-2 procedure (POL-1/2) is used, and additionally the fine-dust contribution to the total fine mode (fine dust plus non-dust aerosols) is estimated. The high dust impact before 12:00UTC yields a mean mass loading of 0.6±0.1gm-2 due to the prevalence of Saharan coarse-dust particles. After that time, the mean mass loading is reduced by two-thirds, showing a rather weak dust incidence. In the smoke case, the arrival of fine biomass-burning particles is detected at altitudes as high as 7km. The smoke particles, probably mixed with less depolarizing non-smoke aerosols, are observed in air masses, having their origin from either North American fires or the Arctic area, as reported by HYSPLIT back-trajectory analysis. The particle linear depolarization ratio for smoke shows values in the 0.10–0.15 range and even higher at given times, and the daily mean smoke mass loading is 0.017±0.008gm-2, around 3% of that found for the dust event. Pollen particles are detected up to 1.5km in height from 10:00UTC during an intense pollination event with a particle linear depolarization ratio ranging between 0.10 and 0.15. The maximal mass loading of Platanus pollen particles is 0.011±0.003gm-2, representing around 2% of the dust loading during the higher dust incidence. Regarding the MEE derived for each aerosol component, their values are in agreement with others referenced in the literature for the specific aerosol types examined in this work: 0.5±0.1 and 1.7±0.2m2g-1 are found for coarse and fine dust particles, 4.5±1.4m2g-1 is derived for smoke and 2.4±0.5m2g-1 for non-smoke aerosols with Arctic origin, and a MEE of 2.4±0.8m2g-1 is obtained for pollen particles, though it can reach higher or lower values depending on predominantly smaller or larger pollen grain sizes. Results reveal the high potential of the P-MPL system, a simple polarization-sensitive elastic backscatter lidar working in a 24/7 operation mode, to retrieve the relative optical and mass contributions of each aerosol component throughout the day, reflecting the daily variability of their properties. In fact, this procedure can be simply implemented in other P-MPLs that also operate within the worldwide Micro-Pulse Lidar Network (MPLNET), thus extending the aerosol discrimination at a global scale. Moreover, the method has the advantage of also being relatively easily applicable to space-borne lidars with an equivalent configuration such as the ongoing Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP) on board NASA CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) and the forthcoming Atmospheric Lidar (ATLID) on board the ESA EarthCARE mission.Peer ReviewedPostprint (published version

    Separation of the optical and mass features of particle components in different aerosol mixtures by using POLIPHON retrievals in synergy with continuous polarized Micro-Pulse Lidar (P-MPL) measurements

    Get PDF
    The application of the POLIPHON (POlarization-LIdar PHOtometer Networking) method is presented for the first time in synergy with continuous 24/7 polarized Micro-Pulse Lidar (P-MPL) measurements to derive the vertical separation of two or three particle components in different aerosol mixtures, and the retrieval of their particular optical properties. The procedure of extinction-to-mass conversion, together with an analysis of the mass extinction efficiency (MEE) parameter, is described, and the relative mass contribution of each aerosol component is also derived in a further step. The general POLIPHON algorithm is based on the specific particle linear depolarization ratio given for different types of aerosols and can be run in either 1-step (POL-1) or 2 steps (POL-2) versions with dependence on either the 2- or 3-component separation. In order to illustrate this procedure, aerosol mixing cases observed over Barcelona (NE Spain) are selected: a dust event on 5 July 2016, smoke plumes detected on 23 May 2016 and a pollination episode observed on 23 March 2016. In particular, the 3-component separation is just applied for the dust case: a combined POL-1 with POL-2 procedure (POL-1/2) is used, and additionally the fine-dust contribution to the total fine mode (fine dust plus non-dust aerosols) is estimated. The high dust impact before 12:00 UTC yields a mean mass loading of 0.6±0.1 g m'2 due to the prevalence of Saharan coarse-dust particles. After that time, the mean mass loading is reduced by two-thirds, showing a rather weak dust incidence. In the smoke case, the arrival of fine biomass-burning particles is detected at altitudes as high as 7 km. The smoke particles, probably mixed with less depolarizing non-smoke aerosols, are observed in air masses, having their origin from either North American fires or the Arctic area, as reported by HYSPLIT back-trajectory analysis. The particle linear depolarization ratio for smoke shows values in the 0.10-0.15 range and even higher at given times, and the daily mean smoke mass loading is 0.017±0.008 g m'2, around 3 % of that found for the dust event. Pollen particles are detected up to 1.5 km in height from 10:00 UTC during an intense pollination event with a particle linear depolarization ratio ranging between 0.10 and 0.15. The maximal mass loading of Platanus pollen particles is 0.011±0.003 g m'2, representing around 2 % of the dust loading during the higher dust incidence. Regarding the MEE derived for each aerosol component, their values are in agreement with others referenced in the literature for the specific aerosol types examined in this work: 0.5±0.1 and 1.7±0.2 m2 g'1 are found for coarse and fine dust particles, 4.5±1.4 m2 g'1 is derived for smoke and 2.4±0.5 m2 g'1 for non-smoke aerosols with Arctic origin, and a MEE of 2.4±0.8 m2 g'1 is obtained for pollen particles, though it can reach higher or lower values depending on predominantly smaller or larger pollen grain sizes. Results reveal the high potential of the P-MPL system, a simple polarization-sensitive elastic backscatter lidar working in a 24/7 operation mode, to retrieve the relative optical and mass contributions of each aerosol component throughout the day, reflecting the daily variability of their properties. In fact, this procedure can be simply implemented in other P-MPLs that also operate within the worldwide Micro-Pulse Lidar Network (MPLNET), thus extending the aerosol discrimination at a global scale. Moreover, the method has the advantage of also being relatively easily applicable to space-borne lidars with an equivalent configuration such as the ongoing Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP) on board NASA CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) and the forthcoming Atmospheric Lidar (ATLID) on board the ESA EarthCARE mission

    On the Sensitivity of a Ground-Based Tropospheric Lidar to Aitken Mode Particles in the Upper Troposphere

    Get PDF
    Airborne observations have shown high concentrations of ultrafine aerosols in the Amazon upper troposphere (UT), which are key for replenishing the planetary boundary layer (PBL) with cloud condensation nuclei that sustain the “green ocean” clouds. Given their climatic relevance, longterm observations are needed, but aircraft measurements are only available in short-term campaigns. Alternatively, continuous observations of the aerosol vertical structure could be performed by a lidar (acronym for “light detection and ranging”) system in long-term campaigns. Here we assess whether a ground-based tropospheric lidar system could detect these ultrafine UT aerosols. To this aim, we simulated the lidar signal of a real instrument and then varied the instrument’s efficiency and the UT-particle concentration to determine under which conditions the detection is possible. Optical properties were computed with a Mie code based on the size distributions and numerical concentration profiles measured by the aircraft, and on the refractive indexes inverted from AERONET measurements. The aerosol optical depth (AOD) was retrieved by inverting the elastic lidar signal, and a statistical test was applied to evaluate the detection of the UT-aerosol layer. Our results indicate that, for the instrument we simulated, a 55-fold increase in the signal-to-noise ratio (SNR) is required for a 100% detection rate. This could be achieved by simultaneously time averaging over 30 min and spatially averaging to vertical bin lengths of 375 m, or by modifying the hardware. We repeated the analysis for under- and overestimated aerosol lidar ratio (Laer), and found that possible systematic errors did not affect the detection rate. Further studies are necessary to assess whether such longtime averages are feasible in the Amazon region (given the very high cloud cover), and to design a hardware upgrade. Although simulations and analyses here were based on a particular instrument and for the presence of new organic particles in the Amazonian upper troposphere, our methodology and results are general and applicable to other instruments and sites.Conselho Nacional de Desenvolvimento Cientifico e Tecnologico (CNPQ) 132402/2020-3 308682/2017-

    Design, implementation, and characterisation of a novel lidar ceilometer

    Get PDF
    A novel lidar ceilometer prototype based on divided lens optics has been designed, built, characterised, and tested. The primary applications for this manufacturable ground-based sensor are the determination of cloud base height and the measurement of vertical visibility. First, the design, which was developed in order to achieve superior performance at a low cost, is described in detail, along with the process used to develop it. The primary design considerations of optical signal to noise ratio, range-dependent overlap of the transmitter and receiver channels, and manufacturability, were balanced to develop an instrument with good signal to noise ratio, fast turn-on of overlap for detection of close range returns, and a minimised number of optical components and simplicity of assembly for cost control purposes. Second, a novel imaging method for characterisation of transmitter-receiver overlap as a function of range is described and applied to the instrument. The method is validated by an alternative experimental method and a geometric calculation that is specific to the unique geometry of the instrument. These techniques allow the calibration of close range detection sensitivity in order to acquire information prior to full overlap. Finally, signal processing methods used to automate the detection process are described. A novel two-part cloud base detection algorithm has been developed which combines extinction-derived visibility thresholds in the inverted cloud return signal with feature detection on the raw signal. In addition, standard approaches for determination of visibility based on an iterative far boundary inversion method, and calibration of attenuated backscatter profile using returns from a fully-attenuating water cloud, have been applied to the prototype. The prototype design, characterisation, and signal processing have been shown to be appropriate for implementation into a commercial instrument. The work that has been carried out provides a platform upon which a wide range of further work can be built
    corecore