Volcanic emissions from the Tajogaite volcano, located on the Cumbre Vieja edifice on the island of La Palma (Canary Islands, Spain), caused significant public health and aviation disruptions throughout the eruption (19 September–13 December 2021, officially declared over on 25 December). Nonetheless, it is considered the most significant volcanic event in Europe over the past 75 years due to the substantial amount of SO2 released into the atmosphere. The Instituto Geográfico Nacional (IGN), the authority responsible for volcano surveillance in Spain, implemented extensive operational monitoring to track volcanic activity and to provide a robust estimation of the volcanic plume height using a video-surveillance network. In parallel, the State Meteorological Agency of Spain (AEMET), in partnership with other Spanish ACTRIS (Aerosol, Clouds, and Trace Gases Research Infrastructure) members and collaborating institutions, conducted an unprecedented instrumental deployment to evaluate the impacts of this volcanic event on atmospheric composition. This effort included a network of aerosol profilers surrounding the volcano. A total of four profiling instruments were installed on La Palma: one MPL-4B lidar and three ceilometers. Additionally, a pre-existing Raman lidar on the island contributed valuable data to this study. These efforts are undertaken due to the importance of monitoring volcanic plume height in terms of air quality (necessary for the implementation of effective civil protection policies), volcanic activity surveillance (for tracking and forecasting eruptive behaviour), and, from a scientific perspective, for improving our understanding of the climatic and radiative impacts of this type of aerosol.
In this study, the eruptive process was characterised in terms of the altitude of the dispersive volcanic plume (hd), measured by both IGN and AEMET-ACTRIS, and the altitude of the eruptive column (hec), measured by IGN. Modulating factors such as seismicity and meteorological conditions were also analysed. The consistency between the two independent and complementary datasets (hd,IGN and hd,AEMET) was assessed throughout the eruption (mean difference of 258.6 m).
Our results confirmed the existence of three distinct eruptive phases, encompassing a range of styles from Strombolian explosive to effusive activity. While these phases have been characterised in previous studies, the results of the present work provide complementary information and novel insights from an alternative observational approach, which may be of use in future volcanic crises and will be applied to operational surveillance during such events.
A subsequent comparison of hd,AEMET with the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) aerosol layer height product (ALHCALIOP) revealed a systematic underestimation by the satellite product, with a mean difference of 392.2 m.
Finally, the impact of using hec in estimating SO2 emissions from the NASA MSVOLSO2L4 satellite-based product was evaluated. When a fixed (standard) plume altitude of 8 km was used instead of the observed hec, the total SO2 emission was significantly underestimated by an average of 56.2 %, and by up to 84.7 %. These findings underscore the importance of accurately determining the volcanic plume height when deriving SO2 emissions from satellite data.
Introduction
Anthropogenic aerosols, as one of the major drivers of climate change, have attracted significant scientific attention for decades. Natural aerosols, such as dust, sea salt, and volcanic particles, have also been thoroughly studied in recent decades, particularly in relation to their radiative forcing . Volcanic aerosols, both tropospheric and stratospheric, have a primary impact on atmospheric chemistry, climate, and the radiative budget. Moreover, they play a critical role in public health, civil aviation, the economy, and ecosystems .
Volcanoes primarily emit gaseous species, including water vapour (H2O), carbon dioxide (CO2), and sulfur dioxide (SO2), the latter being the most abundant gas released during volcanic activity . Secondary sulfate aerosols are formed through in-plume conversion of SO2 into sulfuric acid droplets via gas-phase oxidation , along with other sulfate-bearing compounds . Taking into account that the lifetime of sulfate aerosols is significantly longer than that of primary volcanic emissions (typically around 1.3 weeks in the troposphere and several years when injected into the stratosphere), secondary volcanic aerosols can be transported over long distances and may exert a global impact on the climate system , such as increasing the Earth’s albedo or contributing to ozone depletion in the stratosphere . Such aerosols can also play a key role as cloud condensation nuclei (CCN) and ice-nucleating particles (INPs), influencing cloud formation and their radiative properties . The amount of incoming solar radiation scattered from these emitted atmospheric components, and therefore the radiative effect and impact of the volcanic aerosols, depend on the location of the volcanic layer, the size of the emitted particles , and also the nature of the volcanic aerosol.
Although stratospheric volcanic aerosols have a more significant and prolonged impact on climate, tropospheric volcanic aerosols originating from smaller eruptions and sustained magmatic or hydrothermal degassing, despite their relatively modest scale, can also substantially affect air quality on local scales. This is because they may be injected into either the planetary boundary layer or the free troposphere, depending on the interplay between plume height and the atmospheric vertical structure . Such volcanic activity remains largely understudied, primarily due to the lack of detailed “near-source” characterisations . Consequently, the atmospheric and climatic effects of weak volcanic activity are still poorly represented in current climate models .
Monitoring either tropospheric or stratospheric volcanic aerosol emissions is inherently challenging, as volcanic eruptions are natural events that are highly unpredictable, short-lived, and often occur in remote or inaccessible locations. These factors, together with the difficulties of deploying scientific instrumentation in hazardous environments, make near-source observations particularly scarce. The eruptions of the Eyjafjallajökull volcano in Iceland during April and May 2010 provided the scientific community with an unprecedented opportunity to study, in near-real time, the role of volcanic aerosols in the climate system. The event also had a major impact on air traffic and the global economy, as the volcanic plume was transported over Europe and references therein. As noted by , this eruption triggered significant advancements in lidar technology for monitoring the height of volcanic plumes and promoted its broader deployment by research and environmental institutions, given the suitability of lidar systems for detecting transported aerosol plumes at high altitudes and their ability to assess aerosol type through depolarisation measurements. Assimilation of ground-based lidar profiles into aerosol transport models has also been shown to be beneficial for volcanic monitoring and forecasting .
Other ground-based techniques, including video-surveillance systems, weather radars, infrasound, and lightning detection, have also been employed to estimate volcanic plume-top height . Among these, video-surveillance systems offer a particularly promising and easily automatable solution for near-real-time estimation of hec during emergency events .
Despite the increasing use of lidar technology for volcanic plume monitoring, no comprehensive comparison has yet been carried out with other observation techniques more typically employed within the volcanological community. This limits our ability to assess the consistency and complementarity of the different approaches.
Another key parameter used to monitor volcanic activity and to estimate the radiative and climatic impacts of volcanic eruptions is the SO2 mass loading (or emitted flux) . Sulfur emission rates during eruptions can be quantified through direct observations from various platforms , including satellite remote sensing and in situ measurements . Satellite-based hyperspectral spectrometers operating in the ultraviolet (UV) have provided frequent and increasingly accurate observations of global SO2 levels. These observations rely on retrieval algorithms that process backscattered radiance measurements to estimate both SO2 mass loadings and plume height . According to and , the greatest source of uncertainty in estimating SO2 mass loadings in the lower troposphere arises from the limited a priori knowledge of volcanic plume altitude. SO2 retrievals are more accurate for plumes located above the SO2 cloud and snow/ice layers, with overall uncertainties typically ranging from 20 %–30 % .
The 2021 volcanic eruption on La Palma (Canary Islands, Spain) served as a natural laboratory for studying the eruptive process and its impacts from multiple perspectives e.g.,. The evolution and long-range transport of the volcanic plume have been the subject of numerous studies, including its dispersion toward the Iberian Peninsula , southern France , and Cape Verde . In this context, this study presents the unprecedented instrumental coverage deployed during the Tajogaite (Cumbre Vieja) eruption (19 September to 25 December 2021) by the Instituto Geográfico Nacional (IGN), the Spanish State Meteorological Agency (AEMET), and other Spanish members of ACTRIS (Aerosol, Clouds and Trace Gases Research Infrastructure) to monitor the atmospheric impact of this rare event. As the institution responsible for volcano surveillance in Spain, IGN implemented extensive monitoring from the beginning of the eruption, including robust estimations of the volcanic plume height using a video-surveillance system. This quantification was incorporated into the PEVOLCA (Steering Committee of the Special Plan for Civil Protection and Attention to Emergencies due to Volcanic Risk) reports , the VONA (Volcano Observatory Notice for Aviation) alerts, and the regular reports submitted to the Toulouse VAAC (Volcanic Ash Advisory Centre) during the crisis . In parallel, AEMET and ACTRIS (hereafter AEMET-ACTRIS) established a network of aerosol profilers around the volcano to provide a complementary and reliable estimation of the dispersive plume altitude.
The information presented in this study has been used to shed light on two key parameters in volcanic activity monitoring. First, we aim to provide a comprehensive comparison of independent methods for measuring volcanic plume height. By understanding the intrinsic limitations of each approach, these methods could be employed as complementary tools in future volcanic crises. Additionally, this study seeks to evaluate the impact of discrepancies in plume height measurements on the quantification of volcanic emissions from satellite observations, thereby quantifying the importance of accurate plume height determination for reliable emission estimates.
The manuscript is organised as follows. Section describes the main features of the volcanic eruption. Section details the IGN and AEMET-ACTRIS monitoring networks, as well as the auxiliary information used in this study. The main results are presented in Sect. , including the comparison between the two monitoring networks (Sect. ) and the description of the volcanic event from volcanic plume altitudes and modulating factors (Sect. ). The evaluation of satellite-based products is included in Sect. : CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations) CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization) aerosol height (Sect. ), and the estimation of SO2 volcanic emissions using multi-satellite UV-based observations (Sect. ). Finally, Sect. presents the summary and conclusions that are extracted from this work.
Tajogaite 2021 volcanic eruption
The eruption began on 19 September 2021 at the Cumbre Vieja volcanic edifice. The initial phase was a fissural eruption, and over the 85 d of eruption duration, significant changes were detected in the activity. PEVOLCA Committee declared officially the end of the eruption on 25 December (which ended on 13 December from a seismic perspective) . In September 2022, nearly one year after the eruption began, the volcano on Cumbre Vieja was officially named Tajogaite, a traditional toponym of Guanche origin used by locals to refer to the upper part of the affected area.
This eruption is considered the most significant volcanic event in Europe over the past 75 years due to the substantial amount of SO2 released into the atmosphere, and the extensive damage caused by lava flows . It is also regarded as the longest-lasting historical eruption on the island. The eruption led to the evacuation of more than 7000 residents, with severe consequences for public health and a profound impact on the island's economy. The eruption caused long-lasting impacts due to extensive lava fields (1219 ha), widespread damage to infrastructure, homes, and farmland, and ongoing volcanic gas emissions in a region heavily dependent on tourism. Aviation was severely affected, with 26 % of scheduled flights at La Palma Airport cancelled – 34 % due to airport closures from ash accumulation and 66 % due to ash in the airspace. Initial estimates place total economic losses at over USD 1025 million .
Multiple summit vents emerged during the nearly three months of volcanic activity, exhibiting a range of eruptive styles – from Strombolian explosions to effusive phases – resulting in violent ash-rich eruptions, lava fountains, gas emissions, and extensive lava flows. The main volcanic edifice, a cinder cone that reached an altitude of 1121 m above sea level (a.s.l.), was formed in the northwestern sector of the Cumbre Vieja rift .
Although the eruption was predominantly characterised by Strombolian basaltic activity, different eruptive phases were observed throughout the entire eruption period . As noted by , these phases were distinguishable based on variations in surface and satellite-derived SO2 emission rates. According to their results, the initial phase of the eruption exhibited the highest SO2 and ash emissions, with an accumulated SO2 output of 1.59 Mt. This was followed by a marked decline in emissions during a second phase, beginning around 7 November, with an accumulated SO2 emission of 0.25 Mt. According to , nine different eruptive phases can be identified in terms of seismological data (epicenters, hypocenters and magnitude/depth of earthquakes). attended to the variation in the number of events with magnitudes between 2 and 3 mbLg, as well as changes in the Real-time Seismic Amplitude Measurement (RSAM) values using the seismic pattern recorded by the IGN monitoring network (https://www.ign.es/web/ign/portal/sis-catalogo-terremotos, last access: 15 March 2025). Such variations can be interpreted as upwards migration of fresh, SO2-rich magma batches and related degassing processes through the eruptive vents .
Instruments and methodsIGN video-surveillance monitoring network
The IGN, as the competent institution for volcano surveillance in Spain, monitored volcanic events in terms of seismic activity, ground deformation, geochemical and gravimetric parameters, and volcanic plume behaviour . In this context, a monitoring network based on video-surveillance cameras was implemented to track the altitude of the volcanic plume. This methodology provided the quantitative characterisation of the volcanic plume required by the PEVOLCA Committees during the crisis. The results were also incorporated into the VONA alerts and the regular reports submitted to the Toulouse VAAC .
The monitoring network consisted of calibrated webcams that provided scaled images of the volcanic plume. A pre-existing webcam from the E.U. EELabs project (https://www.eelabs.eu/, last access: 1 January 2024), coordinated by the Instituto de Astrofísica de Canarias (IAC, https://www.iac.es/, last access: 1 January 2024) and located at the Roque de Los Muchachos Observatory (28.76° N, 17.88° W, 2423 m a.s.l.; approximately 16 km north of the eruption site), was selected as the reference calibrated webcam. This camera had to meet several criteria to be suitable for measuring the altitude of the volcanic plume. A favourable location was essential: the camera needed to be situated at a high altitude above sea level and far enough from the eruption to capture the full extent of the plume, while minimising the effect of plume inclination and the obstruction by volcanic material. The camera's position and internal settings must remain unchanged during the eruptive period. Furthermore, it was necessary to ensure continuous data recording and easy access to the image archive. The IAC camera met all these requirements and had already been recording prior to the eruption. Additionally, three more webcams were installed by IGN personnel at different locations surrounding the volcanic vent: Time (28.66° N, 17.94° W), Los Llanos de Aridane (28.67° N, 17.92° W), and Fuencaliente (28.50° N, 17.85° W). These cameras were used to measure the height of the volcanic plume from different angles, allowing for a better understanding of the influence of wind on plume dispersion.
In order to achieve the objective of detecting the altitude of the volcanic plume, it is essential to determine the scale of the images based on a specific distance, particularly the distance from the camera to the volcanic vent. By employing geodetic techniques such as trigonometric levelling, we obtained different altitudes above sea level along the eruptive column at specific moments and selected locations. Once these altitudes were retrieved, they were compared with the corresponding points identified in the images.
To perform direct measurements on the images, a calibrated scale needed to be established (see Fig. ). The calibration process involved multiple simultaneous measurements of the same reference points using both image analysis and geodetic methods. A trigonometrical station, located approximately 30 m from the camera, was used for this purpose. The resulting scale was used to generate a grid that was automatically superimposed on the images, enabling real-time measurements of the volcanic plume height on demand. This grid also included wind direction and speed at different altitudes, as presented in Fig. , which was crucial for accurate estimation of the column height.
Calibrated image taken from the IGN-IAC webcam on 13 December 2021 with grid meteorological information from AEMET (wind speed, WS in m s−1) and wind direction (in terms of Cardinal Points, PC).
This procedure was also applied to the camera installed at the Time station. For the other two webcams, located in Los Llanos de Aridane and Fuencaliente, calibration was performed by comparing their images to those from the calibrated reference camera at Roque de Los Muchachos (RMO), taken at the same moments. Unfortunately, the Fuencaliente camera was installed near the end of the eruption and provided limited useful data.
This technique has its limitations. Since it relies on visual observations from cameras, it is generally not effective during nighttime. However, during the most intense phases of an eruption, it was possible to observe the eruptive column even at night. On days with heavy cloud cover or fog, identifying the top of the eruptive column proved challenging. Furthermore, during certain stages of the eruption, particularly in the final weeks when ash emissions were weak and the volcano displayed a pulsating behaviour, producing discrete ash puffs every few seconds, the visibility was significantly compromised.
Using this methodology, the height reached by the eruptive column (hec) and the dispersive plume (hd) can be calculated, the latter referring to the altitude attained by the volcanic plume as it is transported and modified by thermodynamic processes, including atmospheric stratification, wind direction, and turbulence.
The IGN altitude database (hec,IGN and hd,IGN) consists of 344 individual altitude measurements corresponding to the eruptive column and dispersive volcanic plume heights recorded between 20 September and 13 December 2021. Each observation was performed at the time when a significant change in the height of the eruptive column was detected. Therefore, it can be assumed that both series remained relatively stable between consecutive measurements.
AEMET-ACTRIS profiling monitoring network
Four stations were strategically deployed in record time (between 11 and 40 d after the eruption) around the Tajogaite-Cumbre Vieja volcano by AEMET-ACTRIS members in Spain in collaboration with other institutions (Fig. ). These stations, equipped with various types of cloud and aerosol profilers, enabled continuous tracking of the volcanic plume under variable wind conditions throughout the entire eruption, as well as providing ground-level air quality information and plume height for PEVOLCA reports. In addition to these, one station that was already installed in La Palma prior to the eruption was also added to the profiling network.
The characteristics of the five sites, as shown in Fig. b, are as follows:
Roque de los Muchachos (RMO, 28.75° N, 17.88° W, 2423 m a.s.l.): this high-altitude observatory, part of the Instituto de Astrofísica de Canarias (IAC), is dedicated to high-quality astrophysical observations. It hosts telescopes and other astronomical instruments from nineteen countries, in addition to an AERONET (Aerosol Robotic Network, https://aeronet.gsfc.nasa.gov, last access: 15 March 2025) station, which provides valuable aerosol properties information for this site. RMO is located 16.5 km from the eruptive vents.
Tazacorte (TAZ, 28.64° N, 17.93° W, 140 m a.s.l.): situated approximately 4 km from the eruptive vents in the town of Tazacorte (population ∼ 4500) on the eastern flank of La Palma, this station offers an excellent opportunity to study the aerosol layer near ground level when the predominant trade winds favor the dispersion of the volcanic plume westward.
Fuencaliente (FUE, 28.49° N, 17.85° W, 630 m a.s.l.): located on the southern side of La Palma, with a population of fewer than 2000 people, this station is about 10 km from the volcanic vents. Positioned in a rural background site, it provides a valuable opportunity to detect the volcanic plume under the prevailing trade wind regime.
La Palma Airport (LPA, 28.62° N, 17.75° W, 56 m a.s.l.): this station is crucial for providing civil aviation with reliable information regarding the potential presence of volcanic ash, which can severely impact aircraft fuselages, sensors, and engines.
El Paso (EPA, 28.65° N, 17.87° W, 700 m): located in the central-northern part of La Palma, the municipality of El Paso has a population of approximately 7700. The station here is the closest to the volcano, at about 3 km from the volcanic vents, allowing for early detection of the volcanic plume during its initial stages. However, it also impacts the instruments in terms of the quality of the measurement, mainly in the most intense stages of the eruption.
By combining information from these five sites, a single hd,AEMET time series was constructed, given that the proximity of the stations to the volcano and to each other ensures that no significant dispersive processes occur either horizontally or vertically, nor are intra-plume chemical reactions expected. This assumption is supported by previous findings on intra-plume interactions e.g.,. A total of 137 altitudes of the dispersive column (hd,AEMET) have been included in the comparison analysis, coincident with those values measured by IGN. hd,AEMET observations from these five profilers have been selected for the comparison with IGN's database as those lying in a time interval of ±5 min around IGN observation time. In cases where multiple AEMET heights from different profilers and sites coincide with the IGN height, the one most influenced by the volcanic plume, based on the prevailing wind direction, is selected. Taking into account the inherent problems in retrieving hd,AEMET from this on-site profiler network (signal blocking by the dense volcanic plume, high signal-to-noise ratio (SNR), cloud contamination, or technical problems), we have considered hd,AEMET also valid for the intercomparison provided that they lie between two IGN observations, since hd,IGN is not expected to change appreciably in the interval between two consecutive IGN measurements (see Sect. ).
Detailed information on each system that is part of the network will be presented below, including the different methodologies used in this monitoring network. The approaches used in this study, that is, the Gradient and the Wavelet Covariance Transform methods, were applied to all AEMET-ACTRIS profilers except the RMO station. These methodologies can be considered equivalent and reliable according to .
The ARCADE Raman Lidar at Roque de los Muchachos
This was the only existing profiler in La Palma at the time of the eruption, providing the only information about plume height until the next profiler at Fuencaliente was installed, 11 d after the eruption began. ARCADE is a Raman Lidar system (Fig. a) installed at the Cherenkov Telescope Array Observatory (CTAO) at RMO in October 2018, with the aim of providing nightly and seasonal aerosol attenuation profiles of UV light measured for the CTA two times a day, at sunrise and sunset, in automatic and unattended mode. Since the production of Cherenkov light depends on the molecular profile, aerosol profiles are required to estimate the attenuation and propagation in CTA observations . ARCADE is a joint project between CTA members from the National Institute for Nuclear Physics (INFN) and the University of L'Aquila.
View of the different instruments belonging to the AEMET-ACTRIS profiling network: (a) ARCADE Raman lidar at Roque de los Muchachos, (b) MPL-4B at Tazacorte, (c) CHM15k at Fuencaliente, (d) CL61 at El Paso and (e) CL51 at La Palma Airport.
This system is designed to measure elastic and Raman backscattered photons using a 20 cm diameter Newtonian telescope and a Nd:YAG laser source operating at 355 nm. The laser delivers pulses with a maximum energy of 5–6 mJ and a variable repetition rate ranging from 1 to 100 Hz . A second Raman channel for water vapour was added to retrieve atmospheric water vapour profiles; however, due to the COVID-19 outbreak, the Raman channels have been offline since 2021. The system was dismantled in October 2023 for future deployment at CTAO-South.
hd,AEMET was qualitatively estimated with this instrument from 19 September to 27 October at two fixed times per day (06:00 and 19:00 UTC), resulting in a total of 18 hd,AEMET values retrieved from this site. The estimation was based on visual inspection of the aerosol backscatter profiles, when available. Although the approach used at RMO is qualitative, the plume altitude measurements provided by this instrument have been included in this study, as they represent the only available information on the volcanic plume during the initial phase of the eruption and were the sole profiler data from the northern part of La Palma until the CL61 at El Paso station became operational on 25 October. This information proved to be extremely valuable during the eruptive process, with some of the ARCADE-derived altitudes reportedly used by PEVOLCA. From 27 October 2021 onwards, a technical issue prevented the ARCADE team from continuing plume monitoring, and the instrument's operation was definitively discontinued.
MPL-4B micropulse lidar at Tazacorte
The MPL-4B is a micropulse lidar installed at Tazacorte station on 15 October 2021 (Fig. b), 25 d after the eruption started. This instrument is jointly owned by the Universitat Politècnica de Catalunya (UPC) and the Instituto Nacional de Técnica Aeroespacial (INTA). MPL operates continuously (24/7) with a relatively high frequency (2500 Hz) and low-energy (∼ 7 µJ) using a Nd:YLF laser centred at 532 nm. As a part of NASA/MPLNET (https://mplnet.gsfc.nasa.gov, last access: 15 March 2025), this system provides vertical profiles of clouds and aerosols with a 1 min temporal resolution and 75 m of vertical resolution. MPL-4B also includes polarisation capabilities, which are particularly useful to distinguish between spherical and non-spherical particles .
The methodology used to retrieve hd,AEMET from this instrument is based on the Gradient Method , as described in , , and . In this case, the Range-Square-Corrected lidar Signal (RSCS) was integrated over a ±5 min interval around the IGN observation time in order to retrieve the height corresponding to the peak of the uppermost volcanic plume. A reference RSCS value at 8 km was used, and a simple threshold method was applied to the RSCS profiles to detect cases of signal blocking or strong attenuation. Further details on the algorithms employed in this study can be found in and .
A total of 24 altitudes retrieved with this MPL-4B instrument constitute the final hd,AEMET data series.
Lufft CHM15k ceilometer at Fuencaliente
This Lufft CHM15k ceilometer , belonging to the Group of Atmospheric Optics of the University of Valladolid (GOA-UVa), is part of the Iberian Ceilometer Network . The CHM15k operates continuously (24/7) with a pulse energy of 8.4 µJ and a repetition frequency of 5–7 kHz (Fig. c). The pulsed Nd:YAG laser emits at 1064 nm, and the backscattered signal with a vertical resolution of 15 m is collected using a telescope with a field of view of 0.45 mrad. The use of near-infrared (nIR) laser emission minimises uncertainty due to avoiding water vapour absorption, compared to ceilometers operating at lower wavelengths. According to , the full overlap of the system is reached at approximately 1500 m a.s.l. Other studies report that the maximum detection height of this system is 15 360 m a.s.l., and that the overlap-corrected signal begins at 80 m a.s.l. In this study, overlap has been corrected according to the overlap function provided by the manufacturer. More technical details about the CHM15k can be found in .
This ceilometer, installed at the Fuencaliente station on 30 September, provided RSCS profiles to this study from all products available in the CHM15k Nimbus internal software. These profiles were integrated over a ±5 min interval around the IGN observation time. As in the MPL procedure, 8 km was considered the maximum altitude with significant aerosol influence. The Gradient Method was used to detect the uppermost aerosol layer based on sharp negative gradients in the backscattered signal, which are sensitive to aerosol concentration. An empirical SNR threshold of 3 was applied to limit the maximum height at which data are considered reliable .
Using this methodology, a total of 63 altitudes were included in the final hd,AEMET dataset.
Vaisala CL51 ceilometer at La Palma Airport
The CL51 ceilometer (Fig. d), belonging to AEMET, operates continuously (24/7), using a pulsed laser emitter centred at 910 nm to measure backscattered radiation via a telescope in a coaxial configuration with an avalanche photodiode detector. This configuration allows the system to achieve an overlap height as low as 50 m, according to the manufacturer . The spatial and temporal resolutions of the CL51 are 10 m and 15 s, respectively. The capability of this instrument for aerosol retrieval has been validated in several studies and references therein, and was specifically applied to the La Palma volcanic eruption in . The methodology used to retrieve the aerosol layer top is described in the latter study and is based on the Continuous Wavelet Covariance Transform (WCT) method .
This instrument was installed at La Palma Airport on 7 October and contributed a total of 26 altitudes to the final hd,AEMET dataset.
Vaisala CL61 lidar ceilometer at El Paso
The Vaisala CL61 ceilometer (Fig. e) is a novel, high-performance instrument capable of depolarisation measurements, designed for unattended 24/7 operation under all weather conditions. It uses a single avalanche photodiode detector along a coaxial optical path. The CL61 operates at a wavelength of 910.5 nm, with a vertical resolution of 4.8 m, a temporal resolution of up to 5 s, and a maximum range of up to 15 400 m. Its patented optical design enables the use of a narrow spectral bandwidth, which minimises water vapour absorption and ensures wavelength stability across varying temperatures.
The transmitter operates at a pulsing frequency of 9.5 kHz, with an average laser power of 40 mW and a pulse length of 160 ns (Full Width at Half Maximum, FWHM). The CL61 also incorporates an enhanced single-lens optical system with a coaxial configuration for both transmitted and received signals, featuring individual overlap functions calibrated for each instrument. This optimised overlap function enables the reliable detection of atmospheric layers at low altitudes.
The attenuated backscatter profiles are pre-calibrated at the Vaisala facilities using liquid water clouds, following the methodology described in . Since the CL61 uses the same receiver module for both cross-polarised (x-pol) and parallel-polarised (p-pol) signals, no additional receiver sensitivity calibration is required.
A similar methodology to that described in the previous section for the CL51 was applied to the CL61 to retrieve attenuated backscatter profiles, following .
Six altitude levels retrieved by the CL61 at the El Paso station were included in the final hd,AEMET dataset, contributing low statistical significance to the data from this station.
Seismic records and ground deformation analyses are critically important for monitoring eruptive events . The Real-time Seismic Amplitude Measurement (RSAM), based on the vertical component of broadband seismic data, is widely recognised as a reliable indicator of magma transport from depth to the surface vent .
In this study, RSAM data were extracted from the IGN permanent seismic network at the CENR station (28.63° N; 17.85° W; 1208 m a.s.l.), the closest station to the Tajogaite eruption site. This station is equipped with a three-component broadband seismic sensor and has been operating continuously since 2017. Following , RSAM was calculated using dedicated software (ThomasLecocq/ssxm: RSAM/RSEM – SSAM/SSEM easy code), applying 10 min time windows. The resulting time series was normalised by subtracting the median and dividing by the standard deviation.
Meteorological information: Height of the Trade Wind Inversion (TWI) and wind speed/direction
Radiosonde vertical profiles are launched daily at 00:00 and 12:00 UTC from AEMET’s station at Güímar (28.32° N, 16.38° W; 105 m a.s.l.). This site is a World Meteorological Organization (WMO) Global Climate Observing System (GCOS) Upper-Air Network (GUAN) station (no. 60018). This GCOS-GUAN station is located on the eastern coast of Tenerife, leeward of the prevailing trade wind flow due to the presence of a central ridge in the Güímar Valley, which reaches elevations above 2100 m a.s.l. This station has been providing upper-air profiles since 2003, measuring temperature, pressure, and humidity using Vaisala RS92 radiosondes, along with wind speed and direction obtained via a Global Positioning System (GPS) wind-finding system .
Atmospheric profiling was also conducted on La Palma (28.63° N, 17.91° W, 295 m a.s.l.) during the volcanic eruption, from 6 November to 26 December 2021 at 11:00 UTC, with the collaboration of AEMET and Unidad Militar de Emergencias (UME). At this temporary station, the same Vaisala RS92 radiosondes were deployed to characterise the atmosphere near the volcano. This site shares similar characteristics with the Güímar station, as it is also located on the eastern coast of La Palma, leeward of the prevailing trade winds and downstream of the island's central ridge (Cumbre Vieja), which reaches a maximum elevation of 1949 m a.s.l.
The Trade Wind Inversion height was estimated following the methodology described in , based on the analysis of the temperature lapse rate, water vapour mixing ratio, and wind components on both AEMET stations (at Tenerife and La Palma).
In addition, daily wind speed and direction profiles from the AEMET numerical weather prediction (NWP) model HARMONIE-AROME (HIRLAM–ALADIN Regional/Meso-scale Operational NWP in Euromed) Version 43h2.1.1 were used. Data were extracted from the grid point closest to the eruptive centre. This model has a spatial resolution of 2.5 × 2.5 km2, a temporal resolution of 3 h (with forecasts up to 72 h), and includes 65 vertical levels. Further details about the HARMONIE-AROME model can be found in .
Satellite derived SO<sub>2</sub> and Aerosol Layer Height (ALH)
TROPOMI (TROPOspheric Monitoring Instrument) is a sensor onboard the polar, low-Earth orbiting Copernicus Sentinel-5 Precursor (S-5P) platform, capable of providing SO2 altitude layer information as a near-real-time operational product. This product is generated by the German Aerospace Center (DLR, Deutsches Zentrum für Luft- und Raumfahrt) using the Full-Physics Inverse Learning Machine (FP-ILM) algorithm. The SO2 layer height (SO2LH) product was developed by with an expected accuracy of less than 2 km for SO2 vertical column densities (VCD) greater than 20 DU (Dobson Units) .
The hd,AEMET dataset was used by to validate TROPOMI SO2LH during the La Palma volcanic eruption. Altitude differences, based on median TROPOMI SO2LH, were found to range between 1 and 3 km, decreasing to around 1 km for cloud fractions below 0.5. TROPOMI LH slightly underestimates the hd,AEMET measurements, but the two datasets are correlated, with a Pearson correlation coefficient of r= 0.74 .
Another lidar sensor used in this context to provide information on the aerosol LH (ALH) is CALIOP. CALIOP was part of the CALIPSO mission providing global atmospheric profiles from June 2006 to August 2023 . This instrument orbited the Earth in a sun-synchronous orbit, acquiring lidar backscatter profiles at 532 and 1064 nm, including polarisation capabilities at 532 nm . Level 2 version 4.2 aerosol layer information product (https://subset.larc.nasa.gov/calipso/, last access: 15 March 2025) with a 5 km of spatial resolution was used to provide information on the aerosol layer top altitude. The CALIOP layer detection algorithm is described in detail in .
MSVOLSO2L4 is a multi-satellite, Level 4 version 4 volcanic SO2 long-term global database product that estimates the total SO2 mass (in metric tons) within a volcanic cloud by combining SO2 column data with information on plume altitude . The latter is a critical parameter, typically obtained from direct observations (e.g., pilot reports or ground-based measurements). This independent database contains all significant volcanic eruptions detected from space, offering daily estimates of volcanic SO2 emissions. Since 2018, the product has relied on three spaceborne UV sensors to retrieve SO2 VCD: the Ozone Monitoring Instrument (OMI) , the Ozone Mapping and Profiler Suite (OMPS) , and the TROPOspheric Monitoring Instrument (TROPOMI) . Infrared (IR) satellite observations have also been used to detect SO2 emissions not captured by UV sensors, thereby refining the UV-based product . These IR sensors include the Atmospheric Infrared Sounder (AIRS) , the Moderate Resolution Imaging Spectroradiometer (MODIS) , and the Infrared Atmospheric Sounding Interferometer (IASI) , which are used for both refinement and validation purposes.
In the absence of such observations, a default plume altitude is assumed based on the eruption style (effusive or explosive) and magnitude, inferred from the Volcanic Explosivity Index (VEI). These default altitudes can range from approximately 5 km for effusive events to 10 km for explosive eruptions.
For the initial stages of the La Palma eruption, the NASA team assumed a default injection height of 8 km, which was used in the default SO2 mass loading estimates. Once a profiling network was established on La Palma, a revised product was generated for this study using the observed plume-top height (hec), derived on-site using the IGN video-surveillance method.
A comparative analysis between the default 8 km-based retrievals and the adjusted SO2 emissions, based on ground-based plume height measurements, provides valuable insight into the impact of plume altitude assumptions on the derived SO2 emissions.
ResultsComparison of volcanic plume height from IGN and AEMET-ACTRIS approaches
The eruptive column height measured by the IGN (hec,IGN) is shown in Fig. , along with the dispersive plume height (hd) obtained from both the IGN and AEMET-ACTRIS networks (hd,IGN and hd,AEMET, respectively). hec,IGN and hd,IGN correspond to 344 measurements made by the IGN at times when a significant change in column height was observed. hd,AEMET (137 altitudes) were retrieved concurrently with the IGN observations using the methodology described in Sect. for each of the five profiling stations in the network (integrated over a ±5 min interval around the IGN observation time).
Altitude of the eruptive column in m a.s.l. measured by IGN (hec,IGN) and dispersion plume (hd) measured by IGN (hd,IGN) and by AEMET-ACTRIS profiling network (hd,AEMET). Broken lines represent the corresponding lowess smoothing curves for hd,IGN and hd,AEMET datasets.
Figure and Table shows a good agreement between the IGN and AEMET-ACTRIS datasets. The average absolute difference between hd,IGN and hd,AEMET throughout the eruptive period was 258.6 m (standard deviation of 620.4 m). This high consistency (Person coefficient, r, of 0.81) is noteworthy considering the different methodologies involved: one based on visual inspection from video-surveillance and the other one from a multi-instrumental approach including different profiling systems. These differences were further analysed by station to assess their spatial dependence. The highest absolute difference was observed at ELP (941.7 m), while the lowest was recorded at RMO (139.8 m). However, the minimum r value was observed at this station (0.18), while FUE and TAZ exhibited the maximum r values (0.86 and 0.89, respectively). In terms of relative differences, the mean deviation between the two datasets was 12.9 %. The maximum relative difference was again found at ELP (42.7 %), with the lowest at RMO (2.6 %). Intermediate relative differences were observed at FUE (8.1 %), TAZ (19.8 %), and LPA (22.3 %). The larger differences and the marked worsening of the skill scores observed at ELP may be attributed to its proximity to the source (potential signal attenuation and significant impact of sudden changes in volcanic activity) or to complex topographic effects affecting the site. However, the limited number of measurements at this site (12) prevents us from establishing any robust error metrics for comparison with other instruments in the database. In light of the comparison results shown in Table , it can be observed that the methodology used to derive the volcanic plume height from lidar data does not appear to play a dominant role in the comparison outcomes, with the best statistics obtained for the FUE and TAZ cases (Gradient Method in both cases).
Main statistical skill scores (in m) for the comparison between hd,IGN and hd,AEMET differences, including also the multi-instrument and inter-method comparison between hd,IGN and the height of the volcanic plume measured by AEMET at the five different stations. The methodology for retrieving the altitude of the volcanic plume is also included. GM stands for the Gradient Method, and WCT for the Wavelet Covariance Transform method.
An additional insight from Fig. is the relatively small difference between hec and hd, both measured by IGN, with an average difference of only 236.2 m. Larger differences were found in the early stages of the eruption and during the final paroxysmal event on 13 December , likely reflecting more vigorous eruptive dynamics and increased column buoyancy during those phases. As a preliminary conclusion, it can be stated that both techniques provide comparable measurements and can be combined into a single data series.
Description of the eruptive event in terms of tracers of the volcanic plumes (<inline-formula><mml:math id="M118" display="inline"><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi mathvariant="normal">d</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) and modulating factors (RSAM and TWI)
The height reached by the eruptive column and its subsequent dispersion, as explained in Sect. , depends on both volcanological and meteorological parameters, as well as their interaction. According to , key factors influencing the rise of an eruptive volcanic column include the total grain size distribution of emitted volcanic particles, the amount of gas released during magma fragmentation, and atmospheric crosswinds. Other studies e.g., also highlight the importance of parameters such as vent radius, plume temperature at the vent, upward velocity, and the wind entrainment coefficient.
An important modulating factor of the altitude of the eruptive column is the eruptive dynamics, which can be parameterised using proxies such as SO2 emissions or seismic activity . The evolution of the eruptive process during the Tajogaite eruption is shown in Fig. b, represented by the RSAM signal.
(a) Altitude of the eruptive column in m a.s.l. measured by IGN (hd,IGN) and by AEMET-ACTRIS profiling network (hd,AEMET). Yellow bands indicate the presence of the SAL. (b) Time series of the raw real-time seismic amplitude measurement (RSAM) measured in Cumbre Vieja by the IGN seismic network. (c) Trade Wind Inversion (TWI) height (m a.s.l.) estimated at 00:00 and 12:00 UTC from meteorological sondes launched at Güímar (Tenerife) station and at 12:00 UTC from La Palma station during the volcanic eruption (from 6 November 2021 onwards).
Following , another key factor that modulates the injection height of volcanic emissions is the prevailing meteorological conditions. Wind speed and direction are fundamental variables that influence the direct transport and dispersion of the volcanic plume. In addition, the presence of thermal inversions has a critical effect on the vertical extent, mixing, and dispersal of volcanic emissions. As demonstrated by , a shallower temperature inversion (TWI in our case) confines the volcanic plume to lower levels of the troposphere, with significant implications for air quality and civil aviation safety.
Yellow bands in Fig. indicate days affected by mineral dust conditions over La Palma, as defined by (26 September–3 October; 7–8 October; 19–21 October). The presence of the Saharan Air Layer (SAL) can influence the altitude of the aerosol layer as a result of its effect on vertical stratification and atmospheric stability. According to , the SAL appears as a well-stratified layer capped by a temperature inversion located between ∼ 6.3 km in summer (when the dusty air moves westwards over the subtropical North Atlantic as an elevated layer over the Marine Boundary Layer, MBL). The SAL, therefore, compresses the MBL when transported as an elevated layer.
A comprehensive view of the eruptive process is presented in Fig. , which combines information on plume altitude with the evolution of the two modulating factors previously described: eruptive dynamics and meteorological conditions. This multidisciplinary approach allows for a detailed interpretation of the eruption's progression. Based on the trend in the RSAM signal (Fig. b), three distinct eruptive phases can be identified:
Phase I (19–27 September) is marked by alternating Strombolian and Vulcanian activity, including two strong explosive episodes on 23 and 26 September , which led to the highest RSAM value recorded during the entire eruption (4.4 a.u.). This period coincides with peak values of both eruptive column height (hec) and dispersive plume altitude (hd), reaching up to 6500 m a.s.l. and 5500 m a.s.l., respectively. Simultaneously, a downward trend in the TWI height was observed (Fig. c), with values decreasing from 2074.7 to 601.4 m a.s.l. by the end of the phase. Mean TWI heights during this period were 1576.9 m a.s.l. (12:00 UTC) and 1498.4 m a.s.l. (00:00 UTC), based on radiosonde data from Güímar. These conditions suggest a limited surface-level impact of volcanic emissions, as corroborated by surface SO2 concentration and particulate matter measurements reported in .
Phase II (28 September–2 November) is characterised by predominantly effusive activity with intermittent explosive events. The RSAM signal in this phase remains mostly between 0 and 1 a.u., with a notable increase from 24 October onwards, reaching up to 3.1 a.u. . This intensification coincides with the second-largest increase in coarse-mode aerosol fractions during the eruption (25 October–2 November), detected via ground-based photometry . The TWI shows significant variability, peaking at 2780 m a.s.l. on 25 October and dropping to minimum values of 263 m a.s.l. (4 October) and 220–365 m a.s.l. (16–20 October). Average TWI heights were 1039.0 and 955.9 m a.s.l. at 12:00 and 00:00 UTC, respectively. During this phase, plume heights remained relatively stable, with mean hd values of 2935.6 m a.s.l. (AEMET-ACTRIS) and 3374.4 m a.s.l. (IGN). The maximum hec and hd values (5300 and 4800 m a.s.l., respectively) occurred between 17 and 18 October.
Phase III (3 November–13 December) was dominated by low-intensity effusive activity, reflected in RSAM values generally below 0.5 a.u., with a few sharp peaks above 2 a.u. These peaks correspond to discrete events, including large ash emissions on 17 November, the opening of a new vent on 28 November, intense explosive activity with shockwaves on 1–2 December, and the final paroxysmal phase with Vulcanian episodes on 13 December . Despite the overall lower eruptive energy, this phase exhibited a stable aerosol layer and notable increases in plume heights, particularly on 12–13 December, when hec reached 8500 m a.s.l. TWI values during this period were consistently high, with mean heights of 1712.5 m a.s.l. (12:00 UTC), 1695.1 m a.s.l. (00:00 UTC) at Güímar, and 1819.3 m a.s.l. at La Palma (12:00 UTC). The intense seismic activity on 1–2 December also elevated both hec and hd to approximately 5000 m a.s.l. Observations by confirmed high aerosol optical depth (AOD) coarse-mode values (0.15) during this time.
A key inflexion point between phases II and III is observed as a reduction in RSAM, also identified by and , the latter through daily SO2 emissions from TROPOMI. However, reported that surface-level SO2 concentrations did not decrease proportionally to emission rates or seismic signals. This discrepancy is attributed to the maximum TWI extension during this phase, as seen in Fig. c, with mean heights ranging from 1695.1 to 1819.3 m a.s.l. This factor, combined with the low volcanic intensity, is expected to confine the volcanic plume within the lower layers of the troposphere. Consequently, this final phase of the eruption had the greatest impact on surface air quality despite exhibiting lower overall volcanic activity.
hd,IGN and hd,AEMET mean differences for the 35 d affected by the SAL are reduced to 33.4 m (1.1 %), highlighting the influence of the SAL on vertical atmospheric stratification and the role of the strong temperature inversion at the top of the SAL. This stable layer acts as a physical barrier, hindering the vertical dispersion of the volcanic plume and confining aerosols within the Planetary Boundary Layer (PBL). This situation results from the balance between the upward velocity of the volcanic plume and the strength of the inversion layer.
Another important aspect shown in Fig. a is the spatial representativity of the hd,AEMET measurements, which are derived from the station most affected by the volcanic plume depending on prevailing wind direction. The spatial distribution of these data reveals a significant bias toward the western side of La Palma Island, as expected, given the predominant trade wind regime (more than 63.5 % of the data in the hd,AEMET time series were recorded at the Tazacorte and Fuencaliente stations). Specifically, the number of plume height measurements per station is as follows: 18 from ORM (mainly restricted to Phase I), 24 from TAZ, 63 from FUE, 6 from ELP, and 26 from LPA. This distribution highlights the role of wind direction in plume transport and vertical dispersion, further evidenced in Fig. , which displays the wind rose diagram for the study period. The trade wind regime (originating from the NE) remains predominant while the volcanic plume height stays within the marine boundary layer (up to 2000 m a.s.l.). Above this level, a change in prevailing wind direction occurs, as expected from the vertical balance of forces at this latitude. Up to 4000 m a.s.l., two predominant components are observed (NE and W–WSW). It must be noted that this analysis reflects the wind structure under the specific meteorological conditions of the 85 d eruption period and cannot be interpreted as a climatological pattern.
Wind rose diagram for the volcanic eruptive column hec (in m a.s.l.) and HARMONIE-AROME wind vectors at the pixel of the volcanic edifice in Cumbre Vieja. Contour lines are represented with the frequency of occurrence in %.
Comparison between <inline-formula><mml:math id="M139" display="inline"><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mrow><mml:mi mathvariant="normal">d</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="normal">AEMET</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> with CALIOP aerosol altitude (ALH<sub>CALIOP</sub>)
A total of 12 CALIOP overpasses during the volcanic eruption were used in this study in the evaluation of the CALIOP level 2 layer information product. The different altitudes (ALHCALIOP) corresponding to the highest aerosol layer detected by CALIOP are shown in Fig. together with the hd,AEMET series. A mean difference of 615.0 m was retrieved with a clear underestimation of the CALIOP altitudes with respect to the ground-based measurements. These overpasses were done at an average distance of 61.7 km from the volcanic vent, and the AEMET-ACTRIS versus CALIOP differences do not seem to be dependent on the overpass distance. The maximum difference was observed on 6 October (a discrepancy of 3065.5 m, with an overpass distance of 28 km). The aerosol layer altitude difference on this particular day is attributed to the presence of a thin aerosol layer located between approximately 2 and 4 km, as observed by the CHM15k ceilometer at FUE, with the PBL situated below (not shown for the sake of brevity). It appears that CALIOP was unable to detect weak backscatter signals, presumably because they were below its detection threshold. Excluding this day, the mean difference is reduced to 392.2 m. These differences are consistent with previous studies using ground-based lidar measurements , as well as with results obtained using other techniques .
Altitude of the dispersion plume hd,AEMET (in m) and aerosol layer height determined from CALIOP (ALHCALIOP) measured during the entire volcanic eruption, with the distance of the overpass from the volcanic vent (in km) in the colorbar.
Impact of the altitude of the eruptive plume on satellite-based SO<sub>2</sub> emission estimation
As already stated in the first part of this paper, the height of the volcanic plume is an important factor to take into account when it comes to estimating the SO2 emission loadings using satellite observations. Figure shows the evolution of the eruptive event in terms of daily SO2 emission (in kt d−1) assuming a reference hec of 8 km (standard MSVOLSO2L4 NASA product) for effusive volcanoes such as Tajojaite . In this figure, SO2 estimations using hec,IGN observations are also included. The change in the eruptive activity is apparent from this figure, showing a strong decrease in the measured volcanic SO2 emissions during the first days of November. Average SO2 during the first phase was 23.8 kt d−1 (42.9 kt d−1 in the case of the corrected series using the hec,IGN), in the second phase was 17.3 (44.3) kt d−1, while this average in the case of the last phase is reduced to 4.7 (10.1) kt d−1. This decrease is attributed to the change in the eruptive activity observed from the RSAM series in Sect. . Consistent daily emissions for the Tajogaite volcano were calculated by by means of a novel forward trajectory approach using PlumeTraj, a pixel-based trajectory analysis of an SO2 cloud . Ground-based miniDOAS observations of SO2 emissions were also carried out during the Tajogaite eruption. This is exemplified by the studies of , which reported SO2 emissions ranging from 670 to 17 t d−1 and observed a clear decreasing trend in SO2 during the post-eruptive phase.
NASA daily SO2 emission (in kilotons per day, kt d−1) calculated as the average of OMI, OMPS and TROPOMI UV backscatter radiances considering a standard columnar injection height of 8 km (blue circles) and the real hec measured by the IGN (in red triangles). Dotted vertical lines represent the three eruptive phases derived from RSAM data series.
Regarding the impact of using a real hec instead of an a priori value (8 km in the case of La Palma eruption), we observe significant underestimation of the SO2 emission mass estimated from the satellite. An average SO2 mass emission difference of 17.3 kt d−1 has been found, with maximum values up to 64.2 kt d−1 (found on 22 October). Average relative differences of 56.2 % have been found in the whole eruptive process, with maximum values found on 13 October (84.7 %). According to , if the actual plume altitude is lower than the assumed value – as is the case in this study – the retrieved SO2 column amount and the SO2 emission mass will be underestimated, as it is observed here.
Comparable daily SO2 emission values to those retrieved using hec,IGN were reported by for the same eruption, based on the PlumeTraj analysis, with lower mean relative differences of 21.6 %.
Summary and conclusions
This work provides a comprehensive description of the eruptive process that took place on La Palma Island (19 September–December 2021) from a different scientific perspective, offering complementary information and novel insights to those of previous studies. Regarded as the most significant eruption in Europe over the past 75 years in terms of SO2 emissions and societal impact, the Tajogaite eruption provided a unique natural laboratory to study explosive volcanic activity from multiple scientific perspectives. In our case, this was made possible through an unprecedented collaborative effort among public institutions, research groups, and private entities, coordinated by IGN and AEMET, the latter within the framework of ACTRIS and ACTRIS-Spain .
The characterisation of the eruption presented in this study focused on key plume tracers, including the height of the eruptive column (hec) and the dispersive volcanic plume (hd), as well as modulating factors such as seismic activity and meteorological conditions. Based on the temporal evolution of the RSAM signal, three distinct eruptive phases were identified. These phases reflect changes in the eruptive dynamics, plume altitude, and associated seismic and atmospheric conditions throughout the event. In this regard, our results reveal a complex eruptive process characterised by a range of eruptive styles, from explosive Strombolian phases to predominantly effusive activity.
Two independent and complementary observational techniques were employed to monitor the vertical evolution of the plume: the video-surveillance system operated by IGN and the profiling network deployed by AEMET-ACTRIS, which includes different profiling systems positioned strategically around the volcanic vent. This key information was incorporated into the PEVOLCA reports (in the case of both IGN and AEMET-ACTRIS datasets) during the whole volcanic crisis. The VONA alerts and the regular reports submitted to the Toulouse VAAC were informed by the IGN dataset.
These two techniques present intrinsic limitations. The first method, based on visual observations from the IGN camera network, is generally not applicable during nighttime hours. Nevertheless, during the most intense eruptive phases, nighttime observations were feasible due to the high luminosity of the eruptive column. Additionally, under conditions of dense cloud cover or fog, the top of the plume was often obscured, making accurate estimations difficult. The second method, based on vertical profiling (AEMET-ACTRIS network), only measures the atmosphere directly above each station. Consequently, plume detection is limited to instances when the volcanic layer is present in the station's vertical column. Moreover, the AEMET-ACTRIS dataset integrates data from multiple instruments and locations, which are differently influenced by the prevailing East-to-West trade wind pattern – resulting in a spatial sampling pattern, with over 63.5 % of data originating from Tazacorte and Fuencaliente on the western flank.
The relatively low overall discrepancies found in this study (258.6 m) confirm the coherence between the two retrieval techniques. This agreement is noteworthy given the different methodologies involved, and it demonstrates that both techniques can yield comparable volcanic plume height measurements, which may be combined into a single data series. This approach may be of use in future volcanic crises and could support operational surveillance during such events.
The hd,AEMET dataset was also used to evaluate satellite-derived plume altitudes near the source. Comparisons with the CALIOP aerosol layer height (ALHCALIOP) showed a mean underestimation of 615.0 m. This value can be reduced to 392.2 m if one specific day, characterised by the presence of a thin aerosol layer, is excluded from the comparison. In addition to their inherent operational use in volcanic monitoring, we have demonstrated in this work that the hd and hec products also provide an opportunity to evaluate different satellite-derived datasets, which is crucial for improving their applicability in future volcanic crises.
The impact of assuming a fixed plume height in the estimation of satellite-based SO2 emissions was also assessed. Two versions of NASA’s MSVOLSO2L4 product were analysed: (1) an initial estimation assuming a fixed plume height of 8 km above the vent, and (2) a final version using actual hec values provided by IGN. Our results indicate that assuming a fixed plume height leads to significant underestimations of SO2 mass loadings, with average discrepancies of 56.2 % over the eruption and peak deviations reaching 84.7 %. These findings underscore the critical importance of accurate plume height estimation when quantifying volcanic emissions using satellite data.
Data availability
Data from MPLNET used in the present study can be obtained from https://mplnet.gsfc.nasa.gov/download_tool/ (last access: 21 February 2025). The vertical meteorological soundings can be downloaded from http://weather.uwyo.edu/upperair/sounding.html (last access: 14 March 2025). MSVOLSO2L4 is accessible at https://disc.gsfc.nasa.gov/datasets/MSVOLSO2L4_4/summary (last access: 15 April 2025). Level 2 version 4.2 CALIOP aerosol layer information product is available at https://subset.larc.nasa.gov/calipso/(last access: 14 March 2025). All the information about the height of the volcanic plume will be available upon request to the contributing authors.
Author contributions
AB and OG provided the volcanic plume height from the AEMET-ACTRIS network. FQ and JP-P provided the height from the IGN network. DG did the data analysis concerning CALIOP and wrote and revised the paper. AB-V provided the volcanic height from CL51 and CL61 profilers; MS and CC-J provided the altitudes from MPL-4B; MI and VR provided the height of the volcanic plume from ARCADE; and NK and SC provided the SO2 mass emission estimates from satellites. The remaining authors, listed in alphabetical order, actively contributed scientifically to the development of the paper, either through fieldwork during the eruption or by supporting the retrieval of the volcanic plume height series throughout the three-month eruptive period. Furthermore, EW provided calibrated data of the MPL lidar; AF, SM, SC and NK also contributed to the writing and the discussion on the data. PH and DL provided insight on TROPOMI data. MH was in charge of vertical soundings at La Palma during the eruption. All authors were involved in helpful discussions and contributed to the manuscript.
Competing interests
At least one of the (co-)authors is a member of the editorial board of Atmospheric Measurement Techniques. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.
Disclaimer
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.
Special issue statement
This article is part of the special issue “Sun-photometric measurements of aerosols: harmonization, comparisons, synergies, effects, and applications”. It is not associated with a conference.
Acknowledgements
This work is part of the activities of the WMO-Measurement Lead Centre for aerosols and water vapour remote sensing instruments (MLC). We gratefully acknowledge the data provided by the MPLNet network. The MPLNET project is funded by the NASA Radiation Sciences Program and the Earth Observing System. The authors also wish to acknowledge the ACTRIS (Aerosol, Clouds, and Trace Gases Research Infrastructure) for its essential role in providing access to high-quality data, facilities, and technical support that were integral to the execution of this research. The authors gratefully acknowledge the extraordinary effort carried out by the AEMET staff (both in La Palma and in support of the activities in La Palma) during the volcanic eruption, from the Izaña Atmospheric Research Center and the Delegation of AEMET in the Canary Islands. We also gratefully acknowledge the dedication of Dr. David Suárez and all the information provided by the PEVOLCA Scientific Committee, in addition to the support received from the insular and local governments (Cabildo Insular de La Palma and the Ayuntamientos de Tazacorte, Los Llanos de Aridane and El Paso). Finally, the authors would like to thank Dr. Emilio Cuevas for his support and leadership in this project, making the necessary institutional collaboration possible to carry out the measurements presented in this paper.
Financial support
Scientific activities in this paper were done under the ACTRIS grant (agreement no. 871115), ACTRIS ERIC Spain Grant RED2024-153891-E funded by MICIU/AEI/10.13039/501100011033 and other grants supported by the Ministerio de Ciencia e Innovación (MICINN) (PID2023-151666NB-I00 and PID2021-127588OB-I00). Financial support of the Department of Education, Junta de Castilla y Leon, and FEDER Funds is gratefully acknowledged (Reference: CLU-2023-1-05). Michael Sicard is supported by the Horizon Europe European Research Council (project REALISTIC, grant no. 101086690). This research has also been supported by COST (European Cooperation in Science and Technology) under the HARMONIA (International network for harmonization of atmospheric aerosol retrievals from ground-based photometers) Action CA21119.
Review statement
This paper was edited by Lionel Doppler and reviewed by two anonymous referees.
ReferencesACTRIS 2021: ACTRIS-Spain Coordinating Unprecedented Actions for the Cumbre Vieja Volcanic Emergency | ACTRIS, https://www.actris.eu/news-events/news/actris-spain-coordina-ting-unprecedented-actions-cumbre-vieja-volcanic-emergency, last access: 15 March 2025.Albertos, V. T., Recio, G., Alonso, M., Amonte, C., Rodríguez, F., Rodríguez, C., Pitti, L., Leal, V., Cervigón, G., González, J., Przeor, M., Santana-León, J. M., Barrancos, J., Hernández, P. A., Padilla, G. D., Melián, G. V., Padrón, E., Asensio-Ramos, M., and Pérez, N. M.: Sulphur dioxide (SO2) emissions by means of miniDOAS measurements during the 2021 eruption of Cumbre Vieja volcano, La Palma, Canary Islands, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-5603, 10.5194/egusphere-egu22-5603, 2022.Ansmann, A., Tesche, M., Groß, S., Freudenthaler, V., Seifert, P., Hiebsch, A., Schmidt, J., Wandinger, U., Mattis, I., Müller, 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, 10.1029/2010GL043809, 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 fine-mode particle mass profiles from EARLINET-AERONET observations over central Europe after the eruptions of the Eyjafjallajökull volcano in 2010, J. Geophys. Res.-Atmos., 116, D00U02, 10.1029/2010JD015567, 2011.Ansmann, A., Seifert, P., Tesche, M., and Wandinger, U.: Profiling of fine and coarse particle mass: case studies of Saharan dust and Eyjafjallajökull/Grimsvötn volcanic plumes, Atmos. Chem. Phys., 12, 9399–9415, 10.5194/acp-12-9399-2012, 2012.Aubry, T. J., Staunton-Sykes, J., Marshall, L. R., Haywood, J., Abraham, N. L., and Schmidt, A.: Climate change modulates the stratospheric volcanic sulfate aerosol lifecycle and radiative forcing from tropical eruptions, Nat. Commun., 12, 4708, 10.1038/s41467-021-24943-7, 2021.Baars, H., Ansmann, A., Engelmann, R., and Althausen, D.: Continuous monitoring of the boundary-layer top with lidar, Atmos. Chem. Phys., 8, 7281–7296, 10.5194/acp-8-7281-2008, 2008.Barreto, Á., Cuevas, E., García, R. D., Carrillo, J., Prospero, J. M., Ilić, L., Basart, S., Berjón, A. J., Marrero, C. L., Hernández, Y., Bustos, J. J., Ničković, S., and Yela, M.: Long-term characterisation of the vertical structure of the Saharan Air Layer over the Canary Islands using lidar and radiosonde profiles: implications for radiative and cloud processes over the subtropical Atlantic Ocean, Atmos. Chem. Phys., 22, 739–763, 10.5194/acp-22-739-2022, 2022.Bartolini, S., López, C., Becerril, L., Sobradelo, R., and Martí, J.: A retrospective study of the pre-eruptive unrest on El Hierro (Canary Islands): implications of seismicity and deformation in the short-term volcanic hazard assessment, Nat. Hazards Earth Syst. Sci., 18, 1759–1770, 10.5194/nhess-18-1759-2018, 2018.Bedoya-Velásquez, A., Herreras-Giralda, M., Román, R., Wiegner, M., Lefebvre, S., Toledano, C., Huet, T., and Ceolato, R.: Ceilometer inversion method using water-vapor correction from co-located microwave radiometer for aerosol retrievals, Atmos. Res., 250, 105379, 10.1016/j.atmosres.2020.105379, 2021.Bedoya-Velásquez, A. E., Hoyos-Restrepo, M., Barreto, A., García, R. D., Romero-Campos, P. M., García, O., Ramos, R., Roininen, R., Toledano, C., Sicard, M., and Ceolato, R.: Estimation of the Mass Concentration of Volcanic Ash Using Ceilometers: Study of Fresh and Transported Plumes from La Palma Volcano, Remote Sensing, 14, 5680, 10.3390/rs14225680, 2022.Bengtsson, L., Andrae, U., Aspelien, T., Batrak, Y., Calvo, J., de Rooy, W., Gleeson, E., Hansen-Sass, B., Homleid, M., Hortal, M., Ivarsson, K.-I., Lenderink, G., Niemelä, S., Nielsen, K. P., Onvlee, J., Rontu, L., Samuelsson, P., Muñoz, D. S., Subias, A., Tijm, S., Toll, V., Yang, X., and Køltzow, M. Ø.: The HARMONIE-AROME Model Configuration in the ALADIN-HIRLAM NWP System, Mon. Weather Rev., 145, 1919–1935, 10.1175/MWR-D-16-0417.1, 2017.Benito, M. B., Alvarado, G. E., Marchamalo, M., Rejas, J. G., Murphy, P., Franco, R., Castro, D., Garcia-Lanchares, C., and Sanchez, J.: Temporal and spatial evolution of the 2021 eruption in the Tajogaite volcano (Cumbre Vieja rift zone, La Palma, Canary Islands) from geophysical and geodetic parameter analyses, Nat. Hazards, 118, 2245–2284, 10.1007/s11069-023-06090-y, 2023.Boichu, M., Favez, O., Riffault, V., Petit, J.-E., Zhang, Y., Brogniez, C., Sciare, J., Chiapello, I., Clarisse, L., Zhang, S., Pujol-Söhne, N., Tison, E., Delbarre, H., and Goloub, P.: Large-scale particulate air pollution and chemical fingerprint of volcanic sulfate aerosols from the 2014–2015 Holuhraun flood lava eruption of Bárðarbunga volcano (Iceland), Atmos. Chem. Phys., 19, 14253–14287, 10.5194/acp-19-14253-2019, 2019.Boucher, O., Randall, D., Artaxo, P., Bretherton, C., Feingold, G., Forster, P., Kerminen, V.-M., Kondo, Y., Liao, H., Lohmann, U., Rasch, P., Satheesh, S., Sherwood, S., Stevens, B., and Zhang, X.: Chapter 7 - Clouds and Aerosols, Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 571–658, 10.1017/CBO9781107415324.016, ISBN 978-1-107-66182-0, 2013.Brooks, I.: Finding Boundary Layer Top: Application of a Wavelet Covariance Transform to Lidar Backscatter Profiles, J. Atmos. Ocean. Tech., 20, 10.1175/1520-0426(2003)020<1092:FBLTAO>2.0.CO;2, 2003.Campbell, J. R., Hlavka, D. L., Welton, E. J., Flynn, C. J., Turner, D. D., Spinhirne, J. D., Scott, V. S., and Hwang, I. H.: Full-Time, Eye-Safe Cloud and Aerosol Lidar Observation at Atmospheric Radiation Measurement Program Sites: Instruments and Data Processing, J. Atmos. Ocean. Tech., 19, 431–442, 10.1175/1520-0426(2002)019<0431:FTESCA>2.0.CO;2, 2002.Carn, S.: Multi-Satellite Volcanic Sulfur Dioxide L4 Long-Term Global Database V4 (MSVOLSO2L4 product), Goddard Earth Science Data and Information Services Center (GES DISC) [data set], Greenbelt, MD, USA, 10.5067/MEASURES/SO2/DATA405, 2025.Carn, S., Clarisse, L., and Prata, A.: Multi-decadal satellite measurements of global volcanic degassing, J. Volcanol. Geoth. Res., 311, 99–134, 10.1016/j.jvolgeores.2016.01.002, 2016.Carn, S. A., Yang, K., Prata, A. J., and Krotkov, N. A.: Extending the long-term record of volcanic SO2 emissions with the Ozone Mapping and Profiler Suite nadir mapper, Geophys. Res. Lett., 42, 925–932, 10.1002/2014GL062437, 2015.Carn, S. A., Fioletov, V. E., Mclinden, C. A., Li, C., and Krotkov, N. A.: A decade of global volcanic SO2 emissions measured from space, Scientific Reports, 7, 44095, 10.1038/srep44095, 2017.Carracedo, J. C., Troll, V. R., Day, J. M. D., Geiger, H., Aulinas, M., Soler, V., Deegan, F. M., Perez-Torrado, F. J., Gisbert, G., Gazel, E., Rodriguez-Gonzalez, A., and Albert, H.: The 2021 eruption of the Cumbre Vieja volcanic ridge on La Palma, Canary Islands, Geology Today, 38, 94–107, 10.1111/gto.12388, 2022.Carrillo, J., Guerra, J. C., Cuevas, E., and Barrancos, J.: Characterization of the Marine Boundary Layer and the Trade-Wind Inversion over the Sub-tropical North Atlantic, Bound.-Lay. Meteorol., 158, 311–330, 10.1007/s10546-015-0081-1, 2016.Cazorla, A., Casquero-Vera, J. A., Román, R., Guerrero-Rascado, J. L., Toledano, C., Cachorro, V. E., Orza, J. A. G., Cancillo, M. L., Serrano, A., Titos, G., Pandolfi, M., Alastuey, A., Hanrieder, N., and Alados-Arboledas, L.: Near-real-time processing of a ceilometer network assisted with sun-photometer data: monitoring a dust outbreak over the Iberian Peninsula, Atmos. Chem. Phys., 17, 11861–11876, 10.5194/acp-17-11861-2017, 2017.Chen, S., Tong, B., Dong, C., Wang, F., Chen, B., Cheng, C., Zhang, X., and Liu, D.: Retrievals of aerosol layer height during dust events over the taklimakan and gobi desert, J. Quant. Spectrosc. Ra., 254, 107198, 10.1016/j.jqsrt.2020.107198, 2020.Clarisse, L., Hurtmans, D., Clerbaux, C., Hadji-Lazaro, J., Ngadi, Y., and Coheur, P.-F.: Retrieval of sulphur dioxide from the infrared atmospheric sounding interferometer (IASI), Atmos. Meas. Tech., 5, 581–594, 10.5194/amt-5-581-2012, 2012.Comerón, A., Sicard, M., and Rocadenbosch, F.: Wavelet Correlation Transform Method and Gradient Method to Determine Aerosol Layering from Lidar Returns: Some Comments, J. Atmos. Ocean. Tech., 30, 1189–1193, 10.1175/JTECH-D-12-00233.1, 2013.Córdoba-Jabonero, C., Sicard, M., Ansmann, A., del Águila, A., and Baars, H.: 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, Atmos. Meas. Tech., 11, 4775–4795, 10.5194/amt-11-4775-2018, 2018.Córdoba-Jabonero, C., Sicard, M., Barreto, Á., Toledano, C., López-Cayuela, M. Á., Gil-Díaz, C., García, O., Carvajal-Pérez, C. V., Comerón, A., Ramos, R., Muñoz-Porcar, C., and Rodríguez-Gómez, A.: Fresh volcanic aerosols injected in the atmosphere during the volcano eruptive activity at the Cumbre Vieja area (La Palma, Canary Islands): Temporal evolution and vertical impact, Atmos. Environ., 300, 119667, 10.1016/j.atmosenv.2023.119667, 2023.del Fresno, C., Cesca, S., Domínguez Cerdeña, I., Díaz-Suarez, E., Milkereit, C., Valenzuela, C., López-Díaz, R., Dahm, T., and López, C.: Complex seismicity patterns accompanying the 2021 volcanic eruption at La Palma, Canary Islands, Spain, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-9449, 10.5194/egusphere-egu22-9449, 2022. Endo, E. T. and Murray, T. L.: Real-time Seismic Amplitude Measurement (RSAM): a volcano monitoring and prediction tool, B. Volcanol., 53, 533–545, 1991.Esse, B., Burton, M., Hayer, C., La Spina, G., Cofrades, A. P., Asensio-Ramos, M., Barrancos, J., and Pérez, N.: Forecasting the evolution of the 2021 Tajogaite eruption, La Palma, with TROPOMI/PlumeTraj-derived SO2 emission rates, B. Volcanol., 87, 20, 10.1007/s00445-025-01803-6, 2025a.Esse, B., Smith, J., and Brown, T.: Dataset accompanying: Volcanic SO2 emissions from the Tajogaite eruption (La Palma), Figshare [data set], 10.48420/26161543, 2025b.Fedkin, N. M., Li, C., Krotkov, N. A., Hedelt, P., Loyola, D. G., Dickerson, R. R., and Spurr, R.: Volcanic SO2 effective layer height retrieval for the Ozone Monitoring Instrument (OMI) using a machine-learning approach, Atmos. Meas. Tech., 14, 3673–3691, 10.5194/amt-14-3673-2021, 2021.Felpeto, A., Molina-Arias, A. J., Quirós, F., Pereda, J., García-Cañada, L., and Díaz-Suárez, E. A.: Measuring the height of the eruptive column during the 2021 eruption of Cumbre Vieja (La Palma Island, Canary Islands), EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-9419, 10.5194/egusphere-egu22-9419, 2022.Flamant, C., Pelon, J., Flamant, P., and DURAND, P.: Lidar Determination of the Entrainment Zone Thickness at the Top of the Unstable Marine Atmospheric Boundary Layer, Bound.-Lay. Meteorol., 83, 247–284, 10.1023/A:1000258318944, 1997.Flynn, C. J., Mendoza, A., Zheng, Y., and Mathur, S.: Novel polarization-sensitive micropulse lidar measurement technique, Opt. Express, 15, 2785–2790, 10.1364/OE.15.002785, 2007.Forster, P., Storelvmo, T., Armour, K., Collins, W., Dufresne, J.-L., Frame, D., Lunt, D., Mauritsen, T., Palmer, M., Watanabe, M., Wild, M., and Zhang, H.: The Earth's Energy Budget, Climate Feedbacks, and Climate Sensitivity, Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 923–1054, 10.1017/9781009157896.009, 2021.García, R. D., García, O. E., Cuevas-Agulló, E., Barreto, Á., Cachorro, V. E., Marrero, C., Almansa, F., Ramos, R., and Pó, M.: Spectral Aerosol Radiative Forcing and Efficiency of the La Palma Volcanic Plume over the Izaña Observatory, Remote Sensing, 15, 10.3390/rs15010173, 2023.Gebauer, H., Floutsi, A. A., Haarig, M., Radenz, M., Engelmann, R., Althausen, D., Skupin, A., Ansmann, A., Zenk, C., and Baars, H.: Tropospheric sulfate from Cumbre Vieja (La Palma) observed over Cabo Verde contrasted with background conditions: a lidar case study of aerosol extinction, backscatter, depolarization and lidar ratio profiles at 355, 532 and 1064 nm, Atmos. Chem. Phys., 24, 5047–5067, 10.5194/acp-24-5047-2024, 2024.Girault, F., Carazzo, G., Tait, S., and Kaminski, E.: Combined effects of total grain-size distribution and crosswind on the rise of eruptive volcanic columns, J. Volcanol. Geoth. Res., 326, 103–113, 10.1016/j.jvolgeores.2015.11.007, 2016.Hedelt, P., Efremenko, D. S., Loyola, D. G., Spurr, R., and Clarisse, L.: Sulfur dioxide layer height retrieval from Sentinel-5 Precursor/TROPOMI using FP_ILM, Atmos. Meas. Tech., 12, 5503–5517, 10.5194/amt-12-5503-2019, 2019.Hedelt, P., Reichardt, J., Lauermann, F., Weiß, B., Theys, N., Redondas, A., Barreto, A., Garcia, O., and Loyola, D.: Analysis of the long-range transport of the volcanic plume from the 2021 Tajogaite/Cumbre Vieja eruption to Europe using TROPOMI and ground-based measurements, Atmos. Chem. Phys., 25, 1253–1272, 10.5194/acp-25-1253-2025, 2025.Heese, B., Flentje, H., Althausen, D., Ansmann, A., and Frey, S.: Ceilometer lidar comparison: backscatter coefficient retrieval and signal-to-noise ratio determination, Atmos. Meas. Tech., 3, 1763–1770, 10.5194/amt-3-1763-2010, 2010.Iarlori, M., Pietropaolo, E., Rizi, V., Aramo, C., Valore, L., Silvestri, V., Cirella, A., Dughera, G., Galli, S., and Marengo, M.: The Raman LIDAR for the pre-production phase of Cherenkov Telescope Array, EPJ Web Conf., 197, 02004, 10.1051/epjconf/201919702004, 2019.Jenoptik: CHM15k – Nimbus Ceilometer User Manual, Revision P0, https://www.lufft.com/fileadmin/lufft.com/07_downloads/1_manuals/Manual_CHM15k_EN_R19.pdf (last access: 27 November 2025), 2013.Kampouri, A., Amiridis, V., Solomos, S., Gialitaki, A., Marinou, E., Spyrou, C., Georgoulias, A., Akritidis, D., Papagiannopoulos, N., Mona, L., Scollo, S., Tsichla, M., Tsikoudi, I., Pytharoulis, I., Karacostas, T., and Zanis, P.: Investigation of Volcanic Emissions in the Mediterranean: “The Etna–Antikythera Connection”, Atmosphere, 12, 40, 10.3390/atmos12010040, 2020.Karagulian, F., Clarisse, L., Clerbaux, C., Prata, A. J., Hurtmans, D., and Coheur, P. F.: Detection of volcanic SO2, ash, and H2SO4 using the Infrared Atmospheric Sounding Interferometer (IASI), J. Geophys. Res.-Atmos., 115, 10.1029/2009JD012786, 2010.Kim, S.-W., Berthier, S., Raut, J.-C., Chazette, P., Dulac, F., and Yoon, S.-C.: Validation of aerosol and cloud layer structures from the space-borne lidar CALIOP using a ground-based lidar in Seoul, Korea, Atmos. Chem. Phys., 8, 3705–3720, 10.5194/acp-8-3705-2008, 2008.King, M. D., Kaufman, Y. J., Menzel, W. P., and Tanre, D.: Remote sensing of cloud, aerosol, and water vapor properties from the Moderate Resolution Imaging Spectroradiometer (MODIS), IEEE T. Geosci. Remote, 41, 442–458, 10.1109/TGRS.2002.808226, 2003.Kotthaus, S., O'Connor, E., Münkel, C., Charlton-Perez, C., Haeffelin, M., Gabey, A. M., and Grimmond, C. S. B.: Recommendations for processing atmospheric attenuated backscatter profiles from Vaisala CL31 ceilometers, Atmos. Meas. Tech., 9, 3769–3791, 10.5194/amt-9-3769-2016, 2016.Koukouli, M.-E., Michailidis, K., Hedelt, P., Taylor, I. A., Inness, A., Clarisse, L., Balis, D., Efremenko, D., Loyola, D., Grainger, R. G., and Retscher, C.: Volcanic SO2 layer height by TROPOMI/S5P: evaluation against IASI/MetOp and CALIOP/CALIPSO observations, Atmos. Chem. Phys., 22, 5665–5683, 10.5194/acp-22-5665-2022, 2022.Kremser, S., Thomason, L. W., von Hobe, M., Hermann, M., Deshler, T., Timmreck, C., Toohey, M., Stenke, A., Schwarz, J. P., Weigel, R., Fueglistaler, S., Prata, F. J., Vernier, J.-P., Schlager, H., Barnes, J. E., Antuña-Marrero, J.-C., Fairlie, D., Palm, M., Mahieu, E., Notholt, J., Rex, M., Bingen, C., Vanhellemont, F., Bourassa, A., Plane, J. M. C., Klocke, D., Carn, S. A., Clarisse, L., Trickl, T., Neely, R., James, A. D., Rieger, L., Wilson, J. C., and Meland, B.: Stratospheric aerosol – Observations, processes, and impact on climate, Rev. Geophys., 54, 278–335, 10.1002/2015RG000511, 2016.Krotkov, N., Realmuto, V., Li, C., Seftor, C., Li, J., Brentzel, K., Stuefer, M., Cable, J., Dierking, C., Delamere, J., Schneider, D., Tamminen, J., Hassinen, S., Ryyppö, T., Murray, J., Carn, S., Osiensky, J., Eckstein, N., Layne, G., and Kirkendall, J.: Day–Night Monitoring of Volcanic SO2 and Ash Clouds for Aviation Avoidance at Northern Polar Latitudes, Remote Sensing, 13, 4003, 10.3390/rs13194003, 2021.Lamb, O., De Angelis, S., and Lavallée, Y.: Using infrasound to constrain ash plume rise, Journal of Applied Volcanology, 4, 20, 10.1186/s13617-015-0038-6, 2015.Li, C., Krotkov, N. A., Carn, S., Zhang, Y., Spurr, R. J. D., and Joiner, J.: New-generation NASA Aura Ozone Monitoring Instrument (OMI) volcanic SO2 dataset: algorithm description, initial results, and continuation with the Suomi-NPP Ozone Mapping and Profiler Suite (OMPS), Atmos. Meas. Tech., 10, 445–458, 10.5194/amt-10-445-2017, 2017.Li, C., Krotkov, N. A., Leonard, P. J. T., Carn, S., Joiner, J., Spurr, R. J. D., and Vasilkov, A.: Version 2 Ozone Monitoring Instrument SO2 product (OMSO2 V2): new anthropogenic SO2 vertical column density dataset, Atmos. Meas. Tech., 13, 6175–6191, 10.5194/amt-13-6175-2020, 2020.Li, C., Krotkov, N. A., Leonard, P. J. T., and Joiner, J.: OMI/Aura Sulphur Dioxide (SO2) Total Column 1-orbit L2 Swath 13x24 km V004, Goddard Earth Sciences Data and Information Services Center (GES DISC) [data set], Greenbelt, MD, USA, 10.5067/Aura/OMI/DATA2422, 2025.Lopez, C., Blanco, M. J., and Team, I.: Instituto Geográfico Nacional Volcano Monitoring of the 2021 La Palma eruption (Canary Islands, Spain), EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-11549, 10.5194/egusphere-egu22-11549, 2022.López, C., Benito-Saz, M. A., Martí, J., del-Fresno, C., García-Cañada, L., Albert, H., and Lamolda, H.: Driving magma to the surface: The 2011–2012 El Hierro Volcanic Eruption, Geochem. Geophy. Geosy., 18, 3165–3184, 10.1002/2017GC007023, 2017.Marshall, L. R., Smith, C. J., Forster, P. M., Aubry, T. J., Andrews, T., and Schmidt, A.: Large Variations in Volcanic Aerosol Forcing Efficiency Due to Eruption Source Parameters and Rapid Adjustments, Geophys. Res. Lett., e2020GL090241, 10.1029/2020GL090241, 2020.Masson-Delmotte, V., Zhai, P., Pirani, A., Connors, S., Péan, C., Berger, S., Caud, N., Chen, Y., Goldfarb, L., Gomis, M., Huang, M., Leitzell, K., Lonnoy, E., Matthews, J., Maycock, T., Waterfield, T., Yelekçi, O., Yu, R., and Zhou, B.: Climate Change 2021, Full Report, Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 10.1017/9781009157896, 2023.Mather, T., Pyle, D., and Oppenheimer, C.: Tropospheric Volcanic Aerosol, Washington DC American Geophysical Union Geophysical Monograph Series, 139, 189–212, 10.1029/139GM12, 2003. Measures, R. M.: Laser remote sensing: Fundamentals and applications(Book), New York, Wiley-Interscience, 521 pp., ISBN 10: 0471081930, 1984.Melfi, S., Spinhirne, J., Chou, S.-H., and Palm, S.: Lidar Observations of Vertically Organized Convection in the Planetary Boundary Layer over the Ocean, J. Appl. Meteorol., 24, 806–821, 10.1175/1520-0450(1985)024<0806:LOOVOC>2.0.CO;2, 1985.Milford, C., Cuevas, E., Marrero, C. L., Bustos, J., Gallo, V., Rodríguez, S., Romero-Campos, P. M., and Torres, C.: Impacts of Desert Dust Outbreaks on Air Quality in Urban Areas, Atmosphere, 11, 23, 10.3390/atmos11010023, 2020.Milford, C., Torres, C., Vilches, J., Gossman, A.-K., Weis, F., Suárez-Molina, D., García, O. E., Prats, N., Barreto, Á., García, R. D., Bustos, J. J., Marrero, C. L., Ramos, R., Chinea, N., Boulesteix, T., Taquet, N., Rodríguez, S., López-Darias, J., Sicard, M., Córdoba-Jabonero, C., and Cuevas, E.: Impact of the 2021 La Palma Volcanic Eruption on Air Quality: Insights from a Multidisciplinary Approach, Sci. Total Environ., 869, 161652, 10.1016/j.scitotenv.2023.161652, 2023.Morille, Y., Haeffelin, M., Drobinski, P., and Pelon, J.: STRAT: An Automated Algorithm to Retrieve the Vertical Structure of the Atmosphere from Single-Channel Lidar Data, J. Atmos. Ocean. Tech., 24, 761–775, 10.1175/JTECH2008.1, 2007.Nanda, S., de Graaf, M., Veefkind, J. P., Sneep, M., ter Linden, M., Sun, J., and Levelt, P. F.: A first comparison of TROPOMI aerosol layer height (ALH) to CALIOP data, Atmos. Meas. Tech., 13, 3043–3059, 10.5194/amt-13-3043-2020, 2020.Navas-Guzmán, F., Müller, D., Bravo-Aranda, J. A., Guerrero-Rascado, J. L., Granados-Muñoz, M. J., Pérez-Ramírez, D., Olmo, F. J., and Alados-Arboledas, L.: Eruption of the Eyjafjallajökull Volcano in spring 2010: Multiwavelength Raman lidar measurements of sulphate particles in the lower troposphere, J. Geophys. Res.-Atmos., 118, 1804–1813, 10.1002/jgrd.50116, 2013.Nogales, M., Guerrero-Campos, M., Boulesteix, T., Taquet, N., Beierkuhnlein, C., Campion, R., Fajardo, S., Zurita, N., Arechavaleta, M., García, R., Weiser, F., and Medina, F.: The fate of terrestrial biodiversity during an oceanic island volcanic eruption, Scientific Reports, 12, 19344, 10.1038/s41598-022-22863-0, 2022.O'Connor, E., Illingworth, A., and Hogan, R.: A Technique for Autocalibration of Cloud Lidar, J. Atmos. Ocean. Tech., 21, 10.1175/1520-0426(2004)021<0777:ATFAOC>2.0.CO;2, 2004.Palma, J. L., Calder, E. S., Basualto, D., Blake, S., and Rothery, D. A.: Correlations between SO2 flux, seismicity, and outgassing activity at the open vent of Villarrica volcano, Chile, J. Geophys. Res.-Sol. Ea., 113, 10.1029/2008JB005577, 2008.Pappalardo, G., Amodeo, A., Mona, L., Pandolfi, M., Pergola, N., and Cuomo, V.: Raman lidar observations of aerosol emitted during the 2002 Etna eruption, Geophys. Res. Lett., 31, L05120, 10.1029/2003GL019073, 2004.Pardini, F., Burton, M., de' Michieli Vitturi, M., Corradini, S., Salerno, G., Merucci, L., and Di Grazia, G.: Retrieval and intercomparison of volcanic SO2 injection height and eruption time from satellite maps and ground-based observations, J. Volcanol. Geoth. Res., 331, 79–91, 10.1016/j.jvolgeores.2016.12.008, 2017.Perrone, M., De Tomasi, F., and Burlizzi, P.: Aerosol products by CALIOP at 532 nm and by a ground-based Raman lidar at 355 nm: Intercomparison methodology, Atmos. Res., 101, 438–449, 10.1016/j.atmosres.2011.04.005, 2011.PEVOLCA: Actualización de la actividad volcánica en Cumbre Vieja (La Palma), https://info.igme.es/eventos/la-palma/pevolca/211225_Informe-Comite-Cientifico-PEVOLCA.pdf (last access: 25 December 2021), 2021.Plu, M., Bigeard, G., Sič, B., Emili, E., Bugliaro, L., El Amraoui, L., Guth, J., Josse, B., Mona, L., and Piontek, D.: Modelling the volcanic ash plume from Eyjafjallajökull eruption (May 2010) over Europe: evaluation of the benefit of source term improvements and of the assimilation of aerosol measurements, Nat. Hazards Earth Syst. Sci., 21, 3731–3747, 10.5194/nhess-21-3731-2021, 2021.Prata, A. J. and Bernardo, C.: Retrieval of volcanic SO2 column abundance from Atmospheric Infrared Sounder data, J. Geophys. Res.-Atmos., 112, D20204, 10.1029/2006JD007955, 2007.Realmuto, V. J. and Berk, A.: Plume Tracker: Interactive mapping of volcanic sulfur dioxide emissions with high-performance radiative transfer modeling, J. Volcanol. Geoth. Res., 327, 55–69, 10.1016/j.jvolgeores.2016.07.001, 2016.Rodríguez, F., Pérez, N. M., Amonte, C., Martín-Lorenzo, A., Melián, G. V., Coldwell, B. C., Pankhurst, M. J., Asensio-Ramos, M., Hernández, P. A., and Padrón, E.: Geochemistry of ash leachates during the 2021 eruption of Cumbre Vieja volcano, La Palma, Canary Islands, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-9629, 10.5194/egusphere-egu22-9629, 2022.Rodríguez, O., Barrancos, J., Cutillas, J., Ortega, V., Hernández, P. A., Cabrera, I., and Pérez, N. M.: SO2 emissions during the post-eruptive phase of the Tajogaite eruption (La Palma, Canary Islands) by means of ground-based miniDOAS measurements in transverse mode using a car and UAV, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-3620, 10.5194/egusphere-egu23-3620, 2023.Romero, J. E., Burton, M., Cáceres, F., Taddeucci, J., Civico, R., Ricci, T., Pankhurst, M. J., Hernández, P. A., Bonadonna, C., Llewellin, E. W., Pistolesi, M., Polacci, M., Solana, C., D'Auria, L., Arzilli, F., Andronico, D., Rodríguez, F., Asensio-Ramos, M., Martín-Lorenzo, A., Hayer, C., Scarlato, P., and Perez, N. M.: The Initial Phase of the 2021 Cumbre Vieja Ridge Eruption (Canary Islands): Products and Dynamics Controlling Edifice Growth and Collapse, J. Volcanol. Geoth. Res., 431, 107642, 10.1016/j.jvolgeores.2022.107642, 2022.Román, R., Cazorla, A., Toledano, C., Olmo, F., Cachorro, V., Frutos, A., and Arboledas, L.: Cloud cover detection combining high dynamic range sky images and ceilometer measurements, Atmos. Res., 196, 224–236, 10.1016/j.atmosres.2017.06.006, 2017. Román, R., González, R., Antuña-Sánchez, J.-C., Barreto, A., Martín, P., Toledano, C., Ramos, R., Cazorla, A., Herrero-Anta, S., Mateos, D., García, O. E., González-Fernández, D., Carracedo, R., Herreras-Giralda, M., Carreño, V., Calle, A., Cachorro, V. E., Cuevas, E., and de Frutos, Ã. M.: Vertical profiles of aerosol properties retrieved at La Palma, Canary Islands, during the Cumbre-Vieja volcano eruption in September-October 2021, in: Conference European Lidar Conference 2021 (ELC2021), p. S05P14, 2021.Rossi, E., Bonadonna, C., and Degruyter, W.: A new strategy for the estimation of plume height from clast dispersal in various atmospheric and eruptive conditions, Earth Planet. Sc. Lett., 505, 1–12, 10.1016/j.epsl.2018.10.007, 2019.Salgueiro, V., Guerrero-Rascado, J. L., Costa, M. J., Román, R., Cazorla, A., Serrano, A., Molero, F., Sicard, M., Córdoba-Jabonero, C., Bortoli, D., Comerón, A., Couto, F. T., López-Cayuela, M. Á., Pérez-Ramírez, D., Potes, M., Muñiz-Rosado, J. A., Obregón, M. A., Barragán, R., Oliveira, D. C. F. S., Abril-Gago, J., González, R., Gíl-Díaz, C., Foyo-Moreno, I., Muñoz-Porcar, C., Granados-Muñoz, M. J., Rodríguez-Gómez, A., Herreras-Giralda, M., Bravo-Aranda, J. A., Carvajal-Pérez, C. V., Barreto, A., and Alados-Arboledas, L.: Characterization of Tajogaite Volcanic Plumes Detected over the Iberian Peninsula from a Set of Satellite and Ground-Based Remote Sensing Instrumentation, Remote Sens. Environ., 295, 113684, 10.1016/j.rse.2023.113684, 2023.Scollo, S., Prestifilippo, M., Pecora, E., Corradini, S., Merucci, L., Spata, G., and Coltelli, M.: Eruption Column Height Estimation: the 2011–2013 Etna lava fountains, Ann. Geophys., 57, S0214, 10.4401/ag-6396, 2014.Sellitto, P., di Sarra, A., Corradini, S., Boichu, M., Herbin, H., Dubuisson, P., Sèze, G., Meloni, D., Monteleone, F., Merucci, L., Rusalem, J., Salerno, G., Briole, P., and Legras, B.: Synergistic use of Lagrangian dispersion and radiative transfer modelling with satellite and surface remote sensing measurements for the investigation of volcanic plumes: the Mount Etna eruption of 25–27 October 2013, Atmos. Chem. Phys., 16, 6841–6861, 10.5194/acp-16-6841-2016, 2016.Sicard, M., Guerrero-Rascado, J. L., Navas-Guzmán, F., Preißler, J., Molero, F., Tomás, S., Bravo-Aranda, J. A., Comerón, A., Rocadenbosch, F., Wagner, F., Pujadas, M., and Alados-Arboledas, L.: Monitoring of the Eyjafjallajökull volcanic aerosol plume over the Iberian Peninsula by means of four EARLINET lidar stations, Atmos. Chem. Phys., 12, 3115–3130, 10.5194/acp-12-3115-2012, 2012.Sicard, M., Córdoba-Jabonero, C., Barreto, A., Welton, E. J., Gil-Díaz, C., Carvajal-Pérez, C. V., Comerón, A., García, O., García, R., López-Cayuela, M.-Á., Muñoz-Porcar, C., Prats, N., Ramos, R., Rodríguez-Gómez, A., Toledano, C., and Torres, C.: Volcanic Eruption of Cumbre Vieja, La Palma, Spain: A First Insight to the Particulate Matter Injected in the Troposphere, Remote Sensing, 14, 2470, 10.3390/rs14102470, 2022.Taquet, N., Boulesteix, T., García, O., Campion, R., Stremme, W., Rodríguez, S., López-Darias, J., Marrero, C., González-García, D., Klügel, A., Hase, F., García, M. I., Ramos, R., Rivas-Soriano, P., Léon-Luis, S., Carreño, V., Alcántara, A., Sépulveda, E., Milford, C., González-Sicilia, P., and Torres, C.: Volatile emissions during the 2021 Cumbre Vieja (La Palma) eruption integrating multiplatform atmospheric observations, Atmos. Chem. Phys., 25, 14591–14628, 10.5194/acp-25-14591-2025, 2025.Taylor, I. A., Preston, J., Carboni, E., Mather, T. A., Grainger, R. G., Theys, N., Hidalgo, S., and Kilbride, B. M.: Exploring the Utility of IASI for Monitoring Volcanic SO2 Emissions, J. Geophys. Res.-Atmos., 123, 5588–5606, 10.1002/2017JD027109, 2018.Theys, N., Campion, R., Clarisse, L., Brenot, H., van Gent, J., Dils, B., Corradini, S., Merucci, L., Coheur, P.-F., Van Roozendael, M., Hurtmans, D., Clerbaux, C., Tait, S., and Ferrucci, F.: Volcanic SO2 fluxes derived from satellite data: a survey using OMI, GOME-2, IASI and MODIS, Atmos. Chem. Phys., 13, 5945–5968, 10.5194/acp-13-5945-2013, 2013.Theys, N., De Smedt, I., Yu, H., Danckaert, T., van Gent, J., Hörmann, C., Wagner, T., Hedelt, P., Bauer, H., Romahn, F., Pedergnana, M., Loyola, D., and Van Roozendael, M.: Sulfur dioxide retrievals from TROPOMI onboard Sentinel-5 Precursor: algorithm theoretical basis, Atmos. Meas. Tech., 10, 119–153, 10.5194/amt-10-119-2017, 2017.Theys, N., Fioletov, V., Li, C., De Smedt, I., Lerot, C., McLinden, C., Krotkov, N., Griffin, D., Clarisse, L., Hedelt, P., Loyola, D., Wagner, T., Kumar, V., Innes, A., Ribas, R., Hendrick, F., Vlietinck, J., Brenot, H., and Van Roozendael, M.: A sulfur dioxide Covariance-Based Retrieval Algorithm (COBRA): application to TROPOMI reveals new emission sources, Atmos. Chem. Phys., 21, 16727–16744, 10.5194/acp-21-16727-2021, 2021.Tournigand, P.-Y., Cigala, V., Prata, F., Steiner, A., Kirchengast, G., Brenot, H., Clarisse, L., and Biondi, R.: The 2015 Calbuco Volcanic Cloud Detection Using GNSS Radio Occultation and Satellite Lidar, in: International Geoscience and Remote Sensing (IGARSS) Symposium, 10.1109/IGARSS39084.2020.9323356, 2020.VAAC: Toulouse Volcanic Ash Advisory Center Reports for La Palma, VAAC, http://vaac.meteo.fr/volcanoes/la-palma/ (last access: 1 January 2024), 2022.Valore, L., Aramo, C., Dinelli, B., Di Pierro, F., Dughera, G., Gaug, M., Iarlori, M., Marengo, M., Papandrea, E., Pietropaolo, E., Rizi, V., Vallania, P., and Vigorito, C.: The ARCADE Raman Lidar and atmospheric simulations for the Cherenkov Telescope Array, PoS, ICRC2017, 763 pp., 10.22323/1.301.0763, 2017.Vaughan, M. A., Powell, K. A., Winker, D. M., Hostetler, C. A., Kuehn, R. E., Hunt, W. H., Getzewich, B. J., Young, S. A., Liu, Z., and McGill, M. J.: Fully Automated Detection of Cloud and Aerosol Layers in the CALIPSO Lidar Measurements, J. Atmos. Ocean. Tech., 26, 2034–2050, 10.1175/2009JTECHA1228.1, 2009. Welton, E. J. and Campbell, J. R.: Micropulse Lidar Signals: Uncertainty Analysis, J. Atmos. Ocean. Tech., 19, 2089–2094, 10.1175/1520-0426(2002)019<2089:MLSUA>2.0.CO;2, 2002.Winker, D. M., Vaughan, M. A., Omar, A., Hu, Y., Powell, K. A., Liu, Z., Hunt, W. H., and Young, S. A.: Overview of the CALIPSO Mission and CALIOP Data Processing Algorithms, J. Atmos. Ocean. Tech., 26, 2310–2323, 10.1175/2009JTECHA1281.1, 2009.Winker, D. M., Tackett, J. L., Getzewich, B. J., Liu, Z., Vaughan, M. A., and Rogers, R. R.: The global 3-D distribution of tropospheric aerosols as characterized by CALIOP, Atmos. Chem. Phys., 13, 3345–3361, 10.5194/acp-13-3345-2013, 2013.Woodhouse, M. J., Hogg, A. J., and Phillips, J. C.: A global sensitivity analysis of the PlumeRise model of volcanic plumes, J. Volcanol. Geoth. Res., 326, 54–76, 10.1016/j.jvolgeores.2016.02.019, 2016.Zhu, Y., Toon, O., Jensen, E., Bardeen, C., Mills, M., Tolbert, M., Yu, P., and Woods, S.: Persisting volcanic ash particles impact stratospheric SO2 lifetime and aerosol optical properties, Nat. Commun., 11, 4526, 10.1038/s41467-020-18352-5, 2020.