A spatiotemporal model for Mexico City ozone levels (original) (raw)

Spatio-temporal stochastic modeling of tropospheric ozone concentration in Mexico City

Tropospheric ozone is one of the most damaging atmospheric pollutants for human health. In Mexico City this pollutant systematically exceeds the health official norm. In this city there is a network of atmospheric monitoring stations, which is responsible for watching the air quality. The network reports pollutants concentration with hourly frequency. For knowing the potential risk population is exposed to, the stations records are not enough, but considering pollution in other locations, e.g., in a regular grid, is also necessary, this is why it should be esti- mated. In this contribution it is proposed an ozone pollution model that considers the atmospheric physical and chemical dynamics, expressed by the numerical mete- orology estimation of WRF-Chem, as a component of a spatio-temporal stochastic process. A stationary and separable structure is assumed for the spatio-temporal covariance function, this allows applying the kriging methodology for estimating the ozone concentration in non-observed locations. The efficacy of this approach is evaluated with 6783 estimations in different points. The results confirm that in general combining kriging with WRF-Chem is an improvement with regard to the estimations of the last one alone. A spatial pattern about the efficacy of kriging is also identified, and it is showed to be related to the geographic density of the monitoring stations.

Explicit Modeling of Meteorological Explanatory Variables in Short-Term Forecasting of Maximum Ozone Concentrations via a Multiple Regression Time Series Framework

Atmosphere

Statistical time series forecasting is a useful tool for predicting air pollutant concentrations in urban areas, especially in emerging economies, where the capacity to implement comprehensive air quality models is limited. In this study, a general multiple regression with seasonal autoregressive moving average errors model was estimated and implemented to forecast maximum ozone concentrations with a short time resolution: overnight, morning, afternoon and evening. In contrast to a number of short-term air quality time series forecasting applications, the model was designed to explicitly include the effects of meteorological variables on the ozone level as exogenous variables. As the application location, the model was constructed with data from five monitoring stations in the Monterrey Metropolitan Area of Mexico. The results show that, together with structural stochastic components, meteorological parameters have a significant contribution for obtaining reliable forecasts. The res...

Determination of Spatial-Temporal Correlation Structure of Troposphere Ozone Data in Tehran City

Spatial-temporal modeling of air pollutants, ground-level ozone concentrations in particular, has attracted recent attention because by using spatial-temporal modeling, can analyze, interpolate or predict ozone levels at any location. In this paper we consider daily averages of troposphere ozone over Tehran city. For eliminating the trend of data, a dynamic linear model is used, then some features of correlation structure of de-trended data, such as stationarity, symmetry and separability are considered. Next based on the obtained features, an appropriate model is proposed. This model can be used for future predictions of ozone in Tehran.

A class of covariate-dependent spatiotemporal covariance functions for the analysis of daily ozone concentration

Annals of Applied Statistics, 2011

In geostatistics, it is common to model spatially distributed phenomena through an underlying stationary and isotropic spatial process. However, these assumptions are often untenable in practice because of the influence of local effects in the correlation structure. Therefore, it has been of prolonged interest in the literature to provide flexible and effective ways to model nonstationarity in the spatial effects. Arguably, due to the local nature of the problem, we might envision that the correlation structure would be highly dependent on local characteristics of the domain of study, namely, the latitude, longitude and altitude of the observation sites, as well as other locally defined covariate information. In this work, we provide a flexible and computationally feasible way for allowing the correlation structure of the underlying processes to depend on local covariate information. We discuss the properties of the induced covariance functions and methods to assess its dependence on local covariate information. The proposed method is used to analyze daily ozone in the southeast United States.

Space. time interpolation of daily air temperatures

We propose a model to describe the mean function as well as the spatio-temporal covariance structure of 15 years of both maximum and minimum daily temperature data from 190 stations throughout the region of Catalonia (Spain), with daily data covering the period 1994-2008. Our aim is threefold: (a) estimation of the long-term trend of maximum and minimum temperatures; (b) assessing the spatial and temporal variability of temperatures, and (c) interpolation of the spatial temperatures at any given time. Long-term trend, annual harmonics and winds were considered as explanatory variables of the mean function. The parameters associated with these variables were allowed to vary between stations and within each year. We controlled temporal autocorrelation by means of ARMA models. For the spatial covariance structure we used the Matérn family of covariance functions and a nugget term. Spatio-temporal models were built as Bayesian hierarchical models with two stages following the integrated nested place Laplace approximation (INLA) for Bayesian inference. For the final model estimation we used a two-stage approach, in which we first assumed the stations were spatially independent, and then we modeled the spatio-temporal covariance using the interim posterior from the residuals of the model in the first-stage as prior distributions of replications of a spatial process. We allowed all spatial parameters to also vary with time.

Forecasting Environmental Data: An example to ground-level ozone concentration surfaces

2022

Environmental problems are receiving increasing attention in socio-economic and health studies. This in turn fosters advances in recording and data collection of many related real-life processes. Available tools for data processing are often found too restrictive as they do not account for the rich nature of such data sets. In this paper, we propose a new statistical perspective on forecasting spatial environmental data collected sequentially over time. We treat this data set as a surface (functional) time series with a possibly complicated geographical domain. By employing novel techniques from functional data analysis we develop a new forecasting methodology. Our approach consists of two steps. In the first step, time series of surfaces are reconstructed from measurements sampled over some spatial domain using a finite element spline smoother. In the second step, we adapt the dynamic functional factor model to forecast a surface time series. The advantage of this approach is that ...

Combining measurements and physical model outputs for the spatial prediction of hourly ozone space–time fields

This technical report extends to a spatial setting, an existing temporal two-step linear regression recalibration procedure designed to make the outputs from a deterministic time series simulator comparable with measurements of that series. The result, which Kriges the site specific coefficients of that procedure, enables both temporal forecasting and spatial prediction. Although the extension is somewhat ad hoc, unlike an alternative in work now in preparation, it is computationally simple and fairly transparent. Moreover, in an application where the procedure is used to combing measurements of and simulated model outputs for an hourly ozone field over the eastern and central regions of the United States, we find that it outperforms the model output alone and Kriging, a purely spatial method.

Improved space–time forecasting of next day ozone concentrations in the eastern US

Atmospheric Environment, 2009

There is an urgent need to provide accurate air quality information and forecasts to the general public and environmental health decision-makers. This paper develops a hierarchical space-time model for daily 8-h maximum ozone concentration (O 3 ) data covering much of the eastern United States. The model combines observed data and forecast output from a computer simulation model known as the Eta Community Multi-scale Air Quality (CMAQ) forecast model in a very flexible, yet computationally fast way, so that the next day forecasts can be computed in real-time operational mode. The model adjusts for spatio-temporal biases in the Eta CMAQ forecasts and avoids a change of support problem often encountered in data fusion settings where real data have been observed at point level monitoring sites, but the forecasts from the computer model are provided at grid cell levels. The model is validated with a large amount of set-aside data and is shown to provide much improved forecasts of daily O 3 concentrations in the eastern United States.

Ozone forecasting from meteorological variables

Chemometrics and Intelligent Laboratory Systems, 1998

A multivariate modeling approach is presented for ozone forecasting from meteorological variables. For each prediction, a separate optimized multivariate regression model is constructed. The optimization involves the determination of the size of the training set by an internal validation procedure. A grid search is used in order to determine how many observations the training set should contain. A straightforward ordinary least squares procedure leads to systematic positive or negative devia-Ž . tions between measured and predicted ozone concentrations. Partial least squares regression PLS in combination with training-set selection and variable selection, gave an overall correlation coefficient of 0.83 between observed and measured ozone levels. Appropriate weighting of the observations in the training sets improved the result to give an overall correlation coefficient between measured and predicted ozone levels of 0.86. The dependence of the optimal size of the training set on the number and location of missing data in the data matrix was also investigated. q 1998 Elsevier Science B.V. All rights reserved.

Spatial interpolation of urban air temperatures using satellite-derived predictors

Theoretical and Applied Climatology, 2020

Air temperatures in urban environments are usually obtained from sparse weather stations that provide limited information with regard to spatial patterns. Effective methods that predict air temperatures (T air) in urban areas are based on statistical models which utilize remotely sensed and geographic data. This work aims to compute T air predictions for diurnal and nocturnal time intervals using predictive models that do not exploit information on Land Surface Temperatures. The models are developed based on explanatory variables that describe the urban morphology, land cover and terrain, aggregated at 100 m × 100 m resolution, combined with in situ T air measurements from urban meteorological stations. The case study is the urban and per-urban area of Heraklion, Greece, where a dense meteorological station network is available since 2016. Moran's eigenvector filtering and an autoregressive moving average residual specification are implemented to account for spatial and temporal correlations. The statistical models display satisfactory predictive performance, with mean annual Mean Absolute Error (MAE) equal to 0.36°C, 0.34°C, 0.42°C and 0.54°C,