Abstract: Based on continuous observation data of near-surface ozone from 19 ambient air quality monitoring stations in Zibo City, Shandong Province from 2016 to 2019, the spatial-temporal characteristics of atmospheric O₃ were studied using Pearson correlation analysis and inverse distance weighting method. The results showed that: (1) From 2016 to 2019, the 90^th percentile of ρ(O_{3-8 h}) (maximum 8 h concentration on O₃ day) in Zibo City was 184-203 μg/m³, which is 1.15-1.27 times of the limit value of level Ⅱ (160 μg/m³) in National Ambient Air Quality Standards (GB 3095-2012). (2) The seasonal average of ρ(O_{3-8 h}) was characterized by summer (155 μg/m³) > spring (129 μg/m³) > autumn (104 μg/m³) > winter (60 μg/m³). The monthly variation of ρ(O_{3-8 h}) showed 'double peak' pattern, peaking in June (186 μg/m³) and September (147 μg/m³), respectively. The diurnal variations of O₃ showed a 'single peak' pattern, with the peak value of 117 μg/m³ at around 14:00. (3) The spatial distribution of ρ(O_{3-8 h}) showed 'low in south-north area and high in middle area' pattern, with a uniform trend in the regional. (4) The variation of precursors such as CO, NO, NO₂ and NO_x were negatively correlated with ρ(O_{3-8 h}), but this correlation was weak in summer. Local pollutants at night in winter contributed more to O_x than in daytime. (5) ρ(O_{3-8 h})was significantly positively correlated with temperature and wind speed, and negatively correlated with humidity. High ρ(O_{3-8 h}) was observed under southwesterly dominated wind. The results showed typical photochemical characteristics of O₃ in Zibo, and O₃ pollution gradually increased from 2016 to 2019.