Regional variability and sources of error in mortality risk of total and bushfire-specific PM2.5 (23946)
Background
Bushfire-specific PM2.5 exposure has been shown to increase the risk of all-cause mortality in Australia. However, potential pitfalls exist when quantifying this relationship that may introduce bias and error. For example,
- Modelled PM5 is more valid in areas where observed data (i.e. government monitoring sites) is available, often in more populated urban centres
- Aggregation using meta-regression may bias national estimates towards the most populous areas
- Areas are subject to seasonal variability, and differences in sources of PM5, which may affect the validity of separating bushfire-specific exposure from total PM2.5
- Models are often sensitive to parameterisation.
The current study aimed to explore the impact of such differences on national, and sub-national, estimates of risk.
Methods
Newly developed high resolution (5km² grid) 20-year Australia-wide daily PM2.5 surfaces were decomposed into daily seasonal, trend and remainder components. Bushfire days were identified using a combination of environmental exposure flags. Regression models utilising an underlying quasi-Poisson distribution were fitted. Various sensitivity analyses were conducted, varying the parameterisation of models, the level of smoothing applied to the exposure surfaces, as well as restricting analysis to areas proximal to monitoring stations.
Results
National estimates showed a positive, exposure-dependent risk relationship between PM2.5 and all-cause mortality. Greater Sydney and Melbourne, the two most populous cities, often reflected national estimates. However, other sub-national estimates were considerably more variable and showed greater sensitivity to model specification and rate of smoothing. Areas with greater temperature fluctuation, differential seasonal patterns and lower density of monitoring stations, showed especially high sensitivity.
Discussion
Sub-national results showed sensitivity to exposure definition and model specification. This persisted after adjusting for population size and estimates of socioeconomic status by area. Caution is advised when applying national PM2.5 and mortality risk estimates to sub-national regions.