Observational records are more often than not influenced by residual non‐climatic factors which must be detected and adjusted for prior to their usage. In this work, we present a novel approach, named Radiosounding HARMonization (RHARM), providing a homogenized data set of temperature, humidity and wind profiles along with an estimation of the measurement uncertainties for 697 radiosounding stations globally. The RHARM method has been used to adjust twice daily (0000 and 1200 UTC) radiosonde data holdings at 16 pressure levels in the range 1,000–10 hPa, from 1978 to present, provided by the Integrated Global Radiosonde Archive. Relative humidity (RH) data are limited to 250 hPa. The applied adjustments are interpolated to all reported levels. RHARM is the first data set to provide homogenized time series with an estimation of the observational uncertainty at each sounding pressure level. By construction, RHARM adjusted fields are not affected by cross‐contamination of biases across stations and are fully independent of reanalysis data. Analysis of trends for temperature, RH and winds highlights increased geographical coherency of trends over 1978–2000 globally, but especially in the Northern Hemisphere and South America. RHARM shows warming trends of 0.39 K/decade at 300 hPa in the Northern Hemisphere and of 0.25 K/decade in the tropics. The RHARM adjustments also reduce differences with the European Centre for Medium‐Range Weather Forecast ERA5 reanalysis, with the strongest effect in the Northern Hemisphere for temperature and relative humidity. For wind speed, the comparison indicates a good agreement with ERA5 in the troposphere.

The air in the Lombardy region, Italy, is one of the most polluted in Europe because of limited air circulation and high emission levels. There is a large scientific consensus that the agricultural sector has a significant impact on air quality. To support studies quantifying the role of the agricultural and livestock sectors on the Lombardy air quality, this paper presents a harmonised dataset containing daily values of air quality, weather, emissions, livestock, and land and soil use in the years 2016–2021, for the Lombardy region. The daily scale is obtained by averaging hourly data and interpolating other variables. In fact, the pollutant data come from the European Environmental Agency and the Lombardy Regional Environment Protection Agency, weather and emissions data from the European Copernicus programme, livestock data from the Italian zootechnical registry, and land and soil use data from the CORINE Land Cover project. The resulting dataset is designed to be used as is by those using air quality data for research.

This study presents a comparative analysis of three predictive models with an increasing degree of flexibility: hidden dynamic geostatistical models (HDGM), generalised additive mixed models (GAMM), and the random forest spatiotemporal kriging models (RFSTK). These models are evaluated for their effectiveness in predicting $$\text {PM}_{2.5}$$PM 2.5 concentrations in Lombardy (North Italy) from 2016 to 2020. Despite differing methodologies, all models demonstrate proficient capture of spatiotemporal patterns within air pollution data with similar out-of-sample performance. Furthermore, the study delves into station-specific analyses, revealing variable model performance contingent on localised conditions. Model interpretation, facilitated by parametric coefficient analysis and partial dependence plots, unveils consistent associations between predictor variables and $$\text {PM}_{2.5}$$PM 2.5 concentrations. Despite nuanced variations in modelling spatiotemporal correlations, all models effectively accounted for the underlying dependence. In summary, this study underscores the efficacy of conventional techniques in modelling correlated spatiotemporal data, concurrently highlighting the complementary potential of Machine Learning and classical statistical approaches.