Mortality and morbidity assessment attributed to short- and long-term exposure to fine particles in ambient air of Agadir city, Morocco: The AirQ model approach

Article information

Environ Anal Health Toxicol. 2023;38.e2023009
Publication date (electronic) : 2023 May 12
doi :
1Laboratory of Geosciences, Environment and Geomatics, Faculty of Sciences, Ibn Zohr University, Agadir, Morocco
2High Institute of Nursing Professions and Health Techniques of Agadir, Agadir, Morocco
3High Institute of Nursing Professions and Health Techniques of Agadir, Health Sciences and Environment Laboratory, Health Sciences, Epidemiology and Human Pathologies Research Team (ER-2SEPH), Agadir, Morocco
4Research Laboratory of Innovation in Health Sciences (LARISS), Faculty of Medicine and Pharmacy, Ibn Zohr University, Agadir, Morocco
5Laboratory of Cell Biology and Molecular Genetics, Faculty of Sciences, Ibn Zohr University, Agadir, Morocco
6High Institute of Nursing Professions and Health Techniques, Marrakech, Morocco
7Health Sciences Research Laboratory, Faculty of Medicine and Pharmacy, Ibn Zohr University, Agadir, Morocco
8Clinical Epidemiology and Medico-Surgical Sciences, Faculty of Medicine and Pharmacy, Mohammed V University, Rabat, Morocco
Recommended by: Prof. Sanghyuk Bae
Received 2023 January 22; Accepted 2023 April 18.


It is well established that respiratory mortality and morbidity are associated with high concentrations of fine particles such as PM2.5. The aim of this study was to evaluate the long- and short-term impacts of PM2.5 on the population of Agadir, Morocco, using AirQ 2.1.1 software. The mean PM2.5 values were obtained from data collected at three sites. Baseline incidence data were obtained from the literature, and relative risk (RR) values were referenced from the World Health Organization. This study quantified long-term total mortality (LT-TM), lung cancer mortality (LT-LC), morbidity from acute lower respiratory tract infections (LT-ALRI), and morbidity from chronic obstructive pulmonary disease (LT-COPD), as well as short-term total mortality (ST-TM). The attributable proportions (AP) of LT-TM and LT-LC were estimated to 14.19% and 18.42%, respectively. Their excess deaths were estimated to 279 and 11 persons, respectively, and their RRs to 1.16 (95% CI: 1.10-1.22) and 1.23 (95% CI: 1.12-1.37), respectively. Furthermore, the AP of LT-ALRI and LT-COPD were estimated to 14.36% and 15.68%, respectively, their excess deaths to 33 and 4, and their RRs to 1.17 (95% CI: 1.11-1.31) and 1.19 (95% CI: 1.00-1.02), respectively. In comparison, the AP of ST-TM was estimated to 1.27%, with a 25-person excess death rate. This study was conducted to inform decision-making and to promote local policies on ambient air quality.


The adverse effects of air pollution on health are well established and documented. Currently, about 7 million premature deaths per year are attributable to air pollution [1], and only one in ten people lives in communities that meet World Health Organization air quality guidelines (WHO AQG). Fine particulate matters are among frequently studied air pollutants that constitute major risk factors for comorbidities and have wide range of identified adverse health effects. They are defined as particles less than 2.5 µ m in aerodynamic diameter (PM2.5) [2], and could result from various anthropogenic or natural activities [3].

To date, strong evidence shows causal associations between exposure to PM2.5 air pollution and all-cause mortality, as well as mortality and morbidity from acute lower respiratory tract infections (ALRI), chronic obstructive pulmonary disease (COPD), ischemic heart disease (IHD), lung cancer (LC), and stroke [4-6]. This causality is determined in the long term for all-cause, cardiovascular and respiratory mortality [7,8]. Assuming a linear association, Chen and Hoek (2020) reported a relative risk estimate of 1.08 (95% CI: 1.06-1.09) per 10 µ g/m³ of PM2.5 for nonaccidental mortality [9]. Both authors evoked evidence of a supra-linear relationship implying a greater increase in risk at lower levels of exposure. Long-term effects investigated by cohort studies, revealed also associations between PM2.5 exposure and total, cardiovascular, and lung cancer mortality [10-13]. Sensitive groups, including the elderly, children, and patients with respiratory and cardiovascular diseases exposed to even low concentrations of air pollutants, are likely to be affected [14-17].

The finest particles are more dangerous to human health than the coarsest ones [18]. Particulate matter, primarily from combustion, is deposited throughout the respiratory tract, particularly in the pulmonary alveoli. These particles include a variety of immunogenic substances linked to asthma exacerbation symptom [19]. Particles mostly from combustion, are deposited throughout the respiratory tract in the pulmonary alveoli, and include a variety of immunogenic substances related to asthma symptoms exacerbations [20]. In addition, exposure to particles may increase the risk of neurogenerative and neurodevelopmental disorders [21, 22].

Ambient PM2.5 concentrations vary to a great extent between and within countries, and more than 90% of the world's population in 2019 lived in areas where concentrations exceeded the 2005 WHO air quality guideline (WHO AQG) of 10 µ g/m3 [23]. Likewise, low- and middle-income countries are expected to account for 91% of PM2.5-related mortality in outdoor and indoor air, with a per capita rate three times higher than the one of high-income countries [24]. Furthermore, in their 2017 study on premature death from changes in air pollution attributable to climate change, Silva et al. (2017) reported that PM2.5 would cause 55,600 (-34,300 to 164,000) deaths in 2030 and 215,000 (-76,100 to 595,000) in 2100 (countering by 16% the global decrease in PM2.5-related mortality). If some major pollutants, such as PM2.5 and NO2, persist at present levels, air pollution is expected to cause around 2.5 million cases of noncommunicable diseases also by 2035 [25].

In 2021, the PM2.5 annual WHO AQG level has been lowered from 10 µ g/m3 to 5 µ g/m3. The 24-hour WHO AQG level for PM2.5 changed from 25 µ g/m3 to 15 µ g/m3.

Based on the WHO AQG, Morocco's air quality is moderately affected. In 2014, the economic cost of air pollution in Morocco represented 1.05% of the country's GDP [26]. Moreover, Moroccan cities suffer from high particle air pollution as a result of vehicle traffic and industry, as well as considerable terrigenous inputs related to the aridity of the climate and proximity to the desert. However, because monitoring networks are still in their early stages, data on these cities is limited. Cities in the South are now confronted with pollution from industrial sources as well as urban traffic of polluting vehicles as a result of late industrialization. According to the most recent WHO statistics from 2016, the yearly average concentration of PM2.5 in Morocco is 28.38 µ g/m3, ranging from 23.6 to 35.6 µ g/m3 [27].

So far, only five studies in Morocco have been conducted to quantify the health impact of particulate matter using, in addition to PM10, SO2, NOx, and O3 as exposure markers [28-32].

Despite the adverse health effects associated with ambient exposure to PM2.5, limited information is published on the spatial and temporal variation of PM2.5 concentrations in Morocco, and their health effects in terms of mortality and morbidity, particularly in a metropolis like Agadir City. The objective of this study was thus, to quantify the long-term mortality and short-term morbidities associated to ambient PM2.5 exposure in this area.

Materials and Methods

Study Area

Agadir city is located in southwest Morocco on the Atlantic coast. Its latitudinal and longitudinal coordinates, as well as altitude are respectively, 30°25′12′′ north, 9°35′53′′ west, and 31 meters above the sea level (Figure 1). The city is also the administrative capital of the Souss-Massa Region. According to the General Census of Population and Habitat in 2014 (GCPH 2014), Agadir city had a population of 421,844 persons and constitutes Morocco's first tourist destination and the country's first fishing harbor. Nevertheless, there is currently little evidence on the effects of air pollution on Agadir population. In 2020, the public transport system will have 43 lines and a fleet of 206 vehicles. The system connects the whole of the small communities, while peri-urban lines connect the agglomeration to localities up to 50 kilometers away. The annual rainfall is 250 mm. Between November and March, there are just a few days of rain. There are over 340 days of sunlight every year, however mist and morning dew are typical in the summer. The year-round trade wind front has a considerable impact on temperatures, which fluctuate little between winter and summer. In January, average temperatures vary from 14 to 16°C, while in July, they range from 19 to 22°C. However, upwellings of Saharan air occur on occasion, causing the temperature to climb beyond 40 °C for a few days. Agriculture, along with tourism, are the primary economic activities.

Figure 1.

Location map of the study area: (a) Morocco in Africa; (b) Study area in Agadir Ida Ou Tanan Prefecture; (c) Study area; Locations of the mobile air pollution measuring station in 2016 (green circles)

Air pollution data monitoring and exposure assessment

The average PM2.5 concentrations in 2016 were considered to be as indicators of short- and long-term exposures among Agadir inhabitants. These data were collected from the 2016 monthly report of provided by the Regional Committee for Monitoring and Surveillance of Air Quality in Souss-Massa Region. The automated PM10 sensor measures the amounts of PM suspended in the ambient air instantly. The measurement basic is beta ray absorption, and the PM collecting device is made up of filtering membranes 47 mm in diameter, which can be used for mass measurement control using the gravimetric technique and/or chemical tests.

In the Souss-Massa region, there is only one mobile monitoring unit for air pollutant concentrations. This station has been set for 10 months for the first location (April to November 2016) and 2 months for the second location (May and June 2016) (green circles in Figure 1).

Population and Health Data

The long-term health impact of PM2.5 is assessed considering all-cause death in adults aged 30 and older, as well as LC mortality in adults aged 25 and older. Yearly statistics obtained from the Souss-Massa Regional Hospital Center (SMRHC) on chronic obstructive pulmonary disease (COPD) cases in 2016 were used. Data on acute lower respiratory tract infections (ALRI) in children under the age of five were obtained from examinations of preschool and primary school students, during Ministry of Health campaigns conducted against transmissible ophthalmia in 2015 Health Report. Due to a lack of local data on the incidence of LC, standardized incidence of LC in 2004 was used as a reference from Casablanca Regional Cancer Registry. Finally, to compute baseline cumulative incidences of all-cause mortality and morbidity caused by ALRI and COPD, yearly statistics data from 2015 Health Report and SMRHC were used.

Statistical analysis using AirQ+ software

To assess the burden of disease resulting from air pollution and its effects, AirQ+ software developed by the WHO Regional Office for Europe was implemented [33]. AirQ+ includes methods for assessing the short- and long-term impacts of exposure to ambient air pollution. The main methods make use of scientific data from epidemiological cohort studies that indicate associations between long-term average concentration levels in polluted air and mortality risks in exposed populations.

AirQ+ estimations are based on the application of concentration-response functions, which relate exposure levels to a pollutant with a number of health events across the study area and period. These functions represent the relative health risk associated with a certain change in exposure level. They are developed from epidemiological studies and can be used when the causal nature of the association can be accepted.

The estimated proportion of attributable cases, estimated number of attributable cases, estimated number of attributable cases per 100,000 persons at risk, attributable proportion of cases in each air pollutant concentration category, cumulative distribution by air pollutant concentration, and life years lost are the main results generated by AirQ+.

To quantify long-term and short-term effects, the following information was provided: (1) annual average air quality data for long-term exposure effects, (2) population at-risk data, such as the total number of adults 30 years of age, (3) health data, such as baseline rates of health outcomes in the Greater Agadir population, (4) a threshold value to be considered (10 µ g/m3 as recommended by WHO), and (5) relative risk values (RRs).

For AirQ+, the relative risk due to air pollution is usually modeled with a log-linear function [34], this RR is calculated as RR = expβ×(X-X0), where X denotes the pollutant concentration (µ g/m3), and X0 denotes the threshold or counterfactual value, e.g: the background concentration or lowest possible value (µ g/m3). In this study x0 = 10 µ g/m3 for PM2.5 (annual mean not to be exceeded according to WHO AQG 2005). In the log-linear model, β indicates the change in RR for a one-unit change in concentration x.

The attributable proportion (AP) is the attributed fraction of the health effect caused by exposure in a given population over a given time interval, and it can be calculated as AP = 1 - 1/RR(c), where RR(c) is the relative risk for the health effect in the exposure category (c).

The number of cases per unit population (BE) can be calculated using BE = B × AP and a chosen baseline incidence (B) of the extreme health effect in the population. The number of attributable cases (AC) for a population of size N can be calculated as NE = N × AP [35].

AirQ+ incorporates default values for RRs, or concentration-response functions, for many health outcomes, as well as national conversion factors between PM2.5 and PM10 and solid fuel consumption data. Only deaths reported by the SMRHC that fulfilled the requirements in the AirQ+ software were considered.

Results and Discussion

During the study period, the annual average PM10 was 57.17 µ g/m3, and its concentration was 2.86 times higher than the 2005 WHO AQG (20 µ g/m3). The average PM2.5 concentration was obtained by converting the PM10 concentration according to a factor of 0.62 established by WHO in 2016. The average value of PM2.5 concentration after conversion was then estimated to 35.45 µ g/m3, and was 3.55 times higher than the 2005 WHO AQG (10 µ g/m3), as well as the national average concentration assessed by the WHO to 13.55 µ g/m3 (11.07-16.7 µ g/m3) [36]. The available data on annual mean PM10 concentrations are from 2016. The update of the WHO AQG occurred in 2021; therefore, the WHO AQG of 2005 was adopted.

The relative humidity percentage in Agadir in 2016 varied from 34.75 to 65.88%, with an average of 52.5%. The yearly averages for temperature and wind speed were 21.6 °C and 2.34 m/s, respectively.

Table 1 shows that in 2016 intra-hospital mortality was estimated to 1966 deaths. During the same year, 366 cases of ALRI in children under 5 years of age and 28 cases of COPD in adults 30 years of age were reported. The annual incidence per 100,000 persons of Total mortality, LC, ALRI, and COPD for 2016 was 968.48, 613.91, and 11.63, respectively.

Population characteristics.

Long-term PM2.5 ambient air pollution and adult mortality

Table 2 shows the long-term health effects caused by PM2.5 (total mortality, LC mortality, ALRI morbidity, and COPD morbidity) and the short-term health effects of PM2.5 (total mortality). AirQ+'s default RR for all-cause mortality is estimated to 1.062 (95% CI: 1.040-1.083). According to the results in Table 2, for an annual average PM2.5 concentration of 35.45 µg/m3, 279 deaths caused by long-term exposure to PM2.5 would be theoretically preventable out of the 1,966 deaths in all age groups in 2016. If the PM2.5 concentration did not exceed 10 µg/m3, the AP was 14.19%. The measured RR was estimated to be 1.16 (95% CI: 1.10-1.22). In this case, the estimated RR indicates that for an annual PM2.5 concentration that exceeds 10 µg/m3, the risk of death for 30-year-olds is 1.16 higher than for the other age groups.

PM2.5 health impact assessment results in AirQ+ in 2016.

Recent studies have found that PM has the biggest harmful impact on human health. Although various pollutants have been associated to health effects, PM2.5 has been given the most attention [2].

PM2.5 is recognized as one of the most important environmental health risk factors due to the possibility of particle penetration (deeper into the lung) [18].

Time series studies have revealed that pollution exposure increases risk. It investigates the daily association between air pollution and emergent health impacts such as cardiovascular and respiratory disease death or hospital admissions due to cardiovascular and respiratory disease. In this research, the RR value for the increase in overall mortality induced by PM2.5 was 0.16% for every 10 µg/m3 increment in pollutant concentration.

Given the rise in RR, it is reasonable to assume that PM2.5 had the greatest impact on human health by increasing its concentration. This conclusion was also confirmed in other countries exposed to different degrees and levels of pollution, such as China, the USA, and India [37-39]. In addition, recent studies conducted showed that several deaths have been attributed to long-term PM2.5 exposure [40-43].

Long-term PM2.5 ambient air pollution - integrated exposure-response functions (IER)

The IER function integrates RR information from studies on ambient air pollution, solid home cooking fuel, passive smoking, and active smoking [44]. These functions are used to estimate the disease burden associated with longterm exposure to PM2.5 in ambient air. AirQ+ is used to calculate the effects of PM2.5 on the basis of WHO air quality guidelines. Models such as log-linear or linear-log can be employed. Current IER functions have only been developed for specific adult mortality causes: IHD, stroke, COPD, and LC. Furthermore, IER functions for ALRI in 5-year-old children are estimated.

AirQ+ can also be used to assess changes in impact. Indeed, the IERs have been modified over the years to produce estimates of the Global Burden of Disease (GBD) [45]. In this regard, the counterfactual concentration range of 2.4 to 5.9µg/m3 was chosen to produce the 2015 GBD [6]. The values set by the 2005 WHO AQG were used to explore the benefits of reducing air pollution. Each exposure-response function refers to a specific age group.

Long-term ambient PM2.5 air pollution and adult LC mortality

The main question to be addressed here is: what are the expected health benefits (in terms of reduced LC cases) if the current level of air pollution measured by PM2.5 data, is reduced to the interim targets?

GBD 2015/2016 (integrated function 2016 with WHO interim target 1) was adopted as the calculation method with 2.4 µg/m3 as the threshold value for this model. The standardized incidence per 100,000 people is estimated for all regions in Morocco at 25.53. The population at risk is adults aged 25 years and older estimated at 240,825 according to GCPH 2014 [46].

The impact of PM2.5 exposure is assessed for the 2016 annual average of 35.45 µg/m3. Table 2 shows that lung cancer mortality is modestly decreased due to the expected low incidence, resulting in an estimated 11 avoided deaths in 2016. It's indeed important to note that the AP estimate for the current year is 18.42%.

Pope et al. (2011) calculated the adjusted associated risk of LC mortality for 10 µg/m3 of ambient PM2.5 to 1.14 (95% CI: 1.04-1.23) [47]. In a similar study, the associated risk of LC mortality from PM2.5 exposure was calculated to 1.14 (95% CI: 1.04, 1.23) [41]. Another study found that a 10 µg/m3 rise in PM2.5 was related with an increased incidence of LC (RR = 1.34; 95% CI: 1.10-1.65) [48]. Faridi et al (2018) found that deaths attributed to LC were 83 cases in 2015 [42].

According to Yarahmadi et al. (2018), the average annual number of LC deaths related to long-term PM2.5 exposure in Tehran was 142 cases [41]. Differences in baseline mortality could explain the differences in LC death values of this study compared to those above.

Long-term PM2.5 ambient air pollution and ALRI morbidity in adults

The goal is to estimate the expected health benefits in terms of reduced ALRI cases in children under the age of five if current levels of air pollution, as measured by PM2.5 data, are lowered to the interim targets. The analysis method used is based on the GBD 2015-2016 function (2016 integrated function with WHO 3 interim target), with 15 µg/m3 as the model's threshold value.

In 2015, the cumulative incidence per 1000 persons was estimated to be 613.91. In Table 2, the AP due to PM2.5 exposure over the model target value of 15 µg/m3 was estimated to 14.36% in 2016, with 33 cases of ALRI in children under the age of five potentially avoidable.

These findings are consistent with those of Sonego et al. (2015), who found that air pollution was significantly associated with an increased risk of death from ALRI [49], and Darrow et al. (2014), who reported that PM2.5 was associated with emergency department visits for pneumonia and upper respiratory tract infection in children under the age of five [50].

Long-term ambient PM2.5 air pollution and mortality in adults with COPD

The objective was to estimate the benefits of lowering the number of COPD cases in children and people aged 25 and older if present levels of air pollution, as assessed by PM2.5 data, are lowered to interim targets.

The GBD 2015-2016 function (a 2016 integrated function with the WHO 3 interim target) is used as the calculation method, considering 15 µg/m3 as the model's threshold value. In 2016, the cumulative incidence per 100,000 persons is estimated to 11.63.

The results of Table 2 indicate that in 2016, the AP due to exposure to PM2.5 above the target value taken by the model of 15 µg/m3 is estimated at 15.68%, while 4 of the 29 COPD cases reported by the pulmonary service of the SMRHC for adults aged 25 and over are potentially avoidable.

Short-term PM2.5 ambient air pollution and adult mortality

The WHO does not propose that short-term exposure to PM2.5 be quantified as a substitute to long-term exposure, although it is done for guidance. We are interested in quantifying the number of deaths from non-accidental causes that are attributed to short-term daily exposure to air pollution that exceeds the 2005 WHO AQG for PM2.5 (25 µg/m3).

AirQ+'s default RR for all-cause mortality was estimated to be 1.0123 (95% CI: 1.0045-1.0201). Calculations for mortality in adults aged 30 and older were conducted, using annual PM2.5 exposure averages of 35.45 µg/m3 in 2016. As predicted, short-term impacts are minor in comparison to long-term ones.

In fact, the results indicate that excess deaths would be potentially avoidable if the WHO recommended values for PM2.5 were met. It should be noted that the effects of short-term exposure to PM2.5 are largely included in the estimates of the effects of long-term exposure.

Table 2 shows that in 2016, with an annual average PM2.5 concentration of 35.45 µg/m3, 25 out of the 1,966 deaths across all age groups might have been avoided if the PM2.5 concentration did not exceed 25 µg/m3.

The AirQ software model has some limitations. For instance, it does not consider health impacts induced by exposure to combinations of so many pollutants, their toxicities or their interactions, instead of focusing on the effect of a single pollutant.

Also, the threshold values used to calculate RRs, are relative to studies conducted mainly in Europe and/or the Unites States. Finally, estimations generated by AirQ+ software are subject to uncertainty due to concentration-response functions based on restricted assumptions and levels of evidence.


To our knowledge, no research has been carried out to investigate the disease burdens associated with PM2.5 exposure in Agadir City, Morocco. Furthermore, in a recent systematic review on air quality in Morocco and its health impact on the Moroccan population conducted by Bouchriti et al. (2023), 1230 records were identified, of which 31 were eligible, all of which showed an exceedance of annual concentrations of air pollutants in relation to WHO guidelines [51]. Due to the reduced data from only two monitoring stations, the results of this study suggest that they should be interpreted with caution. Hence, the implementation of the AirQ+ software developed by the WHO, allowed the assessment and quantification of health effects of air pollution in this city for the first time. Based on the findings of the present study, and to help reduce air pollution, and consequently improve air quality, the following might be recommended. First, supporting research on the effects and consequences of air pollution on human health through the use of specific mobile sensors rather than the expensive fixed stations. Second, digitalizing health data in hospitals and primary health care centers to improve the information system and outcomes. Third, encouraging carpooling and developing a sustainable urban public transport system that minimizes greenhouse gas emissions. Finally, bringing national air quality standards up to date with the most recent WHO air quality recommendations.

Despite AirQ+ software limitations, the model provides a helpful and appropriate tool for leaders and health groups to boost collective awareness towards community environmental management and monitoring.


The authors express their gratitude to the person in charge of the monitoring and evaluation unit at the Regional Hospital Center Souss-Massa in Agadir city, as well as the Regional Committee for Monitoring and Surveillance of Air Quality, for providing the data for the study.


The authors declare that no funds, grants, or other support were received during the preparation of this manuscript. The authors have no relevant financial or non-financial interests to disclose.

CRediT author statement

YB: Conceptualization, Methodology, Formal analysis, Writing-Original draft preparation; AK: Validation, Visualization, Writing-Review & Editing; MA & AA: Data curation, Writing-Review & Editing; HS, JR, HS & RA: Investigation; BK: Conceptualization, Methodology, Supervision.


1. World Health Organization (WHO). Air pollution. [cited Jan 22, 2023]. Available from:
2. Dominski FH, Lorenzetti Branco JH, Buonanno G, Stabile L, Gameiro da Silva M, Andrade A. Effects of air pollution on health: A mapping review of systematic reviews and meta-analyses. Environ Res 2021;201:111487.
3. Pöschl U. Atmospheric aerosols: Composition, transformation, climate and health effects. Angew Chem Int Ed Engl 2005;44(46):7520–7540.
4. Anderson HR, Atkinson RW, Peacock JL, Marston L, Konstantinou K, et al. Meta-analysis of time-series studies and panel studies of particulate matter (PM) and ozone (O3): report of a WHO task group WHO Regional Office for Europe, Copenhagen; 2004. Available from:
5. Pope CA, Dockery DW. Health effects of fine particulate air pollution: Lines that connect. Journal of the Air & Waste Management Association 2006;56(6):709–742.
6. Cohen AJ, Brauer M, Burnett R, Anderson HR, Frostad J, Estep K, et al. Estimates and 25-year trends of the global burden of disease attributable to ambient air pollution: an analysis of data from the Global Burden of Diseases Study 2015. Lancet 2017;389(10082):1907–1918.
7. United States, Environmental Protection Agency (US EPA). Integrated Science Assessment (ISA) for particulate matter 2009. [cited Jan 22, 2023]. Available from:
8. Health Canada. Canadian smog science assessment, Volume 2: health effects 2013. [cited Jan 22, 2023]. Available from:
9. Chen J, Hoek G. Long-term exposure to PM and all-cause and cause-specific mortality: A systematic review and metaanalysis. Environ Int 2020;143:105974.
10. Pope 3rd CA, Burnett RT, Thun MJ, Calle EE, Krewski D, Ito K, et al. Lung cancer, cardiopulmonary mortality, and long-term exposure to fine particulate air pollution. JAMA 2002;287(9):1132–1141.
11. Pope 3rd CA, Burnett RT, Thurston GD, Thun MJ, Calle EE, Krewski D, et al. Cardiovascular mortality and long-term exposure to particulate air pollution: epidemiological evidence of general pathophysiological pathways of disease. Circulation 2004;109(1):71–77.
12. Beelen R, Hoek G, van den Brandt PA, Goldbohm RA, Fischer P, Schouten LJ, et al. Long-term effects of traffic-related air pollution on mortality in a Dutch cohort (NLCS-AIR study). Environ Health Perspect 2008;116(2):196–202.
13. Colais P, Serinelli M, Faustini A, Stafoggia M, Randi G, Tessari R, et al. Air pollution and urgent hospital admissions in nine Italian cities. Results of the EpiAir Project. Epidemiol Prev 2009;33(6 Suppl 1):77–94.
14. Cohen AJ, Anderson HR, Ostro B, Pandey KD, Krzyzanowski M, Künzli N, et al. The global burden of disease due to outdoor air pollution. J Toxicol Environ Health A 2005;68(13-14):1301–1307.
15. Künzli N, Tager IB. Air pollution: from lung to heart. Swiss Med Wkly 2005;135(47-48):697–702.
16. Sharma RK, Agrawal M. Biological effects of heavy metals: an overview. J Environ Biol 2005;26(2 Suppl):301–313.
17. Huang YCT, Ghio AJ. Vascular effects of ambient pollutant particles and metals. Current Vascular Pharmacology 2006;4(3):199–203.
18. Manisalidis I, Stavropoulou E, Stavropoulos A, Bezirtzoglou E. Environmental and health impacts of air pollution: A review. Front Public Health 2020;8:14.
19. Kelly FJ, Fussell JC. Size, source and chemical composition as determinants of toxicity attributable to ambient particulate matter. Atmospheric Environment 2012;60:504–526.
20. Riediker M, Cascio WE, Griggs TR, Herbst MC, Bromberg PA, Neas L, et al. Particulate matter exposure in cars is associated with cardiovascular effects in healthy young men. Am J Respir Crit Care Med 2004;169(8):934–940.
21. Lasley SM, Gilbert ME. Glutamatergic components underlying lead-induced impairments in hippocampal synaptic plasticity. Neurotoxicology 2000;21(6):1057–1068.
22. Lasley SM, Green MC, Gilbert ME. Rat hippocampal NMDA receptor binding as a function of chronic lead exposure level. Neurotoxicol Teratol 2001;23(2):185–189.
23. World Health Organization (WHO). WHO global air quality guidelines: particulate matter (PM2.5 and PM10), ozone, nitrogen dioxide, sulfur dioxide and carbon monoxide. [cited Jan 22, 2023]. Available from:
24. World Health Organization (WHO). Burden of disease from ambient air pollution for 2016, version 2. [cited Jan 22, 2023]. Available from:
25. Silva RA, West JJ, Lamarque JF, Shindell DT, Collins WJ, Faluvegi G, et al. Future global mortality from changes in air pollution attributable to climate change. Nat Clim Chang 2017;7(9):647–651.
26. Croitoru L, Sarraf M. Estimating the health cost of air pollution: The case of Morocco. Journal of Environmental Protection 2017;8(10):1087–1099.
27. World Health Organization (WHO). Ambient air pollution data. [cited Jan 22, 2023]. Available from:
28. Houssaini ASS, Messaouri H, Nasri I, Roth MP, Nejjari C, Benchekroun MN. Air pollution as a determinant of asthma among schoolchildren in Mohammedia, Morocco. International Journal of Environmental Health Research 2007;17(4):243–257.
29. Abrouki Y, Anouzla A, Azzallou R, Jioui I, Kastali M, Hatib A, et al. Effects of environmental air pollution on pediatric asthma in Morocco. Green and Applied Chemistry 2021;12:31–39.
30. Nejjari C, Marfak A, Rguig A, Maaroufi A, El Marouani I, El Haloui A, et al. Ambient air pollution and emergency department visits among children and adults in Casablanca, Morocco. AIMS Public Health 2021;8(2):285–302.
31. Ministry of Land Management, Water and Environment, Ministry of Health. Evaluation of the impact of air pollution on the health of populations in the greater Casablanca area: Casa AirPol study 2000 [French].
32. Ministry of Land Management, Water and Environment, Ministry of Health. Evaluation of the impact of air pollution on the health of asthmatic children in Mohammedia: Mohammadia AirPol study, summary of results (French) 2003. [cited Jan 22, 2023]. Available from:
33. WHO Regional Office for Europe, European Centre for Environment and Health. AirQ+: software tool for health risk assessment of air pollution. Bonn: WHO Regional Office for Europe 2019.
34. Ostro B, ; Occupational and Environmental Health Team, ; World Health Organization (WHO). Outdoor air pollution : assessing the environmental burden of disease at national and local levels/Bart Ostro World Health Organization; 2004. [cited Jan 22, 2023]. Available from:
35. Krzyzanowski M. Methods for assessing the extent of exposure and effects of air pollution. Occupational and Environmental Medicine 1997;54(3):145–151.
36. World Health Organization (WHO). WHO SDG Indicator 11.6.2 Concentrations of fine particulate matter (PM2.5). [cited Jan 22, 2023]. Available from:
37. Li Y, Chen Z, Li J. How many people died due to PM2.5 and where the mortality risks increased? A case study in Beijing. In: 7th IEEE International Geoscience and Remote Sensing Symposium (IGARSS). Fort Worth; 2017.
38. Baxter LK, Crooks JL, Sacks JD. Influence of exposure differences on city-to-city heterogeneity in PM2.5-mortality associations in US cities. Environ Health 2017;16(1):1–8.
39. Chowdhury S, Dey S, Smith KR. Ambient PM2.5 exposure and expected premature mortality to 2100 in India under climate change scenarios. Nature Communications 2018;9:318.
40. Ansari M, Ehrampoush MH. Meteorological correlates and AirQ+ health risk assessment of ambient fine particulate matter in Tehran, Iran. Environ Res 2019;170:141–150.
41. Yarahmadi M, Hadei M, Nazari SSH, Conti GO, Alipour MR, Ferrante M, et al. Mortality assessment attributed to longterm exposure to fine particles in ambient air of the megacity of Tehran, Iran. Environ Sci Pollut Res Int 2018;25(14):14254–14262.
42. Faridi S, Shamsipour M, Krzyzanowski M, Künzli N, Amini H, Azimi F, et al. Long-term trends and health impact of PM2.5 and O3 in Tehran, Iran, 2006–2015. Environ Int 2018;114:37–49.
43. Hadei M, Hopke PK, Nazari SSH, Yarahmadi M, Shahsavani A, Alipour MR. Estimation of mortality and hospital admissions attributed to criteria air pollutants in Tehran Metropolis, Iran (2013–2016). Aerosol and Air Quality Research 2017;17(10):2474–2481.
44. Burnett RT, Pope 3rd CA, Ezzati M, Olives C, Lim SS, Mehta S, et al. An integrated risk function for estimating the global burden of disease attributable to ambient fine particulate matter exposure. Environ Health Perspect 2014;122(4):397–403.
45. Ostro B, Spadaro JV, Gumy S, Mudu P, Awe Y, Forastiere F, et al. Assessing the recent estimates of the global burden of disease for ambient air pollution: Methodological changes and implications for low- and middle-income countries. Environ Res 2018;166:713–725.
46. High Commission for Planning (HCP). General Census of Population and Habitat. Demographic and Socio-Economic Characteristics of the Population. National Report 2014 [French].
47. Pope 3rd CA, Burnett RT, Turner MC, Cohen A, Krewski D, Jerrett M, et al. Lung cancer and cardiovascular disease mortality associated with ambient air pollution and cigarette smoke: shape of the exposure-response relationships. Environ Health Perspect 2011;119(11):1616–1621.
48. Tomczak A, Miller AB, Weichenthal SA, To T, Wall C, van Donkelaar A, et al. Long-term exposure to fine particulate matter air pollution and the risk of lung cancer among participants of the Canadian National Breast Screening Study: Fine particulate matter and incident lung cancer. Int J Cancer 2016;139(9):1958–1966.
49. Sonego M, Pellegrin MC, Becker G, Lazzerini M. Risk factors for mortality from acute lower respiratory infections (ALRI) in children under five years of age in low and middle-income countries: a systematic review and meta-analysis of observational studies. PLoS ONE 2015;10(1):2–17.
50. Darrow LA, Klein M, Flanders WD, Mulholland JA, Tolbert PE, Strickland MJ. Air pollution and acute respiratory infections among children 0–4 Years of age: An 18-year Time-Series Study. American Journal of Epidemiology 2014;180(10):968–977.
51. Bouchriti Y, Ait Haddou M, Kabbachi B, Achbani A, Cherrat Z, Rida J et al. Ambient air quality and health impact of exposure to outdoor air pollution in the Moroccan population: A systematic review. Pollution 2023;9(2):660–677.

Article information Continued

Figure 1.

Location map of the study area: (a) Morocco in Africa; (b) Study area in Agadir Ida Ou Tanan Prefecture; (c) Study area; Locations of the mobile air pollution measuring station in 2016 (green circles)

Table 1.

Population characteristics.

Health endpoint Age (years) At-risk population Hospital cases or deaths Annual incidence per 100,000 persons
Total mortality ≥ 30 202,999 1,966 deaths 968.48
LC ≥ 25 240,825 - 25.53 1
ALRI < 5 54,731 336 cases 613.91 2
COPD ≥ 25 240,825 28 cases 3 11.63

Standardized incidence for all Morocco regions (Casablanca region cancer registry results in 2004).


Annual cumulative incidences calculated among students examined in preschool and elementary school as part of the Communicable Ophthalmia Campaign in 2015.


The SMRHC's pneumology service reported cases in 2016.

LC: Lung cancer, ALRI: Acute lower respiratory tract infections, COPD: Chronic obstructive pulmonary disease.

Table 2.

PM2.5 health impact assessment results in AirQ+ in 2016.

Exposure Outcome Health endpoint Age (years) RR (CI 95%) 1 AP (%) 2 NAP 3
Long-term Mortality Total mortality 4 ≥ 30 1.16 (1.10-1.22) 14.19 (9.5-18.37) 279 (187-361)
LC 5 ≥ 25 1.23 (1.12-1.37) 18.42 (11.02-26.26) 11 (7-16)
Morbidity ALRI 6 < 5 1.17 (1.10-1.24) 14.36 (8.78-19.33) 33 (20-44)
COPD 7 ≥ 25 1.19 (1.11-1.31) 15.68 (9.66-23.38) 4 (3-7)
Short-term Mortality Total mortality ≥ 30 1.01 (1.00-1.02) 1.27 (0.47-2.06) 25 (9-40)

Relative risk with its 95% confidence interval for a threshold value of 10 μg/m3.


Estimated attributable proportion (%).


Estimated number of attributable cases.


RR estimation by the log-linear model using formula: RR (x) = eβ×(x-x0), β = 0,006015392281974714 (0,006015392281974714 – 0,007973496801885352).


RR estimation method: GBD 2015/2016 (integrated function 2016 with WHO interim target 1). Threshold value: 2.4μg/m3.


RR estimation method: GBD 2015/2016 (integrated function 2016 vs WHO AGQ value). Threshold value: 10 μg/m3


RR estimation method: GBD 2015/2016 (integrated function 2016 with WHO 3 interim target). Threshold value: 15 μg/m3.


RR estimation method: GBD 2015/2016 (integrated function 2016). Threshold value: 2.4 μg/m3 [6].

LC: Lung cancer, ALRI: Acute lower respiratory tract infections, COPD: Chronic obstructive pulmonary disease, GBD: Global Burden of Disease.