Fractional-Order Viscoelastic Modeling of Seismic P-, S-, and Rayleigh-Wave Attenuation: Marmara Seismotectonic Context, Fault Geometry, and a Reproducible Strong-Motion Framework

Profondamente sotto la superficie della Terra, il suolo non si comporta come un semplice blocco solido di roccia. Quando un terremoto colpisce, l'energia si propaga verso l'esterno sotto forma di onde, ma mentre queste onde si muovono attraverso la crosta, perdono energia e cambiano forma. Questo processo, noto come attenuazione, è il motivo per cui un terreno che trema appare diverso a una certa distanza dalla sorgente rispetto a quanto lo sia proprio accanto alla rottura. Per decenni, gli scienziati hanno cercato di descrivere questa perdita di energia utilizzando le regole standard della fisica, trattando la crosta terrestre come una molla perfettamente elastica o come un semplice fluido viscoso. Tuttavia, il mondo reale è spesso più complesso di quanto questi due estremi permettano. Il terreno è un miscuglio disordinato di roccia, fluido e crepe che ricorda i propri movimenti passati in un modo che la fisica standard fatica a catturare. È qui che una branca della matematica chiamata calcolo frazionario diventa utile. Invece di usare numeri interi per descrivere come i materiali reagiscono, il calcolo frazionario permette di utilizzare valori intermedi, offrendo un modo per modellare la "memoria" della roccia con molte meno assunzioni.

In uno studio recente focalizzato sul Mar di Marmara, una regione di intensa attività geologica appena a sud di Istanbul, i ricercatori Taylan Demir e Atakan Koçyiğit hanno costruito un quadro matematico più preciso per descrivere come le onde sismiche svaniscano. Il Mar di Marmara si trova sopra la Faglia Nord Anatolica, una massiccia frattura nella crosta terrestre dove le placche tettoniche scorrono l'una accanto all'altra. Poiché questa faglia corre direttamente sotto una grande città, comprendere esattamente come le onde sismiche viaggiano e si indeboliscono in quest'area specifica è fondamentale per progettare edifici più sicuri e prevedere l'intensità delle future scosse. I ricercatori si sono posti l'obiettivo di testare se un modello frazionario, che tiene conto del comportamento complesso e dipendente dalla storia della roccia, potesse spiegare lo svanire delle onde sismiche meglio dei modelli tradizionali. Non hanno tentato di prevedere quando avverrà il prossimo grande terremoto; si sono invece concentrati su come ricostruire accuratamente le scosse che si verificano dopo che un evento è stato registrato.

Il team ha iniziato mappando la scena geologica. Hanno esaminato i segmenti specifici della faglia che corrono sotto il mare, inclusi i rami di Tekirdağ, del Bacino Centrale e di Kumburgaz, e hanno analizzato la storia dei terremoti nella regione, come il grande terremoto di İzmit del 1999 e il più piccolo evento di Silivri del 2019. Questo contesto geologico era essenziale perché la matematica che stavano sviluppando doveva essere radicata nella forma e nella struttura reale della zona di faglia. Hanno poi introdotto un nuovo modo di descrivere il comportamento della roccia. Immaginate la crosta terrestre non come una semplice molla che scatta istantaneamente, ma come un materiale che si rilassa lentamente nel tempo, conservando la memoria di come è stato teso o compresso. I ricercatori hanno utilizzato uno strumento matematico chiamato derivata di Caputo per descrivere questo effetto di memoria. Questo strumento ha permesso loro di creare un modello che connette le proprietà fisiche della roccia al modo in cui le onde perdono energia durante il loro percorso.

Utilizzando questo nuovo modello, i ricercatori hanno derivato come dovrebbero comportarsi tre diversi tipi di onde sismiche: le onde P, che sono le più veloci e arrivano per prime; le onde S, che scuotono il terreno lateralmente; e le onde di Rayleigh, che rotolano lungo la superficie. Hanno dimostrato matematicamente che il loro modello produce sempre una perdita di energia, che è una condizione necessaria per qualsiasi descrizione realistica di un terremoto. Se un modello suggerisse che le onde possano guadagnare energia durante il percorso, sarebbe fisicamente impossibile. Le loro equazioni hanno dimostrato che l'approccio frazionario produce una curva di perdita di energia fluida e continua che corrisponde al modo in cui le rocce reali si comportano attraverso un'ampia gamma di frequenze. Questo rappresenta un miglioramento significativo rispetto ai modelli più vecchi, che spesso faticano a adattarsi ai dati attraverso diverse frequenze senza diventare eccessivamente complicati.

Per testare se questa teoria funzioni effettivamente, il team ha condotto un esperimento controllato utilizzando dati generati al computer. Hanno creato uno scenario di terremoto sintetico in cui conoscevano esattamente le regole di come le onde dovrebbero svanire. Hanno poi chiesto al loro modello frazionario di individuare quelle regole basandosi solo sui pattern delle onde. I risultati sono stati sorprendenti. Il modello frazionario è riuscito a recuperare i parametri corretti con un'accuratezza molto elevata, riducendo l'errore nelle sue previsioni di quasi l'87% rispetto a un modello tradizionale che utilizzava regole a numeri interi. Nel modello tradizionale, i ricercatori hanno dovuto forzare i numeri per farli coincidere, ottenendo un adattamento scarso rispetto al reale comportamento delle onde. Il modello frazionario, permettendo quell'effetto di memoria "intermedio", ha catturato la fisica della perdita d'onda in modo molto più naturale.

Lo studio ha anche esaminato come questo framework possa essere applicato ai dati reali della regione del Marmara. I ricercatori hanno identificato il terremoto di Silivri del 2019 come un candidato perfetto per un primo test reale. Hanno delineato un piano chiaro su come gli ingegneri e i sismologi potrebbero scaricare le registrazioni effettive di questo evento e utilizzare le loro nuove equazioni per ricostruire l'intensità delle scosse. Ciò permetterebbe di calcolare misure precise del movimento del suolo, come l'accelerazione o la velocità massima, che sono i numeri utilizzati dagli architetti per progettare edifici resistenti ai terremoti. Gli autori sono stati cauti nel dichiarare che questo strumento serve ad analizzare ciò che è già accaduto, non a prevedere il futuro. Non dice quando avverrà la prossima grande rottura, ma fornisce una lente molto più nitida per comprendere la violenza delle scosse una volta che queste accadono.

Combinando una rigorosa prova matematica con una chiara mappa geologica e un piano per l'analisi dei dati reali, questo articolo sposta la conversazione dalla teoria astratta all'applicazione pratica. Dimostra che l'approccio frazionario non è solo una curiosità matematica, ma un modo valido e più accurato per descrivere la complessa realtà delle onde sismiche nel Mar di Marmara. Il lavoro colma il divario tra la geologia disordinata e stratificata della zona di faglia e le equazioni pulite e prevedibili necessarie per la sicurezza ingegneristica. Sebbene la calibrazione completa rispetto a migliaia di veri record sismici sia un compito per il futuro, questo studio ha gettato le basi necessarie, provando che il metodo funziona in teoria e in simulazione, ed è pronto per essere testato contro il terreno reale sotto Istanbul.