Next Article in Journal
The Impact of Three Communication Channels on the Dissemination of a Serious Game Designed to Enhance COVID-19 Prevention
Previous Article in Journal
Effects of Self-Rated Health Status on Residents’ Social-Benefit Perceptions of Urban Green Space
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Air Pollution in Kosovo: Short Term Effects on Hospital Visits of Children Due to Respiratory Health Diagnoses

by
Zana Shabani Isenaj
1,
Merita Berisha
1,2,
Dragan Gjorgjev
3,
Mirjana Dimovska
3,
Hanns Moshammer
4,5,* and
Antigona Ukëhaxhaj
2,6
1
Medical Faculty, University of Hasan Pristina, George Bush 31, 10000 Pristina, Kosovo
2
National Institute of Public Health, St. Mother Teresa pn, Rrethi i Hospitalit, 10000 Pristina, Kosovo
3
Faculty of Medicine, Ss. Cyril and Methodius University in Skopje, 1000 Skopje, North Macedonia
4
Department of Environmental Health, ZPH, Medical University of Vienna, 1090 Vienna, Austria
5
Department of Hygiene, Medical University of Karakalpakstan, Nukus 230100, Uzbekistan
6
Faculty of Medicine, University Fehmi Agani, Ismail Qemali n.n., 50000 Gjakove, Kosovo
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(16), 10141; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph191610141
Submission received: 14 July 2022 / Revised: 11 August 2022 / Accepted: 15 August 2022 / Published: 16 August 2022
(This article belongs to the Section Environmental Health)

Abstract

:
The Republic of Kosovo is a small country in the Balkans. The capital city of Pristina hosts most of its population and is situated in a mountain basin with poor air exchange, especially during winter. Domestic heating, road transport, industry and coal-fired power plants contribute to high levels of air pollution. We performed a time-series analysis on effects of particulate air pollution (PM2.5) on respiratory health of children and adolescents, using hospital admission and ambulatory visit numbers from the pediatric university clinic. From 2018 until 2020, daily mean concentrations of PM2.5 ranged between 2.41 and 161.03 µg/m³. On average, there were 6.7 ambulatory visits per day with lower numbers on weekends and during the first COVID-19 wave in 2020. An increase in PM2.5 led to an immediate increase in visit numbers that lasted over several days. Averaged over a full week, this amounted to about a 1% increase per 10 µg/m³. There were, on average, 1.7 hospital admissions per day. Two and three days after a rise in air pollution, there was also a rise in admission numbers, followed by a decline during the consecutive days. This might indicate that the wards were overstressed because of high admission numbers and restricted additional admissions.

1. Introduction

Air pollution is the largest contributor to the environmental burden of disease, causing and exacerbating several diseases including cardiovascular diseases such as stroke and ischemic cardiopathy, acute respiratory infections, chronic obstructive pulmonary disease and cancer [1]. According to the World Health Organization (WHO), in 2012, one out of every nine deaths was the result of air pollution-related conditions. Of these deaths, around 3 million are attributable solely to ambient (outdoor) air pollution [2]. Respiratory diseases are among the leading causes of both morbidity and mortality worldwide. Lung infections, chronic obstructive pulmonary diseases, asthma, and bronchial cancer are main contributors of respiratory-related mortality and morbidity [3,4]. By 2030, the WHO estimates that respiratory diseases will be the most important reason for impaired quality of life and death within the EU that will require substantial economic resources for both prevention and treatment [5].
The strongest and most consistent association between air pollution exposure and respiratory morbidities has been seen for ambient particulate matter. Particulate matter (PM) is a mixture of solid particles and liquid droplets that vary in size and origin. Commonly used indicators describing PM that are relevant to health refer to the mass concentration of particles with a diameter of less than 10 µm (PM10) and of particles with a diameter of less than 2.5 µm (PM2.5). Findings from many epidemiological and toxicological studies have highlighted adverse effects of air pollution on both premature mortality and morbidity from respiratory and cardiovascular disease, following both short-term and chronic exposure [6,7]. Apart from other vulnerable groups, children are also susceptible to air pollution. For example, we have shown short-term effects of air pollution on lung function in asthmatic children [8]. We also demonstrated long-term beneficial effects on lung function with declining air pollution [9]. The relationship between air pollutants and respiratory hospital admissions has been reported both in developed countries and in developing countries [10,11]. Studies have also demonstrated a significant impact of air pollution on infants and children, which is manifested primarily as a range of respiratory problems [12].
Several epidemiological studies have found associations between daily concentrations of ambient air pollution, such as particulate matter (PM), ozone (O3), and nitrogen dioxide (NO2), and daily mortality and morbidity outcomes due to respiratory diseases [13,14,15,16,17]. Children are particularly susceptible to ambient air pollution due to their physiological and behavioral characteristics. Previous studies found that children were the most susceptible to hospital admission for respiratory disease associated with PM levels and meteorological factors [18,19,20,21]. Children are considered to be more sensitive to air pollution than adults, and children with asthma are particularly vulnerable to the adverse health effects of air pollution [22]. Ségala and colleagues [23] found PM10 and nitrogen dioxide (NO2) to be positively associated with consultations and hospitalizations for bronchiolitis during the winter. Finally, Rahman and colleagues [24] found higher odds of acute bronchiolitis associated with exposure to PM10, among children who reside in urban areas. A recent study revealed that acute bronchiolitis-related hospitalization among children was associated with temperature and exposure to NO2 and PM10 [25] suggesting a need to adopt sustainable clean air policies to protect young children’s health.
A significant proportion of urban areas in Kosovo suffer from poor air quality, with ambient concentrations of particulate matter with a diameter of 2.5 μm or less (PM2.5) significantly exceeding the national and European Union (EU) standards (20 µg/m³ annual mean [26]) and global air quality guidelines for PM2.5 (5 µg/m³ annual mean [27]) established by the WHO [28]. Air pollution is mainly associated with the burning of solid fuels in homes, obsolete large thermal power plants, industry and exhaust gas from vehicles and has become a serious environmental problem, and its effect on human health is a concern [29].
Kosovo has one of the youngest populations in Europe accompanied by a low life expectancy of 71.6 years, the lowest in the region [30]. The population in Kosovo is growing with now reaching about 1.8 million inhabitants, more than 650,000 of which are aged 19 or younger [31]. The health care system is composed of the public health care network (primary, secondary and tertiary levels) and facilities in private ownership, funded largely from government and municipal tax and direct payments [32]. The main concerns in relation to the health system are the quality of health services and low financing. A 2014 report stated that citizens in Kosovo had to cover expenses for health services with out-of-pocket expenditures ranging from 30 to 50% of total costs [31]. At the same time, the launch of ambitious healthcare reforms, in combination with the fragile structure and governance of the system, have impacted the overall stability and continuity of the service provision [32].
Although the air pollution burden in Kosovo is not well-documented, in part due to not fully functional health information system, a recent report from the World Bank estimates that about 760 people die prematurely every year in Kosovo because of exposure to ambient air pollution, where 11 percent of these cases are from Pristina, representing the most polluted and most populated city (more than 200,000 persons in 2021 [31]) in Kosovo. Of the total number of pollution-related deaths, 90 percent are due to ischemic heart disease and stroke combined [28]. The capital city, Pristina, faces severe smog episodes in the winter, rivaling those of big cities such as Beijing, Mumbai and New Delhi–causing significant health impact [28]. This study aims to examine associations between air pollution and children’s health, thus contributing to local data-based evidence of the environmental burden of disease in Kosovo.

2. Materials and Methods

2.1. Air Pollution Data

Kosovo is a small land-locked country located in the mountainous regions of the Balkan Range. Most of its population lives in the shallow basin in and around the capital city of Pristina and most people living in this area most likely experience a similar time-course of air pollution exposure. The whole population in this basin is served by the hospitals in Pristina with the pediatric university clinic being the primary point of care for children and adolescents.
Data from the monitoring network are collected by the Kosovo Hydro-Meteorological Institute (KHMI) and reported to the European Environment Agency (EEA) [33] and published on the KHMI website [34]. Daily averages were available from both sources starting from 13 May 2018. Over time, the number of monitoring stations and completeness of data increased for all of Kosovo. Daily data for PM2.5, PM10, and NO2, were collected from the EEA website from all available stations until 1 March 2022. Daily average temperature data from Pristina (station at KHMI) were available in annual meteorological reports. The latest report available at the time of data collection early in 2022 was the 2020 report. The US embassy in Pristina measures PM2.5 since several years and daily and hourly values are reported on the AirNow website of the US government [35]. Daily data were extracted from 23 March 2016, until 1 March 2022.

2.2. Health Data

The number of daily ambulatory visits from the district of Pristina were obtained from the Health Information Statistics Platform from Ministry of Health in Kosovo, and also retrieved manually from the protocols of specialist ambulatories of the Pediatric Clinic at the University Clinical Center of Kosovo. Data of the daily hospital admissions were obtained manually from the protocols of the Pediatric Clinic from the University Clinical Center of Kosovo, due to the lack of electronic data. The data included basic information such as sex, age, date of visit or hospital admission and disease diagnosis.
We included data from all children aged 0–18 years admitted to the hospital with respiratory disease (International Classification of Diseases, 10th Revision: J00-J99) between January 2018 and December 2020 and children who visited the ambulatories. Patients who came from locations other than municipality of Pristina were excluded.

2.3. Statistical Analysis

In a first preparatory step, the pairwise Pearson correlation between monitoring stations was assessed for all pollutants and both for the whole observation period and also stratified by month. The pollutant with the strongest correlation between stations was selected as the exposure variable. Daily data were primarily taken from the station that was most closely correlated with the other stations. Missing data at that station were replaced by estimated values based on linear regression between the index station and the other station that was most closely correlated with it.
Next, for each of the two outcomes (number of ambulatory visits and number of hospital admissions), a separate General Additive Model (GAM) was constructed, using the library mgcv in R [36]. To that end, a script was adapted that was originally written by Daniel Rabczenko for the analysis of daily deaths in Vienna, Graz and Linz [37,38]. That script used over-dispersed Poisson GAMs [39], applying a 7-day polynomially distributed lag (pdl) model [40] and stringent convergence criteria [41,42], considering day of the week as “as factor” variable. Seasonal and longer-term variations were modelled using spline functions [43]. In order to find the optimal number of knots (or degrees of freedom) for the non-linear fit, we chose the model with the smallest partial autocorrelation (sum over 30 days). Next, we added temperature (natural spline) to the model. We selected the most relevant lag (lags 0–3) based on the Akaike information criterion (AIC) [44]. Next, we performed the same selection process on the change in temperature between consecutive days. As a sensitivity analysis, splines for the temporal variation were also examined with different number of knots. In addition, effects of single lags were also assessed to check if a third-degree polynomial was a valid approximation of the lag structure of the PM2.5 impact.

3. Results

3.1. Air Pollution Data

PM2.5 displayed the strongest temporal correlation between monitoring stations. The data from the KHMI urban background station were highly correlated with all other stations situated in urban or sub-urban backgrounds. Pearson’s correlation coefficient was generally larger than 0.8, with the exception of the urban background station of Prizren, for which the correlation coefficient with the KHMI station was only 0.765. Details of the monitoring stations in Kosovo and the temporal correlation coefficients between the stations is provided in Appendix A. PM10 and NO2 displayed somewhat weaker correlations, although still quite strong and usually stronger in winter. PM2.5 was chosen as the best measure of exposure, because its time-series reflected best the temporal exposure variation of the total population. Moreover, for this pollutant the US embassy provided an additional series that was strongly correlated with the KHMI urban background station (r = 0.964) that allowed for filling the gaps on days with missing KHMI data using the linear regression formula:
KHMI = 0.8165507 × embassy + 1.152676
The average concentration of PM2.5 was 24.14 µg/m³ and the first, second (median) and third quartile were 11.14, 16.01 and 28.35 µg/m³ (IQR = 12.34). The maximal and the minimal daily concentrations were 2.41 and 161.03 µg/m³.

3.2. Time Course of Ambulatory Visits and Hospital Admissions

Based on the above inclusion and exclusion criteria, a total of 1838 hospital admission records and 7372 ambulatory visits were collected during this study period. On average, there were 6.7 ambulatory visits per day (quartiles: 0, 4, 12, maximum 96). On average, 1.7 children were admitted per day because of respiratory diseases (quartiles: 0, 1, 3, maximum: 18).
Both time-series of health data displayed some seasonal fluctuation with somewhat higher numbers in winter. However, most striking was a steep decline both in ambulatory visits and hospital admissions during the first COVID-19 wave in 2020 (Figure 1). Moreover, during weekends much fewer events occurred.

3.3. Short-Term Effect of Particulate Air Pollution (PM2.5)

An increase in air pollution led to an immediate (same day) increase in ambulatory visits. This increase was significant on the first and second day after the pollution episode. Afterwards, there was a decline in visits. However, added over a whole week, the total effect of an increase in air pollution was significantly positive (Figure 2). An increase in PM2.5 daily mean of 10 µg/m³, added over a full week, resulted in an increase in hospital admissions by 0.99% (95% confidence interval: 0.4; 1.6).
Stationary hospital admissions only increased on the day following the pollution episode This effect was significant on the second and third day after the pollution episode and afterwards declined again and finally turned negative (Figure 2). Summed over the whole week after the pollution episode, no significant effect on hospital admissions could be observed (−0.2%; CI: −2.6; 2.3%).
For both outcomes, changing the number of knots for the temporal/seasonal variation did not substantially affect the main results (Figure A1, Figure A2 and Figure A3). Single lag analyses confirmed that a third-degree polynomial provided a reasonable approximation of the true lag structure (Figure A4).

4. Discussion

A steep decline in the number of ambulatory visits and hospital admissions was observed during the first wave of COVID-19. This phenomenon was attributed to the decrease in the demand for services from the citizens due to the restrictive measures introduced to mitigate the COVID-19 outbreak and the repurposing of health services and health workforce capacities to COVID-19 treatment centers. According to a report from WHO, this situation was also impacted by the lack of Essential Health Service Packages and lack of additional funds to maintain EHS. A rapid assessment of the situation analysis of the impact of COVID-19 in EHS revealed that all non-acute health services in Kosovo were suspended from March till May 2020, with a gradual re-opening afterwards [45]. In addition, individual protective measures such as social distancing and wearing of face masks could also have prevented other respiratory diseases.
This study is only based on three years of observation. This is not very long for a time series analysis. The study may therefore lack sufficient power to detect small effects. However, a longer time series was not feasible due to limitations in data availability. Temperature data (a relevant confounder) were only available until 2020. Air pollution data only became available in 2018. However, our findings are plausible.
We did see an increase in ambulatory visits with increasing levels of PM2.5. This effect was already observed on the same day but reached significance only on the first and second day after the pollution episode. Averaged over 1 week, the effect size per 10 µg/µ³ was about 1%. This order of magnitude has been observed in other similar studies [46,47,48,49].
We expected to find an increase in case numbers soon after the pollution day and a decline in numbers a few days later. Such a shape is well represented by a third-degree polynomial and therefore, we decided in advance to apply such a function in the distributed lag model. Before applying the polynomial, we calculated the effect estimates for single lags to make sure that our assumption was correct. Such a single lag estimate is confounded by the effects from neighboring lags because of strong short-term autocorrelation of pollution concentrations. Therefore, single-day lags cannot inform about the true effect strength but provide valid information about the overall temporal distribution of the effect.
This bi-phasic shape of the temporal exposure-response distribution has been discussed for several decades under the term “harvesting hypothesis” or “harvesting effect” [50,51,52,53,54,55]. Maybe, the term “harvesting” originally really referred to the work of a farmer: During harvest time the farmer will try to harvest the crops during sunny weather, because when the harvest is performed in rainy weather, there is an increased danger of mold. Therefore, in a time series, sunny days will be associated with a higher harvest yield. However, after a sunny day, all the crops on the farmer’s field will already have been harvested. So, on a later day, independent from rain or sunshine, a lower yield is to be expected. However, when this analogy was applied to time-series studies on daily mortality, the term “harvesting” quickly got associated with the mythical concept of the grim reaper. Indeed, some authors argued that an increase in death numbers at or soon after a high pollution day was not a sign of a relevant health effect because death only affected those that otherwise would have died just a few days later. However, further analyses (also applying distributed lag models [52] or compartment models [54]) confirmed, that only a small part of the excess mortality is only due to very short term replacement (“harvesting”).
Of course, every person can only die once. Therefore, when more people out of the “very frail pool” [53] die on one day, there are less people in this pool left to die on a consecutive day. With hospital admissions and ambulatory visits, this is different in the sense that patients can also count as “cases” more often than once. Still, an increase in visit frequency per person would also be a relevant effect of air pollution. Due to anonymous data, we cannot check how many of the visits in our study are first or repeated visits. However, this does not invalidate the results of our study. In comparison, although less well studied, a biphasic temporal shape of the effect estimate is also expected for visits.
Thus, relative risks obtained in time-series analyses refer to the vulnerable group mostly and cannot be applied directly to the general population. Therefore, an estimate of absolute effects is provided here: Per week, there were on average 47 ambulatory visits. An increase per 10 µg/m³ on one day, would therefore lead to about one-half additional visit (0.47%) spread over the coming week. Considering that 10 µg/m³ are nearly the interquartile range of 12.34, achieving a reduction in PM2.5 by 10 µg/m³ would certainly not be an easy task. However, over the whole year, this would reduce the number of visits by about ½ per day, or by about 180 visits. Since the annual mean concentration of PM2.5 was even more than twice as high (24.14 µg/m³), PM2.5 in total was responsible for approximately 1 visit per day. Moreover, this is likely an underestimation, given the fact that air pollution concentrations peak in winter when hospital admissions are highest for respiratory diseases. It is also an underestimation, because part of the study coincided with the 2020 corona lock-down measures. Of course, this is only considering one single outcome (respiratory diseases) and only in the child and adolescent population. Additionally, this impact only considers one although the most important hospital serving this population in Kosovo.
Stationary hospital admissions rose after an increase in air pollution with a latency of one day. This seems plausible, as admissions usually need a consultation with a general practitioner/a pediatrician in primary health care first. This step would allow for a delay of one day. In Austria, we even observed a longer delay for hospital visits because of respiratory diagnoses (2 days for city dwellers and even longer for rural people) [56].
Most likely because of this delay, the biphasic curve produced by the polynomial distributed lag model had not returned to the baseline by lag 6. Therefore, we also examined a longer lag period of 10 days (lag 0–lag 9, Figure A3). In that model, the biphasic curve is clearly complete, but the point estimates for the later days are really imprecise. Data from only 3 years seem insufficient for the analysis of longer lag periods. Thus, we believe that a lag period of 7 days was a good choice.
After this delay, for stationary admissions a significant effect of air pollution was also seen. However, this increase was followed by a decrease on days 5 and 6. Thus, over the whole week, air pollution did not lead to an increase in admissions. If this was to be interpreted as a “perfect harvesting effect” whereby “all crops that could be harvested, had been harvested on the first sunny day”, it would mean that “all children that could be admitted to the hospital have already been admitted during the first days of the pollution episode”. This interpretation is unlikely. It more likely demonstrates that “all hospital beds that could be occupied were already occupied soon after the beginning of each pollution episode, indicating that the hospital ward is already reaching the limit of its resources. When the ward gets overfilled, either due to more admissions or longer duration of stay because of more severe cases, the doctors will be more restrictive in accepting new patients. Therefore, a drop in admissions after the peak triggered by air pollution very likely indicates that the hospital resources are strained to their limits by an air pollution episode. Moreover, in this respect it is noteworthy, that air pollution episodes occur mostly in wintertime, when respiratory diseases are already most frequent. Interpreted this way, a net zero effect over a week in spite of a clear and plausible response on single days, should rather be a cause for concern.

5. Conclusions

The study found plausible and consistent short-term effects of particulate air pollution on respiratory morbidity of children. The effect was most clearly seen in the number of ambulatory visits. It is generally established that short-term effects are not as important as the long-term health consequences of chronic air pollution exposure. However, the study of chronic effects is much more difficult, as it requires longer series of monitoring data. Time-series analyses are an elegant way of demonstrating health effects of air pollution. However, impact assessments solely based on the results of time-series analyses would severely underestimate the magnitude of the impact. Nevertheless, even this time-series study clearly demonstrates the dire need for air pollution abatement measures in Kosovo.

Author Contributions

Conceptualization, Z.S.I.; methodology, H.M.; formal analysis, H.M.; investigation, Z.S.I.; resources, A.U.; data curation A.U.; writing—original draft preparation, H.M. and Z.S.I.; writing—review and editing, D.G. and M.D.; supervision, M.B.; funding acquisition, Z.S.I. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Austrian Development Cooperation and the Government of Kosovo under the HERAS plus (Higher Education, Research and Applied Science plus), grant number K-02-2021.

Institutional Review Board Statement

This is a purely observational study using only anonymous data. Therefore, a review board statement was not necessary.

Informed Consent Statement

Patient consent was waived due to the use of anonymous preexisting data only.

Data Availability Statement

Air pollution and temperature data have been obtained from the internet from freely available sources. The links are provided in the text. Numbers of ambulatory visits and of stationary hospital admissions are available upon request. Send an e-mail to the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Whole-year temporal correlation of PM2.5 between monitoring stations.
Table A1. Whole-year temporal correlation of PM2.5 between monitoring stations.
ABCDEFGHIJKem
A1
B0.92481
C0.90560.85761
D0.87140.86620.86651
E0.85490.84330.84430.86211
F0.7650.72490.77960.75430.87451
G0.82880.82120.77690.85740.83860.68821
H0.82540.80680.74970.76850.83420.72980.75411
I0.91720.84410.90230.85790.86220.73330.81220.74221
J0.92080.86730.92420.87770.87870.78710.83420.82770.94671
K0.90560.83960.91920.84560.8450.74920.79330.80770.95840.96181
em0.90020.9160.87140.84520.85870.77680.79180.83580.91580.92190.91551
Table A2. Summertime (May–October) PM2.5 correlations between monitoring stations.
Table A2. Summertime (May–October) PM2.5 correlations between monitoring stations.
ABCDEFGHIJKem
A1
B0.95111
C0.92670.86431
D0.87710.83840.86261
E0.83870.85490.74220.68881
F0.78690.6830.80240.75940.79371
G0.75210.76550.75210.75440.58570.45861
H0.7930.84830.75570.70420.78690.56890.68051
I0.91670.81060.93420.92040.90280.89510.86410.68441
J0.92180.84310.91320.89790.90990.88360.83360.79790.95611
K0.92270.8240.94360.90120.90830.9040.82820.73470.98180.9641
em0.90260.85610.81270.77420.82680.82850.59520.65640.92910.91650.93141
Table A3. Winter time (November-April) PM2.5 correlations between monitoring stations.
Table A3. Winter time (November-April) PM2.5 correlations between monitoring stations.
ABCDEFGHIJKem
A1
B0.89881
C0.87670.81351
D0.8280.81310.82581
E0.79970.76790.79420.80561
F0.68020.61260.70390.64140.81181
G0.79150.75560.71560.81350.77120.57551
H0.77070.72010.67310.67590.74640.62190.65641
I0.92790.85590.90570.84910.81850.68780.7810.76411
J0.90840.83040.90260.830.82480.70830.76480.79650.94961
K0.88920.80640.89810.79670.77910.66940.71910.79480.95340.95151
em0.88580.89590.83970.7970.80290.69270.75780.79770.92420.91140.90321
Table A4. Description of monitoring stations.
Table A4. Description of monitoring stations.
AbbreviationLocationTypeOperative Since
AAQMS Prishtine–IHMKsuburban backgroundMarch 2018
BAQMS Prishtine, Rilindjeurban backgroundApril 2018
CDrenas-Municipality urban backgroundApril 2018
DMitrovice urban backgroundMay 2018
EAQMS Pejeurban backgroundMay 2018
FAQMS Prizren urban backgroundApril 2018
GAQMS Hani i Elezit suburban backgroundNovember 2019
HAQMS Gjilan urban backgroundApril 2018
IAQMS Palaj, Obiliq suburban backgroundOctober 2019
JAQMS Obiliq, Obiliq urban backgroundSept 2019
KAQMS Dardhishte, Obiliq suburban backgroundOctober 2019
emUS embassy at least 2017
Figure A1. Sensitivity analysis: Effect (relative risk and 95% confidence interval) on ambulatory visits per 10 µg/m³ increase in PM2.5, total effect and daily effects over 7 consecutive days (lag 0–6): Optimal model (25 knots for the time smoother) compared to 35 and 15 knots for the 3-year period.
Figure A1. Sensitivity analysis: Effect (relative risk and 95% confidence interval) on ambulatory visits per 10 µg/m³ increase in PM2.5, total effect and daily effects over 7 consecutive days (lag 0–6): Optimal model (25 knots for the time smoother) compared to 35 and 15 knots for the 3-year period.
Ijerph 19 10141 g0a1
Figure A2. Sensitivity analysis: Effect (relative risk and 95% confidence interval) on stationary hospital admissions per 10 µg/m³ increase in PM2.5, total effect and daily effects over 7 consecutive days (lag 0–6): Optimal model (11 knots for the time smoother) compared to 16 and 6 knots for the 3-year period.
Figure A2. Sensitivity analysis: Effect (relative risk and 95% confidence interval) on stationary hospital admissions per 10 µg/m³ increase in PM2.5, total effect and daily effects over 7 consecutive days (lag 0–6): Optimal model (11 knots for the time smoother) compared to 16 and 6 knots for the 3-year period.
Ijerph 19 10141 g0a2
Figure A3. Sensitivity analysis: Effect (relative risk and 95% confidence interval) on stationary hospital admissions per 10 µg/m³ increase in PM2.5, total effect and daily effects over 10 consecutive days (lag 0–9): Optimal model (16 knots for the time smoother).
Figure A3. Sensitivity analysis: Effect (relative risk and 95% confidence interval) on stationary hospital admissions per 10 µg/m³ increase in PM2.5, total effect and daily effects over 10 consecutive days (lag 0–9): Optimal model (16 knots for the time smoother).
Ijerph 19 10141 g0a3
Figure A4. Effect (relative risk and 95% confidence interval) per 10 µg/m³ increase in PM2.5, daily effects over 7 single consecutive days (lag 0–6): (a) ambulatory visits, (b) stationary hospital admissions.
Figure A4. Effect (relative risk and 95% confidence interval) per 10 µg/m³ increase in PM2.5, daily effects over 7 single consecutive days (lag 0–6): (a) ambulatory visits, (b) stationary hospital admissions.
Ijerph 19 10141 g0a4

References

  1. Dimovska, M.; Gjorgjev, D. Assessing Health Impact of Air Pollution in Macedonian Cities. Biomed. J. Sci. Technol. Res. 2018, 10, 7522–7531. [Google Scholar] [CrossRef]
  2. Smith, K.R.; Bruce, N.; Balakrishnan, K.; Adair-Rohani, H.; Balmes, J.; Chafe, Z.; Dherani, M.; Hosgood, H.D.; Mehta, S.; Pope, D.; et al. HAP CRA Risk Expert Group. Millions dead: How do we know and what does it mean? Methods used in the comparative risk assessment of household air pollution. Annu. Rev. Public Health 2014, 35, 185–206. [Google Scholar] [CrossRef]
  3. Navarro-Torné, A.; Vidal, M.; Trzaska, D.K.; Passante, L.; Crisafulli, A.; Laang, H.; van de Loo, J.W.; Berkouk, K.; Draghia-Akli, R. Chronic respiratory diseases and lung cancer research: A perspective from the European Union. Eur. Respir. J. 2015, 46, 1270–1280. [Google Scholar] [CrossRef]
  4. Taj, T.; Malmqvist, E.; Stroh, E.; Oudin Åström, D.; Jakobsson, K.; Oudin, A. Short-Term Associations between Air Pollution Concentrations and Respiratory Health-Comparing Primary Health Care Visits, Hospital Admissions, and Emergency Department Visits in a Multi-Municipality Study. Int. J. Environ. Res. Public Health 2017, 14, 587. [Google Scholar] [CrossRef]
  5. WHO. Global Status Report on Noncommunicable Diseases 2010. 2015. Available online: http://www.who.int/nmh/publications/ncd_report2010/en/ (accessed on 1 February 2021).
  6. D’Antoni, D.; Smith, L.; Auyeung, V.; Weinman, J. Psychosocial and demographic predictors of adherence and non-adherence to health advice accompanying air quality warning systems: A systematic review. Environ. Health 2017, 16, 100. [Google Scholar] [CrossRef]
  7. Ukëhaxhaj, A.; Gjorgjev, D.; Ramadani, M.; Krasniqi, S.; Gjergji, T.; Zogaj, D. Air pollution in Pristina, influence on cardiovascular hospital morbidity. Med. Arch. 2013, 67, 438–441. [Google Scholar] [CrossRef]
  8. Moshammer, H.; Neuberger, M. The active surface of suspended particles as a predictor of lung function and pulmonary symptoms in Austrian school children. Atmos. Environ. 2003, 37, 1737–1744. [Google Scholar] [CrossRef]
  9. Neuberger, M.; Moshammer, H.; Kundi, M. Declining ambient air pollution and lung function improvement in Austrian children. Atmos. Environ. 2002, 36, 1733–1736. [Google Scholar] [CrossRef]
  10. Linden, B. Air pollution: Outdoor air quality and health. Br. J. Card. Nurs. 2019, 14, 1–4. [Google Scholar] [CrossRef]
  11. COMEAP. Review of the UK Air Quality Index. A Report by the Committee on the Medical Effects of Air Pollutants. 2011. Available online: https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/304633/COMEAP_review_of_the_uk_air_quality_index.pdf (accessed on 9 February 2021).
  12. Jakubiak-Lasocka, J.; Lasocki, J.; Badyda, A.J. The influence of particulate matter on respiratory morbidity and mortality in children and infants. Adv. Exp. Med. Biol. 2015, 849, 39–48. [Google Scholar] [CrossRef]
  13. Atkinson, R.W.; Butland, B.K.; Dimitroulopoulou, C.; Heal, M.R.; Stedman, J.R.; Carslaw, N.; Jarvis, D.; Heaviside, C.; Vardoulakis, S.; Walton, H.; et al. Long-term exposure to ambient ozone and mortality: A quantitative systematic review and meta-analysis of evidence from cohort studies. BMJ Open 2016, 6, e009493. [Google Scholar] [CrossRef] [PubMed]
  14. Bloemsma, L.D.; Hoek, G.; Smit, L.A.M. Panel studies of air pollution in patients with COPD: Systematic review and meta-analysis. Environ. Res. 2016, 151, 458–468. [Google Scholar] [CrossRef]
  15. Xing, Y.F.; Xu, Y.H.; Shi, M.H.; Lian, Y.X. The impact of PM2.5 on the human respiratory system. J. Thorac. Dis. 2016, 8, E69–E74. [Google Scholar] [CrossRef]
  16. Moshammer, H.; Poteser, M.; Kundi, M.; Lemmerer, K.; Weitensfelder, L.; Wallner, P.; Hutter, H.-P. Nitrogen-dioxide remains a valid air quality indicator. Int. J. Environ. Res. Public Health 2020, 17, 3733. [Google Scholar] [CrossRef]
  17. Moshammer, H.; Hutter, H.-P.; Kundi, M. Which metric of ambient ozone to predict daily mortality? Atmos. Environ. 2013, 65, 171–176. [Google Scholar] [CrossRef]
  18. Zhou, H.; Wang, T.; Zhou, F.; Liu, Y.; Zhao, W.; Wang, X.; Chen, H.; Cui, Y. Ambient Air Pollution and Daily Hospital Admissions for Respiratory Disease in Children in Guiyang, China. Front. Pediatr. 2019, 7, 400. [Google Scholar] [CrossRef]
  19. Jo, E.J.; Lee, W.S.; Jo, H.Y.; Kim, C.H.; Eom, J.S.; Mok, J.H.; Kim, M.H.; Lee, K.; Kim, K.U.; Lee, M.K.; et al. Effects of particulate matter on respiratory disease and the impact of meteorological factors in Busan, Korea. Respir. Med. 2017, 124, 79–87. [Google Scholar] [CrossRef]
  20. Zuurbier, M.; Lundqvist, C.; Salines, G.; Stansfeld, S.; Hanke, W.; Babisch, W.; Bistrup, M.L.; Van Den Hazel, P.; Moshammer, H. The Environmental Health of Children: Priorities in Europe. Int. J. Occup. Med. Environ. Health 2007, 20, 291–308. [Google Scholar] [CrossRef]
  21. Ukëhaxhaj, A.; Ramadani, N.; Gjorgjev, D.; Krasniqi, S.; Maloku, T.; Zogaj, D.; Gashi, S.; Loxha, A. Association of particulate air pollution with hospital admissions in pediatric clinics in Pristine, Kosovo. Albanian Med. J. 2013, 21. [Google Scholar]
  22. Yamazaki, S.; Shima, M.; Yoda, Y.; Oka, K.; Kurosaka, F.; Shimizu, S.; Takahashi, H.; Nakatani, Y.; Nishikawa, J.; Fujiwara, K.; et al. Exposure to air pollution and meteorological factors associated with children’s primary care visits at night due to asthma attack: Case-crossover design for 3-year pooled patients. BMJ Open 2015, 5, e005736. [Google Scholar] [CrossRef]
  23. Ségala, C.; Poizeau, D.; Mesbah, M.; Willems, S.; Maidenberg, M. Winter air pollution and infant bronchiolitis in Paris. Environ. Res. 2008, 106, 96–100. [Google Scholar] [CrossRef]
  24. Abdul Rahman, S.; Ismail, S.N.S.; Sahani, M.; Firuz, R.; Latif, M. A case crossover analysis of primary air pollutants association on acute respiratory infection (ARI) among children in urban region of Klang valley, Malaysia. Ann. Trop. Med. Public Health 2017, 10, 44. [Google Scholar] [CrossRef]
  25. Leung, S.Y.; Lau, S.Y.F.; Kwok, K.L.; Mohammad, K.N.; Chan, P.K.S.; Chong, K.C. Short-term association among meteorological variation, outdoor air pollution and acute bronchiolitis in children in a subtropical setting. Thorax 2021, 76, 360–369. [Google Scholar] [CrossRef]
  26. EUR-Lex. Directive 2008/50/EC of the European Parliament and of the Council of 21 May 2008 on Ambient Air Quality and Cleaner Air for Europe. 2008. Available online: https://eur-lex.europa.eu/legal-content/en/ALL/?uri=CELEX:32008L0050 (accessed on 11 August 2022).
  27. World Health Organization. WHO Global Air Quality Guidelines: Particulate Matter (PM2.5 and PM10), Ozone, Nitrogen Dioxide, Sulfur Dioxide and Carbon Monoxide. 2021. Available online: https://apps.who.int/iris/handle/10665/345329 (accessed on 13 July 2022).
  28. World Bank. Document of the World Bank. Report No: AUS0001229. Western Balkans Regional AQM—Western Balkans Report—AQM in Kosovo. 2019. Available online: https://documents1.worldbank.org/curated/en/214511576520047805/pdf/Air-Pollution-Management-in-Kosovo.pdf (accessed on 30 June 2022).
  29. JFE Techno-Research Corporation. Republic of Kosovo, Ministry of Environment and Spatial Planning, Republic of Kosovo Expert for Air Pollution Control. Final Report. Available online: https://openjicareport.jica.go.jp/pdf/12260998.pdf (accessed on 30 June 2022).
  30. UNFPA Kosovo. The UNFPA Kosovo (UNSCR 1244)* Programme Evaluation Report. (Period Covered by the Evaluation: 2013–2018). Available online: https://web2.unfpa.org/public/about/oversight/evaluations/downloadapprsp.unfpa;jsessionid=02A48B3CA4798E5E1B2D64E06B36D52D?docId=136&annex (accessed on 30 June 2022).
  31. Republic of Kosova. Vjetari Statistikor i Republikës së Kosovës 2021. 2021. Available online: https://ask.rks-gov.net/media/6800/vjetari-statistikor-2021f.pdf (accessed on 30 June 2022).
  32. Osmani, A.A.; Marušic, D.; Halimi, R.; Muharremi, R.; Rupel, V.P. Health care reform in Kosovo. Eurohealth 2017, 23, 20–23. [Google Scholar]
  33. European Environment Agency. Air Quality e-Reporting (AQ e-Reporting) EEA. 2021. Available online: https://www.eea.europa.eu/data-and-maps/data/aqereporting-9 (accessed on 30 June 2022).
  34. KHMI. Data from Monitoring Stations. 2022. Available online: https://airqualitykosova.rks-gov.net/en/reports-for-the-monitoring-stations/ (accessed on 30 June 2022).
  35. AirNow. Available online: https://www.airnow.gov/international/us-embassies-and-consulates/#Kosovo$Pristina (accessed on 30 June 2022).
  36. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020. Available online: https://www.R-project.org/ (accessed on 30 June 2022).
  37. Neuberger, M.; Rabczenko, D.; Moshammer, H. Extended effects of air pollution on cardiopulmonary mortality in Vienna. Atmos. Environ. 2007, 41, 8549–8556. [Google Scholar] [CrossRef]
  38. Neuberger, M.; Moshammer, H.; Rabczenko, D. Acute and subacute effects of urban air pollution on cardiopulmonary emergencies and mortality: Time series studies in Austrian cities. Int. J. Environ. Res. Public Health 2013, 10, 4728–4751. [Google Scholar] [CrossRef]
  39. Hastie, T.; Tibshirani, R. Generalized Additive Models, 1st ed.; Chapman & Hall: London, UK, 1990. [Google Scholar]
  40. Schwartz, J. Harvesting and long term exposure effects in the relation between air pollution and mortality. Am. J. Epidemiol. 2000, 151, 440–448. [Google Scholar] [CrossRef]
  41. Dominici, F.; McDermott, A.; Zeger, S.L.; Samet, J.M. On the use of generalized additive models in time-series studies of air pollution and health. Am. J. Epidemiol. 2002, 156, 193–203. [Google Scholar] [CrossRef]
  42. Dominici, F.; Dermott, A.; Hastie, T. Improved Semiparametric Time Series Models of Air Pollution and Mortality. J. Am. Stat. Assoc. 2004, 99, 938–948. [Google Scholar] [CrossRef]
  43. Katsouyanni, K.; Touloumi, G.; Samoli, E.; Gryparis, A.; Le Tertre, A.; Monopolis, Y.; Rossi, G.; Zmirou, D.; Ballester, F.; Boumghar, A.; et al. Confounding and effect modification in the short-term effects of ambient particles on total mortality: Results from 29 European cities within the APHEA2 project. Epidemiology 2001, 12, 521–531. [Google Scholar] [CrossRef]
  44. Akaike, H. Information theory as an extension of the maximum likelihood principle. In Proceedings of the Second International Symposium on Information Theory, Tsahkadsor, Armenia, 2–8 September 1971; Petrov, B.N., Csaksi, F., Eds.; Akademiai Kiado: Budapest, Hungary, 1973. [Google Scholar]
  45. Maintaining Essential Health Services in Kosovo. Available online: https://hlh.who.int/learning-briefs/action-brief-kosovo---deep-dive (accessed on 29 June 2022).
  46. Liu, J.; Li, Y.; Li, J.; Liu, Y.; Tao, N.; Song, W.; Cui, L.; Li, H. Association between ambient PM2.5 and children’s hospital admissions for respiratory diseases in Jinan, China. Environ. Sci Pollut. Res. Int. 2019, 26, 24112–24120. [Google Scholar] [CrossRef] [PubMed]
  47. Jin, L.; Zhou, T.; Fang, S.; Zhou, X.; Bai, Y. Association of air pollutants and hospital admissions for respiratory diseases in Lanzhou, China, 2014–2019. Environ. Geochem. Health 2022, 6, 1–19. [Google Scholar] [CrossRef] [PubMed]
  48. Ren, Z.; Liu, X.; Liu, T.; Chen, D.; Jiao, K.; Wang, X.; Suo, J.; Yang, H.; Liao, J.; Ma, L. Effect of ambient fine particulates (PM2.5) on hospital admissions for respiratory and cardiovascular diseases in Wuhan, China. Respir. Res. 2021, 22, 128. [Google Scholar] [CrossRef] [PubMed]
  49. Albahar, S.; Li, J.; Al-Zoughool, M.; Al-Hemoud, A.; Gasana, J.; Aldashti, H.; Alahmad, B. Air Pollution and Respiratory Hospital Admissions in Kuwait: The Epidemiological Applicability of Predicted PM2.5 in Arid Regions. Int. J. Environ. Res. Public Health 2022, 19, 5998. [Google Scholar] [CrossRef]
  50. Thomas, D.C. Statistical Methods in Environmental Epidemiology; Oxford University Press: New York, NY, USA, 2009; p. 165f. ISBN 10 0199232903. [Google Scholar]
  51. Rabl, A. Air pollution mortality: Harvesting and loss of life expectancy. J. Toxicol. Environ. Health A 2005, 68, 1175–1180. [Google Scholar] [CrossRef]
  52. Schwartz, J. The Distributed Lag between Air Pollution and Daily Deaths. Epidemiology 2000, 11, 320–326. [Google Scholar] [CrossRef]
  53. Zeger, S.L.; Dominici, F.; Samet, J. Harvesting-resistant estimates of air pollution effects on mortality. Epidemiology 1999, 10, 171–175. [Google Scholar] [CrossRef]
  54. Fung, K.; Krewski, D.; Burnett, R.; Dominici, F. Testing the harvesting hypothesis by time-domain regression analysis. I: Baseline analysis. J. Toxicol. Environ. Health A 2005, 68, 1137–1154. [Google Scholar] [CrossRef]
  55. Schimmel, H.; Murawski, T.J. The relation of air pollution to mortality. J. Occup. Med. 1976, 18, 316–333. [Google Scholar] [CrossRef]
  56. Neuberger, M.; Schimek, M.G.; Horak, F.; Moshammer, H.; Kundi, M.; Frischer, T.; Gomiscek, B.; Puxbaum, H.; Hauck, H. Acute effects of particulate matter on respiratory diseases, symptoms and functions: Epidemiological results of the Austrian Project on Health Effects of Particulate Matter (AUPHEP). Atmos. Environ. 2004, 38, 3971–3981. [Google Scholar] [CrossRef]
Figure 1. Smoothed time course of daily ambulatory visits (left) and stationary hospital admissions (right).
Figure 1. Smoothed time course of daily ambulatory visits (left) and stationary hospital admissions (right).
Ijerph 19 10141 g001
Figure 2. Effect (relative risk and 95% confidence interval) per 10 µg/m³ increase in PM2.5, total effect and daily effects over 7 consecutive days (lag 0–6): (a) ambulatory visits, (b) stationary hospital admissions.
Figure 2. Effect (relative risk and 95% confidence interval) per 10 µg/m³ increase in PM2.5, total effect and daily effects over 7 consecutive days (lag 0–6): (a) ambulatory visits, (b) stationary hospital admissions.
Ijerph 19 10141 g002
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shabani Isenaj, Z.; Berisha, M.; Gjorgjev, D.; Dimovska, M.; Moshammer, H.; Ukëhaxhaj, A. Air Pollution in Kosovo: Short Term Effects on Hospital Visits of Children Due to Respiratory Health Diagnoses. Int. J. Environ. Res. Public Health 2022, 19, 10141. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph191610141

AMA Style

Shabani Isenaj Z, Berisha M, Gjorgjev D, Dimovska M, Moshammer H, Ukëhaxhaj A. Air Pollution in Kosovo: Short Term Effects on Hospital Visits of Children Due to Respiratory Health Diagnoses. International Journal of Environmental Research and Public Health. 2022; 19(16):10141. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph191610141

Chicago/Turabian Style

Shabani Isenaj, Zana, Merita Berisha, Dragan Gjorgjev, Mirjana Dimovska, Hanns Moshammer, and Antigona Ukëhaxhaj. 2022. "Air Pollution in Kosovo: Short Term Effects on Hospital Visits of Children Due to Respiratory Health Diagnoses" International Journal of Environmental Research and Public Health 19, no. 16: 10141. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph191610141

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop