Spatial variation of particle number concentration in school microscale environment

There is significant toxicological evidence of the effects of ultrafine particles (<100nm) on human health (WHO 2005). Studies show that the number concentration of particles has been associated with adverse human health effects (Englert 2004). This work is part of a major study called ‘Ultrafine Particles form Traffic Emissions and Children’s Health’ (UPTECH), which seeks to determine the effect of the exposure to traffic related ultrafine particles on children’s health in schools (http://www.ilaqh.qut.edu.au/Misc/UPT ECH%20Home.htm). Quantification of spatial variation of particle number concentration (PNC) in a microscale environment and identification of the main affecting parameters and their contribution levels are the main aims of this analysis.

Spatial variation of particle number concentration in school microscale environments and its impact on exposure assessment

It has not yet been established whether the spatial variation of particle number concentration (PNC) within a microscale environment can have an effect on exposure estimation results. In general, the degree of spatial variation within microscale environments remains unclear, since previous studies have only focused on spatial variation within macroscale environments. The aims of this study were to determine the spatial variation of PNC within microscale school environments, in order to assess the importance of the number of monitoring sites on exposure estimation. Furthermore, this paper aims to identify which parameters have the largest influence on spatial variation, as well as the relationship between those parameters and spatial variation. Air quality measurements were conducted for two consecutive weeks at each of the 25 schools across Brisbane, Australia. PNC was measured at three sites within the grounds of each school, along with the measurement of meteorological and several other air quality parameters. Traffic density was recorded for the busiest road adjacent to the school. Spatial variation at each school was quantified using coefficient of variation (CV). The portion of CV associated with instrument uncertainty was found to be 0.3 and therefore, CV was corrected so that only non-instrument uncertainty was analysed in the data. The median corrected CV (CVc) ranged from 0 to 0.35 across the schools, with 12 schools found to exhibit spatial variation. The study determined the number of required monitoring sites at schools with spatial variability and tested the deviation in exposure estimation arising from using only a single site. Nine schools required two measurement sites and three schools required three sites. Overall, the deviation in exposure estimation from using only one monitoring site was as much as one order of magnitude. The study also tested the association of spatial variation with wind speed/direction and traffic density, using partial correlation coefficients to identify sources of variation and non-parametric function estimation to quantify the level of variability. Traffic density and road to school wind direction were found to have a positive effect on CVc, and therefore, also on spatial variation. Wind speed was found to have a decreasing effect on spatial variation when it exceeded a threshold of 1.5 (m/s), while it had no effect below this threshold. Traffic density had a positive effect on spatial variation and its effect increased until it reached a density of 70 vehicles per five minutes, at which point its effect plateaued and did not increase further as a result of increasing traffic density.

Spatial variation of air pollutants in a multilevel office building

Within-building spatial variability of indoor air quality may influence substantially the reliability of human exposure assessments based on single point samples, but have hitherto been little studied. To investigate and understand the within-building spatial variation of air pollutants, field measurements were conducted in a 7 level office building in Brisbane, Australia. The building consists of 3 sections (A side, Meddler and B side).

Characteristics of spatial variation, size distribution, formation and growth of particles in urban environments

This thesis is the first comprehensive study of important parameters relating to aerosols' impact on climate and human health, namely spatial variation, particle size distribution and new particle formation. We determined the importance of spatial variation of particle number concentration in microscale environments, developed a method for particle size parameterisation and provided knowledge about the chemistry of new particle formation. This is a significant contribution to our understanding of processes behind the transformation and dynamics of urban aerosols. This PhD project included extensive measurements of air quality parameters using state of the art instrumentation at each of the 25 sites within the Brisbane metropolitan area and advanced statistical analysis.

Heat-related morbidity in Brisbane, Australia : spatial variation and area-level predictors

Background: Extreme heat is a leading weather-related cause of illness and death in many locations across the globe, including subtropical Australia. The possibility of increasingly frequent and severe heat waves warrants continued efforts to reduce this health burden, which could be accomplished by targeting intervention measures toward the most vulnerable communities. Objectives: We sought to quantify spatial variability in heat-related morbidity in Brisbane, Australia, to highlight regions of the city with the greatest risk. We also aimed to find area-level social and environmental determinants of high risk within Brisbane. Methods: We used a series of hierarchical Bayesian models to examine city-wide and intracity associations between temperature and morbidity using a 2007–2011 time series of geographically referenced hospital admissions data. The models accounted for long-term time trends, seasonality, and day of week and holiday effects. Results: On average, a 10°C increase in daily maximum temperature during the summer was associated with a 7.2% increase in hospital admissions (95% CI: 4.7, 9.8%) on the following day. Positive statistically significant relationships between admissions and temperature were found for 16 of the city’s 158 areas; negative relationships were found for 5 areas. High-risk areas were associated with a lack of high income earners and higher population density. Conclusions: Geographically targeted public health strategies for extreme heat may be effective in Brisbane, because morbidity risk was found to be spatially variable. Emergency responders, health officials, and city planners could focus on short- and long-term intervention measures that reach communities in the city with lower incomes and higher population densities, including reduction of urban heat island effects.

Space-time variation of malaria incidence in Yunnan province, China

Abstract Background Understanding spatio-temporal variation in malaria incidence provides a basis for effective disease control planning and monitoring. Methods Monthly surveillance data between 1991 and 2006 for Plasmodium vivax and Plasmodium falciparum malaria across 128 counties were assembled for Yunnan, a province of China with one of the highest burdens of malaria. County-level Bayesian Poisson regression models of incidence were constructed, with effects for rainfall, maximum temperature and temporal trend. The model also allowed for spatial variation in county-level incidence and temporal trend, and dependence between incidence in June–September and the preceding January–February. Results Models revealed strong associations between malaria incidence and both rainfall and maximum temperature. There was a significant association between incidence in June–September and the preceding January–February. Raw standardised morbidity ratios showed a high incidence in some counties bordering Myanmar, Laos and Vietnam, and counties in the Red River valley. Clusters of counties in south-western and northern Yunnan were identified that had high incidence not explained by climate. The overall trend in incidence decreased, but there was significant variation between counties. Conclusion Dependence between incidence in summer and the preceding January–February suggests a role of intrinsic host-pathogen dynamics. Incidence during the summer peak might be predictable based on incidence in January–February, facilitating malaria control planning, scaled months in advance to the magnitude of the summer malaria burden. Heterogeneities in county-level temporal trends suggest that reductions in the burden of malaria have been unevenly distributed throughout the province.

Airborne particle concentrations at schools measured at different spatial scales

Potential adverse effects on children health may result from school exposure to airborne particles. To address this issue, measurements in terms of particle number concentration, particle size distribution and black carbon (BC) concentrations were performed in three school buildings in Cassino (Italy) and its suburbs, outside and inside of the classrooms during normal occupancy and use. Additional time resolved information was gathered on ventilation condition, classroom activity, and traffic count data around the schools were obtained using a video camera. Across the three investigated school buildings, the outdoor and indoor particle number concentration monitored down to 4 nm and up to 3 m ranged from 2.8×104 part cm-3 to 4.7×104 part cm-3 and from 2.0×104 part cm-3 to 3.5×104 part cm-3, respectively. The total particle concentrations were usually higher outdoors than indoors, because no indoor sources were detected. I/O measured was less than 1 (varying in a relatively narrow range from 0.63 to 0.74), however one school exhibited indoor concentrations higher than outdoor during the morning rush hours. Particle size distribution at the outdoor site showed high particle concentrations in different size ranges, varying during the day; in relation to the starting and finishing of school time two modes were found. BC concentrations were 5 times higher at the urban school compared with the suburban and suburban-to-urban differences were larger than the relative differences of ultrafine particle concentrations.

Phenotyping cotton ovule fibre initiation with spatial statistics

Yield in cultivated cotton (Gossypium spp.) is affected by the number and distribution of fibres initiated on the seed surface but, apart from simple statistical summaries, little has been done to assess this phenotype quantitatively. Here we use two types of spatial statistics to describe and quantify differences in patterning of cotton ovule fibre initials (FI). The following five different species of Gossypium were analysed: G. hirsutum L., G. barbadense L., G. arboreum, G. raimondii Ulbrich. and G. trilobum (DC.) Skovsted. Scanning electron micrographs of FIs were taken on the day of anthesis. Cell centres for fibre and epidermal cells were digitised and analysed by spatial statistics methods appropriate for marked point processes and tessellations. Results were consistent with previously published reports of fibre number and spacing. However, it was shown that the spatial distributions of FIs in all of species examined exhibit regularity, and are not completely random as previously implied. The regular arrangement indicates FIs do not appear independently of each other and we surmise there may be some form of mutual inhibition specifying fibre-initial development. It is concluded that genetic control of FIs differs from that of stomata, another well studied plant idioblast. Since spatial statistics show clear species differences in the distribution of FIs within this genus, they provide a useful method for phenotyping cotton. © CSIRO 2007.

Regional variation in social isolation amongst older Australians

Regional studies globally have a strong focus on understanding the causes of variation in the economic performance and wellbeing of regions and this emphasis acknowledges that the strength of the local or regional economy plays a determinant role in shaping quality of life. Regional research has been less active in considering spatial variation in other factors that are critical to individual and societal wellbeing. For example, the regional studies community has been absent from the debate on the social determinants of health and how these influences vary spatially. This paper considers the results of a cross sectional survey of Australians aged 65 and over that focussed on social connections and wellbeing. It examines regional variations in the incidence of social isolation within the older population. It finds that while the incidence of self-reported social isolation amongst older persons is broadly consistent with earlier studies, it demonstrates a spatial patterning that is unexpected. The paper considers community-building activities in addressing the impacts of social isolation, including the role of urban design, and suggests that there is a need to supplement the national overview presented there through more detailed studies focussed on individual localities.

Long-term trends in fine particle number concentrations in the urban atmosphere of Brisbane : the relevance of traffic emissions and new particle formation

The measurement of submicrometre (< 1.0 m) and ultrafine particles (diameter < 0.1 m) number concentration have attracted attention since the last decade because the potential health impacts associated with exposure to these particles can be more significant than those due to exposure to larger particles. At present, ultrafine particles are not regularly monitored and they are yet to be incorporated into air quality monitoring programs. As a result, very few studies have analysed their long-term and spatial variations in ultrafine particle concentration, and none have been in Australia. To address this gap in scientific knowledge, the aim of this research was to investigate the long-term trends and seasonal variations in particle number concentrations in Brisbane, Australia. Data collected over a five-year period were analysed using weighted regression models. Monthly mean concentrations in the morning (6:00-10:00) and the afternoon (16:00-19:00) were plotted against time in months, using the monthly variance as the weights. During the five-year period, submicrometre and ultrafine particle concentrations increased in the morning by 105.7% and 81.5% respectively whereas in the afternoon there was no significant trend. The morning concentrations were associated with fresh traffic emissions and the afternoon concentrations with the background. The statistical tests applied to the seasonal models, on the other hand, indicated that there was no seasonal component. The spatial variation in size distribution in a large urban area was investigated using particle number size distribution data collected at nine different locations during different campaigns. The size distributions were represented by the modal structures and cumulative size distributions. Particle number peaked at around 30 nm, except at an isolated site dominated by diesel trucks, where the particle number peaked at around 60 nm. It was found that ultrafine particles contributed to 82%-90% of the total particle number. At the sites dominated by petrol vehicles, nanoparticles (< 50 nm) contributed 60%-70% of the total particle number, and at the site dominated by diesel trucks they contributed 50%. Although the sampling campaigns took place during different seasons and were of varying duration these variations did not have an effect on the particle size distributions. The results suggested that the distributions were rather affected by differences in traffic composition and distance to the road. To investigate the occurrence of nucleation events, that is, secondary particle formation from gaseous precursors, particle size distribution data collected over a 13 month period during 5 different campaigns were analysed. The study area was a complex urban environment influenced by anthropogenic and natural sources. The study introduced a new application of time series differencing for the identification of nucleation events. To evaluate the conditions favourable to nucleation, the meteorological conditions and gaseous concentrations prior to and during nucleation events were recorded. Gaseous concentrations did not exhibit a clear pattern of change in concentration. It was also found that nucleation was associated with sea breeze and long-range transport. The implications of this finding are that whilst vehicles are the most important source of ultrafine particles, sea breeze and aged gaseous emissions play a more important role in secondary particle formation in the study area.

Use of Bayesian geostatistical prediction to estimate local variations in Schistosoma haematobium infection in western Africa

Objective We aimed to predict sub-national spatial variation in numbers of people infected with Schistosoma haematobium, and associated uncertainties, in Burkina Faso, Mali and Niger, prior to implementation of national control programmes. Methods We used national field survey datasets covering a contiguous area 2,750 × 850 km, from 26,790 school-aged children (5–14 years) in 418 schools. Bayesian geostatistical models were used to predict prevalence of high and low intensity infections and associated 95% credible intervals (CrI). Numbers infected were determined by multiplying predicted prevalence by numbers of school-aged children in 1 km2 pixels covering the study area. Findings Numbers of school-aged children with low-intensity infections were: 433,268 in Burkina Faso, 872,328 in Mali and 580,286 in Niger. Numbers with high-intensity infections were: 416,009 in Burkina Faso, 511,845 in Mali and 254,150 in Niger. 95% CrIs (indicative of uncertainty) were wide; e.g. the mean number of boys aged 10–14 years infected in Mali was 140,200 (95% CrI 6200, 512,100). Conclusion National aggregate estimates for numbers infected mask important local variation, e.g. most S. haematobium infections in Niger occur in the Niger River valley. Prevalence of high-intensity infections was strongly clustered in foci in western and central Mali, north-eastern and northwestern Burkina Faso and the Niger River valley in Niger. Populations in these foci are likely to carry the bulk of the urinary schistosomiasis burden and should receive priority for schistosomiasis control. Uncertainties in predicted prevalence and numbers infected should be acknowledged and taken into consideration by control programme planners.

Long-term exposure to gaseous air pollutants and cardio-respiratory mortality in Brisbane, Australia

Air pollution is ranked by the World Health Organisation as one of the top ten contributors to the global burden of disease and injury. Exposure to gaseous air pollutants, even at a low level, has been associated with cardiorespiratory diseases (Vedal, Brauer et al. 2003). Most recent epidemiological studies of air pollution have used time-series analyses to explore the relationship between daily mortality or morbidity and daily ambient air pollution concentrations based on the same day or previous days (Hajat, Armstrong et al. 2007). However, most of the previous studies have examined the association between air pollution and health outcomes using air pollution data from a single monitoring site or average values from a few monitoring sites to represent the whole population of the study area. In fact, for a metropolitan city, ambient air pollution levels may differ significantly among the different areas. There is increasing concern that the relationships between air pollution and mortality may vary with geographical area (Chen, Mengersen et al. 2007). Additionally, some studies have indicated that socio-economic status can act as a confounder when investigating the relation between geographical location and health (Scoggins, Kjellstrom et al. 2004). This study examined the spatial variation in the relationship between long-term exposure to gaseous air pollutants (including nitrogen dioxide (NO2), ozone (O3) and sulphur dioxide (SO2)), and cardiorespiratory mortality in Brisbane, Australia, during the period 1996 - 2004.

Hierarchical models for 2D presence/absence data having ambiguous zeroes: With a biogeographical case study on dingo behaviour

This dissertation is primarily an applied statistical modelling investigation, motivated by a case study comprising real data and real questions. Theoretical questions on modelling and computation of normalization constants arose from pursuit of these data analytic questions. The essence of the thesis can be described as follows. Consider binary data observed on a two-dimensional lattice. A common problem with such data is the ambiguity of zeroes recorded. These may represent zero response given some threshold (presence) or that the threshold has not been triggered (absence). Suppose that the researcher wishes to estimate the effects of covariates on the binary responses, whilst taking into account underlying spatial variation, which is itself of some interest. This situation arises in many contexts and the dingo, cypress and toad case studies described in the motivation chapter are examples of this. Two main approaches to modelling and inference are investigated in this thesis. The first is frequentist and based on generalized linear models, with spatial variation modelled by using a block structure or by smoothing the residuals spatially. The EM algorithm can be used to obtain point estimates, coupled with bootstrapping or asymptotic MLE estimates for standard errors. The second approach is Bayesian and based on a three- or four-tier hierarchical model, comprising a logistic regression with covariates for the data layer, a binary Markov Random field (MRF) for the underlying spatial process, and suitable priors for parameters in these main models. The three-parameter autologistic model is a particular MRF of interest. Markov chain Monte Carlo (MCMC) methods comprising hybrid Metropolis/Gibbs samplers is suitable for computation in this situation. Model performance can be gauged by MCMC diagnostics. Model choice can be assessed by incorporating another tier in the modelling hierarchy. This requires evaluation of a normalization constant, a notoriously difficult problem. Difficulty with estimating the normalization constant for the MRF can be overcome by using a path integral approach, although this is a highly computationally intensive method. Different methods of estimating ratios of normalization constants (N Cs) are investigated, including importance sampling Monte Carlo (ISMC), dependent Monte Carlo based on MCMC simulations (MCMC), and reverse logistic regression (RLR). I develop an idea present though not fully developed in the literature, and propose the Integrated mean canonical statistic (IMCS) method for estimating log NC ratios for binary MRFs. The IMCS method falls within the framework of the newly identified path sampling methods of Gelman & Meng (1998) and outperforms ISMC, MCMC and RLR. It also does not rely on simplifying assumptions, such as ignoring spatio-temporal dependence in the process. A thorough investigation is made of the application of IMCS to the three-parameter Autologistic model. This work introduces background computations required for the full implementation of the four-tier model in Chapter 7. Two different extensions of the three-tier model to a four-tier version are investigated. The first extension incorporates temporal dependence in the underlying spatio-temporal process. The second extensions allows the successes and failures in the data layer to depend on time. The MCMC computational method is extended to incorporate the extra layer. A major contribution of the thesis is the development of a fully Bayesian approach to inference for these hierarchical models for the first time. Note: The author of this thesis has agreed to make it open access but invites people downloading the thesis to send her an email via the 'Contact Author' function.

Identification of area-level influences on regions of high cancer incidence in Queensland, Australia : a classification tree approach

Background: Strategies for cancer reduction and management are targeted at both individual and area levels. Area-level strategies require careful understanding of geographic differences in cancer incidence, in particular the association with factors such as socioeconomic status, ethnicity and accessibility. This study aimed to identify the complex interplay of area-level factors associated with high area-specific incidence of Australian priority cancers using a classification and regression tree (CART) approach. Methods: Area-specific smoothed standardised incidence ratios were estimated for priority-area cancers across 478 statistical local areas in Queensland, Australia (1998-2007, n=186,075). For those cancers with significant spatial variation, CART models were used to identify whether area-level accessibility, socioeconomic status and ethnicity were associated with high area-specific incidence. Results: The accessibility of a person’s residence had the most consistent association with the risk of cancer diagnosis across the specific cancers. Many cancers were likely to have high incidence in more urban areas, although male lung cancer and cervical cancer tended to have high incidence in more remote areas. The impact of socioeconomic status and ethnicity on these associations differed by type of cancer. Conclusions: These results highlight the complex interactions between accessibility, socioeconomic status and ethnicity in determining cancer incidence risk.