Model interakcji trzęsień ziemi napędzany deformacją dla seizmiczności w Campi Flegrei

Dlaczego ten niespokojny wulkan ma znaczenie

Campi Flegrei, rozległy krater wulkaniczny na zachód od Neapolu, leży pod jednym z najgęściej zaludnionych obszarów Europy. Od ponad wieku grunt tam stopniowo się podnosi i opada, czasami o ponad metr, a od 2005 roku zarówno tempo podnoszenia, jak i liczba drobnych trzęsień ziemi przyspieszyły. Ludzie naturalnie pytają: czy to ostrzeżenie przed nadchodzącą erupcją, czy raczej hałaśliwe „oddychanie” długowiecznego systemu wulkanicznego? To badanie zajmuje się tym pytaniem, budując model oparty na fizyce, który łączy deformacje terenu z interakcjami sejsmicznymi, w celu lepszego oszacowania krótkoterminowego zagrożenia sejsmicznego dla tego obszaru.

Niespokojny grunt Campi Flegrei

Zapisy historyczne i współczesne pomiary pokazują, że kaldera Campi Flegrei przechodziła powtarzające się cykle podnoszenia i zapadania się gruntu przynajmniej od 1905 roku. Główne epizody wypiętrzeń miały miejsce na początku lat 50. XX wieku, około 1970 r., na początku lat 80. oraz ponownie od 2005 r. Zwykle trzęsienia skupiają się podczas tych faz podnoszenia, ale nie w prosty, proporcjonalny sposób: liczba wstrząsów rośnie szybciej niż tempo podnoszenia, a znaczna aktywność sejsmiczna często zaczyna się dopiero, gdy poziom gruntu przekracza poprzedni szczyt. Takie zachowanie przypomina zjawisko znane z mechaniki skał jako efekt Kaisera, w którym naprężony materiał pozostaje cichy, dopóki nie zostanie przekroczone wcześniejsze maksimum naprężenia. Jednak obserwacje w Campi Flegrei są bardziej stopniowe niż ścisły próg „włącz‑wyłącz”, co sugeruje, że działają tu bardziej subtelne procesy fizyczne.

Jak poślizgi na uskokach pamiętają przeszłe naprężenia

Aby uchwycić długoterminowy wzorzec, autorzy wykorzystują ramy zwane tarciem rzędu i stanu (rate‑and‑state friction), które opisują, jak uskoki się ślizgają w zależności zarówno od bieżącego naprężenia, jak i historii obciążenia. W ich uproszczonym modelu naprężenie działające na płytkie uskoki jest traktowane jako proporcjonalne do zmierzonego pionowego podniesienia na stacji GPS wewnątrz kaldery. Takie podejście naturalnie wprowadza pamięć: model rejestruje całe podnoszenie od 1905 roku, więc wcześniejsze epizody inflation wpływają na to, jak łatwo uskoki zaczynają się dziś ślizgać. Przy odpowiednich parametrach model rzędu i stanu odtwarza ogólny trend na przestrzeni stulecia, włącznie z opóźnionym początkiem seizmiczności aż do przekroczenia wcześniejszych szczytów podniesienia. Naśladuje on czasowanie sugerowane przez efekt Kaisera, jednak generuje łagodniejszy, przyspieszający wzrost częstości trzęsień, który lepiej pasuje do obserwacji.

Kiedy trzęsienia wywołują kolejne trzęsienia

Na krótszych skalach czasowych — godzin do dni — zapis sejsmiczny wygląda zupełnie inaczej. Zamiast pojedynczych głównych wstrząsów z uporządkowanymi sekwencjami aftershocków, Campi Flegrei często produkuje gęste stada zdarzeń. Na pierwszy rzut oka te skupiska wydają się pozbawione wyraźnych głównych wstrząsów, ale autorzy pokazują, że wiele z nich zawiera ukryte kaskady aftershocków. Poprzez układanie aktywności wokół największych zdarzeń znajdują, że tempo trzęsień skacze natychmiast po tych wstrząsach, a następnie zanika w czasie w sposób charakterystyczny dla aftershocków. Liczba wywołanych zdarzeń rośnie również szybko wraz z wielkością wstrząsu głównego. Ten wzorzec ujawnia, że interakcje trzęsienie‑trzęsienie są kluczowym składnikiem zachowań stadkowych, nawet jeśli klastry są modulowane przez płyny i inne procesy wulkaniczne.

Hybrydowy obraz naprężeń i klastrowania

Ponieważ sama deformacja nie wyjaśnia intensywnego krótkoterminowego klastrowania, badanie łączy dwa podejścia modelowe. Model rzędu i stanu dostarcza zmiennego w czasie „tła” częstości trzęsień napędzanego podnoszeniem, podczas gdy statystyczny model typu epidemic‑aftershock jest nakładany, by przedstawić, jak każde zdarzenie może wywołać kolejne. Ten hybrydowy model ma siedem parametrów, skalibrowanych przy użyciu tysięcy małych trzęsień zarejestrowanych od 2005 r. Udaje mu się tam, gdzie prostsze alternatywy zawodzą: odtwarza zarówno długoterminowy wzrost seizmiczności, jak i wybuchowe, stadkowe klastry, oraz reprodukuje czasowanie i intensywność minionych epizodów podnoszenia przy symulacjach cofanych do połowy XX wieku. Co istotne, dostarcza też realistycznych oszacowań, jak długo uskoki „pamiętają” wcześniejsze naprężenia.

Co model mówi o ryzyku

Aby przetestować wartość praktyczną, zespół użył modelu hybrydowego w pseudo‑prospektywny sposób: zaczynając od 2020 r., wielokrotnie pytali, jak wyglądałby następny tydzień lub miesiąc pod względem liczby i maksymalnej wielkości trzęsień, używając wyłącznie informacji dostępnych w danym momencie. Tysiące symulowanych scenariuszy dla każdego okna prognozy wygenerowały zakresy prawdopodobieństwa, które w dużej mierze objęły późniejsze obserwacje, w tym zdarzenie o magnitudzie 4,6 w połowie 2025 r. Dla mieszkańców i władz wokół Campi Flegrei nie daje to precyzyjnej prognozy pojedynczego trzęsienia czy erupcji. Zamiast tego oferuje bardziej wiarygodny, oparty na fizyce sposób oszacowania, jak intensywna i jak silna może być aktywność sejsmiczna w nadchodzących tygodniach i miesiącach, poprawiając podstawy krótkoterminowych ocen zagrożenia w tym wrażliwym regionie wulkanicznym.

Cytowanie: Hainzl, S., Dahm, T. & Tramelli, A. A deformation-driven earthquake interaction model for seismicity at Campi Flegrei. Commun Earth Environ 7, 244 (2026). https://doi.org/10.1038/s43247-026-03296-3

Słowa kluczowe: Campi Flegrei, trzęsienia wulkaniczne, podniesienie terenu, prognozowanie sejsmiczne, stadka trzęsień