خلاصة:
پیشینه و هدف دشت آققلا که یکی از دشتهای حاصلخیز استان گلستان است. تجمع رواناب و شکلگیری سیلاب خسارات جانی و مالی فراوان را به همراه دارد. شناخت پراکنشهای زمانی و مکانی رطوبت سطح خاک یک متغیر کلیدی در شبیهسازی رواناب و کاهش مخاطرات ناشی از سیل در محدوده تحقیق است. در این تحقیق میزان رطوبت سطح خاک با استفاده از دادههای ماهوارههای راداری آلوس-2 از پالسار در باند L و سنتینل-1 در باند C مورد بررسی قرار گرفت.مواد و روش ها در این تحقیق، از سه سری داده راداری آلوس-2 از سنجنده پالسار، سنتینل-1 و داده اپتیکی سنتیل-1 استفاده شد. برای محاسبه رطوبت سطح خاک ابتدا به روش وزنی در محل هر خوشه نمونههای خاک در عمق 5 سانتیمتری برداشت و توزین گردید پس از خشک نمودن نمونهها در آون وزن خشک و مرطوب محاسبه گردید. در ادامه دو پارامتر مهم رطوبت خاک شامل ثابت دیالکتریک و زبری سطح خاک اندازهگیری شد. ثابت دیالکتریک با استفاده از TDR و زبری با استفاده از دو دوربین دیجیتال محاسبه گردید. با استفاده از نرمافزار Agisoft photo scan و ArcGIS مدل رقومی ارتفاعی و پروفیل ناهمواریهای سطح استخراج گردید. دادههای راداری با نرمافزار SNAP مربوط به آژانس فضایی اروپا ESA پردازش گردید. مراحل پردازش شامل کالیبراسیون رادیومتریکی، استخراج سیگما نات در واحد دسیبل، فیلترگذاری با فیلتر لی و کرنل سایز 5×5 و تصحیح هندسی و نهایتاً ژئوکدینگ هر سه داده انجام شد. با داده اپتیکی سنتینل-2 ضمن تهیه شاخصهای گیاهی و نمناکی نقشه کاربری اراضی با الگوریتم راندوم فارست تهیه گردید. در ادامه با در دست داشتن مقادیر زبری، ثابت دیالکتریک و شاخصهای گیاهی و نمناکی در محیط نرمافزار متلب 2019 دو مدل Oh و WCM کدنویسی و مقادیر رطوبت سطح خاک محاسبه و خروجی گرفته شد. سپس رطوبت حاصل از مدل باز پراکنش سطحی با رطوبت وزنی در محل و مطالعات پیشین مورد ارزیابی قرار گرفت.نتایج و بحث در بخش حساسیت سیگنال رادار به رطوبت سطح خاک در مدل Oh در پلاریزاسیون HH و HV تصاویر پالسار-2، در با 0.86=R2، با 0.905=R2 و برای تصاویر سنتینل-1، در با 0.91=R2 و با 0.93=R2 به رطوبت سطح خاک حساس است. همچنین حساسیت باز پراکنش سطحی مدل Oh در پالسار-2 در پلاریزه HH با 0.96=R2 و در پلاریزه HV با 0.94=R2 و برای داده سنتینل-1 در پلاریزه VV با 0.72=R2 و در پلاریزه VH با 0.74=R2 به زبری سطح حساس است. رابطۀ بین کل ضرایب باز پراکنش سطحی حاصل از مدل ابر آب و رطوبت اندازهگیری شده در پلاریزاسیون HH و HV تصاویر پالسار-2، در و به ترتیب با 0.73=R2 و 0.399=R2 و برای دادۀ سنتینل-1 در پلاریزاسیون VV و VH، در و بهترتیب با 0.073=R2 و 0.13=R2 به رطوبت سطح خاک حساس است. در بخش برآورد رطوبت خاک مبتنی بر مدل Oh برای دادۀ پالسار-2 و سنتینل-1 در اراضی لخت در پلاریزاسیون HH و HV برای تصاویر پالسار-2، و به ترتیب با 1.924=RMSE، 0.88=R2 و 1.45=RMSE، 0.93=R2 و برای داده سنتینل-1، و بهترتیب با 2.57=RMSE، 0.81=R2 و 0.90=RMSE، 0.97=R2 است. نتایج مدل Oh نشان داد که در پلاریزاسیون HV تصاویر پالسار-2 و پلاریزه VH سنتینل-1 نتایج رطوبت خاک قویتر است. نتایج رطوبت خاک در مدل WCM در اراضی دارای پوشش گیاهی برای تصاویر پالسار-2 در پلاریزاسیون HH و HV در و بهترتیب با 2.44=RMSE، 0.81=R2 و 2.64=RMSE، 0.78=R2 و برای Sentinel-1 در و بهترتیب با 2.74=RMSE، 0.75=R2 و 2.69=RMSE، 0.76=R2 است. در بحث برآورد رطوبت خاک با مدل WCM نتایج پالسار-2 در پلاریزه HH با 0.81=R2 در اراضی با تراکم پوشش گیاهی بیشتر از 0.4 درصد دارای نتیجه بهتری است.نتیجه گیری در مدل Oh نتایج تحقیق نشان داد که حساسیت پلاریزه HV در تصاویر پالسار-2 و پلاریزه VH در سنتینل-1 در مناطقی که دارای پوشش گیاهی با تراکم کمتر از 0.4 است بسیار بیشتر از حساسیت پلاریزه HH آلوس-2 و پلاریزه VV سنتینل-1 است. همچنین حساسیت پلاریزه HH تصاویر پالسار-2 و پلاریزه VH سنتینل-1 به زبری سطح بیشتر از سایر پلاریزاسیونها است. در مدل ابر آب نتایج حاصل از داده سنتینل-1 در حذف تأثیرات پوشش گیاهی و شبیهسازی ضرایب باز پراکنش سطحی بسیار ضعیف است، دلیل این امر در نظر نگرفتن پارامتر زبری توسط مدل ابر آب و همچنین بالا بودن سطح رطوبت در محدوده مطالعاتی بوده که باعث کاهش نفوذ موج گردیده است. نتایج تحقیق نشان داد که حساسیت پلاریزه HH در مناطقی که دارای پوشش گیاهی با تراکم بیش از 0.4 است بسیار بیشتر از حساسیت پلاریزه HV بوده، همچنین حساسیت سیگنال راداری به نویز در پلاریزه HH کم ولی در پلاریزه HV زیاد بود. این موضوع نشان داد که قطبش همسان HH نسبت به پوشش گیاهی حساس است، بنابراین مدل در حضور پوشش گیاهی از ثبات بیشتری برخوردار خواهد بود. مقایسه دو داده SAR در محدوده نشان داد، تصاویر آلوس-2 به دلیل طولموج بلند و عمق نفوذ بیشتر در خاک و پوشش گیاهی، برای برآورد رطوبت سطح خاک در محدوده تحقیق و مناطق مشابه کارایی بیشتری دارد. همسو با موضوع پژوهش تحقیقات دیگری صورت گرفته ازجمله آلیشان و همکاران در پژوهشی اقدام به برآورد رطوبت سطح خاک در اراضی بایر با استفاده از دادههای آلوس-2 و سنتینل-1 کردهاند که به این نتیجه رسیدند که مدل WCM در حذف محتوی آب گیاهی و برآورد رطوبت سطح با استفاده از معکوس مدل در هر دو داده آلوس-2 و سنتینل-1 بالا بوده است. در تحقیقی که توسط زریبی و همکاران صورت گرفته بود نشان دادند که مدل Oh در دو پلاریزه HV و HH با RMSE کمتری قادر به برآورد رطوبت سطح خاک است. این تحقیق برای مدیریت محیطزیست، کشاورزی، منابع طبیعی و مدیریت منابع آب مناسب بوده و کمک شایانی به مدلسازیهای هیدرولوژیکی می نماید.
Background and Objective
Aq Qala plain is one of the fertile plains of Golestan province. Runoff accumulation and flood formation cause a lot of human and financial losses. Understanding the temporal and spatial distributions of soil surface moisture is a key variable in simulating runoff and reducing flood hazards in the research area. In this study, the soil moisture content was investigated using data from ALOS PALSAR-2 radar satellites in the L band and Sentinel 1 in the C band.
Materials and Methods
In this research, three series of ALOS PALSAR-2 radar data, Sentinel 1, and optical data of Sentinel 2A were used. To calculate the soil surface moisture, first soil samples were taken and weighed at a depth of 5 cm in each cluster by weight method. After drying the samples in the oven, the dry and wet weight was calculated. Then two important parameters of soil moisture, including dielectric constant and soil surface roughness were measured. The dielectric constant was calculated using TDR and the roughness was calculated using two digital cameras. Using Agisoft photo scan and Arc GIS software, a digital elevation model and surface roughness profiles were extracted. Radar data were processed with SNAP software from the European Space Agency (ESA). The processing steps included radiometric calibration, sigma note extraction in decibels, filtration of Lee filter and 5x5 kernel size, geometric correction and finally geocoding of all three data. With Sentinel-2 optical data, while preparing vegetation and moisture indices, a land-use map was prepared with the Forrest random algorithm. Then, with the values of roughness, dielectric constant and plant and moisture indices in MATLAB 2019 software environment, two models Oh, and WCM were coded and the soil moisture values were calculated and taken out. Then the moisture obtained from the surface redistribution model with gravity, moisture at the site, and previous studies was evaluated.
Results and Discussion
In terms of radar signal sensitivity to soil surface moisture in Oh model in HH and HV polarization ALOS PALSAR- 2, in with R2 = 0.86, with R2 = 0.905 and for Sentinel-1 images, in with R2 = 0.91 and with R2 = 0.93 are sensitive to soil surface moisture. The surface backscattering coefficients sensitivity of the Oh model in ALOS PALSAR-2 in HH polarization and in HV polarization and for Sentinel-1 data in VV polarization and in VH polarization is sensitive to surface roughness. Also, the surface backscattering coefficients sensitivity of the Oh model in ALOS PALSAR-2 in HH polarization with R2 = 0.96 and in HV polarization with R2 = 0.94 and for Sentinel-1 data in VV polarization with R2 = 0.72 and in VH polarization with R2 = 0.74 is sensitive to surface roughness. Relationship between total surface backscattering coefficients obtained from the Water cloud model and soil surface moisture measured in HH and HV polarization PALSAR-2 images, in and with R2 = 0.73 and R2 = 0.399, respectively, and for Sentinel-1 data in polarization VV and VH, in and with R2 = 0.073 and R2 = 0.13, respectively, are sensitive to soil surface moisture. In the soil moisture estimation section based on Oh model for PALSAR-2 and Sentinel-1 data on bare lands in HH and HV polarization for PALSAR-2, and with R2 = 0.88, RMSE = 1.924 and R2 = 0.93, RMSE = 1.45 for Sentinel-1 data, respectively, and are with R2 = 0.81, RMSE = 2.57 and R2 = 0.97, RMSE = 0.90 , respectively. The results of Oh model showed that in HV polarization, Pulsar-2 images and VH sentinel-1 polarization images have stronger soil moisture results. Soil moisture results in WCM model in vegetated lands for PALSAR-2 images in HH and HV polarization in and , respectively, with R2 = 0.81, RMSE = 2.44 and R2 = 0.78, RMSE = 2.64 , R2 = 0.88, RMSE = 1.924 for Sentinel-1 in and with R2 = 0.75, RMSE = 2.74 and R2 = 0.76, RMSE = 2.69 , respectively. In the discussion of estimating soil moisture with WCM model, the results of PALSAR-2 in polarization HH with R2 = 0.81 in lands with vegetation density of more than 0.4% have better results.
Conclusion
In the Oh model, the results showed that the sensitivity of HV polarization in PALSAR-2 images and VH polarization in Sentinel-1 in an area with vegetation density less than 0.4 is much higher than the polarization sensitivity of HH PALSAR-2 and polarized VV Sentinel-1. Also, the HH polarization sensitivity of PALSAR-2 and VH Sentinel-1 polarized images to surface roughness is higher than other polarizations. In the Water cloud model, the results obtained from Sentinel-1 data by eliminating the effects of vegetation and simulating the surface backscattering coefficients are very poor. Has reduced the penetration of the wave. The results showed that the polarization sensitivity of HH in areas with vegetation with a density of more than 0.4 was much higher than the polarization sensitivity of HV, also the sensitivity of the radar signal to noise was low in polarized HH but high in polarized HV. This showed that the polarization of the HH matrix is sensitive to vegetation, so the model will be more stable in the presence of vegetation. A comparison of two SAR data in the range showed that PALSAR-2 images are more efficient for estimating soil surface moisture in the study area and similar areas due to their long wavelength and depth of penetration into soil and vegetation. In line with the research topic, other researchers, including Alishan et al. A study has attempted to estimate soil surface moisture in barren lands using PALSAR-2 and Sentinel-1 data, which concluded that the WCM model eliminates plant water content. And surface moisture estimation using the inverse model was high in both PALSAR-2 and Sentinel-1 data. A study by Zaribi et al. Showed that the OH model in both polarized HV and HH with lower RMSE is able to estimate soil surface moisture. This research is suitable for environmental management, agriculture, natural resources, and water resources management and contributes significantly to hydrological modeling.