Non-Extensive statistical physics properties of seismic coda waves recorded during the Santorini unrest period 2011-2012

Abstract

Η περίοδος της ηφαιστειακής κρίσης στο νησί της Σαντορίνης κατά τη περίοδο 2011- 2012, οδήγησε αρκετές επιστημονικές μελέτες να εγείρουν το κρίσιμο ερώτημα, εάν η σεισμικότητα της ευρύτερης περιοχής της καλντέρας, σχετίζεται με μια επικείμενη και πιθανή ηφαιστειακή έκρηξη ή αν εκφορτίζει τη συσσωρευμένη τεκτονική ενέργεια. Στην παρούσα μεταπτυχιακή διατριβή, μελετήσαμε τα σεισμικά κύματα ουράς (Coda waves), από τα τοπικά σεισμικά γεγονότα πριν, κατά τη διάρκεια και μετά την περίοδο της ηφαιστειακή και σεισμικής κρίσης της Σαντορίνης. Τα σεισμικά κύματα ουράς, οφείλουν την παρουσία τους σε φαινόμενα σκέδασης και στην πολυπλοκότητα της κατανομής των σκεδαστών του υπεδάφους, στις ανομοιογένειες του εσωτερικού της Γης, όπως είναι οι ρωγμές και η σχιστότητα των πετρωμάτων, τα ρήγματα, στα όρια μετάβασης υλικών διαφορετικής πυκνότητας, και στις πάσης φύσεως γεωφυσικές ανομοιογένειες ως προς τις ιδιότητες των υλικών στα ανώτερα στρώματα του φλοιού της Γης. Πιο συγκεκριμένα, λάβαμε υπόψη και τις τρεις συνιστώσες των σεισμολογικών σταθμών στο νησί της Σαντορίνης και εκτιμήσαμε τη διάρκεια των σεισμικών κυμάτων ουράς, εφαρμόζοντας μια διαδικασία πέντε επιπέδων, η οποία περιλαμβάνει την αναλογία σεισμικού σήματος προς θόρυβο (SNR), τη μέθοδο τύπου βραχέως μέσου χρόνου προς μακροπρόθεσμο μέσο χρόνο (STA/LTA) στο πεδίο του χρόνου αλλά και στο αναλυτικό σήμα, όπως αυτό προκύπτει από την εφαρμογή του μετασχηματισμού Hilbert, καθώς και τον μετασχηματισμό Fourier βραχέως χρόνου (STFT), εφαρμόζοντας παράθυρο Γκαουσιανής συνάρτησης (window function). Η τελική εκτίμηση της διάρκειας των κυμάτων ουράς, επαληθεύτηκε ή και επανεκτιμήθηκε αναλογικά από τον χρήστη, εξαιτίας της ύπαρξης ισχυρού σεισμικού θορύβου. Λόγω της πολυπλοκότητας των κυμάτων ουράς, για τη μελέτη του ηφαιστειακού συμπλέγματος της Σαντορίνης, εφαρμόσαμε την θεωρία της μη εκτατικής στατιστικής φυσικής (NESP). Η προσέγγιση NESP εκφράζει την γενίκευση της στατιστικής μηχανικής των Boltzmann-Gibbs και έχει χρησιμοποιηθεί εκτενώς για την ανάλυση πολύπλοκων συστημάτων, με εμφανή χαρακτηριστικά και ιδιότητες αλληλεπίδρασης μακράς εμβέλειας, μνήμης και μορφοκλασματικότητας. Η ανάλυση και τα αποτελέσματα της παρούσας διατριβής δείχνουν ότι οι διαφορές πλάτους των διαδοχικών σημείων των σεισμικών κυμάτων ουράς (στο πεδίο του χρόνου), αποκλίνουν από το στατιστικό πρότυπο τύπου Gauss, καθώς οι αντίστοιχες συναρτήσεις πυκνότητας πιθανότητας (PDFs), μπορούν να περιγραφούν επαρκώς από κατανομές πιθανότητας που προκύπτουν από τον μη εκτατικό φορμαλισμό και πιο συγκεκριμένα από την κατανομή q-Gaussian. Επιπλέον, για να διερευνήσουμε τη δυναμική εξέλιξη της ηφαιστειακής-σεισμοτεκτονικής δραστηριότητας, εκτιμήσαμε τους εντροπικούς δείκτες σκέδασης qs όπως αυτοί προκύπτουν από την μη γραμμική ανάλυση των PDFs, παρουσιάζοντας τις μεταβολές τους ως συνάρτηση του χρόνου, του χώρου (χωροχρονική μεταβολή), του τοπικού μεγέθους (ML) καθώς και της επικεντρικής απόστασης (km) πριν, κατά τη διάρκεια και μετά την περίοδο ηφαιστειακής κρίσης της καλντέρας, από το 2009 έως το 2014.The unrest period in Santorini's caldera during 2011-2012, led several studies to raise the important question of whether seismicity is associated with an impending and potential volcanic eruption or if it solely relieves the accumulated tectonic energy. In the present work, we study seismic coda waves, generated by local earthquake events prior, during and after the Santorini’s unrest period (volcanic and seismic crisis) that occurred within the caldera area. Coda waves are interpreted as scattered seismic waves generated by heterogeneities within the Earth, i.e. by faults, fractures, fissure structures, microcracks, velocity and/or density boundaries/anomalies, etc. In particular, we utilize all three components of the seismograms recorded by three seismological stations on the Santorini island and estimate the duration of the coda waves by implementing a five-step procedure that includes the signal-to-noise ratio, the STA/LTA like method, in time domain and upon analytic signal obtained after implementing Hilbert transformation, and the short time Fourier transform, using gaussian windowing as a window function. The final estimation was verified or re-estimated manually due to the existent ambient seismic noise. Because of the nature and the path complexity of the coda waves and towards achieving a unified framework for the study of the immerse geo-structural seismotectonic complexity of the Santorini volcanic complex, we use Non-Extensive Statistical Physics (NESP) to study the probability density functions (PDFs) of the increments of seismic coda waves. NESP, forms a generalization of the Boltzmann-Gibbs statistical mechanics that has been extensively used for the analysis of semi-chaotic systems that exhibit long-range interactions, memory effects and multi-fractality. The analysis and results demonstrate that the seismic coda waves increments, deviate from the Gaussian shape and their respective pdfs could adequately be described and processed by the q-Gaussian distribution. Furthermore and in order to investigate the dynamical evolution of the volcanic-tectonic activity, we estimate the qs indices derived from the PDFs of the coda wave time series increments, and present their variations as a function of time, of space (spaciotemporal variation), of local magnitude (ML) and epicentral distance (measured in km) prior, during and after the caldera unrest period from 2009 to 2014

    Similar works