Assetto strutturale e tassi di deformazione lungo il margine esterno dell'Appennino Settentrionale, con implicazioni per la migrazione di fluidi profondi

Riassunto Scopo principale di questa tesi di dottorato è l’analisi dell’assetto strutturale del margine esterno dell’Appennino Settentrionale Emiliano-Romagnolo e Marchigiano, nonché lo studio delle relazioni che intercorrono tra le strutture ivi presenti e la migrazione di fluidi profondi. Tale fenomeno, nell’area di studio comporta la formazione di cosiddetti vulcanelli di fango (mud volcanoes), ben conosciuti ed investigati nel settore Emiliano-Romagnolo, tuttavia meno studiati in quello Marchigiano. Inoltre, tale studio ha visto l’utilizzo di dataset sismici 3D, provenienti da altri contesti tettonici, utili per l’investigazione dettagliata delle strutture, geometrie, meccanismi di messa in posto e processi evolutivi di analoghi dei vulcanelli di fango. Il margine Appenninico Emiliano-Romagnolo è stato investigato con un approccio multidisciplinare, integrando dati strutturali, sedimentologici, pedologici e sismici. In aggiunta, si è proceduto, per due transetti individuati attraverso di esso, ad impostare una modellizzazione numerica -effettuata applicando il metodo deformativo del Trishear- atta a calcolare i tassi di deformazione della principale struttura del margine, il cosiddetto Thrust Pedeappenninico (PAT). L’analisi ha messo in evidenza come il settore studiato (compreso tra le Valli del Fiume Panaro e del Fiume Enza) sia interessato da un’attività tettonica recente. Il PAT infatti, spesso emergente o sub-emergente, deforma depositi Medio Pleistocenici. Tali depositi, investigati da un punto di vista sedimentologico, hanno messo in evidenza come le strutture del margine abbiano condizionato non solo l’evoluzione tettonica del margine stesso, ma anche il reticolo idrografico dei paleo-corsi d’acqua. Si è messo in evidenza infatti come durante il Pleistocene Medio questi fossero orientati parallelamente al margine Appenninico, almeno per alcuni tratti. Tale orientamento si riscontra anche in depositi estremamente recenti, e localizzati a fronte del margine in corrispondenza del Plateau di Ghiardo. Datazioni OSL di suoli ivi deformati hanno permesso di individuare fasi tettoniche relative a circa 60-80 kyr. La modellizzazione numerica ha permesso di calcolare i tassi di deformazione del PAT, ed in particolare gli slip rate, che sono risultati essere compresi nell’ordine di 0,68-0,79 mm/anno per gli ultimi 0,8-1,2 Ma. Nello stesso settore, e in quello immediatamente circostante la città di Bologna, si è dunque provveduto ad investigare le relazioni intercorrenti tra le strutture del margine e i fluidi responsabili della formazione dei mud volcanoes. Sezioni sismiche hanno permesso di mappare le strutture del sottosuolo e constatare la corrispondenza tra anticlinali e vulcanelli. In particolare, uno studio di dettaglio sul vulcanello denominato Dragone di Sassuno, ha messo in evidenza questa dipendenza, nonché l’ulteriore relazione intercorrente tra strutture e fluidi: questi, nella risalita verso la superficie, utilizzano infatti i sistemi di frattura associati all’anticlinale sul quale il vulcanello giace, e che probabilmente svolge in profondità il ruolo di trappola per i fluidi stessi. Come ulteriore approfondimento della tematica, si è provveduto ad indagare possibili correlazioni tra l’accadimento di terremoti storici e l’attivazione di tale vulcanello di fango. È stato possibile modellizzare le variazioni degli stress statici associati a tre eventi sismici di M> 4 occorsi tra il 1779 ed il 1780, in concomitanza di forti eruzioni del vulcanello di Sassuno, descritte dettagliatamente dall’Autore Serafino Calindri al tempo dell’eruzione. In questo caso si è osservato che tale variazione degli stress non può aver avuto un’influenza positiva sulle attivazioni, essendo i valori ottenuti negativi. Tali valori indicano infatti condizioni sfavorevoli all’apertura del complesso di fratture costituenti il feeder dyke system del vulcanello. Si è ipotizzato quindi, data la corrispondenza temporale tra eventi sismici ed attivazione del vulcanello, una possibile dipendenza di tali attivazioni dagli stress dinamici, transitori ed associati al passaggio delle onde sismiche. Vulcanelli di fango sono stati individuati ed investigati anche lungo il Margine Apenninico Marchigiano, con particolare focus su tre aree chiave: S.Paolo di Jesi, Montelone di Fermo e Offida. Anche in questo caso si è provveduto a caratterizzarne l’assetto strutturale con rilievo di campagna ed utilizzo di sezioni sismiche. Come nel caso dei vulcanelli emiliani, si è riscontrata una corrispondenza tra anticlinali -in questo caso sepolte-, strutture ad esse associate e vulcanelli di fango. Per i vulcanelli di Monteleone di Fermo e di Offida, si è provveduto ad indagare la possibilità di attivazione a seguito di triggering sismico. Un caso particolarmente interessante è rappresentato dall’attivazione dei vulcanelli di Monteleone di Fermo a seguito degli eventi sismici dell’Italia Centrale, verificatisi a partire dal 24 Agosto 2016. Si è potuto verificare come i valori degli stress dinamici calcolati in corrispondenza dei vulcanelli attivati siano sufficientemente elevati (>4 bar) da essere considerati la principale causa di attivazione, considerando anche i non sufficienti valori degli stress statici (<0.1 bar). In questo caso, le relazioni intercorrenti tra strutture e migrazione di fluidi è duplice. Una relazione “diretta” tra strutture (anticlinali) e vulcanelli, che ne condiziona il posizionamento, ed una “indiretta” tra attivazione di strutture sismogenetiche distanti dai vulcanelli (>70-80 km) e l’innesco delle eruzioni. Come ultima analisi si è provveduto ad utilizzare due volumi sismici 3D provenienti dal Loyal Field (Scozia) e dal Bacino di Ceará (Brasile). Il primo dataset ha messo in evidenza fondamentali punti di similarità tra l’evoluzione delle cosiddette “fluid escape pipes” e i vulcanelli di fango. Entrambe le tipologie di oggetti geologici sono strettamente dipendenti dall’assetto strutturale delle anticlinali che ne ospitano i reservoir e dalle geometrie dei sistemi di frattura che favoriscono la risalita dei fluidi verso la superficie. In aggiunta, il fluido responsabile della formazione delle fluid escape pipes del Loyal field è risultato essere composto da una mistura contenente fango-acqua-idrocarburi, come nel caso dei vulcanelli di fango italiani. Ulteriore analogia risiede nell’evoluzione dell’attività delle fluid escape pipes e dei vulcanelli, entrambi alternanti periodi di quiescenza a periodi di intensa attività eruttiva. Nel caso del Bacino di Ceará, la presenza sul fondale marino di oggetti assimilabili a pockmarks, ovvero zone depresse indicanti l’emissione di fluidi, è risultata essere uno stimolo per l’analisi del dataset, nell’ottica di un’ulteriore analisi dei processi che guidano la migrazione dei fluidi stessi. Lo studio ha però messo in evidenza come queste depressioni siano in realtà associate a forme di fondo (sediment waves), impilate in maniera tale da formare condotti di migrazione preferenziale, tuttavia passiva, per i fluidi intrappolati nel sottostante reservoir, in maniera indipendente da strutture tettoniche. Abstract The principal aim of this Ph.D. thesis was the investigation of the structural setting of the external Northern Apennine margin, in the Emilia-Romagna and Marche Region, as well as the analysis of possible relations linking the structure of the margin and the migration of deep fluids. This last phenomenon is responsible for the formation of the so-called mud volcanoes, which are of particular interest for our study and that have been investigated with a multidisciplinary approach. First, I investigated the structural setting of the Emilia-Romagna pede-Apennine margin, including possible relations with the mud volcanoes of the area. The same study approach was used for the Marche foothills, and for their associated mud volcanoes. Furthermore, I studied two 3D seismic dataset from different tectonic contexts, to be used as analogues for investigation of mud volcanoes and -more in general- fluid migration processes. The Emilia Romagna Pede-Apennine Margin has been analysed using a multidisciplinary approach, integrating structural, sedimentological, pedological and seismic data. In addition, I carried, for two target transects crossing the margin (at Quattro Castella and Scandiano localities), cinematic 2D cross section numerical modelling of the main structure, the so-called Pede-Apennine Thrust (PAT). The implemented Trishear deformation mechanism, has shown how the studied sector (between the Enza and Panaro rivers) has been tectonically active in recent times. The PAT is poorly exposed and deforms at least Middle Pleistocene deposits. These semi-consolidated sediments, investigated from a sedimentological point of view, revealed that the PAT and its related structures influenced also the river drainage pattern, forcing paleo-rivers to flow parallel to the margin (which roughly strikes NW-SE), at least in some sectors during the Pleistocene. This trend is also visible in extremely recent deposits, outcropping in the vicinity of the margin at the Ghiardo Plateau, that have been dated using the OSL technique applied to deformed palaeosoils. Dating highlighted an extremely recent tectonic phase (~60-80 kyr). Numerical modelling allowed us to calculate PAT deformation rates. Of particular interest are the slip rates, which reach maximum values up to 0.68-0.79 mm/yr for the last 0.8-1.2 Ma. In the same sector, extended toward the west to also include the Bologna area, I performed an analysis of the relationships between fluid migration and structural framework, which is believed to be directly responsible for the formation of the locally outcropping mud volcanoes. Subsurface structures were mapped using seismic sections, confirming a straightforward correspondence between anticlines and mud volcanoes, being generally located on the anticline surface projection (roughly close to the anticline axis or in the anticline’s forelimb/backlimb). In particular, I investigated the Dragone di Sassuno mud volcano, evidencing this correlation and the link between the upward migrating fluids and the fracture array associated to the associated anticline: fluids exploit systematic fractures to reach the surface. An in depth analysis was carried out to investigate possible relationships linking the seismic activity of the PAT and the mud volcano eruptions. I investigated three seismic events with M > 4 occurred between 1779 and 1780 and I modelled the variation of static stresses associated to these earthquakes, that occurred almost contemporaneously to large eruption of the Dragone di Sassuno mud volcano, as described by the Author Serafino Calindri at the time of eruption. I observed that static stress changes have not been responsible for triggering the mud volcano eruption, being the calculated values negative. Negative values imply that the fractures forming the feeder dyke systems were maintained closed by the stresses, inhibiting the fluid flow. Therefore, I propose a correspondence between the dynamic stresses associated with these earthquakes and the mud volcano eruptions: dynamic stresses, shaking the reservoir, are able to increase overpressures and therefore favour the eruptions. Mud volcanoes have been investigated also along the Marche Apennine foothill, with particular focus on the mud volcano systems located in three key areas: S. Paolo di Jesi, Monteleone di Fermo and Offida. As for the Emilia-Romagna mud volcanoes, I characterized the structural setting of these systems using field and seismic data. Also in these cases, a correspondence between mud volcano and anticline structures has been observed. For the Monteloene di Fermo and Offida mud volcanoes eruptions, I investigated the possibility of seismic triggering. An interesting case study was offered by the activation of the Monteloene di Fermo mud volcanoes coincident with three major earthquakes of the 2016-2017 Central Italy seismic sequence. I verified that the dynamic stress associated with the seismic events and calculated at the mud volcano location was high enough (>4 bar) to be considered the main trigger of mud volcano activation. Normal stress changes were too small (less than 0.1 bar) to have had any influence on the triggering. Summarizing, I verified a “direct” relationship between structures (anticlines) and mud volcanoes, which steers their location, and an “indirect” one, linking the activation of seismic structures (also far from the mud volcanoes location, ~70-80 km) and the mud volcano eruptions. As a final investigation, I used two 3D seismic analogues from the Loyal Filed (Scotland) and from the Caerá Basin (Brazil) as seismic analogues for fluid migration mechanisms. The first dataset, highlighted fundamental similarities between the so called “fluid escape pipes” -here widely distributed- and mud volcanoes. Both are linked to anticline structures, and fluids responsible for their formation exploit systematic fractures to migrate upward. In addition, the fluid responsible for the Loyal Field fluid escape pipes formation, revealed to be composed of a mud-water-hydrocarbons mixture, as in the case of the Italian mud volcanoes. A further analogy is represented by the evolution and activity of these geological objects, alternating periods of large eruptions to periods of quiescence. The second dataset evidenced the presence of a very significant amount of depressions resembling pockmarks at the seabottom, generally indicating the presence of a strong fluid flux. Therefore, it has been analysed in order to obtain further information on a seismic analogue on fluid migration mechanisms and possible relation with structures. Nonetheless, this analogue case study from a different tectonic context highlighted the dependence of these fluid migration pathways on large-scale sediment waves, vertically stacking their permeable trough to form preferential paths for fluid “passive” migration. Therefore, in this case I did not highlight any correspondence between tectonic structures and fluid migration.