Coastal Evolution , Hydrothermal Migration Pathways and Soft Deformation along the Campania Continental Shelf ( Southern Tyrrhenian Sea ) : Insights from High-Resolution Seismic Profiles

A closely spaced set of high-resolution Chirp-Sonar and Sparker profiles and swath bathymetric data was acquired in 2013 for the I-AMICA Project off the Volturno River mouth (Southern Tyrrhenian Sea) by the Istituto per l’Ambiente Marino Costiero (IAMC), National Research Council of Italy (CNR). The palaeo-topography of three key surfaces, represented by the bounding surfaces of the post-glacial lithosomes, was mapped by the interpolation of seismically detected reflectors. The morphology of the surface related to the Last Glacial Maximum (LGM) regression revealed the presence of fault linkages which defined a small-scale accommodation zone with an E–W trending interbasinal relative high. The observed set of oppositely dipping faults, NNWand ENE-directed, locally controlled the deposition of the paralic/deltaic bodies during the post-glacial rise in sea level, as testified by their wedge-shaped geometries and shifting depocentres. The deformation may be linked to the Campi Flegrei caldera collapse following the Neapolitan Yellow Tuff (NYT) eruption and aged 15 ka BP. The relevant thickness of the Transgressive System Tract (TST) testifies to an increased sediment yield and intense reworking in coastal areas, probably driven by the high volcanoclastic supply during volcanic paroxysm, almost coeval with the post-glacial transgression. Fluid escape features linked to an E–W striking fluid front at the outer shelf suggest the presence of an hydrothermal system controlled by the predominant direction of normal to oblique Quaternary-active faults and by lithologic discontinuities across the sedimentary pile.


Introduction
The stratigraphical architecture of continental shelf deposits reflects the morphodynamic activity of the offshore deltaic system, constituting an ideal laboratory to investigate the complex interplay among fluvial inputs, glacioeustatic fluctuations, and tectonics.Channel aggradation processes, local uplift, and inherited topographic features are often triggering factors for a delta lobe to avulse [1,2] and may be directly inferred from facies distribution, lateral extent, and thickness of the deltaic deposits.Furthermore, modern continental shelves offer the opportunity to observe and map the presence of gas in shallow water sediments [3][4][5].Biogenic methane is frequently produced in deltaic areas where clay-rich sequences are supported with a constant supply of organic matter.Nevertheless, the presence of thermal fluids could also have an effect in decreasing methane solubility, allowing the formation of gas bubbles and fluid-venting structures [5].
In this paper, we present a high-resolution Chirp-Sonar and Sparker seismic dataset along the Campania continental margin to interpret avulsion processes, fluid venting features, and anomalous thickness of the Upper Pleistocene-Holocene transgressive lithosomes in a more comprehensive regional structural framework.A very large dataset, derived from subsurface exploration of the Campania and Latium margin performed in the past decades [6][7][8][9], may serve as a stratigraphic reference for our purposes.
The strongly asymmetric Middle-Upper Pleistocene glacioeustatic fluctuations shaped the continental margins with slow and gradual sea level falling stages and by short-lived rise and highstand stages [10].This usually resulted in forced regression and lowstand deposits characterized by coarsening upward parasequences stacking at the shelf margin and thin transgressive bodies with a gradual fining-upward trends, finally reaching a climax in a maximum flooding area [11].
Local variations in sediment supply and changing rates in sea level rise resulted in facies distribution and thickness variability [11][12][13][14].Ground flow circulation patterns are related to permeability and deeply influenced by the stratigraphic context; accordingly, the variable grain size of the sedimentary pile at the shelf margin can coincide with sealing bodies or carrier beds.Therefore, the stacking pattern as well as the vertical and lateral variations of sediment textures distributions within the depositional sequences can deeply influence the configuration of migration pathways.Additionally, structural discontinuities can also act as preferred migration pathways, mechanical barriers, or mixed systems [15,16].Thus, such physical discontinuities may affect the diffusion of a geothermal system and control its geometric configuration [17,18].
On the basis of the geometric pattern of the structures composing the system, a recent classification of rising fluid types [19] resolved three different categories of features related to: (1) focused vertical fluid flow; (2) sediment remobilization processes; and (3) laterally extensive fluid flow.The most common seismic expressions of fluid escape features are pockmarks and pipe structures: columnar zones of vertically stacked, downward-bending reflections [20] which could be interpreted as collapse structures induced by the hydraulic fracture and fluid-driven erosion in conditions of high pore fluid pressure.These fluid venting structures, induced by focused vertical fluid flow, are very good timing indicators as they constrain the expulsion process to the reflection on which they occur [21].Laterally extensive fronts of acoustic turbidity may be related to fluid accumulation zones occurring in certain stratigraphic intervals or at certain depths.Fluid front morphologies have also been explained as related to pressure and bottom water temperature [3] or to sedimentary facies association and permeability [22].An accurate description of the fluid front morphologies is a fundamental tool to unravel the origin of the observed venting structures in order to better understand the late Quaternary evolution of the Campania margin.

Geologic Framework
The area under investigation is located along the Northern Campania shelf, off the Volturno coastal plain (Figure 1).The late Quaternary succession is characterized by prograding clinostratified deposits, truncated by the youngest Oxygen Isotopic Stage 2 (OIS 2)-related unconformity, recording multiple events from the subaerial exposure of the shelf to the successive sea level rise [6,23].Records of the latest sea level rise, which occurred from 18 ka BP to 5 ka, as well as the present highstand conditions describe the transgressive and highstand deposits constituting a sedimentary wedge decreasing in thickness towards the shelf edge [24][25][26].The proximity of the Volturno alluvial plain as well as the Campi Flegrei and Roccamonfina volcanic districts, respectively responsible for clastic and volcaniclastic sedimentary inputs, deeply controlled the Plio-Quaternary stratigraphic succession [27,28].The Middle-Upper Pleistocene glacioeustatic fluctuations were characterized by high frequency (20-100 ka) and high amplitude, ranging from +10 m to −120 m [29].Such variations occurred in the same depth range, leading to the erosion of the deposits related to a previous cycle during a sea level fall.The total subsidence enabled the preservation of deposits, acting as a key factor to drive the growth of the margin [6].Beginning in the Early Pleistocene, major NE-SW striking faults affected the Tyrrhenian sector of the Campania region with vertical slip rate values of 4.3 mm/year [30], individuating highly subsiding and asymmetrical half grabens [31][32][33][34][35]. Severe subsidence rates characterized this tectonic scenario in which NW-SE striking left-lateral transtensive faults were also active [35].Records of the Early-Middle Pleistocene horizontal throws have been reported on NW-SE left-lateral strike-slip faults in the Cilento region [34] and in the Sorrento area [35].A Middle Pleistocene (0.7-0.4 Ma) tectonic event has also been reported in the Naples Gulf [30,36,37], with the presence of major NE striking normal listric faults.Since the late Quaternary, the onset of a tectonic event with a NE-SW oriented σ3 [35] produced NW-SE normal faults with vertical offsets of a few tens of meters, as observed in the Campania [38] and in the Sele Plains [39].A post-0.4 Ma tectonic event was individuated in the Gaeta Gulf, reactivating previous NE-and NW-oriented structural planes and producing horizontal offsets on E-W striking faults [37,40].The same striking faults, as observed in the offshore of Naples Gulf, controlled in extension the deposition of the 39 ka Campania Grey Tuff [41,42] and influenced the 12-9 ka transgressive unit [30,43].Volcanic activity, ascribed to different districts and provinces (Roccamonfina, Campi Flegrei, Ischia Island, and Somma-Vesuvius), aligned on a NW-SE direction, has been controlled by the same orthogonal systems of deformations [44,45].Records of the oldest volcanic cycle have been observed in the offshore sectors of the western group of the Pontine islands with products emplaced from 4.2 Ma to 1 Ma ago [46].Volcanism followed an eastward migration with records of activity constrained at 0.8-0.13Ma ago on the island of Ventotene and from 0.63 Ma to 0.13 Ma to the Roccamonfina volcano [47].A southeastward migration of volcanic activity was also observed and revealed by the presence of several buried volcanic structures off the Southern Gaeta Gulf, Ischia, Campi Flegrei, and Naples areas, respectively emplaced from the Early to Upper Pleistocene [48].The Campi Flegrei volcanic district, ascribed to the latest tectonic and magmatic event, deeply influenced the evolution and geomorphological setting of the study area.Its structure is the result of two distinct caldera collapses originating, respectively, from major eruptions and pyroclastic currents: the Campanian Grey Tuff (CGT) at about 39 ka [42,49] and the Neapolitan Yellow Tuff (NYT), dating to ca. 15 ka [41,50]. 1

Methods and Dataset
The geophysical acquisition included a swath bathymetric and shallow high resolution (HR) single-channel seismic (Chirp-Sonar and Sparker) survey performed by the IAMC-CNR of Naples during the oceanic cruise I-AMICA 2013, on board the R/V "Urania" [52].Bathymetric data acquisition was carried out with the use of the Kongsberg Maritime AS (NOR) EM710 multibeam echosounder system.The system was keel mounted and interfaced with a Differential Global Positioning System.Positioning and attitude data were corrected during acquisition.A sound velocity probe was installed near the transducers providing the real-time sound speed for the beam steering, while a sound velocity profiler was used to record the property of the water column every 6-8 h in order to obtain the velocity profile required for accurate beam forming [53].The whole dataset covers a surface area of 270 km 2 in a depth range 16-500 m (Figure 2).The data were acquired using the interface Seafloor Information System (SIS) software and were processed using Hips and Sips 7.0.2 by Teledyne CARIS, Inc. (Fredericton, NB, Canada).Data were then used to produce a high-resolution digital terrain model (DTM) with a 5 m grid cell size (Figure 2).A high-resolution morpho-bathymetric map of the seafloor, covering the continental shelf off the Volturno river mouth and the uppermost sector of the slope, provided a fundamental tool to identify different classes of fluid escape features.A 1-kJ Sparker power supply with a multitip Sparker array by Geo-Resources BV (NLD) was used as a seismic source with a frequency of about 800 Hz, fired at a time interval of 1 s.Data were recorded using a single-channel streamer and Triton SB-L 7.2 software by Smartsea Technology Co., Ltd.(Wuxi, China) and processed by the GeoSuite Allworks software by Geo Marine Survey Systems (NLD), by applying an IIR filter, amplitude gain control, signature deconvolution, swell filter, time-variant gain, and muting of the water column; signal penetration exceeds 300 ms t.w.t. and the vertical resolution is ~2 m just beneath the seafloor.A densely spaced net of ultra-high-resolution seismic data was also recorded using the hull mounted mod.CAP-6600 Chirp-Sonar III Profiling System by Teledyne Benthos Inc. (Shingle Springs, CA, USA), operating with 16 transducers in the 2-7 kHz wide frequency band and with a 10-20 ms pulse length.A number of about 60 Chirp-Sonar profiles were logged at the same time as the swath-bathymetric data (Figure 1b).Given the densely spaced net of profiles, a 3D-like reconstruction of the seismically detected horizons ensued.
The uppermost 5 m of the 7 m-long core C5 recovered on the outer shelf (location in Figure 1b, 40 • 58 24.953" N, 13 • 47 02.514"E, 93 m below sea level, b.s.l.), represent the only tool for calibrating the seismics.Stratigraphic and magnetic susceptibility analysis in the first 5 m [54] enabled seismic calibration of the shallowest seismic horizons (Figure 3).
Data on the Late Quaternary outbuilding of the margin and the analogies of the geometries of deep unconformities, interpreted by previous authors along the Latium and Campania continental shelf [6,27,37,48], provided a guide for the interpretation of Sparker profiles (Table 1).

Seafloor Morphology
A variety of structures were displayed by the DTM obtained by the multibeam bathymetric survey (Figure 2).In the inner shelf, a pronounced break in slope at −25 m marked the transition from topset-to-foreset rollover depth of the subaqueous delta [55].The transition between the shelf, gently deepening seawards, and the continental slope occurred at −150 m depth, where the shelf break was located.The slope area exhibited an uneven topography with several submarine gullies and landslide scars [56].Most of the seabed features were located along a WNW-ESE orientation and were represented by pockmarks and elongated steps, with a direction strikingly parallel to the isobaths.However, two classes of seabed depressions were detected: sub-circular/elliptical and elongated ones.Most of the sub-circular craters had smaller dimensions with the diameter size ranging between 30 m to 70 m and were aligned in a NE-SW trending direction.The elongated craters, displaying a bigger diameter, exceeded 200-300 m along their longest axis, and were aligned in a NW-SE direction.Long alignment of such clusters corresponded to areas of coalescent pockmarks, which created a morphological step, evident both on the DTM and seismic profile (Figure 2b).

Seismic-Stratigraphic Surfaces and Units Derived from Chirp-Sonar Profiles
The densely spaced seismic lines collected during the bathymetric survey provided the high coverage necessary to realize 3D-like maps of the most relevant seismically detected horizons.Three main reflectors, bounding four major seismic units, have been mapped, which, top to bottom, are H1, H2, and H3 (Figure 4).A seismically transparent body constituting Unit 1 (U1), displayed an aggrading stacking pattern and downlap terminations.Its thickness decreases towards the shelf break in an external wedge-shaped configuration.This unit displays cylindrical reflection anomalies, constituted by a series of vertically stacked depressions affecting the unit up to the seabed (Figure 5).A clearly recognizable high-amplitude reflector (H1 horizon) separates U1 from the unit underneath (Unit 2, U2).Medium to low amplitude sub-parallel and aggrading reflections characterize this unit, displaying downlap terminations above a highly irregular reflector (H2 horizon).Occasionally patchy and wavy high-amplitude reflectors, which were enclosed in lenses and arranged in a backstepping external configuration, characterize Unit 3 (U3, Figure 4).
Unit 3 is wedge-shaped, with thickness increasing in the landward direction, reaching values up to 10 ms and revealing a sedimentation controlled, in some cases, by normal faulting (Figure 5).Bounding at the bottom Unit 3 (U3), reflector H3 has a highly irregular trend and is characterized by a strong acoustic impedance contrast with the region beneath it in which the signal is not always penetrating.The recognized pattern of faults, controlling the deposition of Unit 3 and cutting also the uppermost seismic units, is composed by NW-dipping normal faults, observed in the northern side of the studied area, forming a downthrown region in a tilted block configuration (Figure 5).In the southeastern side of the studied area, ENE-WSW trending, SE-dipping normal faults defined a relative morphological low.An interbasinal relative high, with a roughly trending E-W direction, defined by two outward-facing normal faults, separated these two topographic lows.The 3-D-like modelling of H1, H2, and H3 surfaces, evidenced two NE-SW-elongated depocentres of U3 unit, separated by a relative topographic high and a NE-SW-elongated depocentre of U2 unit, slightly shifted southwards (see Section 5.5).

The Stacking Pattern of the Seismic Units by the Sparker Profiles
Nine seismic units have been recognized on Sparker sections (Figure 6), of which the uppermost three correspond to those identified in the Chirp-Sonar profiles.Downlap terminations of the wedge-shaped unit (U1), characterized by low amplitude reflections, rest on the high-amplitude reflector (H1).Below it, low-amplitude subparallel reflections display downlap terminations above the undulating and irregular reflector H2.Wavy and high-amplitude reflectors lie below H2 and describe the lenticular bodies piled in a retrograding (backstepping) external configuration, with a thickness increasing in a landward direction (U3), exceeding 10 ms.A flat and gently seaward deepening reflector (H3) cuts, on the observed profiles, almost all the units underneath, producing toplap terminations and erosive truncations on the mostly below and laterally stacked U4, U5, U6, U7, and U8 seismic units.An upslope thinning unit (U4) was observed beyond the shelf break, exhibiting clear onlap terminations above the H4 reflector.Oblique parallel progradings characterize unit U5, passing laterally across a generally concordant H5 reflector to unit U6, displaying an external sigmoidal geometry.The migration of offlap breaks indicates, at the base, an aggradational pattern, then a marked progradational progression of the roll-over points with a decrease in acoustic transparency moving up the section.The H6 reflector marks the transition to the unit U7 with an external tabular configuration and a stratified aggradational to a progradational growth patterns, moving upwards.Unit U7 onlaps on the H7 reflector which represents a major erosional unconformity, gently dipping in a seaward direction.Below H7, unit U8 displays a prevalent aggrading migration of the offlap breaks and upper toplap terminations.This unit shows onlap terminations against the H8 reflector and a variable thickness among reflections, suggesting possible syndepositional ground movements.The unit U9, the lowermost detected unit, corresponds to a transparent acoustic body with a mound-shaped external form (Figure 6), extending in an E-W trending direction (see Section 5.4).
Laterally extensive pockets of amplitude anomalies have been observed on seismic profiles in the inner and outer shelf.
The seismic expression of the outer shelf pocket was an intermittent acoustic turbidity, alternated to columnar zones of amplitude masking, marked at the top by high amplitude disrupted reflectors, arranged in a saw-toothed top profile (Figures 6 and 7).The top boundary of the acoustic anomaly, ranging from 0.18 to 0.27 ms b.s.l., did not pertain to a specific stratigraphic horizon, but tracked several horizons of units U5, U6, and U7 (Figures 6 and 7).Accordingly, the acoustic turbidity front depicts intermittent and sharp lateral flanks and an irregular and gradual top front.In plan view, it is roughly E-W oriented over an approximate area of 40 km 2 .
In the inner portions of the shelf (Figure 6) a wide zone of acoustic signal wipe-out was also observed.Here, the acoustic blanking was complete, right below a stratigraphic reflector of U1 (Figure 6A).The shape of the front was sub-circular in plan view, almost centred at the Volturno River mouth (Figure 6C).
Distinct sites of vertically focused flows, marked by columnar zones of amplitude masking and disrupted reflectors, reached the seabed at the outer shelf and corresponded to craters (Figures 5a  and 7b).Locally, long alignments of seepages, likely of long existence, have been inhibiting the deposition of sediments of the highstand wedge, creating a non-deposition area.This non-deposition area corresponded to the alignment of elongated pockmarks (Figure 2b) and accounted for the creation of a minor morphological step at the seabed (Figure 7b).

Calibration of Seismic Sections: The Stratigraphic Framework
Literature data [54] reported stratigraphic and magnetic susceptibility analyses performed on the C5 core (Figure 3), recovered on the outer shelf, immediately offshore of the Volturno River mouth at 93 m b.s.l.Top to bottom, a dominant pelitic succession characterizes the first 5 m, interbedded with five major tephritic layers, identified by peaks in the magnetic susceptibility record, geochemical, and micropaleontological analysis performed on the >63 µm fractions [54].Coarser-grained deposits with thin sandy intercalations characterized the lowermost two metres at the base of the C5 core.
The hemipelagic succession recovered by the first 5 m of the core recorded a time span of 5 ky, proven by the correlation of tephritic layers to well-known eruptions from the Campi Flegrei and Somma Vesuvius districts.The C5 recorded part of the succession, interpreted as the present-day highstand muddy sequence, calibrated the topmost unit U1 made of low amplitude reflections.Below it, deposits related to the latest sea level rise, occurring from 18 ka BP to 5 ka BP [23][24][25][26]-or 18 ka to 7 ka according to [57]-correspond to distinctive marine and paralic facies associations of units U2 and U3, respectively.Patchy, high-amplitude, and retrograding acoustic bodies, constituting Unit 3 at the base of the transgressive intervals (Figure 5a), should include coarse-grained deposits [58].Such bodies, recovered on the mid-shelf, should have settled in a coastal environment, consistent with the position of the sea level at about −80 m, some 14 ka BP [59].Above it, a more distal facies association (U2) testified to the landward migration of the coast.This interpretation of units U2 and U3 matches with the well-acknowledged stratigraphic significance of the H3 reflector at the Mediterranean scale as the Last Glacial Maximum (LGM) erosional surface, recording multiple time-transgressive morpho-evolutive processes ranging from the progressive exposure of the shelf, due to the sea level drop, to the following post-glacial sea level rise ( [6] and references therein).
Recent studies [7,28] contributed to the definition of a very thick transgressive systems tract (TST), individuating and calibrating in the upper slope area off the Campi Flegrei area the seismic evidence of the NYT [41,50].This tephra has been found below the maximum flooding surface, confirming that the deposits bounded by this major tephritic layer and the LGM erosional surface can be ascribed to an early stage of the transgressive systems tract.
Lacking any possibility of direct ground truthing, the calibration of deeper horizons identified on the Sparker section can be attempted by analogies in geometrical characteristics of reflectors with other profiles acquired along the central Tyrrhenian margin (Table 1).
Volcanic buried bodies and four different unconformity-bounded units (PP, A, B, and C, respectively), constituting the Plio-Quaternary record, have been recognized off Gaeta Bay [37].
In particular, we refer to a major E-W trending structure recognized by the authors offshore of the Volturno River mouth, a volcanic edifice in place prior to 0.7 Ma [60], showing coastal onlaps of the above lying sequences.We have correlated the U9 unit with the well-known V0 Lower Pleistocene volcano [37].Following the interpretation of the authors, the U8 seismic unit should have corresponded to one of the sequences composing the third-order unconformity-bounded unit B (0.7-0.4 Ma), revealing coastal onlap on the described volcanic edifice at about 0.3 twt.This major unit was constituted by sigmoidal prograding sequences with significant aggradation processes.This was consistent with U8 and the lowermost part of the U7 seismic units, which showed the same geometric pattern.Starting from the uppermost part of U7 seismic unit, an oblique parallel clinoform geometry and the largely lateral stacking of sequences testified to a stricter progradational pattern.We correlated this succession to the unconformity-bounded unit C, i.e., the last major sequence composing the margin [37], recording a post-0.4Ma rapid landward shift of the offlap breaks.Thus, Sparker profiles recorded a Middle-Upper Pleistocene to Holocene succession.

Structural Lineaments
Normal faults cutting the transgressive and highstand units with very small offsets (ranging between 2 and 5 m) were displayed by Chirp-Sonar profiles.The 3D-like modelling of the LGM erosional surface (Figure 5b) outlined a relative high, bounded by two sets of outward facing faults.
Extended domains are usually segmented in fault systems whose terminations and geometrical transitions have been extensively studied [61][62][63][64][65].In this typically asymmetrical segmentation, fault linkages are generally ascribed to two main families of structures responsible for displacement transfer, known as accommodation and transfer zones [61].Accommodation zones correspond to areas in which soft and distributed deformation occur at the overlapping fault terminations [63,66,67].Hard-linkage processes, in response to an increased displacement, should be invoked for the formation of transfer zones, defined as areas of discrete deformations, with a very strong strike-slip component [61,64,[68][69][70].The observed deformations revealed a structural setting constituted by two opposite-facing sets of faults, separated by a relative topographic high located in the overstepping region (Figure 8).No strike-slip component has been seismically detected and the accommodation zone, individuated in the topographic relative high, had a slightly oblique orientation with respect to the NE-SW extension trending direction, revealed by the main NNW-SSE trending set of faults and reported in literature [33,35].
A general classification for accommodation zones is known [67] in which the antithetical ones have been further subdivided both on the basis of inward or outward facing faults and on the orientation of the transfer zone with respect to the regional σ3.Following this classification the observed structure, although small-scaled and with minimal offsets, seems to correspond in the arrangement to an antithetic oblique synclinal accommodation zone.The activity of these faults appears to be constrained by the locally observed growth patterns in the interpreted paralic transgressive bodies (Figure 5a), whose deposition is referred to the initial stages of the last sea level rise which occurred from 18 ka to 5 ka BP [23][24][25][26].Around 15 ka, the Campi Flegrei caldera was collapsing [50] in the onland area (Figure 1a), producing an important volcaniclastic aggradation and great volumes of pyroclastic deposits.Such a dramatic event should have influenced and perturbed the stress field distribution, accommodated in the study area by the reconstructed structure.

Degassing Features
The high-resolution bathymetry depicts the presence of two classes of pockmarks: elliptical and elongated, showing a NE and a NW trending direction, respectively (Figure 2).The NE-SW-oriented small-shaped class of pockmarks seems to match the direction of gullies in the slope area and the maximum dip incline of the shelf (Figure 2c,d) quite well.It is acknowledged that buried palaeo-channels and gullies, infilled with coarse-grained lithology which are more permeable than the muddy bulk of the continental margin, may constitute a preferred migration pathway for fluids [71,72].Although no evidence of buried paleo-channels has emerged from the HR profiles parallel to the coasts, possibly due to a combination of low seismic signal penetration and noisy records, the coinciding orientation between gullies on the upper slope and the NE-SW-oriented pockmark on the shelf, lets us infer the occurrence of an undetected buried paleo-drainage system which may have favoured the vertical migration along such preferred paths.
The mapping of craters and depressions at the outer shelf suggests that the majority of them occur above the laterally extensive zones of acoustic anomalies or in their immediate vicinities.Vertical funnel-shaped and downward-tapering anomalies, interpreted as pipes, located immediately below the main detected pockmark field, are rooted in a deeper fluid pocket and have been individuated in the southeastern area of the interpreted grid.The morphology of the fluids' front depicts sharp and distinct lateral flanks, giving this pocket a WNW-ESE trending direction (Figure 9).Thus, the second class of pockmarks seems to match the overall direction of the fluid front and coincides with distinct sites of vertically focused flows (Figure 7) quite well.Observations of the vertical sections show that the longest alignment of coalescent pockmarks (Figure 2b) among those NW-oriented locally corresponds to an abrupt closure of the wedge-shaped HST and to a 6-km long morphological step at the seabed (the red box in Figure 7b), localized at about 140 m b.s.l.Literature data [54] report a depositional rate of about 0.1 cm/year over the last 4.5 ka at 93 m b.s.l. and the prevalence of fine-grained particles pertaining to an hemipelagic deposition.Therefore, a lower depositional rate can be inferred at the outer shelf where the morphological step occurs.Under these conditions, even a sporadic activity of fluid seeps could be effective in controlling the dispersal of fines and modify the shape of the HST pinch-out over time.Such coinciding features may prove an active venting process up to recent times, inhibiting the deposition of post-glacial fine sediments locally at the shelf edge.
The bulk of the laterally extensive fluid front has been individuated in the Upper Pleistocene prograding units.More coarse-grained sediments, in comparison with the Holocene highstand units, constitute the prograding sequences composing the margin.Such sediments could host fluidsproviding good conditions for their lateral migration [22].The main trending direction of this front is similar to the one observed for the detected roughly E-W striking volcanic body, despite its location seems to have shifted southward (Figure 10).The acoustic masking zone's peculiar facies association, strikingly different in the acoustic response when compared to inner shelf observations yet similar in the trending direction of the volcanic body, suggests that these venting features could be related to buried magmatic bodies.

Biogenic Gas or Hydrothermal Venting?
Amplitude anomalies reveal the boundaries of the detected turbidity pocket along the geologic pile.This suggests abrupt or more gradual variations in acoustic waves velocity, in function of the spatial distribution and concentrations of fluids [73], which is dependent on temperature, pressure, and on the medium permeability [5].Accumulation of fluids may also have resulted from cross-stratal migration processes, which can occur through faults, intrusive bodies, permeable beds and pipe-like structures [74].Presently, a true and unambiguous identification of the type, source, and concentration of fluids based on seismic features is not straightforward.
In this study, two distinct fronts with reflection amplitude anomalies and acoustic blanking are compared.The first was located in the inner shelf and consisted in a total acoustic blanking below a stratigraphic reflector of U1 in the first sub-bottom of the prodelta area, down to 50 m b.s.l.The acoustic blanking was characterized by sub-vertical sharp lateral boundaries, delimiting a sub-circular region centred off the Volturno River mouth.A comparison of this acoustic facies with others HR profiles suggests that such a process could have been induced by the occurrence of subsurface biogenic methane, pervasive to the sedimentary environment (submarine delta lobes) and produced in situ by an almost syndepositional degradation of organic compound [21].
The second turbidity front was located at the shelf-to-slope transition area, where lateral distinct boundaries delimited a region in which a diffusive acoustic attenuation of reflections was observed, marked at the top by enhanced reflections arranged in a saw-toothed profile.Such strong impedance contrast could have been driven by the lithology of the depositional sequences building blocks [6,11,24] and thus was related to permeability: fluids migrating into coarser grained sediments seem to produce bubbles which scatter the sound and produce high amplitude reflectors [21].The saw-toothed profile was revealed by enhanced reflection outlines of the top boundary of several sub-units within U7 up to U5, which were characterized by a pronounced progradational stacking patterns.The lateral and vertical stacking of coarsening-upward sub-units of the outer margin succession [6] accounted for abrupt variations in permeability along the sedimentary pile, driving the formation of fluids-saturated horizons.The mapping of this front revealed a WNW-ESE trending direction, likewise the volcanic structure (U9).A pockmark field was observed above it at the seabed, directly linked to the deep acoustic turbidity pocket (Figure 7a).Hydrothermal fluids usually come from devolatilization processes and have been frequently observed in other volcanic continental margins [75].
Thus, below the laterally extensive gas pocket, a magmatic body should have triggered the venting structures.This suggests that such bodies may have followed an E-W trend, potentially correlating to the main direction of normal to oblique faults active during the Quaternary period [35,40,[76][77][78].A recent study [79] demonstrated the presence of an active hydrothermal system in the onshore sector, corroborating our hypothesis of hydrothermal pumping.

Coastal Evolution
The facies distribution and the spatial organization of the shallow marine depocentres reflect the morpho-dynamicity of coastal sectors along clastic coasts.The elaboration of two isopach maps representing, respectively, the paralic/deltaic and marine lithosomes (U3 and U2) composing the transgressive wedge, enabled observations on the depocentres modifications, which occurred during the latest sea level rise (Figure 11).As reported for similar deposits all over the world, we assisted with a general shore retreat and a main landward shift of the coastal systems [15,57,[80][81][82][83].Furthermore, a southward migration of the river mouth, as suggested by the different distribution of the shallow-marine depocentres (U2) with respect to the paralic/deltaic ones (U3), suggests channel mobility processes at the river mouth.Two main depocentres, NE-SW oriented, bounding a relative high, were observed by the 3D-like modelling of the U3 units.On the contrary, the thickness distribution of the marine sub-units (U2) reflected the presence of one main axis with the same direction perpendicular to the coast; in this case, shifted southwards.Avulsion processes could be induced by different factors, including tectonics, differences in substrate lithologies, and channel progressive infilling [84,85].As mentioned above, the soft and distributed fault linkages might have controlled the distribution of the paralic/deltaic sub-units at the edge of the developing relative topographic high (Figure 5b).An increased sediment yield and intense reworking of the coastal environment, probably driven by high volcanoclastic supply during volcanic paroxysm, almost coeval to the post-glacial initial stage sea level rise, has been reported in the literature [7,28].The combined effect of a relatively uplifting central area with a strong input in sedimentation probably induced a decrease in the accommodation space in that sector, forcing avulsion processes.We also propose an evolutionary model showing the disposition of the main modelled transgressive depocentres (Figure 11).During Stage 1 (Figure 11a), the deformation of the acoustic substratum accounted for the maximal thickness of the paralic/deltaic lobe in the northern sector of the studied area in a NE-SW trending direction perpendicular to the coast.In Stage 2 (Figure 11b), the delta front started to be shaped by erosional processes after the lowering of the tectonic activity and a new site, shifted southward, became the favourable substrate to accommodate the shallow marine deposition.

Conclusions
The late Quaternary evolution of the northern Campania continental shelf was affected by the emplacement of magmatic bodies, degassing features, and minor faulting.Small-scale deformations defined a structure characterized by two opposite-facing sets of faults, separated by a relative topographic high located in the overstepping region, about 11.4 km off the Volturno river mouth.Outward facing NW-SE and ENE-WSW opposite-dipping normal faults identified a structural horst, controlling in extension the transgressive units, as testified by the locally observed growth pattern of the paralic/deltaic transgressive bodies, referring to the initial stages of the latest sea-level rise.The interbasinal ridge, displayed as a belt at the overlapping fault systems terminations, described a roughly E-W oriented zone, accommodating the deformations.We followed the definition proposed in the literature [67] of accommodation zones as structures in regions of overlaps between separated systems of faults as well as accommodating deformations.We also used the geometrical classification proposed by the authors for such zones, interpreting the observed deformational pattern as the result of what should correspond to a small scale antithetic oblique synclinal accommodation zone-like type.Systems of NW-SE and ENE-WSW faults which define the structure should correspond to fault planes reactivated in a different kinematic context, as revealed by their trending direction as well as match the orientation of the most important discontinuities controlling the Tertiary to Quaternary evolution of the area [31,33,35].
Major late Quaternary NW-SE normal faults with quite modest displacements have been reported by several authors [30,33,35].The accommodation zone, individuated by a roughly E-W-oriented relative high, has a slightly oblique orientation with respect to the NE-SW extension trending direction, revealed by the main NNW-SSE trending set of faults and as reported by data from the literature [33,35].The Campi Flegrei caldera collapse, as a major dramatic event occurring at about 15 ka [41,50], should have locally perturbed and rearranged the stress field distribution.Evidence of such a process also could be identified in the renewed location of deposition of the shallow marine transgressive deposits when compared to the position of the paralic transgressive lobes.The combined effect of a relatively uplifting central area with a strong input in sedimentation probably induced a decrease in the accommodation space in that sector, possibly forcing the shallow marine depocentres to shift southwards.
The degassing features, displayed by seismic sections as vertically stacked amplitude anomalies affecting the late Quaternary outer shelf units, coincide with pockmark fields at the seafloor, distributed along an ESE-WNW trending direction and interpreted as the superficial response of a deeper laterally extensive acoustic turbidity front with the same striking direction.This laterally extensive fluids pocket, marked by distinct lateral fronts and an irregular top boundary and arranged in a peculiar saw-toothed configuration, was probably controlled by sharp variations in the permeability of outer shelf succession.The main trending direction of this front is similar to the one observed for the buried volcanic body [37,60], resulting in an E-W-trending flat-topped edifice.Magmatic bodies trigger thermal effects on the hosting sedimentary sequences and produce the expulsion of pore-water volumes [86][87][88].This observation suggests that the observed fluid front is a hydrothermal-related feature.

Figure 2 .
Figure 2. (a) Digital terrain model off the Volturno River: b, c, d boxes refer to the location of different subsets; (b-d) details of the round-shaped or elongated depressions detected at the shelf-edge and marking the sites of fluid seapages; diameter sizes range between 30 and 70 m for the rounded features and 200-300 m for the elongated features.

Figure 3 .
Figure 3. (a) Stratigraphy and magnetic susceptibility curve of the C5 core (redrawn from [54]); (b) seismic facies calibration; and (c) the location on seismic section SP1 of the deepest stratigraphic horizon intercepted by the core.

Figure 4 .
Figure 4. Three main horizons (H1, H2, and H3) were seismically detected on Chirp-Sonar profiles and mapped.They bound Unit 1, Unit 2, and Unit 3, related to the post-glacial deposition on the shelf.

Figure 5 .
Figure 5. (a) Four Chirp-Sonar profiles (1, 2, 3, 4) revealed the presence of NW dipping normal faults on the northern side of the interpreted grid, forming a downthrown region in a tilted block configuration.In the southeastern side of the grid, trending ENE-WSW, SE-dipping normal faults defined a relative morphological low; disrupted reflectors and columnar zones of acoustic masking depicted the sites of vertically focused fluid flows reaching the seabed and correspond to pockmarks (1 and 3); (b) 3D-like modelling of the H3 erosional surface revealed the presence of a relative high, bounded by two sets of outward-facing faults of limited offsets; and (c) the location of the interpreted profiles.

Figure 6 .
Figure 6.(A)The Sp1 Sparker profile shows the stacking pattern of the nine uppermost seismic units of the continental margin and the zone of the inner shelf with shallow gas (total acoustic blanking); Frame i marks the detail of the intermittent acoustic masking zone and the location of its saw-toothed top profile, shown in Figure9; (B) the Sp3 Sparker profile shows nine seismic units and the peculiar top profile of the acoustic turbidity front, which develops for about 10 km, striking the coast; Frame ii marks the detail of the acoustic turbidity zone, shown in Figure9; and (C) the zone of total acoustic blanking is semi-circular in plan view and centred at the Volturno River mouth.

Figure 7 .
Figure 7. (a) The vertically focused zones of rising fluids corresponded to extensive fields of pockmarks at the seabed.The top profile of the pocket developed as a saw-toothed configuration, probably related to vertical and lateral variations in permeability of the sequences, composing the succession; and (b) the detail of the acoustic turbidity front (location in Figure 6, Frame i); the red box evidences the locus of a vertically focused flow, possibly active as a seepage up to recent times, inhibiting the deposition of late Holocene sediments and resulting in a morphological step (see also Figure 2b).

Figure 8 .
Figure 8. Block diagram of the H3 erosional surface deformed by two opposite facing set of faults, separated by a relative topographic high, located in the overstepping region (offsets not to scale).

Figure 9 .
Figure 9. (A) Extension of the acoustic turbidity zone in Sparker lines (Sp1, Sp2 and Sp4) showing a remarkable shape of the top boundary, marked by enhanced acoustic reflections and intermittent acoustic masking; and (B) in dashed yellow is the lateral boundary of biogenic gas at the Volturno River mouth and, in red, the extension along the Sparker profiles of the inferred hydrothermal fluid front.

Figure 10 .
Figure 10.Conduits and pipes located immediately below the outer shelf pockmark field are sourced by a deeper fluid/gas pocket, with a main E-W trending direction.The morphology of the fluids' front depicts sharp and distinct lateral flanks and an irregular top profile arranged in a saw-toothed configuration; (1) volcanics, (2) acoustic turbidity front, (3) deposits related to the latest sea level rise, occurring from 18 ka BP to 5 ka BP, and (4) present-day highstand muddy deposits.

Figure 11 .
Figure 11.(a) Isopach map of the transgressive paralic/deltaic sub-unit; (b) isopach map of the transgressive shallow marine/prodeltaic sub-unit; (c) surface of the LGM erosional surface (H1); (d) H2 surface bounding the paralic from the shallow marine deposits of the Upper Pleistocene-Holocene transgressive systems tract; (e) top of the transgressive systems tract; and (f) two-step model (qualitative) of the depositional systems shift during transgression (Stage 1 is on the left and Stage 2 is on the right of the figure); 1, present-day coastline and plain (the present day location of the Volturno River mouth is at the seaward bulge of the coast); 2, palaeo-delta plain; 3, paralic/deltaic deposits (U3): in dark yellow, the location of the detected paralic depocentres; 4, shallow marine deposits (U2): in dark rose, the location, slightly shifted southward of the main marine depocentres; 5, slope area; 6, fault controlling the settlement of paralic/deltaic deposits; and 7, paleo-coastline location (not to scale).