Reflektometria Interferometriki ya GNSS (Global Navigation Satellite System Interferometric Reflectometry, GNSS-IR) ni mbinu ya reflektometria pasivu inayotumia mwingiliano kati ya ishara ya microwave ya GNSS inayofika moja kwa moja kutoka satelaiti na ishara inayofika kwenye kipokezi kilekile baada ya kuakisiwa na udongo au uoto ili kutoa vipimo vinavyohisi sifa za kijiometri na dielektriki za uso. Utafiti huu ulichunguza uwezo wa GNSS-IR kufuatilia urefu wa tabaka la mimea, msongamano wa uoto na mabadiliko ya unyevu kwa kuchakata data za GPS/GLONASS SNR zilizokusanywa katika kampeni nne za shambani kwenye malisho ya kudumu ya hektari 2 kaskazini-magharibi mwa Italia. Makadirio ya urefu wa tabaka la mimea yalionyesha uhusiano wa r = 0.89 na vipimo vya ardhini, huku kosa la wastani likiwa −0.067 m na MAPE %18. Uhusiano wa r = −0.81 uliripotiwa kati ya amplitudo ya mtikisiko wa dSNR na Sentinel-2 NDVI, na r = 0.86 kati ya sehemu ya awamu ya dielektriki na unyevu wa udongo wa shambani. Hata hivyo, kwa kuwa uchanganuzi ulitegemea tu thamani zilizojumlishwa za tarehe nne, matokeo hayapaswi kufasiriwa kama milinganyo ya urekebishaji inayoweza kujumlishwa kwa mazingira mapana ya kilimo, bali kama ushahidi wa awali wa shambani unaoonyesha kuwa viashiria vinavyoonekana vya GNSS-IR vina unyeti kwa muundo wa mimea na mabadiliko ya dielektriki.

Mbinu ya utafiti haijaribu kutoa sifa zote za shamba kutoka kwa kigezo kimoja cha GNSS-IR. Mzunguko wa mtikisiko na urefu wa kituo cha awamu unaohesabiwa kutokana nao ulihusishwa na muundo wima wa mimea; kupungua kwa amplitudo ya mtikisiko na msongamano wa tabaka la mimea pamoja na udhoofishaji wa sumakuumeme; na sehemu ya awamu ya dielektriki na mabadiliko katika permitiviti changamano ya mazingira ya udongo–mimea, hasa mabadiliko yanayotokana na unyevu. Ufafanuzi huu wa vigezo vingi unaonyesha kuwa ishara ileile ya GNSS inaweza kutoa taarifa zinazokamilishana kuhusu vipengele tofauti vya kimwili. Hata hivyo, utafiti haukuunda mfumo wa inversion unaotenganisha kikamilifu michango ya kimuundo na ya dielektriki.

Uvutio wa kiutendaji wa GNSS-IR ni kwamba, badala ya kuzalisha ishara mpya ya rada amilifu, hutumia ishara zilizopo za satelaiti za GNSS ambazo hutangazwa mfululizo. Hali hii inaweza kuruhusu mbinu kuunganishwa katika vituo vya GNSS vilivyowekwa, mashine za kilimo au, siku zijazo, vipokezi vya gharama nafuu. Kipokezi kilichotumika katika utafiti huu, hata hivyo, ni Topcon HiPer Pro ya kiwango cha utafiti; utendaji wa smartphone au vifaa vya gharama nafuu haukujaribiwa katika jaribio hili.

Tatizo la ufuatiliaji wa kilimo

Urefu wa mimea, msongamano wa tabaka la mimea na unyevu wa udongo ni vigezo muhimu vya kibaiofizikia vinavyoeleza michakato tofauti ya kilimo. Jiometri ya tabaka la mimea inahusiana na kunasa mwanga, transpiration na mwingiliano wa anga–mimea; unyevu wa udongo unahusiana na uotaji, ufyonzwaji wa virutubisho, michakato ya vijidudu na msongo wa maji. Kwa kuwa vigezo hivi vinaweza kubadilika kwa wakati na nafasi ndani ya shamba lilelile, vipimo vya nukta vilivyo haba vinaweza kutotosha kuwakilisha hali ya eneo lote.

Uchunguzi wa satelaiti wa macho unaweza kupiga picha maeneo makubwa mara kwa mara, lakini huathiriwa na mawingu, erosoli na hali ya mwangaza, na haupimi moja kwa moja muundo wa tabaka la mimea wa pande tatu. LiDAR inaweza kubainisha jiometri ya tabaka la mimea kwa usahihi mkubwa, lakini haijaenea kama uchunguzi endelevu wa satelaiti wenye azimio kubwa la wakati. Mifumo ya rada na microwave pasivu ni nyeti kwa unyevu na maudhui ya maji ya mimea, lakini inaweza kuwa na vikwazo vya azimio la wakati/nafasi katika kiwango cha shamba.

Pengo la uchunguzi linalochunguzwa na utafiti huu ni kutafuta mbinu ya kiwango cha shamba ambayo inaweza kuwa nyeti kwa wakati mmoja kwa muundo wa mimea na mabadiliko ya dielektriki yanayotokana na unyevu, huku ikifanya kazi kwa marudio ya juu ya wakati.

GNSS-IR hufanyaje kazi katika ufuatiliaji wa kilimo?

GNSS-IR hupima pamoja ishara inayofika moja kwa moja kwenye antena kutoka satelaiti ya GNSS na ishara inayofika kwenye antena ileile baada ya kuakisiwa na udongo na uoto na kusafiri njia ndefu zaidi; tofauti ya awamu kati ya njia hizi mbili huunda mitikisiko ya kawaida ya mwingiliano katika SNR, na mzunguko, amplitudo na awamu ya mitikisiko hii zinaweza kubeba taarifa kuhusu jiometri ya uso, udhoofishaji wa ishara na sifa za dielektriki mtawalia.

Mwingiliano wa ishara ya moja kwa moja na iliyoakisiwa

Chanzo kinaeleza uwiano wa ishara-kwa-kelele unaoonekana na kipokezi cha GNSS kwa Mlinganyo (1a):

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

\(\theta\) ni pembe ya mwinuko wa satelaiti. \(\widetilde{S}\) inawakilisha ishara changamano ya jumla inayofika kwenye kipokezi, huku \(\widetilde{N}\) ikiwakilisha sehemu ya kelele.

Mwingiliano wa nyuga ya moja kwa moja na iliyoakisiwa katika nguvu ya jumla ya ishara iliyopokelewa umefunguliwa katika chanzo kama ifuatavyo:

\[ |\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\) ni amplitudo ya ishara ya moja kwa moja na \(A_r\) ni amplitudo ya ishara iliyoakisiwa. \(\cos(\Delta\phi)\) katika sehemu ya mwisho inaonyesha kwa nini SNR hupanda na kushuka kwa namna ya vipindi wakati satelaiti inaposonga.

Tofauti ya awamu imetolewa katika chanzo kama:

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

\(\lambda\) ni urefu wa wimbi la carrier, \(h_a\) ni urefu wa antena kutoka ardhini, na \(h_v\) ni urefu wa kituo cha awamu kinachofaa cha uoto. \(R_1-R_2\) inaeleza tofauti ya njia ya kijiometri kati ya njia ya moja kwa moja na iliyoakisiwa.

Hitimisho kuu la mlinganyo ni hili: kadiri tofauti ya urefu kati ya antena na kituo cha awamu kinachoakisi kwa ufanisi inavyoongezeka, ndivyo mabadiliko ya awamu yanayotegemea pembe ya mwinuko wa satelaiti yanavyobadilika. Kwa hiyo, mzunguko wa mtikisiko wa SNR unaweza kubadilishwa kuwa taarifa ya urefu wa kijiometri.

Kwa nini zilitumika tu pembe za mwinuko 5°–30°?

Satelaiti inapopanda angani, tofauti ya kijiometri kati ya njia ya ishara ya moja kwa moja na iliyoakisiwa hupungua. Katika hali hii amplitudo na uwezekano wa kutenganisha mitikisiko ya multipath interference hupungua. Chanzo kilionyesha kuwa alama ya mwingiliano iko wazi zaidi katika masafa ya 5°–30° na kilificha thamani za SNR zilizo nje ya masafa haya.

Katika pembe za chini za mwinuko, eneo la uso linalochangia ishara iliyoakisiwa kufika kwenye kipokezi linawakilishwa takriban na eneo la kwanza la Fresnel. Chanzo kinatoa nusu-mihimili na umbali wa sehemu ya uakisi kutoka antena kwa:

\[ 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\) na \(b\) zinawakilisha takriban nusu-mihimili ya duaradufu ya Fresnel, huku \(d\) ikiwakilisha umbali wa mlalo wa sehemu ya uakisi kutoka antena.

Njia zote zinazofaa za satelaiti katika utafiti zilipotathminiwa kwa pamoja, ilikadiriwa kuwa GNSS-IR ilipokea taarifa ya uakisi kutoka takriban 5000 m² za eneo la shamba katika kila kampeni ya shambani. Footprint hii inawakilisha eneo kubwa zaidi kuliko sensa ya udongo ya nukta au kipimo kimoja cha rula; lakini utofauti wa ndani ya eneo pia huwa sehemu ya kutokuwa na uhakika wa kipimo.

Muktadha wa shamba na fenolojia

Eneo la majaribio ni malisho ya kudumu ya aina nyingi ya hektari 2 katika Pogliola, Piedmont, kaskazini-magharibi mwa Italia. Uoto una aina nyingi za Gramineae na Leguminosae. Mwaka 2025, shamba lilikatwa tarehe 20 Mei, 8 Julai na 13 Agosti na kulishwa tarehe 15 Septemba.

Kwa nini Sentinel-2 ilitumika?

Data ya Sentinel-2 si ingizo la modeli ya GNSS-IR. Utafiti ulitumia NDVI kama uchunguzi huru ili kutathmini kama amplitudo ya mtikisiko wa GNSS-IR ni nyeti kwa mabadiliko ya msimu ya msongamano wa tabaka la mimea.

Jumla ya picha 90 za Sentinel-2 Level-2A Harmonized kutoka kipindi cha 1 Aprili–1 Oktoba zilichakatwa katika Google Earth Engine, na pikseli za wingu/cirrus/kivuli zilitolewa kwa Scene Classification Layer. Kwa kila tarehe halali, NDVI ilihesabiwa kutoka bendi za 10 m NIR na nyekundu na mediani ya eneo ikachukuliwa kwa kampeni nne za GNSS.

Katika utafiti huu NDVI si kipimo cha moja kwa moja cha kiasi cha biomasi au Leaf Area Index. Chanzo kiliitafsiri NDVI tu kama proxy ya kulinganisha inayowakilisha mabadiliko ya msimu ya uoto wa kijani na msongamano wa tabaka la mimea.

Marejeo ya ardhini

Kwa unyevu wa udongo, sensa ya uwezo ya Drill & Drop Bluetooth ilizikwa kwa kina cha 10 cm na maudhui ya maji kwa ujazo yakarekodiwa kila saa. Katika siku za kampeni za GNSS, thamani za SMC kati ya 08:00–14:00 zilifanyiwa wastani.

Urefu wa tabaka la mimea \(h_v\) ulipimwa kwa rula iliyoshikiliwa wima ardhini katika pointi saba kuzunguka antena ya GNSS. Wastani na mkengeuko sanifu wa vipimo vya ardhini zilihesabiwa.

Kutoka SNR hadi dSNR

Mfululizo mbichi wa SNR haujumuishi tu ishara iliyoakisiwa. Kwa sababu ya gain pattern ya antena na athari za anga, ishara ya moja kwa moja huunda mwelekeo laini satelaiti inapopanda. Chanzo kiliondoa mwelekeo huu wa chini wa masafa wenye tabia ya deterministiki kwa modeli ya daraja la pili dhidi ya \(\sin\theta\) kama kigezo huru, na kuunda mfululizo wa SNR uliotolewa mwelekeo (dSNR) wenye wastani karibu na sufuri.

Baada ya mchakato huu, mitikisiko ya kawaida ndani ya dSNR hubeba hasa athari za kijiometri na sumakuumeme za sehemu koherent ya uakisi kutoka uso.

GNSS-IR hukadiriaje urefu wa mimea?

GNSS-IR haipimi urefu wa mimea moja kwa moja kwa kuona; hutambua mzunguko dominant wa mitikisiko ya dSNR dhidi ya \(\sin\theta\) kwa periodogram ya Lomb–Scargle, hubadilisha mzunguko huo kuwa urefu wa reflector kati ya antena na kituo cha awamu kinachoakisi kwa ufanisi, kisha hutoa urefu wa reflector kutoka urefu wa antena ili kupata urefu wa kituo cha awamu cha tabaka la mimea \(h_v\).

Periodogram ya Lomb–Scargle

Chanzo kinachanganua maudhui ya harmoniki ya mfululizo wa dSNR dhidi ya \(\sin\theta\) kwa periodogram ya Lomb–Scargle. Muundo mkuu wa matriki uliotumika ni:

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

na:

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

kwa muundo huo. \(y\) inawakilisha uchunguzi wa dSNR; \(n\) ni idadi ya uchunguzi katika mfululizo wa wakati; na \(\omega=2\pi f\) ni kigezo cha spektra ya pembe.

Hapa \(f\) si masafa ya Hertz kwa wakati. Chanzo kinaifafanua kama mizunguko ya mtikisiko kwa kila kitengo cha \(\sin\theta\). Masafa ya utafutaji ni 0.1–50 cycles per unit of \(\sin\theta\), na hatua ni 0.01.

Kutoka mzunguko dominant hadi urefu wa tabaka la mimea

Chanzo kinatoa Mlinganyo (4a) kwa namna iliyochapishwa kama ifuatavyo:

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

Dokezo la uandishi wa chanzo: Ingawa katika sehemu hiyo hiyo \(f\) inafafanuliwa kama mzunguko wa mtikisiko kwa kila kitengo cha \(\sin\theta\), kigezo cha derivativi katika Mlinganyo (4a) kimechapishwa kama \(t\) katika PDF ya chanzo. Verianla haibadilishi kimya kimya usemi huu na kigezo kingine ambacho hakijaandikwa katika chanzo.

Baada ya mzunguko kupatikana, urefu wa reflector ni:

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

na urefu wa kituo cha awamu cha mimea ni:

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

huhesabiwa kwa namna hiyo.

\(h_a\) ni urefu wa antena kutoka ardhini, ambao ni 2 m. \(RH\) ni umbali kati ya antena na kituo cha awamu kinachoakisi kwa ufanisi; \(h_v\) ni urefu wa kituo hicho cha awamu juu ya ardhi.

Udhibiti wa ubora

Baada ya uchujaji, ingawa mwanzoni kulikuwa na takriban njia 100 kwa kila kampeni, chini ya %8 ya njia zote za satelaiti zilibaki kwa uundaji wa modeli wa mwisho. Hali hii huongeza ubora na uwezekano wa kutafsiri kimwili data iliyobaki, lakini inaweza kuleta athari ya uteuzi katika uwakilishi wa hali changamano za usambazaji.

Modeli isiyo ya mstari ya dSNR

Baada ya mzunguko dominant wa kijiometri na \(h_v\) kubainishwa, chanzo kilitumia modeli ya harmoniki isiyo ya mstari katika hatua ya pili ili kutenganisha amplitudo ya mtikisiko na awamu ya dielektriki:

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

Muundo uliotumika wa modeli ni:

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

kwa namna hiyo.

\(A_0\) ni kipimo cha msingi cha amplitudo, \(m_0\) inawakilisha udhoofishaji wa amplitudo unaotegemea pembe ya mwinuko, na \(\phi_0\) inawakilisha mabadiliko ya awamu takriban thabiti yanayohusishwa na tabia ya dielektriki ya mazingira yanayoakisi. \(h_v\) haikuoptimishwa tena kwa uhuru katika hatua hii; taarifa ya kijiometri iliyotokana na uchanganuzi wa mzunguko wa Lomb–Scargle ilishikiliwa kuwa thabiti.

Fit isiyo ya mstari ilifanywa katika mazingira ya R kwa kutumia kazi `nlsLM` katika kifurushi `minpack.lm`, kwa algoritimu ya bounded least-squares Levenberg–Marquardt. Amplitudo ya awali iliwekwa kwenye thamani ya juu zaidi ya dSNR iliyoonekana, \(m_0=-0.1\), na \(\phi_0=0\). Mipaka ni \(A>0\), \(m\le0\) na \(\phi\in[-\pi,\pi]\).

Modeli ya kimwili ya uga ulioakisiwa

Chanzo kinaeleza uakisi koherent na usio koherent kwa:

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

usemi huo.

\(\Gamma\) ni mgawo wa Fresnel unaotegemea permitiviti changamano inayofaa ya mazingira ya udongo–mimea \(\widetilde{\varepsilon}\). Sehemu \(\widetilde{S}_{incoh}\), inayotokana na mwelekeo nasibu wa majani, mashina, nafasi tupu na utofauti wa unyevu, huongeza sakafu ya kelele badala ya kutoa mtikisiko wa deterministiki.

Kwa hiyo, sehemu ya harmoniki yenye maana ya dSNR ni matokeo yanayoonekana ya mwingiliano wa si uga wote ulioakisiwa, bali sehemu yake inayobaki koherent vya kutosha, na ishara ya moja kwa moja.

Kwa nini kituo cha awamu kinaweza kuwa chini ya kilele halisi cha mimea?

Tabaka la mimea la malisho au shamba si uso mmoja mgumu kwa mtazamo wa microwave. Majani na mashina hutawanya ishara katika urefu tofauti. Chanzo kinaeleza muundo huu uliosambazwa kwa:

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

usemi huo.

\(F_v(z,\theta)\) ni profaili wima ya reflektiviti ya tabaka la mimea; \(z_0\) ni kimo cha ardhi, \(h_{top}\) ni urefu wa kimwili wa tabaka la mimea, na \(k=4\pi/\lambda\) ni namba ya wimbi.

Kituo cha awamu cha usambazaji huu kinafafanuliwa kama:

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

kwa namna hiyo.

Mlinganyo (7) unaeleza kwa nini GNSS-IR si lazima itoe kiwango cha juu kabisa cha jani kinachopimwa kwa rula. Kwa kuwa udhoofishaji na usambazaji huongezeka katika tabaka za ndani zaidi za mimea, kituo cha uzito cha uakisi koherent kinaweza kuwa chini ya kilele cha kimwili cha tabaka la mimea. Bias hasi ya urefu iliyoonekana katika utafiti imehusishwa na tabia hii ya kituo cha awamu.

Mbinu na Matokeo ya Utafiti

Mlolongo wa ukusanyaji na uthibitishaji wa data

Matokeo ya urefu wa tabaka la mimea

Makadirio ya urefu yalifuata hasa hali ya tabaka la mimea la juu zaidi katika Juni–Julai na la chini zaidi baada ya kukatwa katika Agosti–Septemba. Chanzo kinajadili kwamba urefu wa mimea ulikuwa karibu 50 cm katika Juni–Julai, na vipindi vya kukatwa vilivyofuata vilipunguza thamani za urefu na NDVI.

Msongamano wa tabaka la mimea na amplitudo ya dSNR

Kati ya Sentinel-2 NDVI na amplitudo ya modeli ya dSNR iliyohesabiwa kwa θ = 15°:

r = −0.81

ulipatikana uhusiano.

Alama hasi inaonyesha kuwa kadiri tabaka la mimea linavyozidi kuwa zito, amplitudo ya mtikisiko wa dSNR hupungua. Tafsiri ya kimwili ya chanzo ni kwamba kuongezeka kwa eneo la majani, biomasi na maudhui ya maji hudhoofisha uakisi koherent wa L-band.

Matokeo haya hayamaanishi kuwa NDVI ilikadiriwa na GNSS-IR. NDVI ni proxy huru tu inayotathmini kama kigezo cha amplitudo cha GNSS-IR ni nyeti kwa maendeleo ya uoto.

Unyevu wa udongo na awamu ya dielektriki

Kati ya kipimo cha ardhini cha soil moisture content (SMC) na \(|\phi_{die}|\):

r = 0.86

ulipatikana uhusiano chanya.

Kuongezeka kwa maudhui ya maji ya udongo hubadilisha permitiviti changamano ya mazingira ya udongo–mimea, kwa hiyo awamu ya mgawo wa uakisi wa Fresnel pia inaweza kubadilika. Kwa sababu hii, chanzo kinatafsiri kigezo \(\phi_0\) kama kiashiria kinachoonekana ambacho ni nyeti kwa mabadiliko ya dielektriki yanayotokana na unyevu.

Je, GNSS-IR inaweza kutenganisha kikamilifu muundo wa mimea na unyevu wa udongo katika utafiti huu?

Hapana. Utafiti unaunga mkono kwamba urefu wa tabaka la mimea, udhoofishaji wa amplitudo ya dSNR na awamu ya dielektriki huonyesha unyeti tofauti kwa muundo, msongamano wa tabaka la mimea na mabadiliko ya unyevu mtawalia; lakini kwa kuwa ishara ya GNSS iliyoakisiwa hutoka katika mazingira ya sumakuumeme yanayoundwa pamoja na udongo na uoto, michango ya kimuundo na ya dielektriki haikutenganishwa kikamilifu na kwa uhakika katika seti hii ya data.

Maana inayokamilishana ya viashiria vitatu vya GNSS-IR

Jedwali hili halimaanishi kuwa vigezo vitatu ni vipimo vya kimwili vinavyojitegemea kabisa. Katika mchakato uleule wa uakisi, jiometri, msongamano wa tabaka la mimea, maudhui ya maji na muundo wa usambazaji vinaweza kuathiriana. Hitimisho la chanzo ni kwamba vigezo tofauti vya GNSS-IR vinaonyesha unyeti mkuu kwa vipengele tofauti vya kimwili; separation kamili bado haijapatikana.

Kutokulingana kati ya mbinu na matokeo katika matokeo ya unyevu wa udongo

Sehemu ya mbinu haifafanui NDVI na SMC kama malengo ya moja kwa moja ya retrieval na inaeleza kuwa uhusiano na marejeo haya mawili utatathminiwa kwa Pearson correlation pekee. Hata hivyo, katika sehemu ya matokeo pia yameandikwa “ME = −5.3” na “MAPE = %21” kwa SMC.

Chanzo hakielezi katika sehemu ya mbinu jinsi makadirio ya moja kwa moja ya SMC yalivyotolewa kutoka uhusiano kati ya \(\phi_{die}\) na SMC, wala mlolongo wa hesabu uliotumika kupata metriki hizi mbili za kosa. Kwa hiyo, Verianla haitumii namba hizi mbili kama metriki iliyothibitishwa ya “GNSS-IR soil-moisture retrieval error”. Uhusiano wa kiasi unaoungwa mkono kwa uhakika ni Pearson \(r=0.86\).

Kikomo cha kitakwimu cha uhusiano

Mahusiano yote makuu yanategemea tu thamani zilizojumlishwa za kampeni nne za shambani: \(n=4\). Kwa hiyo, chanzo kinaeleza wazi kwamba mahusiano haya yanapaswa kufasiriwa kama viashiria vya uhusiano wa kimwili, si kama jumlisho kali za kitakwimu.

Tarehe nne si marudio ya bahati nasibu ya hali ileile; zinawakilisha mifumo tofauti ya fenolojia kabla na baada ya kukatwa. Utofauti huu husaidia kuonyesha kuwa mbinu inajibu kwa hali tofauti za shamba, lakini hauongezi ukubwa wa sampuli.

Uchujaji na selection bias

Vichujio vya ubora huchagua njia zenye mtikisiko wa multipath koherent unaoweza kutafsiriwa kimwili ndani ya SNR. Mbinu hii hupunguza suluhisho zenye kelele na zisizo thabiti.

Kwa upande mwingine, njia za satelaiti katika hali ya uoto usio sawa, jiometri changamano ya uakisi au hali za unyevu zinazobadilika haraka ambazo hazilingani vizuri na dhana za modeli zinaweza kuondolewa mara nyingi zaidi kutoka kwenye seti ya data. Kwa hiyo, waandishi wanaeleza kuwa njia za mwisho si sampuli isiyo na upendeleo ya hali zote za shamba, bali ni uchunguzi ambao modeli hufanya kazi kwa uwazi zaidi.

Bias hasi katika urefu wa Lomb–Scargle

Periodogram ya Lomb–Scargle ni zana yenye nguvu ya kupata mzunguko dominant katika sampuli za SNR zisizo na vipindi sawa; lakini dirisha finyu la pembe ya satelaiti linaweza kusababisha spectral leakage. Njia za ziada za uakisi pia zinaweza kuathiri kilele kikuu cha periodogram.

Chanzo kinajadili kuwa athari hizi zinaweza kusogeza mzunguko dominant kidogo chini na hilo linaweza kujitokeza kama bias katika makadirio ya urefu wa reflector. Zaidi ya hayo, moja ya sababu kuu za kimwili ni kwamba tabaka la mimea ni kisambazaji cha ujazo: kituo cha awamu kiko chini ya uso wa juu wa kimwili.

Matokeo yanayoungwa mkono na utafiti

Utafiti unaunga mkono kwamba katika hali nne tofauti za fenolojia za shambani, ishara ya GNSS-IR dSNR hubeba mabadiliko ya kimfumo yanayohusishwa na muundo wa tabaka la mimea na unyevu.

Utafiti unaunga mkono kwamba urefu wa reflector wa GNSS-IR unaotegemea Lomb–Scargle unaweza kufuatilia mabadiliko ya urefu wa tabaka la mimea katika malisho na unaonyesha uhusiano mkubwa na rejeleo la rula.

Utafiti unaunga mkono kwamba udhoofishaji wa amplitudo ya dSNR ni nyeti kwa mabadiliko ya msimu ya msongamano wa tabaka la mimea/NDVI, na kigezo cha awamu ya dielektriki ni nyeti kwa mabadiliko ya unyevu wa shambani.

Utafiti unaunga mkono kwamba GNSS-IR ina uwezo wa kutoa njia nyingi za taarifa za kibaiofizikia za kilimo kutoka katika miundombinu moja ya uchunguzi.

Matokeo yasiyoungwa mkono na utafiti

Utafiti hauonyeshi kwamba GNSS-IR itafanya kazi kwa kosa la %18 la urefu wa tabaka la mimea katika aina zote za mazao.

Utafiti hauonyeshi kwamba \(\phi_{die}\) ni kipimo safi cha unyevu wa udongo kilichoondolewa kabisa athari za uoto.

Utafiti hauonyeshi kwamba umwagiliaji wa shamba ulidhibitiwa moja kwa moja na GNSS-IR au kwamba matumizi ya maji yalipunguzwa.

Utafiti hauonyeshi kwamba mifumo ya GNSS ya smartphone ya gharama nafuu hutoa utendaji uleule; matumizi ya smartphone ni mwelekeo wa kazi ya baadaye.

Utafiti hauonyeshi kwamba mahusiano yale yale yamethibitishwa katika maize, wheat, soybean au mazao mengine ya shambani.

Vikwazo vikuu

Mielekeo ya kazi ya baadaye

Chanzo kinapendekeza usampulishaji wa mara nyingi zaidi kwa wakati, maeneo zaidi, miundo tofauti ya mazao na uthibitishaji katika hali zisizo sare. Hasa maize, wheat na soybean zimetajwa miongoni mwa mazao yanayopangwa kujaribiwa siku zijazo.

Modeli zinazonyumbulika zaidi za multi-frequency au time-varying harmonic, mbinu za robust regression na modeli za kimwili za scattering zinazoshughulikia moja kwa moja utofauti wa tabaka la mimea zinaweza kuchunguzwa ili kutoa vigezo bila kupoteza njia nyingi za GNSS.

Vipokezi vya GNSS vya gharama nafuu vya masafa mawili na sensa za smartphone pia zimetajwa kama njia zinazoweza kupanuliwa za matumizi ya baadaye. Uwezo huu si matokeo yaliyothibitishwa ya utafiti huu, bali ni mwelekeo wa maendeleo ya baadaye wa waandishi.

Maelezo ya Chanzo na Mbinu

Kichwa asili: Monitoring Crop Structure and Moisture Using GNSS Interferometric Reflectometry Based on SNR Modeling

Waandishi: Samuele De Petris; Enrico Borgogno-Mondino

Mwandishi wa mawasiliano: Enrico Borgogno-Mondino

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

Jarida: Agronomy

Mchapishaji: MDPI

Juzuu / Toleo / Makala: 16 / 9 / 922

DOI: 10.3390/agronomy16090922

Tarehe: Imepokelewa: 25 Machi 2026; imerekebishwa: 27 Aprili 2026; imekubaliwa: 29 Aprili 2026; imechapishwa: 1 Mei 2026.

Aina ya chanzo: Makala ya utafiti wa majaribio iliyopitiwa na wataalamu; kipimo cha GNSS shambani, reflektometria interferometriki ya microwave, uchanganuzi wa spektra na uundaji wa modeli ya ishara isiyo ya mstari.

Toleo: Version of Record.

Leseni: Creative Commons Attribution (CC BY). Kwa kuwa nambari ya toleo la leseni haijaonyeshwa tofauti katika maandishi yanayoonekana ya leseni kwenye PDF, toleo la leseni halijakisiwa.

Eneo kuu: Pogliola, Piedmont, kaskazini-magharibi mwa Italia; 2 ha za malisho ya kudumu ya aina nyingi.

Kipokezi cha GNSS: Topcon HiPer Pro; GPS na GLONASS; L1/L2; urefu wa antena 2 m; 1 Hz.

Kampeni za GNSS: 1 Juni, 4 Julai, 2 Agosti na 7 Septemba 2025.

Data saidizi ya macho: Copernicus Sentinel-2 Level-2A Harmonized; picha 90; 1 Aprili–1 Oktoba 2025; Google Earth Engine.

Rejeleo la unyevu wa udongo: Sensa ya uwezo ya Drill & Drop Bluetooth; kina cha 10 cm; hourly volumetric water content.

Rejeleo la urefu wa tabaka la mimea: Kipimo cha rula kwa mkono katika pointi saba kuzunguka antena ya GNSS.

Uchakataji mkuu wa ishara: Mgawanyo wa arc za ascending/descending; 5°–30° elevation mask; quadratic detrending; periodogram ya Lomb–Scargle; nonlinear harmonic regression.

Masafa ya mzunguko wa LSP: 0.1–50 cycles per unit of sin θ; hatua = 0.01.

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

Fit isiyo ya mstari: R `nlsLM`, `minpack.lm`; bounded least-squares Levenberg–Marquardt.

Matokeo makuu: Urefu wa tabaka la mimea r = 0.89, ME = −0.067 m, MAPE = %18; NDVI na amplitudo ya dSNR r = −0.81; SMC na awamu kamili ya dielektriki r = 0.86.

Kikomo cha sampuli ya kitakwimu: Ulinganisho mkuu unategemea thamani zilizojumlishwa za kampeni nne, n = 4.

Ufadhili: Utafiti haukupokea ufadhili wa nje.

Upatikanaji wa data: Imeripotiwa kuwa data zitatolewa kwa ombi.

Nyenzo za ziada: Table S1 ina idadi ya njia za satelaiti zinazobaki baada ya kila hatua ya udhibiti wa ubora.

Mgongano wa maslahi: Waandishi hawakuripoti mgongano wa maslahi.

Michango ya waandishi: Samuele De Petris alifanya uundaji wa dhana, mbinu, programu, uthibitishaji, uchanganuzi rasmi, uhifadhi wa data, rasimu ya awali na uwasilishaji wa picha. Utafiti na rasilimali ziligawanywa kati ya waandishi wote wawili. Enrico Borgogno-Mondino alifanya ushauri, usimamizi wa mradi na upatikanaji wa ufadhili; waandishi wote wawili walichangia mapitio na uhariri.

Onyo la notasheni ndani ya chanzo: Ingawa maandishi yanafafanua kigezo cha spektra cha Lomb–Scargle kama mizunguko ya mtikisiko kwa kila kitengo cha sin θ, kigezo cha derivativi katika Mlinganyo (4a) kimechapishwa kama t. Usemi huu haujabadilishwa kimya kimya na Verianla.

Onyo la mbinu/matokeo ndani ya chanzo: Sehemu ya mbinu inaeleza kuwa Pearson r pekee itatumika katika ulinganisho wa soil moisture na NDVI, lakini sehemu ya matokeo pia inaripoti ME = −5.3 na MAPE = %21 kwa soil moisture. Kwa kuwa namna ya kupata metriki hizi za kosa haijaelezwa, hazikutumika kama utendaji huru wa retrieval wa unyevu.

Kizuizi kikuu cha kimetodolojia: Kwa sababu ya eneo moja, aina moja ya tabaka la mimea, kampeni nne, marejeo machache ya ardhini na vichujio vikali vya ubora, matokeo si urekebishaji wa mwisho kwa mfumo wa kilimo wa mazao mengi unaoendeshwa kiutendaji.

Kikomo kikuu cha tafsiri ya kisayansi: Utafiti unaunga mkono kwamba viashiria vinavyoonekana vya GNSS-IR vinaonyesha unyeti tofauti kwa muundo wa mimea na mabadiliko ya dielektriki/unyevu; hauonyeshi kwamba michango hii miwili ya kimwili imetenganishwa kikamilifu.