Spatiotemporal joint analysis of PM and Ozone in California with INLA approach.

The substantial threat of concurrent air pollutants to public health is increasingly severe under climate change. To identify the common drivers and extent of spatiotemporal similarity of PM and ozone (O), this paper proposed a log Gaussian-Gumbel Bayesian hierarchical model allowing for sharing a stochastic partial differential equation and autoregressive model of order one (SPDE-AR(1)) spatiotemporal interaction structure. The proposed model, implemented by the approach of integrated nested Laplace approximation (INLA), outperforms in terms of estimation accuracy and prediction capacity for its increased parsimony and reduced uncertainty, especially for the shared O sub-model. Besides the consistently significant influence of temperature (positive), extreme drought (positive), fire burnt area (positive), gross domestic product (GDP) per capita (positive), and wind speed (negative) on both PM and O, surface pressure and precipitation demonstrate positive associations with PM and O, respectively. While population density relates to neither. In addition, our results demonstrate similar spatiotemporal interactions between PM and O, indicating that the spatial and temporal variations of these pollutants show relatively considerable consistency in California. Finally, with the aid of the excursion function, we see that the areas around the intersection of San Luis Obispo and Santa Barbara counties are likely to exceed the unhealthy O level for USG simultaneously with other areas throughout the year. Our findings provide new insights for regional and seasonal strategies in the co-control of PM and O. Our methodology is expected to be utilized when interest lies in multiple interrelated processes in the fields of environment and epidemiology.