La fisica delle onde sismiche e l'intervallo S−P
La localizzazione dell'epicentro di un terremoto si basa sulle proprietà fisiche di propagazione delle onde sismiche all'interno della Terra. Quando si verifica una rottura lungo una faglia, l'energia si libera sotto forma di onde di volume che si propagano in tutte le direzioni. Queste onde si dividono principalmente in due categorie:
- Onde P (primarie o compressionali): sono onde longitudinali che comprimono e dilatano la roccia nella stessa direzione di propagazione. Essendo le più veloci, vengono registrate per prime dai sismografi.
- Onde S (secondarie o di taglio): sono onde trasversali che fanno oscillare le particelle di roccia perpendicolarmente alla direzione di propagazione. Viaggiano a una velocità inferiore rispetto alle onde P e arrivano alla stazione di registrazione con un ritardo temporale.
Poiché le onde P viaggiano più velocemente delle onde S, l'intervallo di tempo che intercorre tra il loro arrivo (chiamato intervallo S−P, o Δ t) aumenta progressivamente con la distanza dalla sorgente sismica. Sotto l'ipotesi di velocità costanti all'interno di un mezzo uniforme, la relazione tra la distanza della stazione dall'epicentro (D) e l'intervallo temporale soddisfa l'equazione Δ t = D(1/Vs - 1/Vp). Questa formula può essere riscritta per ricavare direttamente la distanza stimata:
D = Δ t ÷ (1/Vs - 1/Vp)
Il Calcolatore triangolazione epicentro terremoto utilizza questa relazione fondamentale per convertire i tempi di arrivo inseriti in raggi di distanza sferica per ciascuna stazione di rilevamento.
La geometria della triangolazione sismica
Una singola stazione sismica che registra un intervallo S−P permette di calcolare esclusivamente la distanza radiale dall'evento, definendo una circonferenza di possibili posizioni epicentrali attorno alla stazione stessa.
- Una stazione: fornisce unicamente un raggio di distanza, lasciando l'epicentro indefinito lungo un intero cerchio.
- Due stazioni: i due cerchi di distanza si intersecano, nella maggior parte dei casi, in due punti distinti sulla superficie terrestre, non consentendo una risoluzione univoca.
- Tre stazioni: l'intersezione dei tre cerchi permette di risolvere l'ambiguità e identificare un unico punto nello spazio bidimensionale.
- Da quattro a otto stazioni: l'aggiunta di ulteriori stazioni introduce una sovraddeterminazione del sistema. Nella realtà, a causa di incertezze strumentali e variazioni geologiche, i cerchi non si intersecheranno mai in un unico punto perfetto. Questa discrepanza rende visibili i residui di distanza e richiede un approccio di ottimizzazione matematica.
La disposizione spaziale delle stazioni, definita come geometria delle stazioni, gioca un ruolo cruciale nella stabilità della soluzione. Se le stazioni sono disposte in modo uniforme attorno all'epicentro, la posizione risulta fortemente vincolata sia in latitudine che in longitudine. Al contrario, se le stazioni sono allineate lungo un'unica linea o concentrate in un gruppo ristretto, la geometria non è in grado di vincolare una posizione bidimensionale stabile.
Ottimizzazione dei minimi quadrati e calcolo sferico
Per risolvere la posizione dell'epicentro in presenza di dati reali e imperfetti, lo strumento applica un trattamento dei minimi quadrati. Questo metodo matematico non cerca un'intersezione geometrica perfetta, ma identifica le coordinate di latitudine e longitudine che minimizzano la discrepanza complessiva tra le distanze calcolate dai tempi di arrivo e le distanze sferiche effettive sulla superficie terrestre.
La funzione obiettivo da minimizzare su tutte le stazioni è:
minimizza Σ (dᵢ - Dᵢ)²
Dove dᵢ rappresenta la distanza superficiale sferica tra l'epicentro stimato e la stazione i-esima, mentre Dᵢ è la distanza teorica dedotta dall'intervallo S−P della medesima stazione.
La distanza superficiale sferica d viene calcolata utilizzando la formula dell'haversine per l'angolo centrale su un modello sferico della Terra con un raggio costante R = 6371 km:
d = R × angolo centrale haversine
Attraverso questo processo iterativo, l'algoritmo calcola la soluzione ottimale e restituisce i parametri di accuratezza dell'adattamento, tra cui il residuo RMS (Root Mean Square) e il residuo di distanza maggiore.
Funzionamento del calcolatore e parametri di input
L'interfaccia del calcolatore richiede l'inserimento di parametri fisici e osservazioni dirette per eseguire la triangolazione. L'utente può inserire da un minimo di 3 a un massimo di 8 stazioni.
Parametri del modello di propagazione
- Velocità P-wave ($Vp$): deve essere maggiore di 0 e non superiore a 20 km/s, e deve essere strettamente maggiore della velocità S-wave.
- Velocità S-wave ($Vs$): deve essere maggiore di 0 e non superiore a 20 km/s, e deve essere inferiore alla velocità P-wave.
- Incertezza temporale (1σ): l'errore standard associato alla lettura dei tempi di arrivo sui sismogrammi. Deve essere maggiore di 0 e non superiore a 60 secondi.
Osservazioni delle stazioni
Per ciascuna stazione inserita, l'utente deve specificare:
- Nome stazione
- Latitudine: compresa tra -90^° e 90^°.
- Longitudine: compresa tra -180^° e 180^°.
- Intervallo S−P: il tempo in secondi tra l'arrivo dell'onda P e dell'onda S. Deve essere maggiore di 0 e non superiore a 2000 secondi.
L'utente può utilizzare il pulsante Carica esempio di Tokyo per popolare istantaneamente i campi con un set di dati predefinito, oppure cliccare su Cancella per azzerare tutti i moduli.
Analisi dei risultati e messaggi di errore
Una volta inseriti i dati, il calcolatore elabora le informazioni localmente nel browser dell'utente; le coordinate delle stazioni, i tempi S−P e i risultati rimangono su questo dispositivo e non vengono caricati su server esterni.
Sezione dei risultati
I dati elaborati vengono mostrati nelle seguenti aree dell'interfaccia:
- Soluzione dell'epicentro: mostra le coordinate dell'epicentro stimato.
- Qualità dell'adattamento: l'adattamento di distanza viene classificato come adattamento di distanza stretto, adattamento di distanza misto o adattamento di distanza scarso.
- Qualità della geometria: la disposizione spaziale viene valutata come geometria delle stazioni forte, geometria delle stazioni moderata o geometria delle stazioni debole.
- Metriche di errore: vengono mostrati il Residuo di distanza RMS, il Residuo di distanza maggiore, l'Incertezza di posizione 1σ approssimativa e la Distanza 1σ per stazione.
- Schema della geometria delle stazioni: una rappresentazione grafica semplificata accompagnata dalla didascalia: "Solo geometria relativa — nessuna mappa, confine o dato sismico in tempo reale."
- Adattamento stazioni: una tabella dettagliata con le colonne Stazione, Distanza S−P, Distanza sferica e Residuo.
Messaggi di stato e validazione
A seconda dello stato dei dati, il sistema mostra uno dei seguenti messaggi nella sezione Metodo, data e ambito (verificato il 6 agosto 2026):
- "Modifica le osservazioni delle stazioni per riadattare l'epicentro."
- "Adattate
3stazioni su un modello sferico della Terra." - "Adattate
3stazioni, ma i residui o la geometria delle stazioni richiedono una revisione."
In caso di inserimenti non validi, l'interfaccia blocca il calcolo e mostra errori specifici:
- Numero di stazioni errato:
Usa tra tre e otto stazioni. - Valore non numerico:
input: "abc" non è un numero. - Latitudine fuori limite:
Stazione 1: la latitudine deve essere compresa tra −90° e 90°. - Longitudine fuori limite:
Stazione 1: la longitudine deve essere compresa tra −180° e 180°. - Intervallo S−P non valido:
Stazione 1: inserisci un intervallo S−P maggiore di 0 e non superiore a 2000 secondi. - Velocità fuori limite:
Vₚ/Vₛ deve essere maggiore di 0 e non superiore a 20 km/s. - Ordine velocità errato:
La velocità P-wave deve essere maggiore della velocità S-wave. - Incertezza temporale fuori limite:
L'incertezza temporale deve essere maggiore di 0 e non superiore a 60 secondi. - Coordinate duplicate:
Le stazioni 1 e 2 hanno le stesse coordinate. Usa posizioni di stazione distinte. - Geometria instabile:
Queste stazioni non vincolano una posizione bidimensionale stabile. Aggiungi una stazione lontana dalla linea o dal gruppo attuale. - Overflow dei calcoli:
Un valore o un risultato intermedio supera l'intervallo numerico supportato.
Semplificazioni del modello didattico rispetto alla sismologia professionale
Questo calcolatore è uno strumento didattico progettato per illustrare i principi geometrici e matematici della localizzazione sismica. Per mantenere il modello accessibile, vengono applicate alcune importanti semplificazioni fisiche rispetto ai sistemi di calcolo professionali utilizzati dai centri di monitoraggio sismico:
| Caratteristica | Modello del Calcolatore | Sismologia Professionale |
|---|---|---|
| Forma della Terra | Sfera perfetta con raggio costante R = 6371 km | Ellissoide di riferimento (es. WGS84) |
| Struttura della velocità | Uniforme e costante ($Vp$ e $Vs$ fisse) | Modelli di velocità stratificati (1D) o tridimensionali (3D) |
| Profondità ipocentrale | Assunta pari a zero (sorgente superficiale) | Calcolo della profondità focale (ipocentro) |
| Traiettoria dei raggi | Lineare lungo la superficie sferica | Raggi curvi rifratti e riflessi attraverso gli strati terrestri |
| Incertezza | Propagazione dell'errore temporale inserito (1σ) | Matrice di covarianza completa comprendente errori di modello e di lettura |
Domande frequenti (FAQ)
Perché sono necessarie tre stazioni?
Una singola misurazione S−P fornisce un cerchio di distanza, non una direzione. Due cerchi possono incontrarsi in due punti. Una terza stazione normalmente risolve la scelta, mentre una quarta o successiva stazione rende visibile la discrepanza attraverso i residui. Le stazioni distribuite intorno alla sorgente vincolano sia la latitudine che la longitudine meglio delle stazioni allineate lungo un'unica linea.
La distanza S−P è l'esatta distanza superficiale?
No. L'equazione a velocità costante stima un percorso di propagazione comune a partire dalla differenza dei tempi di arrivo. Trattare tale raggio come distanza superficiale sferica presuppone una sorgente superficiale e una struttura di velocità semplice. Le curve dei tempi di viaggio reali dipendono dalla profondità, dalla distanza e dal materiale attraversato da ciascun raggio.
Cosa significa un residuo positivo o negativo?
Il residuo è la distanza sferica adattata meno la distanza dedotta da S−P. Un valore positivo significa che il punto adattato è più lontano da quella stazione rispetto al suo raggio S−P; un valore negativo significa che è più vicino. Residui ampi e misti indicano che i cerchi non condividono un'unica intersezione netta con le velocità scelte.
Cosa è incluso nell'incertezza della posizione?
Propaga il singolo valore di incertezza temporale 1σ attraverso il fattore di distanza S−P e la geometria delle stazioni vicino al punto adattato. Non include un modello di velocità errato, la profondità focale, fasi identificate erroneamente, errori sistematici degli orologi o la struttura tridimensionale della Terra, quindi non costituisce una dichiarazione di accuratezza completa.