GNSS İnterferometrik Reflektometri (Global Navigation Satellite System Interferometric Reflectometry, GNSS-IR), uydudan doğrudan gelen GNSS mikrodalga sinyali ile toprak veya bitki örtüsünden yansıdıktan sonra aynı alıcıya ulaşan sinyal arasındaki girişimi kullanarak yüzeyin geometrik ve dielektrik özelliklerine duyarlı ölçümler üreten pasif bir reflektometri yöntemidir. Bu çalışma, kuzeybatı İtalya'daki 2 hektarlık kalıcı çayırda dört saha kampanyasından toplanan GPS/GLONASS SNR verilerini işleyerek GNSS-IR'nin kanopi yüksekliği, bitki örtüsü yoğunluğu ve nem değişimlerini izleme potansiyelini araştırmıştır. Kanopi yüksekliği tahmini yer ölçümleriyle r = 0.89 korelasyon göstermiş, ortalama hata −0.067 m ve MAPE %18 olarak bulunmuştur. dSNR salınım genliği ile Sentinel-2 NDVI arasında r = −0.81, dielektrik faz terimi ile saha toprak nemi arasında r = 0.86 ilişki raporlanmıştır. Bununla birlikte analiz yalnız dört tarihin toplulaştırılmış değerlerine dayandığından sonuçlar geniş tarımsal koşullara genellenmiş kalibrasyon denklemleri olarak değil, GNSS-IR gözlenebilirlerinin bitki yapısı ve dielektrik değişimlere duyarlı olduğunu gösteren ilk saha kanıtı olarak yorumlanmalıdır.

Çalışmanın yaklaşımı tek bir GNSS-IR parametresinden bütün tarla özelliklerini çıkarmaya çalışmaz. Salınım frekansı ve bundan hesaplanan faz merkezi yüksekliği bitkinin dikey yapısıyla; salınım genliğindeki zayıflama kanopi yoğunluğu ve elektromanyetik sönümlemeyle; dielektrik faz terimi ise toprak–bitki ortamının kompleks permitivitesindeki, özellikle nem kaynaklı değişimlerle ilişkilendirilmiştir. Bu çoklu yorum, aynı GNSS sinyalinden farklı fiziksel bileşenler hakkında tamamlayıcı bilgi üretilebileceğini göstermektedir. Ancak çalışma yapısal ve dielektrik katkıları tamamen birbirinden ayıran bir inversion sistemi geliştirmemiştir.

GNSS-IR'nin pratik çekiciliği, yeni aktif radar sinyali üretmek yerine sürekli yayın yapan GNSS uydularının mevcut sinyallerini kullanmasıdır. Bu durum yöntemin sabit GNSS istasyonlarına, tarımsal makinelere veya gelecekte düşük maliyetli alıcılara entegre edilebilmesine olanak verebilir. Çalışmada kullanılan alıcı ise araştırma sınıfında Topcon HiPer Pro'dur; smartphone veya düşük maliyetli cihaz performansı bu deneyde test edilmemiştir.

Tarımsal izleme problemi

Bitki yüksekliği, kanopi yoğunluğu ve toprak nemi farklı tarımsal süreçleri açıklayan temel biyofiziksel değişkenlerdir. Kanopi geometrisi ışık yakalama, transpirasyon ve atmosfer–bitki etkileşimleriyle; toprak nemi ise çimlenme, besin alımı, mikrobiyal süreçler ve su stresiyle ilişkilidir. Bu değişkenler aynı tarlada zaman ve mekân içinde değişebildiğinden seyrek nokta ölçümleri bütün alanın durumunu temsil etmekte yetersiz kalabilir.

Optik uydu gözlemleri geniş alanı düzenli olarak görüntüleyebilmesine rağmen bulut, aerosol ve aydınlatma koşullarından etkilenir ve üç boyutlu kanopi yapısını doğrudan ölçmez. LiDAR kanopi geometrisini yüksek doğrulukla belirleyebilir fakat sürekli yüksek zaman çözünürlüklü uydu gözlemi olarak yaygın değildir. Radar ve pasif mikrodalga sistemleri nem ve bitki su içeriğine duyarlı olmakla birlikte alan ölçeğinde zaman/mekân çözünürlüğü sınırlamaları bulunabilir.

Çalışmanın ele aldığı gözlemsel boşluk, bitki yapısı ile nem kaynaklı dielektrik değişimlere aynı anda duyarlı ve yüksek zamansal sıklıkta çalışabilecek saha ölçekli bir yöntemin araştırılmasıdır.

GNSS-IR Tarımsal İzlemede Nasıl Çalışır?

GNSS-IR, GNSS uydusundan doğrudan antene ulaşan sinyal ile toprak ve bitki örtüsünden yansıyıp daha uzun bir yol izleyerek aynı antene ulaşan sinyali birlikte ölçer; iki yol arasındaki faz farkı SNR'de düzenli girişim salınımları oluşturur ve bu salınımların frekansı, genliği ve fazı sırasıyla yüzey geometrisi, sinyal sönümü ve dielektrik özellikler hakkında bilgi taşıyabilir.

Doğrudan ve yansıyan sinyalin girişimi

Kaynak, GNSS alıcısının gözlediği sinyal-gürültü oranını Denklem (1a) ile ifade eder:

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

\(\theta\) uydu yükseklik açısıdır. \(\widetilde{S}\), alıcıya ulaşan toplam kompleks sinyali; \(\widetilde{N}\) ise gürültü bileşenini temsil eder.

Toplam alınan sinyal gücünde doğrudan ve yansıyan alanların girişimi kaynakta şöyle açılır:

\[ |\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\) doğrudan sinyal genliği, \(A_r\) yansıyan sinyal genliğidir. Son terimde bulunan \(\cos(\Delta\phi)\), SNR'nin uydu hareketi sırasında neden periyodik olarak yükselip düştüğünü gösterir.

Faz farkı kaynakta:

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

şeklinde verilir. \(\lambda\) taşıyıcı dalga boyu, \(h_a\) antenin yerden yüksekliği, \(h_v\) ise bitki örtüsünün etkili faz merkezi yüksekliğidir. \(R_1-R_2\) doğrudan ve yansıyan yollar arasındaki geometrik yol farkını ifade eder.

Denklemin temel sonucu şudur: anten ile etkili yansıtıcı faz merkezi arasındaki yükseklik farkı büyüdükçe uydu yükseklik açısına bağlı faz değişimi de değişir. Bu nedenle SNR salınım frekansı geometrik yükseklik bilgisine dönüştürülebilir.

Neden yalnız 5°–30° yükseklik açıları kullanıldı?

Uydu gökyüzünde yükseldikçe doğrudan ve yansıyan sinyal yollarının geometrik farkı küçülür. Bu durumda multipath girişim salınımlarının genliği ve ayrıştırılabilirliği azalır. Kaynak, 5°–30° aralığında girişim imzasının daha belirgin olduğunu göstermiş ve analiz dışındaki SNR değerlerini maskelemiştir.

Düşük yükseklik açılarında yansıyan sinyalin alıcıya ulaşmasını sağlayan yüzey alanı birinci Fresnel bölgesiyle yaklaşık olarak temsil edilir. Kaynak yarı eksenleri ve yansıma noktasının antene uzaklığını:

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

ilişkileriyle vermektedir. \(a\) ve \(b\) Fresnel elipsinin yaklaşık yarı eksenlerini, \(d\) yansıma noktasının antene yatay uzaklığını ifade eder.

Çalışmadaki bütün uygun uydu izleri birlikte değerlendirildiğinde GNSS-IR'nin her saha kampanyasında yaklaşık 5000 m²'lik tarla alanından yansıma bilgisi aldığı hesaplanmıştır. Bu footprint noktasal toprak probu veya tek cetvel ölçümünden çok daha geniş bir alanı temsil eder; fakat alan içindeki heterojenlik de ölçüm belirsizliğinin bir parçası haline gelir.

Saha ve fenolojik bağlam

Deney alanı Pogliola, Piedmont, kuzeybatı İtalya'da bulunan 2 hektarlık kalıcı bir çok türlü çayırdır. Vejetasyon çoğunlukla Gramineae ve Leguminosae türlerinden oluşmaktadır. 2025 yılında tarla 20 Mayıs, 8 Temmuz ve 13 Ağustos tarihlerinde biçilmiş ve 15 Eylül'de otlatılmıştır.

Sentinel-2 neden kullanıldı?

Sentinel-2 verisi GNSS-IR modelinin girdisi değildir. Çalışma, GNSS-IR salınım genliğinin mevsimsel kanopi yoğunluğu değişimine duyarlı olup olmadığını bağımsız bir gözlemle değerlendirmek için NDVI kullanmıştır.

1 Nisan–1 Ekim döneminden toplam 90 Sentinel-2 Level-2A Harmonized görüntüsü Google Earth Engine'de işlenmiş, bulut/cirrus/gölge pikselleri Scene Classification Layer ile çıkarılmıştır. Geçerli her tarihte 10 m NIR ve kırmızı bantlardan NDVI hesaplanmış ve dört GNSS kampanyası için alan medyanı alınmıştır.

NDVI bu çalışmada biyokütlenin veya Leaf Area Index'in doğrudan nicel ölçümü değildir. Kaynak, NDVI'yi yalnız yeşil bitki örtüsü ve kanopi yoğunluğundaki mevsimsel değişimleri temsil eden göreli bir proxy olarak yorumlamıştır.

Yer referansları

Toprak nemi için Drill & Drop Bluetooth kapasitif sensör 10 cm derinliğe gömülmüş ve hacimsel su içeriği saatlik olarak kaydedilmiştir. GNSS kampanya günlerinde 08:00–14:00 arasındaki SMC değerleri ortalanmıştır.

Kanopi yüksekliği \(h_v\) ise GNSS anteni çevresindeki yedi noktada yere dik tutulan cetvelle ölçülmüştür. Yer ölçümlerinin ortalaması ve standart sapması hesaplanmıştır.

SNR'den dSNR'ye

Ham SNR zaman serisi yalnız yansıma sinyalini içermez. Antenin gain pattern'i ve atmosferik etkiler nedeniyle doğrudan sinyal uydu yükseldikçe yumuşak bir eğilim oluşturur. Kaynak bu düşük frekanslı deterministik eğilimi \(\sin\theta\) bağımsız değişkenine göre ikinci dereceden bir modelle çıkararak ortalaması sıfıra yakın detrended SNR (dSNR) serisi oluşturmuştur.

Bu işlemden sonra dSNR içindeki düzenli salınımlar ağırlıklı olarak yüzeyden yansıyan koherent bileşenin geometrik ve elektromanyetik etkilerini taşır.

GNSS-IR Bitki Boyunu Nasıl Tahmin Eder?

GNSS-IR bitki boyunu doğrudan görsel olarak ölçmez; dSNR salınımlarının \(\sin\theta\)'ya göre dominant frekansını Lomb–Scargle periodogramıyla belirler, bu frekansı anten ile etkili yansıtıcı faz merkezi arasındaki reflektör yüksekliğine dönüştürür ve reflektör yüksekliğini anten yüksekliğinden çıkararak kanopi faz merkezi yüksekliği \(h_v\)'yi elde eder.

Lomb–Scargle periodogramı

Kaynak, dSNR serisinin \(\sin\theta\)'ya göre harmonik içeriğini Lomb–Scargle periodogramıyla analiz etmektedir. Kullanılan temel matris biçimi:

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

ve:

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

şeklindedir. \(y\), dSNR gözlemlerini; \(n\), zaman serisindeki gözlem sayısını; \(\omega=2\pi f\) ise açısal spektral değişkeni ifade eder.

Buradaki \(f\), zamana göre Hertz cinsinden frekans değildir. Kaynak bunu \(\sin\theta\) birimi başına salınım çevrimi olarak tanımlamaktadır. Arama aralığı 0.1–50 cycles per unit of \(\sin\theta\), adım ise 0.01'dir.

Dominant frekanstan kanopi yüksekliğine

Kaynak Denklem (4a)'yı basılı haliyle şöyle vermektedir:

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

Kaynak notasyon notu: Aynı bölümde \(f\), \(\sin\theta\) birimi başına salınım frekansı olarak tanımlanmasına rağmen Denklem (4a)'daki türev değişkeni kaynak PDF'de \(t\) olarak basılmıştır. Verianla bu ifadeyi kaynakta yazılmayan başka bir değişkenle sessizce değiştirmemektedir.

Frekans elde edildiğinde reflektör yüksekliği:

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

ve bitki faz merkezi yüksekliği:

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

olarak hesaplanır.

\(h_a\), antenin yerden yüksekliği olan 2 m'dir. \(RH\), anten ile etkili yansıtıcı faz merkezi arasındaki mesafeyi; \(h_v\) ise bu faz merkezinin yer üzerindeki yüksekliğini temsil eder.

Kalite kontrolü

Filtreleme sonucunda başlangıçta kampanya başına yaklaşık 100 iz olmasına rağmen toplam uydu izlerinin %8'inden azı nihai modellemeye kalmıştır. Bu durum kalan verinin kalite ve fiziksel yorumlanabilirliğini artırırken karmaşık saçılma koşullarının temsili açısından seçim etkisi oluşturabilir.

Doğrusal olmayan dSNR modeli

Dominant geometrik frekans ve \(h_v\) belirlendikten sonra kaynak, salınım genliği ile dielektrik fazı ayırabilmek için ikinci aşamada doğrusal olmayan bir harmonik model kullanmıştır:

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

Modelin uygulanan biçimi:

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

şeklindedir.

\(A_0\) temel genlik ölçeği, \(m_0\) yükseklik açısına bağlı genlik sönümünü, \(\phi_0\) ise yansıtıcı ortamın dielektrik davranışıyla ilişkili yaklaşık sabit faz kaymasını temsil eder. \(h_v\) bu aşamada serbestçe yeniden optimize edilmemiş; önceki Lomb–Scargle frekans analizinden türetilen geometrik bilgi sabit tutulmuştur.

Doğrusal olmayan fit R ortamında `minpack.lm` paketindeki `nlsLM` fonksiyonuyla, bounded least-squares Levenberg–Marquardt algoritması kullanılarak yapılmıştır. Başlangıç genliği gözlenen maksimum dSNR değerine, \(m_0=-0.1\)'e, \(\phi_0=0\)'a ayarlanmıştır. Sınırlar \(A>0\), \(m\le0\) ve \(\phi\in[-\pi,\pi]\) biçimindedir.

Yansıyan alanın fiziksel modeli

Kaynak koherent ve inkoherent yansımayı:

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

ile açıklamaktadır.

\(\Gamma\), toprak–bitki ortamının etkili kompleks permitivitesine \(\widetilde{\varepsilon}\) bağlı Fresnel katsayısıdır. Rastgele yönlenmiş yaprak, sap, boşluk ve nem heterojenliklerinin oluşturduğu \(\widetilde{S}_{incoh}\) bileşeni deterministik salınım üretmek yerine gürültü tabanını yükseltir.

Bu nedenle dSNR'nin anlamlı harmonik bileşeni, yansıyan alanın tamamen değil, yeterince koherent kalan bölümünün doğrudan sinyalle girişiminin gözlenebilir sonucudur.

Faz merkezi neden gerçek bitki tepesinden daha aşağıda olabilir?

Bir çayır veya tarla kanopisi mikrodalga açısından tek bir sert yüzey değildir. Yaprak ve saplar farklı yüksekliklerde saçılma oluşturur. Kaynak bu dağıtık yapıyı:

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

ile ifade eder.

\(F_v(z,\theta)\), kanopinin düşey reflektivite profilidir; \(z_0\) zemin yüksekliği, \(h_{top}\) fiziksel kanopi yüksekliği ve \(k=4\pi/\lambda\) dalga sayısıdır.

Bu dağılımın faz merkezi:

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

olarak tanımlanır.

Denklem (7), GNSS-IR'nin neden cetvelle ölçülen en üst yaprak seviyesini zorunlu olarak vermediğini açıklar. Derindeki bitki katmanlarında sönüm ve saçılma arttığından koherent yansımanın ağırlık merkezi fiziksel kanopi tepesinin altında bulunabilir. Çalışmada gözlenen negatif yükseklik bias'ı bu faz-merkezi davranışıyla ilişkilendirilmiştir.

Çalışmanın Yöntemi ve Bulguları

Veri toplama ve doğrulama zinciri

Kanopi yüksekliği sonucu

Yükseklik tahminleri özellikle Haziran–Temmuz döneminde daha yüksek, Ağustos–Eylül döneminde ise biçim sonrası daha düşük kanopi koşullarını izlemiştir. Kaynak, Haziran–Temmuz döneminde vejetasyon yüksekliğinin yaklaşık 50 cm civarında olduğunu; sonraki biçim dönemlerinin hem yükseklik hem NDVI değerlerini düşürdüğünü tartışmaktadır.

Kanopi yoğunluğu ve dSNR genliği

Sentinel-2 NDVI ile θ = 15° için hesaplanan dSNR model genliği arasında:

r = −0.81

ilişkisi bulunmuştur.

Negatif işaret, bitki örtüsü yoğunlaştıkça dSNR salınım genliğinin azaldığını gösterir. Kaynağın fiziksel yorumu, artan yaprak alanı, biyokütle ve su içeriğinin L-band koherent yansımasını zayıflatmasıdır.

Bu bulgu NDVI'nin GNSS-IR tarafından tahmin edildiği anlamına gelmez. NDVI yalnız GNSS-IR genlik parametresinin bitki örtüsü gelişimine duyarlı olup olmadığını değerlendiren bağımsız bir proxy'dir.

Toprak nemi ve dielektrik faz

Yer ölçümü soil moisture content (SMC) ile \(|\phi_{die}|\) arasında:

r = 0.86

pozitif ilişki bulunmuştur.

Artan toprak su içeriği toprak–bitki ortamının kompleks permitivitesini değiştirdiği için Fresnel yansıma katsayısının fazı da değişebilir. Kaynak bu nedenle \(\phi_0\) parametresini nem kaynaklı dielektrik değişimlere duyarlı bir gözlenebilir olarak yorumlamaktadır.

GNSS-IR Bu Çalışmada Bitki Yapısı ile Toprak Nemini Tam Olarak Birbirinden Ayırabiliyor mu?

Hayır. Çalışma kanopi yüksekliği, dSNR genlik sönümü ve dielektrik fazın sırasıyla yapı, kanopi yoğunluğu ve nem değişimlerine farklı duyarlılıklar gösterdiğini desteklemektedir; ancak yansıyan GNSS sinyali toprak ile bitki örtüsünün birlikte oluşturduğu elektromanyetik ortamdan geldiği için yapısal ve dielektrik katkılar bu veri setinde tam ve kesin biçimde ayrıştırılmamıştır.

Üç GNSS-IR gözlenebilirinin tamamlayıcı anlamı

Bu tablo üç parametrenin tamamen bağımsız fiziksel ölçümler olduğu anlamına gelmez. Aynı yansıma sürecinde geometri, kanopi yoğunluğu, su içeriği ve saçılma yapısı birbirini etkileyebilir. Kaynağın sonucu, farklı GNSS-IR parametrelerinin farklı fiziksel bileşenlere baskın duyarlılık gösterdiğidir; tam separation henüz sağlanmamıştır.

Toprak nemi sonuçlarında yöntem–sonuç tutarsızlığı

Yöntem bölümü NDVI ve SMC'yi doğrudan retrieval hedefi olarak tanımlamaz ve bu iki referansla ilişkinin yalnız Pearson korelasyonu üzerinden değerlendirileceğini belirtir. Buna rağmen sonuç bölümünde SMC için ayrıca “ME = −5.3” ve “MAPE = %21” yazılmıştır.

Kaynak, \(\phi_{die}\) ile SMC arasında nasıl doğrudan bir SMC tahmini üretildiğini ve bu iki hata metriğinin hangi hesaplama zincirinden geldiğini yöntem bölümünde açıklamamaktadır. Bu nedenle Verianla, bu iki sayıyı doğrulanmış bir “GNSS-IR soil-moisture retrieval error” metriği olarak kullanmamaktadır. Güvenilir biçimde desteklenen nicel ilişki Pearson \(r=0.86\)'dır.

Korelasyonların istatistiksel sınırı

Bütün temel korelasyonlar yalnız dört saha kampanyasının toplulaştırılmış değerlerine dayanır: \(n=4\). Kaynak bu nedenle korelasyonların katı istatistiksel genellemeler olarak değil, fiziksel ilişkilerin göstergesi olarak yorumlanması gerektiğini açıkça belirtmektedir.

Dört tarih rastgele seçilmiş aynı koşulların tekrarları değildir; biçme öncesi ve biçme sonrası farklı fenolojik rejimleri temsil etmektedir. Bu çeşitlilik yöntemin farklı saha durumlarına tepki verdiğini göstermeye yardımcı olsa da örneklem sayısını büyütmez.

Filtreleme ve selection bias

Kalite filtreleri, SNR içinde fiziksel olarak yorumlanabilir koherent multipath salınımı bulunan izleri seçmektedir. Bu yaklaşım gürültülü ve kararsız çözümleri azaltır.

Bunun ters tarafı, heterojen kanopi, karmaşık yansıma geometrisi veya daha hızlı değişen nem koşullarının model varsayımlarına daha az uyan uydu izlerinin veri setinden daha sık çıkarılabilmesidir. Yazarlar bu nedenle nihai izlerin bütün saha koşullarının tarafsız örneklemi değil, modelin en açık çalıştığı gözlemler olduğunu belirtmektedir.

Lomb–Scargle yüksekliğindeki negatif bias

Lomb–Scargle periodogramı eşit aralıklı olmayan SNR örnekleminde dominant salınım frekansını bulmak için güçlü bir araçtır; ancak sonlu uydu açısı penceresi spectral leakage oluşturabilir. Ek yansıma yolları da ana periodogram tepesini etkileyebilir.

Kaynak bu etkilerin dominant frekansı bir miktar aşağı kaydırabileceğini ve bunun reflektör yüksekliği kestirimine bias olarak yansıyabileceğini tartışmaktadır. Bunun yanında asıl fiziksel nedenlerden biri kanopinin bir hacim saçıcı olmasıdır: faz merkezi fiziksel üst yüzeyden aşağıdadır.

Çalışmanın desteklediği sonuçlar

Çalışma, dört farklı fenolojik saha koşulunda GNSS-IR dSNR sinyalinin kanopi yapısı ve nemle ilişkili sistematik değişimler taşıdığını desteklemektedir.

Çalışma, Lomb–Scargle tabanlı GNSS-IR reflektör yüksekliğinin çayır kanopisindeki yükseklik değişimini izleyebildiğini ve cetvel referansıyla güçlü ilişki gösterdiğini desteklemektedir.

Çalışma, dSNR genlik sönümünün mevsimsel kanopi yoğunluğu/NDVI değişimine ve dielektrik faz parametresinin saha nem değişimlerine duyarlı olduğunu desteklemektedir.

Çalışma, GNSS-IR'nin tek bir gözlem altyapısından birden fazla tarımsal biyofiziksel bilgi kanalı çıkarma potansiyeline sahip olduğunu desteklemektedir.

Çalışmanın desteklemediği sonuçlar

Çalışma GNSS-IR'nin bütün ürün türlerinde %18 kanopi yüksekliği hatasıyla çalışacağını göstermemektedir.

Çalışma \(\phi_{die}\)'nin bitki örtüsü etkisinden tamamen arındırılmış saf bir toprak nemi ölçümü olduğunu göstermemektedir.

Çalışma tarla sulamasının otomatik olarak GNSS-IR ile kontrol edildiğini veya su tüketiminin azaltıldığını göstermemektedir.

Çalışma düşük maliyetli smartphone GNSS sistemlerinin aynı performansı verdiğini göstermemektedir; smartphone kullanımı gelecek çalışma yönüdür.

Çalışma mısır, buğday, soya veya başka tarla bitkilerinde aynı ilişkilerin doğrulandığını göstermemektedir.

Başlıca sınırlılıklar

Gelecek çalışma yönleri

Kaynak daha sık zaman örneklemesi, daha fazla saha, farklı ürün mimarileri ve heterojen koşullarda doğrulama önermektedir. Özellikle maize, wheat ve soybean gelecekte test edilmesi planlanan ürünler arasında belirtilmektedir.

Daha esnek multi-frequency veya time-varying harmonic modeller, robust regression yaklaşımları ve kanopi heterojenliğini doğrudan ele alan fiziksel scattering modelleri daha fazla GNSS izini kaybetmeden parametre üretmek için araştırılabilir.

Düşük maliyetli çift frekanslı GNSS alıcıları ve smartphone sensörleri de gelecekteki ölçeklenebilir uygulama yolları arasında gösterilmiştir. Bu potansiyel, mevcut çalışmanın test edilmiş sonucu değil, yazarların gelecek geliştirme yönüdür.

Kaynak ve Yöntem Notu

Özgün başlık: Monitoring Crop Structure and Moisture Using GNSS Interferometric Reflectometry Based on SNR Modeling

Yazarlar: Samuele De Petris; Enrico Borgogno-Mondino

Yazışma yazarı: Enrico Borgogno-Mondino

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

Dergi: Agronomy

Yayınevi: MDPI

Cilt / Sayı / Makale: 16 / 9 / 922

DOI: 10.3390/agronomy16090922

Tarihler: Alındı: 25 Mart 2026; revize edildi: 27 Nisan 2026; kabul edildi: 29 Nisan 2026; yayımlandı: 1 Mayıs 2026.

Kaynak türü: Hakemli deneysel araştırma makalesi; tarla GNSS ölçümü, mikrodalga interferometrik reflektometri, spektral analiz ve doğrusal olmayan sinyal modelleme.

Sürüm: Version of Record.

Lisans: Creative Commons Attribution (CC BY). PDF'nin görünür lisans metninde sürüm numarası ayrıca belirtilmediğinden lisans sürümü tahmin edilmemiştir.

Ana saha: Pogliola, Piedmont, kuzeybatı İtalya; 2 ha kalıcı çok türlü çayır.

GNSS alıcısı: Topcon HiPer Pro; GPS ve GLONASS; L1/L2; anten yüksekliği 2 m; 1 Hz.

GNSS kampanyaları: 1 Haziran, 4 Temmuz, 2 Ağustos ve 7 Eylül 2025.

Optik yardımcı veri: Copernicus Sentinel-2 Level-2A Harmonized; 90 görüntü; 1 Nisan–1 Ekim 2025; Google Earth Engine.

Toprak nemi referansı: Drill & Drop Bluetooth kapasitif sensör; 10 cm derinlik; saatlik volumetric water content.

Kanopi yüksekliği referansı: GNSS anteni çevresinde yedi noktada manuel cetvel ölçümü.

Temel sinyal işleme: Ascending/descending uydu yayı ayrımı; 5°–30° elevation mask; quadratic detrending; Lomb–Scargle periodogramı; nonlinear harmonic regression.

LSP frekans aralığı: 0.1–50 cycles per unit of sin θ; adım = 0.01.

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

Doğrusal olmayan fit: R `nlsLM`, `minpack.lm`; bounded least-squares Levenberg–Marquardt.

Ana sonuçlar: Kanopi yüksekliği r = 0.89, ME = −0.067 m, MAPE = %18; NDVI ile dSNR genliği r = −0.81; SMC ile mutlak dielektrik faz r = 0.86.

İstatistiksel örneklem sınırı: Temel karşılaştırmalar dört kampanyanın toplulaştırılmış değerlerine dayanır, n = 4.

Finansman: Araştırma dış finansman almamıştır.

Veri erişilebilirliği: Verilerin talep üzerine sağlanacağı bildirilmiştir.

Ek materyal: Table S1, her kalite kontrol aşamasından sonra kalan uydu izi sayılarını içermektedir.

Çıkar çatışması: Yazarlar çıkar çatışması olmadığını bildirmektedir.

Yazar katkıları: Samuele De Petris kavramsallaştırma, yöntem, yazılım, doğrulama, formal analiz, veri kürasyonu, özgün taslak ve görselleştirme rollerini yürütmüştür. Araştırma ve kaynaklar iki yazar tarafından paylaşılmıştır. Enrico Borgogno-Mondino danışmanlık, proje yönetimi ve finansman edinimi rollerini üstlenmiş; her iki yazar gözden geçirme ve düzenlemeye katkıda bulunmuştur.

Kaynak içi notasyon uyarısı: Metin Lomb–Scargle spektral değişkenini sin θ başına salınım çevrimi olarak tanımlamasına rağmen Denklem (4a)'da türev değişkeni t olarak basılmıştır. Bu ifade Verianla tarafından sessizce değiştirilmemiştir.

Kaynak içi yöntem/sonuç uyarısı: Yöntem bölümü soil moisture ve NDVI karşılaştırmalarında yalnız Pearson r kullanılacağını belirtirken sonuç bölümünde soil moisture için ayrıca ME = −5.3 ve MAPE = %21 raporlamaktadır. Bu hata metriklerinin türetilme yöntemi açıklanmadığından bağımsız nem retrieval performansı olarak kullanılmamıştır.

Temel yöntemsel sınırlılık: Tek saha, tek kanopi tipi, dört kampanya, seyrek yer referansları ve yoğun kalite filtreleri nedeniyle sonuçlar operasyonel çok ürünlü tarım sistemi için nihai kalibrasyon değildir.

Ana bilimsel yorum sınırı: Çalışma GNSS-IR gözlenebilirlerinin bitki yapısı ve dielektrik/nem değişimlerine farklı duyarlılık gösterdiğini destekler; bu iki fiziksel katkının tam olarak birbirinden ayrıldığını göstermez.