Experimental study of the seismic properties of the Eastern Alps (Italy) along the Aurina–Tures–Badia Valleys transect

In order to better constrain the interpretation and the nature of the seismic reflectors, experimental measurements at high confining pressure (up to 300 MPa) and room temperature of the compressional wave velocity (Vp) on 10 samples representative of the most common lithologies along the Aurina (Ahrntal), Tures (Tauferer Tal), and Badia (Abtei Tal) Valleys profile (Eastern Alps, Italy) have been performed. For each sample, the speed of ultrasonic waves was measured in three mutually perpendicular directions, parallel and normal to the rock foliation and lineation.The main results are:(a) Good agreement between the calculated vs. measured modal compositions of the considered rocks, indicating that they were presumably equilibrated at the estimated P–T conditions; therefore, the seismic properties are representative of the crustal level indicated by the thermobarometry.(b) Measured and calculated average Vp are in good agreement, and are typical of mid-crustal level (6.0–6.5 km/s). Only the amphibolites show Vp typical of the lower crust (7.2 km/s).(c) The seismic anisotropy of metapelites is very high (12–27%), both with orthorhombic and transverse isotropy symmetry; amphibolites are transversely isotropic with an anisotropy of 8%; orthogneisses and granitoids are isotropic or weakly anisotropic.(d) The contacts between amphibolites and all other rock types may generate good reflections, provided they are not steeply inclined. Although the metamorphic foliation remains steeply inclined, discordant buried sub-horizontal igneous contacts may be detected.     

Note on the seismicity of northeastern Italy (Friuli area)

The kinematics and mechanics of the seismicity of the Friuli area in northeastern Italy are discussed, especially the focal mechanism of the catastrophic earthquake of May 6, 1976. Alternative solutions are analyzed. A N-S section, down to 60 km below sea-level, across the Alps from the Bavarian foreland to the Adriatic Sea illustrates these possibilities.     

Physical properties of ultrahigh-pressure metamorphic rocks from the Sulu terrain, eastern central China: implications for the seismic structure at the Donghai (CCSD) drilling site

About 30 samples representing major lithologies of Sulu ultrahigh-pressure (UHP) metamorphic rocks were collected from surface exposures and exploration wells, and compressional (Vp) and shear wave (Vs) velocities and their directional dependence (anisotropy) were determined over a range of constant confining pressures up to 600 MPa and temperatures ranging from 20 to 600 °C. Samples range in composition from acidic to ultramafic. P- and S-wave velocities measured at 600 MPa vary from 5.08 to 8.64 km/s and 2.34 to 4.93 km/s, respectively. Densities are in the range from 2.60 to 3.68 g/cm³. To make a direct tie between seismic measurements (refraction and reflection) and subsurface lithologies, the experimental velocity data (corresponding to shallow depths) were used to calculate velocity profiles for the different lithologies and profiles of reflection coefficients at possible lithologic interfaces across the projected 5000-m Chinese Continental Scientific Drilling Program (CCSD) crustal segment. Comparison of calculated in situ velocities with respective intrinsic velocities suggests that the in situ velocities at shallow depths are lowered by an increased abundance of open microcracks. The strongly reflective zone beneath the Donghai drill site can be explained by the impedance contrasts between the different lithologies. Contacts between eclogite/peridotite and felsic rocks (gt-gneiss, granitic gneiss), in particular, may give rise to strong seismic reflections. In addition, shear-induced (lattice preferred orientation (LPO)-related) seismic anisotropy can increase reflectivity. For the explanation of the high velocity bodies (>6.4 km/s) around 1000 m and below 3200-m depth, large proportions of eclogite/peridotite (about 40 and 30 vol.%, respectively) are needed.     

Anisotropic seismic properties of the upper mantle beneath the Torre Alfina area (Northern Apennines, Central Italy)

Central Italy is an active tectonic area that has been recently studied by several regional mantle, Pn and SKS, studies which revealed the presence of a strong regional anisotropy. In this paper, we present the first petrophysical results on the only mantle xenoliths from Central Italy, which place new constraints on the upper mantle structures of this region. The Torre Alfina mantle xenoliths are very small in size, from few millimetres to about 1.5 cm. They are mainly dunites and harzburgites, with subordinate lherzolites and wehrlites. Since olivine and spinel are always present, they should have crystallised in the spinel-bearing lherzolite field. Their mineralogical composition is ol+spl±opx±cpx. Both olivines and pyroxenes are present as porphyroclasts and as neoblasts. The xenoliths show different degrees of recrystallization. Geothermobarometry on these xenoliths give a temperature range of 1040±40 °C and a pressure estimate of about 1.5 GPa, corresponding to 50 to 60 km depth. Previous seismic studies have estimated the Moho to be at 20 to 25 km in this region, hence the xenoliths come from a hot mantle, probably asthenospheric, below a lithosphere of about 25 to 40 km in thickness below the Moho. We measure the crystallographic preferred orientation (CPO) of olivines and pyroxenes using a SEM and the Electron Back Scattered Diffraction (EBSD) technique. The CPO shows all three axes of olivine are tightly clustered: [100] axis is typically more tightly clustered than [010] and [001] is the most widely distributed axis. The fabric strength expressed by the integral J index, varies from 4.5 to 25.9, and decreases with the degree of recrystallization. We use CPO data to calculate anisotropic seismic properties of the xenoliths. They are very homogenous and probably statistically representative of the mantle below the Torre Alfina area. Vp ranges from 8.4 to 9.1 km/s, Vs₁ from 4.8 to 5.0 km/s. The seismic anisotropy is more variable; AVp ranges from 9.8% to 19.3% and AVs from 7.3% to 13.4%. The majority of the xenoliths display an orthorhombic seismic symmetry, but xenoliths with a transverse isotropic behaviour have also been observed.

We consider four geodynamic models for the source region of the xenoliths (extension, shear, upwelling, slab tilted), defined by different orientations of the structural reference frame, and we calculated for each model the variation of the seismic properties with temperature, pressure and volume fraction of orthopyroxene. After comparing this variation of calculated seismic parameters with seismic observations from the region, we form the hypothesis that the xenoliths come from either an extensional tectonic zone (lineation X and foliation plane XY horizontal) or transcurrent shear zone (lineation X horizontal and foliation plane XY vertical) and that the mantle beneath Torre Alfina is composed by 70% olivine and 30% orthopyroxene forming an anisotropic layer of about 160 or 110 km in thickness, respectively.     

Seismic reflectivity of detachment faults of the Iberian and Tethyan distal continental margins based on geological and petrophysical data

Low-angle detachment faults are key to our understanding of the tectonic evolution of magma-poor rifted continental margins. In seismic images of present-day rifted margins the identification and interpretation of such features is, however, notoriously difficult and ambiguous. We address this problem by studying the structure and seismic response of such faults through a synoptic interpretation of petrophysical data and geological evidence from the distal segments of the present-day West Iberian and the ancient Tethyan margins. On the basis of the geologically well-constrained remnants of the Tethyan margins, which are spectacularly preserved and exposed in the Alps of Eastern Switzerland, vertical profiles at four key geological settings of a typical magma-poor rifted margin are constructed and their synthetic seismic responses are compared to the observed seismic data from corresponding locations in the present-day Iberian margin. The seismic structure of these profiles is considered as the sum of deterministic large-scale and the stochastic small-scale components. Both components are analyzed for all pertinent lithologies. The large-scale structures are derived from laboratory measurements on samples from both, the West Iberian and Tethyan margins, whereas the small-scale fluctuations are constrained predominantly on the basis of well-log data from the Iberian margin. Different realizations of the simulated stochastic small-scale velocity fluctuations illustrate the potential variability of impedance contrasts and its impact on the seismic response from lithological interfaces and fault structures. Our results indicate that the nature of the seismic response from low-angle detachment faults is largely determined through the fracture-healing behavior of the surrounding rocks. Geological evidence from the exposed fragments of the Tethyan margins indicate that fracture-healing is generally well developed in crustal lithologies, but largely absent in mantle lithologies. It is for this reason that low-angle, intra-crustal detachment faults tend to be seismically undetectable. Conversely, crust–mantle detachments have a complex and variable seismic response, depending on the nature of the damaged zone and on the frequency content of the seismic data. These model-based inferences are consistent with the available evidence from the present-day Iberian passive margin and thus open new perspectives for the interpretation of the corresponding seismic images.     

Seismic hazard in the Sannio-Matese area of Italy

The assessment of seismic hazard at five selected sites in the Sannio-Matese region is based on the computer program SRAMSC. Owing to the extensive historical data base for the output parameter, the MSK intensiy is chosen. The seismicity model is made up of five narrow area seismic sources. Circular or elliptical macroseismic fields are assigned to individual sources. A generalized Kövesligethy equation is used for this purpose as the attenuation relationship. The study reveals similar and a rather high hazard at the sites at Benevento, Boiano, and Melfi, which are located in the zone of highest seismic activity. At the Pomigliano and Lucera sites, the assessed hazard is much lower.     

Seismic properties of the upper mantle beneath Lanzarote (Canary Islands): Model predictions based on texture measurements by EBSD

We present a petrophysical analysis of upper mantle xenoliths, collected in the Quaternary alkali basalt fields (Series III and IV) from the island of Lanzarote. The samples consist of eight harzburgite and four dunite nodules, 5 to 15 cm in size, and exhibit a typical protogranular to porphyroclastic texture. An anomalous foliation resulting from strong recovery processes is observed in half of the specimens. The lattice preferred orientations (LPO) of olivine, orthopyroxene and clinopyroxene were measured using electron backscatter diffraction (EBSD). In most samples, olivine exhibits LPOs intermediate between the typical single crystal texture and the [100] fiber texture. Occasionally, the [010] fiber texture was also observed. Simultaneous occurrence of both types of fiber textures suggests the existence of more than one deformation regime, probably dominated by a simple shear component under low strain rate and moderate to high temperature. Orthopyroxene and clinopyroxene display a weaker but significant texture. The LPO data were used to calculate the seismic properties of the xenoliths at PT conditions obtained from geothermobarometry, and were compared to field geophysical data reported from the literature. The velocity of P-waves (7.9 km/s) obtained for a direction corresponding to the existing seismic transect is in good agreement with the most recent geophysical interpretation. Our results are consistent with a roughly W–E oriented fastest P-wave propagation direction in the uppermost mantle beneath the Canary Islands, and with the lithosphere structure proposed by previous authors involving a crust–mantle boundary at around 18 km in depth, overlaid by intermediate material between 11 and 18 km.     

Seismic velocity structure of a large mafic intrusion in the crust of central Denmark from project ESTRID

The origin of regional sedimentary basins is being investigated by the ESTRID project (Explosion Seismic Transects around a Rift In Denmark). This project investigates the mechanisms of the formation of wide, regional basins and their interrelation to previous rifting processes in the Danish–Norwegian Basin in the North Sea region. In May 2004 a 143 km long refraction seismic profile was acquired along the strike direction of a suspected major mafic intrusion in the crust in central Denmark. The data confirms the presence of a body with high seismic velocity (> 6.5 km/s) extending from a depth of  10–12 km depth into the lower crust. There is a remarkable Moho relief between 27 and 34 km depth along this new along-strike profile as based on ray-tracing modelling of PmP reflections. The lack of PmP reflections at a zone of very high velocity in the lowest crust (7.3–7.5 km/s) suggests a possible location of a feeder channel to the batholith. The presence of volcanic rocks of Carboniferous–Permian age above the intrusion (mafic batholith) suggests a similar age of the intrusion. An older obliquely crossing profile and two new fan profiles deployed perpendicular to the main ESTRID profile, show that the batholith is about 30–40 km wide. The existence of this large mafic batholith supports the hypothesis that the origin of the Danish–Norwegian Basin is related to cooling and contraction after intrusion of large amounts of mafic melts into the crust during the late Carboniferous and early Permian. The data and interpretations from project ESTRID will form the basis for subsidence modelling. Tentatively, we interpret the formation of the Danish–Norwegian Basin as a thermal subsidence basin, which developed after widespread rifting of the region.     

An integrated geophysical investigation of the upper crust in the epicentral area of the 1980, Ms=6.9, Irpinia earthquake (Southern Italy)

In this paper, we investigate the upper crustal structure of the Irpinia region, Southern Apennines thrust belt, Italy, through analysis and joint interpretation of gravity data, seismic reflection lines and subsurface information from many deep wells. The investigated region includes the epicentral area of the 1980 (M_s=6.9) Irpinia earthquake and is one of the Italian regions with the highest seismic hazard. The upper crustal structure is imaged by modeling a series of 15 SW-trending gravity profiles, spaced about 5 km apart, plentifully constrained by seismic reflection lines and wells, thus reducing the inherent ambiguity of the gravity modeling. Despite of the complexity of the modeled Bouguer anomalies, the application of a calibrating procedure to constrain the range of variability of the density values, as well as the use of geometric constraints, results in a good level of stability in the final density cross-sections, which in fact appear coherent both in the density values and in the geometrical features. The inferred model shows important lateral density variations that can be mostly related to NW-trending geologic structures. High-density bodies delineate carbonate platform thrust sheets and broad antiforms involving Mesozoic basinal rocks, while low-density shallow bodies are associated with Pliocene basins. In addition, important density (i.e. lithological) variations are evident along the strike of the range, the most relevant being an abrupt deepening of the Apulia Carbonate Platform in the southeastern part of the investigated region. In the epicentral region of the 1980 event, we find that the geometry of the high-density, high-velocity carbonates of the Apulia Platform appears correlated with the distribution of the aftershocks and with the P-wave velocity anomaly pattern as inferred from a previous local earthquake tomography. The structural highs of the Apulia Platform correspond to high-velocity regions, where aftershocks and coseismic slip of the mainshock are concentrated. This correlation suggests that the Apulia Carbonate Platform geometry played an important role in the rupture propagation and in the aftershock distribution.     

Seismological, geological and geophysical constraints for the Gualdo Tadino fault, Umbria–Marche Apennines (central Italy)

The April 3, 1998 Mw = 5.1 Gualdo Tadino earthquake (central Italy) was the last significant event in the 6-month-long Umbria–Marche seismic crisis. This event and its aftershocks occurred in an area where active faulting produces no striking geological and geomorphological effects. In this study, we investigated the ruptured fault using detailed seismological data and a re-processed and re-interpreted seismic reflection profile. Aftershock location and focal mechanisms were used to constrain the geometry and kinematics of the ruptured fault and a comparison was made with the subsurface image provided by the seismic profile. We found that the 1998 Gualdo Tadino earthquake occurred on a WSW-dipping, normal fault, with a length of about 8 km and a relatively gentle dip (30°–40°), confined between 3.5 and 7 km in depth. Kinematics of the mainshock and aftershocks revealed a NE-trending extension, in agreement with the regional stress field active in the Northern Apennines belt. The Mw = 5.1 earthquake originated above the top of the basement and ruptured within the sedimentary cover, which consists of an evaporites–carbonates multilayer. We hypothesised that the active fault does not reach the surface (blind normal fault).     

Miocene low-angle detachments and upper crust megaboudinage in the Mt. Amiata geothermal area (Northern Apennines,Italy)

Information from surface and subsurface geology (boreholes and seismic reflection lines) are used to depict the geometry of the extensional structures (low-angle normal faults and related Tuscan Nappe megaboudins) affecting the Mt. Amiata geothermal area and developed during the early stage of the extensional tectonics which affected the inner Northern Apennines and Tyrrhenian Sea from the Early-Middle Miocene. Normal faulting involved the thickened middle-upper crust after the collisional stage and, in the Mt. Amiata region, took place over relatively short periods (5-7 Ma) characterised by rapid extensional strain rates. Normal faults showing articulated geometry (flat-ramp-flat) characterised by subhorizontal detachments (flats) and synthetic ramps, caused widespread megaboudinage mainly in the sedimentary tectonic units and particularly in the Tuscan Nappe. Evaporites occurring at the base of the Tuscan Nappe, the deepest sedimentary tectonic unit of the Northern Apennines, controlled the geometry of the faults, and rift-raft tectonics may be the style of this first extensional phase. Three Tuscan Nappe extensional horses (megaboudins) have been recognised in the subsurface of the Mt. Amiata area. They are characterised, in map view, by elliptical shapes and show a mean NNW-SSE lengthening. They are delimited at the base and at the top by east-dipping flats, while their western and eastern margins coincide with east-dipping ramps. On the whole, considering their geometrical features, these megaboudins correspond to extensional horses belonging to an asymmetrical east-dipping extensional duplex system.

Rollover anticlines deformed the western ramp of the megaboudins and rotated the uppermost flat as well as all the structures previously developed, which became steeply-dipping to the west.     

Regional thermo-rheological field related to granite emplacement in the upper crust: implications for the Larderello area (Tuscany,Italy)

We modelled thermo-rheological perturbations, related to the emplacement of a magmatic body in the upper crust. This approach was considered relevant for the areas characterized by elevated surface heat flow and chiefly for the geothermal fields. The numerical conductive thermal model applied to the Larderello geothermal area in Tuscany, allowed to constrain size, depth and timing of emplacement of the pluton. We inferred that the emplacement of a magmatic body, at a minimum depth of 3 km, having a horizontal extension of 14 km and a maximum thickness of 8 km, can reasonably reproduce the observed regional surface heat flow anomaly of the Larderello area, when 300 (± 100) kyr are elapsed from the magma emplacement. Even assuming an incremental growth, the first magma injection should not be older than 1 ± 0.3 Ma.

Results of the thermal model were used to set up a rheological model and to simulate the drifting of the brittle-ductile transition during the cooling of the pluton. A comparison with the K-horizon profile, a prominent seismic reflector in the Larderello area, was then performed. It was found that the K-horizon approximately corresponds with the pluton roof and with the current location of the brittle-ductile transition.     

Three-dimensional earthquake locations and upper crustal structure of the Sannio-Matese region (southern Italy)

The Sannio-Matese region is one of the most seismically active regions of Italy and has been struck by large historical earthquakes. At present, the area is characterized by low magnitude background seismicity and small seismic sequences following M4 main events. In this paper, we show Vp and Vp/Vs models and 3D locations for a complete set of earthquakes occurring in the period 1991–2001. We observe a significant crustal heterogeneity, with large scale east-verging high Vp fault-related-folds, stacked by the Pliocene compression. The relocated earthquakes cluster along a 70° east-dipping, NW-striking plane located at the border of the high Vp thrust units. Normal fault earthquakes related to the young and active extension occur within these high Vp zones, interpreted as high strength material. We expect large future earthquakes to occur within these high Vp zones actually characterized by low magnitude seismicity at their borders.     

Elastic modeling of the Alto Tiberina normal fault (central Italy): geometry and lithological stratification influences on the local stress field

Two-dimensional, elastic, plane-strain, finite element models (FEMs) are generated to study the state of stress and failure induced by a low-angle normal fault, the Alto Tiberina Fault of the Northern Apennines (central Italy): it is beyond the scope of the present work to show that slip can occur on such a fault. The numerical study was performed to evaluate the influences on the local stress field of the litho-mechanical stratification of formations surrounding the fault, and those induced by the geometry of the fault. The performed models have shown the important role played by evaporites and basement formations of the Umbria-Marche succession as seismogenetic layers. The model results have also underlined that the flat-ramp geometry of the fault induces high relative concentration of stresses in correspondence with the low-angle, east-dipping, synthetic normal faults observed today in shallow depth near the Alto Tiberina surface trace. The stress regime predicted by the final model, in which the behavior of the Alto Tiberina together with the antithetic normal fault of Gubbio was simulated, reconciles available geological and geophysical observable to a greater extent. The numerical results can assist interpretation of the tectonic evolution of the region.     

Seismic anisotropy of the lithosphere around the Trans-European Suture Zone (TESZ) based on teleseismic body-wave data of the TOR experiment

A passive teleseismic experiment (TOR), traversing the northern part of the Trans-European Suture Zone (TESZ) in Germany, Denmark and Sweden, recorded data for tomography of the upper mantle with a lateral resolution of few tens of kilometers as well as for a detailed study of seismic anisotropy. A joint inversion of teleseismic P-residual spheres and shear-wave splitting parameters allows us to retrieve the 3D orientation of dipping anisotropic structures in different domains of the sub-crustal lithosphere. We distinguish three major domains of different large-scale fabric divided by first-order sutures cutting the whole lithosphere thickness. The Baltic Shield north of the Sorgenfrei–Tornquist Zone (STZ) is characterised by lithosphere thickness around 175 km and the anisotropy is modelled by olivine aggregate of hexagonal symmetry with the high-velocity (ac) foliation plane striking NW–SE and dipping to NE. Southward of the STZ, beneath the Norwegian–Danish Basin, the lithosphere thins abruptly to about 75 km. In this domain, between the STZ and the so-called Caledonian Deformation Front (CDF), the anisotropic structures strike NE–SW and the high-velocity (ac) foliation dips to NW. To the south of the CDF, beneath northern Germany, we observe a heterogeneous lithosphere with variable thickness and anisotropic structures with high velocity dipping predominantly to SW. Most of the anisotropy observed at TOR stations can be explained by a preferred olivine orientation frozen in the sub-crustal lithosphere. Beneath northern Germany, a part of the shear-wave splitting is probably caused by a present-day flow in the asthenosphere.     

Petrological constraints on seismic properties of the Slave upper mantle (Northern Canada)

Modes and compositions of minerals in Slave mantle xenoliths, together with their pressures and temperatures of equilibrium were used to derive model depth profiles of P- and S-wave velocities (Vp, Vs) for composites equivalent to peridotite, pyroxenite and eclogite. The rocks were modeled as isotropic aggregates with uniform distribution of crystal orientations, based on single-crystal elastic moduli and volume fractions of constituent minerals. Calculated seismic wave velocities are adjusted for in situ pressure and temperature conditions using (1) experimental P- and T- derivatives for bulk rocks' Vp and Vs, and (2) calculated P- and T- derivatives for bulk rocks' elastic moduli and densities. The peridotite seismic profiles match well with the globally averaged IASP91 model and with seismic tomography results for the Slave mantle. In peridotite, an observed increase of seismic wave velocities with depth is controlled by lower degrees of chemical depletion in the deeper upper mantle. In eclogite, seismic velocities increase more rapidly with depth than in peridotite. This follows from contrasting first-order pressure derivatives of bulk isotropic moduli for eclogite and peridotite, and from the lower compressibility of eclogite at high pressures. Our calculations suggest that depletion in cratonic mantle has a distinct seismic signature compared to non-cratonic mantle. Depleted mantle on cratons should have slower Vp, faster Vs and should show lower Poisson's ratios due to an orthopyroxene enrichment. For the modelled Slave craton xenoliths, the predicted effect on seismic wave velocities would be up to 0.05 km/s.     

Statistical occurrence analysis and spatio-temporal distribution of earthquakes in the Apennines (Italy)

We present two examples of statistical analysis of seismicity conducted by integrating geological, geophysical and seismological data with the aim to characterize the active stress field and to define the spatio-temporal distribution of large earthquakes. Moreover, our data will help to improve the knowledge of the “seismogenic behavior” of the areas and to provide useful information for seismic hazard evaluation.The earthquakes are described by two non-parametric statistical procedures integrating also tectonic-physical parameters to study the spatio-temporal variability.The results show that the areas are characterized by: 1) a stress regime with mainly extensional kinematics; 2) tectonic structures mainly oriented with the active stress field (S_hmin = N44° ± 18° in the southern Apennines and S_hmin = N50° ± 17° in the central Apennines); 3) cluster distribution of seismicity and 4) a high probability of earthquake occurrence (M > 5.5) in the next 10 years.     

Piperno from Campi Flegrei: a relevant stone in the historical and monumental heritage of Naples (Italy)

This paper reports novel data concerning the volcanic formation of the Piperno, a stone that has been used for a long time in the monumental architecture of Naples and the Campanian region (Southern Italy). A careful field survey of the outcrops or underground sites so far accessible allowed us to draw a lithostratigraphic reconstruction of the formation and enabled the recognition of the three main layers exploited in past centuries and which provided dimension stones for architectural and structural purposes. Petrophysical and mineralogical characterization also enabled us to discriminate the same layers previously identified. Different scoriae/matrix ratio, scoriae dimension and welding degree are the most peculiar features of these horizons. The combination of these three parameters, even though they did not affect the use of Piperno as dimension stone, somehow determined its variable susceptibility to the weathering agents. However, a linear relationship between the exposure time and the decay of the stone itself, at a parity of weathering intensity, was not observed. This assumption was proved by verifying that the specific layers of the volcanic sequence controlled the weathering degree of Piperno from monuments disregarding their age of emplacement.     

Impact of groundwater withdrawals on the interaction of multi-layered aquifers in the Viterbo geothermal area (central Italy)

The impact of groundwater withdrawals on the interaction between multi-layered aquifers with different water qualities in the Viterbo geothermal area (central Italy) was studied. In this area, deep thermal waters are used to supply thermal spas and public pools. A shallow overlying aquifer carries cold and fresh water, used for irrigation and the local drinking-water supply. Starting with a conceptual hydrogeological model, two simplified numerical models were implemented: a steady-state flow model of the entire groundwater system, and a steady-state flow and heat transport model of a representative area, which included complex interactions between the aquifers. The impact of increased withdrawals associated with potential future development of the thermal aquifer must be considered in terms of the water temperature of the existing thermal sources. However, withdrawals from the shallow aquifer might also influence the discharge of thermal sources and quality of the water withdrawn from the shallow wells. The exploitation of the two aquifers is dependent on the hydraulic conductivity and thickness of the intervening aquitard, which maintains the delicate hydrogeological equilibrium. Effective methods to control this equilibrium include monitoring the vertical gradient between the two aquifers and the residual discharge of natural thermal springs.     

Structural setting of the Bay of Naples (Italy) seismic reflection data: implications for Campanian volcanism

This paper focuses on the recent tectonic evolution of the Bay of Naples with the aim of exploring the connection between local tectonics and volcanism. We reprocessed the seismic reflection dataset acquired in the area in the late 1973. The new processing was highly successful in obtaining a decisive strong reduction of random noise, removal of coherent noise and reduction of spatial aliasing. Classical interpretative schemes and complex attributes of seismic traces were used to reconstruct fault kinematics and reflector patterns. The results show that the faults affecting the Bay of Naples exhibit prevailing NE structural strikes, with the exception of the Pozzuoli Caldera where NW patterns are also common. Many faults are subvertical and show seismic evidence of volcanic activity along them. A main alignment of conjugate NE–SW faults, named here as “Magnaghi–Sebeto line”, intersects several submarine volcanic banks and separates the bay into two sectors, characterized by important geological, geophysical and petrochemical differences. The structural configuration of the bay may reflect the occurrence of either oblique extension or a transfer zone of the NW–SE fault system, along which, in the Campanian–Lucanian Apennine chain, great vertical displacements occur.