Bir Saniyelik Düşey Kanalımızda Gerçekte Ne Var

Her GeoShake istasyonu kesintisiz tek bir kanal yayınlar: LNZ, saniyede bir örnek. 10 Eylül 2026'ya kadar bu kanaldaki örnekler, kanal kodunun söylediği şey değildi. Bu yazı örneklerin eskiden ne olduğunu, şimdi ne olduğunu, değerin nasıl hesaplandığını ve tanımın hâlâ nerelerde yumuşak kaldığını anlatıyor.

Adın vaat ettiği ve eskiden verdiği

SEED kanal adlandırmasında L ~1 örnek/s bandı, N ivmeölçer, Z düşey eksendir. Harfiyen okunduğunda LNZ, 1 sps'lik bir düşey ivme serisi vaat eder.

2026-09-10T14:24Z'deki filo geçişine kadar her örnek, ortalaması çıkarılmış ivmenin üç eksenli vektör büyüklüğünün, √(x²+y²+z²), bir saniyelik tepesiydi. Bu niceliğin yönü yoktur: yer yukarı, yana ya da çapraz hareket etse de aynıdır. Z, verinin tutmadığı bir iddiaydı.

Sebep mekanik. T1 zemine de duvara da monte edilebilir; dolayısıyla sensör paketinin kendi z ekseni bir kurulumda yukarı, ötekinde yana bakar: figürdeki istasyon G274Y duvara monteli ve gövde z ekseni düşeyden 89,74° eğik. Gövde eksenlerine dayanan bir tanım her montajda farklı bir nicelik verirdi. Vektör büyüklüğü yönü tümüyle atarak bundan kaçınıyordu; bu yüzden o da bir düşey kanal değildi.

Hangisi olduğunu cihaz kendisi bildirir. Aynı pencere-ortalama yerçekimi vektörü her sağlık ve tetik paketiyle yayımlanır; gövde z ekseniyle yaptığı açıdan montaj sınıflanır: eksen düşeyden 20° içindeyse düz, 70°–110° arasındaysa duvar, diğer hâllerde eğik. Sınıf her pakette yeniden hesaplanır; cihazın nasıl kurulduğunu kimsenin söylemesi gerekmez, taşınırsa bir sonraki paket onu söyler. Jiroskopun bunda payı yok: LSM6DSO'da jiroskop var, firmware onu kapatıyor.

Algoritma, adım adım

Firmware 0.9.28'den beri cihaz değeri her bir saniyelik pencere için şöyle hesaplıyor:

1. Örnekleme. Kart üzerindeki dört LSM6DSO ivmeölçer yaklaşık 109 Hz'de okunur ve örnek örnek ortalanır: örnek başına tek bir üç eksenli ivme vektörü; döndürme ya da filtre yok.
2. Yerçekimi tahmini. Pencerenin ortalama ivme vektörü, o pencerenin yerçekimi yönü ĝ olarak alınır. Uzun dönemli filtre yoktur; her saniye kendi ĝ'sini kendi yüz küsur örneğinden tahmin eder.
3. Ortalama çıkarma. Pencere ortalaması her örnekten çıkarılır.
4. İzdüşüm. Ortalaması çıkarılmış her örnek ĝ üzerine izdüşürülür: örnek başına işaretli tek bir skaler, yerçekimi çerçevesinde yer ivmesinin düşey bileşeni.
5. Tepe. Penceredeki izdüşümün en büyük mutlak değeri, raporlanan örnektir. İşaret atılır.
6. Çıktı. Saniyede bir örnek, tam sayı count olarak; nominal 1.000.000 count/(m/s²), yani bir count 1 µm/s².

Eksen muhafazadan değil yerçekiminden geldiği için sonuç montaja bağlı değildir: duvardaki cihaz ile raftaki cihaz aynı niceliği hesaplar. Z artık bunu ifade ediyor: her saniye pencere-ortalama yerçekimi vektöründen tahmin edilen, yerçekimi referanslı düşey.

Buna artık "zarf" demiyoruz; kendi belgelerimiz yakın zamana kadar diyordu. Sinyal işleme anlamında zarf — Hilbert ya da RMS — bir dalga formunu sınırlayan pürüzsüz bir fonksiyondur. Buradaki ise saniyede tek bir sayı, o saniyede görülen en büyük |düşey ivme|: bir saniyelik tepe mutlak düşey ivme serisi. Negatif olmaz, fazı yoktur, sıfırı hiç kesmez.

Figür ve okunuşu

Yukarıdaki figür gerçek bir kayıt. Yeni Zelanda, Fiordland'daki GW.G274Y istasyonu 2 Eylül 2026'da 06:13:40 UTC'de, 77 km uzaklıktaki 5,3 büyüklüğünde bir depremle tetiklendi. Camgöbeği iz tetik kaydı ENZ: saniyede yaklaşık 108 örnek, tetikten altı saniye önce, altı saniye sonra, yerçekimi referanslı düşey, işaretli. Beyaz basamak ise aynı izden yukarıdaki algoritmayla hesaplanan bir saniyelik tepe |düşey|.

Köken hakkında iki not. 2 Eylül'de canlı LNZ hâlâ vektör büyüklüğüydü; yani beyaz seri arşivin o dakika için tuttuğu şey değil, sonradan tetik kaydından türetildi — tetik kaydı 5 Ağustos 2026'dan beri yerçekimi çerçevesinde yayınlanıyor. Ayrıca tetik kaydı on iki saniyelik bir pencere; dolayısıyla bu yalnızca P girişi: katalog kaynak zamanı tetikten 13,0 s önce, P artığı +0,05 s, P'den yaklaşık dokuz saniye sonra beklenen S dalgası ise kaydın dışında.

Basamaktan okunacak iki şey. Kayıttaki tepe 0,039 m/s². Ve basamak tetik çizgisinde değil, ondan önceki pencerede yükseliyor: dedektörün ateşlemesi için eşiğin üstünde birkaç örnek gerekir; bu yüzden ilk P enerjisini taşıyan pencere çoktan yüksekken tetik zaman damgası biraz sonraya düşer. LNZ'den bir varış zamanlayacaksanız, tetiğin bir pencereye kadar önüne geçmesini ve bir saniyelik kuantalamayı bekleyin.

Sınırlamalar, açıkça

• Ölçek nominaldir. Count'lar m/s²'ye datasheet duyarlılığıyla çevrilir; cihaz başına kalibrasyon, sarsma masası testi, kutup ve sıfır yoktur, level=response da o nominal değeri döndürür. Genlikler göstergedir, metrolojik olarak izlenebilir değildir.
• Yerçekimi tahmini güçlü sarsıntıda kayabilir. ĝ pencere ortalamasıdır, başka bir şey değil. Zayıf harekette bu iyi bir yerçekimi tahminidir; güçlü harekette, pencere içindeki sarsıntı asimetrikse ortalama yerçekiminden sapar ve izdüşüm ekseni onunla birlikte kayar. Bunu henüz ölçmedik.
• Zamanlama NTP'dir, GNSS değil. Cihaz saati pool.ntp.org'dan SNTP ile saatte bir kurulur; GNSS alıcısı ve PPS yoktur. Alım sunucusu cihaz zaman damgasının yanına kendi alış zamanını da kaydeder.
• Bazı istasyonlar hâlâ eski tanımda. Bu yazının yazıldığı 12 Eylül 2026 itibarıyla listelenen on dört istasyonun on biri yeni seriyi üretiyor. Üçü — 0.9.27'deki QA2HF ile 0.9.14'teki BZAZP ve G8M8U — hâlâ vektör büyüklüğü üretiyor; LNZ kanalları havadan güncellenene kadar StationXML'de TRANSITIONAL yorumu taşıyor.
• 10 Eylül öncesi arşiv eski tanımdır. 2026-09-10T14:24Z'den önce damgalanmış her LNZ örneğini vektör büyüklüğü olarak ele alın. Kanal kodu bu sınırda değişmedi; StationXML'deki geçmiş yorumu tarihi belirtir.

Nasıl okunur

İşlemenizi koda değil StationXML'deki kanal <Description>'ına bağlayın. Sözleşme açıklamadır; bugünkü ifadesi bu yazıdan eski ama tarif ettiği algoritma yukarıdaki; altındaki yorumlar da tarih sınırını ve yerine göre TRANSITIONAL işaretini taşıyor.

from obspy.clients.fdsn import Client

client = Client(base_url="https://api.geoshake.org")
inv = client.get_stations(network="GW", station="G274Y", level="channel")
for cha in inv[0][0]:
    print(cha.code, "-", cha.description)

Canlı akış için:

from obspy.clients.seedlink.easyseedlink import create_client

def on_data(trace):
    print(trace.id, trace.stats.starttime, trace.data.max())

client = create_client("seedlink.geoshake.org:18000", on_data)
client.select_stream("GW", "G274Y", "LNZ")
client.run()

LNZ'yi yataylarla döndürmeyin ve işaretli bir bileşen gibi işlemeyin. İşaretli üç bileşenli küme ENZ, EN1, EN2'dir ve yalnızca istasyon tetiklendiğinde yayınlanır.

Ad üzerine

SEED'e göre üç harf doğru: 1 sps bandı, ivmeölçer, düşey. Ama nicelik türev bir ürün — örneklenmiş bir dalga formu değil, bir pencere üzerindeki tepe — ve SEED'in bunun için harfi yok. Bu kanalı alan bir arşiv kendi geleneğine göre başka bir kodu tercih edebilir; öyle olursa kod değişir, açıklama değişmez.

Ağ GW; kalıcı ve FDSN tescilli, DOI 10.7914/ft1b-4x43; veri merkezi GEOSHAKE. Veriyi kullanıyorsanız ağa FDSN'in verdiği metinle atıf verin:

GeoShake (2026): GeoShake Global Seismic Network. International Federation of Digital Seismograph Networks. Dataset/Seismic Network. https://doi.org/10.7914/ft1b-4x43