Integrated Morpho-Bathymetric, Seismic-Stratigraphic, and Sedimentological Data on the Dohrn Canyon (Naples Bay, Southern Tyrrhenian Sea): Relationships with Volcanism and Tectonics

: Submarine canyons are geomorphologic lineaments engraving the slope / outer shelf of continental margins. These features are often associated with signiﬁcant geologic hazard when they develop close to densely populated coastal zones. The seaﬂoor of Naples Bay is deeply cut by two incisions characterized by a dense network of gullies, namely the Dohrn and Magnaghi canyons, which develop from the shelf break of the Campania margin, down to the peripheral rise of the Eastern Tyrrhenian bathyal plain. Seismic-stratigraphic interpretation of multichannel seismic reﬂection proﬁles has shown that quaternary tectonics and recent to active volcanism have exerted a signiﬁcant control on the morphological evolution and source-to sink depositional processes of the Dohrn and Magnaghi submarine canyons. The Dohrn canyon is characterized by relatively steep walls hundreds of meters high, which incise a Middle-Late Pleistocene prograding wedge, formed by clastic and volcaniclastic deposits associated with the paleo-Sarno river system during the Mid-Late Pleistocene. The formation of the Dohrn canyon predates the onset of the volcanic eruption of the Neapolitan Yellow Tu ﬀ (NYT), an ignimbrite deposit of ca. 15 ka that represents the bedrock on which the town of Napoli is built. Integrated stratigraphic analysis of high-resolution seismic proﬁles and marine gravity core data (C74_12) collected along the ﬂanks of the eastern bifurcation of the head of Dohrn Canyon suggests that depositional processes along the canyon ﬂanks are dominated by gravity ﬂows (e.g., ﬁne-grained turbidites, debris ﬂows) and sediment mass transport associated with slope instability and failure. feeds the eastern branch the shelf break of the Sorrento A remarkable rectilinear shape of the channels in plain suggests a structural control on the morphology of the drainage


Introduction
The continental slope and outer shelf of the Naples Bay are deeply cut by two submarine incisions-the Dohrn and the Magnaghi canyons-that form a large drainage system of this coastal marine area during the Late Quaternary. The Bay is located on the Eastern margin of the Tyrrhenian Sea, a Neogene-Quaternary extensional domain that includes small areas floored with oceanic-type crust. The Tyrrhenian Sea encompasses a number of peripheral ("peri-Tyrrhenian") basins, along its borders, like Naples Bay [1][2][3][4], that exhibit remarkable differences in size, thickness, and stratigraphic architecture.
Canyons and gullies are major components of continental slope morphology, as they represent a critical link between coastal waters and abyssal depths, by transferring sediments, nutrients and even litter and pollutants into deep-water settings. The aim of this research is to provide a contribution toward the understanding of depositional and erosional processes along a major canyon system developing on a tectonically active continental margin, associated with active volcanism, during the

Material and Methods
The acquisition of geological and geophysical data presented in this study has been carried out within the frame of a national project (CARG) for the geological mapping of Italy at 1:50,000 scale. The morphobathymetry of the Naples Bay has been obtained as a mosaic of Multibeam Ecosounder (MBES) surveys conducted by the Institute of Marine Sciences of National Research Council of Naples, Italy, mostly between 2000 and 2014, using various MBES systems. The processing of multibeam data was carried out with PDS2000 software (Teledyne Marine, Cypress, TX, USA) following the International Hydrographic Organization standards [83]. Processing included the removal of navigation errors, noise reduction (i.e., despiking), removal of poor-quality beams, and tidal and sound velocity corrections. The final digital terrain model (DTM) derived from the MBES data covers an area of 900 km 2 with a cell size of 10 m. Onshore topographic data were obtained from the official topographic grid of the Military Institute for Geography (20 m grid cell), acquired between 1985 and 1990 by aerial photogrammetry [84] (Figure 1).
The geological framework of the shallow structures has been investigated through seismic and sequence stratigraphic interpretation of a grid of multichannel and single-channel, high resolution (Sparker 1 kJ and 4 kJ) reflection seismic profiles. Selected sub-bottom Chirp profiles were also used to document the detailed stratigraphic architecture of channel-levee complexes and gravity flow deposits. The geological framework of the shallow structures has been investigated through seismic and sequence stratigraphic interpretation of a grid of multichannel and single-channel, high resolution (Sparker 1 kJ and 4 kJ) reflection seismic profiles. Selected sub-bottom Chirp profiles were also used to document the detailed stratigraphic architecture of channel-levee complexes and gravity flow deposits.
The analysis of seismic reflection profiles was conducted according with the methods and procedures of seismic and sequence stratigraphic interpretation [85][86][87][88][89][90][91]. Seismic facies of stratigraphic units identified on seismic profiles were described particularly based on the internal and external geometries of individual units and on the continuity, amplitude, and frequency of reflectors [92][93][94][95][96]. Relative sea-level fluctuations were analyzed by constructing chronostratigraphic diagrams and curves of relative sea-level cycles [90,91].
The interpretation of reflection seismic data was calibrated with integrated stratigraphic analysis of gravity core C74_12, collected at a water depth of 586.6 m, along the Dohrn canyon ( Figure 1). Core samples were half-split, photographed, and described at cm-scale using a 10× lens. Microscope observation was conducted on selected samples of sieved sand fraction to identify the major lithologic components (lithoclasts, bioclasts, volcaniclasts) and determine the shape and surface texture of particles. Sedimentological analysis included the recognition of major lithofacies associations, sedimentary structures, and grain-size analysis of the fraction <63 μ by laser diffractometry (Sympatec). A total number of 81 samples have been processed for the grain-size analysis. Statistical parameters have been calculated according to the classic graphical equations of Folk and Ward (1957) [97].
The analysis of sediment samples has been complemented by high-resolution (2 cm step) measurements of physical properties of core C74_12 to assist stratigraphic correlation. Core logging was conducted on half-split cores, using a fully automated multi sensor core logger (MSCL-Geotek) (Geotek, Northamptonshire, UK) at the CNR-ISMAR petrophysical laboratory. The MSCL system includes (i) a Bartington MS2E Point sensor, to measure the low-field volume magnetic susceptibility; (ii) a gamma ray attenuation porosity evaluator (GRAPE) sensor to determine the bulk density; (iii) an ultrasonic P-wave (PWTs) system to analyze the P-wave velocity through the core; and (iv) a The analysis of seismic reflection profiles was conducted according with the methods and procedures of seismic and sequence stratigraphic interpretation [85][86][87][88][89][90][91]. Seismic facies of stratigraphic units identified on seismic profiles were described particularly based on the internal and external geometries of individual units and on the continuity, amplitude, and frequency of reflectors [92][93][94][95][96]. Relative sea-level fluctuations were analyzed by constructing chronostratigraphic diagrams and curves of relative sea-level cycles [90,91].
The interpretation of reflection seismic data was calibrated with integrated stratigraphic analysis of gravity core C74_12, collected at a water depth of 586.6 m, along the Dohrn canyon ( Figure 1). Core samples were half-split, photographed, and described at cm-scale using a 10× lens. Microscope observation was conducted on selected samples of sieved sand fraction to identify the major lithologic components (lithoclasts, bioclasts, volcaniclasts) and determine the shape and surface texture of particles. Sedimentological analysis included the recognition of major lithofacies associations, sedimentary structures, and grain-size analysis of the fraction <63 µ by laser diffractometry (Sympatec). A total number of 81 samples have been processed for the grain-size analysis. Statistical parameters have been calculated according to the classic graphical equations of Folk and Ward (1957) [97].
The analysis of sediment samples has been complemented by high-resolution (2 cm step) measurements of physical properties of core C74_12 to assist stratigraphic correlation. Core logging was conducted on half-split cores, using a fully automated multi sensor core logger (MSCL-Geotek) (Geotek, Northamptonshire, UK) at the CNR-ISMAR petrophysical laboratory. The MSCL system includes (i) a Bartington MS2E Point sensor, to measure the low-field volume magnetic susceptibility; (ii) a gamma ray attenuation porosity evaluator (GRAPE) sensor to determine the bulk density; (iii) an ultrasonic P-wave (PWTs) system to analyze the P-wave velocity through the core; and (iv) a Minolta Spectrophotometer CM 2002 to detect the reflectance (L%,) acquired within 1 h after the core splitting to minimize redox-associated color changes.

Geological Setting
The Gulf of Naples stretches in the southern part of a tectonic depression, namely the Campania Plain, whose formation was controlled by the Neogene back-arc extension of the Tyrrhenian Basin that accompanied the northeastward-verging accretion of the Southern Apennines during the roll-back of the subducting Adriatic plate [98][99][100][101][102][103][104][105]. The Early Pleistocene marks the onset of NW-SE extension and tectonic subsidence of the Campania Plain-Napoli Bay basin, whereas a following phase (Middle Pleistocene in age) mainly resulted in NE-SW extension. Normal faults associated with these two phases were later reactivated during the Late Pleistocene-Holocene and caused either subsidence or tectonic uplift, especially in areas affected by volcano-tectonic deformation [106][107][108][109][110][111][112].
The regional geological structure of the Naples Bay is characterized by prevailing NE-SW and ENE-WSW trends of normal faults, namely expressed in the directions of the Dohrn canyon and the Capri-Sorrento lineament ( Figure 2) [113]. These faults have often accommodated significant offsets within the Meso-Cenozoic sedimentary succession, representing the acoustic basement on the seismic profiles in the area.
The regional geological structure of the Naples Bay is characterized by prevailing NE-SW and ENE-WSW trends of normal faults, namely expressed in the directions of the Dohrn canyon and the Capri-Sorrento lineament ( Figure 2) [113]. These faults have often accommodated significant offsets within the Meso-Cenozoic sedimentary succession, representing the acoustic basement on the seismic profiles in the area.
The present-day morpho-structural setting of the Campania continental margin and the Naples Bay is the result of significant interaction between volcanic and sedimentary processes during the Late Quaternary [114]. The Dohrn canyon fault bounds the south-western flank of the Banco di Fuori, with a total offset between the top of the structural high NW of the fault and the top of the downthrown Meso-Cenozoic carbonates below the Dohrn canyon, SE of the fault, in the order of one thousand meters. The enhanced thickness of the Middle-Late Pleistocene prograding wedge associated with the Dohrn canyon fault suggests a coeval age for tectonic activity along this structure, which is then buried by relatively undeformed Holocene deposits.  The present-day morpho-structural setting of the Campania continental margin and the Naples Bay is the result of significant interaction between volcanic and sedimentary processes during the Late Quaternary [114]. The Dohrn canyon fault bounds the south-western flank of the Banco di Fuori, with a total offset between the top of the structural high NW of the fault and the top of the downthrown Meso-Cenozoic carbonates below the Dohrn canyon, SE of the fault, in the order of one thousand meters. The enhanced thickness of the Middle-Late Pleistocene prograding wedge associated with the Dohrn canyon fault suggests a coeval age for tectonic activity along this structure, which is then buried by relatively undeformed Holocene deposits.
Based on previous seismic data, three regional geological sections have been constructed across the Campania continental margin (Figure 3) [113], showing the main morpho-structures occurring at a regional scale ("Banco di Fuori" morpho-structural high, Dohrn canyon, Magnaghi canyon, Capri Basin, Salerno Valley, Volturno basin). branches, which converge into a NE-SW thalweg, bounded by the Capri Basin to the South. The Magnaghi canyon, located to the West of the Dohrn canyon, is characterized by diffused erosional processes along its flanks accompanied by the occurrence of volcaniclastic gravity flow deposits along the main axis of the valley, mostly fed by the products of the eruptive activity of the Ischia and Procida islands during the Late Quaternary. The Capri Basin forms a deep basin located in the Tyrrhenian bathyal plain south of the Naples Bay and filled by a several hundred meters of Pleistocene-Holocene sediments overlying the Meso-Cenozoic carbonate acoustic basement. The Salerno Valley is a half-graben basin characterized by three seismic units composed of Quaternary marine and continental deposits, laterally resting above chaotic deposits of the Cilento Group. Off the northern sector of the Campania Plain, the Volturno Basin is filled up by a thick succession of Quaternary marine siliciclastic, volcaniclastic and deltaic units, overlying a substrate of Meso-Cenozoic carbonates and Miocene deposits (sand and shale).  [113], showing the main morpho-structures on a regional scale (Banco di Fuori morpho-structural high, Dohrn canyon, Magnaghi canyon, Capri Basin, Salerno Valley, Volturno Basin).
Campi Flegrei is volcanic district along the coastline of the Naples Bay, where volcanism has been active since ~60 ka B.P. [18,[21][22][23][24]32,33,[115][116][117][118][119][120][121][122]. The area was characterized by a large caldera forming eruption, which caused the emplacement the Campanian Ignimbrite, (ca. 39 ka B.P.). The regressive peak of the last glacial maximum (isotopic stage 2) resulted in ca. 30 km progradation of the shoreline that reached the present-day location of the shelf break in the Gulf of Naples, close to the 140 m isobath. The glacio-eustatic emersion was accompanied by significant volcaniclastic input  [113], showing the main morpho-structures on a regional scale (Banco di Fuori morpho-structural high, Dohrn canyon, Magnaghi canyon, Capri Basin, Salerno Valley, Volturno Basin). "Banco di Fuori" is a morpho-structural high composed of Meso-Cenozoic carbonates bounding the Naples Bay to the south. The flanks and top of this high are overlain by a thick succession of Pleistocene deposits. The Dohrn canyon is a major erosional lineament that broadly divides the Naples Bay into an eastern sector, mostly characterized by siliciclastic deposits, and a western sector where volcanic and mixed volcaniclastic-siliciclastic units prevail. The canyon is composed of two branches, which converge into a NE-SW thalweg, bounded by the Capri Basin to the South. The Magnaghi canyon, located to the West of the Dohrn canyon, is characterized by diffused erosional processes along its flanks accompanied by the occurrence of volcaniclastic gravity flow deposits along the main axis of the valley, mostly fed by the products of the eruptive activity of the Ischia and Procida islands during the Late Quaternary. The Capri Basin forms a deep basin located in the Tyrrhenian bathyal plain south of the Naples Bay and filled by a several hundred meters of Pleistocene-Holocene sediments overlying the Meso-Cenozoic carbonate acoustic basement. The Salerno Valley is a half-graben basin characterized by three seismic units composed of Quaternary marine and continental deposits, laterally resting above chaotic deposits of the Cilento Group. Off the northern sector of the Campania Plain, the Volturno Basin is filled up by a thick succession of Quaternary marine siliciclastic, volcaniclastic and deltaic units, overlying a substrate of Meso-Cenozoic carbonates and Miocene deposits (sand and shale).
Previous interpretations have suggested that the Campanian Ignimbrite possibly originated along a NW-SE trending normal fault (Napoli-Vico Equense normal fault) [124] that caused the tilting of coastal plain toward the NE. The tilting block was bounded by NE-SW trending normal faults, located along the northern coast of the Sorrento Peninsula and along the opposite side of the Gulf of Naples. A significant lineament pertaining to this structural trend is represented by the Acerra-Dohrn canyon fault, which extends onshore from the Acerra graben to the Dohrn canyon offshore [124][125][126][127]. During the last decades a series of detailed studies on the Campania Ignimbrite deposits and the crustal mechanisms controlling the emplacement of the ignimbrites of the Campania Volcanic Zone (CVZ) have been presented [32,[123][124][125][128][129][130][131].
Absolute age determinations on the products of the CI yielded ca. 37 ka BP [131,132], 39 ky [130], and 39.85 ± 0.14 ka B.P. [133], whereas the NYT has been dated to ca. 15 ka B.P. [75]. In the northern CVZ, the stratigraphic data from more than 600 boreholes [19,134,135] allowed the construction of geological sections showing the stratigraphic architecture of the ignimbrite deposits [123] (Figure 4). Figure 4 shows that the lowermost unit is composed of marine sediments (Tyrrhenian Layer M1), whose top has been dated at 126 ka B.P. [134]. A volcanic unit overlies the marine deposits and is composed of massive ignimbrite deposits and lava flows. This unit is overlain by another succession of marine deposits, dated at 55-50 ka (Tyrrhenian unit M2) [134], which is covered in turn by the CI deposits. The Giugliano Ignimbrite (GI) follows above the CI deposits. The intermediate CVZ, whose stratigraphy is described in Figure 4 [123], is characterized by the occurrence of the NYT caldera, in its western sector, whereas its eastern sector occupies the area close to the town of Naples. whose top has been dated at 126 ka B.P. [134]. A volcanic unit overlies the marine deposits and is composed of massive ignimbrite deposits and lava flows. This unit is overlain by another succession of marine deposits, dated at 55-50 ka (Tyrrhenian unit M2) [134], which is covered in turn by the CI deposits. The Giugliano Ignimbrite (GI) follows above the CI deposits. The intermediate CVZ, whose stratigraphy is described in Figure 4 [123], is characterized by the occurrence of the NYT caldera, in its western sector, whereas its eastern sector occupies the area close to the town of Naples.

Morpho-Bathymetry and Gravity Instabilities of the Naples Bay Canyons
An extensive, high-resolution bathymetric survey of the continental shelf/slope system of the Campania region (southern Italy) has been carried out in the frame of research programs of marine geological mapping. The relative bathymetric data were acquired between 1997 and 2017, using multibeam systems with an average vertical resolution of <0.25% water depth and a position accuracy of <10 m. The survey data were successively merged with a Digital Terrain Model (DTM) created from topographic maps of the Naples Bay onshore coastal area and islands, to produce a Digital Elevation Model (DEM) based on a homogeneous grid with cell-spacing of 20 m ( Figure 5) [78]. Major morphological elements of the Naples Bay include: (1) the Dohrn-Magnaghi canyon system, carving the continental slope at water depths between −250 m and −1100 m [77][78][79][80][136][137][138][139]; (2) the continental slope system of the Ischia island and related submarine canyons [74]; (3) the onshore and offshore volcanoes of the Campi Flegrei (Pentapalummo, Nisida, and Miseno banks) [18,21,28,[30][31][32][33]; (4) the rough sea bottom area of the outer shelf of the city of Naples (Banco della Montagna) [78,81,[137][138][139]; (5) the sediment wave field of the inner continental shelf off the Vesuvius [34,35,78,81,[137][138][139][140].
the continental slope at water depths between −250 m and −1100 m [77][78][79][80][136][137][138][139]; (2) the continental slope system of the Ischia island and related submarine canyons [74]; (3) the onshore and offshore volcanoes of the Campi Flegrei (Pentapalummo, Nisida, and Miseno banks) [18,21,28,[30][31][32][33]; (4) the rough sea bottom area of the outer shelf of the city of Naples (Banco della Montagna) [78,81,[137][138][139]; (5) the sediment wave field of the inner continental shelf off the Vesuvius [34,35,78,81,[137][138][139][140]. The Dohrn canyon represents the most important morphological feature of the Naples Bay and develops into two branches that converge into a main valley. The thalweg of the canyon is relatively flat and displays a variable width (200 m-1 km). The canyon slopes are relatively steep, sometimes terraced, and are characterized by several instability phenomena, such as slide scars, slumps, and slides. On the southeastern flanks of the canyon system, at least four rectilinear, sub-parallel incisions occur, likely because of a structural control by ENE-WSW trending faults ( Figure 5). The Dohrn canyon represents the most important morphological feature of the Naples Bay and develops into two branches that converge into a main valley. The thalweg of the canyon is relatively flat and displays a variable width (200 m-1 km). The canyon slopes are relatively steep, sometimes terraced, and are characterized by several instability phenomena, such as slide scars, slumps, and slides. On the southeastern flanks of the canyon system, at least four rectilinear, sub-parallel incisions occur, likely because of a structural control by ENE-WSW trending faults ( Figure 5). The Naples canyon system originates from the shelf break off the Campi Flegrei volcanic district along the 140 m isobath. The general morphological setting of the canyon system is controlled by two major morpho-structural lineaments of the Naples offshore, namely the Banco di Fuori structural high, bounding the canyon system to the South and the Capri structural high, bounding the Dohrn canyon to the East, next to its confluence with the bathyal plain. Both the Magnaghi canyon and the western branch of the Dohrn canyon show morphological evidence of retreatment of the canyon heads. Particularly, the Dohrn shows two retrogressive head areas, as documented by the multibeam bathymetry and morpho-bathymetric profiles.

Geomorphologic Map of Naples Canyons
Based on marine DTM derived by MBES data, a geomorphologic interpretation of the Dohrn canyon and of the continental slope of the Naples bay was obtained. Major morphological lineaments of the canyon system include the main drainage axes, slide scars, and associated gravity instability areas, relict volcanic edifices, shelf break, terrace rims, and gravity flow deposits (e.g., turbidite slope fans; Figure 6). The Dohrn canyon is fed by several tributary channels, which can be readily identified from morpho-bathymetric interpretation. Two main tributary channels originate from the shelf break of the Naples Bay, located approximately along the 140 m isobath and cross the entire continental slope between the Dohrn and Magnaghi canyons. A system of tributary channels, developing at water depths between 200 m and 500 m, feeds the Dohrn eastern branch from the shelf break of the Sorrento Peninsula. A remarkable rectilinear shape of the channels in plain view suggests a structural control on the morphology of the drainage system. areas, relict volcanic edifices, shelf break, terrace rims, and gravity flow deposits (e.g., turbidite slope fans; Figure 6). The Dohrn canyon is fed by several tributary channels, which can be readily identified from morpho-bathymetric interpretation. Two main tributary channels originate from the shelf break of the Naples Bay, located approximately along the 140 m isobath and cross the entire continental slope between the Dohrn and Magnaghi canyons. A system of tributary channels, developing at water depths between 200 m and 500 m, feeds the Dohrn eastern branch from the shelf break of the Sorrento Peninsula. A remarkable rectilinear shape of the channels in plain view suggests a structural control on the morphology of the drainage system.

Bathymetric Profiles
Two bathymetric profiles perpendicular to the main axis of the Dohrn canyon have been constructed (Figures 7 and 8). At the canyon head, the Dohrn drainage system is characterized by two main tributary channels. The western tributary is deeper and more incised than the eastern one (P1; Figure 8). Profile-P2 represents the sector of the canyon system where the two tributary channels merge to form a single valley involved by severe erosional processes, particularly to the NW and W of the Capri Basin [139].
Two bathymetric profiles perpendicular to the main axis of the Dohrn canyon have been constructed (Figures 7 and 8). At the canyon head, the Dohrn drainage system is characterized by two main tributary channels. The western tributary is deeper and more incised than the eastern one (P1; Figure 8). Profile-P2 represents the sector of the canyon system where the two tributary channels merge to form a single valley involved by severe erosional processes, particularly to the NW and W of the Capri Basin [139].  Two bathymetric profiles perpendicular to the main axis of the Dohrn canyon have been constructed (Figures 7 and 8). At the canyon head, the Dohrn drainage system is characterized by two main tributary channels. The western tributary is deeper and more incised than the eastern one (P1; Figure 8). Profile-P2 represents the sector of the canyon system where the two tributary channels merge to form a single valley involved by severe erosional processes, particularly to the NW and W of the Capri Basin [139].  The central-eastern part of the Naples Bay is characterized by the occurrence of the Ammontatura channel, a relatively minor incision of the Dohrn canyon system characterized by a curved morphology in plain view, a flat thalweg and asymmetrical (locally irregular) slopes. Ammontatura channel originates from a narrow area between the Nisida island and the Nisida bank and terminates at the head of the western tributary channel of the Dohrn canyon (Figures 9 and 10). It also separates the "Banco della Montagna" structure from the volcanic margins of the Campi Flegrei Bathymetric profiles indicate that the western slope of the Ammontatura Channel is generally steeper than the western slope, where evidence of a series of terraced surfaces is also documented. The valley is 2.5 km wide and ca 20-40 m deep. In its northernmost part, the axis of the channel bends toward the north-west and terminates abruptly N of the Nisida Bank. The central-eastern part of the Naples Bay is characterized by the occurrence of the Ammontatura channel, a relatively minor incision of the Dohrn canyon system characterized by a curved morphology in plain view, a flat thalweg and asymmetrical (locally irregular) slopes. Ammontatura channel originates from a narrow area between the Nisida island and the Nisida bank and terminates at the head of the western tributary channel of the Dohrn canyon (Figures 9 and 10). It also separates the "Banco della Montagna" structure from the volcanic margins of the Campi Flegrei Bathymetric profiles indicate that the western slope of the Ammontatura Channel is generally steeper than the western slope, where evidence of a series of terraced surfaces is also documented. The valley is 2.5 km wide and ca 20-40 m deep. In its northernmost part, the axis of the channel bends toward the north-west and terminates abruptly N of the Nisida Bank.   The central-eastern part of the Naples Bay is characterized by the occurrence of the Ammontatura channel, a relatively minor incision of the Dohrn canyon system characterized by a curved morphology in plain view, a flat thalweg and asymmetrical (locally irregular) slopes. Ammontatura channel originates from a narrow area between the Nisida island and the Nisida bank and terminates at the head of the western tributary channel of the Dohrn canyon (Figures 9 and 10). It also separates the "Banco della Montagna" structure from the volcanic margins of the Campi Flegrei Bathymetric profiles indicate that the western slope of the Ammontatura Channel is generally steeper than the western slope, where evidence of a series of terraced surfaces is also documented. The valley is 2.5 km wide and ca 20-40 m deep. In its northernmost part, the axis of the channel bends toward the north-west and terminates abruptly N of the Nisida Bank.

Channel-Levee Complex
Seismic facies interpreted as channel-levee complexes have been identified on subbottom Chirp profiles crossing the Ammontatura channel ( Figure 11). This facies is characterized by alternating sub-parallel seismic reflectors and acoustically transparent intervals, separating an incision, lacking filling, and corresponding with the channel. This facies is composed of gravity flow (e.g., turbidite, debris flow) deposits, characterized by alternating sands and muds, that can be interpreted as channel-levee complexes in the areas surrounding the Dohrn canyon and the related tributary channels.

Channel-Levee Complex
Seismic facies interpreted as channel-levee complexes have been identified on subbottom Chirp profiles crossing the Ammontatura channel ( Figure 11). This facies is characterized by alternating sub-parallel seismic reflectors and acoustically transparent intervals, separating an incision, lacking filling, and corresponding with the channel. This facies is composed of gravity flow (e.g., turbidite, debris flow) deposits, characterized by alternating sands and muds, that can be interpreted as channel-levee complexes in the areas surrounding the Dohrn canyon and the related tributary channels.  Channel-levee complexes have been also documented by Chirp profiles NA18 ( Figure 12) and GP23 ( Figure 13). Profile NA18 shows the continental shelf of the Naples Bay, the Ammontatura channel, a slope morphological high and the Pentapalummo Bank, with its SE slope cut by a tributary channel ( Figure 12). The Naples continental shelf is characterized by wavy, sub-parallel and continuous seismic reflectors, whereas the Ammontatura channel shows a sequence of chaotic seismic reflectors, overlain by parallel seismic reflectors. The slope morphological high is characterized by wavy and sub-parallel continuous seismic reflectors and the Pentapalummo bank shows an acoustically transparent seismic facies ( Figure 12). Sub-parallel and continuous seismic reflectors on the continental shelf of the Naples Bay have been interpreted as Holocene marine deposits, whereas the chaotic sequence recognized in the Ammontatura channel has been interpreted as corresponding to slump deposits, overlain by sub-parallel Holocene marine strata ( Figure 12). Wavy and sub-parallel reflectors forming most of the slope morphological high have been interpreted as channel-levee complexes, whereas the acoustically transparent seismic facies on the Pentapalummo bank have been regarded as representing the volcanic acoustic basement (Figure 12). reflectors on the continental shelf of the Naples Bay have been interpreted as Holocene marine deposits, whereas the chaotic sequence recognized in the Ammontatura channel has been interpreted as corresponding to slump deposits, overlain by sub-parallel Holocene marine strata ( Figure 12). Wavy and sub-parallel reflectors forming most of the slope morphological high have been interpreted as channel-levee complexes, whereas the acoustically transparent seismic facies on the Pentapalummo bank have been regarded as representing the volcanic acoustic basement ( Figure 12). Seismic profile GP23 crosses a relict morphology of the Late Pleistocene continental shelf, located in the central slope of the Naples Bay and the continental slope towards the thalweg of the Dohrn eastern branch (Figure 13). The relict morphological high is characterized by seismic facies showing prograding clinoforms systems, truncated at the top by an angular unconformity and overlain by parallel seismic reflectors ( Figure 13). The continental slope shows inclined seismic reflectors and significant detachment scours (Figure 13), whereas at the base-of-the-slope truncated reflectors and thin sediment drapes draping on the canyon walls are also documented ( Figure 13; refer to the B rectangle on the figure). The prograding clinoforms detected on the relict morphological high have been interpreted as Middle-Late Pleistocene continental shelf deposits, truncated by an erosional unconformity and overlain by Holocene marine layers ( Figure 13). Thin Holocene marine deposits also drape older strata truncated along the channel walls ( Figure 13). Seismic profile GP23 crosses a relict morphology of the Late Pleistocene continental shelf, located in the central slope of the Naples Bay and the continental slope towards the thalweg of the Dohrn eastern branch (Figure 13). The relict morphological high is characterized by seismic facies showing prograding clinoforms systems, truncated at the top by an angular unconformity and overlain by parallel seismic reflectors ( Figure 13). The continental slope shows inclined seismic reflectors and significant detachment scours (Figure 13), whereas at the base-of-the-slope truncated reflectors and thin sediment drapes draping on the canyon walls are also documented ( Figure 13; refer to the B rectangle on the figure). The prograding clinoforms detected on the relict morphological high have been interpreted as Middle-Late Pleistocene continental shelf deposits, truncated by an erosional unconformity and overlain by Holocene marine layers ( Figure 13). Thin Holocene marine deposits also drape older strata truncated along the channel walls ( Figure 13).

Gravity Flow Deposits
Seismic units represented by lenticular bodies with chaotic reflections, poor lateral continuity of seismic reflectors and a turbid acoustic facies ( Figure 11) are often found sandwiched between acoustically sharp, wavy inclined reflectors, and have been interpreted as corresponding to gravity flow (e.g., debris flow/turbidite) deposits formed along the continental slope.
A sketch table was constructed in order to summarize the most significant seismic facies detected in the Naples Bay on sub-bottom Chirp profiles and their corresponding seismic stratigraphic interpretation ( Figure 14). It should be noted that the table is based on previously published data [142] and includes a wider spectrum of seismic facies than that specifically discussed in this work. Chaotic reflectors with poor lateral continuity and transparent to turbid facies have been interpreted as volcaniclastic and/or pyroclastic deposits (mud flows, lahars, pyroclastic flows) and are genetically related to the Somma-Vesuvius volcano (A in Figure 14). Convex and acoustically transparent units, partly exposed at the sea floor, alternating with wavy or parallel reflectors with a clear acoustic appearance have been interpreted as volcaniclastic mounds composed of alternating coarse-grained pumice, volcanic sand and mud (B in Figure 14). These features occur as a field of small mounds (Banco della Montagna) that develop over the continental shelf of the Naples Bay and have been recently interpreted as the result of dragging and uplift of unconsolidated volcaniclastic deposits from deeper stratigraphic layers due to the degassing of overpressured volcanic fluids originated from deeper structural levels [81]. Convex units with fair lateral continuity of with a transparent acoustic facies, resting above high amplitude reflectors have been recently interpreted as Holocene marine sediments involved by gravitational instability and slow creep over a basal sliding surface likely corresponding to a maximum flooding surface (MFS) (C in Figure 14). Slightly inclined acoustically clear, parallel reflections, passing upwards to an area with a turbid acoustic appearance have been interpreted as outcrops of relict sands, deposited during the last sea level lowstand (18-20 ky B.P.), which is genetically related to the last 4th order eustatic cycle (D in Figure 14). Acoustically transparent bodies (E in Figure 14), respectively, correspond to the dome-shaped volcanic units genetically related to (i) Campanian Ignimbrite (E1 in Figure 14); (ii) volcanic acoustic basement detected at Pentapalummo, Miseno, and Nisida banks (E2 in Figure 14); and (iii) volcanic basement interpreted as lava flows of parasitic vents genetically related to the Somma-Vesuvius volcano (E3 in Figure 14). Gravity flow (debris flows/turbidite) deposits have been documented as lenticular bodies with chaotic laterally discontinuous reflectors with a turbid acoustic appearance (F in Figure 14). Channel-levee complexes have been characterized as a succession of convex and sub-parallel or inclined reflectors with good lateral continuity, with channel and embankment morphologies (G in Figure 14). A succession of inclined and wavy reflectors with a good lateral continuity and sharp to turbid acoustic characters have been interpreted as a relict morphology of the Middle-Late Pleistocene continental shelf, which is genetically related to the Sarno prograding wedge (H in Figure 14). Geosciences 2020, 10, x FOR PEER REVIEW 16 of 37

Seismic Stratigraphy
Erosion and the transport of sediments in the Bay of Naples are chiefly driven by the Dohrn and Magnaghi canyons, which expose along its walls hundreds of meters of a Middle-Late Pleistocene prograding wedge represented by siliciclastic and volcaniclastic deposits, genetically related to the paleo-Sarno fluvial system. The Dohrn canyon is separated from the Magnaghi canyon by a major morpho-structural high ("Banco di Bocca Grande" or "Banco di Fuori"), likely formed as a result of

Seismic Stratigraphy
Erosion and the transport of sediments in the Bay of Naples are chiefly driven by the Dohrn and Magnaghi canyons, which expose along its walls hundreds of meters of a Middle-Late Pleistocene prograding wedge represented by siliciclastic and volcaniclastic deposits, genetically related to the paleo-Sarno fluvial system. The Dohrn canyon is separated from the Magnaghi canyon by a major morpho-structural high ("Banco di Bocca Grande" or "Banco di Fuori"), likely formed as a result of the local uplift and tilting of the basement rocks (Mesozoic carbonates) that are widely exposed in the Sorrento Peninsula-Capri Island structural high.
The location of the DEB (Dohrn eastern branch) and several lines of evidence deriving from seismic stratigraphic interpretation have suggested a genetic linkage between the activity of the Dohrn canyon and the Sarno paleo-drainage system during sea-level lowstands. A 30-40-m-thick seismic sequence overlying older, undisturbed reflectors and including channel-levee complexes, crops out at the sea floor on both the sides of the DEB. These deposits do not seem related to the main branches, which are deeply incised, but rather than the occurrence of tributary channels, controlling the overflow of sediment fluxes in the surrounding areas.
The deltaic system of the paleo-Sarno river has directly fed the continental slope areas during the lowstand phases of the Middle-Late Pleistocene, producing a large and thick prograding wedge, deeply incised by the Dohrn canyon ( Figure 15). We infer that the DEB was directly linked with the paleo-Sarno river mouth during lowstand phases. This is also suggested by the occurrence of two non-volcanic morphobathymetric highs detected on Chirp profiles close to the heads of the Dohrn canyon, (Figure 13). Seismic-stratigraphic analysis shows they represent deltaic sediments arranged with clinoform patterns and suggest they are related to the distal part of the prograding wedge fed by the paleo-Sarno river mouth. These features reach an elevation of 20-30 m with respect to the surrounding seafloor water depth and may be interpreted as relict morphologies of the Middle-Late Pleistocene continental shelf ( Figure 16). The origin of these morpho-bathymetric high was likely the result of differential erosion during the subaerial exposure of the continental shelf following sea level lowstand. Seismic interpretation also suggests this erosional phase predates the eruption of the Campanian Ignimbrite age (39 ka B.P.).
Geosciences 2020, 10, x FOR PEER REVIEW 17 of 37 the local uplift and tilting of the basement rocks (Mesozoic carbonates) that are widely exposed in the Sorrento Peninsula-Capri Island structural high. The location of the DEB (Dohrn eastern branch) and several lines of evidence deriving from seismic stratigraphic interpretation have suggested a genetic linkage between the activity of the Dohrn canyon and the Sarno paleo-drainage system during sea-level lowstands. A 30-40-m-thick seismic sequence overlying older, undisturbed reflectors and including channel-levee complexes, crops out at the sea floor on both the sides of the DEB. These deposits do not seem related to the main branches, which are deeply incised, but rather than the occurrence of tributary channels, controlling the overflow of sediment fluxes in the surrounding areas.
The deltaic system of the paleo-Sarno river has directly fed the continental slope areas during the lowstand phases of the Middle-Late Pleistocene, producing a large and thick prograding wedge, deeply incised by the Dohrn canyon ( Figure 15). We infer that the DEB was directly linked with the paleo-Sarno river mouth during lowstand phases. This is also suggested by the occurrence of two non-volcanic morphobathymetric highs detected on Chirp profiles close to the heads of the Dohrn canyon, (Figure 13). Seismic-stratigraphic analysis shows they represent deltaic sediments arranged with clinoform patterns and suggest they are related to the distal part of the prograding wedge fed by the paleo-Sarno river mouth. These features reach an elevation of 20-30 m with respect to the surrounding seafloor water depth and may be interpreted as relict morphologies of the Middle-Late Pleistocene continental shelf ( Figure 16). The origin of these morpho-bathymetric high was likely the result of differential erosion during the subaerial exposure of the continental shelf following sea level lowstand. Seismic interpretation also suggests this erosional phase predates the eruption of the Campanian Ignimbrite age (39 ka B.P.). Figure 15. Geologic interpretation of the Sparker seismic profile SPK8. Note that the Mesozoic carbonate acoustic basement of the Naples Bay is overlain by several stratigraphic sequences, including the pre and syn-rift prograding wedge, Early-Middle Pleistocene in age and the post-rift prograding wedge, Middle Pleistocene in age. C and D represent main regional unconformities.  Selected single-channel (Sparker) profiles have been interpreted in order to illustrate the main seismic-stratigraphic units recognized in the Naples Bay (SP20, SP20_1, SP2, SP3; . Seismic line SP20 documents the stratigraphic setting of the continental shelf are of the Naples Bay, that is characterized by three main seismic units ( Figure 17). The lowest unit is represented by truncated prograding clinoforms and has been interpreted as a relict prograding wedge fed by the Sarno river mouth, probably Middle-Late Pleistocene in age. This succession is overlain by a thick acoustically transparent seismic unit, likely corresponding with the Campanian Ignimbrite, [15][16][17][18][19][20][21][22][23][24][25][26][27]29], whose submarine occurrence in the Naples Bay has been recently mapped by Aiello et al., 2016 [31]. The upper seismic unit is characterized by parallel-continuous seismic reflectors, that are interpreted as a succession of Pleistocene-Holocene deposits covering the erosional surface at the top of the relict prograding wedges and the Campanian Ignimbrite units (Figure 17). Seismic line SP20 documents the stratigraphic setting of the continental shelf are of the Naples Bay, that is characterized by three main seismic units ( Figure 17). The lowest unit is represented by truncated prograding clinoforms and has been interpreted as a relict prograding wedge fed by the Sarno river mouth, probably Middle-Late Pleistocene in age. This succession is overlain by a thick acoustically transparent seismic unit, likely corresponding with the Campanian Ignimbrite, [15][16][17][18][19][20][21][22][23][24][25][26][27]29], whose submarine occurrence in the Naples Bay has been recently mapped by Aiello et al., 2016 [31]. The upper seismic unit is characterized by parallel-continuous seismic reflectors, that are interpreted as a succession of Pleistocene-Holocene deposits covering the erosional surface at the top of the relict prograding wedges and the Campanian Ignimbrite units (Figure 17). Seismic line SP20_1 illustrates the seismic stratigraphic architecture of the Naples Bay from the Banco di Fuori morpho-structural high towards the tributary channels of the Dohrn canyon ( Figure  18). The seismic-stratigraphic setting of the Banco di Fuori high is characterized by the occurrence of turbiditic slope systems with mostly chaotic seismic facies, overlain by marine deposits, displaying sub-parallel reflectors, probably Late Pleistocene-Holocene in age. Tributary channels of the Dohrn canyon are incised in thick turbiditic slope systems, which form the bulk of the stratigraphic succession of the continental slope area, along with channel levee complexes ( Figure 18). Seismic line SP20_1 illustrates the seismic stratigraphic architecture of the Naples Bay from the Banco di Fuori morpho-structural high towards the tributary channels of the Dohrn canyon ( Figure 18). The seismic-stratigraphic setting of the Banco di Fuori high is characterized by the occurrence of turbiditic slope systems with mostly chaotic seismic facies, overlain by marine deposits, displaying sub-parallel reflectors, probably Late Pleistocene-Holocene in age. Tributary channels of the Dohrn canyon are incised in thick turbiditic slope systems, which form the bulk of the stratigraphic succession of the continental slope area, along with channel levee complexes ( Figure 18). Geosciences 2020, 10, x FOR PEER REVIEW 20 of 37 Figure 18. Geologic interpretation of the Sparker seismic profile SP20_1, exhibiting the seismostratigraphic setting of the Dohrn canyon next to the Banco di Fuori morpho-structural high.
Seismic line SP2 shows the stratigraphic relationships between the volcanic seismic units of the Naples Bay; namely, the Campanian Ignimbrite (CI) and the Neapolitan Yellow Tuff (NYT).

Figure 18.
Geologic interpretation of the Sparker seismic profile SP20_1, exhibiting the seismo-stratigraphic setting of the Dohrn canyon next to the Banco di Fuori morpho-structural high.
during the isotopic stage 3 [8], whereas the Neapolitan Yellow Tuff was deposited during the upper part of the isotopic stage 2 [8] corresponding to the development of the Transgressive Systems Tract (TST) in the Naples Bay [77,110,125]. The marine sequence partly formed during the lower part of the isotopic stage 2 [8], corresponding with the Lowstand Systems Tract (LST) in the Naples Bay ( Figure  19) [77,110,125,136]. Nevertheless, this unit also includes a significant volcaniclastic component, as indicated by the occurrence of volcanic domes interbedded in the LST sequence ( Figure 19). Additional information on the stratigraphic relationships between the major volcaniclastic units of the Naples Bay (namely CI and NYT) have been obtained by the Sparker seismic line SP3 ( Figure  20). In this area, the lowstand deposits are missing, and the NYT unit directly overlies the CI unit. The NYT seismic unit is rooted in the Banco della Montagna structure, a wide volcanic field occurring in the Naples Bay continental shelf [81,137]. In this area the NYT unit is overlain by a succession of marine and coastal sediments, Late Pleistocene-Holocene (<15 ka) in age ( Figure 20).

Sedimentological Data
Gravity core C74_12 has been sampled collected at water depths of 586.6 m in the lower southern part of the Dohrn western branch, close to the confluence with the Dohrn eastern branch (Figures 1  and 21). The recovered core sample displays a total length of ca 7.2 m. The core description was performed as follows, from the top (A) to the bottom by a homogeneous, bioturbated, hemipelagic mud (clayey silt and silty clay 65%) (H) of the core (Figure 21). Starting from the sea bottom, considered as a zero level a 1.5 m thick silty interbedded with coarse to very fine silty sand layers Seismic line SP2 shows the stratigraphic relationships between the volcanic seismic units of the Naples Bay; namely, the Campanian Ignimbrite (CI) and the Neapolitan Yellow Tuff (NYT). Particularly, this seismic profile highlights the occurrence of a thick marine sequence intercalated between the CI and the NYT seismic units (Figure 19). The Campanian Ignimbrite was deposited during the isotopic stage 3 [8], whereas the Neapolitan Yellow Tuff was deposited during the upper part of the isotopic stage 2 [8] corresponding to the development of the Transgressive Systems Tract (TST) in the Naples Bay [77,110,125]. The marine sequence partly formed during the lower part of the isotopic stage 2 [8], corresponding with the Lowstand Systems Tract (LST) in the Naples Bay ( Figure 19) [77,110,125,136]. Nevertheless, this unit also includes a significant volcaniclastic component, as indicated by the occurrence of volcanic domes interbedded in the LST sequence ( Figure 19).
Additional information on the stratigraphic relationships between the major volcaniclastic units of the Naples Bay (namely CI and NYT) have been obtained by the Sparker seismic line SP3 (Figure 20). In this area, the lowstand deposits are missing, and the NYT unit directly overlies the CI unit. The NYT seismic unit is rooted in the Banco della Montagna structure, a wide volcanic field occurring in the Naples Bay continental shelf [81,137]. In this area the NYT unit is overlain by a succession of marine and coastal sediments, Late Pleistocene-Holocene (<15 ka) in age ( Figure 20).

Sedimentological Data
Gravity core C74_12 has been sampled collected at water depths of 586.6 m in the lower southern part of the Dohrn western branch, close to the confluence with the Dohrn eastern branch (Figures 1  and 21). The recovered core sample displays a total length of ca 7.2 m. The core description was performed as follows, from the top (A) to the bottom by a homogeneous, bioturbated, hemipelagic mud (clayey silt and silty clay 65%) (H) of the core (Figure 21). Starting from the sea bottom, considered as a zero level a 1.5 m thick silty interbedded with coarse to very fine silty sand layers even a few decimeters thick, made of mostly volcaniclasts (pumice, scoriae, glass shards, and fragments of crystals) that probably indicate tephra layers (e), and bioclasts (forams, pteropods, ostracods, mollusks fragments) (sandy/gravelly sediments 35%) ( Figure 21) bioturbated interval has been found ( Figure 21). Based on a sedimentological study of the cored succession and the results of statistical grain-size analysis of the sediment, three main lithofacies associations (units) have been recognized (Figures 21  and 22). These are represented, from bottom to top by the following units [143].
-Facies A (720-530 cm)-From base to top, it consists of semi-consolidated, silty sand including lens of small, mostly rounded pumice and bioclasts with evidence of reworking. At 625 cm, a decimeter thick level of fining upward medium coarse sand, with an erosive base, consisting mainly of volcaniclasts (scoriae, minerals, and glass shards); -Facies B (530-240 cm)-It consists of clayey silt with bioclasts (planktonic and benthonic foraminifera) interbedded with decimetric to centimetric layers and lenses of parallel and cross laminated, medium and fine sandy silt, often with erosive base and normal grading. The coarser fraction is mostly represented by gray pumice, scoriae, glass shard, well preserved minerals (e.g., sanidine, peroxide, biotite). At a depth of 530 cm, a 55 cm-thick layer of coarse pumice with sub-angular to sub-rounded clasts occurs with reverse grading (e10). Coarse-grained components of these deposits are formed by moderately sorted, sub-angular volcaniclastic fragments (pumice, scoria, K-feldspar, pyroxene, volcanic glass). These deposits pass upward to 30 cm of silty sandy pumice with rare shallow water faunas (benthic foramifers, shell fragments). At a depth of 425 cm follows a 20 cm-thick layer, with erosional base, including sub-angular to sub-rounded pumice and occasional iron-oxidized patinas, rare pelitic pebble, and reworked fragments of Mollusk shells (e9). Toward the top follow centimetric lenses (403 cm, 378 cm, 252 cm, and 249 cm) and fining upward and/or parallel and crossed laminated decimetric levels (438 cm, 300 cm) of volcaniclasts, (e4-e8); -Facies C (240-0 cm)-It is composed of homogeneous grey mud, represented by bioturbated clayey silt and clayey sandy silt with a few interbedded layers rich in volcaniclasts (170-160 cm, 160-153 cm, 43-37 cm) (e1-e3) and/or bioclasts (153-135 cm, 60-50 cm).

Physical Properties and Stratigraphic Correlation
The major petrophysical properties of the sampled sediments along the core logs have been analyzed in order to identify macroscopic variations and general patterns for each analyzed parameter and assist stratigraphic correlation at decimetric scale. Particular emphasis has been placed on Magnetic Susceptibility which is one of the most important parameters for the recognition to correspond to a fair correlation with GRAPE density maxima. Continuous high-resolution logs based on physical properties measurements of the cores have been acquired on a split-core by a fully automated Multi Sensor Core Logger (Geotek) at the CNR-ISMAR of Naples (Italy) petrophysical laboratory. The measured parameters are the magnetic susceptibility, the GRAPE density, the compressional P-wave velocity, and the spectral reflectance.

Discussion
Submarine canyons represent important erosive lineaments of the seafloor morphology. They deeply cut long belts of the continental slopes and provide a main mode for the transfer of coastal sediments mode through which both siliciclastic and volcaniclastic materials coming from the coastal zones are transported towards the deep oceanic basins. Previous research on Neogene evolution of the main continental margins of the Mediterranean, North Atlantic and Pacific areas has confirmed the common occurrence of large palaeo-channels, often inter-layered in the Quaternary deposits of the sedimentary basins, suggesting that the origin of these kinds of canyons could be genetically related with the subaerial environment [56,[144][145][146][147][148][149][150].
The main parameters in the study of submarine canyons include (a) the depth range, (b) the slope gradient, (c) the relationships between canyon morphology and tectonic lineaments, (d) the

Physical Properties and Stratigraphic Correlation
The major petrophysical properties of the sampled sediments along the core logs have been analyzed in order to identify macroscopic variations and general patterns for each analyzed parameter and assist stratigraphic correlation at decimetric scale. Particular emphasis has been placed on Magnetic Susceptibility which is one of the most important parameters for the recognition of tephras and cryptotephras in mixed siliciclastic-volcaniclastic environments like the Naples Bay. Measured parameters are (1) magnetic susceptibility, (2) Gamma-ray attenuation porosity evaluator (GRAPE) density, (3) spectral reflectance, and (4) P-wave velocity. Corresponding core logs and correlation are presented in Figure 22. Continuous high-resolution (2-cm step) logs, based on measurements of physical properties of the cored sediments, have been acquired on a split-core by a fully automated Multi Sensor Core Logger (Geotek) at the petrophysical laboratory of CNR-ISMAR, Naples.
Measured parameters for the core C74_12 are: (a) magnetic susceptibility, (b) GRAPE density, (c) spectral reflectance (L%), and (d) P-wave velocity are presented in Figure 22. As a first step, the general patterns for each parameter were correlated with the lithostratigraphic log in order to identify the macroscopic facies changes at centimeter scale. Successively, the most reliable parameters for bedded siliciclastic and volcaniclastic successions (e.g., magnetic susceptibility) were analyzed in detail, trying to pinpoint single prominent peaks, as well as characteristic patterns. Finally, the preliminary parameter-based correlation was checked versus the sedimentological logs and other measured parameters.
The analysis of plot diagrams indicates that the values of volume magnetic susceptibilities display values between 2 and 1036 *10-5 SI, with six and show six ranges (indicated by peaks 1-6) where values reach maxima between 350 and 1036 *10-5 SI. All six peaks correlate to stratigraphic levels in which there are lithological variations with respect to embedding sediments. Particularly, peaks 1, 3, and 5 correspond to sandy silt with fine-grained pumice and scoria, with evidence of basal erosion, and representing partly reworked pyroclastic fall volcaniclastic deposits. On the other hand, peaks 4 and 6, mark normally graded, undisturbed, fine-grained pumice, and scoria sandy levels with scoria overlying muddy intervals and showing absence of erosional surface which could be interpreted as probable undisturbed volcanic sandy levels corresponding with primary pyroclastic fall deposits. Peak 2 finally appears to correlate with the base of a gravity flow deposit (e.g., grain flow/upper turbidite). Texturally homogeneous deposits (mainly composed of silts), display almost constant values of magnetic susceptibility values of volume, while the interval between 458-524 cm, which is characterized by the lowest values, consists almost exclusively of coarse pumice level of massive.
GRAPE density values swear in a range of 1.1 g/cc to 2.1 g/cm 3 g/cc all along the wel corel. The trend pattern of variations of this parameter seems appears to follow mirror that of the trend of the magnetic susceptibility values log, including a with the main values fall remarkable minimum of (form 11.6 gr/cc to 1.1 g/cm 3 gr/cc) occurring again at 458-524 cm that correlates with so far the density drop corresponds to the beginning of the above mentioned coarse pumice interval, and corresponding to the magnetic susceptibility minimum, recorded between 458 cm and 524 cm and extending all over it.
Peaks 1-6 are also distinguishable in the density record, and may yield either apositive correlation with respect to magnetic volume susceptibilities, (in the case of there is a positive correlation for peaks 1, 3, and 5), while there is an anti-parallel or anti-correlation (e.g., for peaks 2, 4, and 6) with magnetic volume susceptibility.
Spectral reflectance indicates a general trend in the variation of lower and higher values of the change in sediment's lightness, and particularly interesting is to note that it shows a tendency to correlation with the variation of lower and higher trend values in correspondence of magnetic susceptibility peaks, which could indicate a relative abundance between darker or/and lighter volcanic sediments depending by their different contents (e.g., (pumices/sanidine or scoriaceous/lithics enrichness).
Finally, it should be noted that the quality of the P-wave velocity log data is not fully recorded, not due to the partial loss of signal (the signal is lost, perhaps as a consequence of the presence of voids, in the coarse-grained volcaniclastic interval). On the other hand, where the signal was well recorded within intervals 0-400 cm and 540-720 cm, it is interesting to note that where the trend of the variations is similar to that one of the density log, particularly between 180 cm and 340 cm the P-wave velocity peaks 3 (1538 m/s) and 4 (maximum speed 1538 m/s and 1694 m/s respectively) seem to correspond to a fair correlation with GRAPE density maxima. Continuous high-resolution logs based on physical properties measurements of the cores have been acquired on a split-core by a fully automated Multi Sensor Core Logger (Geotek) at the CNR-ISMAR of Naples (Italy) petrophysical laboratory. The measured parameters are the magnetic susceptibility, the GRAPE density, the compressional P-wave velocity, and the spectral reflectance.

Discussion
Submarine canyons represent important erosive lineaments of the seafloor morphology. They deeply cut long belts of the continental slopes and provide a main mode for the transfer of coastal sediments mode through which both siliciclastic and volcaniclastic materials coming from the coastal zones are transported towards the deep oceanic basins. Previous research on Neogene evolution of the main continental margins of the Mediterranean, North Atlantic and Pacific areas has confirmed the common occurrence of large palaeo-channels, often inter-layered in the Quaternary deposits of the sedimentary basins, suggesting that the origin of these kinds of canyons could be genetically related with the subaerial environment [56,[144][145][146][147][148][149][150].
The main parameters in the study of submarine canyons include (a) the depth range, (b) the slope gradient, (c) the relationships between canyon morphology and tectonic lineaments, (d) the shape of canyon section, i.e., U-shape (depositional) and V-shape (erosive), and (e) the canyon profile (rectilinear or meandering) [56].
Several studies on the Dohrn canyon have been carried out during the last decades, mostly regarding the distribution of volcanic centers of the Naples Bay area and regional tectonic lineaments, as well as the seismic stratigraphic architecture of the continental margin of the Campania region and associated depositional environments offshore [77,78,138,151,152]. Other research has focused on the morpho-bathymetric setting of the Dohrn and Magnaghi canyons the and submarine gravity instabilities [81,137] and the potential run-up and tsunami reconstruction from fossil submarine landslides [79].
Previous studies have suggested that, in the submarine sector of the Naples Baym the deposition of thick volcanic units (CI and NYT) was coeval with the formation of the Dohrn canyon, suggesting this process was acting during the eustatic sea level fall and the tectonic uplift of the outer shelf domains during the Late Quaternary [77]. Other authors have proposed that the Dohrn and the Magnaghi canyons may be related to a much older, inherited drainage system genetically linked to the latest Neogene-Pleistocene tectonic evolution of the Campania continental margin [151,152].
The Magnaghi canyon and the Dohrn western branch show morphological indication of retreating canyon heads. This evidence is essential for the understanding of the submarine canyon formation, since the sediment flows may have repeatedly triggered retrogressive failures, according to the commonly accepted models on submarine canyon evolution [54]. The process of canyon formation can be interpreted by reconciling morphological evidence for headward canyon erosion by mass wasting with the stratigraphic evidence for canyon inception by downslope-eroding sediment flows. The canyon morphology and evolution could be summarized by hypothesizing a scenario from "youthful" to "mature" canyon morphology [64]. The youthful stage of submarine canyon evolution begins with extensive slope failures. Retrogressive mass wasting of continental slope sediments along the failure headwall leads to the formation and upslope extension of relatively straight, steep-walled chute. If this headward migrating chute breaches the shelf break, the canyon taps into a new sediment source of outer shelf sands and enters into a mature phase of canyon evolution. The failure of shelf sediments in the vicinity of the canyon head initiates coarse-grained gravity flow processes (e.g., turbidity currents, debris/grain flows), which become an important agent in canyon erosion by downcutting the canyon and creating a sinuous thalweg [54,64].
Previous interpretations have reported that there is no evidence of a high sediment supply along the shelf edge of the Naples Bay during the Late Quaternary, as the sediment originating from the emerged areas was possibly deposited on the continental shelf, as a consequence of tectonic subsidence [77]. The results of our study indicate, on the contrary, that during the late Pleistocene sediment supply (particularly volcaniclastic input) was extremely high in the area, due to the eruption of the Campanian Ignimbrite (CI) and the Neapolitan Yellow Tuff (NYT). Moreover, the considerable amount of available sediment also played a major role in the building of the stratigraphic architecture of the Naples Bay and the formation of the canyon system.
The stratigraphic analysis of gravity core C74_12 documents a succession of primary and mostly reworked tephra layers interbebbed with hemipelagic deposits. Tephra layers ostensibly represent periods of significant volcaniclastic input and sediment transport along the axis of the Dohrn canyon, whereas hemipelagic intervals (i.e., bioturbated mud) denote phases of quiescence of sediment transport throughout the canyon. Volcaniclastic input mainly consists of sediment gravity flows (e.g., turbidite and debris/grain flow), likely originated in the upper part of the canyon from massive accumulation of primary pyroclastic deposits on the outer shelf and upper slope of the Naples Bay that become mixed with hemipelagic material during the sedimentary transport downslope (e.g., facies C, layers e2 and e3; facies B, layers e9 and e10) (Figures 21 and 22). In some cases, gravity flow deposits stack up one upon another, suggesting the occurrence of periods of enhanced activity in sediment transport and a high frequency of re-sedimentation processes.
A thin drape of hemipelagic sediments, detected on seismic profiles and confirmed by sedimentological analysis, appears to cover the whole Dohrn canyon system, suggesting a substantial present-day inactivity of the Dohrn canyon, likely as a consequence of the recent decrease of volcaniclastic input into the Naples Bay and the concomitant high sea level stillstand. In core C74_12 the recent phase of quiescence of the Dohrn canyon is evidenced by ca. 1.5 m thick interval of bioturbated silt (hemipelagite) immediately beneath the seafloor ( Figure 21).
Our interpretation confirms the suggestion that the Dohrn canyon is genetically related with the main eruptive events involving the Naples Bay during the latest Pleistocene, including the CI; (ca. 39 ka B.P.) [21] and the NYT (ca. 15 ky B.P.) [75], recently correlated with main tephra levels occurring in the Northern Campi Flegrei offshore [76]. The significant volcaniclastic supply during the eruptive phases occurring in the Neapolitan area during this period has probably controlled the individuation of the Dohrn canyon western branch. Similarly, we infer that the individuation of the Dohrn eastern branch occurred by the joining of a paleo-valley located offshore the Sorrento Peninsula, likely as a consequence of high sedimentary supply from the Schiazzano river mouth [77].
The results of our research are in agreement with the geological model recently presented for the Pozzuoli Bay by Steinmann et al. (2018) [33], which includes a scenario for the formation of the Dohrn canyon. Based on this model a possible kinematic reconstruction of the northern part of the Naples Bay during the last 15 ka, can be summarized into five main stages (Figures 23 and 24). The paleo-coastlines for each stage have been inferred through an integrated analysis of seismic and bathymetric data [78], onshore outcrops [20], and borehole stratigraphic columns [20,21].
-Stage 1 (15 ka) During stage 1 the sea level stationed at ca. −93 m below the present-day level and the main volcano-tectonic event was represented by the caldera collapse related to the NYT eruption (Figures 23 and 24). The central part of the NYT caldera collapse was in the submarine domain, whereas outside of the caldera depression, subaerial conditions prevailed in the coastal area. Unconformity U2 correlates with the top of the NYT deposits and represents a partly submarine and partly subaerial feature. This unconformity also marks a channel-like topographic low, genetically related with the erosional truncation of the M1 seismic sequence (Figure 24), corresponding to the actual location of the Dohrn Canyon and its continuation toward the north-west. The Miseno bank was partly emerged, and formed a peninsula connected to the mainland, whereas the Pentapalummo bank was an island in northern Naples Bay. During this stage, the Dohrn canyon extended toward the north-west and drained significant amounts of volcaniclastic materials originated during the NYT eruption toward the outer shelf and slope.  (Figures 23 and 24) and at least six magmatic to phreato-magmatic eruptions occurred. The sea level rise resulted in additional landward retreat of the shoreline outside of the NYT caldera structure. Early during this stage the La Starza terrace (Pozzuoli Bay) progressively emerged due to the ground uplift associated with the NYT caldera resurgence and, at same time, marine unconformities U3 and U4 formed as a response to the consequent decrease in the accommodation space over the continental shelf (Figures 23 and 24).  (Figures 23 and 24). Since the beginning of this stage a significant phase of subsidence, affected the inner sector of the Pozzuoli Bay as testified by dramatic drowning of coastal settlements and infrastructures of Roman age. Subsidence in the range of 20 m to 8 m, with average rates of 10-4 mm/year also led to the drowning of the coastal shallow-water terrace associated with unconformity U6 and related prograding clinoforms. This created significant accommodation space, resulting in the deposition of the present day coastal prograding wedge ( Figure 23).    [33], showing the main marine units and the corresponding regional unconformities (U2-U6), that correlate with the main volcano-tectonic events (see Figure 23).

Conclusions
Two deeply incised canyons drain the continental shelf and the upper slope of the Bay of Naples (Southern Italy), an active volcanic area since Middle-Late Pleistocene. Erosion and transport of volcaniclastic sediments from the Campi Flegrei and Procida Island have acted mainly along the Magnaghi canyon and the Dohrn western branch, whereas a consistent fluvial supply coming from the Sarno river fed the Dohrn eastern branch. The Magnaghi trilobate canyon head displays retrogressive failures engraving the continental shelf. Sedimentary processes are dominated by gravity flows (turbidite and debris/grain flows) mostly moving downslope toward the southwest due to the local control of a morpho-structural high represented by Meso-Cenozoic basement rocks.
Based on the results of seismic stratigraphic interpretation, the Dohrn canyon is definitely older than the volcanic deposits of the Neapolitan Yellow Tuff (NYT) (ca. 15 ka B.P.) that form the backbone of the town of Naples [21,75].
The walls of Dohrn canyon were affected by slope failures mostly along the western flanks, where large slide scars occur. Material displaced by slope failure processes was effectively removed and redistributed toward the southeast by dominant gravity flows (turbidite and debris/grain flows) Sedimentological analysis of gravity core C74_12 revealed the occurrence, along the axis of Dohrn canyon of reworked, K-feldspar-rich pumiceous tephra layers and volcaniclastic gravity flow deposits interbedded with bioturbated muddy intervals (hemipelagite) marking phases of quiescence in the canyon activity. Gravity flow processes were accompanied by the overflow of sediment spills over levees of tributary channels, producing overbank deposits along the canyon sides. Fossil tributary channels suspended over the main branches testify for phases of rapid re-incision switching  [33], showing the main marine units and the corresponding regional unconformities (U2-U6), that correlate with the main volcano-tectonic events (see Figure 23).

Conclusions
Two deeply incised canyons drain the continental shelf and the upper slope of the Bay of Naples (Southern Italy), an active volcanic area since Middle-Late Pleistocene. Erosion and transport of volcaniclastic sediments from the Campi Flegrei and Procida Island have acted mainly along the Magnaghi canyon and the Dohrn western branch, whereas a consistent fluvial supply coming from the Sarno river fed the Dohrn eastern branch. The Magnaghi trilobate canyon head displays retrogressive failures engraving the continental shelf. Sedimentary processes are dominated by gravity flows (turbidite and debris/grain flows) mostly moving downslope toward the southwest due to the local control of a morpho-structural high represented by Meso-Cenozoic basement rocks.
Based on the results of seismic stratigraphic interpretation, the Dohrn canyon is definitely older than the volcanic deposits of the Neapolitan Yellow Tuff (NYT) (ca. 15 ka B.P.) that form the backbone of the town of Naples [21,75].
The walls of Dohrn canyon were affected by slope failures mostly along the western flanks, where large slide scars occur. Material displaced by slope failure processes was effectively removed and redistributed toward the southeast by dominant gravity flows (turbidite and debris/grain flows) Sedimentological analysis of gravity core C74_12 revealed the occurrence, along the axis of Dohrn canyon of reworked, K-feldspar-rich pumiceous tephra layers and volcaniclastic gravity flow deposits interbedded with bioturbated muddy intervals (hemipelagite) marking phases of quiescence in the canyon activity. Gravity flow processes were accompanied by the overflow of sediment spills over levees of tributary channels, producing overbank deposits along the canyon sides. Fossil tributary channels suspended over the main branches testify for phases of rapid re-incision switching off the feeding from lateral sources. A thin drape of hemipelagic muddy sediments covering the whole canyon system, suggests its present-day inactivity, probably due to the decreasing volcaniclastic supply and to the high sea-level stillstand.
The instability phenomena documented along the Dohrn canyon slopes may be interpreted as the result of a phases of rapid accumulation of volcaniclastic deposits on the continental shelf due to (a) climatically induced enhanced erosion and associated high sediment input in the Bay of Naples or (b) the deposition of large volumes of volcaniclastic sediments as a consequence of massive volcaniclastic input into the coastal sedimentary environment of the Naples Bay during the Late Pleistocene. Seismic stratigraphic interpretation calibrated with the sedimentological analysis of gravity core samples suggests an upper fan valley-levee depositional setting for this stratigraphic interval.