Time Series Study of Air Pollution Health Effects in COPSAC Children

2 Data

2.1 The COPSAC Cohort
2.2 Air Pollutants
2.3 Meteorological Data

Data for this project include the COPSAC cohort, which provided data on the incidence of respiratory symptoms in small children, air pollution data and meteorological data provided by the Danish National Environmental Institute (DMU), and Influenza epidemics data provided by Staten's Serum Institute (SSI).

2.1 The COPSAC Cohort

The COPSAC (COpenhagen Prospective Study on Atopy in Children) cohort was designed to study childhood asthma, eczema and allergy with focus on epidemiological, clinical, cellular and molecular research with consideration to environment and lifestyle factors. COPSAC is unique as it consists of a large group of high-risk children and has long-term continuous data registration. Clinical data are collected prospectively, and biological materials collected into a data bank, which enables par clinical analysis during, before, and after disease development. COPSAC is therefore a unique resource for studying causes and development (course) of asthma, eczema and allergy in children. Details about design of the COPSAC cohort are published elsewhere (Bisgaard, 2004).

The COPSAC cohort consists of 411 children genetically predisposed to atopic illnesses (children of mothers with asthma) living in Denmark. COPSAC children were enrolled in the study at birth and followed until 3 years of age. The first child was enrolled in the cohort on 02.08.1998 and the latest on 29.12.2001. Daily registration of airways symptoms via dairies is carried out by parents for all three years. The presence of the registered symptoms (ICD-10 diagnosis R068) is the outcome in this study. R068 diagnosis includes airways symptoms such as wheezing, apnoeic episode, breath-holding attacks, etc. R068 diagnoses were defined by their starting and finishing date. Most children had R068 diagnoses that lasted more than one day.

For this project, data on incidence of respiratory disease in COPSAC children up to 18 months of age are available. Due to study design where air pollution data are available from a single measuring station located in the centre of Copenhagen (HCØ Institute), three populations of COPSAC children were created according to their home vicinity to the HCØ. By this design it is assumed that air pollution exposure from central monitoring station is most representative exposure for children living in the centre of the city defined by postal code 2450. This, Population 1, is thus the primary study population. Further, Populations 2 and 3 are created, representing those who live in Copenhagen suburbs (2450 < and 2930) and those living in rest of Sealand (2930 < and 5000), as illustrated in Figure 2.1.1. These three populations allow us to study a geographical gradient of centrally measured air pollution's effect on the development of respiratory symptoms in small children, with a hypothesis that this exposure has strongest effect on children living in the centre of the city and this effect gradually diminishing further away from the Copenhagen city centre into the rural areas of Sealand. Geographical distribution of COPSAC children and definition of three populations according to children's' home address postal code can be seen in Table 2.1.1 below. A number of children changed their home address during the first 18 months of life, due to which those children fall into two or all three populations at the different relevant time periods, and the totals do not add up to 411.

Table 2.1.1: Definition of study populations according to postal code of COPSAC children's home address

  Population Definition Postal Code Sample Size
Study Populations Population 1 Inner City Copenhagen 2450 115
Population 2 Copenhagen Suburbs 2450 < and 2930 134
Population 3 Rest of Sealand 2930 < and 5000 186
Not Included in the Study   Odense C 5000 1
  Odense M 5230 1
  Svendborg 5700 1
  Esbjerg 6700 1
  Vejle 7100 2
  Silkeborg 8600 1

Figure 2.1.1: Postal code area definition of the three study populations

Figure 2.1.1: Postal code area definition of the three study populations

Population 1 consists of 115 children who during some of or the whole period of their first 18 months of life lived in Inner City of Copenhagen (postal code 2450). The follow-up period starts with the earliest child's birthrate 16.08.1998 and ends with the latest 18-month follow-up date on 29.06.2003 totaling in 1.779 days. During this period total of 346 new (incident) cases of R068 diagnoses were observed and 2.332 total (prevalent) cases. Table 2.1.2 and Figures 2.1.2 and 2.1.3 below describe Population 1.

Table 2.1.2: Population 1 - Definition of the outcome

  Study Period
16.08.1998 – 29.06.2003
1.779 days
R068 Symptom Duration Incident Cases Prevalent Cases
One Day (startdate = finishdate) 67 67
More than One Day 271 2.265
Unknown Duration (missing finishdate) 8 Missing
Total number of days with R068 346 (mean 0.2) 2.332 (mean=1.3)

Figure 2.1.2: Daily count of incident cases and number of children in COPSAC Cohort in Population 1 during study period (16.08.1998-29.06.2003)

Figure 2.1.2: Daily count of incident cases and number of children in COPSAC Cohort in Population 1 during study period (16.08.1998-29.06.2003)

Figure 2.1.3: Daily count of prevalent cases and number of children in COPSAC Cohort in Population 1 during study period (16.08.1998-29.06.2003)

Figure 2.1.3: Daily count of prevalent cases and number of children in COPSAC Cohort in Population 1 during study period (16.08.1998-29.06.2003)

Population 2 consists of 134 children who during some of or the whole period of their first 18 months of life lived in the suburbs of Copenhagen (2450 < postal code 2930). The follow-up period starts with the earliest child's birthrate 06.10.1998 and ends with the latest 18-month follow-up date on 28.06.2003 totaling in 1.727 days. During this period total of 691 new (incident) cases of R068 diagnoses were observed and 5.663 total (prevalent) cases. Table 2.1.3 and Figures 2.1.4 and 2.1.5 describe Population 2.

Table 2.1.3: Population 2 - Definition of the outcome

  Study Period
06.10.1998 – 28.06.2003
1.727 days
R068 Symptom Duration Incident Cases Prevalent Cases
One Day (startdate = finishdate) 139 139
More than One Day 541 5.524
Unknown Duration (missing finishdate) 11 Missing
Total number of days with R068 691 (mean 0.4) 5.663 (mean=3.3)

Figure 2.1.4: Daily count of incident cases and number of children in COPSAC Cohort in Population 2 during study period (06.10.1998-28.06.2003)

Figure 2.1.4: Daily count of incident cases and number of children in COPSAC Cohort in Population 2 during study period (06.10.1998-28.06.2003)

Figure 2.1.5: Daily count of prevalent cases and number of children in COPSAC Cohort in Population 1 during study period (06.10.1998-28.06.2003)

Figure 2.1.5: Daily count of prevalent cases and number of children in COPSAC Cohort in Population 1 during study period (06.10.1998-28.06.2003)

Population 3 consists of 186 children who during some of or the whole period of their first 18 months of life lived in the rest of Sealand (2930 < postal code 5000). The follow-up period starts with the earliest child's birthrate 02.08.1998 and ends with the latest 18 month follow-up date on 22.06.2003 totaling in 1.786 days. During this period total of 914 new (incident) cases of R068 diagnoses were observed and 7.098 total (prevalent) cases. Table 2.1.4 and Figures 2.1.6 and 2.1.7 describe Population 3.

Table 2.1.4: Population 3 – Definition of the outcome

  Study Period
02.08.1998 – 22.06.2003
1.786 days
R068 Symptom Duration Incident Cases Prevalent Cases
One Day (startdate = finishdate) 177 177
More than One Day 732 6.921
Unknown Duration (missing finishdate) 5 missing
Total number of days with R068 914 (mean 0.5) 7.098 (mean=4.0)

Figure 2.1.5: Daily count of incident cases and number of children in COPSAC Cohort in Population 3 during study period (02.08.1998-22.06.2003)

Figure 2.1.5: Daily count of incident cases and number of children in COPSAC Cohort in Population 3 during study period (02.08.1998-22.06.2003)

Figure 2.1.6: Daily count of prevalent cases and number of children in COPSAC Cohort in Population 3 during study period (02.08.1998-22.06.2003)

Figure 2.1.6: Daily count of prevalent cases and number of children in COPSAC Cohort in Population 3 during study period (02.08.1998-22.06.2003)

2.2 Air Pollutants

Air pollutant data are available as daily averages of air pollutant levels from following measuring stations in Copenhagen: HCØ Institute (city background pollution levels), Jagtvej and H.C.Andersen's Boulevard (HCAB) (street level pollution levels). From all measuring stations measurements are available for the following pollutants: CO (ppm), NOx (ppb), NO2 (ppb), O3 (ppb), PM10 (g/m³), and TON (part./m³). PM10 at street level measuring stations Jagtvej and HCAB was measured by SM200 gravimetric method. PM10 city background level variable (HCØ) is combined with PM10 measurements from HCØ (SM200 gravimetric method) and values extrapolated from PM10 measurements at Jagtvej (beta method), corrected for traffic contribution using NOx measurements from HCØ and Jegtvej and the ratio 0.144*(PM10 /NOx). The extrapolation of PM10 values was done for the study periods where HCØ measurements were missing, to obtain more complete variable. Furthermore, PM2.5 (µg/m³) measurements are available from HCAB measuring station, and PM10 (g/m³) measurements (SM200 gravimetric method) from Lille Valby (rural Zealand pollution levels). Three high PM10 values measured at HCØ were excluded (176.6 g/m³ on 01.01.2000, 248.5 g/m³ on 01.01.2001, and 283.7 g/m³ on 10.08.2000). Furthermore, a PM10 value of 215.5 g/m³ measured at Jagtvej on 01.01.2001 was excluded. Description of air pollutants can be seen in Table 2.2.1 below, and Figures 2.2.1-2.2.7. Pearson correlation coefficients between pollutants measured at the same station can be seen in Tables A.1 - A.3 in Appendix A.

Table 2.2.1: Air pollutant levels during study period

  Study Period (02.08.1998-29.06.2003)
n = 1.793
  Mean SD n Median Range
CO (ppm)
HCØ
Jagtvej
HCAB
0.29 ± 0.11
0.99 ± 0.43
0.81 ± 0.3
1.736
1.737
808
0.28
0.96
0.76
0.09 - 0.97
0.14 - 3.32
0.28 - 2.31
NOx (ppb)
HCØ
Jagtvej
HCAB
15.5 ± 9.3
61.1 ± 31.6
87.3 ± 35.6
1.658
1.736
800
12.9
58.3
81.2
2.9 - 78.9
5.8 - 229.6
14.0 - 263.9
NO2 (ppb)
HCØ
Jagtvej
HCAB
11.9 ± 5.2
23.5 ± 8.2
30.9 ± 8.6
1.658
1.736
800
11.0
23.1
30.2
2.6 - 41.1
3.9 - 62.3
7.6 - 78.0
O3 (ppb)
HCØ
Jagtvej
HCAB
24.6 ± 10.3
16.5 ± 8.3
17.4 ± 8.9
338
1.748
714
25.4
16.1
17.3
2.6 - 48.9
0.8 - 50.3
1.1 - 48.2
PM10 (g/m³)
HCØ
Jagtvej
Lille Valby
27.2 ± 15.7
34.5 ± 16.3
25.0 ± 13.9
1.313
831
754
23.8
31.3
22.3
0.9 - 129.0
3.5 - 88.4
1.1 - 106.1
PM2,5 (g/m³)
HCAB 18.0 ± 9.6 399 15.2 4.6 - 72.2
TON (part./m³) per 100
HCØ
Jagtvej
HCAB
88.0 ± 38.5
234.3 ± 110.2
364.5 ± 176.7
427
192
312
78.4
217.2
320.1
23.0 - 280.8
63.2 - 547.1
60.9 - 1199.6

Figure 2.2.1: CO (ppm) daily levels during study period (02.08.1998-29.06.2003)

Figure 2.2.1: CO (ppm) daily levels during study period (02.08.1998-29.06.2003)

Figure 2.2.2: NOx (ppb) daily levels during study period (02.08.1998-29.06.2003)

Figure 2.2.2: NO<sub>x</sub> (ppb) daily levels during study period (02.08.1998-29.06.2003)

Figure 2.2.3: NO2 (ppb) daily levels during study period (02.08.1998-29.06.2003)

Figure 2.2.3: NO<sub>2</sub> (ppb) daily levels during study period (02.08.1998-29.06.2003)

Figure 2.2.4: O3 (ppb) daily levels during study period (02.08.1998-29.06.2003)

Figure 2.2.4: O<sub>3</sub> (ppb) daily levels during study period (02.08.1998-29.06.2003)

Figure 2.2.5: PM10 (g/m³) daily levels during study period (02.08.1998-29.06.2003)

Figure 2.2.5: PM<sub>10</sub> (g/m³) daily levels during study period (02.08.1998-29.06.2003)

Figure 2.2.6: TON(part./m³)/100 daily levels during study period (02.08.1998-29.06.2003)

Figure 2.2.6: TON(part./m³)/100 daily levels during study period (02.08.1998-29.06.2003)

Figure 2.2.7: PM2.5 (g/m³) daily levels measured at HCAB during study period (02.08.1998-29.06.2003)

Figure 2.2.7: PM<sub>2.5</sub> (g/m³) daily levels measured at HCAB during study period (02.08.1998-29.06.2003)

2.3 Meteorological Data

Meteorological data are available as daily hour averages measured at HCØ and include: wind speed (m/s), temperature (C), relative humidity (%), and global radiation (W/m²). Description of meteorological data can be seen in the Table 2.3.1 and Figures 2.3.1 – 2.3.4 below. Pearson correlation coefficients between meteorological variables can be seen in Table 2.3.2. Pearson correlation coefficients between air pollutants measured at HCØ and weather variables can be seen in Table A.4 in Appendix A.

Gennemsnit m.h.t. vindhastighed, temp og fugtighed er vel uinteressant, her er fx kvartiler vel mere beskrivende.

Table 2.3.1: Meteorological data level during study period

  Study Period (02.08.1998-29.06.2003)
n = 1.793
  Percentilles
Mean SD n 25th 50th 75th
Wind Speed (m/s)
Temperature (C)
Relative Humidity (%)
Global Radiation (W/m²)
4.2 ± 1.5
9.0 ± 6.6
76.0 ± 11.6
111.6 ± 95.0
1.668
1.734
1.733
1.720
3.08
3.62
68.42
25.52
3.94
8.68
77.21
87.91
5.16
14.3
84.58
185.12

Figure 2.3.1: Wind Speed daily averages during the study period (02.08.1998-29.06.2003)

Figure 2.3.1: Wind Speed daily averages during the study period (02.08.1998-29.06.2003)

Figure 2.3.2: Temperature daily averages during the study period (02.08.1998-29.06.2003)

Figure 2.3.2: Temperature daily averages during the study period (02.08.1998-29.06.2003)

Figure 2.3.3: Relative Humidity daily averages during the study period (02.08.1998-29.06.2003)

Figure 2.3.3: Relative Humidity daily averages during the study period (02.08.1998-29.06.2003)

Figure 2.3.4: Global Radiation daily averages during the study period

Figure 2.3.4: Global Radiation daily averages during the study period

Table 2.3.2: Correlation of meteorological variables during the study period in Copenhagen (02.08.1998 - 29.06.2003)

  Wind Speed Temperature Global
Radiation
Relative
Humidity
Wind Speed 1.00 -0.16* -0.25* 0.11*
Temperature   1.00 0.67* -0.35*
Global Radiation     1.00 -0.75*
Relative Humidity       1.00

* p < 0.01 - significance level for the Pearson correlation coefficients

In the Table 2.3.2 above it can be seen that there is statistically significant positive association between wind speed and relative humidity, and temperature and global radiation. There is negative and statistically significant correlation between wind speed and temperature, wind speed and global radiation, temperature and relative humidity and global radiation and relative humidity. This table shows that all meteroloigcal variables are mutually correlated, which implies that fitting them in the model together could cause some colinearity. Therefore, we chose only temperature in the final model, as the strongest predictor of the incident respiratory cases in small children.

 



Version 1.0 Maj 2005, © Danish Environmental Protection Agency