Meteorological associations of hourly PM₂.₅ in Lhasa, Tibetan Plateau: antecedent weather, trend specification and out-of-sample evaluation
Abstract
Statistical models that relate PM₂.₅ to meteorology are widely used to separate weather from emission effects, but they are rarely evaluated out of sample in clean, high-altitude cities, and the role of their trend term is seldom examined. This study fitted penalized additive spline models to hourly city-average PM₂.₅ in Lhasa, Tibetan Plateau (about 3,650 m a.s.l.), for 2019–2024 (51,955 valid hours), using ERA5 temperature, relative humidity, wind speed and precipitation as current-hour and 24-hour antecedent predictors with seasonal, diurnal and long-term controls, and benchmarked them against calendar, linear and random forest models. Three results emerge. First, the point-estimated wind and humidity responses are carried mainly by antecedent conditions: in perturbation decompositions, the 24-hour terms carry most of the fitted wind and humidity sensitivity. Second, the trend specification governs both the weather-normalized series and out-of-sample bias: a linear trend forces a near-monotonic normalized decline (10.98 to 9.66 µg/m³), although observed annual means fell from 12.16 to 8.55 µg/m³ by 2022 and rebounded to 10.36 µg/m³ in 2024. Third, a model frozen after 2019–2024 explained 42% of hourly and 60% of daily variance in 2025 (RMSE 5.35 and 2.80 µg/m³; calendar-only benchmark: 32% and 41%) and modestly outperformed the linear benchmark, but errors were concentrated in winter and in a few high-concentration hours; Gutor-evening residuals exceeded the same-year, same-hour winter background in five of the six years of the 2019–2024 record. Weather-normalized trends for Plateau cities should therefore be reported with alternative trend specifications and event-level diagnostics.
Related
- Early warning of road icing from antecedent meteorological conditions without consistent critical slowing down
- Long-term evolution and climatic associations of dry and warm snow droughts across the Tibetan plateau
- Evaluation of turbulent flux parameterizations over a continental glacier on the Tibetan Plateau
- Characteristics and simulation of ozone pollution in Haikou City under combined effects of meteorological factors and pollutants
- A support vector machine classifier of PM₂.₅-based air quality index categories from co-pollutants and meteorological variables in Ahvaz, Iran
BibTeX
@misc{ma2026meteorological,
title = {{Meteorological associations of hourly PM₂.₅ in Lhasa, Tibetan Plateau: antecedent weather, trend specification and out-of-sample evaluation}},
author = {Zhongjie Ma},
howpublished = {EarthArXiv},
year = {2026},
month = nov,
doi = {10.31223/x5z817},
url = {https://eartharxiv.org/repository/view/15941/},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Meteorological associations of hourly PM₂.₅ in Lhasa, Tibetan Plateau: antecedent weather, trend specification and out-of-sample evaluation
AU - Zhongjie Ma
PY - 2026
DA - 2026/11/03
PB - EarthArXiv
DO - 10.31223/x5z817
UR - https://eartharxiv.org/repository/view/15941/
AB - Statistical models that relate PM₂.₅ to meteorology are widely used to separate weather from emission effects, but they are rarely evaluated out of sample in clean, high-altitude cities, and the role of their trend term is seldom examined. This study fitted penalized additive spline models to hourly city-average PM₂.₅ in Lhasa, Tibetan Plateau (about 3,650 m a.s.l.), for 2019–2024 (51,955 valid hours), using ERA5 temperature, relative humidity, wind speed and precipitation as current-hour and 24-hour antecedent predictors with seasonal, diurnal and long-term controls, and benchmarked them against calendar, linear and random forest models. Three results emerge. First, the point-estimated wind and humidity responses are carried mainly by antecedent conditions: in perturbation decompositions, the 24-hour terms carry most of the fitted wind and humidity sensitivity. Second, the trend specification governs both the weather-normalized series and out-of-sample bias: a linear trend forces a near-monotonic normalized decline (10.98 to 9.66 µg/m³), although observed annual means fell from 12.16 to 8.55 µg/m³ by 2022 and rebounded to 10.36 µg/m³ in 2024. Third, a model frozen after 2019–2024 explained 42% of hourly and 60% of daily variance in 2025 (RMSE 5.35 and 2.80 µg/m³; calendar-only benchmark: 32% and 41%) and modestly outperformed the linear benchmark, but errors were concentrated in winter and in a few high-concentration hours; Gutor-evening residuals exceeded the same-year, same-hour winter background in five of the six years of the 2019–2024 record. Weather-normalized trends for Plateau cities should therefore be reported with alternative trend specifications and event-level diagnostics.
ER -