1. معرفی
رطوبت خاک یک پارامتر مهم در تعادل انرژی جهانی و چرخه آب است. علاوه بر این، نقش مهمی در مدیریت منابع آب، رشد پوشش گیاهی، پایش سیل و پیشبینی آب و هوا دارد [ 1 ، 2 ، 3]]. سنجش از دور یکی از موثرترین ابزارهای مورد استفاده در پایش رطوبت خاک به دلیل طیف وسیع مشاهدات و بسامد بالای اندازه گیری های مکرر است. در حال حاضر، سنجش از دور مایکروویو روشی موثر برای اندازهگیری رطوبت خاک در نظر گرفته میشود، زیرا میتواند به ابرها و پوشش گیاهی نفوذ کند و در هر شرایط آب و هوایی، به ویژه در باند L، هم در روز و هم در شب کار میکند. رادیومترهای مایکروویو باند L پتانسیل بالایی برای نظارت بر رطوبت خاک نشان می دهند زیرا تضعیف کمی در جو را تجربه می کنند، قابلیت نفوذ پوشش گیاهی قوی دارند و حساسیت بالایی به رطوبت خاک دارند. بنابراین، بسیاری از ماموریت های ماهواره ای از پرتو سنج های باند L برای پایش رطوبت خاک استفاده می کنند. ماموریت رطوبت خاک و شوری اقیانوس (SMOS) از یک رادیومتر باند L چند زاویه ای استفاده می کند. 1]]، در حالی که ماموریت های Soil Moisture Active and Passive (SMAP) و Aquarius هر دو دارای رادار و رادیومترهای باند L هستند [ 2 ، 4 ]. ماموریت مشاهده چرخه آب برنامه ریزی شده (WCOM) از یک رادیومتر سه فرکانسی که شامل باند L [ 5 ] است استفاده می کند.
الگوریتم های زیادی برای پرتو سنج های باند L ایجاد شده است. مدل انتشار سطحی جزء حیاتی این الگوریتمها است. عمدتاً زبری سطح را بر اساس تابش سطح در نظر می گیرد. از آنجایی که اندازهگیری زبری سطح در مقیاسهای بزرگ دشوار است، مدلهای متعددی برای تخمین اثر زبری سطح، از جمله مدلهای فیزیکی و نیمه تجربی، پیشنهاد شدهاند. مدل های فیزیکی مانند مدل معادلات انتگرال (IEM) و مدل معادلات انتگرال پیشرفته (AIEM) [ 6] بسیار پیچیده هستند و شامل بسیاری از پارامترهای غیر قابل اندازه گیری هستند. اگرچه این مدلها توضیحات مفصلی از فرآیندهای تابشی ارائه میدهند، پیچیدگیهای آنها کاربردهای آنها را محدود میکند. در مقابل، مدلهای نیمه تجربی از پارامترهای سادهسازی شدهای استفاده میکنند که به راحتی قابل بازیابی هستند. تقریباً تمام الگوریتمهای مرتبط با بازیابی رطوبت خاک از مدلهای نیمه تجربی توسعه یافتهاند. را Hpاچپمدل توسط Choudhury [ 7 ] پیشنهاد شد که مؤلفه غیر منسجم را در فرکانس پایین نسبتاً کوچک در نظر گرفت. این مدل معمولاً در بازیابی رطوبت خاک باند L استفاده می شود و در بسیاری از مدل های نیمه تجربی دیگر توسعه یافته است. با این حال، پارامترهای زبری سطح در آن مدلها عمدتاً به صورت مقادیر ثابت با توجه به مناطق مختلف بیان میشوند و عدم قطعیتها را به دلیل ناهمگونی سطح قوی افزایش میدهند. Shi [ 8 ] یک مدل پارامتری برای باند L بر اساس پایگاه داده شبیه سازی شده IEM ارائه کرد. با این حال، این مدل به اندازه کافی برای کاربردهای عملی ساده نیست. اخیراً برخی مدلها با تجزیه و تحلیل بازتاب سطح دو قطبی با استفاده از پایگاه داده شبیهسازیشده AIEM ایجاد شدهاند [ 9 ، 10 ،11 ، 12 ]. این روشهای جدید از قطبش دوگانه برای کاهش وابستگی به پارامترهای زبری سطح استفاده میکنند و الگوریتم را قابل اعتمادتر میکنند. مدل انتشار سطحی ارائه شده در این مطالعه بر اساس این روش ها است.
پس از در نظر گرفتن اثرات زبری سطح، رطوبت خاک را می توان با استفاده از سه روش اساسی محاسبه کرد: یک الگوریتم تکرار معکوس، یک الگوریتم عددی رو به جلو و یک الگوریتم تجربی. الگوریتم تکرار معکوس معمولاً در بسیاری از الگوریتم های تک کانال استفاده می شود [ 13 ]. ابتدا رطوبت خاک با توجه به مناطق مختلف به مقدار اولیه تنظیم می شود. سپس بازتاب فرنل بر اساس معادله بازتاب فرنل به صورت معکوس محاسبه می شود. در نهایت، مقدار شبیه سازی شده بر اساس مقدار مشاهده شده با استفاده از الگوریتم حداقل مربعات تا همگرایی تکرار می شود. در رادیومترهای دو کاناله و چند کاناله، الگوریتم عددی رو به جلو به طور گسترده ای مورد استفاده قرار می گیرد. رطوبت خاک با به حداقل رساندن یک تابع هزینه با استفاده از الگوریتم تکراری حداقل مربعات بازیابی می شود.14 ]. این روش در بهترین استفاده از داده های چند کاناله به خوبی عمل می کند. الگوریتم تجربی شامل چندین معادله تجربی است که روابط نامشخص پارامترها را توصیف می کند. با توجه به داده های آموزشی کافی، این الگوریتم دقت خوبی را در مناطق خاص ارائه می دهد [ 15 ، 16 ، 17 ، 18 ]. با این حال، نقاط ضعف این مدل ها هنوز کاملاً مشهود است. در الگوریتم تکرار معکوس، اجتناب از پدیده شبه همگرایی دشوار است. در مورد الگوریتم عددی رو به جلو، تابع هزینه همیشه بر اساس چندین پارامتر بازیابی شده (رطوبت خاک، عمق نوری در نادر و پارامتر زبری خاک) است [ 14] .]. همه این پارامترها در به حداقل رساندن تفاوت بین مقادیر اندازه گیری شده و شبیه سازی شده تأثیر دارند. بنابراین، الگوریتم تکراری نمی تواند توضیح دهد که کدام پارامتر مسئول تغییر اندازه گیری ماهواره در مکانیسم است [ 19 ]. محدودیت ها و مقادیر اولیه تعیین شده در تابع هزینه نیز عوامل عدم قطعیت را به الگوریتم اضافه می کنند. هم در تکرار معکوس و هم در الگوریتم عددی رو به جلو، مدل دی الکتریک خاک برای تخمین اثرات خاک اتخاذ میشود که دادههای جانبی قابلتوجهی را به ارمغان میآورد و عدم قطعیتهای بیشتری را به الگوریتم وارد میکند. الگوریتم تجربی فاقد اهمیت فیزیکی واضح است و کاربرد آن را در مقیاس های بزرگ دشوار می کند. با توجه به محدودیت های این مدل ها، الگوریتم جدیدی در این مقاله پیشنهاد شده است.
هدف از این مطالعه توسعه یک الگوریتم جدید است که به راحتی می تواند برای بازیابی رطوبت خاک باند L استفاده شود. الگوریتم مبتنی بر فیزیکی دو الگوریتم فرعی را ترکیب می کند: یک مدل انتشار سطحی و یک مدل بازیابی رطوبت خاک. مدل انتشار سطح به اثرات زبری سطح می پردازد. این مدل بر اساس بازتاب سطح دو قطبی بدون نیاز به پارامترهای زبری سطح است. مدل بازیابی رطوبت خاک بر خواص دی الکتریک خاک متمرکز است. این مدل بازیابی از یک راه حل تحلیلی استفاده می کند که بر اساس مدل فیزیکی است و تنها به محتویات ماسه و خاک رس در طول فرآیند بازیابی نیاز دارد. پیشرفت اصلی این الگوریتم این است که پارامترهای نامشخص کمتری دارد و از نظر فیزیکی مبتنی است و وابستگی منطقهای ندارد.
بقیه این مقاله به شرح زیر سازماندهی شده است. بخش 2 منطقه مورد مطالعه و مجموعه داده های میدانی مورد استفاده در این مطالعه را توصیف می کند. بخش 3 الگوریتم پیشنهادی در این مطالعه را معرفی می کند. بخش 4 نتایج تجربی و تحلیل حساسیت مدل جدید را ارائه میکند. نتیجه گیری مقاله در بخش 5 خلاصه شده است .
2. منطقه مطالعه و مجموعه داده ها
2.1. منطقه مطالعه
از سال 1979 تا 1982، مرکز پرواز فضایی ناسا/گودارد (GSFC) و وزارت کشاورزی ایالات متحده (USDA)/مرکز تحقیقات کشاورزی Beltsville (BARC) آزمایشهای چهار ساله را در سایتهای آزمایش BARC انجام دادند [20 ، 21 ، 22 ، 23 . ]. هدف از این آزمایش ها بررسی اثرات بافت خاک، ناهمواری سطح و پوشش گیاهی بر تخمین از راه دور میزان رطوبت خاک با استفاده از رادیومترهای مایکروویو [ 21] بود.]. این مطالعه در چندین سایت آزمایش کشاورزی محلی در بلتسویل، مریلند، که شامل انواع خاک و پوشش گیاهی مختلف بود، انجام شد. مشاهدات رطوبت خاک، دمای خاک، زیست توده گیاهی و سایر پارامترهای خاک و تاج پوشش همزمان با اندازهگیریهای مایکروویو [ 23 ] انجام شد. این داده ها از آن زمان به طور گسترده ای برای اعتبارسنجی مدل های مختلف بازیابی رطوبت خاک استفاده شده اند [ 9 ، 10 ]. در این مطالعه، محل های آزمایش سطح برهنه به عنوان منطقه مورد مطالعه انتخاب شدند.
2.2. داده های مایکروویو
داده های مایکروویو با رادیومترهای باند L (1.4 گیگاهرتز) نصب شده بر روی یک کامیون بوم به دست آمد. سنسورها رادیومترهای دیک دو قطبی بودند که انتشار امواج مایکروویو حرارتی را در زوایای فرود مختلف از 10 درجه تا 70 درجه در مراحل 10 درجه اندازهگیری میکنند. کالیبراسیون سیستم مایکروویو هر روز با استفاده از یک بازتابنده آسمان سرد و یک ماده جاذب مایکروویو (Eccosorb) تأیید شد. دقت تخمین زده شده رادیومتر باند L تقریباً 3± K [ 23 ] است.
2.3. داده های میدانی
این آزمایش چهار ساله در سایت های آزمایش BARC انجام شد. در سالهای 1979 و 1980، مزرعهای با خاک لوم شنی در Elinsboro انتخاب شد. دو سایت آزمایشی جدید در سالهای 1981 و 1982 انتخاب شدند: یکی در ادمونستون و دیگری در مزرعه جنوبی. در سال 1982، سه سایت انتخاب شد: دو سایت قبلی به اضافه یک سایت در گیش. داده های صحرایی از سطوح لخت اندازه گیری شد و شامل چگالی ظاهری خاک، رطوبت خاک و دمای خاک بود. اندازهگیریهای دقیق چگالی ظاهری و بافت خاک چندین بار در هر محل آزمایش در طول آزمایشها انجام شد [ 23 ]. توضيحات تفصيلي تمام سايت هاي سطح برهنه در جدول 1 فهرست شده است. رطوبت خاک در عمق های 0-0.5 سانتی متر، 0-2.5 سانتی متر، 2.5-5 سانتی متر و 5-10 سانتی متر اندازه گیری شد. پروفیل دمای خاک در چهار عمق 1.25 سانتی متر، 2.5 سانتی متر، 7.5 سانتی متر و 15 سانتی متر نمونه برداری شد. دادههای رطوبت و دما خاک در محل به طور همزمان از اندازهگیریهای رادیومتر مایکروویو در یک میدان مشخص بهدست آمد [ 20 ].
3. روش شناسی
3.1. مدل انتشار سطحی
مدل انتشار سطحی، رابطه بین بازتاب سطح موثر و بازتاب فرنل را با پارامترسازی اثرات مختلف از زبری سطح توصیف میکند. با توجه به اینکه جزء غیر منسجم در فرکانس های پایین ناچیز است، Hp�پمدل توسعه داده شد. این مدل به طور گسترده در الگوریتم های مختلف رادیومتر باند L استفاده شده است [ 14 ]. با این حال، پارامترهای زبری سطح در این مدل به طور تجربی با توجه به مناطق مختلف اختصاص داده شده است. این مقادیر نامشخص خطاهای قابل توجهی را در سطوح ناهمگن ایجاد می کند. اخیراً، چن [ 9 ] یک مدل انتشار سطحی ساده شده با توجه به پایگاه داده شبیه سازی شده AIEM توسعه داده است:
جایی که آرپ��و آرq��بازتاب های موثر هستند. rپ��و rq��بازتاب های فرنل هستند. θ�زاویه بروز است. س�ارتفاع ریشه میانگین مربع (RMS) است. ل�طول همبستگی است. و ضرایب آ�، ب�، و سی�فقط به قطبش و زاویه تابش وابسته هستند، که می توان از رابطه (3) با استفاده از پایگاه داده شبیه سازی شده AIEM در مطالعه قبلی به دست آورد.
این مدل نسبت را انتخاب می کند s / l�/�به عنوان پارامتر زبری سطح بر اساس اچپ��مدل. پارامتر سی�برای اصلاح حساسیت انتشار سطحی به رطوبت خاک استفاده می شود [ 9 ]. با اينكه س�و ل�نمی توان به راحتی از سطح اندازه گیری کرد، این شکل از مدل پارامترهای زبری سطح را در قطبش V- و H یکی می کند، که تا حدی اثرات زبری سطح را از طریق قطبش دوگانه حذف می کند. با حذف متغیر زبری سطح s / l�/�از معادلات (1) و (2)، رابطه بین بازتاب موثر و بازتاب فرنل را می توان به سادگی به دست آورد:
در حالی که تنها چهار پارامتر در این معادله وجود دارد، حل آن به صورت مستقیم آسان نیست. سمت چپ معادله (4) شامل بازتاب موثر است آرپ��و آرq��، که می تواند از دمای روشنایی دو قطبی محاسبه شود. بازتاب های فرنل rپ��و rq��دو کمیت مجهولی هستند که این معادله را نادیده میگیرند.
در سالهای اخیر، با توجه به دادههای شبیهسازیشده، ثابت شده است که رابطه بین بازتابپذیریهای فرنل در قطبش دوگانه یک تابع توان است. هونگ یک مشتق تقریبی را از طریق یک سری تیلور ارائه کرد که توسط مطالعات حساسیت و استخراج تحلیلی تأیید شد [ 24 ]:
تقریب هنگ در کاربردهای مختلف سنجش از دور اعمال شده است [ 25 ، 26 ، 27 ]. در این تحقیق می توان رابطه (4) را با این رابطه تقریبی ساده کرد. با وارد کردن معادلات (5) و (6) در معادله (4)، rپ��را می توان جایگزین کرد rq��. سپس معادله (4) را می توان به صورت زیر بیان کرد:
هفت ضریب را می توان در سه ضریب جدید ترکیب کرد و یک رابطه واضح بین بازتاب موثر دو قطبی و بازتاب فرنل ایجاد کرد. rq��:
روابط بین ضرایب در معادلات (7) و (8) به شرح زیر است:
جایی که آ�، ب�، و ج�می توان از تجزیه و تحلیل رگرسیون با توجه به پایگاه داده شبیه سازی شده AIEM تعیین کرد.
3.2. مدل بازیابی رطوبت خاک
برای سطوح خالی، معادله انتقال تابش را می توان به صورت زیر ساده کرد:
جایی که تیبپ���و تیبq���درجه حرارت روشنایی مایکروویو قطبش دوگانه و تیه�هدمای سطح موثر است که می تواند از مدل های نیمه تجربی استخراج شود یا با دمای سطح جایگزین شود.
ترکیب معادلات (12) و (13) در معادله (8)، بازتاب فرنل rq��را می توان به صورت زیر بیان کرد:
شکل دیگری از معادله بازتاب فرنل توسط لیو [ 28 ] پیشنهاد شد. در این معادله قسمت واقعی دی الکتریک خاک ε�با ضریب شکست واقعی تنظیم شده جایگزین می شود نrن�. ضریب شکست واقعی تعدیل شده دارای معنای فیزیکی مشابه دی الکتریک خاک است که به خواص خاک نیز بستگی دارد. بنابراین، این شکل از معادله بازتاب فرنل را می توان برای بازیابی رطوبت خاک اعمال کرد. شکل کامل معادله به صورت زیر نوشته می شود:
جایی که نrن�ضریب شکست واقعی و ضریب شکست تنظیم شده است n =nr– منnمن�=��–من�من. nr��. اینجا، nمن�منبه ترتیب قسمت های واقعی و موهومی ضریب شکست هستند که می توان از ثابت دی الکتریک پیچیده بدست آورد:
جایی که ε“�“و ε” “�““به ترتیب اجزای واقعی و خیالی ثابت دی الکتریک مختلط هستند.
حل تحلیلی معادله (15) برای یک ماده غیر جاذب توسط Sohn [ 29] ارائه شد. ] ارائه شد. در این مطالعه، شکل مشابهی از محلول به محیط های جاذب تعمیم داده می شود.
با گرفتن جذر معادله (15)، r1/2 _ _q��1/2را می توان به صورت زیر نوشت:
جایی که نr= گناهθمن/ گناهθتین�=گناه�من/گناه�تی، θمن�منزاویه بروز است و θتی�تیزاویه شکست است. طبق فرآیند انتقال تابشی، موج از طریق هوا به خاک یا آب نفوذ می کند. هوا در مقایسه با خاک و آب از نظر نوری محیطی رقیقتر است. بر اساس قانون اسنل، نr> 1ن�>1. با حذف قدر مطلق، سمت راست معادله (19) را می توان به صورت زیر ساده کرد:
پس از گرفتن جذر دو طرف معادله (20)، می توان آن را به شکل زیر بازنویسی کرد:
با گرفتن مربع های دو طرف رابطه (21) رابطه (22) بدست می آید:
موارد مشابه را می توان به صورت زیر ادغام کرد:
با گرفتن مربع دو طرف رابطه (23) معادله (24) به دست می آید:
ترتیب معادله (24) معادله (25) را به دست می دهد:
جمله دوم سمت چپ صفر است زیرا نr> 1ن�>1. بنابراین، معادله (25) را می توان به صورت زیر بازنویسی کرد:
پس از ساده سازی، ضریب شکست واقعی تنظیم شده است نrن�را می توان به صورت زیر بیان کرد:
در این مطالعه، رابطه بین ضریب شکست واقعی تعدیل شده و رطوبت خاک پیشنهاد شده است. شکل مدل شبیه به مدل ثابت دی الکتریک هالیکاینن [ 30 ] است که به شرح زیر است:
جایی که S�و C�به ترتیب درصد محتوای ماسه و خاک رس هستند، ε′�′بخش واقعی ثابت دی الکتریک مختلط است و mv��رطوبت حجمی خاک است.
مدل جدید بازیابی رطوبت خاک با جایگزینی سمت چپ مدل هالیکاینن توسعه یافته است. ε′�′به ضریب شکست واقعی تعدیل شده Nr��:
ضرایب این مدل ( a0�0، a1�1، a2�2، b0�0، b1�1، b2�2، c0�0، c1�1، و c2�2) را می توان از طریق تحلیل رگرسیون با توجه به پایگاه داده شبیه سازی شده ثابت دی الکتریک تعیین کرد.
معادله درجه دوم با یک مجهول در معادله (29) دو راه حل دارد که یکی از آنها رطوبت خاک است. تنها داده های کمکی مورد نیاز این مدل درصد محتوای شن و رس است که می تواند از تحقیقات میدانی به دست آید. سپس می توان رطوبت خاک را با استفاده از مدل دی الکتریک خاک از طریق یک کمیت ناشناخته محاسبه کرد Nr��. پس از ساده سازی، رطوبت حجمی خاک را می توان به صورت زیر بیان کرد. طرح انتگرال این روش در شکل 1 نشان داده شده است :
4. نتایج و بحث
4.1. تخمین رطوبت خاک با داده های شبیه سازی شده
برای کالیبره کردن ضرایب مدل انتشار سطحی (معادله (8))، ما یک پایگاه داده انتشار سطحی شبیهسازی شده تحت پارامترهای حسگر فرکانس SMOS، از جمله 1.41 گیگاهرتز و قطبش V و H، با استفاده از مدل AIEM تولید کردیم. این پایگاه داده طیف وسیعی از زوایای فرود، رطوبت حجمی خاک و پارامترهای زبری، از جمله ارتفاع RMS، طول همبستگی و توابع همبستگی را پوشش می دهد. پارامترهای دقیق مورد استفاده در پایگاه داده شبیه سازی شده در جدول 2 خلاصه شده است خلاصه شده است .
سه ضریب ( a�، b�، و c�) در زوایای فرود مختلف در معادله (8) در جدول 3 آورده شده است . شکل 2 بازتاب فرنل را مقایسه می کند rq��با استفاده از مدل AIEM با مدل های تخمین زده شده با استفاده از معادله (6) به عنوان تابعی از زوایای فرود از 5 درجه تا 60 درجه محاسبه شده است. ریشه میانگین مربعات خطا (RMSEs) در زوایای بروز بزرگ، به ویژه در 45 درجه، کمتر است. این تمایل ناشی از سبک مدل ساده شده و تقریب شیب برای زبری است. پارامتر s/l�/�استفاده شده در معادلات (1) و (2) زبری سطح را در زوایای فرود بزرگ، به ویژه در 45 درجه، بهتر تخمین می زند. در زوایای کوچک، هر دو قطبش H و V دقت پایینی دارند. هنگامی که زاویه تابش بیشتر از 50 درجه باشد، قطبش H خطاهای قابل توجهی ایجاد می کند [ 9 ]. مدل ساده شده بر اساس این فرضیه است که مقادیر مشابهی به دست آورد. به طور کلی، خطای مدل انتشار سطحی جدید در تمام زوایای بروز کوچکتر از 0.02 است، که به وضوح نشان می دهد که این مدل انتشار ساده شده می تواند جایگزین مدل AIEM در طیف گسترده ای از شرایط سطحی بدون خطاهای قابل توجه شود.
برای تأیید مدل جدید بازیابی رطوبت خاک (معادله (30))، یک پایگاه داده شبیه سازی شده ثابت دی الکتریک تولید شد. مطالعات اخیر نشان داده است که پنج مدل متداول دی الکتریک خاک در شرایط مختلف برتری های متفاوتی دارند. مدل دابسون دقت خوبی از خود نشان می دهد و در بسیاری از الگوریتم های بازیابی رطوبت خاک استفاده شده است [ 31 ، 32 ]. از این رو مدل دابسون به عنوان مدل استاندارد در این پژوهش انتخاب شد. پایگاه داده شبیه سازی شده ثابت دی الکتریک حاوی بیش از 1.3 میلیون گروه از داده های بافت خاک، از جمله طیف وسیعی از رطوبت حجمی خاک، چگالی ظاهری، دمای خاک و درصد محتوای شن و خاک بود. پارامترهای دقیق مورد استفاده در پایگاه داده شبیه سازی شده در جدول 4 آمده است . سپس از این پایگاه داده برای واسنجی ضرایب مدل بازیابی رطوبت خاک در رابطه (29) استفاده شد. نتایج در جدول 5 فهرست شده است .
مقایسه بین رطوبت خاک برآورد شده با استفاده از مدل ثابت دی الکتریک جدید توسعه یافته و رطوبت ورودی خاک محاسبه شده با استفاده از مدل دابسون در شکل 3 نشان داده شده است . RMSE کلی رطوبت حجمی خاک 0.014 است که دقت بسیار خوبی را نشان می دهد. نتیجه نشان می دهد که رابطه بین ضریب شکست واقعی تعدیل شده وجود دارد Nr��و رطوبت حجمی خاک تطابق خوبی را نشان می دهد. ضریب شکست واقعی تنظیم شده Nr��اثرات مشابه ثابت دی الکتریک خاک را نشان می دهد ε�. هر دو می توانند خواص دی الکتریک خاک را منعکس کنند. بنابراین، ضریب شکست واقعی تعدیل شده می تواند جایگزین خواص دی الکتریک خاک در بازیابی رطوبت خاک شود. تنها دو ورودی پارامتر خاک و برخی از پارامترهای تجربی در مدل جدید توسعهیافته دخیل هستند که حجم کار و عدم قطعیتهای مرتبط با اندازهگیریهای درجا را کاهش میدهند. انحراف اندک در رطوبت خاک به دلیل ساده شدن جرم ظاهری و دمای خاک است. به طور کلی مدل ساده است و نسبت به شبیه سازی های مدل دابسون دارای خطای ناچیز است.
شکل 4 مقایسه ای بین رطوبت خاک برآورد شده با استفاده از مدل جدید توسعه یافته و رطوبت ورودی خاک محاسبه شده با استفاده از مدل AIEM را نشان می دهد. RMSE کلی در تمام زوایای بروز از 5 درجه تا 60 درجه کمتر از 0.03 است. دقت مدل جدید توسعهیافته در زوایای بزرگ بهتر است و همان تمایل مدل انتشار سطح قبلی را نشان میدهد. این نتایج به وضوح نشان میدهد که مدل جدید توسعهیافته دارای دقت بالایی در ارتباط با بازیابی رطوبت خاک، بدون استفاده از پارامترهای زبری سطح است.
4.2. تجزیه و تحلیل میزان حساسیت
با توجه به خطاهای اجتناب ناپذیر در اندازه گیری ها و مدل ها، تحلیل حساسیت برای ارزیابی عملکرد این مدل انجام شد. پارامترهای ورودی دمای روشنایی مایکروویو و دمای سطح از پایگاه داده شبیه سازی شده AIEM بودند. برای رادیومتر باند L، عدم قطعیت رادیومتری 0.8 تا 3 K برای SMOS [ 1 ] است، در حالی که SMAP تازه راه اندازی شده 1.3 K [ 2 ] است. حداکثر خطای طیفسنج تصویربرداری با وضوح متوسط (MODIS) در بازیابی دمای سطح بین 2 تا 3 کلوین است. 33]. بنابراین، سطوح مختلف نویز توزیع شده گاوسی به طور تصادفی به پایگاه داده شبیه سازی شده در این مطالعه اضافه شد. محدوده نویز اختصاص داده شده به دمای روشنایی و دمای سطح از 1 K تا 5 K است.
شکل 5 دقت مدل بازیابی رطوبت خاک را در سطوح مختلف نویز نشان می دهد. مدل جدید توسعه یافته به نویز در زوایای بروز متفاوت پاسخ متفاوتی می دهد. زوایای بزرگ نتایج بهتری را به همراه دارد و روند مشابهی را با آنچه در شکل 2 و شکل 3 نشان داده شده است نشان می دهد . شکل 5 a نشان می دهد که زوایای فرود کم، به ویژه آنهایی که در 5 درجه، 10 درجه و 15 درجه هستند، در هنگام ایجاد نویز عملکرد ضعیفی از خود نشان می دهند. برخی از این مقادیر بیش از 10 درصد هستند که در صورت عدم کیفیت داده های ورودی، عدم اطمینان زیادی ایجاد می کند. شکل 5b نشان می دهد که از آنجایی که زوایای فرود بزرگ از 20 درجه تا 60 درجه متفاوت است، تقریباً همه RMSE ها زیر 5٪ هستند. تاثیرگذارترین زاویه 45 درجه است که کمترین RMSE را در هر سطح نویز نشان می دهد. برای برآورده کردن آستانه دقت 0.04 متر مکعب بر متر مکعب ، که هدف محصولات رطوبت خاک SMOS و SMAP است، خطای اندازهگیری باید کمتر از 5 K باشد و زاویه برخورد انتخابی باید از 35 درجه تا 60 درجه باشد. .
4.3. اعتبارسنجی رطوبت خاک با داده های درجا
برای تأیید و آزمایش بیشتر مدل بازیابی رطوبت خاک توسعهیافته در این مطالعه، از مجموعه دادههای تجربی چهار ساله BARC از سال 1979 تا 1982 استفاده شد که شامل هر دو قطبش V و H بود. با توجه به توزیع ناهمگن دمای خاک و رطوبت خاک در جهت عمودی، داده های اعتبارسنجی بر اساس مقادیر میانگین کالیبره می شوند. دمای خاک بر اساس مقادیر متوسط در اعماق 0-2.5 سانتی متر و 2.5-5 سانتی متر انتخاب شد و رطوبت خاک با میانگین اندازه گیری رطوبت خاک از اعماق 2.5 سانتی متر و 7.5 سانتی متر به دست آمد که مطابق با اعماق نفوذ مربوط به باند L. شکل 6اندازهگیریهای حجمی رطوبت خاک را با اندازهگیریهای تخمین زده شده با استفاده از مدل جدید توسعهیافته مقایسه میکند. RMSE رطوبت حجمی خاک به ترتیب 7.0، 5.4، 4.3 و 3.4 درصد در زوایای 20 درجه، 30 درجه، 40 درجه و 50 درجه است. به طور خاص، دقت در 40 درجه و 50 درجه به طور کامل الزامات ماموریت های SMOS و SMAP را برآورده می کند، که نشان می دهد مدل جدید توسعه یافته قابل اعتماد است و می تواند به طور موثر برای بازیابی رطوبت خاک اعمال شود.
5. نتیجه گیری ها
رادیومتری باند L پتانسیل بالایی در پایش رطوبت خاک نشان می دهد. مدل سازی دقیق در باند L می تواند به طور قابل توجهی بازیابی رطوبت خاک را بهبود بخشد. در این مطالعه، ما یک الگوریتم جدید برای بازیابی رطوبت خاک باند L ایجاد کردیم که دو الگوریتم فرعی را ترکیب میکند. ابتدا، یک مدل انتشار سطح جدید برای ارزیابی اثرات زبری سطح پیشنهاد شد. این مدل اثرات زبری سطح را با استفاده از بازتاب سطح دو قطبی کاهش می دهد. ضرایب در این مدل با استفاده از پایگاه داده شبیه سازی شده AIEM کالیبره شده است که به طور دقیق در طیف گسترده ای از شرایط سطح ثابت شده است. روش جدید می تواند به طور مستقیم و بدون پارامترهای تجربی مورد استفاده قرار گیرد. سپس این مدل برای بازیابی رطوبت خاک سطحی با استفاده از مدل بازیابی رطوبت خاک اعمال شد. Nr��و رطوبت حجمی خاک ضریب شکست واقعی تنظیم شده Nr��که از نظر فیزیکی از معادله بازتاب فرنل استنتاج می شود، معنای فیزیکی مشابه دی الکتریک خاک را نشان می دهد. ε�. به طور کلی، الگوریتم ارائه شده در این مقاله در مقایسه با روش های قبلی دارای مزایای بسیاری است. این مدل از روش حل تحلیلی استفاده می کند که پایداری و کارایی بیشتری را نشان می دهد. اشتقاق دقیق بر اساس مدلهای فیزیکی است که الگوریتم را قابل اعتمادتر میکند. این الگوریتم به تعداد کمی از ورودی داده های خاک کمکی نیاز دارد که آن را در مقایسه با روش های قبلی کاربردی تر می کند. تنها پارامترهای ورودی الگوریتم پیشنهادی، دمای روشنایی مایکروویو با پلاریزاسیون دوگانه، دمای سطح، و محتویات ماسه و خاک رس است. به هیچ پارامتر زبری سطحی نیاز ندارد. در مقایسه با مقادیر رطوبت ورودی خاک شبیهسازی شده توسط مدل AIEM، رطوبت خاک برآورد شده با استفاده از این مدل دقت بسیار خوبی را با RMSEهای زیر 3 درصد در همه زوایای وقوع نشان داد. مدل پیشنهادی در زوایای بروز بالاتر، به ویژه در 40 درجه، 45 درجه و 50 درجه به خوبی عمل کرد. تجزیه و تحلیل حساسیت نشان داد که 45 درجه بهترین زاویه تابش برای این مدل است. اعتبار داده های اندازه گیری شده بر اساس آزمایش چهار ساله انجام شده در بلتسویل، مریلند بود. نتایج به ویژه برای زوایای فرود 20 درجه، 30 درجه، 40 درجه و 50 درجه خوب بود. خطا تقریباً 4٪ در زوایای فرود 40 درجه و 50 درجه بود که الزامات مأموریت SMOS و SMAP را برآورده می کند. نتایج نشان میدهد که مدل پیشنهادی در این مقاله پتانسیل عالی برای بازیابی رطوبت خاک را نشان میدهد و برای پرتوسنجهای باند L با زاویه برخورد بزرگ مانند SMAP مناسب است. اگرچه دقت کلی کافی بود، اما تمام دادههای تایید شده بر اساس سطوح برهنه بودند.
بدون نظر