Urban forests provide important ecosystem services, such as urban air quality improvement by removing pollutants. While robust evidence exists that plant physiology, abundance, and distribution within cities are basic parameters affecting the magnitude and efficiency of air pollution removal, little is known about effects of plant diversity on the stability of this ecosystem service. Here, by means of a spatial analysis integrating system dynamic modeling and geostatistics, we assessed the effects of tree diversity on the removal of tropospheric ozone (O3) in Rome, Italy, in two years (2003 and 2004) that were very different for climatic conditions and ozone levels. Different tree functional groups showed complementary uptake patterns, related to tree physiology and phenology, maintaining a stable community function across different climatic conditions. Our results, although depending on the city-specific conditions of the studied area, suggest a higher function stability at increasing diversity levels in urban ecosystems. In Rome, such ecosystem services, based on published unitary costs of externalities and of mortality associated with O3, can be prudently valued to roughly US$2 and $3 million/year, respectively.