space space   space
space

Environmental Health and 
Biostatistics and Computing Groups
Department of Community Medicine
The University of Hong Kong

Environmental Health and
Biostatistics and Computing Groups
Department of Community Medicine
The University of Hong Kong

Project Team

Dr CM Wong
(Data analysis and report writing)

Mr Stefan Ma
(Computation and statistical advice)

   
Professor AJ Hedley
(Head of Department)
Professor TH Lam
(Epidemiological advice)

CONTENTS

Executive Summary

1.0 Background and introduction
   
2.0 Scope and objectives
   
3.0

Materials and methods

3.1 Study design
3.2 Databases
3.3 Statistical modelling
   
4.0

Findings

4.1 Descriptive statistics
4.2 Statistical modelling
   
5.0

Discussion

5.1 Validity and reliability of the models
5.2 Summary of findings
5.3 Comparison with APHEA studies
5.4 Comparison with CUHK study
5.5 Limitations of the Hong Kong study
   
6.0

Conclusions

6.1 Critique of CUHK study

References

Figures

Basic Tables

  Appendix A Basic Tables
   
  Appendix B Tables with estimates per 100 ug/m3 changes

Operation Manual


Executive Summary

Background and objectives

Valuable indicators of the possible benefits of environmental management and control can be obtained by extrapolation from analysis carried out in other locations; but governments and local regulatory agencies are usually unable to draft or implement effective legislation without relevant local information to support their proposals. A study had been commissioned to the Chinese University of Hong Kong (CUHK) by the Environmental Protection Department, to evaluate the acute health effects of air pollution, using data for 1994-96 as a first attempt towards utilizing local intelligence. This study is a follow-up of the first study, aiming to validate the methods and results in the first study.

Methods

A series of daily hospital admissions and deaths in 1995, 1996 and the first half of 1997, due to respiratory and circulatory diseases, were obtained from routinely collected data and were analysed using Poisson regression with adjustment for overdispersion and long term effects of covariates

(including trend, seasonality, weekdays, holidays, after holidays, temperature and humidity). The health effects due to daily pollutant concentrations of sulphur dioxide (SO2), nitrogen dioxide (NO2), ozone (O3) and respirable suspended particulates (RSP) were then estimated and compared with those obtained from other similar studies.

Findings

(a)

For hospital admissions

Air pollution was found to have an effect on circulatory and respiratory diseases combined and separately (relative risk, RR=1.03-1.10; p<0.084) for all ages; on circulatory admissions (RR=1.05-1.10; p<0.001) for the 65 or above age group; on respiratory admissions, the effects of which appeared to be j-shaped from the younger to the older age groups; on asthma (RR=1.10-1.16; p<0.018, except SO2); on chronic obstructive pulmonary disease (RR=1.08-1.14; p<0.0001); and on ischaemic heart disease (RR=1.04-1.09; p<0.051, except SO2 and RSP).

   
(b)

For hospital deaths

Both NO2 and O3 were positively associated with circulatory and respiratory diseases combined and separately (for NO2: RR= 1.10-1.14, p<0.038; and for O3: RR=1.07-1.22, p<0.010).

   
(c)

Validation and composite score

The above estimates were consistent with and in between those obtained from similar studies using the European (APHEA) approach overseas and in the CUHK. But in addition, a composite score was derived from the four pollutants and was found to provide consistent estimates for all the health outcomes under study in all ages (RR=1.04-1.11, p<0.098, except hospital deaths due to circulatory diseases).

Conclusions

Routine hospital morbidity and mortality data, air pollution and meteorological observations can be utilized to provide information for the estimation of acute health effects of air pollution. Environmental management and control should and could take into account health effects of air pollution based on locally derived information. However the processes are errors prompted (as large and complex data sets are involved) and are vulnerable to the misuse of health and other parameters and to misinterpretation of the results as it involves knowledge from several different fields. A team approach with expertise from epidemiological, environmental, statistical and computational professionals is required.

Back to Contents


1.0 Background and introduction

1.1 The contribution of epidemiological studies to the process of environmental management and control of health hazards is well established world-wide. Valuable indicators of possible benefits can be obtained by extrapolation from analyses carried out in other locations but, in general, governments and local regulatory agencies are unable to draft or implement effective legislation without relevant local information to support their proposals.
   
1.2

Several epidemiological studies have now shown an association between particulate air pollution and exacerbations of illness in individuals with respiratory disease and also increases in the numbers of deaths from cardiovascular and respiratory disease, particularly in the elderly. New hypotheses have been advanced to postulate mechanisms underlying these observed effects (Seaton et al 1995).1

Respirable particulates (RSP) with an aerodynamic diameter of <10 um(or particles measured as black smoke by the smoke stain method) comprise the principal pollutant associated with these findings. In addition RSP sulphate and sulphur dioxide concentrations, are reported to be associated with all causes mortality and respiratory mortality in recent studies from the USA (Pope et al 1995)2 and the UK (Anderson et al 1996)3 respectively. A recent review in the United Kingdom concludes that the associations between daily concentrations of particles and acute health effects principally reflect a causal relationship (Committee on the Medical Effects of Air Pollutants 1995).4 After a lengthy scientific review, the USEPA determined that new standard should be added for particulates less than 2.5 ? of aerodynamic in size and the welfare-base standards were also revised by making them identical to the health-based standards.5

In Hong Kong we have found that SO2, RSP and SO4 concentrations are associated with excess risks for symptoms of cough, phlegm and wheeze and also bronchial hyper-responsiveness (by histamine challenge test) in primary school children (Hedley et al 1993; Peters et al 1996; Tam et al 1994; Wong et al 1998).6,7,8,9 However in the London study the strongest association with daily mortality was for ozone. The effects of ozone and black smoke were independent of the effects of other pollutants.

Evidence on the risks associated with other pollutants is variable and less consistent. In the recent London study (Anderson et al 1996)3 the NO2 (1 hour maximum) was associated with all causes mortality and cardiovascular mortality; however a negative effect was seen for respiratory mortality. A significant positive effect on mortality was seen for SO2 and all cause mortality in the warm season period.

Several authors point to the complex between-season covariation of several pollutants, which is sometimes negative and at other times positive.

   
1.3 The published literature in this field is growing rapidly. The Department of Community Medicine is monitoring this through different databases and will aim to carry out the analysis using a state-of-the-art approach. This will enhance the utility of the outputs and ensure comparability as far as possible with studies in other countries (Katsouyanni, Schwartz, Spix et al 1996).10

Back to Contents


2.0 Scope and objectives

2.1 to examine the variation of daily air pollution data i.e. 24 hour average for sulphur dioxide, nitrogen dioxide and respirable suspended particulates and 8 hour average for ozone, among the various monitoring stations in Hong Kong for the years 1995-1996, as available;
   
2.2 to investigate the availability and the use of the various health outcome measures including data on hospital admissions and hospital deaths due to respiratory and circulatory problems collected routinely in the Hong Kong hospitals;
   
2.3 to investigate the short-term effects of the air pollutants considered in 2.1 (in the same day and one or more days lagged) individually and compositely on some of the health outcome measures considered in 2.2 above, with adjustment for seasonal variations, secular trends as well as meteorological conditions including temperature and humidity;
   
2.4 to validate and update models developed earlier and to develop a mechanism for the use and maintenance of the model for continuous study by the Environmental Protection Department.

Back to Contents


3.0 Materials and methods

3.1

Study design

It was an ecological study utilizing routinely collected hospital admission data, air pollutant concentration data and weather data by the Hospital Authority, Environmental Protection Department and the Observatory respectively. Variations in the daily number of hospital admissions due to circulatory and respiratory diseases were studied, and their relationships with each of the pollutants were modelled to assess the effects of air pollution on health after adjustment for time trends, seasonality, weather conditions and some other factors including days of week, holidays and days after holidays.

This study follows a previous one performed by the Chinese University of Hong Kong (CUHK)11 which followed the general approach of the protocol of the APHEA (a European approach using epidemiological time series data), developed within the frame of the EC Environment 1991-94 Programme. However data for some disease categories (ICD Rubrics) included in the APHEA protocol were not analysed in the CUHK study. The data sets and disease categories used for this new study are shown in Table 1 below. Those categories which were excluded or missing from the CUHK study are indicated in the table.

Table 1: Number of hospital admissions by disease groups

Disease groups 1995 1996 1997*
I. Diseases of the Circulatory System (ICD9 390-459)12:
  Acute rheumatic fever (390-392)** 32 - 12 11
  Chronic rheumatic heart disease (393-398)** 1,572 - 1,626 672
  Hypertensive disease (401-405)*** 4,396 (0) 4,319 2,000
  Ischaemic heart disease (410-414) 12,281 (11,884) 13,741 6,560
  Disease of pulmonary circulation (415-417) 291 (289) 274 144
  Other forms of heart disease (420-429) 15,494 (15,549) 17,567 8,932
  Cerebrovascular disease (430-438) 12,474 (10,224) 13,326 6,826
  Diseases of arteries, arterioles and capillaries (440-448)** 1,909 (1,300) 2,144 944
  Diseases of veins and lymphatics, and other diseases of circulatory system (451-459)** 5,591 - 6,448 3,240
  Sub-total: 54,040 (39,246) 59,457 29,329
II. Diseases of the Respiratory System (ICD9 460-519):
  Acute respiratory infections (460-466) 13,845 (13,871) 16,906 9,968
  Other diseases of upper respiratory tract (470-478)** 3,005 (2,228) 3,650 1,641
  Pneumonia and influenza (480-487) 12,567 (12,574) 14,944 7,648
  Chronic obstructive pulmonary disease and allied conditions (490-496) 25,330 (25,357) 28,344 13,567
  Pneumoconioses and other lung diseases due to external agents (500-508)*** 386 (0) 471 275
  Other diseases of respiratory system (510-519)** 11,358 - 13,478 7,733
  Sub-total: 66,491 (54,030) 77,793 40,832
I & II Total: 120,531 (93,276) 137,250 70,161

* first half year
** data for ICD9 390-392, 393-398, 446-448, 451-459, 470 and 510-519 were not included for analysis in the CUHK Final Report
*** data for ICD9 401-405 and 500-508 were missing in the CUHK data files and were not analysed in the CUHK Final Report
( ) CUHK data in brackets
- not included for analysis in the CUHK study

In order not to exclude categories which might show effects from air pollution (e.g. hospital admissions for cardiovascular diseases under ICD9 390 - 429 which were found to be related to respirable suspended particulates and carbon monoxide by Schwartz (1997)13 but not all included in the CUHK study) and to ensure comparability to the results of those studies which follow the APHEA protocol (Bacharova 1996; Schouten 1996),14,15 we decided to follow strictly the categories recommend by the APHEA protocol for this study.

Data for 1995-1996 (generated and cleaned by the Department of Community Medicine, the University of Hong Kong in this study) were used in establishing the statistical models and in estimating the effects; and data for first half year of 1997 were used for validation of the established models. It was also decided that data for the year 1994, which were used by the CUHK group, should not be used in this study because the data quality is not consistent with that of the 1995-96 data set to be used in this study. In 1994 only 3 of the 12 hospitals under study had already adopted the MRAS database; but the number increased to 7 in 1995. Besides, the percentages of valid daily data for air pollutant concentrations from the monitoring stations were lower in 1994 than those in the other years.

   
3.2

Databases

Hospital admission data: Hospitals included for generation of hospital admissions were the publicly funded hospitals (accounting for 90% of hospital beds in Hong Kong) which either had an accident and emergency department, or was a referral base from the accident and emergency department of another nearby hospital (9) or had a 24-hour outpatient department (2). One other hospital which was the only hospital in the most polluted district in Hong Kong, was also included. All hospitals should have a computerized system for inputting and retrieval of patient data. The hospitals included in the study, together with the type of information system they were using and an indication as to whether they had an A & E department, were listed in Table 2 below:

Table 2: List of HA hospitals included in the study

Hospital Information system* Whether having A&E
1. Kong Wah Hospital (KWH) MRAS Yes
2. Our Lady of Maryknoll Hospital (OLM) IPAS No
3. Princess Margaret Hospital (PMH) IPAS/MRAS Yes
4. Pok Oi Hospital (POH) IPAS Yes
5. Prince of Wales Hospital (PWH) MRAS Yes
6. Pamela Youde Nethersole Hospital (PYN) MRAS Yes
7. Queen Elizabeth Hospital (QEH) MRAS Yes
8. Queen Mary Hospital (QMH) IPAS/MRAS Yes
9. Ruttonjee Hospital (RH) MRAS No
10. Tuen Mun Hospital (TMH) MRAS Yes
11. United Christian Hospital (UCH) MRAS Yes
12. Yan Chai Hospital (YCH) IPAS/MRAS Yes

* IPAS - Integrated Patient Administrative System
MRAS - Medical Records Abstracting System

Patients admitted to hospitals between 1.1.1995 and 30.6.1997 with data on: dates of admission and discharge, socio-demographic information (age, gender, marital status, ethnic group, district of residence, pseudo identifier of patient), admission source, discharge diagnosis in ICD9 codes and discharge status were retrieved from the databases for each of the hospitals under study. In order to validate the completeness of the retrieved data, the total number of inpatients for the period 1.4.1995 to 31.3.1996, for each hospital, were compared with those reported in the Hospital Authority Statistical Report 1995/96. When there were big differences between the two, the Hospital Authority Information Technology Department was ask for an explanation and we then revised the databases if necessary until they were reasonably close to each other. This process took almost a half year to complete. The data are shown in Table 3 below:

Table 3: Comparison of total hospital discharges between HKU data obtained from the HA IPAS/MRAS and those from HA Statistical Report 95/96

 

1.4.1995 - 31.3.1996#

 

HKU

HA

Kong Wah Hospital (KWH)

63,583

63,583

Our Lady of Maryknoll Hospital (OLM)

8,564

8,564

Princess Margaret Hospital (PMH)

104,149

104,149

Pok Oi Hospital (POH)

11,432

11,432

Prince of Wales Hospital (PWH)

104,679

104,679

Pamela Youde Nethersole Hospital (PYN)

66,672@

61,263

Queen Elizabeth Hospital (QEH)

112,557

112,559

Queen Mary Hospital (QMH)

87,308

87,313

Ruttonjee Hospital (RH)

16,727

16,727

Tuen Mun Hospital (TMH)

83,117

83,118

United Christian Hospital (UCH)

51,266

51,293

Yan Chai Hospital (YCH)

39,508

39,508

Total:

749,562

744,188

# total hospital patient discharges
@ The excess 5409 cases in the HKU data set were day cases which could not be excluded when the Information Technology Department generated the data set due to missing of the identifier.

The total data sets were then extracted for circulatory (ICD9 390-459) and respiratory (ICD9 460-519) diseases. The numbers of hospital admissions by disease groups in the three years were as shown in the previous Table 1 and subset of specific disease categories are shown in Table 4 below:

Table 4: Number of hospital admissions by specific diseases

Disease

1995

1996

1997*

Asthma (ICD-9 493)

8,682

9,672

3,803

Chronic obstructive pulmonary disease (ICD-9 490-496)

25,330

28,344

13,567

Ischaemic heart disease (ICD-9 410-414)

12,281

13,741

6,560

* first half year

Acute myocardial infarction was not analysed, as had been done by the CUHK, because diagnosis for the disease has been changing and subject to misclassification over the past years.

Pollutant concentration data: Pollutant concentration data in CD-ROM were made available by the Air Services Group of the Environmental Protection Department with hourly data from all monitoring stations in Hong Kong. The following stations in various urban, suburban and industrial areas were included in the study:

Table 5: List of air pollution monitoring stations

Station

Sampling height

Above ground

Date start operation

1. Central / West (C/W)

78m

18m (4th floor)

11/83

2. Kwai Chung (KC)

82m

25m (6th floor)

7/88

3. Kwun Tong (KT)

34m

25m (6th floor)

7/83

4. Sham Shui Po (SSPO)

21m

17m (4th floor)

7/84

5. Shatin (ST)

27m

21m (5th floor)

7/91

6. Tai Po (TP)

31m

25m (6th floor)

2/90

7. Tsuen Wan (TW)

21m

17m (4th floor)

8/88


Two stations, one in Yuen Long (in suburban area) due to the extent of missing data and one in Mong Kok (in urban area) which provided concentrations measured on ground level, were excluded from the study. The pollutants included in this study are in Table 6 below:

Table 6: List of pollutants used in the study

Pollutant Unit
1. Nitrogen Dioxide {24-hr} (NO2)  
2. Sulphur Dioxide {24-hr} (SO2) All in micrograms/
3. Respirable Suspended Particulates {24-hr} (TEOM)* cubic metre
4. Ozone {9.00 am - 5.00 pm} (O3)  

* TEOM - Tapered Element Oscillating Microbalance, an instrument for the continuous measurement of particulates matter in air

In order to maintain consistency and similar quality standards in the data, missing data were defined and replaced in accordance with the APHEA recommendations. The guidelines were slightly modified to suit local situations and these are described in Table 7 below:

Table 7: Definitions of and methods for and replacement of missing daily data

Procedure Computation
(a) Define non-missing daily data on a particular day (i)

For SO2, NO2 and RSP, if number of non-missing hourly data in that day 318, it will be defined as non-missing.

(ii) For O3 (9 am - 5 pm), if number of non-missing hourly data in that day (during the 8 hour interval) 36, it will be defined as non-missing
(b) Exclude pollutant from a station for further analysis (i)

For each of SO2, NO2 and O3, if the proportion of non-missing daily data in a station over the study period <75%, it will be excluded from the analysis.

(ii) For RSP, if the proportion of non-missing daily data in a station over the study period <67%, it will be excluded from the analysis.#
(c) Compute non-missing daily data   Mean of non-missing hourly data in that day
(d)

Define seasonal i

(i = 1, 2, 3, 4) in a particular year (December include in the next coming year)

 

w(i)=Mean for non-missing daily data in the station in season i / Mean for non-missing daily data in all the other stations in season i

1 = December - February

2 = March - May

3 = June - August

4 = September - November

(e) Define a weight w(i) for a station in a particular season i of the year  

i = 1, 2, 3, 4

w(i) will be missing if the proportion of non-missing daily data in either the above numerator or denominator are less than 75% for SO2, NO2 and O3 and less than 67% for RSP.#

(f) Define missing daily data in a particular day   Data in the day not regarded as non-missing according to (a) were missing data
(g) Replace missing daily data in a particular day during a particular season i for a particular station   Mean of all non-missing daily data over all the other stations multiply by a seasonal weight w(i) defined in (e)

Note: RSP measured by TEOM.

# A minimum of 67% non-missing daily data was used as criteria for inclusion of a pollutant in the analysis instead of 75%. This was set in the computer programme at the beginning of the study when the HKU was trying to adopt a similar procedure as that used by the CUHK. However this might not be necessary for data after 1994 as the data were more complete.

The data in monthly averages were comparable to those in the EPD 1995 and 1996 statistical reports. The percentage of data valid after replacement are shown in the following Table 8:

Table 8: Percentage of valid daily measures of air pollutants by stations in the study (1995-96 and first half of 1997)

  Station#
Pollutant C/W KC KT SSPO ST TP TW
NO2 (24-hr average)
1995 92.60

92.33

95.89 89.04 97.26 96.71 95.34
1996 93.44 95.90 94.81 90.16 95.63 94.54 93.17
1997* 94.48 93.37 87.85 97.24 96.69 95.03 86.74
SO2 (24-hr average)
1995 98.08 96.71 97.26 93.42 99.45 - 95.07
1996 97.54 95.08 98.36 89.07 98.91 - 97.54
1997* 96.69 98.90 98.34 98.90 98.90 - 85.08
TEOM (24-hr average)
1995 85.21 94.79 79.18 - 87.40 - 94.25
1996 98.63 97.81 91.80 - 87.43 - 86.07
1997* 93.92 99.45 92.82 - 96.13 - 86.19
O3 (8-hr average)
1995 94.25 93.15 - - - - -
1996 92.35 94.54 - - - - -
1997* 91.16 96.69 - - - - -

* the data are available only for first half year of 1997
- not available
# Abbreviations referred to Table 5.

A daily mean concentration