Urban–Rural PM2.5 Dynamics in Kraków, Poland: Patterns and Source Attribution

Hourly PM2.5 concentrations were measured from February to May 2025 by a network of low-cost sensors located in urban Kraków and its surrounding municipalities. Temporal variability associated with the transition from the heating period to the spring months, together with spatial contrasts, were assessed with principal component analysis (PCA), urban–rural difference curves, and a detailed examination of the most severe smog episode (12–13 February). Particle trajectories generated with the HYSPLIT dispersion model, run in a coarse-grained, 36-task parallel configuration, were combined with kernel density mapping to trace emission pathways. The results show that peak concentrations coincide with the heating season; rural sites recorded higher amplitudes and led the urban signal by up to several hours, implicating external sources. Time-series patterns, PCA loadings, and HYSPLIT density fields provided mutually consistent evidence of pollutant advection toward the city. Parallelizing HYSPLIT on nine central processing unit (CPU) cores reduced the runtime from more than 600 s to about 100 s (speed-up ≈ 6.5), demonstrating that routine episode-scale analyses are feasible even on modest hardware. The findings underline the need to extend monitoring and mitigation beyond Kraków’s administrative boundary and confirm that coarse-grained parallel HYSPLIT modeling, combined with low-cost sensor data and relatively basic statistics, offers a practical framework for rapid source attribution.