Multi-elemental isotopic approach associated with a land-use characteristic sampling strategy may be relevant for conducting biomonitoring studies to determine the spatial extent of atmospheric contamination sources. In this work, we investigated how the combined isotopic signatures in epiphytic lichens of two major metallic pollutants, lead (206Pb/207Pb) and mercury (¿202Hg, ¿199Hg), together with the isotopic composition of nitrogen and carbon (¿15N, ¿13C), can be used to better constrain atmospheric contamination inputs. To this end, an intensive and integrated sampling strategy based on land-use characteristics (Geographic information system, GIS) over a meso-scale area (Pyrénées-Atlantiques, SW France) was applied to more than 90 sampling stations. To depict potential relationships between such multi-elemental isotopic fingerprint and land-use characteristics, multivariate analysis was carried out. Combined Pb and Hg isotopic signatures resolved spatially the contribution of background atmospheric inputs from long range transport, from local legacy contamination (i.e. Pb) or actual industrial inputs (i.e. Pb and Hg from steel industry). Application of clustering multivariate analysis to all studied isotopes provided a new assessment of the region in accordance with the land-use characteristics and anthropogenic pressures.