Skip to main content

Health impact of the Tajogaite volcano eruption in La Palma population (ISVOLCAN study): rationale, design, and preliminary results from the first 1002 participants

Abstract

Background

The eruption of the Tajogaite volcano began on the island of La Palma on September 19, 2021, lasting for 85 days. This study aims to present the design and methodology of the ISVOLCAN (Health Impact on the Population of La Palma due to the Volcanic Eruption) cohort, as well as the preliminary findings from the first 1002 enrolled participants.

Methods

A prospective cohort study was conducted with random selection of adult participants from the general population, with an estimated sample size of 2600 individuals. The results of the first 857 participants are presented, along with a group of 145 voluntary participants who served as interveners during the eruption. Data on epidemiology and volcano exposure were collected, and participants underwent physical examinations, including anthropometry, blood pressure measurement, spirometry, and venous blood extraction for toxicological assessment.

Results

In the general population (n = 857), descriptive analysis revealed that the participants were mostly middle-aged individuals (50.8 ± 16.4), with a predominance of females. Before the eruption, the participants resided at a median distance of 6.7 km from the volcano in the Western region and 10.9 km in the Eastern region. Approximately 15.4% of the sample required evacuation, whose 34.8% returning to their homes on average after 3 months. A significant number of participants reported engaging in daily tasks involving cleaning of volcanic ash both indoors and outdoors. The most reported acute symptoms included ocular irritation, insomnia, mood disorders (anxiety-depression), and respiratory symptoms. Multivariate analysis results show that participants in the western region had a higher likelihood of lower respiratory tract symptoms (OR 1.99; 95% CI:1.33–2.99), depression and anxiety (OR 1.95; 95% CI:1.30–2.93), and insomnia (OR 2.03; 95% CI:1.33–3.09), compared to those in the eastern region.

Conclusion

The ongoing follow-up of the ISVOLCAN cohort will provide valuable insights into the short, medium, and long-term health impact related to the material emitted during the Tajogaite eruption, based on the level of exposure suffered by the affected population.

Peer Review reports

Background

Approximately one billion people worldwide live within the influence zone of an active volcano, at about 100 km [1], and thus could be affected by the effects of an eruption at some point. On La Palma Island (Canary Islands, Spain), a volcanic eruption began on September 19, 2021, in the Valle de Aridane, lasting for 85 days and resulting in the formation of a new volcano named Tajogaite. The eruption generated a significant expulsion of volcanic ash and gas emissions, leading to days of highly unfavourable air quality with elevated toxicity levels in the breathable air [2].

Volcanic eruptions can have a wide range of deleterious effects on human health. Despite their often-short duration, the emission of toxic gases, particles, and ash deposits can persist in the local environment for years or even decades, being mobilized and redistributed by climatic factors or human activities [3]. Gases emitted during volcanic activity, such as CO and CO2, SO2, HCl, HF, H2S, radon, and permanent degassing, have the potential to impact human health [4]. There are several causes for degassing, both natural and anthropogenic sources whose contribute to air pollution. In La Palma, several years before the Tajogaite eruption, the concentration of CO2 were recorded in Cumbre Vieja, beeing related, in much more amount, to anthropogenic and other natural sources than magmatic emissions, accounting for only 4% [5]. Acute and prolongated exposure of individuals to high concentrations of CO2 is uncommon and incompatible with life. Nevertheless, long-term exposure occurs at low concentrations and outdoor is more frequent. In this context, physiological adaptation mechanisms are activated, promote oxidation followed by the release of proinflammatory cytokines; these mechanisms entail the development of a pro-inflammatory status and, consequently, the onset of diseases related to such conditions [6]. Particularly harmful to health is the atmospheric transformation of SO2 and other gases into particulate matter (PM) [7]. This transformation process is influenced by various factors, including the emission plume’s characteristics and maturity, as well as meteorological variables such as humidity, solar radiation, and temperature [7].

Previous studies have demonstrated an increase in acute symptoms of ocular irritation and upper respiratory tract issues [8], as well as elevated respiratory morbidity and visits to hospital or primary care services due to exacerbation of respiratory pathologies associated with peaks in airborne emissions of these types of toxic gases or particles [3, 9, 10]. Many of these associations are independent of age, sex, education level, and smoking habits, exhibiting a dose-response gradient [11]. Furthermore, exposure has been linked to increased cardiovascular morbidity and all-cause mortality [12]. Few studies have assessed long-term chronic health effects, and the scarce longitudinal studies have methodological limitations due to the analysis of samples from hospitalized patients, low reliability of data sources in some countries, and short follow-up periods, often not exceeding six months. Consequently, longitudinal studies with extended follow-up periods in the general population are needed to analyse the occurrence of deleterious medium to long-term effects.

In January 2022, the ISVOLCAN study (Health Impact on the Population of La Palma caused by the Tajogaite Volcano Eruption) was started. This study involves the recruitment and follow-up of a cohort from the general adult population to assess the impact of the Tajogaite volcano eruption on the health of the population of the island. The purpose of this paper is to present the methodology of ISVOLCAN study and provide a preliminary analysis of the data obtained from the first 1002 enrolled participants.

Methods

Study design

This is an observational epidemiological study, using a prospective cohort design, targeting the general adult population residing in multiple municipalities on La Palma Island. Additionally, a group of volunteers from professional personnel with access to the exclusion zone or operations centre during the eruption (including civil protection workers, Spanish Security Forces, Emergency Services, scientists, etc.) was included. The study consists of two different stages: the first stage involved recruitment and baseline assessments conducted from 2022 to 2023, while the second phase will involve follow-up of the cohort at 2, 5, and 10-year intervals.

The study has obtained authorization from the health authorities and received a favorable decision from the Provincial Ethics and Medicines Committee (ref. CHUNSC_2021_88). Participants were required to provide written consent before being included in the study.

Settings

La Palma Island is a volcanic island located in the Atlantic Ocean, within the Canary Islands archipelago, Spain. Geographically, it is positioned at 28° 26’ N latitude and 14° 01’ W longitude from Madrid. Covering an approximate area of just over 700 km2, it ranks as the fifth-largest island in this archipelago [13]. The island counts with a population of 83,439 inhabitants [14]. Los Llanos de Aridane in the west side, is the city with highest population density, followed of Santa Cruz de La Palma. The Canarian Public Health System provides healthcare to the entire population through a hospital and a network of primary care health centres throughout the island.

On September 19, 2021, the eruption started in the western region of La Palma Island, in Cabeza de Vaca area, on the western flank of the Cumbre Vieja ridge, belowing to the municipality of El Paso. This volcanic process persisted for 85 days until its conclusion on December 13 of the same year. It consisted in a long-lasting, hybrid eruption associated with multiple eruptive styles (effusive, lava fountains, ash emissions, strombolian explosions) with the formation of cones of various heights, widespread tephra blankets and extensive lava-flow fields and was characterized by simultaneous effusive and explosive activity [15]. The eruption affected the Valle de Aridane, which was greatly impacted by the lava flows, gases, and particulate matter emitted during the eruption. The newly formed volcano, named Tajogaite, reached a maximum altitude of 1131 m above sea level and extended 200 m from the pre-eruptive topography, with its base situated at 1080 m above sea level [16].

Subjects, sampling, inclusion and exclusion criteria

Participants from the general population

The sample selection was conducted using a random, stratified approach based on age and gender groups, according to the 2020 municipal census data of the population residing in the western region (El Paso, Los Llanos, Tazacorte, and Puntagorda) and eastern region (Mazo, Santa Cruz de La Palma, and San Andrés y Sauces) of the island.

The sample was drawn from the health card registry of the Canarian Health Service, which is continuously updated and includes all individuals above the age of 18 who were residents on the island during the eruption and provided informed consent to participate in the study.

To ensure the achievement of the intended objectives, the population of the western region, closer to the eruption, was oversampled. The sample sizes for each municipality in the western region were as follows: Los Llanos de Aridane: 820; El Paso: 405; Tazacorte: 305; Puntagorda: 205. In contrast, for the eastern region, the sample sizes were 505 for Santa Cruz de La Palma, 205 for Mazo, and 155 for San Andrés y Sauces.

Highly exposed participants (intervening personnel)

Using a non-probabilistic convenience sampling method, participants in the study also included members of various professional and volunteer groups involved in different tasks related to the eruption and who had access to the volcano’s exclusion zone or operations centre. Although access to these areas was controlled and followed safety and protection measures, we expected that these participants from the different groups were highly exposed during their workdays throughout the nearly four-month duration of the eruption.

Study size

An initial sample size of 1207 persons was estimated (precision 3%, confidence level 95%) based on an expected prevalence of acute respiratory symptoms (the most frequently associated with such phenomena). Considering an anticipated participation rate in this type of study of 60–70% and a dropout rate during follow-up exceeding 30%, the sample was increased to 2600 individuals.

Recruitment and baseline assessment

In January 2022, telephone contact with the selected sample started. Those who agreed to participate in the study were administered an epidemiological questionnaire specifically designed for this purpose. The questionnaire was completed by Primary Care professionals, including both physicians and nurses, who were trained for this study. Subsequently, participants were scheduled to visit the health centres in the two regions for physical examinations, pulmonary function tests, and venous blood extraction, all conducted by qualified nursing staff.

An electronic questionnaire was designed in accordance with the recommendations of the International Volcanic Health Hazard Network (IVHHN), an organization under the World Health Organization (WHO), aimed at standardizing epidemiological protocols for assessing health effects in volcanic eruptions [17].

The questionnaire is available on the website of the study (www.estudioisvolcan.com) and included sociodemographic data (age, gender, employment status, occupation type, educational level), variables related to the level of exposure to the volcano (residence before and during the eruption, need for evacuation and subsequent return to the usual residence, access to exclusion zones, involvement in activities related to volcanic ash cleaning, daily hours spent in outdoor environments, and use of masks and eyeglasses for protection), pre-existing comorbidities (lung diseases, cardiovascular diseases, type 2 diabetes, blood hypertension, etc.), acute symptoms (cough, sneezing, wheezing, headache, fatigue, tearing, ocular irritation, etc.), suffering from any respiratory infection (flu, COVID 19 or cold) and visits to emergency services during the eruption, lifestyle factors (smoking habits and leisure-time physical activity). Additionally, the questionnaire included a shortened version of the scale for assessing post-traumatic stress disorder, adapted for the Spanish population [18].

During the visit to the health centre, measurements of weight, height, waist circumference, heart rate, and two separate blood pressure (separately by 10 min) were recorded. Additionally, a venous blood sample of approximately 20 mL was collected, divided into 4 tubes (2 tubes for complete blood count and 2 tubes for biochemistry), for the toxicological determination of persistent contaminants in whole blood and serum.

The tubes for complete blood were stored in a refrigerator at 4 °C, while the biochemistry tubes were centrifuged at 3000 rpm for 10–15 min, and then allowed to rest for 20–25 min until clot retraction. Daily, the samples were transported to the Laboratory of the University Hospital of La Palma and finally stored at the Research Unit of the Hospital Nuestra Señora de Candelaria in Tenerife at -80 °C for the sera and − 20 °C for the whole blood.

In the blood samples, organic contaminants, primarily polycyclic aromatic hydrocarbons (PAHs), will be quantitatively determined due to their possible formation in eruptive processes and their known carcinogenic and teratogenic properties. Among them, the following will be determined: naphthalene, acenaphthene, acenaphthylene, fluorene, anthracene, phenanthrene, pyrene, fluoranthene, benzo(a)anthracene, chrysene, benzo(b)fluoranthene, benzo(k)fluoranthene, benzo(a)pyrene, benzo(ghi)perylene, indene(1,2,3,cd)pyrene, and dibenzo(ah)anthracene. Additionally, inorganic contaminants that may have been emitted in these eruptive processes will be quantified in whole blood. This includes: (a) trace elements (Co, Cr, Cu, Fe, Mn, Ni, Se, and Zn); (b) toxic elements listed in the Agency for Toxic Substances and Disease Registry (ATSDR) inventory, such as Ag, Al, As, Be, Cd, Hg, Pb, Pd, Sb, Sr, Th, Ti, Tl, U, and V; (c) rare earth elements and other minor elements (Au, Bi, Ce, Dy, Eu, Er, Ga, Gd, Ho, In, La, Lu, Nb, Nd, Os, Pr, Pt, Ru, Sm, Sn, Tb, Ta, Tm, Y, and Yb). All these analyses will be performed using gas chromatography coupled with triple quadrupole mass spectrometry (GC-MS/MS) for organic contaminants and inductively coupled plasma mass spectrometry (ICP-MS) for inorganic contaminants. The determinations will be carried out in the Toxicology laboratories of the two public universities of the Canary Islands.

Additionally, each participant underwent forced spirometry to measure lung function following the recommendations by the American Thoracic Society and the European Respiratory Society during the current SARS-CoV-2 pandemic. Forced expiratory volume in the first second (FEV-1) and forced vital capacity (FVC), among other parameters, were measured. Spirometry tests were conducted using a portable spirometer acquired specifically for this study (Sibelmed, model Datospir Touch 3000).

Statistical analysis

Categorical variables will be presented with their distribution of absolute and relative frequencies. Quantitative variables that follow a normal distribution will be summarized using the mean and standard deviation (± SD), while those that do not follow this distribution will be presented with the median and interquartile range (IQR). To calculate the distance to the volcano, participant home coordinates during the eruption were obtained using the geodist command in STATA, and elevation was obtained using the elevatr Statistical package in R.

A comparison of the distribution of sociodemographic characteristics, variables related to the level of exposure during the eruption and previous comorbidities of the participants in the general population between the two regions (west and east) was performed. For categorical variables, the Chi-square test were used. Comparisons of means between two regions were performed by Student’s t-test if the variables followed a normal distribution, or by the nonparametric Mann-Whitney U test for asymmetric variables. Finally, multivariate logistic regression models were performed to evaluate the independent effect of the place of residence (west vs. east) on acute symptomatology during the eruption. Those variables considered to be of interest were introduced as adjustment variables. The crude and adjusted odds ratios (OR) are presented together with their 95% confidence intervals (CI). Statistical significance was assumed as p < 0.05. Analyses were performed using the statistical package SPSS 26.0® (SPSS Inc., Chicago, IL, USA).

Results

Preliminary results of the descriptive analysis are presented for the first 1002 participants: 857 participants from the ISVOLCAN cohort, representing the general adult population of La Palma Island, and 145 intervening personnel who accessed the exclusion zone during the eruption.

Figure 1 shows the flowchart of the study sample. As of December 31, 2022, a total of 2355 phone calls were made to randomly selected individuals from the general population, and 857 participants were included (36.4% of those initially selected). In addition to the general population sample, the interveners (n = 145) were mainly composed of members of State Security Forces, Emergency Services, and cleaning workers.

Fig. 1
figure 1

Flowchart of the ISVOLCAN study cohort until December 31, 2022

Table 1 describes the sociodemographic characteristics of the analysed sample from the general population. The mean age was 50.8 years (± 16.4), with a higher proportion of females. The majority had secondary education, and 20.8% of the sample were unemployed before the eruption; a similar situation was found in the two regions. During the eruption 662 (77.2%) resided in the western region and 198 (22.8%) in the eastern region. The group of participants from the western region presented a higher percentage of women and a higher percentage of unemployed people significantly.

Table 1 Sociodemographic characteristics of the general population according to the place of residence during the eruption (western region/eastern region)

In the interveners, the mean age was slightly younger (45.7 years (± 11.8)) with a predominance of males (supplemental Table 1).

Figure 2 shows the geolocation of the ISVOLCAN cohort based on the coordinates of participants addresses before and during the eruption in the general population. It can be observed that during the eruption, there was a displacement of residents from the Valle de Aridane area to other parts of the island.

Fig. 2
figure 2

Place of residence of ISVOLCAN cohort participants: (a) before and (b) during Tajogaite volcano eruption

Characteristics related to exposure during the volcanic eruption in general population are described in Table 2. The median distance from participants residence to the volcano during the eruption was 7.1 km (IQR:6.1–9.3); for the western region, it was 6.7 km (IQR: 4.9–7.3), while for the eastern region, it was 10.9 km (IQR: 9.3–12.7). Most of the population in the sample engaged in cleaning up volcanic ash, both inside and outside their homes, using tools with a high capacity for particle projection, such as brooms and blowers. During the eruption, 85% of the general population always used masks when outdoors, with FFP2 masks being the most used. In the bivariate analysis, it found that the location of the usual residence with less distance to the volcano and higher altitude, a more frequent cleaning of volcanic ash both inside and outside the homes and a more daily hours spent in outdoor environments were registered between participants from western region compared to those from the eastern one. The frequency of use of face masks and protective eyeglasses in outdoor environments did not differ between the two regions.

Table 2 Aspects related to the level of exposure on the general population according to the place of residence during the eruption (western region/eastern region)

The intervining group showed a similar distribution to the general population regarding variables related to volcanic ash cleaning (location, tools used, and cleaning frequency), as well as mask usage frequency and type in outdoor environments (supplemental Table 2).

Table 3 shows baseline characteristics in general population related to lifestyle and pre-existing comorbidities before the eruption and use of healthcare resources and acute symptoms reported by participants during the eruption. The most prevalent pre-eruption comorbidities included blood hypertension (24.3%), depression and anxiety. The most frequently reported acute symptoms by the general population were eye irritation (45.9%), insomnia (44.9%), anxiety and depression (44.7%), and respiratory symptoms. In addition, 12.1% of the sample reported having an emergency visit at a hospital or primary care centres. The main reason for primary care visits was anxiety or depression, while hospital emergency visits were mainly due to osteomuscular traumas. Only 1.8% of the participants reported being hospitalized during the eruption, with surgical intervention being the primary reason. Participants from the western region compared to those from the eastern one were, significantly, more current smokers. Regarding acute symptomatology, western participants showed, in a statistically significant way, higher prevalence of nausea and vomiting, headache, lower respiratory tract symptoms (cough, dyspnea or wheezing), chest pain, insomnia, depression and anxiety, ocular, nasal and ear symptoms.

Table 3 Lifestyle factors, prevalence of previous comorbidities, use of healthcare resources and acute symptoms on the general population according to the place of residence during the eruption (western region/eastern region)

In the interveners, the acute symptoms reported during the eruption were like those of the general population, as well as the utilization of healthcare services. However, the percentage of hospitalizations was lower in this group (supplemental Table 3).

Table 4 shows the adjusted and unadjusted effect of region of residence during the volcano eruption (west/east) in general population on each of the most prevalent acute symptoms that showed statistically significant differences between the two regions. Age, gender, education level, employment, distance to the volcano, ash cleaning, type of cleaning tool used, daily hours in outdoor environments and type of smoker were entered as adjustment variables in all multivariate models. In addition, for the acute symptom lower respiratory tract symptoms, we adjusted for having suffered from any respiratory infection (influenza, COVID 19 or cold) during the months of the volcano eruption. Adjusted multivariate analysis results show that participants in the western region had a higher likelihood of lower respiratory tract symptoms (OR 1.99; 95% CI:1.33–2.99), depression and anxiety (OR 1.95; 95% CI:1.30–2.93) and insomnia (OR 2.03; 95% CI:1.33–3.09), compared to those in the eastern region.

Table 4 Association of place of residence on the general population (western vs. eastern region) in acute symptomatology during the eruption of the volcano. Logistic regression

Discussion

This article presents the methodology of the ISVOLCAN study, as well as a descriptive analysis of the baseline characteristics of the first 1002 participants (857 participants from the general adult population of La Palma Island, and 145 interveners, potentially highly exposed).

After the initial telephone contact was established with the selected individuals from the general population of the island, an initial response rate of 36.4% was observed. Although a higher participation rate was expected, the conditions of uncertainty and vulnerability experienced by the population immediately after the eruption was extinguished and during the subsequent months, generated certain limitations. At the beginning of the ISVOLCAN study, part of the evacuated population was still displaced or involved in bureaucratic and administrative procedures related to the disaster.

As mentioned previously, epidemiological data for each participant were collected through a health questionnaire. Analysis of this data revealed that the participants had a mean age within the working-age range, with a predominance of women and most individuals who had completed secondary education. The recruited population mainly resided in the municipalities affected by the volcano, with the highest number of displacements during the eruption occurring among the inhabitants of Los Llanos de Aridane, which coincided with the movement of the lava flows. Regarding the intervining group, it was observed that they were younger and predominantly male, reflecting the male dominance in certain professions related to the field of public safety.

Factors related to the level of exposure of the participants were also considered in the analysis. It was observed that the proximity to the volcano was about 7 km, even less for the residents of Valle de Aridane. This proximity is unusual compared to other volcanic phenomena documented in scientific literature. For instance, in the case of Holuhraun, population centres were located at least 100 km away from the volcano, with only a few isolated farms found at a closer distance, approximately 70 km [19]. Another recent example concerns the Nyragongo or Nyamulagira volcanoes in the Republic of Congo, which affected a population of nearly one million people around the volcano, at approximately 15–30 km [20]. Therefore, in La Palma Island, the local population resided much closer to the eruption at the time compared to other mentioned populations.

Various health risks associated with the size of PM and their potential environmental impact on agriculture and water reservoirs have been reported [4]. Indeed, the deposition of several heavy metals, such as chromium and arsenic, in soils near volcanic eruptions has been documented, both of which have carcinogenic effects at certain levels [19]. In line with this, a very recent publication shows the chemical characterization of ash samples from Tajogaite eruption, founding that the most of the water-soluble compounds were SO4, F, Cl, Na, Ca, Ba, Mg and Zn; worryingly, the authors conclude that F and Cl concentration may exceed both the recommended levels for irrigation purpose and for health [21].

Moreover, the size of PM is of critical importance; particles smaller than 10 μm (PM10) can penetrate and reach the alveolar region of the lungs [3], while those smaller than 2.5 μm (PM2.5) may even cross the lung barrier and enter the bloodstream. There is an extensive body of evidence in relation to the health effects of the long-term exposure to PM 2.5 or lesser. The main reported effects are on all-causes and cause-specific mortality [22], incidence of cardiovascular or respiratory diseases [19, 23], incidence of endocrine and metabolic disorders such as type 2 diabetes [24] and incidence of lung cancer among others, even at concentrations below current EU limit values and possibly WHO Air Quality Guidelines [25]. However, to the best of our knowledge, there are no published studies that analyze the potential effects of degassing exposure on the population of La Palma, neither before nor during the eruption.

During the Tajogaite eruption, daily air quality monitoring was carried out through eight stations located in different points of Valle de Aridane and the eastern region of the island. Based on these records, the average levels of SO2 concentration in the island were recently published, and it was observed that the threshold recommended as safe by the European Commission was exceeded in the Valle area during 1 to 4% of the eruption duration. Furthermore, during the first month of the eruption, the threshold of 400 μm-3 was frequently exceeded, especially in the later stages of the phenomenon, in contrast to the emissions of particulate matter [2, 26].

It is noteworthy to mention that, due to the recommendations of authorities and scientists, as well as the activation of volcanic emergency protocols, the integrity of the population was successfully safeguarded. However, it is reasonable to assume that the displacements of the evacuated population during the eruption could have had an impact on their health. Throughout the volcanic event, the island’s population received daily information about the necessary preventive measures in each municipality, based on air quality and the evolution of volcanic ash. In the case of our sample from the general population, 15.4% were evacuated during the eruption, and less than half of the evacuated individuals returned to their usual homes after an average of approximately 3 months.

On the other hand, exposure to volcanic gases and ash has been widely associated with increased respiratory morbidity and short-term irritation in the respiratory tract, ocular mucosa, and skin due to their chemical and mechanical irritant effects [3, 19, 27]. In the case of ISVOLCAN cohort participants, ocular and upper respiratory tract irritation were the most frequent acute symptoms. These findings are consistent with epidemiological studies conducted in the general population, both during the acute phase [28] and 6–9 months after exposure [11], as well as in highly exposed professionals [29]. Other studies evaluating the reasons and number of visits to hospital emergency departments have detected an increase in visits due to respiratory diseases and ocular disorders [30, 31].

During the volcanic eruption, a significant proportion of the participants carried out ash cleaning tasks both indoors and outdoors, thereby increasing their exposure to the emitted material. As the eruption coincided with the second year of the SARS-CoV-2 pandemic, the population already had access to masks and was used to wearing them; the majority of the participants stated using masks when outdoors, with FPP2 masks being the most commonly used in these environments during the eruption, as they have demonstrated effectiveness in protecting against the inhalation of volcanic ash [32]. Certainly, it is imperative to maintain a surveillance over this excessive exposure in the coming years to comprehensively gauge potential medium and long-term repercussions. In the aftermath of the Tajogaite volcanic eruption, numerous supplementary investigations have been instigated, in addition to ISVOLCAN, with the aim of enhancing the monitoring of the health of the local population. Notably, the ASHES study is among these initiatives, with its principal focus being the assessment of respiratory health outcomes associated with exposure to volcanic emissions [33].

Moreover, prior investigations following volcanic eruptions have demonstrated a notable rise in the occurrence of psychiatric disorders within the general population [34]. Evacuated individuals, in particular, exhibited a pronounced prevalence of post-traumatic stress and depressive symptoms [35]. During the eruption period, nearly half of the individuals reported insomnia and symptoms indicative of mood disorders, such as anxiety or depression. Notably, those who had to undergo evacuation displayed a higher incidence of these symptoms. The eruption caused significant disruptions in the daily routines of the population in specific municipalities, especially those directly affected by evacuation orders.

The elevated prevalence of anxiety and depression can be related to several factors, including increased work demands during the eruption and the uncertainty concerning personal health, the well-being of others, property, and crop security, as well as the outlook for the future. Furthermore, given the substantial number of seismic events and the explosive nature of the eruption, it is plausible that these anxiety-related symptoms contributed to the substantial percentage of reported insomnia among the affected population.

Adjusted multivariate analysis results show that participants in the western region compared to those in the eastern region had a higher likelihood of lower respiratory tract symptoms, depression and anxiety, and insomnia. These results are similar to those found in the few epidemiological studies conducted in the general population that evaluate symptomatology, acute or short-term, during the eruption according to the level of exposure. These results are in concordance to previous evidence [11, 36].

Furthermore, the recognition of volcanic eruptions as sources of toxic elements underscores the environmental exposure faced by populations residing in close proximity to these emission sites. Environmental studies conducted worldwide, including the Canary Islands, have consistently identified volcanic eruptions as significant contributors of inorganic elements known to be toxic to humans, such as Se, Cd, Pb and Hg [37, 38]. Notably, recent findings from the ISVOLCAN study have documented elevated levels of Fe, Al, Ti, V, Ba, Pb, Mo, Co, and Rare Earths in banana crops on the island during the eruption period [39].

However, studies focused on monitoring toxin levels in populations affected by eruptions are limited, primarily due to the challenge of simultaneously quantifying these inorganic toxins in blood samples collected from affected individuals. Furthermore, the necessary analytical methods are mostly expensive, limiting their inclusion into epidemiological studies. In this context, our research team, as experts in toxicological analysis of both major inorganic and organic pollutants, is presently conducting determinations using venous blood samples from study participants, although results are pending.

The main limitation of the ISVOLCAN study, as is common in cohort studies, is its high cost, which is exacerbated in our case by logistical difficulties inherent in a fragmented territory like the Canary Islands, limiting the transfer of biological samples and human or material resources between islands. Additionally, while the participants were randomly selected from the general population, there may exist a selection bias if those who chose not to participate had some differential characteristics (e.g., older age, pre-existing health issues, etc.) compared to the participants, which could limit the detection of certain relevant associations.

Furthermore, the epidemiological data relies on self-reporting by the participants, which could introduce information biases affecting the validity of the results. Additionally, the high percentage of losses during follow-up, related to this type of design, could generate a survival bias. To address these concerns, several methodological strategies have been implemented. The sample size was increased to more than double the initial estimate, that is why recruitment and inclusion of participants are ongoing at this moment. Moreover, as a strength of the study, data collection started as soon as possible after the eruption was finished, carried out by personnel specially trained to ensure rigor and thoroughness in the process, following the recommendations of the IVHHN regarding epidemiological data records for such phenomena. Additionally, prior to analysis, the data undergo rigorous quality control and verification processes.

Given that the data come from a randomized sample of the general population of the island, followed over several years, this study will allow for the detection of causal associations. It is worth noting that the inclusion of interveners in the ISVOLCAN cohort provides a significant area of study since they can be considered as individuals with high prior exposure.

Conclusion

The ISVOLCAN study has been meticulously designed as a 10-year follow-up study aimed at assessing the medium to long-term health impact on the adult general population of La Palma Island following the recent eruption of the Tajogaite volcano. Despite currently being in a recruitment phase, the study has successfully completed several stages of biological sample collection and biomedical data gathering. Once the baseline measurements are finalized and toxicological determinations are conducted, data from over 2000 individuals with varying levels of exposure during the eruption are expected to be obtained. Lastly, in our knowledge this study is the first to publish data related to the short-term health impact on the population of La Palma following the eruption of the Tajogaite volcano.

Data availability

The data that support the findings of this study are not openly available due to reasons of sensitivity and are available from the corresponding author upon reasonable request.

Abbreviations

CO:

Carbon monoxide

CO2 :

Carbon dioxide

SO2 :

Sulfur dioxide

HCl:

Hydrogen chloride

HF:

Hydrogen fluoride

H2S:

Hydrogen sulfide

PM:

Particulate matter

IVHHN:

International Volcanic Health Hazard Network

WHO:

World Health Organization

PAHs:

Primarily polycyclic aromatic hydrocarbons

ATSDR:

Agency for Toxic Substances and Disease Registry

SD:

Standard deviation

IQR:

Interquartile range

PM2.5 :

Particulate matter ≤ 2.5 μm

PM10 :

Particulate matter < 10 μm

References

  1. Freire S, Florczyk AJ, Pesaresi M, Sliuzas R. An Improved Global Analysis of Population Distribution in proximity to active volcanoes, 1975–2015. ISPRS Int J Geo-Inf. 2019;8(8):341.

    Article  Google Scholar 

  2. Milford C, Torres C, Vilches J, Gossman AK, Weis F, Suárez-Molina D, et al. Impact of the 2021 La Palma volcanic eruption on air quality: insights from a multidisciplinary approach. Sci Total Environ. 2023;869:161652.

    Article  PubMed  ADS  CAS  Google Scholar 

  3. Horwell CJ, Baxter PJ. The respiratory health hazards of volcanic ash: a review for volcanic risk mitigation. Bull Volcanol. 2006;69(1):1–24.

    Article  ADS  Google Scholar 

  4. Hansell AL, Horwell CJ, Oppenheimer C. The health hazards of volcanoes and geothermal areas. Occup Environ Med. 2006;63(2):149–56.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  5. Padrón E, Pérez N, Rodríguez F, Melián G, Hernández PA, Sumino H. Dynamics of diffuse carbon dioxide emissions from Cumbre Vieja volcano, La Palma, Canary Islands. Bull Volcanol. 2015;77:28.

    Article  ADS  Google Scholar 

  6. Zappulla D. Environmental stress, erythrocyte dysfunctions, inflammation, and metabolic syndrome: adaptations to CO2 increases? J Cardiometab Syndr. 2008;3(1):30–4.

    Article  PubMed  Google Scholar 

  7. Stewart C, Damby DE, Horwell CJ, Elias T, Ilyinskaya E, Tomašek I, et al. Volcanic air pollution and human health: recent advances and future directions. Bull Volcanol. 2021;84(1):11.

    Article  ADS  Google Scholar 

  8. Gudmundsson G. Respiratory health effects of volcanic ash with special reference to Iceland. A review. Clin Respir J. 2011;5(1):2–9.

    Article  PubMed  Google Scholar 

  9. Higuchi K, Koriyama C, Akiba S. Increased mortality of respiratory diseases, including lung cancer, in the area with large amount of ashfall from Mount Sakurajima volcano. J Environ Public Health. 2012;2012:257831.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Longo BM, Yang W. Acute bronchitis and volcanic air pollution: a community-based cohort study at Kilauea Volcano, Hawai’i, USA. J Toxicol Environ Health A. 2008;71(24):1565–71.

    Article  PubMed  CAS  Google Scholar 

  11. Carlsen HK, Hauksdottir A, Valdimarsdottir UA, Gíslason T, Einarsdottir G, Runolfsson H, et al. Health effects following the Eyjafjallajökull volcanic eruption: a cohort study. BMJ Open. 2012;2(6):e001851.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Oudin A, Carlsen HK, Forsberg B, Johansson C. Volcanic Ash and Daily Mortality in Sweden after the Icelandic volcano eruption of May 2011. Int J Environ Res Public Health. 2013;10(12):6909–19.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Instituto Geográfico Nacional. Instituto Geográfico Nacional - Servicio de Documentación. 2022 [cited 2023 Jul 26]. La Palma (Isla). Mapas topográficos. 1996. Available from: https://www.ign.es/web/catalogo-cartoteca/resources/html/017034.html.

  14. Instituto Canario de Estadística ISTAC. Official Population Figures Referring to Revision of Municipal Register 1 January, 2022. Canary Islands Statistic Institute, Spain. [Internet]. 2022 [cited 2023 Sep 17]. Available from: http://www.gobiernodecanarias.org/istac/estadisticas/.

  15. Bonadonna C, Pistolesi M, Dominguez L, Freret-Lorgeril V, Rossi E, Fries A, et al. Tephra sedimentation and grainsize associated with pulsatory activity: the 2021 Tajogaite eruption of Cumbre Vieja (La Palma, Canary Islands, Spain). Front Earth Sci. 2023;11:1166073.

    Article  ADS  Google Scholar 

  16. Gobierno de Canarias. PEVOLCA - Comité científico [Internet]. 2021 [cited 2023 Jul 28]. Available from: https://www.gobiernodecanarias.org/infovolcanlapalma/pevolca/.

  17. The International Volcanic Health Hazard Network (IVHHN). Standardized Epidemiological Study Protocol to Assess Short-term Respiratory and Other Health Impacts in Volcanic Eruptions [Internet]. 2019 [cited 2023 Jul 19]. Available from: https://www.ivhhn.org/uploads/IVHHN%20Basic%20Protocol.pdf.

  18. Bobes J, Calcedo-Barba A, García M, François M, Rico-Villademoros F, González MP, et al. [Evaluation of the psychometric properties of the Spanish version of 5 questionnaires for the evaluation of post-traumatic stress syndrome]. Actas Esp Psiquiatr. 2000;28(4):207–18.

    PubMed  CAS  Google Scholar 

  19. Carlsen HK, Ilyinskaya E, Baxter PJ, Schmidt A, Thorsteinsson T, Pfeffer MA, et al. Increased respiratory morbidity associated with exposure to a mature volcanic plume from a large Icelandic fissure eruption. Nat Commun. 2021;12(1):2161.

    Article  PubMed  PubMed Central  ADS  CAS  Google Scholar 

  20. Michellier C, Katoto P, de Dramaix MC, Nemery M, Kervyn B. Respiratory health and eruptions of the Nyiragongo and Nyamulagira volcanoes in the Democratic Republic of Congo: a time-series analysis. Environ Health. 2020;19(1):62.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. Ruggieri F, Forte G, Bocca B, Casentini B, Bruna Petrangeli A, Salatino A, Gimeno D. Potentially harmful elements released by volcanic ash of the 2021 Tajogaite eruption (Cumbre Vieja, La Palma Island, Spain): implications for human health. Sci Total Environ. 2023;905:167103.

    Article  PubMed  ADS  Google Scholar 

  22. Chen J, Hoek G. Long-term exposure to PM and all-cause and cause-specific mortality: a systematic review and meta-analysis. Environ Int. 2020;143:105974.

    Article  PubMed  CAS  Google Scholar 

  23. Wolf K, Hoffmann B, Andersen ZJ, Atkinson RW, Bauwelinck M, Bellander T, et al. Long-term exposure to low-level ambient air pollution and incidence of stroke and coronary heart disease: a pooled analysis of six European cohorts within the ELAPSE project. Lancet Planet Health. 2021;5(9):e620–32.

    Article  PubMed  Google Scholar 

  24. Yang Z, Mahendran R, Yu P, Xu R, Yu W, Godellawattage S, Li S, Guo Y. Health effects of Long-Term exposure to ambient PM2.5 in Asia-Pacific: a systematic review of Cohort studies. Curr Environ Health Rep. 2022;9(2):130–51.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Hvidtfeldt UA, Severi G, Andersen ZJ, Atkinson R, Bauwelinck M, Bellander T. Long-term low-level ambient air pollution exposure and risk of lung cancer - A pooled analysis of 7 European cohorts. Environ Int. 2021;146:106249.

    Article  PubMed  CAS  Google Scholar 

  26. Filonchyk M, Peterson MP, Gusev A, Hu F, Yan H, Zhou L. Measuring air pollution from the 2021 Canary Islands volcanic eruption. Sci Total Environ. 2022;849:157827.

    Article  PubMed  ADS  CAS  Google Scholar 

  27. Hansell A, Oppenheimer C. Health hazards from volcanic gases: a systematic literature review. Arch Environ Health. 2004;59(12):628–39.

    Article  PubMed  CAS  Google Scholar 

  28. Carlsen HK, Gislason T, Benediktsdottir B, Kolbeinsson TB, Hauksdottir A, Thorsteinsson T, et al. A survey of early health effects of the Eyjafjallajokull 2010 eruption in Iceland: a population-based study. BMJ Open. 2012;2(2):e000343.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Carlsen HK, Aspelund T, Briem H, Gislason T, Jóhannsson T, Valdimarsdóttir U, et al. Respiratory health among professionals exposed to extreme SO2 levels from a volcanic eruption. Scand J Work Environ Health. 2019;45(3):312–5.

    Article  PubMed  Google Scholar 

  30. Baxter PJ, Ing R, Falk H, French J, Stein GF, Bernstein RS, et al. Mount St Helens eruptions, May 18 to June 12, 1980. An overview of the acute health impact. JAMA. 1981;246(22):2585–9.

    Article  PubMed  CAS  Google Scholar 

  31. Lombardo D, Ciancio N, Campisi R, Di Maria A, Bivona L, Poletti V, et al. A retrospective study on acute health effects due to volcanic ash exposure during the eruption of Mount Etna (Sicily) in 2002. Multidiscip Respir Med. 2013;8(1):51.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Steinle S, Sleeuwenhoek A, Mueller W, Horwell CJ, Apsley A, Davis A, et al. The effectiveness of respiratory protection worn by communities to protect from volcanic ash inhalation. Part II: total inward leakage tests. Int J Hyg Environ Health. 2018;221(6):977–84.

    Article  PubMed  Google Scholar 

  33. Ruano-Ravina A, Acosta O, Díaz Pérez D, Casanova C, Velasco V, Peces-Barba G, et al. A longitudinal and multidesign epidemiological study to analyze the effect of the volcanic eruption of Tajogaite volcano (La Palma, Canary Islands). The ASHES study protocol. Environ Res. 2023;216(Pt 2):114486.

    Article  PubMed  CAS  Google Scholar 

  34. Shore JH, Tatum EL, Vollmer WM. Evaluation of mental effects of disaster, Mount St. Helens eruption. Am J Public Health. 1986;76(3 Suppl):76–83.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  35. Goto T, Wilson JP, Kahana B, Slane S. The Miyake Island Volcano Disaster in Japan: loss, uncertainty, and Relocation as predictors of PTSD and Depression. J App Soc Psychol. 2006;36(8):2001–26.

    Article  Google Scholar 

  36. Shimizu Y, Dobashi K, Hisada T, Ono A, Todokoro M, Iijima H, et al. Acute impact of volcanic ash on asthma symptoms and treatment. Int J Immunopathol Pharmacol. 2007;20(2 Suppl 2):9–14.

    Article  PubMed  CAS  Google Scholar 

  37. Ilyinskaya E, Mason E, Wieser PE, Holland L, Liu EJ, Mather TA, et al. Rapid metal pollutant deposition from the volcanic plume of Kīlauea, Hawai’i. Commun Earth Environ. 2021;2(1):1–15.

    Google Scholar 

  38. Rodríguez Martín JA, Nanos N, Miranda JC, Carbonell G, Gil L. Volcanic mercury in Pinus canariensis. Naturwissenschaften. 2013;100(8):739–47.

    Article  PubMed  ADS  Google Scholar 

  39. Rodríguez-Hernández A, Díaz-Díaz R, Zumbado M, Bernal-Suárez M, del Acosta-Dacal M, Macías-Montes A. Impact of chemical elements released by the volcanic eruption of La Palma (Canary Islands, Spain) on banana agriculture and European consumers. Chemosphere. 2022;293:133508.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

The authors give special thanks to Marta Rodríguez Pérez for her invaluable contribution in the technical support of this study. She has managed the contacts with selected persons and scheduled participants. She also designed all the electronic documents to recording data from the participants and has carried out the data quality control of the database of the cohort. Many thanks too to the Primary Care health staff of La Palma, nurses and family physicians, around the island, for their help to disseminate the study and to administer epidemiological questionnaires to the participants.

Funding

Fundación Canaria Instituto de Investigación Sanitaria de Canarias (FIISC: ST22/07). Instituto de Salud Carlos III, Madrid (PI22/00395).

Author information

Authors and Affiliations

Authors

Contributions

Conception and design: M.C.R.P, L.D.B, A.H.T, A.C.L, I.G.T. Study recruitment and sample processing: A.D.A.P, M.C.D.A, J.F.F.J, L.V.G. Samples management and analysis: I.G.T, L.D.B, A.H.T, L.V.G. Acquisition of epidemiological data: M.C.R.P, A.D.A.P, M.C.D.A, J.F.F.J. Complete data curation and analysed: M.C.R.P, M.E.F.F. Interpretation of the data: M.C.R.P, M.E.F.F, A.C.L, L.D.B, K.S.R. Draft the article: M.C.R.P, M.E.F.F, A.C.L, L.D.B, K.S.R. All authors revised and approved the final manuscript.

Corresponding author

Correspondence to María Cristo Rodríguez-Pérez.

Ethics declarations

Ethics approval and consent to participate

The study has obtained authorization from the health authorities and received a favorable decision from the Provincial Ethics and Medicines Committee (ref. CHUNSC_2021_88). Participants were required to provide written consent before being included in the study.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Rodríguez-Pérez, M.C., Ferrer, M.E.F., Boada, L.D. et al. Health impact of the Tajogaite volcano eruption in La Palma population (ISVOLCAN study): rationale, design, and preliminary results from the first 1002 participants. Environ Health 23, 19 (2024). https://doi.org/10.1186/s12940-024-01056-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12940-024-01056-4

Keywords