Seismic images and rock properties of the very shallow structure of Campi Flegrei caldera (southern Italy)

In September 2001, an extensive active-seismic investigation (Serapis experiment) was carried out in the Gulfs of Naples and Pozzuoli, with the aim of investigating and reconstructing the shallow crustal structure of the Campi Flegrei caldera, and possibly identifying its feeding system at depth. The present study provides a joint analysis of the very shallow seismic reflection data and tomographic images based on the Serapis dataset. This is achieved by reflection seismic sections obtained by the 3D data gathering and through refined P-velocity images of the shallowest layer of Pozzuoli Gulf (z < 1,000 m). From the refined Vp model, the overall picture of the velocity distribution confirms the presence of a complex arc-shaped anomaly that borders the bay offshore. The deeper part of the anomaly (beneath 700 m, with Vp > 3,500 m/s) correlates with units made up of agglomerate tuffs and interbedded lava, which form the southern edge of the caldera, which was probably formed following the two large ignimbritic eruptions that marked the evolutionary history of the area under study. The upper part of the anomaly that tends to split into two parallel arcs is correlated with dikes, volcanic mounds and hydrothermal alteration zones noted in previous shallow reflection seismic analyses. The depth of the transition between the upper and lower parts of the anomaly is characterized by an abrupt Vp increase on the one-dimensional (1D) profiles extracted from the 3D tomographic model and by the presence of a strong reflector located at about 0.6/0.7 s Two Way Time (TWT) on Common Mid Point gathers. The move-out velocity analysis and stack of the P–P and P–S reflections at the layer bottom allowed to estimate relatively high Vp/Vs values (3.7 ± 0.9). This hypothesis has been tested by a theoretical rock physical modeling of the Vp/Vs ratio as a function of porosity suggesting that the shallow layer is likely formed by incoherent, water saturated, volcanic and marine sediments that filled Pozzuoli Bay during the post-caldera activity.     

Magma chamber evolution prior to the Campanian Ignimbrite and Neapolitan Yellow Tuff eruptions (Campi Flegrei,Italy)

The Campi Flegrei (Campanian Region, Italy) experienced two cataclysmic caldera-forming eruptions which produced the Campanian Ignimbrite (39 ka, CI) and the Neapolitan Yellow Tuff (15 ka, NYT). We studied the minor eruptions before both these large events to understand magma chamber evolution leading towards such catastrophic eruptions. Major, trace element, and Sr and Nd isotope compositions of pre-Campanian Ignimbrite and pre-Neapolitan Yellow Tuff products define distinct geochemical groups, which are here interpreted as distinct magma batches. These batches do not show any transitional trend towards the CI and NYT eruptions. The CI and NYT systems are decoupled geochemically and isotopically. At least one of the pre-CI and one of the pre-NYT erupted magma batches qualifies as mixing endmembers for the large CI and NYT eruptions, and thus, must have been stored in reservoirs for some time to remain available for the CI and NYT eruptions. The least evolved, isotopically distinct magma compositions that are typical of the last phases of the NYT and CI eruptions did not occur before caldera-forming events. Based on the new data, we propose the following scenario: Multiple magma chambers with distinct compositions existed below the Campi Flegrei before the CI and NYT eruptions and remained generally separated for some time unless new magma was recharged. In each case, one of the residing magma reservoirs was recharged by a new large-volume magma input of intermediate composition from a deeper differentiating magma reservoir. This may have triggered the coalescence of the previously separated reservoirs into one large chamber which fed the cataclysmic caldera-forming eruption. Large magma chambers in the Campi Flegrei may therefore be ephemeral features, interrupted by periods of evolution in individual, separated magma reservoirs.     

Structure and evolution of the Bay of Pozzuoli (Italy) using marine seismic reflection data: implications for collapse of the Campi Flegrei caldera

Digital marine seismic reflection data acquired in 1973 in the Bay of Pozzuoli, and recently reprocessed, were used to study the volcanological evolution of the marine sector of Campi Flegrei Caldera during the last 37 ka. In order to gain more information, interpretation also involved estimation of the "pseudo-velocity" and the "pseudo-density" from the resistivity logs of two onshore deep exploration wells. The main results are: (1) discovery of ancient pre-18 ka and post-37 ka submarine and mainly effusive volcanic activity, along coeval emission centers located at the edges of Campi Flegrei Caldera; (2) confirmation that the caldera collapse in the marine sector of Campi Flegrei seems strongly controlled by regional NE–SW and NW–SE structural discontinuities; (3) the finding of at least two episodes of collapse in the bay; and (4) identification of a post-18 ka volcanic deflation phase that has caused about 150–200 m of subsidence in the central sector of the Bay of Pozzuoli in the last 18 ka.Editorial responsibilty: T. Druitt     

Shallow Seismic Velocity Structure of the Betic Cordillera (Southern Spain) from Modelling of Rayleigh Wave Dispersion

A detailed dispersion analysis of Rayleigh waves generated by local earthquakes and occasionally by blasts that occurred in southern Spain, was undertaken to obtain the shear-wave velocity structure of the region at shallow depth. Our database includes seismograms generated by 35 seismic events that were recorded by 15 single-component short-period stations from 1990 to 1995. All these events have focal depths less than 10 km and body-wave magnitudes between 3.0 and 4.0, and they were all recorded at distances between 40 and 300 km from the epicentre. We analysed a total of 90 source-station Rayleigh-wave paths. The collected data were processed by standard digital filtering techniques to obtain Rayleigh-wave group-velocity dispersion measurements. The path-averaged group velocities vary from 1.12 to 2.25 km/s within the 1.0-6.0 s period interval. Then, using a stochastic inversion approach we obtained 1-D shear-wave velocity–depth models across the study area, which were resolved to a depth of circa 5 km. The inverted shear-wave velocities range approximately between 1.0 and 3.8 km/s with a standard deviation range of 0.05–0.16 km/s, and show significant variations from region to region. These results were combined to produce 3-D images via volumetric modelling and data visualization. We present images that show different shear velocity patterns for the Betic Cordillera. Looking at the velocity distribution at various depths and at vertical sections, we discuss of the study area in terms of subsurface structure and S-wave velocity distribution (low velocity channels, basement depth, etc.) at very shallow depths (0–5 km). Our results characterize the region sufficiently and lead to a correlation of shear-wave velocity with the different geological units features.     

Vent area and depositional mechanisms of the Upper Member of the Neapolitan Yellow Tuff (Campi Flegrei, Italy): new insights from directional fabric through image analysis

In order to provide new information about the source area and depositional mechanisms of the Upper Member of the Neapolitan Yellow Tuff (NYT), a prominent pyroclastic deposit of the Campi Flegrei Volcanic District (southern Italy), statistics on directional fabric, by means of computer-assisted image analysis on 32 rock samples, were compiled. Seventeen samples were collected along vertical direction on two selected exposures and fifteen were taken from outcrops widely distributed all around the Campi Flegrei Volcanic District. Fabric measurements within the investigated successions reveal a vertically homogeneous direction of the mean particle iso-orientation, with considerable variability in the strength of particle iso-orientation even at cm-scale. The existence of particle iso-orientation can be related to continuous sedimentation from a concentrated bedload region beneath suspension currents, producing massive or inversely graded beds by traction carpet sedimentation. The considerable vertical variability in the strength of iso-orientation is the result of very unstable flow regimes, up to the extreme condition of discrete depositional events, with a variable combination of traction carpet and/or direct suspension sedimentation. The vertical homogeneity in the mean orientation values, found in the investigated sections, may derive from the sequential deposition of laminae to thin beds, whose relatively flat upper surfaces were unable to significantly deflect the depositional system of the following currents. According to the observed homogeneous mean particle orientation values along the investigated vertical profiles, samples collected through areal distribution are considered representative of the local paleo-flow directions of the whole deposit. The mean directions of the samples collected areally show two different coherent patterns which point to the existence of two different source areas. The first, which includes all samples from the northern outcrops, appears to converge in a narrow area about 2 km NE of the town of Pozzuoli, largely in coincidence with the inferred area on the basis of the pumice fall distribution. The second, which includes samples from Capo Miseno and Posillipo areas, points to the central part of the Pozzuoli Bay, about 4 km offshore the town of Pozzuoli.     

Chemical and isotopic signatures of the basement rocks from the Campi Flegrei geothermal field (Naples, southern Italy): inferences about the origin and evolution of its hydrothermal fluids

The Campi Flegrei (Naples, Campanian Plain, southern Italy) geothermal system is hosted by Quaternary volcanic rocks erupted before, during and after the formation of the caldera that represents one of the major structural features in the Neapolitan area. The volcanic products rest on a Mesozoic carbonate basement, cropping out north, east and south of the area. Chemical (major, minor and trace elements) and stable isotope (C, H, O) analyses were conducted on drill-core samples recovered from geothermal wells MF-1, MF-5, SV-1 and SV-3, at depths of ˜ 1100 to 2900 m. The study was complemented by petrographic and SEM examination of thin sections. The water which feeds the system is both marine and meteoric in origin. Mineral zonation typical of a high-temperature geothermal system exists in all the geothermal wells; measured temperatures in wells are as high as ˜ 400 °C. The chemical composition of the waters suggests the existence of two reservoirs: a shallow reservoir (depth < 2000 m) fed by seawater that boiled at 320 °C and became progressively diluted by steam-heated local meteoric water during its ascent; and a deeper reservoir (depth > 2000 m) of hypersaline water. The drill-cores are mainly hydrothermally altered volcanics of trachy-latitic affinity, but some altered pelites and limestones are also present. Published Na, Mg and K concentrations of selected geothermal waters indicate that the hydrothermal fluids are in equilibrium with their host rocks, with respect to K-feldspar, albite, sericite and chlorite. The measured δ¹⁸O_(SMOW) values of rocks range from +4.3 to + 16.5%. The measured δD_(SMOW) values range from − 79 to − 46%. The calculated isotopic composition of the fluids at equilibrium with the samples vary from + 1 to + 8.3%. δ¹⁸O and from − 52 to + 1%. δD. The estimated isotopic composition of the waters at equilibrium with the studied samples confirmed the existence of two distinct fluid types circulating in the geothermal system. The shallower has a marine water signature, while the deeper water has a signature consistent both with magmatic and meteoric origins. In the latter case, the recharge of this aquifer likely occurs at the outcrop of the Mesozoic Limestones surrounding the Campanian Plain; after infiltration, the water percolates through evaporitic layers, becoming hypersaline and D-depleted.     

Simultaneous inversion of deformation and gravity changes in a horizontally layered half-space: Evidences for magma intrusion during the 1982–1984 unrest at Campi Flegrei caldera (Italy)

A very large uplift (about 1.8 m) occurred in the period 1982–1984 at Campi Flegrei caldera, Italy, without culminating in an eruption. A still-standing controversy accompanies the interpretation of deformation and gravity changes recorded during the unrest, which were interpreted to result from the sub-surface magmatic reservoir by some authors and from the hydrothermal system or hybrid sources by others. Here for the first time we take into account crustal layering while inverting leveling, EDM, and gravity data using uniformly-pressurized sources, namely small vertical spheroids and finite horizontal penny-shaped sources. The weight of EDM data in the misfit function is chosen from a trade-off curve in order to balance the compromise between fitting the leveling and the EDM data well. Models using a homogeneous medium cannot give a good simultaneous fit to leveling and EDM deformation data of the 1982–1984 unrest, whereas incorporating a layered structure (determined from seismically derived estimates of the P wave speed for the crust, and not adjusted to improve the fit in any of the inversions) allows a significantly better fit. Also, layering affects the sub-surface mass redistribution effects on gravity changes, and we show that the retrieved intrusion density is in full agreement with densities for highly evolved magmas expected at the Campi Flegrei caldera for depths of 3 to 4 km, ruling out hydrothermal fluids as the primary cause of the 1982–1984 unrest. The source of the 1982–1984 CF unrest was probably a shallow (about 3-km deep) penny-shaped magma intrusion fed by a deeper magma chamber; source overpressure was few MPa.     

Monte Carlo Inversion of DInSAR Data for Dislocation Modeling: Application to the 1997 Umbria-Marche Seismic Sequence (Central Italy)

— The study of surface deformation due to seismic activity is often made using dislocations with uniform slip and simple geometries. A better modeling of coseismic and postseismic surface displacements can be obtained by using dislocations with variable slip and nonregular shapes. This is consistent with the asperity model of fault surfaces, assuming a friction distribution on faults made of locked zones with much higher friction than surrounding zones. In this paper we consider the 1997–1998 Colfiorito seismic sequence. The coseismic surface displacements in the Colfiorito zone are used in order to infer the slip distribution on the fault surface at different stages of the sequence. The displacement field has been modeled varying the slip distribution on the fault, and comparing the deformation observed by SAR and GPS techniques with model results. The slip distribution is calculated by Monte Carlo simulations on a normal fault with the dip angle equal to 40°. A good approximation is obtained by using square asperity units of 1.5×1.5 km². In the first stage, we employed a simplified model with uniform slip, in which each asperity unit is allowed to slip a constant amount or not to slip at all, and in the second stage, we evaluate the slip distribution in the dislocation area determined by the Monte Carlo inversion: in this case we allow unit cells to undergo different values of slip in order to refine the initial dislocation model. The results show that the 1997 seismic events of the sequence can be modeled by irregular dislocations, obtaining a good fit to the DInSAR and GPS observations. The model also confirms the results of previous studies by a different methodology, defining the distribution of asperities on the fault plane using the fault geometry, the geodetic data and the seismic moment of the 1997–1998 Colfiorito seismic sequence. Furthermore, the analysis of 1997 aftershocks in the seismogenic region shows a strong correlation between most events and the asperity distribution, which can be considered as an independent test of the validity of the model.