Die Physik seismischer Wellen und die S−P-Ankunftszeit
Bei einem Erdbeben breiten sich verschiedene Arten von elastischen Wellen im Erdkörper aus. Die beiden wichtigsten Raumwellen für die Lokalisierung des Epizentrums sind die P-Wellen (Primär- oder Kompressionswellen) und die S-Wellen (Sekundär- oder Scherwellen). Da P-Wellen eine höhere Ausbreitungsgeschwindigkeit als S-Wellen besitzen, treffen sie an einer Registrierstation stets zuerst ein.
Mit zunehmender Entfernung vom Erdbebenherd wächst der zeitliche Abstand zwischen dem Eintreffen der P-Welle und dem Eintreffen der S-Welle kontinuierlich an. Dieses Zeitintervall wird als S−P-Intervall bezeichnet. Unter der Annahme konstanter Wellengeschwindigkeiten in einem homogenen Medium lässt sich die Distanz D zwischen dem Erdbebenherd und einer Messstation direkt aus dieser Zeitdifferenz ableiten. Die mathematische Beziehung lautet:
Δ t = D × (1 / Vₛ - 1 / Vₚ)
Durch Umstellung nach der Distanz ergibt sich die im Rechner verwendete Formel für die Distanznäherung:
D = Δ t ÷ (1 / Vₛ - 1 / Vₚ)
Hierbei stellt Vₚ die P-Wellen-Geschwindigkeit und Vₛ die S-Wellen-Geschwindigkeit dar.
Die Geometrie der Triangulation auf einer Kugeloberfläche
Eine einzelne seismische Station liefert über das S−P-Intervall lediglich eine Information: die Entfernung zum Erdbebenherd. Geometrisch beschreibt diese Entfernung einen Kreis auf der Erdoberfläche, auf dem das Epizentrum liegen muss.
- Eine Station: Der Aufpunkt liegt irgendwo auf einem Kreis mit dem berechneten Radius (Berechneter Radius) um die Station.
- Zwei Stationen: Die zwei Kreise schneiden sich im Normalfall an zwei diskreten Punkten, sodass zwei mögliche Orte für das Epizentrum verbleiben.
- Drei Stationen: Ein dritter Kreis löst diese Zweideutigkeit auf und schneidet die anderen Kreise im Idealfall in einem einzigen, eindeutigen Punkt.
- Vier bis acht Stationen: Zusätzliche Stationen überbestimmen das System. Da reale Messdaten immer mit kleinen Fehlern behaftet sind, schneiden sich die Kreise in der Praxis nicht mehr in einem exakten Punkt. Hier hilft die mathematische Ausgleichsrechnung, um die bestmögliche Position zu bestimmen.
Der Rechner nutzt ein vereinfachtes Modell, bei dem die Erde als Kugel mit einem konstanten Radius R = 6371 km angenommen wird. Die Berechnung der tatsächlichen Oberflächendistanz d zwischen zwei Koordinatenpunkten erfolgt über die Haversine-Formel zur Bestimmung des Zentralwinkels:
Oberflächendistanz: d = R × Haversine-Zentralwinkel, mit R = 6371 km
Methode der kleinsten Quadrate zur Epizentrum-Anpassung
Wenn reale Beobachtungsdaten von drei oder mehr Stationen vorliegen, weichen die berechneten S−P-Distanzen (Dᵢ) meist geringfügig von den tatsächlichen geometrischen Oberflächendistanzen (dᵢ) zum geschätzten Epizentrum ab. Um die optimale Position des Epizentrums zu finden, nutzt das Tool eine mathematische Optimierung nach der Methode der kleinsten Quadrate.
Das Ziel dieser Optimierung ist es, die Summe der quadrierten Abweichungen (Residuen) über alle beteiligten Stationen hinweg zu minimieren:
Epizentrum-Anpassung: Minimierung von Σ (dᵢ - Dᵢ)² über alle Stationen
Das Residuum einer einzelnen Station ist definiert als die Differenz zwischen der angepassten sphärischen Distanz und der aus dem S−P-Intervall abgeleiteten Distanz:
Residuum = Sphärische Distanz - S−P-Distanz
- Ein positives Residuum bedeutet, dass das berechnete Epizentrum weiter von der Station entfernt liegt, als es das gemessene S−P-Intervall nahelegt.
- Ein negatives Residuum bedeutet, dass das Epizentrum näher an der Station liegt, als der berechnete Radius vorgibt.
Die Qualität der Gesamtanpassung wird über das RMS-Distanzresiduum (Root Mean Square) und das größte Distanzresiduum quantifiziert.
Fehlerfortpflanzung und geometrische Einflüsse
Die Genauigkeit der Epizentrumsbestimmung hängt nicht nur von der Präzision der Zeitmessung ab, sondern maßgeblich von der räumlichen Anordnung der Messstationen (Stationsgeometrie).
Zeitliche Unsicherheit und Distanz-1σ
Jede manuelle Bestimmung der Ankunftszeiten von P- und S-Wellen in einem Seismogramm unterliegt einer gewissen Ungenauigkeit. Diese wird im Rechner als zeitliche Unsicherheit (1σ) erfasst. Multipliziert mit dem Distanzfaktor der Wellengeschwindigkeiten ergibt sich daraus die Distanz-1σ pro Station:
Distanz-1σ = zeitliche Unsicherheit × (Vₚ × Vₛ ÷ (Vₚ - Vₛ))
Einfluss der Stationsgeometrie
Die räumliche Verteilung der Stationen bestimmt, wie stark sich diese Distanzunsicherheit auf die berechnete Position auswirkt.
- Starke Stationsgeometrie: Die Stationen umschließen das Epizentrum gleichmäßig aus verschiedenen Himmelsrichtungen. Die Unsicherheit der Position ist gering und in alle Richtungen ähnlich verteilt.
- Schwache Stationsgeometrie: Die Stationen liegen alle in einer Richtung vom Epizentrum aus oder sind nahezu linear angeordnet. In diesem Fall ist die Position entlang der Querachse extrem schlecht eingegrenzt. Das Tool gibt in solchen Fällen eine Warnung aus und berechnet die Unsicherheit entlang der horizontalen Hauptachse.
Vereinfachungen des Lehrmodells gegenüber der Realität
Dieses Tool ist als interaktives Lehrmodell konzipiert, um die grundlegenden geometrischen und mathematischen Prinzipien der Seismologie zu verdeutlichen. Es unterscheidet sich in wesentlichen Punkten von professionellen Systemen zur Erdbebenlokalisierung:
- Homogene Geschwindigkeiten: Das Modell nimmt konstante Geschwindigkeiten für die P- und S-Wellen an. In der realen Erde nehmen die Geschwindigkeiten mit der Tiefe zu, was zu gekrümmten Strahlwegen führt.
- Flache Herdtiefe: Es wird angenommen, dass sich das Erdbeben direkt an der Oberfläche (Herdtiefe ≈ 0 km) ereignet. Reale Erdbeben besitzen eine Herdtiefe (Hypozentrum), die bei der Triangulation als zusätzliche Dimension berücksichtigt werden muss.
- Sphärische Geometrie ohne Schichtung: Professionelle seismische Kataloge nutzen komplexe, dreidimensionale Geschwindigkeitsmodelle der Erdkruste und des Erdmantels, um Laufzeiten präzise zu berechnen.
Bedienung und Dateneingabe
Das Tool erlaubt die Eingabe von 3 bis maximal 8 Stationen.
Eingabeparameter und Grenzwerte
- Breite: Muss für jede Station im Bereich von -90^° bis 90^° liegen.
- Länge: Muss im Bereich von -180^° bis 180^° liegen.
- S−P-Intervall: Muss größer als 0 und darf maximal 2000 s betragen.
- P-Wellen-Geschwindigkeit: Muss größer als 0 und maximal 20 km/s sein sowie über der S-Wellen-Geschwindigkeit liegen.
- S-Wellen-Geschwindigkeit: Muss größer als 0 und maximal 20 km/s sein sowie unter der P-Wellen-Geschwindigkeit liegen.
- Zeitliche Unsicherheit (1σ): Muss zwischen 0 und maximal 60 s liegen.
Über die Schaltfläche Beispiel Tokio laden können vordefinierte Beispieldaten geladen werden, während Leeren alle Eingabefelder zurücksetzt.
Fehlermeldungen bei der Validierung
Sollten Eingaben außerhalb der zulässigen Bereiche liegen oder die Geometrie instabil sein, gibt das Tool präzise Fehlermeldungen aus:
- Bei ungültigen Zeichen:
input: „abc“ ist keine Zahl. - Bei unzureichender Stationsanzahl:
Verwenden Sie zwischen drei und acht Stationen. - Bei identischen Koordinaten:
Die Stationen 1 und 2 haben dieselben Koordinaten. Verwenden Sie unterschiedliche Stationsstandorte. - Bei linearer Anordnung:
Diese Stationen lassen keine stabile zweidimensionale Lokalisierung zu. Fügen Sie eine Station abseits der aktuellen Linie oder Gruppe hinzu.
Datenschutz und Datenverarbeitung
Sämtliche Berechnungen werden direkt im Webbrowser des Nutzers ausgeführt. Es findet keine Übertragung der eingegebenen Stationskoordinaten, S−P-Zeiten oder Berechnungsergebnisse an externe Server statt.
Häufig gestellte Fragen (FAQ)
Ist die S−P-Distanz die exakte Oberflächendistanz?
Nein. Die Gleichung mit konstanter Geschwindigkeit schätzt einen gemeinsamen Laufweg aus der Ankunftszeitdifferenz. Diesen Radius als sphärische Oberflächendistanz zu behandeln, setzt eine geringe Herdtiefe und eine einfache Geschwindigkeitsstruktur voraus. Reale Laufzeitkurven hängen von der Tiefe, der Distanz und dem von jedem Strahl durchquerten Material ab.
Warum werden drei Stationen benötigt?
Eine einzelne S−P-Messung liefert einen Distanzkreis, keine Richtung. Zwei Kreise können sich an zwei Stellen schneiden. Eine dritte Station löst diese Zweideutigkeit im Normalfall auf, während eine vierte oder weitere Station die Abweichung durch Residuen sichtbar macht. Stationen, die um die Quelle herum verteilt sind, grenzen sowohl die Breite als auch die Länge besser ein als Stationen entlang einer einzelnen Linie.
Was bedeutet ein positives oder negatives Residuum?
Das Residuum ist die angepasste sphärische Distanz minus der aus S−P abgeleiteten Distanz. Ein positiver Wert bedeutet, dass der angepasste Punkt weiter von dieser Station entfernt ist als sein S−P-Radius; ein negativer bedeutet, dass er näher liegt. Große, gemischte Residuen zeigen, dass die Kreise bei den gewählten Geschwindigkeiten keinen eindeutigen gemeinsamen Schnittpunkt haben.
Was ist in der Positionsunsicherheit enthalten?
Sie pflanzt den einzelnen zeitlichen 1σ-Wert durch den S−P-Distanzfaktor und die Stationsgeometrie nahe dem angepassten Punkt fort. Sie berücksichtigt kein fehlerhaftes Geschwindigkeitsmodell, keine Herdtiefe, keine falsch zugeordneten Phasen, keinen Uhrzeitversatz und keine dreidimensionale Erdstruktur, weshalb sie keine vollständige Aussage über die Genauigkeit darstellt.