Nesibe KÖSE, Ünal AKKEMİK, Nüzhet DALFES

1. ANADOLU’NUN İKLİM TARİHİNİN SON 500 YILI: ESKİŞEHİR VE ÇEVRESİ İÇİN BİR DENDROKLİMATOLOJİK RÖKONSTRÜKSİYON Nesibe KÖSE, Ünal AKKEMİK, Nüzhet DALFES

2. Projenin amaçları • Anadolu için iklime duyarlı ana kronolojiler oluşturmak; • Sıcaklık ve yağış ile bu bölgedeki ağaçların yıllık halka gelişimi arasındaki ilişkileri ortaya koymak; • Sıcaklık ve yağışın rekonstrüksiyonunu yapmak; Anadolu’da geçmişte yaşanan kurak ve yağışlı yılları ve bunların sıklığını belirlemek; • NAOI (Kuzey Atlantik Salınımı İndisi) ile bu bölgedeki ağaçların yıllık halka gelişimi arasındaki ilişkiyi ortaya koymak ve yüksek ilişki bulunduğu taktirde NAOI nin rekonstrüksiyonunu yapmak; • İklim parametrelerinin uzayda ve zamanda değişkenliğini belirlemektir.

3. Proje kapsamında Ankara, Eskişehir ve Çorum’da arazi gezileri yapılmış; toplam 17 araştırma alanından örnekler alınmıştır. Örnek alınan ağaçlar çoğunlukla 200-400 yaşlarında olup, en yaşlılarına (500-700 yaşında) Çorum-Kargı’da rastlanmıştır. Buna karşın, önümüzdeki dönem gidilmesi planlanan Muğla-Sandras dağındaki ağaçların ortalama yaşı yaklaşık 750 yıldır.

4. Projenden elde edilen ilk sonuç olarak Ankara ve Eskişehir civarındaki 4 yöreden alınan örneklere dayanarak yapılan rökonstrüksiyonu sunulacaktır.

5. MATERYAL VE YÖNTEMYıllık Halka Verisi Örneklerin alınması için yöre seçiminde, kuraklık etkilerini yansıtacak ve yaşlı bireylerden oluşan yörelerin seçimine özen gösterilmiştir. Bu amaçla, güneye bakan eğimli alanlar üzerinde bulunan yaşlı meşcereler tercih edilmiştir.

6. Örnek alınan yöreler (GOV, CAT, KAS, INC) ve Eskişehir Meteoroloji İstasyonu (E.M.İ)

7. Yıllık Halka Ölçümleri ve Kronoloji Oluşturma Kalemlere ait yıllık halka genişlikleri, LINTAB-TSAP ölçüm sistemiyle, 0,001 mm duyarlılıkta ölçülmüştür. Eksik ve yalancı halkalar COFECHA programı kullanılarak giderilmiştir

8. Yıllık halka genişliği serileri, başta yaşa bağlı etkiler olmak üzere iklimden kaynaklanmayan etkileri gidermek amacıyla standart hale dönüştürülmüştür. Standartlaştırılma işlemiyle indis değerlere dönüştürülmüş kronolojilere, otokorelasyonu ortadan kaldırmak amacıyla AR1 otoregressif modeli uygulanmıştır. Bireysel kalemlerden bu şekilde elde edilen indis değerleri, biweight robust ortalama yöntemi kullanılarak her bir yöre için yöre kronolojileri elde edilmiştir

9. Yöre kronolojilerinin iklime verdiği tepkiyi ölçmek amacıyla, her bir kronoloji için “ortalama duyarlılık katsayıları” ve “işaretin gürültüye oranı” hesaplanmıştır. Elde edilen dört yöre kronolojisinin birbiriyle uyumu eğrilerin uyum yüzdesi değerleri hesaplanarak denetlenmiştir.

10. İklim Verisi İklim verisi olarak Eskişehir Meteoroloji İstasyonu aylık toplam yağış ve aylık ortalama sıcaklık verileri kullanılmıştır. İklim verileri, takvim yılı değil, biyolojik yıl temel alınarak düzenlenmiştir. Buna göre, aylık toplam yağış ve aylık ortalama sıcaklık değerleri, önceki yılın ekim ayından halka oluşum yılının eylül ayına kadar olan dönem için 1930-2002 aralığında sıralanmıştır.

11. Kalibrasyon ve Rökonstrüksiyon Yöntemleri Yöre kronolojileri ile aylık toplam yağış ve aylık ortalama sıcaklık değerleri arasındaki ilişki, basit korelasyon katsayısı yöntemi ile belirlenmiş ve t testiyle ilişkinin anlamlılığı denetlenmiştir. Buna dayanarak rekonstrüksiyonda kullanılacak sıcaklık ve yağış değişkenlerine karar verilmiştir. Dört yöre kronolojisi, ortak zaman aralığında sıralanarak bir yıllık halka veri seti elde edilmiştir. Bu veri seti, Ana Bileşen Analizi (PCA) ile birbirine dik bağımsız değişkenlere dönüştürülmüştür.

12. Elde edilen dört ana bileşenden (PC) veri setini %64 oranında ifade eden birinci ana bileşen, rekonstrüksiyonda bağımsız değişken olarak kullanılmıştır. 1930-2002 aralığındaki iklim verisi ikiye bölünerek rekonstrüksiyonda iki ayrı model oluşturulmuştur. Birinci modelde, 1930-1966 dönemi kalibrasyonda, 1967-2002 dönemi verifikasyonda; ikinci modelde ise 1967-2002 dönemi kalibrasyonda, 1930-1966 dönemi verifikasyonda kullanılmıştır.

13. Kalibrasyonda iklim verisi bağımlı değişken, PC1 bağımsız değişen olarak girilerek regresyon katsayıları elde edilmiş, daha sonra bu regresyon katsayıları denklemde yerine konularak iklim değişkenleri tahmin edilmiştir. Verifikasyon döneminde ise oluşturulan modelin anlamlılığı denetlenmiştir. Modelin denetlenmesinde pearson korelasyon katsayısı, ayarlı R² ve hata indirgeme değerleri (RE) hesaplanmıştır.

14. Tahmini yağış değerlerinin 1 standart sapma üstündeki yıllar yağışlı, 2 standart sapma üstündeki yıllar çok yağışlı; 1 standart sapma altındaki yıllar kurak, 2 standart sapma altındaki yıllar ise çok kurak yıllar olarak belirlenmiştir

15. BULGULARYöre Kronolojileri GOV 1710-2002 (293 yıl), INC 1674-2002 (329 yıl), CAT 1734-2002 (264 yıl) ve KAS 1757-2002 (246 yıl) yıllarını kapsamaktadır. Buna göre en uzun kronoloji 329 yıl uzunluğunda INC, en kısa kronoloji ise 246 yıl uzunluğundaki KAS yöre kronolojisidir.

16. Tablo 2: Yöre kronolojilerine ait ARSTAN programından elde edilen istatistik sonuçlar.

18. Yöre Kronolojileri arasındaki eğrilerin uyum yüzdesi (EUY) değerleri ‘***’ işareti 0.999 güven düzeyinde anlamlı olduğunu göstermektedir.

19. Aylık toplam yağış değerleri ile yöre kronolojileri arasındaki ilişki

20. Aylık toplam sıcaklık değerleri ile yöre kronolojileri arasındaki ilişki

21. 1757-2002 yıllarını kapsayan yıllık halka veri seti, Ana Bileşen Analizi ile (PCA) birbirine dik bağımsız değişkenlere dönüştürülmüştür. Birinci ana bileşen (PC1) toplam varyansı %64 oranında, ikici ana bileşen (PC2) %17, üçüncü ana bileşen (PC3) %11 ve dördüncü ana bileşen (PC4) %8 oranında açıklamaktadır. Rekonstrüksiyonda, toplam varyansı %64 oranında açıklayan birinci ana bileşen, bağımsız değişen olarak kullanılmıştır.

22. Tablo 4. Kalibrasyon ve verifikasyon istatistik sonuçları Y (mm) = -136,243 + 116,2161 * PC1

23. Gerçek ve tahmini Mayıs-Ağustos toplam yağışı

24. Tahmini Mayıs-Ağustos toplam yağışı. Yatay düz çizgi ortalamayı, yatay kesikli çizgiler ±1 standart sapma (içtekiler) ve ±2 standart sapmayı (dıştakiler) göstermektedir.

25. TARTIŞMA VE SONUÇ H) Hughes ve diğ. (2001); D) D’Arrigo ve Cullen (2001); T) Touchan ve diğ. (2003); A) Akkemik ve Aras (2005); Ü) Akkemik ve diğ. (2005); 1) Osmanlı Arşivleri; 2) Inalcik (1997); 3) Erler (2000)

26. Projeden elde edilen bu sonuçlar başlangıç niteliğinde olup, çalışmalar halen devam etmektedir. Ülkemiz ormanlarındaki yaşlı meşcerelerden alınan/alınacak örneklerle rekonstrüksiyonların sayısı arttırılacaktır. Ülkemiz iklim tarihinin ortaya konmasında büyük ölçekli iklim olaylarının da etkisi bulunduğundan, proje kapsamında uzay-zaman analizleri de yapılarak, iklim etkileri ile yağış arasındaki ilişkilerin bölgeler arasındaki benzerlik ve farklılık gösterdiği alanlar tanımlanacaktır.