AFAD Event Web Service Earthquake Data Download Package: AFADEarthQuakeData¹²

Bu paketin amacı; Acil ve Afet Durum Yönetim Başkanlığının (AFAD) sunmakta olduğu Deprem Web Servisi aracılığı ile dağıtımı yapılan deprem katalog verilerinin istatistik hesaplama ve grafik yazılımı R'ye aktarınını kolaylaştırmak ve açıklayıcı veri analizine yönelik bazı temel istatistiklerin grafiklerin hızlıca elde edilebilmesini sağlamaktır. Paket içeriği an itibari 10 ögeden oluşmaktadır:

Fonksiyon	Açıklama
calc.Havershine	Verilen bir nokta koordinatlarına göre deprem olaylarının Havershine uzaklıklarını hesaplar.
calc.meanCenter	Ortalama ve ağırlıklı ortalama merkez değerlerini hesaplar.
calc.StDevEllipse	Standart elipsi hesaplar.
calc.StDistance	Standard uzaklığı hesaplar.
histog.depth	İndirilen deprem verisi için hızlı deprem derinliği histogramı çizer.
histog.magni	İndirilen deprem verisi için hızlı deprem büyüklüğü histogramı çizer.
map.AFAD	maps paketi ile indirilen deprem verisinin haritalar.
opensmap.AFAD	OpenStreetMap paketi ile indirilen deprem verisini haritalar.
read.AFAD	AFAD Deprem Web Servisinden deprem verisi çeker.
turkiye_and_neighborhood_gazetteer	Türkiye ve Çevresi Coğrafi Ad Dizini (CAD)

Örneğin, eventID'si bilinen bir depreme ilişkin özellikler, Web Servisinden çekilebilir.

  equake<-read.AFAD(eventid=512359)
# View(equake)
# equake
#  country                   date depth district eventid latitude        location longitude magnitude neighborhood  province  rms type iseventupdate lastupdatedate
#1 Türkiye 2021-09-16T09:40:27.44  7.07      Çan  512359  40.0066 Çan (Çanakkale)   27.0065       1.6        Hurma Çanakkale 0.71   ML         false               

Diğer bir örnek de belli bir zaman aralığını kapsayacak şeklinde deprem verisi çekmek üzerine verilebilir. 01 Eylül - 30 Eylül 2023 tarihleri arasında Türkiye genelinde meydana gelmiş deprem verisi çekilmek istendiğinde aşağıdaki şekilde bir veri çerçevesi elde edilecektir.

  equake<-read.AFAD(start="2023-09-01T00:00:00",end="2023-09-30T23:59:59")
  head(equake)
#  country                date depth   district eventid latitude             location longitude magnitude neighborhood province  rms type iseventupdate             lastupdatedate
#1    Irak 2023-09-10T21:43:48  7.00             602422   35.024 Daquq, Kerkük (Irak)    44.411       4.8                       0.68   MW          true 2023-09-10T22:12:31.523207
#2         2023-09-16T18:56:31 33.96             603132   34.436              Akdeniz    32.675       1.9                       0.44   ML         false                           
#3 Türkiye 2023-09-01T00:02:01 12.15   Sungurlu  601270   40.120     Sungurlu (Çorum)    34.467       2.0  Büyükincesu    Çorum 0.38   ML         false                           
#4 Türkiye 2023-09-01T00:04:50  7.00    Pütürge  601271   38.135    Pütürge (Malatya)    38.637       1.5       Üçyaka  Malatya 0.39   ML         false                           
#5 Türkiye 2023-09-01T00:06:26  7.00 Bahçesaray  601272   38.109     Bahçesaray (Van)    42.687       2.6      Çiçekli      Van 0.46   ML         false                           
#6 Türkiye 2023-09-01T00:54:51  3.34 Bahçesaray  601273   38.104     Bahçesaray (Van)    42.702       2.1     Ulubeyli      Van 0.63   ML         false    

AFAD Deprem Web Servisinden indirilen veri hızlı bir şekilde haritalandırılmak ve deprem derinlik ve büyüklük dağılımları hızlı bir şekilde görselleştirilmek istenir ise sırası ile map.AFAD ve/veya opensmap.AFAD, histog.depth ve histog.magni fonksiyonlarından yararlanılabilir.

opensmap.AFAD(equake,gap=.5,title="Standart Arama Deprem Haritası",subtitle="Sadece Tarih-Zaman Filtresi: Eylül 2023" ,xlab="Boylam",ylab="Enlem",minnumtiles=10)
histog.depth(equake,title="Eylül 2023 - Deprem Derinlik Dağılımı")
histog.magni(equake,title="Eylül 2023 - Deprem Büyüklük Dağılımı")

Dikdörtgen Arama Örneği

Bu örnekte, map.AFAD fonksiyonu ile kullanılmakta olan RDS harita dosyası için Old GADM data linki üzerinden versiyon 3.6'yı seçerek ulaşılabilir.

equake<-read.AFAD(start="1990-01-01T12:00:00",end="2023-01-02T12:00:00",minlat=39,maxlat=42,minlon=26,maxlon=42,orderby ="magnitude",minmag=4)
#RDS dosyası için https://gadm.org/old_versions.html linki üzerinden versiyon 3.6'yı seçerek ulaşılabilir.
map.AFAD(equake,"gadm36_TUR_0_sp.rds")
opensmap.AFAD(equake,gap=.5,title="Box Search Plot",xlab="Boylam",ylab="Enlem",minnumtiles = 10); 
histog.depth(equake,col="grey",fill = "black")
histog.magni(equake,col="blue",fill = "black")
c.mean<-calc.meanCenter(equake,type="mean")
stdis<-calc.StDistance(equake)
stdev<-calc.StDevEllipse(equake,addplt=TRUE,plt=p)
# Takip eden grafik çizimi için ggforce paketi gerekmektedir.
library(ggforce)
p<-opensmap.AFAD(equake,gap=.5,minnumtiles=10)
p+geom_point(color="red",aes(x=c.mean[,1],y=c.mean[,2]))
+geom_circle(aes(x0=c.mean[,1], y0=c.mean[,2], r=calc.StDistance(equake)), inherit.aes=FALSE)
+geom_ellipse(color="red",aes(x0=c.mean[,1],y0=c.mean[,2],angle=stdev$rot,a=stdev$A,b=stdev$B))
#Deprem büyüklüğü 5'in üzerinde olan depremlerin şiddetlerinin harita üzerinde metin şeklinde gösterilmesi
eqmag<-equake[which(equake$magnitude>5),]
opensmap.AFAD(equake,gap=.5,title="Box Search Plot",xlab="Boylam",ylab="Enlem",minnumtiles = 10)+ geom_text(data=eqmag,aes(longitude,latitude,label=magnitude),color="red")

Github'dan Paket Yükleme

devtools paketi, GitHub'dan paketlerin yüklenmesini sağlayan install_github() işlevini sağlar.

library(devtools)
install_github("lterlemez/AFAD_Package")

Footnotes

1. İçöz, C., & Terlemez, L. (2024). AFAD Deprem Veri ve Görselleştirmelerine İlişkin Bir R Paketi: AFADEarthQuakeData. Türk Deprem Araştırma Dergisi. https://doi.org/10.46464/tdad.1375464 ↩

2. AFAD tarafından yapılan son değişiklik nedeni ile maalesef eski deprem verilerine artık ulaşılamamaktadır. Şu an itibari ile Deprem Web Servis sadece son 5 gün'e ayıp deprem olaylarına ilişkin veriyi sağlamaktadır. Örnek hata sonucunu görmek için https://deprem.afad.gov.tr/apiv2/event/filter?start=2020-09-14%2010:00:00&end=2021-09-16%2010:00:00&format=xml linkinden yararlanabilirsiniz. ↩