Kalkylator för triangulering av jordbävningsepicentrum

Beräkning av jordbävningsepicentrum genom triangulering

Kalkylatorn för triangulering av jordbävningsepicentrum är ett kostnadsfritt pedagogiskt verktyg som gör det möjligt att uppskatta epicentrum för en jordbävning på en sfärisk jordmodell med hjälp av observationer från tre till åtta registreringsstationer. Genom att mata in koordinater och ankomsttidsintervall för S- och P-vågor (S−P-intervall) för varje station, tillsammans med användardefinierade våghastigheter och tidsosäkerhet, beräknas den bäst anpassade positionen för epicentrum.

Verktyget använder en minsta kvadrat-metod för att minimera avvikelsen i avstånd över samtliga stationer. Resultatet visar det beräknade epicentrumet, styrkan i stationsgeometrin samt individuella residualer för varje station för att illustrera hur väl observationerna stämmer överens.

Fysiken bakom seismiska vågor och triangulering

När en jordbävning inträffar frigörs energi i form av seismiska vågor som sprider sig genom jorden. P-vågor (primära kompressionsvågor) färdas snabbare än S-vågor (sekundära skjuvningsvågor). Detta hastighetsskillnadsförhållande innebär att tidsskillnaden mellan vågornas ankomst till en mätstation ökar ju längre bort från källan stationen befinner sig.

Under antagandet om konstanta hastigheter uppfyller skillnaden i ankomsttid ekvationen:

Δt = D × (1/Vs − 1/Vp)

Vilket kan skrivas om för att beräkna avståndet (D) från stationen till jordbävningen:

D = Δt ÷ (1/Vs − 1/Vp)

En enskild station kan endast fastställa ett avstånd, vilket skapar en cirkel av möjliga positioner. Två stationer begränsar de möjliga positionerna till två skärningspunkter. För att entydigt bestämma ett unikt epicentrum krävs observationer från minst tre stationer. När verkliga mätdata används, som ofta innehåller små felmarginaler, skär cirklarna sällan i en exakt punkt. Därför används en minsta kvadrat-anpassning för att hitta den punkt som minimerar skillnaderna mellan de beräknade avstånden och de faktiska observationerna.

Indata och parametrar i kalkylatorn

Varje station behöver latitud från −90° till 90°, longitud från −180° till 180° och ett S−P-intervall som är större än 0 och högst 2000 sekunder.

För att utföra beräkningen krävs att användaren anger specifika värden för stationerna och vågutbredningen. Följande gränser och valideringsregler tillämpas:

*   **Hastighet för P-wave**: Måste vara större än 0 och högst 20 km/s, samt större än hastigheten för S-vågor. Vid felaktigt värde visas: `Hastighet för P-wave måste vara större än 0 och högst 20 km/s.`
*   **Hastighet för S-wave**: Måste vara större än 0 och högst 20 km/s, samt mindre än hastigheten för P-vågor. Vid felaktigt värde visas: `Hastighet för S-wave måste vara större än 0 och högst 20 km/s.` Om ordningen är felaktig visas: `Hastigheten för P-vågor måste vara större än hastigheten för S-vågor.`
*   **Tidsosäkerhet (1σ)**: Måste vara större än 0 och högst 60 sekunder. Vid felaktigt värde visas: `Tidsosäkerheten måste vara större än 0 och högst 60 sekunder.`

Användaren kan klicka på Läs in exempel för Tokyo för att fylla i kalkylatorn med exempeldata, eller klicka på Rensa för att återställa alla fält.

Beräkningsmetod och matematiska formler

Kalkylatorn modellerar jorden som en perfekt sfär med en konstant radie (R) på 6 371 km. Beräkningen av det sfäriska ytavståndet (d) mellan en station och en punkt på ytan görs med haversine-formeln:

d = R × haversine-centralvinkel, med R = 6371 km

För att hitta det bäst anpassade epicentrumet minimeras summan av de kvadrerade avvikelsena (residualerna) för alla stationer:

Epicentrumanpassning: minimera Σ(dᵢ − Dᵢ)² över alla stationer

Residualen för varje enskild station beräknas som det anpassade sfäriska avståndet minus det avstånd som härletts från S−P-intervallet. En positiv residual innebär att det beräknade epicentrumet ligger längre bort från stationen än dess beräknade S−P-radie, medan en negativ residual innebär att det ligger närmare.

Tolkning av resultat och kvalitet

När beräkningen har genomförts presenterar verktyget resultaten under följande rubriker och etiketter:

• Hero-etikett: Beräknat epicentrum
  • Beräknat epicentrum: Visar koordinaterna för latitud och longitud.
  • tät avståndsanpassning, blandad avståndsanpassning eller dålig avståndsanpassning.
  • stark stationsgeometri, måttlig stationsgeometri eller svag stationsgeometri.
• Schematisk bild: En grafisk representation av stationernas inbördes förhållande visas under rubriken Schematisk bild över stationsgeometri med bildtexten: "Endast relativ geometri — ingen karta, gränser eller live-jordbävningsdata.".
• Stationsanpassning (tabell): Innehåller kolumnerna Station, S−P-avstånd, Sfäriskt avstånd och Residual.
• Formel och insättning: Visar de matematiska stegen, inklusive beräkningen av Avståndsfaktor = 6 × 3.5 ÷ (6 − 3.5) = 8.4 km per sekund av S−P samt de enskilda stationernas härledda avstånd.
• Metod, datum och omfattning: Innehåller texten: "S−P-avståndsekvationen, trestationscirkelmetoden och minsta kvadrat-behandlingen följer de länkade pedagogiska referenserna från USGS. Anpassningen använder sferiskt ytavstånd och en enhetlig hastighetsmodell som ställs in av användaren. Källor och metod kontrollerade 6 augusti 2026.".

Integritet och lokal databehandling

Alla beräkningar sker lokalt i webbläsaren på användarens enhet. Inga koordinater, tider eller beräknade resultat laddas upp till externa servrar.

Vanliga frågor (FAQ)

Varför behövs tre stationer?
En S−P-mätning ger en avståndscirkel, inte en riktning. Två cirklar kan mötas på två ställen. En tredje station avgör i regel valet, medan en fjärde eller senare station gör avvikelsen synlig genom residualer. Stationer utspridda runt källan begränsar både latitud och longitud bättre än stationer längs en enda linje.

Är S−P-avståndet det exakta ytavståndet?
Nej. Ekvationen med konstant hastighet uppskattar en gemensam löpväg utifrån skillnaden i ankomsttid. Att behandla den radien som sferiskt ytavstånd förutsätter en grund källa och en enkel hastighetsstruktur. Verkliga löptidskurvor beror på djup, avstånd och det material som varje stråle passerar.

Vad betyder en positiv eller negativ residual?
Residualen är det anpassade sferiska avståndet minus det avstånd som härletts från S−P. Positivt betyder att den anpassade punkten ligger längre bort från den stationen än dess S−P-radie; negativt betyder att den ligger närmare. Stora blandade residualer visar att cirklarna inte har en ren gemensam skärningspunkt med de valda hastigheterna.

Vad ingår i positionsosäkerheten?
Den propagerar det enskilda tidsosäkerhetsvärdet 1σ genom S−P-avståndsfaktorn och stationsgeometrin nära den anpassade punkten. Den inkluderar inte felaktig hastighetsmodell, fokalpunktens djup, felaktigt identifierade faser, klockavvikelse eller tredimensionell jordstruktur, så det är inte ett fullständigt anspråk på noggrannhet.

Vetenskapliga källor och referenser

• USGS: Jordbävningarnas vetenskap – Grundläggande principer för seismiska vågor och triangulering med tre stationer. Tillgänglig på USGS: Jordbävningarnas vetenskap.
• USGS: Seismiska vågor och lokalisering av jordbävningar – Utbildningsmaterial om S−P-tidsskillnader och minsta kvadrat-anpassning. Tillgänglig på USGS: Seismiska vågor och lokalisering av jordbävningar.
• USGS: Gångtider för seismiska vågor – Information om hur löptider påverkas av jordens skiktade struktur och källans djup. Tillgänglig på USGS: Gångtider för seismiska vågor.