آگاهی از توزیع مکانی غلظت فلزات سنگین جهت پایش آلودگی خاک و حفظ کیفیت محیط زیست ضروری است. این تحقیق با هدف تهیه نقشه توزیع مکانی آلودگی غلظت فلزات سنگین منگنز، مس، روی و آهن در زمین های کشاورزی شهرستان هریس واقع در استان آذربایجان شرقی انجام شده است. بدین منظور با استفاده از روش نمونه برداری سیستماتیک تصادفی 370 نمونه خاک سطحی در عمق صفر تا 30 سانتیمتری جمع آوری گردید و غلظت کل این فلزات، در آن ها اندازه گیری شد. به منظور مدل سازی تغییرات مکانی غلظت فلزات سنگین در خاک های منطقه، از روش های کریجینگ معمولی و توابع پایه شعاعی در محیط ArcGIS بهره گرفته شد. در بین مدل های مختلف کریجینگ معمولی و توابع پایه شعاعی، پایین ترین مقادیر RMSE و MAE و نزدیکی بسیار زیاد آماره ی MBE به صفر و بالاترین ضریب همبستگی R به منزله ملاک ارزیابی بهترین روش برای مدل سازی الگوی توزیع و پراکنش فلزات سنگین در نظر گرفته شد. هر دو روش از دقت خوبی برخوردار بودند ولی بر اساس ملاک ارزیابی، برای فلزات منیزیم و مس، روش کریجینگ با مدل نمایی، برای فلز روی، روش کریجینگ با مدل کروی و برای فلز آهن، روش RBF با تابع مالتی کوادریک انتخاب شد. نقشه ی پهنه بندی بعد از انتخاب مدل برتر برای هر یک از عناصر، ترسیم و با درنظر گرفتن استاندارد مؤسسه حفاظت خاک و آب ، کلاسه بندی شد و عناصر در چهار کلاس کمبود، کفایت، زیاد و آلودگی قرار گرفتند. تجزیه و تحلیل نقشه توزیع مکانی آلودگی فلزات سنگین مورد مطالعه نشان داد که برای عنصر منگنز بیشتر اراضی در محدوده ی کفایت قرار داشتند و برای عنصر مس مشخص شد که حدود 92 درصد در کلاس زیاد قرار داشته و حدود 8 درصد از محدوده ی مورد مطالعه دچار آلودگی بود. برای عناصر روی و آهن نه تنها آلودگی مشاهده نشد بلکه به ترتیب حدود 96 و 80 درصد دچار کمبود بود .
