GNSS Interferometrik Reflektometriyasi (Global Navigation Satellite System Interferometric Reflectometry, GNSS-IR) — sun’iy yo‘ldoshdan bevosita keladigan GNSS mikroto‘lqinli signali bilan tuproq yoki o‘simlik qoplamasidan qaytgandan so‘ng ayni qabul qilgichga yetib keladigan signal o‘rtasidagi interferensiyadan foydalanib, sirtning geometrik va dielektrik xususiyatlariga sezgir o‘lchovlarni hosil qiladigan passiv reflektometriya usulidir. Ushbu tadqiqot Italiyaning shimoli-g‘arbida joylashgan 2 gektarlik doimiy o‘tloqda to‘rtta dala kampaniyasida yig‘ilgan GPS/GLONASS SNR ma’lumotlarini qayta ishlab, GNSS-IRning qoplama balandligi, o‘simlik qoplami zichligi va namlik o‘zgarishlarini kuzatish salohiyatini o‘rgandi. Qoplama balandligini baholash yer usti o‘lchovlari bilan r = 0.89 korrelyatsiyani ko‘rsatdi, o‘rtacha xato −0.067 m va MAPE %18 deb topildi. dSNR tebranish amplitudasi bilan Sentinel-2 NDVI o‘rtasida r = −0.81, dielektrik faza hadi bilan dala tuproq namligi o‘rtasida r = 0.86 munosabat qayd etildi. Biroq tahlil faqat to‘rtta sananing birlashtirilgan qiymatlariga asoslanganligi sababli natijalar keng qishloq xo‘jaligi sharoitlariga umumlashtirilgan kalibrlash tenglamalari sifatida emas, balki GNSS-IR kuzatiladigan kattaliklari o‘simlik tuzilishi va dielektrik o‘zgarishlarga sezgirligini ko‘rsatuvchi dastlabki dala dalili sifatida talqin qilinishi kerak.

Tadqiqot yondashuvi butun dala xususiyatlarini bitta GNSS-IR parametridan chiqarishga urinmaydi. Tebranish chastotasi va undan hisoblangan faza markazi balandligi o‘simlikning vertikal tuzilishi bilan; tebranish amplitudasidagi susayish qoplama zichligi va elektromagnit so‘nish bilan; dielektrik faza hadi esa tuproq–o‘simlik muhitining kompleks dielektrik singdiruvchanligidagi, ayniqsa namlik sababli yuz beradigan o‘zgarishlar bilan bog‘langan. Bunday ko‘p parametrli talqin bir xil GNSS signalidan turli fizik tarkibiy qismlar haqida bir-birini to‘ldiruvchi ma’lumot olish mumkinligini ko‘rsatadi. Biroq tadqiqot strukturaviy va dielektrik hissalarni to‘liq bir-biridan ajratadigan inversiya tizimini ishlab chiqmagan.

GNSS-IRning amaliy jozibadorligi yangi faol radar signali yaratish o‘rniga uzluksiz uzatiladigan GNSS sun’iy yo‘ldoshlarining mavjud signallaridan foydalanishidadir. Bu usulni doimiy GNSS stansiyalariga, qishloq xo‘jaligi texnikasiga yoki kelajakda arzon qabul qilgichlarga integratsiya qilish imkonini berishi mumkin. Tadqiqotda ishlatilgan qabul qilgich esa ilmiy tadqiqot sinfidagi Topcon HiPer Pro bo‘lgan; smartfon yoki arzon qurilmalarning ishlashi ushbu tajribada sinovdan o‘tkazilmagan.

Qishloq xo‘jaligi monitoringi muammosi

O‘simlik balandligi, qoplama zichligi va tuproq namligi turli qishloq xo‘jaligi jarayonlarini tushuntiruvchi asosiy biofizik o‘zgaruvchilardir. Qoplama geometriyasi yorug‘likni tutib qolish, transpiratsiya va atmosfera–o‘simlik o‘zaro ta’sirlari bilan; tuproq namligi esa unib chiqish, oziqa moddalarning o‘zlashtirilishi, mikrobiologik jarayonlar va suv stressi bilan bog‘liq. Bu o‘zgaruvchilar bir dalaning o‘zida vaqt va makon bo‘yicha o‘zgarishi mumkinligi sababli siyrak nuqtaviy o‘lchovlar butun maydon holatini ifodalash uchun yetarli bo‘lmasligi mumkin.

Optik sun’iy yo‘ldosh kuzatuvlari katta maydonlarni muntazam tasvirlay olsa-da, bulut, aerozol va yoritish sharoitlaridan ta’sirlanadi hamda uch o‘lchamli qoplama tuzilishini bevosita o‘lchamaydi. LiDAR qoplama geometriyasini yuqori aniqlikda aniqlashi mumkin, ammo doimiy yuqori vaqt aniqligidagi sun’iy yo‘ldosh kuzatuvi sifatida keng tarqalmagan. Radar va passiv mikroto‘lqin tizimlari namlik hamda o‘simlik suv tarkibiga sezgir bo‘lsa-da, dala miqyosida vaqt/makon aniqligi bo‘yicha cheklovlarga ega bo‘lishi mumkin.

Tadqiqot ko‘rib chiqayotgan kuzatuv bo‘shlig‘i — o‘simlik tuzilishi va namlik tufayli yuz beradigan dielektrik o‘zgarishlarga bir vaqtning o‘zida sezgir, yuqori vaqt chastotasida ishlay oladigan dala miqyosidagi usulni o‘rganishdir.

GNSS-IR qishloq xo‘jaligi monitoringida qanday ishlaydi?

GNSS-IR GNSS sun’iy yo‘ldoshidan bevosita antennaga keladigan signal bilan tuproq va o‘simlik qoplamasidan qaytib, uzoqroq yo‘l bosib ayni antennaga yetib keladigan signalni birgalikda o‘lchaydi; ikki yo‘l orasidagi faza farqi SNRda muntazam interferensiya tebranishlarini hosil qiladi va bu tebranishlarning chastotasi, amplitudasi hamda fazasi mos ravishda sirt geometriyasi, signal susayishi va dielektrik xususiyatlar haqida ma’lumot berishi mumkin.

Bevosita va qaytgan signal interferensiyasi

Manba GNSS qabul qilgichi kuzatadigan signal-shovqin nisbatini Tenglama (1a) bilan ifodalaydi:

\[ \mathrm{SNR}(\theta)= \frac{|\widetilde{S}(\theta)|^2} {|\widetilde{N}(\theta)|^2} \tag{1a} \]

\(\theta\) — sun’iy yo‘ldoshning balandlik burchagi. \(\widetilde{S}\) qabul qilgichga yetib keladigan umumiy kompleks signalni, \(\widetilde{N}\) esa shovqin komponentini ifodalaydi.

Qabul qilingan umumiy signal quvvatida bevosita va qaytgan maydonlarning interferensiyasi manbada quyidagicha ochiladi:

\[ |\widetilde{S}(\theta)|^2 = |\widetilde{S}_d(\theta)+\widetilde{S}_r(\theta)|^2 = A_d^2(\theta)+A_r^2(\theta) + 2A_d(\theta)A_r(\theta)\cos(\Delta\phi) \tag{1b} \]

\(A_d\) bevosita signal amplitudasi, \(A_r\) qaytgan signal amplitudasidir. Oxirgi had tarkibidagi \(\cos(\Delta\phi)\) sun’iy yo‘ldosh harakati davomida SNR nima uchun davriy ravishda ko‘tarilib-pasayishini ko‘rsatadi.

Faza farqi manbada:

\[ \Delta\phi = \frac{2\pi}{\lambda}(R_1-R_2) = \frac{2\pi}{\lambda}\,2(h_a-h_v)\sin\theta \tag{1c} \]

ko‘rinishida beriladi. \(\lambda\) tashuvchi to‘lqin uzunligi, \(h_a\) antennaning yerdan balandligi, \(h_v\) esa o‘simlik qoplamasining samarali faza markazi balandligidir. \(R_1-R_2\) bevosita va qaytgan yo‘llar orasidagi geometrik yo‘l farqini anglatadi.

Tenglamaning asosiy natijasi shuki: antenna bilan samarali qaytaruvchi faza markazi orasidagi balandlik farqi ortgani sari sun’iy yo‘ldosh balandlik burchagiga bog‘liq faza o‘zgarishi ham o‘zgaradi. Shu sababli SNR tebranish chastotasi geometrik balandlik ma’lumotiga aylantirilishi mumkin.

Nega faqat 5°–30° balandlik burchaklari ishlatildi?

Sun’iy yo‘ldosh osmonda ko‘tarilgani sari bevosita va qaytgan signal yo‘llari orasidagi geometrik farq kamayadi. Bunda multipath interferensiya tebranishlarining amplitudasi va ajratib olish imkoniyati pasayadi. Manba 5°–30° oralig‘ida interferensiya izi aniqroq bo‘lishini ko‘rsatgan va shu oraliqdan tashqaridagi SNR qiymatlarini niqoblagan.

Past balandlik burchaklarida qaytgan signalning qabul qilgichga yetib kelishiga hissa qo‘shadigan sirt maydoni birinchi Fresnel zonasi bilan taxminan ifodalanadi. Manba yarim o‘qlarni va qaytish nuqtasining antennaga masofasini:

\[ a= \frac{\sqrt{\lambda h_a\sin\theta}} {(\sin\theta)^2} \tag{2a} \]

\[ b= \frac{\sqrt{\lambda h_a\sin\theta}} {\sin\theta} \tag{2b} \]

\[ d= \frac{h_a}{\tan\theta} \tag{2c} \]

munosabatlari bilan beradi. \(a\) va \(b\) Fresnel ellipsining taxminiy yarim o‘qlarini, \(d\) esa qaytish nuqtasining antennaga gorizontal masofasini ifodalaydi.

Tadqiqotdagi barcha mos sun’iy yo‘ldosh izlari birgalikda baholanganda, GNSS-IR har bir dala kampaniyasida taxminan 5000 m² dala maydonidan qaytish ma’lumotini olgani hisoblangan. Bu footprint nuqtaviy tuproq zondi yoki bitta lineyka o‘lchovidan ancha katta maydonni ifodalaydi; biroq maydon ichidagi geterogenlik ham o‘lchov noaniqligining bir qismiga aylanadi.

Dala va fenologik kontekst

Tajriba maydoni Italiyaning shimoli-g‘arbida, Piedmont hududidagi Pogliola joyida joylashgan 2 gektarlik doimiy ko‘p turlik o‘tloqdir. O‘simlik qoplami asosan Gramineae va Leguminosae turlaridan tashkil topgan. 2025 yilda dala 20 may, 8 iyul va 13 avgust kunlari o‘rilgan, 15 sentyabrda esa boqilgan.

Sentinel-2 nima uchun ishlatildi?

Sentinel-2 ma’lumoti GNSS-IR modelining kirishi emas. Tadqiqot GNSS-IR tebranish amplitudasi mavsumiy qoplama zichligi o‘zgarishiga sezgir yoki sezgir emasligini mustaqil kuzatuv bilan baholash uchun NDVIdan foydalangan.

1 aprel–1 oktyabr davridagi jami 90 Sentinel-2 Level-2A Harmonized tasviri Google Earth Engineda qayta ishlangan, bulut/cirrus/soya piksellari Scene Classification Layer yordamida chiqarib tashlangan. Har bir yaroqli sanada 10 m NIR va qizil diapazonlardan NDVI hisoblangan hamda to‘rtta GNSS kampaniyasi uchun maydon mediani olingan.

NDVI ushbu tadqiqotda biomassa yoki Leaf Area Indexning bevosita miqdoriy o‘lchovi emas. Manba NDVIni faqat yashil o‘simlik qoplami va qoplama zichligidagi mavsumiy o‘zgarishlarni ifodalovchi nisbiy proxy sifatida talqin qilgan.

Yer usti etalonlari

Tuproq namligi uchun Drill & Drop Bluetooth sig‘imli sensori 10 cm chuqurlikka ko‘milgan va hajmiy suv tarkibi har soatda qayd etilgan. GNSS kampaniya kunlarida 08:00–14:00 oralig‘idagi SMC qiymatlari o‘rtachasi olingan.

Qoplama balandligi \(h_v\) esa GNSS antennasi atrofidagi yetti nuqtada yerga tik tutilgan lineyka bilan o‘lchangan. Yer o‘lchovlarining o‘rtacha qiymati va standart og‘ishi hisoblangan.

SNRdan dSNRga

Xom SNR vaqt qatori faqat qaytgan signalni o‘z ichiga olmaydi. Antennaning gain patterni va atmosfera ta’sirlari sababli bevosita signal sun’iy yo‘ldosh ko‘tarilgan sari silliq trend hosil qiladi. Manba ushbu past chastotali deterministik trendni \(\sin\theta\) mustaqil o‘zgaruvchisiga nisbatan ikkinchi darajali model bilan olib tashlab, o‘rtachasi nolga yaqin trenddan tozalangan SNR (dSNR) qatorini hosil qilgan.

Shundan keyin dSNR ichidagi muntazam tebranishlar asosan sirtdan qaytgan koherent komponentning geometrik va elektromagnit ta’sirlarini tashiydi.

GNSS-IR o‘simlik bo‘yini qanday baholaydi?

GNSS-IR o‘simlik bo‘yini bevosita ko‘rish orqali o‘lchamaydi; dSNR tebranishlarining \(\sin\theta\) bo‘yicha dominant chastotasini Lomb–Scargle periodogrami orqali aniqlaydi, ushbu chastotani antenna bilan samarali qaytaruvchi faza markazi orasidagi reflektor balandligiga aylantiradi va reflektor balandligini antenna balandligidan ayirib qoplama faza markazi balandligi \(h_v\)ni oladi.

Lomb–Scargle periodogrami

Manba dSNR qatorining \(\sin\theta\) bo‘yicha garmonik tarkibini Lomb–Scargle periodogrami bilan tahlil qiladi. Ishlatilgan asosiy matritsa ko‘rinishi:

\[ P(\omega)=0.5\,\beta(\omega)^{T}\beta(\omega) \tag{3a} \]

\[ \beta(\omega) = \left[X(\omega)^TX(\omega)\right]^{-1} X(\omega)^Ty \tag{3b} \]

va:

\[ X(\omega)= \begin{bmatrix} \cos(\omega\sin\theta_1) & \sin(\omega\sin\theta_1)\\ \vdots & \vdots\\ \cos(\omega\sin\theta_n) & \sin(\omega\sin\theta_n) \end{bmatrix} \tag{3c} \]

ko‘rinishidadir. \(y\) dSNR kuzatuvlarini, \(n\) vaqt qatoridagi kuzatuvlar sonini, \(\omega=2\pi f\) esa burchak spektral o‘zgaruvchisini ifodalaydi.

Bu yerdagi \(f\) vaqt bo‘yicha Hertz birligidagi chastota emas. Manba uni \(\sin\theta\) birligiga to‘g‘ri keladigan tebranish sikllari deb ta’riflaydi. Qidiruv oralig‘i 0.1–50 cycles per unit of \(\sin\theta\), qadam esa 0.01.

Dominant chastotadan qoplama balandligiga

Manba Tenglama (4a)ni bosma ko‘rinishda quyidagicha beradi:

\[ f= \frac{\partial\Delta\phi}{\partial t} \frac{1}{2\pi} \simeq \frac{2RH}{\lambda} \tag{4a} \]

Manba notatsiyasi bo‘yicha eslatma: Shu bo‘limning o‘zida \(f\) \(\sin\theta\) birligiga to‘g‘ri keladigan tebranish chastotasi sifatida ta’riflansa-da, Tenglama (4a)dagi hosila o‘zgaruvchisi manba PDFida \(t\) sifatida bosilgan. Verianla ushbu ifodani manbada yozilmagan boshqa o‘zgaruvchi bilan jim tarzda almashtirmaydi.

Chastota olingach, reflektor balandligi:

\[ RH=\frac{f\lambda}{2} \tag{4b} \]

va o‘simlik faza markazi balandligi:

\[ h_v=h_a-RH \tag{4c} \]

tarzida hisoblanadi.

\(h_a\) antennaning yerdan balandligi bo‘lgan 2 mdir. \(RH\) antenna bilan samarali qaytaruvchi faza markazi orasidagi masofani, \(h_v\) esa ushbu faza markazining yer ustidagi balandligini ifodalaydi.

Sifat nazorati

Filtrlash natijasida dastlab har bir kampaniyada taxminan 100 iz bo‘lsa-da, umumiy sun’iy yo‘ldosh izlarining %8dan kamrog‘i yakuniy modellashtirishda qolgan. Bu qolgan ma’lumotlarning sifati va fizik talqin qilinishini oshiradi, biroq murakkab sochilish sharoitlarini ifodalash nuqtai nazaridan tanlanish ta’sirini keltirib chiqarishi mumkin.

Nolinear dSNR modeli

Dominant geometrik chastota va \(h_v\) aniqlangach, manba tebranish amplitudasi bilan dielektrik fazani ajratish uchun ikkinchi bosqichda nolinear garmonik modeldan foydalangan:

\[ dSNR(\theta) = A(\theta) \cos \left[ \phi_{geom}(\theta)+\phi_{die}(\theta) \right] \tag{5a} \]

Modelning qo‘llangan ko‘rinishi:

\[ dSNR(\theta) = A_0e^{m_0\sin\theta} \cos \left[ \frac{4\pi}{\lambda}(h_a-h_v)\sin\theta+\phi_0 \right] \tag{5b} \]

ko‘rinishidadir.

\(A_0\) asosiy amplituda masshtabi, \(m_0\) balandlik burchagiga bog‘liq amplituda so‘nishini, \(\phi_0\) esa qaytaruvchi muhitning dielektrik xatti-harakati bilan bog‘liq taxminan doimiy faza siljishini ifodalaydi. \(h_v\) bu bosqichda erkin ravishda qayta optimallashtirilmagan; oldingi Lomb–Scargle chastota tahlilidan olingan geometrik ma’lumot sobit saqlangan.

Nolinear moslashtirish R muhitida `minpack.lm` paketidagi `nlsLM` funksiyasi bilan, bounded least-squares Levenberg–Marquardt algoritmi yordamida bajarilgan. Boshlang‘ich amplituda kuzatilgan maksimal dSNR qiymatiga, \(m_0=-0.1\)ga, \(\phi_0=0\)ga o‘rnatilgan. Chegaralar \(A>0\), \(m\le0\) va \(\phi\in[-\pi,\pi]\) ko‘rinishidadir.

Qaytgan maydonning fizik modeli

Manba koherent va nokoherent qaytishni:

\[ \widetilde{S}_r(\theta) = |\Gamma(\theta,\widetilde{\varepsilon})| e^{j\,\angle[\Gamma(\theta,\widetilde{\varepsilon})]} + \widetilde{S}_{incoh} \tag{5c} \]

bilan tushuntiradi.

\(\Gamma\) tuproq–o‘simlik muhitining samarali kompleks dielektrik singdiruvchanligi \(\widetilde{\varepsilon}\)ga bog‘liq Fresnel koeffitsiyentidir. Tasodifiy yo‘nalgan barg, poya, bo‘shliq va namlik geterogenliklari hosil qiladigan \(\widetilde{S}_{incoh}\) komponenti deterministik tebranish hosil qilish o‘rniga shovqin bazasini oshiradi.

Shu sababli dSNRning mazmunli garmonik komponenti qaytgan maydonning barchasi emas, balki yetarli darajada koherent bo‘lib qolgan qismi bevosita signal bilan interferensiyaga kirishishining kuzatiladigan natijasidir.

Nega faza markazi haqiqiy o‘simlik cho‘qqisidan past bo‘lishi mumkin?

O‘tloq yoki dala qoplami mikroto‘lqin nuqtai nazaridan yagona qattiq sirt emas. Barg va poyalar turli balandliklarda sochilish hosil qiladi. Manba ushbu taqsimlangan tuzilishni:

\[ \widetilde{S}_r(\theta) = \int_{z_0}^{z_0+h_{top}} F_v(z,\theta)e^{ikz}\,dz \tag{6} \]

bilan ifodalaydi.

\(F_v(z,\theta)\) qoplamaning vertikal reflektivlik profilidir; \(z_0\) yer sathi balandligi, \(h_{top}\) fizik qoplama balandligi va \(k=4\pi/\lambda\) to‘lqin sonidir.

Ushbu taqsimotning faza markazi:

\[ PCH(\theta)= \frac{ \int_{z_0}^{z_0+h_{top}} z\,F_v(z,\theta)\,dz }{ \int_{z_0}^{z_0+h_{top}} F_v(z,\theta)\,dz } \tag{7} \]

tarzida ta’riflanadi.

Tenglama (7) GNSS-IR nima uchun lineyka bilan o‘lchangan eng yuqori barg sathini majburiy ravishda bermasligini tushuntiradi. Chuqurroq o‘simlik qatlamlarida so‘nish va sochilish ortgani sababli koherent qaytishning og‘irlik markazi fizik qoplama cho‘qqisidan pastda joylashishi mumkin. Tadqiqotda kuzatilgan manfiy balandlik biasi ana shu faza-markazi xatti-harakati bilan bog‘langan.

Tadqiqot Usuli va Natijalari

Ma’lumot yig‘ish va tekshirish zanjiri

Qoplama balandligi natijasi

Balandlik baholari, ayniqsa, iyun–iyul davrida yuqoriroq, avgust–sentyabrda esa o‘rimdan keyingi pastroq qoplama sharoitlarini kuzatgan. Manba iyun–iyul davrida o‘simlik balandligi taxminan 50 cm atrofida bo‘lganini, keyingi o‘rim davrlari esa balandlik va NDVI qiymatlarini pasaytirganini muhokama qiladi.

Qoplama zichligi va dSNR amplitudasi

Sentinel-2 NDVI bilan θ = 15° uchun hisoblangan dSNR model amplitudasi o‘rtasida:

r = −0.81

munosabat aniqlangan.

Manfiy ishora o‘simlik qoplami zichlashgan sari dSNR tebranish amplitudasi kamayishini ko‘rsatadi. Manbaning fizik talqiniga ko‘ra, barg maydoni, biomassa va suv tarkibining ortishi L-band koherent qaytishini susaytiradi.

Ushbu natija NDVI GNSS-IR tomonidan baholanganini anglatmaydi. NDVI faqat GNSS-IR amplituda parametri o‘simlik qoplami rivojlanishiga sezgir yoki sezgir emasligini baholovchi mustaqil proxydir.

Tuproq namligi va dielektrik faza

Yer o‘lchovi soil moisture content (SMC) bilan \(|\phi_{die}|\) o‘rtasida:

r = 0.86

musbat munosabat aniqlangan.

Tuproqdagi suv miqdori ortishi tuproq–o‘simlik muhitining kompleks dielektrik singdiruvchanligini o‘zgartirgani uchun Fresnel qaytish koeffitsiyentining fazasi ham o‘zgarishi mumkin. Shu sababli manba \(\phi_0\) parametrini namlikdan kelib chiqadigan dielektrik o‘zgarishlarga sezgir kuzatiladigan kattalik sifatida talqin qiladi.

GNSS-IR ushbu tadqiqotda o‘simlik tuzilishi bilan tuproq namligini to‘liq bir-biridan ajrata oladimi?

Yo‘q. Tadqiqot qoplama balandligi, dSNR amplituda so‘nishi va dielektrik faza mos ravishda tuzilish, qoplama zichligi va namlik o‘zgarishlariga turlicha sezgirlik ko‘rsatishini qo‘llab-quvvatlaydi; biroq qaytgan GNSS signali tuproq va o‘simlik qoplami birgalikda hosil qiladigan elektromagnit muhitdan kelganligi sababli strukturaviy va dielektrik hissalar ushbu ma’lumotlar to‘plamida to‘liq va aniq ajratilmagan.

Uchta GNSS-IR kuzatiladigan kattaligining bir-birini to‘ldiruvchi ma’nosi

Ushbu jadval uchta parametr mutlaqo mustaqil fizik o‘lchov ekanini anglatmaydi. Bir xil qaytish jarayonida geometriya, qoplama zichligi, suv miqdori va sochilish tuzilishi bir-biriga ta’sir qilishi mumkin. Manbaning xulosasi turli GNSS-IR parametrlari turli fizik komponentlarga ustun sezgirlik ko‘rsatishidir; to‘liq separation hali ta’minlanmagan.

Tuproq namligi natijalaridagi usul–natija nomuvofiqligi

Usul bo‘limi NDVI va SMCni bevosita retrieval maqsadi sifatida belgilamaydi va ushbu ikki etalon bilan munosabat faqat Pearson korrelyatsiyasi orqali baholanishini bildiradi. Shunga qaramay, natijalar bo‘limida SMC uchun qo‘shimcha ravishda “ME = −5.3” va “MAPE = %21” yozilgan.

Manba \(\phi_{die}\) bilan SMC o‘rtasida qanday qilib bevosita SMC bahosi hosil qilinganini va ushbu ikki xato metrikasi qaysi hisoblash zanjiridan kelib chiqqanini usul bo‘limida tushuntirmaydi. Shu sababli Verianla bu ikki raqamni tasdiqlangan “GNSS-IR soil-moisture retrieval error” metrikasi sifatida ishlatmaydi. Ishonchli ravishda qo‘llab-quvvatlangan miqdoriy munosabat Pearson \(r=0.86\)dir.

Korrelyatsiyalarning statistik chegarasi

Barcha asosiy korrelyatsiyalar faqat to‘rtta dala kampaniyasining birlashtirilgan qiymatlariga asoslanadi: \(n=4\). Shu sababli manba korrelyatsiyalar qat’iy statistik umumlashtirish sifatida emas, fizik munosabatlarning ko‘rsatkichi sifatida talqin qilinishi kerakligini aniq ta’kidlaydi.

To‘rtta sana ayni sharoitlarning tasodifiy tanlangan takrorlari emas; o‘rimdan oldingi va o‘rimdan keyingi turli fenologik rejimlarni ifodalaydi. Bu xilma-xillik usulning turli dala holatlariga javob berishini ko‘rsatishga yordam beradi, ammo namuna sonini oshirmaydi.

Filtrlash va selection bias

Sifat filtrlari SNR ichida fizik jihatdan talqin qilish mumkin bo‘lgan koherent multipath tebranishi mavjud izlarni tanlaydi. Bu yondashuv shovqinli va beqaror yechimlarni kamaytiradi.

Buning boshqa tomoni shuki, geterogen qoplama, murakkab qaytish geometriyasi yoki tezroq o‘zgaruvchi namlik sharoitlarida model farazlariga kamroq mos keladigan sun’iy yo‘ldosh izlari ma’lumotlar to‘plamidan ko‘proq chiqarib tashlanishi mumkin. Mualliflar shu sababli yakuniy izlar barcha dala sharoitlarining xolis namunasi emas, balki model eng aniq ishlaydigan kuzatuvlar ekanini qayd etadi.

Lomb–Scargle balandligidagi manfiy bias

Lomb–Scargle periodogrami teng oraliqli bo‘lmagan SNR namunalarida dominant tebranish chastotasini topish uchun kuchli vositadir; biroq chekli sun’iy yo‘ldosh burchagi oynasi spectral leakage hosil qilishi mumkin. Qo‘shimcha qaytish yo‘llari ham asosiy periodogram cho‘qqisiga ta’sir qilishi mumkin.

Manba bu ta’sirlar dominant chastotani biroz pastga siljitishi va bu reflektor balandligi bahosida bias sifatida aks etishi mumkinligini muhokama qiladi. Bundan tashqari, asosiy fizik sabablardan biri qoplamaning hajmiy sochuvchi muhit ekanidir: faza markazi fizik yuqori sathdan pastda bo‘ladi.

Tadqiqot qo‘llab-quvvatlaydigan natijalar

Tadqiqot to‘rtta turli fenologik dala sharoitida GNSS-IR dSNR signali qoplama tuzilishi va namlik bilan bog‘liq tizimli o‘zgarishlarni tashishini qo‘llab-quvvatlaydi.

Tadqiqot Lomb–Scargle asosidagi GNSS-IR reflektor balandligi o‘tloq qoplamasidagi balandlik o‘zgarishini kuzata olishini va lineyka etaloni bilan kuchli munosabat ko‘rsatishini qo‘llab-quvvatlaydi.

Tadqiqot dSNR amplituda so‘nishi mavsumiy qoplama zichligi/NDVI o‘zgarishiga, dielektrik faza parametri esa dala namligi o‘zgarishlariga sezgir ekanini qo‘llab-quvvatlaydi.

Tadqiqot GNSS-IR bitta kuzatuv infratuzilmasidan bir nechta qishloq xo‘jaligi biofizik ma’lumot kanallarini chiqarish salohiyatiga ega ekanini qo‘llab-quvvatlaydi.

Tadqiqot qo‘llab-quvvatlamaydigan natijalar

Tadqiqot GNSS-IR barcha ekin turlarida %18 qoplama balandligi xatosi bilan ishlashini ko‘rsatmaydi.

Tadqiqot \(\phi_{die}\) o‘simlik qoplami ta’siridan to‘liq tozalangan sof tuproq namligi o‘lchovi ekanini ko‘rsatmaydi.

Tadqiqot dala sug‘orishi avtomatik ravishda GNSS-IR bilan boshqarilganini yoki suv iste’moli kamaytirilganini ko‘rsatmaydi.

Tadqiqot arzon smartfon GNSS tizimlari ayni ishlash natijasini berishini ko‘rsatmaydi; smartfonlardan foydalanish kelajakdagi tadqiqot yo‘nalishidir.

Tadqiqot makkajo‘xori, bug‘doy, soya yoki boshqa dala ekinlarida ayni munosabatlar tasdiqlanganini ko‘rsatmaydi.

Asosiy cheklovlar

Kelajakdagi tadqiqot yo‘nalishlari

Manba vaqt bo‘yicha tez-tez namuna olish, ko‘proq maydon, turli ekin arxitekturalari va geterogen sharoitlarda tekshirishni tavsiya qiladi. Xususan maize, wheat va soybean kelajakda sinash rejalashtirilgan ekinlar qatorida qayd etilgan.

Yanada moslashuvchan multi-frequency yoki time-varying harmonic modellar, robust regression yondashuvlari va qoplama geterogenligini bevosita hisobga oladigan fizik scattering modellar ko‘proq GNSS izini yo‘qotmasdan parametr olish uchun o‘rganilishi mumkin.

Arzon ikki chastotali GNSS qabul qilgichlari va smartfon sensorlari ham kelajakdagi kengaytiriladigan qo‘llash yo‘llari qatorida ko‘rsatilgan. Bu salohiyat joriy tadqiqotning sinovdan o‘tgan natijasi emas, mualliflarning kelajakdagi rivojlantirish yo‘nalishidir.

Manba va Usul Haqida Izoh

Asl sarlavha: Monitoring Crop Structure and Moisture Using GNSS Interferometric Reflectometry Based on SNR Modeling

Mualliflar: Samuele De Petris; Enrico Borgogno-Mondino

Yozishmalar muallifi: Enrico Borgogno-Mondino

Muassasa: Department of Agriculture, Forestry and Food Science, University of Turin, Grugliasco, Italy.

Jurnal: Agronomy

Nashriyot: MDPI

Jild / Son / Maqola: 16 / 9 / 922

DOI: 10.3390/agronomy16090922

Sanalar: Qabul qilindi: 25 mart 2026; qayta ko‘rib chiqildi: 27 aprel 2026; ma’qullandi: 29 aprel 2026; nashr qilindi: 1 may 2026.

Manba turi: Taqrizdan o‘tgan eksperimental tadqiqot maqolasi; dala GNSS o‘lchovi, mikroto‘lqinli interferometrik reflektometriya, spektral tahlil va nolinear signal modellashtirish.

Versiya: Version of Record.

Litsenziya: Creative Commons Attribution (CC BY). PDFning ko‘rinadigan litsenziya matnida versiya raqami alohida ko‘rsatilmagani uchun litsenziya versiyasi taxmin qilinmagan.

Asosiy maydon: Pogliola, Piedmont, Italiyaning shimoli-g‘arbi; 2 ha doimiy ko‘p turlik o‘tloq.

GNSS qabul qilgichi: Topcon HiPer Pro; GPS va GLONASS; L1/L2; antenna balandligi 2 m; 1 Hz.

GNSS kampaniyalari: 1 iyun, 4 iyul, 2 avgust va 7 sentyabr 2025.

Optik yordamchi ma’lumot: Copernicus Sentinel-2 Level-2A Harmonized; 90 tasvir; 1 aprel–1 oktyabr 2025; Google Earth Engine.

Tuproq namligi etaloni: Drill & Drop Bluetooth sig‘imli sensori; 10 cm chuqurlik; soatlik volumetric water content.

Qoplama balandligi etaloni: GNSS antennasi atrofidagi yetti nuqtada qo‘lda lineyka o‘lchovi.

Asosiy signalga ishlov berish: Ascending/descending sun’iy yo‘ldosh yoyi ajratish; 5°–30° elevation mask; quadratic detrending; Lomb–Scargle periodogrami; nonlinear harmonic regression.

LSP chastota oralig‘i: 0.1–50 cycles per unit of sin θ; qadam = 0.01.

Sifat filtrlari: P2N > 3 dB-Hz; Δθ > 15°; 0.05 m < h_v < 2 m; nonlinear model R² > 0.6.

Nolinear fit: R `nlsLM`, `minpack.lm`; bounded least-squares Levenberg–Marquardt.

Asosiy natijalar: Qoplama balandligi r = 0.89, ME = −0.067 m, MAPE = %18; NDVI bilan dSNR amplitudasi r = −0.81; SMC bilan mutlaq dielektrik faza r = 0.86.

Statistik namuna cheklovi: Asosiy taqqoslashlar to‘rtta kampaniyaning birlashtirilgan qiymatlariga asoslanadi, n = 4.

Moliyalashtirish: Tadqiqot tashqi moliyalashtirish olmagan.

Ma’lumotlardan foydalanish: Ma’lumotlar so‘rov bo‘yicha taqdim etilishi bildirilgan.

Qo‘shimcha material: Table S1 har bir sifat nazorati bosqichidan keyin qolgan sun’iy yo‘ldosh izlari sonini o‘z ichiga oladi.

Manfaatlar to‘qnashuvi: Mualliflar manfaatlar to‘qnashuvi yo‘qligini bildirgan.

Mualliflar hissasi: Samuele De Petris konseptuallashtirish, usul, dasturiy ta’minot, tekshirish, formal tahlil, ma’lumotlar kuratsiyasi, dastlabki matn va vizualizatsiya rollarini bajargan. Tadqiqot va resurslar ikki muallif o‘rtasida taqsimlangan. Enrico Borgogno-Mondino maslahat, loyiha boshqaruvi va moliyalashtirishni jalb qilish rollarini bajargan; har ikki muallif ko‘rib chiqish va tahrirlashga hissa qo‘shgan.

Manba ichidagi notatsiya ogohlantirishi: Matn Lomb–Scargle spektral o‘zgaruvchisini sin θ birligiga to‘g‘ri keladigan tebranish sikli sifatida ta’riflagan bo‘lsa-da, Tenglama (4a)da hosila o‘zgaruvchisi t sifatida bosilgan. Bu ifoda Verianla tomonidan jim tarzda o‘zgartirilmagan.

Manba ichidagi usul/natija ogohlantirishi: Usul bo‘limi soil moisture va NDVI taqqoslashlarida faqat Pearson r ishlatilishini bildirsa-da, natijalar bo‘limida soil moisture uchun qo‘shimcha ravishda ME = −5.3 va MAPE = %21 berilgan. Ushbu xato metrikalarining qanday hosil qilingani tushuntirilmagani sababli mustaqil namlik retrieval ishlashi sifatida foydalanilmagan.

Asosiy metodologik cheklov: Bitta maydon, bitta qoplama turi, to‘rtta kampaniya, siyrak yer etalonlari va qat’iy sifat filtrlari sababli natijalar amaliy ko‘p ekinli qishloq xo‘jaligi tizimi uchun yakuniy kalibrlash emas.

Asosiy ilmiy talqin chegarasi: Tadqiqot GNSS-IR kuzatiladigan kattaliklari o‘simlik tuzilishi va dielektrik/namlik o‘zgarishlariga turlicha sezgirlik ko‘rsatishini qo‘llab-quvvatlaydi; bu ikki fizik hissa to‘liq bir-biridan ajratilganini ko‘rsatmaydi.