The Commission of the European Union, as well the United States Environmental Protection Agency, have set limit values for some pollutants in the ambient air that have been shown to have adverse effects on human and environmental health. It is therefore important to identify regions where the probability of exceeding those limits is high. We propose a two-step procedure for estimating the probability of exceeding the legal limits that combines smoothing in the time domain with spatial interpolation. For illustration, we show an application to particulate matter with diameter less than 10 microns (PM₁₀) in the North-Italian region Piemonte.