هشتمين كنگره علوم خاك ايران، ۱۲ - ۹ شهريور ۱۳۸۲، رشت
بصورت سخنرانی
تعيين شار سايش و انتشار ذرات خاك با استفاده از خصوصيات خاك و شبيه سازي سرعت وجهت باد
سيد فخرالدين افضلي
دانشجوي دكتري خاكشناسي دانشگاه صنعتي اصفهان
مقدمه
پديده فرسايش بادي يا به عبارت ديگر حساسيت اراضي نسبت به باد مانند بسياري از پديده هاي طبيعت مورد توجه دانشمندان جهان قرار گرفته، و تا كنون مخصوصا" در دهه هاي اخير سعي و تلاش زيادي در زمينه كمي كردن معيارهاي شناسائي اين پديده بعمل آمده است. توانائي پيش بيني دقيق فرسايش خاك براي بسياري منظورها، از جمله برنامه هاي حفاظتي، منابع طبيعي و كاهش آلودگي هوائي ناشي از طوفان ضروري است (6) .چپيل (4) با استفاده از تونل باد نشان داد كه فرسايش پذيري خاك كاملا" تحت تأثير توزيع اندازه خاكدانه هاي خشك قرار مي گيرد. هاگن (5) چگونگي فرسايش بادي و سايش خاكدانه را در تونلهاي باد مورد مطالعه قرار داده و اظهار داشت كه سايش وانتشاربه عنوان منابع حركت خاك عمل مي كنند. اسكيدمور و ليتون (9) تأثير پايداري خاكدانه ها در حالت خشك را برميزان فرسايش بادي در ده خاك ايالت كانزاس آمريكا مورد ارزيابي قرار داده، و مدلي را جهت پيش بيني پايداري خاكدانه هابا استفاده از خصوصيات خاك در شرايط فرسايش بادي ارائه دادند. هاگن(8) دريافت كه پايداري خاكدانه ها يكي از عوامل پيش بيني كننده سايش سطحي و فرسايش است. وي اظهار ميدارد كه فقط از طريق پيش بيني پايداري خاكدانه و سله ها مي توان به پيشرفتهاي خوبي در زمينه معادلات شار سائيدگي دست يافت. از طرف ديگر پيش بيني ساعتي سرعت و جهت باد، با استفاده از شبيه سازي و مدلهاي استوكاستيك ويبل (2) با تخمين بسيار عالي ( باد، بعنوان عامل اصلي فرساينده ) در تركيب با عامل خاك (فرسايش پذير) بعنوان پيش بيني فرسايش در منطقه بسيار مهم ميباشد. لذا در اين تحقيق شار سايش و شار انتشار ذرات خاك در خاكهاي گوناگون در منطقه يزد، كه به نوعي پيش بيني فرسايش بادي با استفاده از آمار باد تاريخي و شبيه سازي شده اين منطقه نيز مي باشد، محاسبه گرديد.
مواد و روشها
نمونه برداري از محدوده اراضي واقع در تپه هاي شني واقع در شمال شهرستان يزد و همچنين اراضي كشاورزي در نقاطي با كيفيت مناسب خاك انجام گرفت(2). خصوصيات فيزيكوشيميائي اين خاكها نظير توزيع اندازه خاكدانه ها ، درصد ماده آلي، كربنات كلسيم معادل، بافت خاك، حداكثر ارتفاع پستي و بلندي خاك يا ( Zmax ) در هر كدام از اراضي و ميزان شكستگي اوليه خاكدانه در هريك از خاكها اندازه گيري شد(2). چنانچه شرايط فرسايش بادي فراهم باشد و سرعت باد نيز از 8 متر بر ثانيه بيشتر شود، ذرات خاك شروع به حركت مي كنند (7). در اين مرحله با توجه به نتايج بدست آمده از پيش بيني سرعت و جهت باد با استفاده از مدل ويبل و شبيه سازي استوكاستيك( كه با استفاده از آمار سي ساله ايستگاه سينوپتيك يزد انجام داديم ) و با در نظر گرفتن سرعتهاي بيشتر از 8 متر بر ثانيه( باد هاي فرساينده ) (1) ميزان شار انتشار(Gen) و شار سايش (Gan) محاسبه و مجموع آنها تعيين گرديد(1).
نتايج و بحث
ضرايب همبستگي بين شارها ي محاسبه شده با استفاده از آمار داده هاي شبيه سازي شده سالهاي 1980 - 1966 و همچنين با داده هاي تاريخي سالهاي 1981 تا 1995 در ماههاي پر باد و در كل سال، نشان دهنده همبستگي خوبي بين آنها ميباشد. خـاك Sandy Loam، مساعدترين خاك از نظرفرسايش پذيري در ميان خاكهاي مورد آزمايش در اين تحقيق بود. از نظر مقايسه ميان خواص مختلف خاكهاي مورد آزمايش، تغييرات بافت خاك و از همه مهمتر درصد رس در اين زمينه نقش تعيين كننده دارد. بنابراين مي بينيم خاك Sandy Loamاز جهات دو خصوصيت ذكر شده، وضعيت نا مناسبي را نسبت به خاكهاي ديگر دارد. بيشتر بودن مقدار Gen نسبت به Gan در خاك Loamy Sand به اين دليل است كه درصد ذرات آزاد كوچكتر از1/0 ميليمتر اين خاك بيشتر و ميزان خاكدانه كمتر است( تپه هاي شني) . اين محاسبات و نظاير آن در ماههاي پربادتر نتايج بهتري را نسبت به ماههاي كم باد، براي داده هاي شبيه سازي شده و تاريخي نشان مي دهد. با توجه به ضرايب همبستگي تعيين شده و با وجود اينكه ماه آگوست، يكي از ماههاي كم باد به حساب مي آيد، همبستگي نشان داده شده بين اين داده ها متوسط است، ولي در ماههاي پر باد و در كل سال ضرايب همبستگي، نشان دهنده همبستگي خوبي بين داده هاي ذكر شده هستند و نشان دهنده اعتبار مدل شبيه سازي باد و محاسبه فرسايش ميباشند. بين سرعت و جهت باد تاريخي و شييه سازي شده اختلاف كوچكي وجود دارد ولي داده هاي ذكر شده، تطابق خوبي را نشان مي دهند. در ماههاي كم بادي همچون آگوست، تطابق داده هاي سرعت و جهت باد تاريخي و شبيه سازي شده نسبت به ماههاي پر باد كمتر است كه با توجه به اينكه فرسايش بادي با توان سوم سرعت باد متناسب است، بنابراين تفاوت هاي كوچك در سرعت هاي باد، در محاسبه شار فرسايش بادي، به اختلاف بزرگي تبديل شده و موجب مي شوند همبستگي نشان داده شده بين اين داده ها، در ماههاي كم بادي همچون آگوست، نسبت به ماههاي پر باد كمتر باشد. بنابراين با توجه به نتايج حاصل از تركيب اين دو مدل، ميتوان به نتايج ارزشمندي در پيش بيني فرسايش بادي منطقه دست يافت.
منابع مورد استفاده
1- افضلي، س. ف. 1380. مطالعه بادهاي فرساينده در منطقه يزد. پايان نامه كارشناسي ارشد بخش خاكشناسي . دانشكده كشاورزي. دانشگاه شيراز.
2- افضلي، س. ف. و ن. ميرزامصطفي. 1381. مدل احتمال تجربي ( استوكاستيك ) ويبل بعنوان ابزار لازم جهت پيش بيني پتانسيل فرسايندگي باد ( مطالعه موردي: منطقه يزد ). مجموعه مقالات همايش ملي مديريت اراضي – فرسايش خاك و توسعه پايدار. اراك . صفحه 516.
3-حكيم زاده، م. ع. 1379. تخمين شدت فرسايش بادي در خاكهاي مختلف دشت يزد - اردكان با استفاده از نحوه توزيع اندازه و پايداري خاكدانه ها.گزارش مسئله مخصوص دوره دانشوري دانشگاه شيراز.
4- Chepil, W. S. 1950. Properties of soil which influence wind erosion: II. Dry aggregate structure as an index of erodibility. Soil Sci. 69: 403-414.
5- Hagen, L. J. 1991. Wind erosion mechanics: Abrasion of aggregated soil. Trans. ASAE. 34:831-837.
6- Hagen, L. J. 1995. Preface. P. X. In Proc. WEPP/WEPS Symp. Des Moines, IA. 10-11 Aug. 1995. Soil Water Conserv. Soc. Ankeny, IA. , U. S. A.
7- Hagen, L. J., L. E. Wagner, and J. Tatarko. 1995. Wind erosion prediction system (WEPS). P. I1-I8. In Proc. WEPP/WEPS Symp. Des Moines, IA. 10-11 Aug. 1995. Soil Water Conserv. Soc. Ankeny, IA. , U. S. A.
http://www.fafzali.com/farsi/publications/ps2-fa.htm : منبع مقاله