Bu çalışma coğrafi veri, uzaktan algılama ve saha doğrulama verilerini içermekte olup; insan katılımcıları, tanımlanabilir kişisel verileri, hayvanları veya omurgalı dokularını kapsamamaktadır; bu nedenle, kurumsal insan veya hayvan etiği onayı gerekmemiştir.
Çalışma bölgesi
Panjkora nehir havzası, Pakistan'ın kuzeyindeki doğu Hindu Kuş Dağları'nda yer alan önemli bir fizyografik bölgedir (Şekil 1). Panjkora nehri, havzanın ana nehri olup (113 km uzunluğunda ve 5758,27 km2 havza alanına sahiptir) Hindu Kuş'un buzlarla kaplı dağlarından bir akarsu olarak başlar. Dir Lower'daki Chakdara yakınlarında Swat Nehri ile birleşir20. Panjkora Nehri'ne Barawal, Dir, Gawaldai, Jandol ve Kohistan dahil olmak üzere beş önemli akarsu veya dere katılır. Nehir, 34°39′30′′ ile 35°46′1′′ Kuzey enlemleri ve 71°13′08′′ ile 72°22′13′′ Doğu boylamları arasında uzanır. Bölgenin konumu ve engebeli topografyası, iklimini (dağlık ve ılıman) önemli ölçüde etkilemektedir. Havzanın Üst (Kumrat, Thal) bölgesinde kış mevsimi daha uzun, yaz mevsimi ise daha soğuk geçer. Kasım ayından itibaren sıcaklıklar hızla düşer. Bununla birlikte, Dir Lower'da (Timergara, Talaash, Maidan, Samarbagh) sıcaklıklar genellikle Aralık'tan Şubat'a kadar donma noktasının üzerindedir. Timergara'da en sıcak aylar, ortalama maksimum sıcaklıkların 35 °C'nin üzerinde olduğu Haziran ile Ağustos arasıdır; Dir Kasabası'nda ise en sıcak aylar Haziran ve Temmuz'dur (maksimum sıcaklıklar sırasıyla 32,4 °C ve 31,5 °C). Musonlar yaz yağışlarının kaynağıyken, Batı Alçak Basınç sistemi kış yağışlarını getirir. Çalışma alanı, yıl boyunca yüksek bağıl nem ile karakterize edilir. Özellikle Wari'nin yukarı ve aşağı akış bölgelerinde (nerdeyse) her yıl nehir taşkınları ve ani seller21 meydana gelir. Bölgede yetiştirilen başlıca ürünler arasında pirinç, buğday, mısır, patates ve soğan yer alırken, çalışma alanında yetiştirilen önemli meyveler hurma, portakal, elma, ceviz, kayısı, erik, malta eriği ve duttur.

Şekil 1: Kuzey Pakistan'daki Panjkora Nehri Havzası'nın çalışma alanı haritası. (A) Khyber Pakhtunkhwa'nın Pakistan içindeki konumu; (B) Panjkora Nehri Havzası'nın Khyber Pakhtunkhwa içindeki konumu; ve (C) havza sınırını, yükselti dağılımını, nehir ağını ve çalışma alanındaki önemli konumları gösteren Panjkora Nehri Havzası. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.
Veri toplama ve hazırlama
Bu çalışma için veriler farklı kaynaklardan toplanmıştır. Yağış verileri, 2014 ile 2023 yılları arasını kapsayacak şekilde Küresel Yağış Ölçümü (GPM) Ulusal Havacılık ve Uzay Dairesi'nden (NASA) https://gpm.nasa.gov/missions/GPM adresinden indirilmiştir. Toprak dokusu verileri, Pakistan Khyber Pakhtunkhwa Toprak Etüt Müdürlüğü'nden (www.soilconservation.kp.org) toplanmıştır. Jeolojik veriler, Pakistan Jeolojik Etüt bölgesel ofisinden (https://gsp.gov.pk/) temin edilmiştir. Arazi Senaryolarının (Arazi kullanımı/Arazi Örtüsü) toplanması ve hesaplanması için Sentinel 2 görüntüleri Avrupa Uzay Ajansı (ESA) Copernicus Open Access Hub'dan (https://scihub.copernicus.eu/) alınmıştır. Arazi kullanımı/arazi örtüsü (LULC) haritalaması için 10 Eylül 2025 tarihinde elde edilen Sentinel-2B görüntüleri kullanılmıştır. Görüntü, Maksimum Olabilirlik Sınıflandırması (MLC) algoritması kullanılarak işlenmiş ve sınıflandırılmıştır. Çalışma alanı genelinde; su kütleleri, orman, tarım arazisi, kentsel alanlar, çıplak toprak, kar/buz ve mera olmak üzere yedi LULC sınıfını temsil eden toplam 65 eğitim örneği toplanmıştır. Hazırlanan eğitim örnekleri, denetimli sınıflandırma yapmak ve nihai LULC haritasını oluşturmak için kullanılmıştır. Sınıflandırma doğruluğu, genel doğruluk ve Kappa katsayısını içeren, doğrulama örneklerine dayalı bir doğruluk değerlendirme yaklaşımı kullanılarak analiz edilmiştir. 12.5 uzamsal çözünürlüğe sahip dijital yükseklik modeli (DEM), 2/12/2023 tarihinde Alaska Uydu Tesisi'nden (ASF) (https://asf.alaska.edu/) temin edilmiştir. DEM modeli daha sonra eğim, drenaj ağı, drenaj yoğunluğu ve yükseklik katmanlarını oluşturmak için kullanılmıştır. Mevcut yağmur suyu hasadı yapılarına ait veriler, çapraz doğrulama amacıyla ilgili departmanlardan toplanmıştır.
Tüm mekansal veri setleri, coğrafi bilgi sistemi (CBS) yazılımı kullanılarak işlenmiş ve analiz edilmiştir (bakınız Materyaller Tablosu). CBS tematik katman verileri Ek Dosya 1'de sunulmuştur. Mekansal tutarlılığı sağlamak ve doğru alan hesaplamaları yapmak için tüm girdi veri setleri ortak bir projeksiyon koordinat referans sistemine (CRS) (WGS 1984 UTM Zon 42N) projekte edilmiştir. Farklı mekansal çözünürlükteki raster veri setleri, kategorik veri setlerinin orijinal mekansal özellikleri korunarak, en yakın komşu yeniden örnekleme yöntemiyle ortak bir ızgaraya göre yeniden örneklenmiş ve hizalanmıştır. Mekansal hizalama için 12.5 m mekansal çözünürlüğe sahip Sayısal Yükseklik Modeli (DEM) referans raster olarak kullanılmış ve tüm tematik katmanlar aynı hücre boyutu ve kapsama sahip raster formatına dönüştürülmüştür. Tüm girdi katmanlarını dışarı aktarmak ve analiz için tutarlı bir mekansal kapsam sağlamak amacıyla Panjkora Nehri Havzası'nın çalışma alanı sınırı maske olarak kullanılmıştır. Eksik pikseller ve havza sınırı dışındaki alanlar analizden çıkarılmış ve NoData değerleri olarak kabul edilmiştir. Tematik katmanlar (yağış, eğim, drenaj yoğunluğu, lineament yoğunluğu, toprak, jeoloji ve arazi kullanımı/arazi örtüsü), Jenks Doğal Kırılmalar sınıflandırma yöntemi kullanılarak uygunluk sınıflarına yeniden sınıflandırılmış ve MIF ve AHP yaklaşımlarına göre ilgili rütbe/ağırlıklar atanmıştır. Veri kaynakları Tablo 1'de gösterilmektedir.
Tablo 1: Yağmur suyu hasadı uygunluk değerlendirmesi için kullanılan veri kaynakları ve özellikleri. Bu tabloyu indirmek için lütfen buraya tıklayın.
MIF uygunluk modellemesi
Başlangıçta, çeşitli parametrelerin seçimi literatür taramasına dayanarak gerçekleştirilmiştir12. Yağmur Suyu Hasadı (RWH) için uygun konumları belirlemek amacıyla; yağış, litoloji, lineament yoğunluğu, drenaj yoğunluğu, toprak dokusu, eğim ve arazi kullanımı/arazi örtüsü ayrı etkileyici faktörler olarak dikkate alınmıştır. Bu amaç doğrultusunda, parametrelerin etki ölçeğini oluşturmak için parametrelerin ön işlemesi yapılmış; ardından veriler RWH açısından önem derecelerine göre kategorize edilmiş ve çoklu etkileyici faktör formülü (Denklem 1) kullanılarak ana ve ikincil önem dereceleri belirlenmiştir. Tablo 2, farklı faktörlerin ana ve ikincil önem derecelerini göstermektedir22 (Bkz. Ek Dosya 2)
Tablo 2: Çoklu Etki Faktörü (MIF) modelinde kullanılan seçilmiş etkileyici faktörler ile bunların ana ve ikincil etki puanları. Bu tabloyu indirmek için lütfen buraya tıklayın.
Seçilen faktörler, şu ilişki kullanılarak sıralanmıştır,
[(X+Y) ÷ ∑(X+Y)] × 100 (1)
burada Y faktörlerin ikincil etkisini, X ise ana etkisini temsil eder. Her faktörün ana ve ikincil etkileri Denklem 1 kullanılarak hesaplanır.
Majör (X) ve minör (Y) etki puanları, önceki çalışmalara ve her bir faktörün akış üretimi, sızma ve yağmur suyu hasadı potansiyelini kontrol etmedeki hidrolojik önemine göre atanmıştır12. RWH uygunluğu üzerinde doğrudan etkisi olan faktörlere majör etki atanırken, minör etki, kontrol parametreleri arasındaki dolaylı ilişkileri temsil etmiştir. Faktör ağırlıkları, birleşik majör ve minör etki puanlarının normalize edilmesiyle Denklem (1) kullanılarak hesaplanmıştır. Alt sınıf ağırlıkları; akış birikimine, sızma kapasitesine, su tutma kapasitesine ve RWH yapılarına olan uygunluğa yaptıkları göreceli katkıya göre atanmıştır. Bu yaklaşım, CBS tabanlı uygunluk analizi için şeffaf ve tekrarlanabilir bir ağırlıklandırma çerçevesi sağlamıştır.
Saaty ölçeğine dayalı göreceli önem Tablo 3'te gösterilmiştir.
Tablo 3: Analitik Hiyerarşi Süreci (AHP) analizi için kullanılan Saaty'nin göreceli önem ölçeği. Bu tabloyu indirmek için lütfen buraya tıklayın.
Tüm parametrelerin tematik seviye puanları birleştirilmiş olup, MIF parametrelerinin her bir alt sınıf puanı Tablo 4'te listelenmiştir. Yeniden sınıflandırma tekniği kullanılarak, MIF çıktısı yağmur suyu hasadı için beş kategoriye ayrılmıştır. Son olarak, setler, çiftlik göletleri, şev dolguları ve diğer koruma odaklı yapılar gibi farklı YSH yapılarının kurulması için önerilen nihai konumların haritaları oluşturulmuş ve analiz edilmiştir. Şekil 2 metodoloji çerçevesini göstermektedir.

Şekil 2: CBS tabanlı MIF ve AHP modelleri kullanılarak yağmur suyu hasadı saha uygunluk değerlendirmesi için metodolojik çerçeve. Çerçeve; jeoloji, toprak, eğim, drenaj yoğunluğu, lineament yoğunluğu, arazi kullanımı/arazi örtüsü (LULC) ve yağış gibi tematik katmanlar elde etmek için saha araştırması, jeoloji, toprak, ALOS PALSAR DEM, ESA ve GPM verilerinin toplama ve işleme sürecini göstermektedir. Bu katmanlar, RWH uygunluk haritasını oluşturmak için Çoklu Etkileyici Faktör (MIF) yaklaşımı kullanılarak entegre edilmiş ve ardından nihai doğrulanmış haritaları üretmek için saha tabanlı doğrulama yapılmıştır Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.
Tablo 4: Yağmur suyu hasadı uygunluk haritalaması için tematik faktörlerin ve alt sınıfların Çoklu Etki Faktörü (MIF) tabanlı sıralamaları ve ağırlıkları. Bu Tabloyu indirmek için lütfen buraya tıklayın.
AHP uygunluk modellemesi
Analitik Hiyerarşi Süreci (AHP), karar vericinin öncelikleri belirlemesine ve en iyi seçeneği seçmesine yardımcı olan, zor karar verme durumlarını yönetmek için etkili bir tekniktir23. AHP tekniği, matematik ve uzman bilgisinin uygulanması yoluyla karmaşık kararları organize etmek ve değerlendirmek için kullanılan sistematik bir çerçevedir24. AHP, karmaşık yargıları ikili karşılaştırmalar yoluyla basitleştirip ardından sonuçları değerlendirerek, bir kararın hem öznel hem de nesnel yönlerinin belirlenmesine yardımcı olur25. Karşılaştırmalar öznel veya bireysel bakış açılarına dayandığı için kaçınılmaz olarak bazı farklılıklar olacaktır. Tutarlılık oranını hesaplayarak ve karar verme yanlılığını ortadan kaldırarak AHP tekniği, karar vericinin yargılarının tutarlılığını değerlendirmek için yararlı bir araç sunarak algıların tutarlılığını sağlar. AHP'nin temel avantajlarından biri, farklı kriterlerin ikili karşılaştırmaları arasındaki tutarlılık derecesini nicelleştiren tutarlılık oranıdır26,27,28,29. Coğrafi veri girişleri, AHP tarafından birleştirilerek bir karar çıktısına dönüştürülür. Saaty'nin ölçeği (Tablo 3) kullanılarak, ikili karşılaştırma matrisi oluşturulması yoluyla çeşitli temalar ve nitelikler hakkındaki kalitatif veriler kantitatif verilere dönüştürülür30,31. Temel süreç; hedefin belirlenmesini, nihai kararı etkileyen faktörlerin veya standartların değerlendirilmesini ve her bir kritere puan atamak için Saaty'nin ölçeğinin kullanılmasını içerir. Atanan ağırlıkların tutarlılığını kontrol etmek için, Saaty tarafından önerilen23 tutarlılık oranı (CR), Denklem 2 ve 3 kullanılarak hesaplanmıştır:
CR = CI/RCI (2)
burada CI tutarlılık indeksi, RCI ise rastgele tutarlılık indeksidir.
Tutarlılık indeksi (CI) şu denklemle verilir:
(3)
burada n kriter sayısını ve λmax en büyük özdeğeri temsil eder. Tutarlılık indeksinin ortalama değeri, rastgele indeks ile tahmin edilir.
RWH yapı seçimi
Arazi örtüsü ve arazi kullanımı (LULC)
Arazi kullanımı, arazinin nasıl kullanıldığını tanımlarken; arazi örtüsü, arazinin doğal özelliklerini tanımlar. Akış yayılımı hakkında önemli bilgiler LULC'de32 yer almaktadır. Bitki örtüsüyle kaplı alanlarda, daha yüksek emilim ve sızma oranları daha az akışla ilişkilendirilirken; çıplak araziler ve yapılaşmış alanlar yüksek akış oluşumunu teşvik eder33,34. Çalışma alanının arazi kullanımı/arazi örtüsü modellerini hazırlamak için Sentinel 2b uydu verileri kullanılmıştır. Panjkora nehir havzasının arazi kullanımı; su kütleleri, orman, ürün ve tarım arazisi, kentsel alan, çıplak arazi, kar/buz ve mera olmak üzere yedi sınıfa ayrılmıştır. Farklı arazi kullanımı/arazi örtüsü sınıflarına atanan uygunluk ağırlıkları, bunların akış üretimi, sızma ve yağmur suyu depolama potansiyeli üzerindeki etkilerine göre belirlenmiştir. Tarım arazileri, genellikle orta düzeyde akış üretmesi ve sulama için hasat edilen sudan doğrudan faydalanması nedeniyle en yüksek uygunluk puanını almıştır. Çıplak arazilere de nispeten yüksek bir ağırlık verilmiştir çünkü seyrek bitki örtüsü ve açık toprak yüzeyleri, yoğun bitki örtüsü olan alanlara kıyasla daha fazla yüzey akışını teşvik eder. Buna karşılık, ormanlık alanlara daha düşük ağırlıklar verilmiştir çünkü yoğun bitki örtüsü yağışı engeller, geniş kök sistemleri aracılığıyla sızmayı artırır ve yüzey akışını azaltır. Kentsel alanlara ve mevcut su kütlelerine, ya ek RWH yapıları inşa etmek için sınırlı imkanlara sahip oldukları ya da halihazırda geçirimsiz yüzeyler veya mevcut su kütleleri tarafından kaplı oldukları için daha düşük uygunluk atanmıştır (Şekil 3A).
Drenaj yoğunluğu
Bir bölgenin yer altı suyu sızması ve su yüzey akışı, drenaj yoğunluğu ile tanımlanır. Drenaj yoğunluğu, hem yer altı hidrolojik oluşumlarını hem de yüzey özelliklerini yansıtır. Kanal aralıklarının sıkılığını ve yüzey malzemesinin özelliklerini gösterir. Drenaj yoğunluğu azaldıkça yüzey akışı azalır ve tersi durum geçerlidir12. Düşük drenaj yoğunluğuna sahip bölgelerde genellikle daha düşük sızma ve daha düşük yüzey akışı görülür ve tersi durum geçerlidir. Yoğun drenaj ağları, yağmur suyu toplama için temel öneme sahiptir. RWH, suyun akmasına ve bir toplama noktasına hızla iletilmesine olanak tanıyan bir sistem sağladıkları için yüksek drenaj yoğunluğuna sahip bölgelere daha uygundur34,35. Panjkora Nehri havzasının drenaj yoğunluğu, Jenks Doğal Kırılmalar sınıflandırmasına göre beş sınıfa ayrılmıştır: 0–9,4907, 9,4907–27,207, 27,207–48,219, 48,219–79,089 ve 79,089–161,34 km/km2. Düşükten orta düzeyde drenaj yoğunluğuna sahip bölgeler, yağmur suyu hasadı için ideal konumlar olarak kabul edildikleri için daha yüksek bir ağırlık değerine atanmıştır (Şekil 3B).
Lineament yoğunluğu
Lineamentler, tipik olarak jeolojik haritalardan türetilen ve uydu görüntülerinde de görülebilen doğrusal yüzey altı özelliklerdir. (Lokalize veya yapısal ayrışma bölgelerinin altında gömülü olan) lineamentler, artmış gözeneklilik ve geçirgenlik sergiler12. Lineamentler, uzaktan algılama görüntü işleme yazılımı kullanılarak Landsat 8 görüntüsünden çıkarılmıştır. Lineament raster katmanını oluşturmak için çizgi yoğunluğu (line-density) aracı kullanılmıştır. Lineament yoğunluğu, Jenks Doğal Kırılımlar (Natural Breaks) sınıflandırma yöntemi kullanılarak beş sınıfa ayrılmıştır: 0.0072-0.406 km/km2, 0.406-0.664 km/km2, 0.664-0.921 km/km2, 0.921-1.33 km/km2 ve 1.33-2.13 km/km2 (Şekil 3C).
Toprak
Toprak tekstürü, RWH planlaması ve saha seçimi açısından önemli bir faktördür. Toprağın infiltrasyon kapasitesi, tekstürü tarafından belirlenir. Genel olarak, kumlu topraklar killi topraklara kıyasla düşük yüzey akışı oluşturur36. Silt, kum ve kil yüzdeleri toprağın tekstürel sınıfını tanımlar. Killi toprak düşük geçirgenliğe sahiptir ve toplanan suyu tutabilir, bu nedenle yağmur suyu toplama için genellikle orta ve ince taneli toprağa sahip alanlar tercih edilmiştir8,37. Çalışma bölgesi beş toprak tekstürü ile karakterize edilmiştir: Buzullar ve kar örtüleri, tınlı, kireçsiz killi toprak, Tınlı sığ kireçsiz toprak, Tınlı çok sığ toprak ve kaya çıkıntıları (Şekil 3D).
Eğim
Sızma ve yüzeysel akış, topografyadan önemli ölçüde etkilenir8. Havzadaki eğim varyasyonu, sağanak yağış sırasında ve sonrasında suyun nasıl aktığı üzerinde belirgin bir etkiye sahiptir. Dik eğimli alanlarda, ihtiyaç duyulan yoğun toprak işleri nedeniyle yağmur suyu hasadı (RWH) yapılarının inşa edilmesi maliyet etkin değildir38. Yüksek RWH potansiyeli için en uygun konum, hafif eğimli alanlardır. Dik eğimli bölgelerde (yüzde 5'ten büyük eğimler) RWH yapıları dayanıklı değildir39. Dik eğimli alanlarda erozyon kontrol önlemleri de dikkate alınır40. Eğim derece cinsinden hesaplanmış ve çalışma alanı Jenks Natural Breaks sınıflandırması kullanılarak beş sınıfa ayrılmıştır: 0°–11.9°, 12°–22.5°, 22.6°–31.8°, 31.9°–42.4° ve 42.5°–82° (Şekil 3E).
Yağış
Yüzey akışını oluşturan temel bileşen yağışlardır. 2014-2023 yılları arasındaki yağış/presipitasyon verileri, Küresel Yağış Ölçümü (GPM) verileri olarak NASA'nın https://gpm.nasa.gov/missions/GPM adresinden indirilmiştir41. 2014-2023 dönemine ait GPM yağış verileri ve Jenks Doğal Kırılmalar sınıflandırması kullanılarak çalışma alanı beş yağış (mm) sınıfına ayrılmıştır: 49,93–57,014, 57,014–61,773, 61,773–65,262, 65,262–68,646 ve 68,646–76,894 (Şekil 3F).
Jeoloji
Bir havzanın fiziksel bileşimi ve ürettiği toprak miktarı, bölgenin jeolojisinden büyük ölçüde etkilenir. Jeolojik özellikler, suyun yer altı akiferlerine akışını kontrol eder40. Mevcut çalışma alanında bulunan iki ana kayaç türü sedimanter ve metamorfik kayaçlardır. Litoloji genel olarak Alt Paleozoyik kayaçlar, Karbonifer sedimanter kayaçlar, Kretase sedimanter kayaçlar, Mezozoyik intrüzif ve metamorfik kayaçlar, Triyas kayaçlar, ayrılmamış Paleozoyik kayaçlar, ayrılmamış Paleozoyik kayaçlar ve ayrılmamış Prekambriyen kayaçlar ile ayrılmamış Silüryen kayaçlar olarak ayrılmıştır. Kullanılabilirlik ve depolama kapasitesi, litoloji türünden önemli ölçüde etkilenir; bazı kayaçlar yüzey sularını süzme ve akiferi besleme yeteneğine sahiptir41. Öte yandan, bazı kayaçlar suyun geçişine izin vererek yer altı sularının yeniden beslenmesine yardımcı olur. Litoloji; geçirgenlik, gözeneklilik ve infiltrasyon kapasitesi üzerindeki etkileri aracılığıyla yüzeysel akış oluşumunu güçlü bir şekilde kontrol eder. Panjkora Havzası'nda, kompakt metamorfik kayaçlar genellikle konsolide olmamış veya yüksek gözenekli sedimanter çökeltilere göre daha düşük birincil gözeneklilik ve geçirgenlik sergiler. Sonuç olarak, yağışların sızma olasılığı daha düşük, yüzeysel akış oluşturma olasılığı ise daha yüksektir; bu da bu formasyonları yüzey yağmur suyu hasadı yapıları için daha uygun hale getirir. Buna karşılık, iri taneli veya kumlu malzemeler içeren sedimanter formasyonlar genellikle daha fazla infiltrasyona izin vererek depolama için mevcut olan yüzeysel akışı azaltır. Bu nedenle, metamorfik kayaçlara daha yüksek uygunluk ağırlıkları atanırken, sedimanter formasyonlara nispeten daha düşük ağırlıklar atanmıştır. Şekil 3G, Panjkora nehir havzasının jeolojik haritasını göstermektedir. Tüm veriler Ek Dosyalar 1 ve 3'te mevcuttur.

Şekil 3: Panjkora Nehri Havzası'nda yağmur suyu hasadı alanı değerlendirmesi için kullanılan tematik faktörlerin mekansal dağılımı. (A) Arazi kullanımı/arazi örtüsü, (B) drenaj yoğunluğu, (C) lineament yoğunluğu, (D) toprak dokusu, (E) eğim, (F) yağış ve (G) jeoloji. Farklı renkler, her bir tematik faktörün ilgili sınıflarını temsil etmektedir. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayınız.