Abstract:
هدایت هیدرولیکی اشباع خاک (Ks) در اکثر مدل های شبیه سازی فرسایش و فرایندهای هیدرولوژیکی خاک در آبخیزها، نقش مهمی را ایفا می کند. اندازه گیری مستقیم هدایت هیدرولیکی اشباع (Ks) خاک، کاری وقت گیر، دشوار و پرهزینه است. هدف از این پژوهش، مقایسه ی دقت توابع انتقالی (PTFs) رگرسیونی، شبکه عصبی مصنوعی (ANN) و نروفازی در برآورد Ks خاک های متاثر از نمک واقع در شمال غرب دریاچه ی ارومیه بود. برای تعیین برخی متغیرهای فیزیکی و شیمیایی زودیافت خاک، تعداد 100 نمونه خاک دست خورده و دست نخورده از عمق 0 تا 10 سانتیمتری اراضی کشاورزی و بایر بخش شندآباد در منطقهی شبستر برداشته شد. متغیر Ks، در آزمایشگاه به روش بار ثابت یا افتان اندازه گیری شد. برای اشتقاق توابع رگرسیونی از نرمافزار SPSS استفاده شد و برای توابع ANN و نروفازی از نرم افزار MATLAB. هشتاد درصد داده ها برای آموزش و بیست درصد آن برای آزمون توابع به کار رفت. نتایج توابع رگرسیونی، ANN و نروفازی نشان داد که تابع انتقالی با دو متغیر سیلت و جرم مخصوص ظاهری، بهترین تابع برای برآورد Ks خاک در منطقه ی مورد مطالعه است. مقادیر ضریب تبیین (R2)، مجذور میانگین مربعات خطا (RMSE) و میانگین خطا (ME) بهترتیب 0/65، cm/min 0/119 وcm/min 059/0- و 0/73، 00/87 cm/min وcm/min 0/006 و 0/69،cm/min 0/127 و cm/min 0/051- به ترتیب برای بهترین تابع رگرسیونی، ANN و نروفازی به دست آمد. بنابراین، توابع ANN به دلیل داشتن R2 بالا و RMSE پایین در مقایسه با توابع رگرسیونی و نروفازی، دقت بیشتری برای برآورد Ks خاک در منطقه ی مورد مطالعه دارد.
Extended abstract
1- Introduction
Soil saturated hydraulic conductivity (Ks) is an important factor in the estimation of water, solute transport models and erosion processes. Direct measurement of soil saturated hydraulic conductivity (Ks) in field and laboratory is time-consuming, laborious and expensive because of high temporal and spatial variability; especially in salt-affected soils around Urmia Lake, Ks measurement is difficult because of high sodium concentration and consequently poor stability of soil aggregates. Therefore, many different regressions, artificial neural network (ANN) and the neuro-fuzzy pedotransfer functions (PTF) have been developed to estimate Ks from readily available soil variables such as sand, silt, clay, bulk density (BD), particle density (PD), electrical conductivity (EC), pH, and organic carbon (OC). The objectives of this study were to derive pedotransfer functions by using regression, artificial neural network, and neuro-fuzzy methods to estimate Ks from some soil variables in the salt-affected soils selected from the northwest of Urmia Lake and to compare the performance of the neuro-fuzzy, artificial neural network and regression models.
2- Methodology
Disturbed and undisturbed (steel cylinders with 5 cm diameter and height) soil samples (n= 100) were systematically taken from 0-10 cm soil depth of bare and agricultural lands of Shend Abad region located at the 15 km of Shabestar city, northwest of Urmia Lake, Iran (45° 36ʹ 34ʺ E and 38° 6ʹ 37ʺ N). The values of sand, silt, and clay (hydrometer method), CaCO3 (titration method), bulk density (cylinder method), particle density (pycnometer method), organic carbon (wet oxidation method), and total porosity (calculating from BD and PD) were measured in the laboratory. The mean geometric diameter (dg) of soil particles was computed using the percentages of sand, silt, and clay. The EC and sodium adsorption ratio (SAR) were measured in 1:2.5 (soil: distilled water) extra. The pHe was determined in a/the saturated paste. The soil saturated hydraulic conductivity (Ks) was measured by constant (agricultural lands) and falling (bare lands) head method using steel cylinders in the laboratory. The data were divided into two series as 80 data for training and 20 data for testing. The SPSS 18 software with a/the stepwise method to derive the regression PTFs and MATLAB software to derive the artificial neural network and neuro-fuzzy PTFs were used. A three-layer perceptron network and the tangent sigmoid transfer function were used for the artificial neural network modeling. In estimating soil saturated hydraulic conductivity, the accuracy of neuro-fuzzy, artificial neural network and regression pedotransfer functions were evaluated by the coefficient of determination (R2), root mean square error (RMSE) and mean error (ME) criteria.
3- Results & Discussion
Most of the studied soil variables had good distribution for developing and evaluating regression, ANN, and neuro-fuzzy PTFs. The high values of the coefficient of variation (CV) were found for SAR (167.86%), Ks (130.36%), EC (117.05%), dg (88.44%), clay (73.23%), OC (58.46%) and sand (51.47%) in the studied area. The textural classes of studied soils were loamy sand (n= 3), sandy loam (n= 39), loam (n= 20), silt loam (n= 22), silty clay loam (n= 7), silty clay (n= 7) and clay (n= 2). There were found significant correlations between soil saturated hydraulic conductivity (Ks) and sand (r= 0.60**), silt (r= -0.60**), clay (r= -0.43**), organic carbon (r= 0.36**), bulk density (r= -0.52**), particle density (r= -0.53**), total porosity (r= 0.31**), CaCO3 (r= -0.58**), mean geometric diameter (r= 0.57**), SAR (r= -0.35**), EC (r= -0.22*) and pHe (r= -0.44**). These results are in line with the findings of the former studies that reported direct relation of Ks with OC, sand, and inverse relation of Ks with silt, clay, BD, and SAR. Generally, 8 regression, artificial neural network, and neuro-fuzzy pedotransfer functions were constructed to estimate soil saturated hydraulic conductivity (Ks) from measured readily available soil variables. The results of the best regression, artificial neural network and neuro-fuzzy pedotransfer functions indicated that the most suitable input variables to estimate soil saturated hydraulic conductivity (Ks) were bulk density and silt in the studied region. The values of R2, RMSE and ME were obtained equal to 0.65, 0.119 cm min-1, 0.059 cm min-1 and 0.73, 0.087 cm min-1, 0.006 cm min-1 and 0.69, 0.127 cm min-1, -0.051 cm min-1 for the best regression, artificial neural network, and neuro-fuzzy Ks pedotransfer functions, respectively. According to these results, the ANN PTF was the best in estimating Ks because of having high R2 and low RMSE compared with regression and neuro-fuzzy PTFs. The former researchers also obtained bulk density and silt as the best input variables for estimating soil saturated hydraulic conductivity (Ks) in different soils and regions.
4- Conclusions
The results showed that bulk density and silt are the most suitable readily available soil variables to estimate soil saturated hydraulic conductivity (Ks) in the studied salt-affected soils. According to the RMSE criterion, the precision of an/the artificial neural networks in estimating Ks was more than regression and neuro-fuzzy pedotransfer functions in this research. Also, regression and neuro-fuzzy PTFs have not an/the observable difference in estimating Ks.
Machine summary:
1 Pedotransfer Functions (PTF) 2 Artificial Neural Network (ANN) 3 Neuro-Fuzzy 193 مدل هاي مختلفي براي تخمين Ks ارائه شده است ، ولي با عنايت به اينکه توابع انتقالي خاک به مکان وابسته هسـتند؛ از اين رو اين پژوهش به منظور اشتقاق توابع رگرسيوني، شبکه عصبي مصنوعي و نروفازي، براي تخمين Ks خاکهاي متأثر از نمک منتخب از شمال غرب درياچه ي اروميه و مقايسه ي دقت برآورد مدل هاي مذکور انجام شد.
119 جدول ٥: معيارهاي ارزيابي توابع شبکه عصبي مصنوعي براي برآورد هدايت هيدروليکي اشباع (Ks) خاک داده هاي آموزشي داده هاي آزموني شماره تابع متغيرهاي ورودي R2 (cRmM/mSiEn) (cmM/mEin) (cmM/mEin) (cRmM/mSiEn R 2 0/71 0/529 0/436 0/47 0/531 0/454 sand, BD 1 0/73 0/087 0/006 0/48 0/183 0/013 silt, BD 2 0/69 0/099 -0/010 0/47 0/189 -0/013 sand, BD, pH 3 0/76 0/097 -0/008 0/49 0/184 0/010 silt, BD, pH, EC, OC 4 0/74 0/120 0/037 0/49 0/208 0/056 silt, BD, pH, EC 5 0/71 0/098 -0/028 0/56 0/168 0/002 silt, BD, pH, EC, SAR 6 0/72 0/096 -0/026 0/47 0/190 -0/044 CaCO3, BD 7 0/74 0/094 -0/008 0/47 0/189 -0/014 BD, CaCO3, d 8 g توصيف متغيرها و علائم در جداول ١ و ٤ آورده شده است .
Therefore, many different regressions, artificial neural network (ANN) and the neuro-fuzzy pedotransfer functions (PTF) have been developed to estimate Ks from readily available soil variables such as sand, silt, clay, bulk density (BD), particle density (PD), electrical conductivity (EC), pH, and organic carbon (OC).