📘🌍 EQ-LiveTurkey Teknik Raporu

🌍 SeismoReport v3.1 · Otomatik Çeviri

EQ-LiveTürkiye · Sismik İzleme Platformu · Teknik Metodoloji Raporu · SeismoReport 2026

📡 ISC 1900–1997 + EMSC 1998–2026 📐 Scordilis (2006) Mw Homojenizasyonu 📊 Gutenberg-Richter b-değeri Izgarası ⚡ EMSC Gerçek Zamanlı Veri 🛠️ OpenQuake · Leaflet · Python 🇹🇷 83 918 Olay · 19 719 Grid Noktası

🌐 EQ-LiveTürkiye

Sismik İzleme ve Tehlike Görselleştirme Platformu

Prof. Dr. Ali Osman Öncel

İstanbul Üniversitesi–Cerrahpaşa · Jeofizik Mühendisliği Bölümü

🗺️ Canlı Harita ⭐ EQ-LiveTurkiye 🔶 GEM / oq-engine 📡 EMSC FDSN API 🌐 ISC FDSN API

83 918

Katalog Olayı

126 yıl

Zaman Aralığı

5

Tamamlılık Dönemi

1.031

Türkiye b-değeri

19 719

Grid Noktası

100 km

Sabit Alan Yarıçapı

M ≥ 3.0

EMSC Tamamlılık Eşiği

Gerçek Zamanlı

EMSC Canlı Veri

ÖZET EQ-LiveTürkiye, Türkiye için gerçek zamanlı deprem akışını bilimsel açıdan sağlam bir Gutenberg-Richter b-değeri tehlike katmanıyla birleştiren açık kaynaklı bir sismik izleme platformudur. Platform; ISC (1900–1997) ve EMSC (1998–2026) verilerinden dönem bazlı tamamlılık eşikleriyle (merdiven yöntemi) derlenen 126 yıllık bileşik bir katalog kullanmaktadır. Büyüklük homojenizasyonu Scordilis (2006) ampirik ilişkileriyle (ML/MD/mb → Mw) gerçekleştirilmiştir. b-değeri ızgarası, Öncel & Wyss (2000) Sabit Alan Yöntemi (R=100 km, 10 km kaydırma adımı, N≥50) çerçevesinde Aki (1965) MLE ve Utsu (1966) binom düzeltmesiyle hesaplanmıştır. Tüm Türkiye için elde edilen b=1.031 değeri literatürle uyumludur.

Deprem kataloğuGutenberg-Richterb-değeri Türkiye sismisitesiMerdiven tamamlılıkScordilis (2006) Aki (1965) MLEÖncel & Wyss (2000)PSHAOpenQuake

🌐 1. Platforma Genel Bakış — EQ-LiveTürkiye Ne Gösteriyor?

Amaç

EQ-LiveTürkiye, Türkiye için tarayıcı tabanlı bir sismik izleme ve tehlike görselleştirme platformudur. Birbirini tamamlayan iki veri akışını bir arada sunar: (1) son depremleri otomatik olarak güncelleyen gerçek zamanlı EMSC deprem akışı ve (2) 126 yıllık bileşik bir katalogdan önceden hesaplanmış Gutenberg-Richter b-değeri ızgarası. Her iki katman da etkileşimli Leaflet.js haritası üzerinde çakıştırılmakta; böylece süregelen sismisitenin ve uzun vadeli sismik tehlike örüntülerinin eş zamanlı değerlendirilmesi mümkün kılınmaktadır.

Harita Katmanları

🔴 Son Depremler
EMSC canlı

+

🟡 b-değeri Izgarası
CircleMarker

+

🗺️ Altlık Harita
CartoDB Dark

=

🌐 Etkileşimli Harita
Leaflet.js

Açık Kaynak Araçlar ve Veri Kaynakları

Tablo 1. EQ-LiveTürkiye'de kullanılan araçlar ve veri kaynakları

Bileşen	Araç / Kaynak	Amaç
Canlı harita	Leaflet.js v1.9	Etkileşimli harita, CircleMarker, LayerGroup
Deprem akışı	EMSC FDSN API	Gerçek zamanlı M≥2.0 olaylar, Türkiye sınır kutusu
Tarihsel katalog (1900–1997)	ISC FDSN API (web servis)	Merdiven dönemleri M≥4.5–5.5
Tarihsel katalog (1998–2026)	EMSC FDSN API	Homojen dönem M≥3.0
Katalog derleyici	fetch_combined_catalog.py	ISC+EMSC birleştirme, tekrar temizleme, Scordilis Mw
b-değeri ızgarası	build_bvalue_grid.py	Aki MLE, Sabit Alan Yöntemi, JSON çıktı
PSHA çerçevesi (referans)	GEM / oq-engine	b-değeri OpenQuake kaynak modeline girdi olarak kullanılır
Yayın	GitHub Pages	index.html + veri JSON statik yayın

🧭 2. Haritayı Nasıl Kullanırım? — Hızlı Kılavuz

Harita açılınca ilk bakışta sade görünür ama aslında altında bir kontrol paneli saklıdır. Aşağıdakileri bilirseniz haritayı çok daha etkili kullanırsınız:

💡 Bilmeniz gereken 6 şey

• Kontrol paneli haritanın ALTINDA gizli. Ekranın en altındaki ▲ ▲ ▲ ok işaretine tıklayın — tüm butonlar ve ayarlar açılır.
• Aktif olan buton içi mavi (camgöbeği) dolu ve yazısı koyu olandır. Pasif butonlar şeffaftır.
• Açılışta "Isı Haritası" görünür — bu, depremlerin yoğunluğunu soluk bir sis gibi gösterir. Depremleri tek tek, büyüklüğüne göre renkli noktalar hâlinde görmek isterseniz "Seismicity" butonuna basın.
• TR / MAR / EGE / DAF / KAF butonları yalnızca haritayı o bölgeye yakınlaştırır. Analizi o bölgeyle sınırlamaz — alttaki sayılar (Toplam, Mw≥3, Mw≥4, Maks…) her zaman o an haritada yüklü olan tüm veriyi gösterir, sadece baktığınız bölgeyi değil.
• Depremler ortalama 7-10 dakikada bir geliyor, bazen bir saate kadar sessiz kalabilir — "hareket yok" görmeniz çoğu zaman normaldir. Veri saatte bir otomatik yenilenir; sağ alttaki "EMSC Canlı ●" ibaresi bağlantının çalıştığını gösterir.
• b-değeri, T_L Asperite gibi bilimsel katmanlar önceden hesaplanmış Türkiye-geneli haritalardır — bir bölgeye yakınlaşsanız bile hep ülke genelini gösterirler, o bölgeye özel yeniden hesaplanmazlar.

Canlı haritayı doğrudan açmak için: aliosmanoncel.github.io/EQ-LiveTurkiye

🔄 3. Veri Boru Hattı — Ham Katalogdan Izgara Üretimine

📡 ISC FDSN
1900–1997

+

📡 EMSC FDSN
1998–2026

→

🔀 Merdiven Birleştirme
dönem M_c

→

🔁 Tekrar Temizleme
±60s, <50km

→

⚖️ Mw Homojenizasyon
Scordilis 2006

→

📊 b-değeri Izgarası
Aki+Utsu

→

🌐 JSON → Harita
GitHub Pages

Adım 1 — FDSN Veri Çekimi

Hem ISC hem EMSC, FDSN uyumlu web servis API'leri (fdsnws/event/1/query) sunmaktadır. Python betiği fetch_combined_catalog.py, API sonuç sınırlarını aşmamak için her servisi 5–10 yıllık dilimler hâlinde sorgular (sorgu başına 20 000 olay sınırı). Coğrafi sınır kutusu 33–45°K, 23–48°D olarak belirlenmiştir.

⚙️ Çekim Parametreleri

• Biçim: text (boru karakteriyle ayrılmış, hafif)
• Sıralama: dilim başına time-asc
• Hız sınırı: istekler arası 1.5 sn bekleme
• Yeniden deneme: ağ hatasında 2× gecikme
• Dilim boyutu: 10 yıl (ISC) / 5 yıl (EMSC)

Adım 2 — Merdiven (Staircase) Tamamlılık Birleştirmesi

Sismik ağ yoğunluğu tarihsel süreçte evrilmiştir; 126 yıllık kaydın tamamına tek bir büyüklük eşiği uygulamak, ya büyük sayıda küçük olayı dışarıda bırakır ya da tamamlanmamış erken dönemleri dahil eder. Merdiven yaklaşımı, her zaman dilimine uygun minimum tamamlılık büyüklüğünü (M_c) uygulayarak ve ardından tüm kümelerin bileşimini (küme birleşimi) alarak bu sorunu çözer.

Tablo 2. Merdiven tamamlılık eşikleri (+0.5 muhafazakâr düzeltme uygulandı)

Dönem	Kaynak	Mc	Olay Sayısı	Gerekçe
1900–1965	ISC	5.5	173	Pre-WWSSN; yalnız büyük olaylar güvenilir
1965–1980	ISC	4.5	1 058	WWSSN tam operasyonel; ref 4.0+0.5
1980–1990	ISC	4.0	1 318	Global ağ yoğunlaşması; ref 3.0+1.0
1990–1998	ISC	3.5	4 855	Dijital geniş bant geçişi; ref 2.5+1.0
1998–2026	EMSC	3.0	76 506	Modern geniş bant; ref 2.5+0.5
Toplam (tekrar temizleme sonrası)			83 918

Adım 3 — Tekrar Temizleme

Hem ISC hem EMSC kataloglarında yer alan olaylar zamansal ve konumsal yakınlıkla belirlenir. Birbirinden 60 saniye ve 50 km'den az uzakta bulunan çiftler birleştirilerek daha yüksek Mw değerine sahip olan tutulur. Bu işlem 396 yinelenen çifti kaldırarak 83 918 benzersiz olay elde edilmesini sağlar.

⚖️ 4. Büyüklük Homojenizasyonu — Scordilis (2006)

Homojenizasyon Neden Gereklidir?

Bileşik ISC+EMSC kataloğu, en az sekiz farklı büyüklük ölçeğinde (MD, ML, mb, Ms, Mw ve türevleri) rapor edilmiş olaylar içermektedir. Ölçeklerin karıştırılması sistematik önyargıya yol açar: mb=5.0, Mw=5.0 ile eşdeğer değildir. PSHA için tüm büyüklüklerin tek ve fiziksel olarak anlamlı bir ölçekte — moment büyüklüğü Mw — ifade edilmesi gerekmektedir.

Tablo 3. Bileşik katalogdaki büyüklük türü dağılımı

Ölçek	Olay Sayısı	Birincil Kaynak
MD (koda süresi)	32 686	EMSC bölgesel ağlar
ML (yerel)	25 991	EMSC bölgesel ağlar
mb (cisim dalgası)	9 847	ISC küresel istasyonlar
ML (büyük harf)	9 232	EMSC
md	2 545	EMSC
mw / Mw	1 294	ISC/EMSC moment tensör çözümleri
mL	707	EMSC

Büyüklük Dönüşüm Bağıntıları

mb → Mw (cisim dalgası büyüklüğü)

Mw = 0.85·mb + 1.03   (3.5 ≤ mb ≤ 6.2)

Scordilis (2006). ISC küresel istasyon ağından elde edilen teleseismik cisim dalgası büyüklükleri için geçerlidir.

ML / MD → Mw (yerel / koda büyüklüğü)

Mw = 0.0376·M² + 0.646·M − 0.269   (1.0 ≤ M ≤ 6.5)

İkinci dereceden bağıntı, yerel büyüklüğün yüksek değerlerdeki doğrusal olmayan doygunluğunu hesaba katar; Scordilis (2006) ve büyüklük türlerine uygun diğer ampirik dönüşümlerden derlenmiştir. ML, MD, md ve mL'ye uygulanır.

Ms → Mw (yüzey dalgası büyüklüğü)

Mw = 0.646·Ms + 2.079   (3.0 ≤ Ms ≤ 6.1)
Mw = 0.994·Ms + 0.115   (Ms > 6.1)

İki parçalı bağıntı, büyük depremlerde Ms'nin doygunluğunu yansıtır.

⚠️ b-Değeri Hesabında Dönüştürülmüş-Mw Önyargısı

• ML→Mw dönüşümü: ML=3.0 için → Mw = 0.0376×9 + 0.646×3 − 0.269 = 2.007
• Dönüştürülmüş Mw≥3.0 filtresi ~66 000 EMSC olayını dışarıda bırakır (MD/ML ≈ 3.0–4.5)
• Bu durum ortalama(Mw)'yi ~2.5'ten ~4.2'ye yükseltir → yapay düşük b değerleri (~0.28) üretir
• Çözüm: b-değeri ızgarası dönüştürülmüş Mw değil, orijinal mag kullanılarak hesaplanır

Dönüştürülmüş Mw alanı eq_historical.json dosyasında diğer kullanımlar için (PSHA kaynak modeli, moment hız hesapları) korunmaktadır.

📊 5. b-Değeri Metodolojisi — Aki (1965) + Öncel & Wyss (2000)

Gutenberg-Richter Frekans-Büyüklük Bağıntısı

log₁₀(N) = a − b·M

b-değeri eğimi, büyük ve küçük depremler arasındaki göreli oranı tanımlar. Küresel ortalama b≈1, büyüklükteki her bir birim artışın frekansta 10 katlık bir düşüşe karşılık geldiği anlamına gelir. b değerinin yanal değişimleri stres durumundaki, malzeme heterojenliğindeki ve sıcaklıktaki farklılıkları yansıtır.

Vatandaş için özet: b=1 ise M5'ten büyük her 10 depremin karşısında M6'dan büyük yaklaşık 1 deprem beklenir. b küçüldükçe büyük depremler göreceli olarak daha sık; b büyüdükçe küçük depremler daha baskın hâle gelir.

Aki (1965) Maksimum Olabilirlik Tahmini (MLE)

b = log₁₀(e) / (M̄ − M_c,eff)

MLE, en küçük kareler regresyonuna kıyasla istatistiksel olarak üstündür; bağımlı değişkenin dağılımı hakkında varsayımda bulunmaz ve yansız tahminler üretir. M̄, M_c üzerindeki örneklem ortalama büyüklüğüdür.

Utsu (1966) Binom Düzeltmesi

M_c,eff = M_c − ΔM = 3.0 − 0.05 = 2.95

Deprem katalogları büyüklükleri genellikle en yakın 0.1 birimine yuvarlayarak raporlar. Utsu (1966), bu ayrıklaştırmanın b'yi sistematik olarak küçümsediğini göstermiştir; düzeltme ΔM = aralık genişliği/2 = 0.05 olarak uygulanır.

Öncel & Wyss (2000) Sabit Alan Yöntemi

Kaynak: Öncel, A. O. & Wyss, M. (2000). The major asperities of the 1999 Mw 7.4 İzmit earthquake defined by the microseismicity of the two decades before it. Geophysical Journal International, 143(2), 501–506.

Sabit Alan Yöntemi'nde, her grid düğümü etrafında sabit yarıçaplı bir çember içindeki olaylar örneklenerek b-değeri hesaplanır. Bu yaklaşım, hareketli pencere yöntemine kıyasla alansal örneklem büyüklüğünü tutarlı kılar ve alan kaynaklı önyargıyı azaltır. Öncel & Wyss (2000), İzmit kırık zonu için R=20 km ve 5 km adım kullanmıştır; bu çalışmada Türkiye geneli analiz için R=100 km ve 10 km adım kullanılmıştır.

Tablo 4. Sabit Alan Yöntemi parametreleri

Parametre	Değer	Gerekçe
Arama yarıçapı R	100 km	Bölgesel ölçek; Türkiye genelinde N≥50'yi sağlar
Grid adımı	0.09° (~10 km)	Öncel & Wyss: 5 km; Türkiye geneline ölçeklendirildi
Min. olay sayısı	Nmin = 50	Öncel & Wyss (2000) kriteri
Tamamlılık eşiği	Mc = 3.0	EMSC 1998–2026 homojen dönemi
Efektif eşik	Mc,eff = 2.95	Utsu düzeltmesi ΔM=0.05
Grid alanı	34–43°K, 25–45°D	Türkiye kara alanı + yakın kıyı
Grid boyutu	101×223 = 22 523 aday	0.09° adım
Haritalanan nokta	19 719 (%87.5)	2 804 nokta N≥50 kriterini karşılayamadı

🆕 Zaman İçinde b-Değeri — Kayan Pencere Analizi (Time–Magnitude Panel)

Yukarıdaki b-değeri ızgarası tek bir sabit değer verir — 126 yıllık kataloğun tamamına ait, zamanla değişmeyen bir özet. Canlı haritaya eklenen Zaman–Magnitude (Time–Magnitude) paneli ise farklı bir soruya cevap arar: seçili bölge ve pencerede b-değeri zaman içinde nasıl değişiyor?

Kaynak: Öncel, A. O. & Wilson, T. (2007). Anomalous seismicity preceding the 1999 Izmit event, NW Turkey. Geophysical Journal International, 169(1), 259–270.

Yöntem, bu makalenin §3.1'inde uygulanan tekniğin birebir aynısıdır: katalog zaman sırasına dizilir, sabit N=100 olaylık bir pencere tanımlanır, pencere içindeki olaylara Aki (1965) MLE uygulanarak bir b-değeri hesaplanır; ardından pencere 10 olay ileri kaydırılıp işlem tekrarlanır. Sonuç, zaman ekseninde ilerleyen bir b-değeri dizisidir.

⚙️ Neden N=100, neden adım=10?

• İstatistiksel güç: Aki MLE'nin %95 güven aralığı ±1,966/√n ile daralır; n=100 için ±0,196 — kabul edilebilir bir kesinlik. N çok küçük olsaydı (ör. 20-30), tek bir aykırı deprem b-değerini büyük ölçüde oynatabilirdi.
• Zamansal çözünürlük: N büyüdükçe her pencere daha uzun bir zaman aralığını temsil eder, kısa vadeli değişimler bulanıklaşır. N=100, Öncel & Wilson (2007)'da İzmit öncesi 8,5 yıllık dönemi ~70 anlamlı pencereye ayırmaya yetmiştir.
• Adım=10 (pencere N'sinin %10'u): ardışık pencereler arasında güçlü örtüşme sağlar, ani sıçramalar yerine yumuşak bir zaman serisi üretir; komşu pencereler arası korelasyon yüksek kalır.

Global b-değeri ile Kayan (Rolling) b-değeri farkı
Global b: Tüm uygun katalog (N=76 514, M≥Mc) üzerinden hesaplanan tek sayı — "Türkiye'nin ortalama sismisite karakteri nedir?" sorusuna cevap verir.
Rolling b: N=100'lük kayan pencerelerden üretilen bir dizi — "sismisite karakteri şu anda nasıl değişiyor?" sorusuna cevap verir. İkisi farklı sorulara cevap verir, biri diğerinin yerine geçmez.

Mc (tamamlılık eşiği), kayan pencerede pencere başına değil, görüş alanı için global/sabit olarak hesaplanır (MAXC yöntemi — Wiemer & Wyss 2000). Pencere başına ayrı Mc hesaplamak, küçük örneklemlerde gürültülü ve tutarsız sonuçlar üretirdi.

⚠️ Bu bir deprem tahmin aracı değildir

Scholz (2019, The Mechanics of Earthquakes and Faulting, §4.4.2)'a göre swarm dizileri genellikle yüksek b-değeriyle ilişkilendirilir (pore-fluid diffüzyonu, kendini sınırlayan bir mekanizma), düşük b ise kilitlenme/mainshock-aftershock tipik davranışına işaret eder. Ancak aynı kaynak, Mogi (1963)'nin dizi tiplerini kesin sınıflandırma iddiasını "fazla genellenmiş" bulur — aynı bölgede, aynı fay mekanizmasıyla, birkaç km fark ile tamamen farklı dizi tipleri gözlenebilir. Bu nedenle platform hiçbir zaman "bu swarm'dır" ya da "foreshock tespit edildi" gibi kesin bir sınıflandırma üretmez; yalnızca "swarm-benzeri" veya "artçı-şok-benzeri" gibi geriye-dönük-doğrulanabilir ifadeler kullanır. Foreshock, tanımı gereği ancak daha büyük bir olay gerçekleştikten sonra geriye dönük olarak teşhis edilebilir.

Panelin "Forecast" ve "Marmara FC" katmanları da aynı temkinle okunmalıdır: bunlar deprem zamanı veya yerinin öngörüsü değil, geçmiş dönem verisiyle eğitilip (1998–2019) sonraki dönemde (2020–2030) test edilen istatistiksel bir doğrulama denemesidir (Wiemer & Wyss 2002 yaklaşımı) — b-değeri haritasının öngörü gücünü retrospektif olarak sınamak içindir, gelecekteki bir depremin kesin konumunu vermez.

🗺️ 6. Etkileşimli b-Değeri Haritası

Aşağıdaki harita 19 719 grid noktasından oluşan b-değeri dağılımını göstermektedir. Her daire, merkezdeki 100 km yarıçaplı alan içindeki depremlere dayanarak Aki (1965) MLE yöntemiyle hesaplanan b-değerini temsil eder. Fareyi bir noktanın üzerine getirerek o bölgenin b-değerini ve kullanılan deprem sayısını (N) görebilirsiniz. Yakınlaştırdıkça ızgara daha yüksek çözünürlüğe geçer.

■ b≈0.5 (yüksek stres)

■ b≈1.5 (düşük stres) Aki (1965) MLE · Utsu (1966) · M_c=2.95 · R=100 km · N≥50 · EMSC 1998–2026

⚠️ Bu harita bir deprem tahmin haritası değildir

b-değeri haritası tek başına deprem tehlikesini veya gelecekteki büyük bir depremin yerini göstermez; sismisite içindeki büyüklük-frekans dağılımının mekânsal değişimini gösterir. Düşük b bölgeleri stres birikimiyle ilişkilendirilen istatistiksel bir örüntüdür, kesin bir öngörü aracı değildir.

Bölgesel Yorumlama

🔴 Düşük b-değeri bölgeleri (b < 0.8) — olası yüksek gerilim/stres ile ilişkili anomaliler

Doğu Anadolu Fay Zonu (DAFZ) boyunca, özellikle Kahramanmaraş üçlü kavşağı çevresinde belirginleşmektedir. Ancak bu anomalilerin Şubat 2023 Mw7.8–7.7 deprem dizisiyle nedensel olarak ilişkilendirilmesi için artçı kümelerinin giderilmesi (declustering) ve yerel M_c/belirsizlik analizlerinin yapılması gerekmektedir — bkz. Kısıtlar bölümü.

🟢 Yüksek b-değeri bölgeleri (b > 1.3) — Heterojen / Jeotermal

Batı Anadolu genişleme kuşağı (Gediz, Büyük Menderes grabenleri), Orta Anadolu volkanik yayı ve Karadeniz kıyı kuşağının bazı kesimleri. Bu alanlardaki yüksek b-değerleri kabuk heterojenliğini, yüksek gözenek sıvısı basınçlarını veya ısıl olarak zayıflamış fay zonlarını yansıtır.

✅ Doğrulama

• Tüm Türkiye b = 1.031 (N=76 514) — literatürle uyumlu (Kalafat vd. 2007: b≈1.0–1.1)
• b aralığı 0.460–2.249, ortalama 1.167 — tektonik açıdan çeşitli bölge için fiziksel olarak makul
• 2 804 haritalanmamış nokta (N<50) Karadeniz açığı ve düşük sismisiteli Orta Anadolu'da yoğunlaşıyor

Neden 83 918 değil 76 514? b-değeri ızgarası yalnızca homojen EMSC dönemini kullanır (1998–2026, M_c=3.0 → Tablo 2'de 76 506 olay). 83 918, ISC (1900–1997) dönemini de içeren 126 yıllık bileşik kataloğun tamamıdır; farklı dönemler farklı tamamlılık eşikleri taşıdığından (merdiven yöntemi), tutarlı tek bir M_c gerektiren b-değeri hesabında yalnızca en güncel homojen dönem kullanılmıştır.

🔶 7. PSHA Entegrasyonu — GEM OpenQuake Engine

Olasılıksal Sismik Tehlike Analizi (PSHA), belirli bir zaman penceresinde bir alanda verilen zemin hareketi yoğunluk ölçütünü (örn. PGA, Sa) aşma olasılığını nicel olarak ortaya koyar. GEM Vakfı'nın OpenQuake Engine'i, Küresel Deprem Modeli'nde ve ulusal tehlike çalışmalarında kullanılan önde gelen açık kaynaklı PSHA platformudur. OpenQuake kaynak modelinde her alan kaynağı veya nokta kaynağı, bir Gutenberg-Richter MFD (Büyüklük-Frekans Dağılımı) gerektirmektedir.

OpenQuake source_model.xml — b-değeri girdisi örneği:
<areaSource id="TR_01" name="Marmara Bölgesi">
  <truncGutenbergRichterMFD aValue="5.2" bValue="0.95"
    minMag="5.0" maxMag="7.5"/>
</areaSource>
bValue alanı, ilgili kaynak zonuyla çakışan b-değeri grid hücrelerinin uzamsal ortalamasından alınır.

bvalue_grid.json
19 719 nokta

→

Zon ortalaması
konumsal birleştirme

→

source_model.xml
zon başına b

→

oq engine
PSHA

→

Tehlike eğrileri
PGA, Sa haritaları

⚠️ 8. Kısıtlar ve Gelecek Çalışmalar

⚠️ Mevcut Kısıtlar

• Kümeleme giderme yok: Artçı şok sürgünleri (örn. 2023 Kahramanmaraş: ~10 000 olay) etkilenen bölgelerde b-değerlerini şişirir. Gardner-Knopoff (1974) veya Reasenberg (1985) henüz uygulanmadı.
• Sabit M_c=3.0: Yerel M_c, Türkiye genelinde değişkenlik gösterebilir. Wiemer & Wyss (2000) MAXC yöntemi konumsal değişken M_c sağlayabilir.
• Belirsizlik sınırı yok: Shi & Bolt (1982) b-değeri güven aralıkları henüz haritada gösterilmemektedir.
• Statik ızgara: b-değeri ızgarası önceden hesaplanmıştır; yeni EMSC verileri geldiğinde otomatik güncelleme yapılmamaktadır.

✅ Planlanan İyileştirmeler

• Gardner-Knopoff kümeleme giderme (yalnız ana şok kataloğu)
• MAXC konumsal değişken M_c haritalama
• Shi & Bolt (1982) b-değeri belirsizlik gösterimi
• GitHub Actions ile otomatik katalog güncelleme (aylık yenileme)
• OpenQuake kaynak modeli XML dışa aktarım betiği

📖 9. Kaynaklar

Aki, K. (1965). Maximum likelihood estimate of b in the formula log N = a − bM and its confidence limits. Bulletin of the Earthquake Research Institute, 43, 237–239.

Gardner, J. K. & Knopoff, L. (1974). Is the sequence of earthquakes in Southern California, with aftershocks removed, Poissonian? Bulletin of the Seismological Society of America, 64(5), 1363–1367.

GEM Foundation (2023). OpenQuake Engine. github.com/gem/oq-engine

Gutenberg, B. & Richter, C. F. (1944). Frequency of earthquakes in California. Bulletin of the Seismological Society of America, 34(4), 185–188.

ISC — International Seismological Centre. FDSN Web Service. isc.ac.uk/fdsnws/event/1/

EMSC — European Mediterranean Seismological Centre. FDSN Web Service. seismicportal.eu/fdsnws/event/1/

Mogi, K. (1963). Some discussions on aftershocks, foreshocks and earthquake swarms. Bulletin of the Earthquake Research Institute, 41, 615–658.

Öncel, A. O. & Wilson, T. (2007). Anomalous seismicity preceding the 1999 Izmit event, NW Turkey. Geophysical Journal International, 169(1), 259–270.

Öncel, A. O. & Wyss, M. (2000). The major asperities of the 1999 Mw 7.4 İzmit earthquake defined by the microseismicity of the two decades before it. Geophysical Journal International, 143(2), 501–506.

Reasenberg, P. (1985). Second-order moment of central California seismicity, 1969–1982. Journal of Geophysical Research, 90(B7), 5479–5495.

Scholz, C. H. (2019). The Mechanics of Earthquakes and Faulting (3rd ed.). Cambridge University Press.

Scordilis, E. M. (2006). Empirical global relations converting Ms and mb to moment magnitude. Journal of Seismology, 10(2), 225–236.

Shi, Y. & Bolt, B. A. (1982). The standard error of the magnitude-frequency b value. Bulletin of the Seismological Society of America, 72(5), 1677–1687.

Utsu, T. (1966). A statistical significance test of the difference in b-value between two earthquake groups. Journal of Physics of the Earth, 14(2), 37–40.

Wiemer, S. & Wyss, M. (2000). Minimum magnitude of completeness in earthquake catalogs. Bulletin of the Seismological Society of America, 90(4), 859–869.

Woessner, J. & Wiemer, S. (2005). Assessing the quality of earthquake catalogues. Bulletin of the Seismological Society of America, 95(2), 684–698.

SeismoReport 2026 · Prof. Dr. Ali Osman Öncel · İstanbul Üniversitesi–Cerrahpaşa · Jeofizik Mühendisliği
aliosmanoncel.github.io/EQ-LiveTurkiye · github.com/aliosmanoncel/EQ-LiveTurkiye