GNSS интерферометриялык рефлектометриясы (Global Navigation Satellite System Interferometric Reflectometry, GNSS-IR) — спутниктен түз келген GNSS микротолкундуу сигналы менен топурактан же өсүмдүк катмарынан чагылып, ошол эле кабылдагычка жеткен сигналдын ортосундагы интерференцияны колдонуп, беттин геометриялык жана диэлектрдик касиеттерине сезимтал өлчөөлөрдү түзгөн пассивдүү рефлектометрия ыкмасы. Бул изилдөө Италиянын түндүк-батышындагы 2 гектар туруктуу шалбаада төрт талаа кампаниясында чогултулган GPS/GLONASS SNR маалыматтарын иштетип, GNSS-IRдин өсүмдүк жамынтыгынын бийиктигин, өсүмдүк катмарынын тыгыздыгын жана нымдуулук өзгөрүүлөрүн көзөмөлдөө мүмкүнчүлүгүн изилдеген. Өсүмдүк жамынтыгынын бийиктигин баалоо жердеги өлчөөлөр менен r = 0.89 корреляция көрсөткөн, орточо ката −0.067 m жана MAPE %18 болгон. dSNR термелүү амплитудасы менен Sentinel-2 NDVI ортосунда r = −0.81, диэлектрдик фаза мүчөсү менен талаадагы топурак нымдуулугунун ортосунда r = 0.86 байланыш билдирилген. Бирок талдоо болгону төрт датанын агрегатталган маанилерине негизделгендиктен, натыйжаларды кеңири айыл чарба шарттарына жалпыланган калибрлөө теңдемелери катары эмес, GNSS-IR байкалуучу чоңдуктары өсүмдүк түзүлүшүнө жана диэлектрдик өзгөрүүлөргө сезимтал экенин көрсөткөн алгачкы талаа далили катары чечмелөө керек.

Изилдөөнүн мамилеси бардык талаа касиеттерин бир GNSS-IR параметринен чыгарууга аракет кылбайт. Термелүү жыштыгы жана андан эсептелген фаза борборунун бийиктиги өсүмдүктүн тик түзүлүшү менен; термелүү амплитудасынын басаңдашы өсүмдүк катмарынын тыгыздыгы жана электромагниттик басаңдоо менен; диэлектрдик фаза мүчөсү болсо топурак–өсүмдүк чөйрөсүнүн комплекстүү диэлектрдик өтүмдүүлүгүндөгү, өзгөчө нымдуулуктан келип чыккан өзгөрүүлөр менен байланышкан. Бул көп параметрлүү түшүндүрмө бир эле GNSS сигналынан ар башка физикалык компоненттер жөнүндө бири-бирин толуктаган маалымат алууга болорун көрсөтөт. Бирок изилдөө структуралык жана диэлектрдик салымдарды толук ажыраткан инверсия системасын иштеп чыккан эмес.

GNSS-IRдин практикалык артыкчылыгы жаңы активдүү радар сигналын түзүүнүн ордуна тынымсыз берилүүчү GNSS спутниктеринин бар сигналдарын колдонушунда. Бул ыкманы туруктуу GNSS станцияларына, айыл чарба машиналарына же келечекте арзан кабылдагычтарга интеграциялоого мүмкүнчүлүк бериши ыктымал. Изилдөөдө колдонулган кабылдагыч болсо илимий класстагы Topcon HiPer Pro болгон; смартфондордун же арзан түзмөктөрдүн иштеши бул экспериментте текшерилген эмес.

Айыл чарба мониторингинин көйгөйү

Өсүмдүктүн бийиктиги, өсүмдүк катмарынын тыгыздыгы жана топурак нымдуулугу ар кандай айыл чарба процесстерин түшүндүргөн негизги биофизикалык өзгөрмөлөр. Өсүмдүк катмарынын геометриясы жарыкты кармоо, транспирация жана атмосфера–өсүмдүк өз ара аракеттери менен; топурак нымдуулугу болсо өнүп чыгуу, азык заттарды сиңирүү, микробдук процесстер жана суу стресси менен байланыштуу. Бул өзгөрмөлөр бир эле талаада убакыт жана мейкиндик боюнча өзгөрүшү мүмкүн болгондуктан, сейрек чекиттик өлчөөлөр бүт аянттын абалын жетиштүү чагылдырбай калышы мүмкүн.

Оптикалык спутниктик байкоолор чоң аймактарды үзгүлтүксүз сүрөткө тарта алганы менен, булут, аэрозоль жана жарыктандыруу шарттарынан таасирленет жана өсүмдүк катмарынын үч өлчөмдүү түзүлүшүн түздөн-түз өлчөбөйт. LiDAR өсүмдүк катмарынын геометриясын жогорку тактыкта аныктай алат, бирок үзгүлтүксүз жогорку убакыттык чечилишке ээ спутниктик байкоо катары кеңири колдонулбайт. Радар жана пассивдүү микротолкун системалары нымдуулукка жана өсүмдүктүн суу курамына сезимтал болгону менен талаа масштабында убакыт/мейкиндик чечилиши боюнча чектөөлөр болушу мүмкүн.

Изилдөө караган байкоо боштугу — өсүмдүк түзүлүшү менен нымдуулуктан келип чыккан диэлектрдик өзгөрүүлөргө бир убакта сезимтал жана жогорку убакыттык жыштыкта иштей ала турган талаа масштабындагы ыкманы изилдөө.

GNSS-IR айыл чарба мониторингинде кантип иштейт?

GNSS-IR GNSS спутнигинен түз антеннага келген сигнал менен топурак жана өсүмдүк катмарынан чагылып, узагыраак жол өтүп ошол эле антеннага жеткен сигналды бирге өлчөйт; эки жолдун ортосундагы фаза айырмасы SNRде үзгүлтүксүз интерференциялык термелүүлөрдү жаратат жана бул термелүүлөрдүн жыштыгы, амплитудасы жана фазасы тиешелүүлүгүнө жараша беттин геометриясы, сигналдын басаңдашы жана диэлектрдик касиеттер жөнүндө маалымат бере алат.

Түз жана чагылган сигналдардын интерференциясы

Булак GNSS кабылдагычы байкаган сигнал-ызы-чуу катышын Теңдеме (1a) менен берет:

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

\(\theta\) спутниктин көтөрүлүү бурчу. \(\widetilde{S}\) кабылдагычка жеткен жалпы комплекстүү сигналды, \(\widetilde{N}\) ызы-чуу компонентин билдирет.

Жалпы кабыл алынган сигнал кубатында түз жана чагылган талаалардын интерференциясы булакта төмөнкүдөй ачылат:

\[ |\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\) түз сигналдын амплитудасы, \(A_r\) чагылган сигналдын амплитудасы. Акыркы мүчөдөгү \(\cos(\Delta\phi)\) спутник кыймылдаганда SNR эмне үчүн мезгил-мезгили менен көтөрүлүп-түшөрүн көрсөтөт.

Фаза айырмасы булакта:

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

түрүндө берилет. \(\lambda\) ташуучу толкун узундугу, \(h_a\) антеннанын жерден бийиктиги, \(h_v\) өсүмдүк катмарынын эффективдүү фаза борборунун бийиктиги. \(R_1-R_2\) түз жана чагылган жолдордун ортосундагы геометриялык жол айырмасын билдирет.

Теңдеменин негизги жыйынтыгы: антенна менен эффективдүү чагылдыруучу фаза борборунун ортосундагы бийиктик айырмасы чоңойгон сайын спутниктин көтөрүлүү бурчуна байланышкан фаза өзгөрүшү да өзгөрөт. Ошондуктан SNR термелүү жыштыгын геометриялык бийиктик маалыматына айландырууга болот.

Эмне үчүн 5°–30° көтөрүлүү бурчтары гана колдонулду?

Спутник асманда көтөрүлгөн сайын түз жана чагылган сигнал жолдорунун геометриялык айырмасы азаят. Бул учурда multipath интерференциялык термелүүлөрдүн амплитудасы жана айырмалануу мүмкүнчүлүгү төмөндөйт. Булак 5°–30° аралыгында интерференция белгиси айкыныраак болорун көрсөткөн жана бул аралыктан тышкаркы SNR маанилерин маскалаган.

Төмөн көтөрүлүү бурчтарында чагылган сигналдын кабылдагычка жетишине салым кошкон бет аянты биринчи Fresnel зонасы менен болжолдуу көрсөтүлөт. Булак жарым окторду жана чагылуу чекитинин антеннага чейинки аралыкты:

\[ 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} \]

катыштары менен берет. \(a\) жана \(b\) Fresnel эллипсинин болжолдуу жарым окторун, \(d\) чагылуу чекитинин антеннага горизонталдык аралыкты билдирет.

Изилдөөдөгү бардык ылайыктуу спутник издери бирге бааланганда, GNSS-IR ар бир талаа кампаниясында болжол менен 5000 m² талаа аянтынан чагылуу маалыматын алган деп эсептелген. Бул footprint чекиттик топурак сенсорунан же бир сызгыч өлчөөсүнөн кыйла кең аянтты билдирет; бирок аянт ичиндеги гетерогендүүлүк да өлчөө белгисиздигинин бир бөлүгүнө айланат.

Талаа жана фенологиялык контекст

Эксперименттик аянт Италиянын түндүк-батышындагы Piedmont аймагында, Pogliola жергесинде жайгашкан 2 гектар туруктуу көп түрдүү шалбаа. Өсүмдүк катмары негизинен Gramineae жана Leguminosae түрлөрүнөн турат. 2025-жылы талаа 20 май, 8 июль жана 13 августта чабылып, 15 сентябрда жайыт катары пайдаланылган.

Sentinel-2 эмне үчүн колдонулду?

Sentinel-2 маалыматы GNSS-IR моделинин киргизүүсү эмес. Изилдөө GNSS-IR термелүү амплитудасы сезондук өсүмдүк катмарынын тыгыздыгынын өзгөрүшүнө сезимталбы же жокпу деген суроону көз карандысыз байкоо менен баалоо үчүн NDVI колдонгон.

1 апрель–1 октябрь мезгилинен жалпы 90 Sentinel-2 Level-2A Harmonized сүрөтү Google Earth Engineде иштетилип, булут/cirrus/көлөкө пикселдери Scene Classification Layer аркылуу чыгарылган. Ар бир жарактуу датада 10 m NIR жана кызыл диапазондордон NDVI эсептелип, төрт GNSS кампаниясы үчүн аянт медианасы алынган.

NDVI бул изилдөөдө биомассанын же Leaf Area Indexтин түз сандык өлчөөсү эмес. Булак NDVIни жашыл өсүмдүк катмары жана катмар тыгыздыгынын сезондук өзгөрүүлөрүн билдирген салыштырмалуу proxy катары гана чечмелеген.

Жердеги шилтемелер

Топурак нымдуулугу үчүн Drill & Drop Bluetooth сыйымдуулук сенсору 10 cm тереңдикке көмүлүп, көлөмдүк суу курамы саат сайын жазылып турган. GNSS кампания күндөрүндө 08:00–14:00 аралыгындагы SMC маанилеринин орточосу алынган.

Өсүмдүк катмарынын бийиктиги \(h_v\) болсо GNSS антеннасынын айланасындагы жети чекитте жерге тик кармалган сызгыч менен өлчөнгөн. Жер өлчөөлөрүнүн орточо мааниси жана стандарттык четтөөсү эсептелген.

SNRден dSNRге

Чийки SNR убакыт катары чагылган сигналды гана камтыбайт. Антеннанын gain pattern-и жана атмосфералык таасирлерден улам түз сигнал спутник көтөрүлгөн сайын жылмакай тренд түзөт. Булак бул төмөн жыштыктагы детерминисттик трендди \(\sin\theta\) көз карандысыз өзгөрмөсүнө карата экинчи даражадагы модель менен алып салып, орточо мааниси нөлгө жакын трендден тазаланган SNR (dSNR) катарын түзгөн.

Бул процесстен кийин dSNR ичиндеги үзгүлтүксүз термелүүлөр негизинен беттен чагылган когеренттүү компоненттин геометриялык жана электромагниттик таасирлерин алып жүрөт.

GNSS-IR өсүмдүктүн бийиктигин кантип баалайт?

GNSS-IR өсүмдүктүн бийиктигин түз визуалдык өлчөбөйт; dSNR термелүүлөрүнүн \(\sin\theta\) боюнча доминанттык жыштыгын Lomb–Scargle периодограммасы менен аныктайт, бул жыштыкты антенна менен эффективдүү чагылдыруучу фаза борборунун ортосундагы рефлектор бийиктигине айлантат жана рефлектор бийиктигин антенна бийиктигинен алып, өсүмдүк катмарынын фаза борборунун бийиктиги \(h_v\) алынат.

Lomb–Scargle периодограммасы

Булак dSNR катарынын \(\sin\theta\) боюнча гармоникалык курамын Lomb–Scargle периодограммасы менен талдайт. Колдонулган негизги матрицалык түр:

\[ 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} \]

жана:

\[ 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} \]

түрүндө. \(y\) dSNR байкоолорун, \(n\) убакыт катарындагы байкоолордун санын, \(\omega=2\pi f\) бурчтук спектралдык өзгөрмөнү билдирет.

Бул жердеги \(f\) убакыт боюнча Hertz бирдигиндеги жыштык эмес. Булак муну \(\sin\theta\) бирдигине туура келген термелүү циклдери катары аныктайт. Издөө аралыгы 0.1–50 cycles per unit of \(\sin\theta\), кадам 0.01.

Доминанттык жыштыктан өсүмдүк катмарынын бийиктигине

Булак Теңдеме (4a)ны басылган түрүндө төмөнкүдөй берет:

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

Булактын нотациясы боюнча эскертүү: Ошол эле бөлүмдө \(f\) \(\sin\theta\) бирдигине туура келген термелүү жыштыгы катары аныкталганы менен, Теңдеме (4a)дагы туунду өзгөрмөсү булак PDFинде \(t\) болуп басылган. Verianla бул туюнтманы булакта жазылбаган башка өзгөрмө менен үнсүз алмаштырбайт.

Жыштык алынгандан кийин рефлектор бийиктиги:

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

жана өсүмдүк фаза борборунун бийиктиги:

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

түрүндө эсептелет.

\(h_a\) антеннанын жерден бийиктиги болгон 2 m. \(RH\) антенна менен эффективдүү чагылдыруучу фаза борборунун ортосундагы аралыкты, \(h_v\) бул фаза борборунун жерден бийиктигин билдирет.

Сапатты көзөмөлдөө

Фильтрлөөдөн кийин башында ар бир кампанияда болжол менен 100 из болгондугуна карабай, жалпы спутник издеринин %8ден азы гана акыркы моделдөөгө калган. Бул калган маалыматтын сапатын жана физикалык чечмеленишин жогорулатат, бирок татаал чачыроо шарттарын чагылдыруу жагынан тандоо таасирин жаратышы мүмкүн.

Сызыктуу эмес dSNR модели

Доминанттык геометриялык жыштык жана \(h_v\) аныкталгандан кийин булак термелүү амплитудасы менен диэлектрдик фазаны ажыратуу үчүн экинчи этапта сызыктуу эмес гармоникалык моделди колдонгон:

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

Моделдин колдонулган түрү:

\[ 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} \]

түрүндө.

\(A_0\) негизги амплитуда масштабы, \(m_0\) көтөрүлүү бурчуна көз каранды амплитуда басаңдоосун, \(\phi_0\) чагылдыруучу чөйрөнүн диэлектрдик жүрүм-туруму менен байланышкан болжолдуу туруктуу фаза жылышын билдирет. \(h_v\) бул этапта эркин кайра оптималдаштырылган эмес; мурдагы Lomb–Scargle жыштык талдоосунан алынган геометриялык маалымат туруктуу калтырылган.

Сызыктуу эмес fit R чөйрөсүндө `minpack.lm` пакетиндеги `nlsLM` функциясы аркылуу bounded least-squares Levenberg–Marquardt алгоритми менен жүргүзүлгөн. Баштапкы амплитуда байкалган максималдуу dSNR маанисине, \(m_0=-0.1\)ге, \(\phi_0=0\)гө коюлган. Чектери \(A>0\), \(m\le0\) жана \(\phi\in[-\pi,\pi]\) түрүндө.

Чагылган талаанын физикалык модели

Булак когеренттүү жана когерентсиз чагылууну:

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

менен түшүндүрөт.

\(\Gamma\) топурак–өсүмдүк чөйрөсүнүн эффективдүү комплекстүү диэлектрдик өтүмдүүлүгүнө \(\widetilde{\varepsilon}\) көз каранды Fresnel коэффициенти. Туш келди багытталган жалбырак, сабак, боштук жана нымдуулук гетерогендүүлүктөрү жараткан \(\widetilde{S}_{incoh}\) компоненти детерминисттик термелүү жаратпастан, ызы-чуу фонун көтөрөт.

Ошондуктан dSNRдин маанилүү гармоникалык компоненти чагылган талаанын бүтүндөй эмес, жетиштүү когеренттүү калган бөлүгүнүн түз сигнал менен интерференциясынын байкалуучу натыйжасы.

Фаза борбору эмне үчүн чыныгы өсүмдүк чокусунан төмөн болушу мүмкүн?

Шалбаа же талаа өсүмдүк катмары микротолкун жагынан бир катуу бет эмес. Жалбырак жана сабактар ар башка бийиктиктерде чачыроо жаратат. Булак бул бөлүштүрүлгөн түзүлүштү:

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

менен берет.

\(F_v(z,\theta)\) өсүмдүк катмарынын тик рефлективдүүлүк профили; \(z_0\) жер деңгээли, \(h_{top}\) физикалык өсүмдүк катмарынын бийиктиги жана \(k=4\pi/\lambda\) толкун саны.

Бул бөлүштүрүүнүн фаза борбору:

\[ 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} \]

түрүндө аныкталат.

Теңдеме (7) GNSS-IR эмне үчүн сызгыч менен өлчөнгөн эң жогорку жалбырак деңгээлин сөзсүз бербей турганын түшүндүрөт. Өсүмдүктүн терең катмарларында басаңдоо жана чачыроо көбөйгөндүктөн когеренттүү чагылуунун салмак борбору физикалык өсүмдүк чокусунан төмөн болушу мүмкүн. Изилдөөдө байкалган терс бийиктик biasы ушул фаза-борбор жүрүм-туруму менен байланышкан.

Изилдөөнүн Ыкмасы жана Жыйынтыктары

Маалымат чогултуу жана текшерүү чынжыры

Өсүмдүк катмарынын бийиктиги боюнча жыйынтык

Бийиктик баалоолору айрыкча июнь–июль мезгилинде жогорураак, август–сентябрь мезгилинде болсо чабуудан кийинки төмөнүрөөк өсүмдүк катмарынын шарттарын көзөмөлдөгөн. Булак июнь–июль мезгилинде өсүмдүк бийиктиги болжол менен 50 cm болгонун, кийинки чабуу мезгилдери бийиктик жана NDVI маанилерин азайтканын талкуулайт.

Өсүмдүк катмарынын тыгыздыгы жана dSNR амплитудасы

Sentinel-2 NDVI менен θ = 15° үчүн эсептелген dSNR моделинин амплитудасынын ортосунда:

r = −0.81

байланыш табылган.

Терс белги өсүмдүк катмары тыгыздаган сайын dSNR термелүү амплитудасы азаярын көрсөтөт. Булактын физикалык түшүндүрмөсү: жалбырак аянты, биомасса жана суу курамы көбөйгөн сайын L-band когеренттүү чагылуу алсырайт.

Бул табылга NDVI GNSS-IR тарабынан бааланган дегенди билдирбейт. NDVI GNSS-IR амплитуда параметри өсүмдүк катмарынын өнүгүшүнө сезимталбы же жокпу дегенди баалоочу көз карандысыз proxy гана.

Топурак нымдуулугу жана диэлектрдик фаза

Жердеги soil moisture content (SMC) өлчөөсү менен \(|\phi_{die}|\) ортосунда:

r = 0.86

оң байланыш табылган.

Топурактын суу курамы көбөйгөндө топурак–өсүмдүк чөйрөсүнүн комплекстүү диэлектрдик өтүмдүүлүгү өзгөргөндүктөн Fresnel чагылуу коэффициентинин фазасы да өзгөрүшү мүмкүн. Ошондуктан булак \(\phi_0\) параметрин нымдуулуктан келип чыккан диэлектрдик өзгөрүүлөргө сезимтал байкалуучу чоңдук катары чечмелейт.

GNSS-IR бул изилдөөдө өсүмдүк түзүлүшү менен топурак нымдуулугун толук ажырата алабы?

Жок. Изилдөө өсүмдүк катмарынын бийиктиги, dSNR амплитуда басаңдоосу жана диэлектрдик фаза тиешелүүлүгүнө жараша түзүлүш, катмар тыгыздыгы жана нымдуулук өзгөрүүлөрүнө ар башка сезимталдык көрсөткөнүн колдойт; бирок чагылган GNSS сигналы топурак менен өсүмдүк катмары чогуу түзгөн электромагниттик чөйрөдөн келгендиктен, структуралык жана диэлектрдик салымдар бул маалымат топтомунда толук жана так ажыратылган эмес.

Үч GNSS-IR байкалуучу чоңдугунун толуктоочу мааниси

Бул таблица үч параметр толугу менен көз карандысыз физикалык өлчөөлөр дегенди билдирбейт. Ошол эле чагылуу процессинде геометрия, өсүмдүк катмарынын тыгыздыгы, суу курамы жана чачыроо түзүлүшү бири-бирине таасир этиши мүмкүн. Булактын жыйынтыгы — ар башка GNSS-IR параметрлери ар башка физикалык компоненттерге басымдуу сезимталдык көрсөтөт; толук separation азырынча камсыздалган эмес.

Топурак нымдуулугу жыйынтыктарындагы ыкма–жыйынтык шайкешсиздиги

Ыкма бөлүмү NDVI жана SMCни түз retrieval максаты катары аныктабайт жана бул эки шилтеме менен байланышты Pearson корреляциясы аркылуу гана баалоо керектигин көрсөтөт. Ошого карабастан жыйынтык бөлүмүндө SMC үчүн кошумча “ME = −5.3” жана “MAPE = %21” жазылган.

Булак \(\phi_{die}\) менен SMCнин ортосунда түз SMC баалоосу кандай алынганын жана бул эки ката метрикасы кайсы эсептөө чынжырынан чыкканын ыкма бөлүмүндө түшүндүрбөйт. Ошондуктан Verianla бул эки санды тастыкталган “GNSS-IR soil-moisture retrieval error” метрикасы катары колдонбойт. Ишенимдүү түрдө колдоого алынган сандык байланыш Pearson \(r=0.86\).

Корреляциялардын статистикалык чеги

Бардык негизги корреляциялар болгону төрт талаа кампаниясынын агрегатталган маанилерине негизделет: \(n=4\). Ошондуктан булак корреляцияларды катаал статистикалык жалпылоо катары эмес, физикалык байланыштардын көрсөткүчү катары чечмелөө керектигин ачык белгилейт.

Төрт дата бирдей шарттардын туш келди тандалган кайталоолору эмес; чабууга чейинки жана чабуудан кийинки ар башка фенологиялык режимдерди билдирет. Бул ар түрдүүлүк ыкманын ар башка талаа абалдарына реакциясын көрсөтүүгө жардам берет, бирок үлгү санын көбөйтпөйт.

Фильтрлөө жана selection bias

Сапат фильтрлери SNR ичинде физикалык жактан чечмелене турган когеренттүү multipath термелүүсү бар издерди тандайт. Бул мамиле ызы-чуулуу жана туруксуз чечимдерди азайтат.

Мунун экинчи жагы — гетерогендүү өсүмдүк катмары, татаал чагылуу геометриясы же тезирээк өзгөрүүчү нымдуулук шарттарында модель божомолдоруна азыраак туура келген спутник издери маалымат топтомунан көбүрөөк чыгарылышы мүмкүн. Авторлор ошондуктан акыркы издер бардык талаа шарттарынын калыс үлгүсү эмес, модель эң ачык иштеген байкоолор экенин белгилейт.

Lomb–Scargle бийиктигиндеги терс bias

Lomb–Scargle периодограммасы тең аралыкта эмес SNR үлгүлөрүндө доминанттык термелүү жыштыгын табуу үчүн күчтүү курал; бирок чектелген спутник бурчу терезеси spectral leakage жаратышы мүмкүн. Кошумча чагылуу жолдору да негизги периодограмма чокусуна таасир этиши мүмкүн.

Булак бул таасирлер доминанттык жыштыкты бир аз төмөн жылдырып, бул рефлектор бийиктигин баалоодо bias болуп чагылышы мүмкүн экенин талкуулайт. Мындан тышкары негизги физикалык себептердин бири өсүмдүк катмары көлөмдүк чачыраткыч болушу: фаза борбору физикалык үстүңкү беттен төмөн турат.

Изилдөө колдогон жыйынтыктар

Изилдөө төрт ар башка фенологиялык талаа шартында GNSS-IR dSNR сигналы өсүмдүк катмарынын түзүлүшү жана нымдуулук менен байланышкан системалуу өзгөрүүлөрдү камтый турганын колдойт.

Изилдөө Lomb–Scargle негизиндеги GNSS-IR рефлектор бийиктиги шалбаа өсүмдүк катмарындагы бийиктик өзгөрүшүн көзөмөлдөй аларын жана сызгыч шилтемеси менен күчтүү байланыш көрсөткөнүн колдойт.

Изилдөө dSNR амплитуда басаңдоосу сезондук өсүмдүк катмарынын тыгыздыгы/NDVI өзгөрүшүнө, диэлектрдик фаза параметри талаа нымдуулугунун өзгөрүшүнө сезимтал экенин колдойт.

Изилдөө GNSS-IR бир байкоо инфраструктурасынан бир нече айыл чарба биофизикалык маалымат каналдарын алуу мүмкүнчүлүгүнө ээ экенин колдойт.

Изилдөө колдобогон жыйынтыктар

Изилдөө GNSS-IR бардык өсүмдүк түрлөрүндө %18 өсүмдүк катмарынын бийиктик катасы менен иштей турганын көрсөтпөйт.

Изилдөө \(\phi_{die}\) өсүмдүк катмарынын таасиринен толук тазаланган таза топурак нымдуулугу өлчөөсү экенин көрсөтпөйт.

Изилдөө талааны сугаруу автоматтык түрдө GNSS-IR менен башкарылганын же суу керектөөсү азайтылганын көрсөтпөйт.

Изилдөө арзан смартфон GNSS системалары ошол эле натыйжалуулукту бергенин көрсөтпөйт; смартфон колдонуу келечектеги изилдөө багыты.

Изилдөө жүгөрү, буудай, соя же башка талаа өсүмдүктөрүндө ошол эле байланыштар тастыкталганын көрсөтпөйт.

Негизги чектөөлөр

Келечектеги изилдөө багыттары

Булак убакыт боюнча тез-тез үлгүлөөнү, көбүрөөк аянтты, ар башка өсүмдүк архитектураларын жана гетерогендүү шарттарда текшерүүнү сунуштайт. Өзгөчө maize, wheat жана soybean келечекте сыналышы пландалган өсүмдүктөрдүн арасында көрсөтүлгөн.

Ийкемдүү multi-frequency же time-varying harmonic моделдер, robust regression мамилелери жана өсүмдүк катмарынын гетерогендүүлүгүн түздөн-түз эсепке алган физикалык scattering моделдер көбүрөөк GNSS изин жоготпой параметр алуу үчүн изилдениши мүмкүн.

Арзан кош жыштыктуу GNSS кабылдагычтары жана смартфон сенсорлору да келечектеги масштабдалуучу колдонуу жолдорунун арасында көрсөтүлгөн. Бул мүмкүнчүлүк учурдагы изилдөөнүн текшерилген жыйынтыгы эмес, авторлордун келечектеги өнүктүрүү багыты.

Булак жана Ыкма Жөнүндө Эскертүү

Түпнуска аталыш: Monitoring Crop Structure and Moisture Using GNSS Interferometric Reflectometry Based on SNR Modeling

Авторлор: Samuele De Petris; Enrico Borgogno-Mondino

Кат алышуучу автор: Enrico Borgogno-Mondino

Мекеме: Department of Agriculture, Forestry and Food Science, University of Turin, Grugliasco, Italy.

Журнал: Agronomy

Басмакана: MDPI

Том / Сан / Макала: 16 / 9 / 922

DOI: 10.3390/agronomy16090922

Даталар: Алынды: 25 март 2026; кайра каралды: 27 апрель 2026; кабыл алынды: 29 апрель 2026; жарыяланды: 1 май 2026.

Булак түрү: Рецензияланган эксперименттик изилдөө макаласы; талаа GNSS өлчөөсү, микротолкундуу интерферометриялык рефлектометрия, спектралдык талдоо жана сызыктуу эмес сигнал моделдөө.

Версия: Version of Record.

Лицензия: Creative Commons Attribution (CC BY). PDFдеги көрүнгөн лицензия текстинде версия номери өзүнчө көрсөтүлбөгөндүктөн лицензия версиясы божомолдонгон эмес.

Негизги аянт: Pogliola, Piedmont, Италиянын түндүк-батышы; 2 ha туруктуу көп түрдүү шалбаа.

GNSS кабылдагычы: Topcon HiPer Pro; GPS жана GLONASS; L1/L2; антенна бийиктиги 2 m; 1 Hz.

GNSS кампаниялары: 1 июнь, 4 июль, 2 август жана 7 сентябрь 2025.

Оптикалык жардамчы маалымат: Copernicus Sentinel-2 Level-2A Harmonized; 90 сүрөт; 1 апрель–1 октябрь 2025; Google Earth Engine.

Топурак нымдуулугу шилтемеси: Drill & Drop Bluetooth сыйымдуулук сенсору; 10 cm тереңдик; сааттык volumetric water content.

Өсүмдүк катмарынын бийиктиги шилтемеси: GNSS антеннасынын айланасындагы жети чекитте кол менен сызгыч өлчөөсү.

Негизги сигнал иштетүү: Ascending/descending спутник жаасын ажыратуу; 5°–30° elevation mask; quadratic detrending; Lomb–Scargle периодограммасы; nonlinear harmonic regression.

LSP жыштык аралыгы: 0.1–50 cycles per unit of sin θ; кадам = 0.01.

Сапат фильтрлери: P2N > 3 dB-Hz; Δθ > 15°; 0.05 m < h_v < 2 m; nonlinear model R² > 0.6.

Сызыктуу эмес fit: R `nlsLM`, `minpack.lm`; bounded least-squares Levenberg–Marquardt.

Негизги жыйынтыктар: Өсүмдүк катмарынын бийиктиги r = 0.89, ME = −0.067 m, MAPE = %18; NDVI менен dSNR амплитудасы r = −0.81; SMC менен абсолюттук диэлектрдик фаза r = 0.86.

Статистикалык үлгү чектөөсү: Негизги салыштыруулар төрт кампаниянын агрегатталган маанилерине негизделет, n = 4.

Каржылоо: Изилдөө тышкы каржылоо алган эмес.

Маалымат жеткиликтүүлүгү: Маалыматтар суроо боюнча берилери билдирилген.

Кошумча материал: Table S1 ар бир сапатты көзөмөлдөө этабынан кийин калган спутник издеринин санын камтыйт.

Кызыкчылыктардын кагылышы: Авторлор кызыкчылыктардын кагылышы жок экенин билдиришкен.

Авторлордун салымы: Samuele De Petris концептуалдаштыруу, ыкма, программалык камсыздоо, текшерүү, формалдык талдоо, маалымат курациясы, баштапкы текст жана визуалдаштыруу ролдорун аткарган. Изилдөө жана ресурстар эки автор тарабынан бөлүшүлгөн. Enrico Borgogno-Mondino кеңеш берүү, долбоорду башкаруу жана каржылоону алуу ролдорун аткарган; эки автор тең карап чыгуу жана редакциялоого салым кошкон.

Булак ичиндеги нотация эскертүүсү: Текст Lomb–Scargle спектралдык өзгөрмөсүн sin θ бирдигине туура келген термелүү цикли катары аныктаганы менен Теңдеме (4a)да туунду өзгөрмөсү t болуп басылган. Бул туюнтма Verianla тарабынан үнсүз өзгөртүлгөн эмес.

Булак ичиндеги ыкма/жыйынтык эскертүүсү: Ыкма бөлүмү soil moisture жана NDVI салыштырууларында Pearson r гана колдонуларын айтса да, жыйынтык бөлүмүндө soil moisture үчүн кошумча ME = −5.3 жана MAPE = %21 берилген. Бул ката метрикаларынын чыгарылуу ыкмасы түшүндүрүлбөгөндүктөн, өз алдынча нымдуулук retrieval көрсөткүчү катары колдонулган эмес.

Негизги методологиялык чектөө: Бир аянт, бир өсүмдүк катмарынын тиби, төрт кампания, сейрек жер шилтемелери жана катуу сапат фильтрлери себептүү натыйжалар операциялык көп өсүмдүктүү айыл чарба системасы үчүн акыркы калибрлөө эмес.

Негизги илимий чечмелөө чеги: Изилдөө GNSS-IR байкалуучу чоңдуктары өсүмдүк түзүлүшү жана диэлектрдик/нымдуулук өзгөрүүлөрүнө ар башка сезимталдык көрсөтөрүн колдойт; бул эки физикалык салым толук ажыратылганын көрсөтпөйт.