Impact of the 2021 La Palma volcanic eruption on air quality: Insights from a multidisciplinary approach.

The La Palma 2021 volcanic eruption was the first subaerial eruption in a 50-year period in the Canary Islands (Spain), emitting ~1.8 Tg of sulphur dioxide (SO2) into the troposphere over nearly 3 months (19 September-13 December 2021), exceeding the total anthropogenic SO2 emitted from the 27 European Union countries in 2019. We conducted a comprehensive evaluation of the impact of the 2021 volcanic eruption on air quality (SO2, PM10 and PM2.5 concentrations) utilising a multidisciplinary approach, combining ground and satellite-based measurements with height-resolved aerosol and meteorological information. High concentrations of SO2, PM10 and PM2.5 were observed in La Palma (hourly mean SO2 up to ~2600 µg m-3 and also sporadically at ~140 km distance on the island of Tenerife (> 7700 µg m-3) in the free troposphere. PM10 and PM2.5 daily mean concentrations in La Palma peaked at ~380 and 60 µg m-3. Volcanic aerosols and desert dust both impacted the lower troposphere in a similar height range (~ 0-6 km) during the eruption, providing a unique opportunity to study the combined effect of both natural phenomena. The impact of the 2021 volcanic eruption on SO2 and PM concentrations was strongly influenced by the magnitude of the volcanic emissions, the injection height, the vertical stratification of the atmosphere and its seasonal dynamics. Mean daily SO2 concentrations increased during the eruption, from 38 µg m-3 (Phase I) to 92 µg m-3 (Phase II), showing an opposite temporal trend to mean daily SO2 emissions, which decreased from 34 kt (Phase I) to 7 kt (Phase II). The results of this study are relevant for emergency preparedness in all international areas at risk of volcanic eruptions; a multidisciplinary approach is key to understand the processes by which volcanic eruptions affect air quality and to mitigate and minimise impacts on the population.

The need for Pan-European automatic pollen and fungal spore monitoring: A stakeholder workshop position paper.

BACKGROUND: Information about airborne pollen concentrations is required by a range of end users, particularly from the health sector who use both observations and forecasts to diagnose and treat allergic patients. Manual methods are the standard for such measurements but, despite the range of pollen taxa that can be identified, these techniques suffer from a range of drawbacks. This includes being available at low temporal resolution (usually daily averages) and with a delay (usually 3-9 days from the measurement). Recent technological developments have made possible automatic pollen measurements, which are available at high temporal resolution and in real time, although currently only scattered in a few locations across Europe. MATERIALS & METHODS: To promote the development of an extensive network across Europe and to ensure that this network will respond to end user needs, a stakeholder workshop was organised under the auspices of the EUMETNET AutoPollen Programme. Participants discussed requirements for the groups they represented, ranging from the need for information at various spatial scales, at high temporal resolution, and for targeted services to be developed. RESULTS: The provision of real-time information is likely to lead to a notable decrease in the direct and indirect health costs associated with allergy in Europe, currently estimated between €50-150 billion/year.1 DISCUSSION & CONCLUSION: A European measurement network to meet end user requirements would thus more than pay for itself in terms of potential annual savings and provide significant impetus to research across a range of disciplines from climate science and public health to agriculture and environmental management.

A Comparative Analysis of Aerosol Optical Coefficients and Their Associated Errors Retrieved from Pure-Rotational and Vibro-Rotational Raman Lidar Signals.

This paper aims to quantify the improvement obtained with a purely rotational Raman (PRR) channel over a vibro-rotational Raman (VRR) channel, used in an aerosol lidar with elastic and Raman channels, in terms of signal-to-noise ratio (SNR), effective vertical resolution, and absolute and relative uncertainties associated to the retrieved aerosol optical (extinction and backscatter) coefficients. Measurements were made with the European Aerosol Research Lidar Network/Universitat Politècnica de Catalunya (EARLINET/UPC) multi-wavelength lidar system enabling a PRR channel at 353.9 nm, together with an already existing VRR (386.7 nm) and an elastic (354.7 nm) channels. Inversions were performed with the EARLINET Single Calculus Chain (SCC). When using PRR instead of VRR, the measurements show a gain in SNR of a factor 2.8 and about 7.6 for 3-h nighttime and daytime measurements, respectively. For 3-h nighttime (daytime) measurements the effective vertical resolution is reduced by 17% (20%), the absolute uncertainty (associated to the extinction) is divided by 2 (10) and the relative uncertainty is divided by 3 (7). During daytime, VRR extinction coefficient is retrieved in a limited height range (<2.2 km) preventing the SCC from finding a suitable calibration range in the search height range. So the advantage of using PRR instead of VRR is particularly evidenced in daytime conditions. For nighttime measurements, decreasing the time resolution from 3 to 1 h has nearly no effect on the relative performances of PRR vs. VRR.

Calculation of the Overlap Function and Associated Error of an Elastic Lidar or a Ceilometer: Cross-Comparison with a Cooperative Overlap-Corrected System.

This paper establishes the relationship between the signal of a lidar system corrected for the incomplete overlap effect and the signal of another lidar system or a ceilometer for which the overlap function is unknown. Simple mathematical relationships permit the estimation of the overlap function of the second system as well as the associated error. Several overlap functions have been retrieved with this method over a period of 1.5 years with two lidar systems of the Universitat Politècnica de Catalunya, Barcelona, Spain. The error when the overlap function reaches 1 is usually less than 7%. The temporal variability estimated over a period of 1.5 years is less than 11% in the first 1.5 km from the surface and peaks at 18% at heights between 1.7 and 2.4 km. The use of a non-appropriate overlap function in the retrieval of the backscatter coefficient yield errors up to 60% in the first 0.5 km and up to 20% above.

Considerations about the Determination of the Depolarization Calibration Profile of a Two-Telescope Lidar and Its Implications for Volume Depolarization Ratio Retrieval.

We propose a new method for calculating the volume depolarization ratio of light backscattered by the atmosphere and a lidar system that employs an auxiliary telescope to detect the depolarized component. It takes into account the possible error in the positioning of the polarizer used in the auxiliary telescope. The theory of operation is presented and then applied to a few cases for which the actual position of the polarizer is estimated, and the improvement of the volume depolarization ratio in the molecular region is quantified. In comparison to the method used before, i.e., without correction, the agreement between the volume depolarization ratio with correction and the theoretical value in the molecular region is improved by a factor of 2⁻2.5.

An Architecture Providing Depolarization Ratio Capability for a Multi-Wavelength Raman Lidar: Implementation and First Measurements.

A new architecture for the measurement of depolarization produced by atmospheric aerosols with a Raman lidar is presented. The system uses two different telescopes: one for depolarization measurements and another for total-power measurements. The system architecture and principle of operation are described. The first experimental results are also presented, corresponding to a collection of atmospheric conditions over the city of Barcelona.

Current Research in Lidar Technology Used for the Remote Sensing of Atmospheric Aerosols.

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.

Dimensionless parameterization of lidar for laser remote sensing of the atmosphere and its application to systems with SiPM and PMT detectors.

In this paper, we show a renewed approach to the generalized methodology for atmospheric lidar assessment, which uses the dimensionless parameterization as a core component. It is based on a series of our previous works where the problem of universal parameterization over many lidar technologies were described and analyzed from different points of view. The modernized dimensionless parameterization concept applied to relatively new silicon photomultiplier detectors (SiPMs) and traditional photomultiplier (PMT) detectors for remote-sensing instruments allowed predicting the lidar receiver performance with sky background available. The renewed approach can be widely used to evaluate a broad range of lidar system capabilities for a variety of lidar remote-sensing applications as well as to serve as a basis for selection of appropriate lidar system parameters for a specific application. Such a modernized methodology provides a generalized, uniform, and objective approach for evaluation of a broad range of lidar types and systems (aerosol, Raman, DIAL) operating on different targets (backscatter or topographic) and under intense sky background conditions. It can be used within the lidar community to compare different lidar instruments.

Silicon photomultiplier detector for atmospheric lidar applications.

The viability and performance of using a silicon photomultiplier (SiPM) in atmospheric lidar applications is experimentally compared against the well-established use of photomultiplier tubes. By using a modified lidar setup for simultaneous data acquisition of both types of sensors, we demonstrate that a SiPM can offer appropriate qualities for this specific application where the detection of fast, extremely low light pulses and large dynamic range signals are essential capabilities. The experimental results show that the SiPM has an appropriate behaviour offering suitable capabilities for elastic, backscatter aerosol lidars. To the best of our knowledge, this is the first study showing SiPM for atmospheric lidar applications.

Use of a field lens for improving the overlap function of a lidar system employing an optical fiber in the receiver assembly.

This paper presents a method to compute the overlap function of a lidar system in which a step-index optical fiber (or a bundle of such fibers) is used to carry the light collected by the telescope to the photoreceiver and a field lens is placed between the telescope and the optical fiber to increase the receiver field of view (FOV). The use of field lenses is a classical way to increase the FOV of radiometric systems (such as the receiving part of a lidar) when there is no numerical aperture (NA) limitation after the lens. However, when such a limitation exists, as in the case studied here, it will place a limit on the maximum attainable FOV. In the case of lidars, which have range-resolution capabilities, the limited FOV has an effect on the fraction of power coming from scattering volumes at different ranges that actually reaches the photodetector. This fraction is a function (the so-called overlap function) of the range of the scattering volume and its behavior has an impact on the accuracy of the retrievals. The application of the method developed in this paper shows that, in spite of the fiber NA limit, in practical situations the goal is attained of making the overlap function steeper and reaching higher values by using a field lens.

Practical analytical backscatter error bars for elastic one-component lidar inversion algorithm.

We present an analytical formulation to compute the total-backscatter range-dependent error bars from the well-known Klett's elastic-lidar inversion algorithm. A combined error-propagation and statistical formulation approach is used to assess inversion errors in response to the following error sources: observation noise (i.e., signal-to-noise ratio) in the reception channel, the user's uncertainty in the backscatter calibration, and in the (range-dependent) total extinction-to-backscatter ratio provided. The method is validated using a Monte Carlo procedure, where the error bars are computed by inversion of a large population of noisy generated lidar signals, for total optical depths tau < or = 5 and typical user uncertainties, all of which yield a practical tool to compute the sought-after error bars.

Quasi-analytical determination of noise-induced error limits in lidar retrieval of aerosol backscatter coefficient by the elastic, two-component algorithm.

The elastic, two-component algorithm is the most common inversion method for retrieving the aerosol backscatter coefficient from ground- or space-based backscatter lidar systems. A quasi-analytical formulation of the statistical error associated to the aerosol backscatter coefficient caused by the use of real, noise-corrupted lidar signals in the two-component algorithm is presented. The error expression depends on the signal-to-noise ratio along the inversion path and takes into account "instantaneous" effects, the effect of the signal-to-noise ratio at the range where the aerosol backscatter coefficient is being computed, as well as "memory" effects, namely, both the effect of the signal-to-noise ratio in the cell where the inversion is started and the cumulative effect of the noise between that cell and the actual cell where the aerosol backscatter coefficient is evaluated. An example is shown to illustrate how the "instantaneous" effect is reduced when averaging the noise-contaminated signal over a number of cells around the range where the inversion is started.

Effects of noise on lidar data inversion with the backward algorithm.

The lidar data-inversion algorithm widely known as the Klett method (and its more elaborate variants) has long been used to invert elastic-lidar data obtained from atmospheric sounding systems. The Klett backward algorithm has also been shown to be robust in the face of uncertainties concerning the boundary condition. Nevertheless electrical noise at the photoreceiver output unavoidably has an impact on the data-inversion process, and describing in an explicit way how it affects retrieval of the atmospheric optical coefficients can contribute to improvement in inversion quality. We examine formally the way noise disturbs backscatter-coefficient retrievals done with the Klett backward algorithm, derive a mathematical expression for the retrieved backscatter coefficient in the presence of noise affecting the signal, and assess the noise impact and suggest ways to limit it.

Variational method for the retrieval of the optical thickness and the backscatter coefficient from multiangle lidar profiles.

A variational method for retrieving the aerosol optical thickness and backscatter coefficient profiles from multiangle lidar measurements is presented and discussed. A monostatic single-wavelength low-energy lidar system was operated at different zenith angles during the Indian Ocean Experiment (INDOEX) campaign in 1999 to characterize the aerosol plumes in the Indian monsoon. The variational method was applied to lidar data to retrieve profiles of optical thickness and the backscatter coefficient for nighttime and daytime measurements. Results are obtained with an uncertainty of 10% below 3 km (nighttime) and 2.8 km (daytime) and a bias of less than 0.01. During daytime the retrieval of optical parameters is indeed limited to a lower altitude owing to the sky background signal and the atmospheric inhomogeneity. In both cases the total aerosol optical thickness is consistent (+/- 10%) with the integrated value derived from sunphotometer measurements. Backscatter-to-extinction ratios estimated in different regions by two distinct methods compared well, which proves the capability of the method to assess optical measurements and account for the altitude dependence of the phase function.