From Siliciclastic to Bioclastic Deposits in the Gulf of Naples: New Highlights from Offshore Ischia and Procida–Pozzuoli Based on Sedimentological and Seismo-Stratigraphic Data

: This study discusses the siliciclastic to bioclastic deposits (in particular, the rhodolith deposits) in the Gulf of Naples based on sedimentological and seismo-stratigraphic data. The selected areas are offshore Ischia Island (offshore Casamicciola, Ischia Channel), where a dense network of sea-bottom samples has been collected, coupled with Sparker Multi-tip seismic lines, and offshore Procida–Pozzuoli (Procida Channel), where sea-bottom samples are available, in addition to Sparker seismic proﬁles. The basic methods applied in this research include sedimentological analysis, processing sedimentological data, and assessing seismo-stratigraphic criteria and techniques. In the Gulf of Naples, and particularly offshore Ischia, bioclastic sedimentation has been controlled by seaﬂoor topography coupled with the oceanographic setting. Wide seismo-stratigraphic units include the bioclastic deposits in their uppermost part. Offshore Procida–Pozzuoli, siliciclastic deposits appear to prevail, coupled with pyroclastic units, and no signiﬁcant bioclastic or rhodolith deposits have been outlined based on sedimentological and seismo-stratigraphic data. The occurrence of mixed siliciclastic–carbonate depositional systems is highlighted in this section of the Gulf of Naples based on the obtained results, which can be compared with similar systems recognized in the central Tyrrhenian Sea (Pontine Islands).


Introduction
Rhodolith and the maërl deposits are genetically related to coralline algae, representing important tools for assessing the paleoclimatic and paleoenvironmental conditions [1][2][3][4]. Analysis of the structure of these algae is an important technique in order to investigate the paleoenvironmental conditions in different regions, including tropical, temperate, and polar regions. Coralline algae have some characteristics which make them important paleoenvironmental proxies, comprising the growth of annual and sub-annual growth bands, consisting of high-Mg calcite, and the progressive burial of new thalli from sediments, when new thalli grow, composing part of the dead deposit [1][2][3][4]. Moreover, their widespread distribution, extending from polar to tropical zones, makes them a key stratigraphic tool in order to reconstruct the environmental conditions during the geological past, from middle to high latitudes, where they can be compared with a few other recorders, to low latitudes, where many other recorders exist [1][2][3][4]. However, it has been highlighted that the geographical distribution of coralline algae on a global scale is scattered, if compared with the geographical distribution of the rhodolith and maërl deposits [1][2][3][4].
In carbonate platform settings, rhodolith and the maërl deposits are important facies [5][6][7][8]. On the Pontian shelf (Tyrrhenian Sea), the patterns of carbonate facies have been reconstructed in order to quantify the contribution of red coralline algae to carbonate sedimentation [5]. The studied deposits are mainly composed of gravelly sands. The maximum diameter of the algal particles, representing the rhodoliths or the unattached branches and constituting these deposits, ranges between 1 and 3 cm. The analysis of the thin sections of the collected samples has shown five main taxa of coralline algae, including Lithotamnion minervae, Lithothamnion sp., Titanoderma sp., Mesophyllum sp., and Lithophyllum sp. [5]. Two main carbonate facies have been identified, namely, calcareous algae facies and carbonate matrix facies [5]. Calcareous algae facies is characterized by a carbonate content higher than 80% and is mainly composed of red algae. Carbonate matrix facies has a carbonate content ranging between 42% and 84%, and individual carbonate grains are not distinguishable. These two carbonate facies are typical of distinct bathymetric intervals, the first one ranging between water depths of 40 m and 70 m, and the second one between 70 m and 100 m. Facies interpretation has highlighted that this area of the Pontian shelf can be interpreted as a mobile substratum of the circa-littoral zone ("Détritique Cotier"). In this area, red coralline algae play a major role in controlling carbonate production and sedimentation.
In the Gulf of Naples, the occurrence of bioclastic deposits controlled by red coralline algae has been highlighted by previous authors [14], as well as the distribution of deep rhodolith layers offshore the Campania region [15]. Previous studies identifying red coralline algae in the Gulf of Naples are the studies of Funk, 1927;1955 [16,17], Bacci, 1947 [18], Lewalle, 1961 [19] and Segre, 1967;1972 [20,21]. Although the studies by Segre have been realized in the frame of the construction of the first geological map of Italy (geological sheet "Napoli-Isola d'Ischia), the other papers have been published by the Stazione Zoologica "Anton Dohrn" of Naples (Italy).
One key aim of this paper is to study the bioclastic deposits (rhodolith deposits) in the northern Gulf of Naples (offshore Ischia) and to highlight that a transition from bioclastic to siliciclastic deposits exists towards offshore Procida-Pozzuoli based on sedimentological and seismo-stratigraphic data. In the study area of the Gulf of Naples, mixed carbonate-siliciclastic sedimentation exists. Both offshore Ischia [22][23][24][25][26][27][28][29] and offshore Procida-Pozzuoli [30][31][32][33][34] represent areas which were actively volcanic during the Late Quaternary period. The sedimentological and seismo-stratigraphic data have been analyzed in the general framework of the geological maps of Ischia and Procida Islands at a 1:10,000 scale [35,36].
The distribution and the characterization of rhodolith layers offshore Campania, and particularly in the Gulf of Naples (Capri, Punta Campanella, Bocca Piccola, and Ischia), have recently been discussed based on the results obtained during a research project performed within the Marine Strategy Framework Directive [15]. ROV investigations and grab samples collected at water depths ranging between 42 m (Capri) and 73 m (Ischia) have been coupled with multibeam and side-scan sonar acoustic imagery of the studied sites. The obtained results show that pralines are the most important morphological types between Capri, Punta Campanella, and Bocca Piccola, whereas at Ischia, there is a prevalence of pralines, but boxwork rhodoliths also occur [15]. All the selected sites are characterized by a dominance of the genus Lithothamnion, but at Ischia, the specimen Lithophyllum racemus has been documented [15].
In the Gulf of Naples, rhodolith deposits have been reported offshore Nisida-Posillipo, closely to the city of Naples, where they overlie rocky bottoms and surrounding areas of the mobile sea bottom in an inner shelf domain [14]. These deposits are composed of bioclastic gravels, bioclastic gravelly sands, and bioclastic sands. The pyroclastic component is locally abundant. The dimensions of the rhodoliths range between 15 and 20 mm. The red coralline algae grow up around the pyroclastics, producing the rhodoliths. In this area, the rhodolith deposits are mainly composed of dead thalli, which have been transported from the rocky outcrops and consequently deposited in the surrounding areas after reworking [14].
In the Gulf of Pozzuoli, rhodolith deposits have been found on the Miseno volcanic bank [14]. Here, living crustose thalli have been found, grading into palimpsest deposits composed of bioclastic sands. Moreover, maërl deposits have been found at water depths ranging between 35 and 58 m.
At Ischia, only one sample has been previously studied [14], while a dense network of samples located offshore Ischia has been recently analyzed [37] and has been re-interpreted in this paper and compared with the data from offshore Procida-Pozzuoli. The Sparker seismic profiles interpreted at Ischia [37] have been herein revised, re-interpreted, and compared with new seismo-acoustic data, including a Sparker seismic profile located in the Procida Channel.

Geologic and Oceanographic Setting
The Gulf of Naples is a Neogene-Quaternary basin of the peri-Tyrrhenian area, located between the Southern Apennines and the Tyrrhenian basin. In this area the extensional processes were active simultaneously with the compressional processes, leading to the formation of the Apenninic chain and the migration of the chain-foredeep system towards the Apulian foreland [38][39][40][41][42].
The Gulf of Naples is bounded southwards by the Sorrento Peninsula, eastwards by the southern sector of the Campania Plain, and north-westwards by the Phlegraean Fields ( Figure 1). The Sorrento Peninsula and the Capri island represent a WSW-ENE trending morpho-structural high, interposed between the Gulf of Naples-Campania Plain northwards and the Gulf of Salerno-Sele Plain southwards [55]. The morpho-structural high is mainly composed of Mesozoic carbonate rocks overlain by Miocene flysch deposits. The Campania Plain is a depression filled with alluvial deposits and volcanic products, whose thickness exceeds 3000 m [56][57][58].  The volcanologic setting of the Ischia and Procida-Pozzuoli volcanic districts is outlined. The Phlegraean volcanic district embraces the Ischia, Procida, and Vivara islands and the submarine volcanic conduits of the NW Naples and Pozzuoli Gulfs. Wide explosive eruptions, both Plinian and ignimbritic, occurred at Ischia, whose deposits constitute the stratigraphic unit known as the "Sintema del Rifugio di S. Nicola" [59]. The oldest eruptive activity (150-75 ky B.P.) is documented in outcrop in the southern sector of the island, in correspondence to the promontories located between Punta Imperatore and Carta Romana [60][61][62]. At Ischia, other important eruptive periods ranging in age between 98 and 73 ky B.P. and between 73 and 56 ky B.P. are responsible for the formation of the stratigraphic units of the Lower Ischia Channel and of the Pignatiello pyroclastics, after which the Ischia caldera started to form [62].
The geologic setting of the Procida island is characterized by volcanic deposits, respectively ascribed to the Vivara, Terra Murata, and Punta Serra volcanoes (age < 75 ky B.P.), dated back after radiometric absolute dating [35].
From 30 to 10 ky B.P. a volcanic activity characterized by a strong area diffusion has controlled the construction of the Solchiaro volcano (average age 22 ky B.P.) and the coeval deposition of white tuffs, outcropping in Soccavo and in the Naples town. Other main volcanological events were the construction of the Torregaveta volcano and the eruption of the Neapolitan Yellow Tuff, 15 ky ago [63]. Starting from 10 ky B.P. in the surrounding sectors, lava domes, scoria cones, and monogenetic pyroclastic volcanoes were established (Monte Nuovo, Astroni, Agnano, Capo Miseno, Bacoli, Montagna Spaccata, Monte Ruscello [64,65]. Based on the interpretation of seismic and core data, a collapse of the Phlegraean caldera in the Gulf of Pozzuoli has been hypothesized as controlled by a ring-fault system [33,34,66]. This system is coincident with a faulted area reaching a width of 2 km, surrounding an inner region of the caldera, reaching a width of 5 km [34]. A nested caldera complex, formed during and after the Campanian Ignimbrite and Neapolitan Yellow tuff eruptions, was recognized based on the interpretation of seismic profiles [33,34,66]. The main outcome is that the ring-fault, formed during the collapse of the caldera triggered by the eruption of the Campanian Ignimbrite, was re-activated during the collapse of the caldera and enhanced by the eruption of the Neapolitan Yellow Tuff [33,34]. A differential volcano-tectonic subsidence has been observed during these geological processes. Five main regional unconformities (U2, U3, U4, U5, and U6) have been distinguished in the Gulf of Pozzuoli based on an updated seismo-stratigraphic framework of the area [34], which is summarized in Figure 2. These unconformities have been correlated with significant geological processes, including volcanic eruptions (Neapolitan Yellow Tuff), temporary and ultimate emersion of marine terraces (La Starza), and coastal subsidence ( Figure 2) [34]. The unconformity U2 is located at the top of the seismo-stratigraphic unit M1 and is a transgressive surface, whose age ranges between 22 ky B.P. and 15 ky B.P. Its formation is genetically related to erosional phases that accompanied the lowstand and early stages of sea level rise ( Figure 2) [34].
The unconformities U3 and U4, whose regional extent is notable, were developed when the rate of tectonic uplift controlled by the caldera resurgence surpassed the rate of sea level rise ( Figure 2) [34]. As a result, the "La Starza" marine terrace emerged (Gulf of Pozzuoli), and a resurgent dome was formed in the caldera depression. While the unconformity U3 has been correlated with the temporary emersion of the "La Starza" marine terrace, the unconformity U4 has been correlated with the ultimate emersion of the same terrace ( Figure 2) [34] significant geological processes, including volcanic eruptions (Neapolitan Yellow Tuff), temporary and ultimate emersion of marine terraces (La Starza), and coastal subsidence ( Figure 2) [34]. The unconformity U2 is located at the top of the seismo-stratigraphic unit M1 and is a transgressive surface, whose age ranges between 22 ky B.P. and 15 ky B.P. Its formation is genetically related to erosional phases that accompanied the lowstand and early stages of sea level rise ( Figure 2) [34]. The unconformity U5 is located at the top of the volcanic sequences V2 and V4, represented by volcaniclastic deposits genetically related to the volcanic eruptions of the Nisida Island and of the Capo Miseno volcanic edifice ( Figure 2) [34]. The age of these deposits ranges between 4.8 and 3.7 ky, although the exact age is strongly debated in the scientific literature [66][67][68][69].
The unconformity U6 is a local erosional unconformity, separating the two seismic units, the seismic unit M4.2a and the seismic unit M4.2b, and was formed during several phases of tectonic uplift preceding the coastal subsidence ( Figure 2) [34]. The onset of the coastal subsidence was dated back at about 2 ky B.P., when the progradational unit M4.2a, coupled with the former unconformity U6, started to subside, thus creating the accommodation space, then available for the deposition of the wedged-shaped seismic unit M4.2b ( Figure 2) [34].
The Gulf of Naples has been deeply studied with respect to its oceanographic setting [70][71][72][73][74]. This area is distinguished by peculiar topographic aspects, which influence the dynamics of the wind and of the sea. The Vesuvius volcano and the hills surrounding Naples City, protecting the basin from the winds coming from north-east, which dominate in the winter, create currents triggering rapid processes of transport from the coastline towards the basin [75][76][77][78][79].
The Gulf of Naples is characterized by the occurrence of two main water masses, which are typical of the central-southern Tyrrhenian sea [75]. The first one is the Atlantic Water, which, being modified along its pathway (MAV), reaches water depths ranging between 50 m and 100 m. Its salinity is 37.5%, and the temperature follows the seasonal trend, decreasing accordingly with the increase in water depths. During the winter, the temperature assumes a quite constant value of about 14 • , due to the whole remixing of the water column. In the deepest stations, located in correspondence with the Bocca Grande high, the Levantine Intermediate Water (LIW) occurs. Moreover, other types of water masses occur, depending on the seasonal trend. The winter mixing controls the formation of the Tyrrhenian Intermediate Water (TIW). After the summer heating, being located at water depths of about 75, the TIW turns into the TSW (Tyrrhenian Surface Water), which is characterized by a temperature of about 26 • C and by a salinity of 38%. Being controlled by the input of freshwaters supplied by the Sarno and Volturno rivers in the Gulf of Naples, another water mass, similar to the TSW, but warmer and less salty, can be identified [75].
Five sketches (Figure 3) of oceanographic circulation have been developed, depending on the seasonal and dynamic conditions [79]. The first two sketches ( Figure 3) take into account the Tyrrhenian offshore currents, which, being oriented towards the north-west, control the entrance of waters in the Gulf of Naples from Bocca Piccola (Capri) and create a regional meander, moving counter-clockwise. After the reversal of the Tyrrhenian current, a cyclonic gyre separates the offshore areas from the coastal ones. In the other three sketches (Figure 3), the wind is the main control factor. During the winter, the north-eastern winds prevail and the Vesuvius volcano, due to its topographic location, channels the wind and produces a jet directed towards the basin. In the presence of south-western winds, frequently during the onset of air-pressure lows, cyclonic and anti-cyclonic structures form in the Gulf of Naples, triggering the formation of coastal currents and the retaining of the surface waters in the littoral area. During the summer, diurnal alternating land and sea breezes occur in the coastal areas, controlling the clockwise rotation of the field of the surface currents in the Gulf. Strong along-shore currents are controlled by the sea breeze, while the land breeze, mainly occurring during the night, controls relatively weaker currents.

Materials and Methods
The methods used for the interpretation of the seismic profiles include the seismic stratigraphy, the sequence stratigraphy, and the volcanic seismic stratigraphy. The seismic stratigraphy was born as a key methodology in the geological interpretation of seismic profiles at the end of the 70s [80][81][82]. The key idea of the seismo-stratigraphic approach is that the horizons are comparable to the strata plans and their geometry reflects the geometry of deposition. Due to this concept, seismic stratigraphy conjugates geological methods with geophysical ones [83]. Catuneanu (2002) [84] has revised the concepts of sequence stratigraphy, regarding the base level changes, transgressions and regressions, the stratigraphic surfaces (types of stratal terminations and sequence stratigraphic surfaces) and their time attributes, the systems tracts, the models of the depositional sequences and the importance of several control factors. The LST is limited in its lower part by a subaerial unconformity and its correlative conformity in a marine environment and in its upper part by the maximum regressive surface. The corresponding deposits may include incised valley fills, amalgamated fluvial channels, and shelf margin deltas, whose preservation potential is controlled by the consequent ravinement surface [84][85][86][87]. The TST is limited in its lower part by the maximum regressive surface and in its upper part by the maximum flooding surface. The corresponding deposits include coastal marine deposits and correlative condensed sections and healing-phase deposits, deposited in the shoreface, more than backstepping beach deposits [84]. The HST is limited in its lower part by the MFS (maximum flooding surface) and by the regressive surface of marine erosion in its upper part. The corresponding deposits include, in their marine part, progradational and aggradational sequences [84,85], whose preservation potential is controlled by the interaction between erosional processes, both in subaerial and submarine environment [88,89]. The FST includes the deposits corresponding to base-level fall in the marine part, contemporaneously with the evolution of the subaerial unconformity towards the land, represented by progradational strata [90,91].
Another important research method is represented by the volcanic seismic stratigraphy, which is based on the interpretation of the volcano-seismic facies [92], where key volcanic facies were defined based on the volcanic facies' architecture. Planke et al. [92] have pointed out the occurrence of six volcanic seismic facies, including the landward flows, the lava delta, the inner flows, the Inner Seaward Dipping Reflectors (Inner SDR), the Outer High, and the Outer SDR. These facies have been put in the geological framework of the five-stage evolution of a continental margin, including the explosive volcanism in a broad basin, the subaerial effusive volcanism, forming Gilbert-type lava deltas, the subaerial effusive volcanism filling a narrow rift basin and, finally, the deep marine sheet flow or pillow-basalt volcanism [92]. Additional studies have been carried out, which link good onshore volcanic facies outcrops with offshore interpretation to build up a more robust interpretation of the seismic data, where no well data exist to confirm the interpretations [93,94]. More recently, the stratigraphic architecture of the volcanic seismic facies at rifted margins was summarized by clarifying the volcanic key facies, including lava flows, prograding hyaloclastites and other important facies, which can be recognized on the seismic sections in volcanic settings [95]. The concepts of the volcanic seismic stratigraphy were recently applied offshore northern Ischia, where a correlation between volcanic seismic units and outcropping volcanic facies has been established [96].
The Sparker seismic profiles were acquired during the technical activities aimed at the realization of the CARG (CARta Geologica) project at the 1:10,000 scale by using a Sparker Multi-tip seismic source [36]. During the seismic survey, several seismic profiles perpendicular to the shoreline were recorded in the southern, eastern, and offshore northern Ischia [36,97]. The seismic acquisition was realized through an EG&G Trigger Capacitor Bank mod. 231 and PSU mod. 232, and a SAM96 Sparker seismic source. The range of frequencies included in the seismic signal extends from 100 Hz to 3000 Hz. Moreover, the acquisition system D-Seismic was used [36]. The recorded data were stored in a DDA (Digital Data Array) format on a hard disk. In a playback mode, the D-Seismic software generated the images of the seismic profiles in a jpg or bmp format. The seismic profiles, in a jpg format, were imported in the Coreldraw Graphic Suite in order to perform the geological interpretation. In this paper, particular attention was paid to the detailed interpretation of seismic lines located in the northern and eastern offshore of the island that are significant for the identification of the bioclastic deposits and the interpretation of a regional seismic line located in the Procida Channel, a morphological threshold located between the Ischia and Procida islands. The location of the seismic profiles analyzed in this paper is shown in Figure 1.
At Ischia and Procida-Pozzuoli, the sea bottom samples analyzed in this paper were collected during several oceanographic cruises carried out on the R/V Urania (National Research Council, Italy) in the frame of the CARG project of the Campania region (geological map n. 464 "Island of Ischia" scale 1:10,000) [36]; (geological map n. 465 "Island of Procida", scale 1:10,000; scale 1:50,000) [35]. The Ischia samples were collected during the oceanographic cruises GMS02_01 (2002) Figure 4 shows the DEM of Ischia-Procida, superimposed with the location of sea bottom samples. The sedimentological data have been plotted in ternary diagrams. While we can refer to the previous study for the Ischia samples and their elaboration ( Figure 5) [37], at Procida-Pozzuoli, the samples were plotted by dividing them into two groups, from sample B1A2 to sample B1054 (Figure 6,c) and from sample BXMIS1 to sample PM27 (Figure 6b,d).   [37], at Procida-Pozzuoli, the samples have been plotted by dividing them into two groups, from sample B1A2 to sample B1054 (Figure 6a,c) and from sample BXMIS1 to sample PM27 (Figure 6b,d).

Sedimentological Data
Gravelly muds have been found in the northern Ischia offshore (Ischia, Ischia harbor and Casamicciola). The gravelly fraction is composed of volcaniclasts (pumice and tuff) and bioclasts (lamellibranchs, red coralline algae and single corals). These facies can be interpreted as the biocenosis of the "Détritique Cotier". Moreover, offshore northern Ischia, middle-grained bioclastic silty sands with Posidonia Oceanica sea-grass leaves and silts with bioclastic fragments (M.te Vico) and bioclastic sands immersed in a silty matrix (Punta La Scrofa) have been found. The bioclastic fragments consist of gastropods, bivalves, calcareous algae, and bryozoans.

Sedimentological Data
Gravelly muds have been found in the northern Ischia offshore (Ischia, Ischia harbor and Casamicciola). The gravelly fraction is composed of volcaniclasts (pumice and tuff) and bioclasts (lamellibranchs, red coralline algae and single corals). These facies can be interpreted as the biocenosis of the "Détritique Cotier". Moreover, offshore northern Ischia, middle-grained bioclastic silty sands with Posidonia Oceanica sea-grass leaves and silts with bioclastic fragments (M.te Vico) and bioclastic sands immersed in a silty matrix (Punta La Scrofa) have been found. The bioclastic fragments consist of gastropods, bivalves, calcareous algae, and bryozoans.
Offshore southwestern Ischia (Citara), volcanic rocky outcrops with superficial algal incrustations colonized by red coralline algae have been found at water depths of 40 m. Here, coarse-grained volcanic sands have also been found. In the same area, offshore the Punta del Chiarito promontory, the bioclastic deposits are composed of medium-to-coarsegrained sands and gravelly sands with red coralline algae, gastropods, and lamellibranch fragments. These deposits disappear at increasing water depths (about 168 m), where they are replaced by sandy silts with reworked Posidonia Oceanica.
On the Ischia Bank, a relict volcanic edifice located offshore south-eastern Ischia, at water depths of 24 m sands, with gravelly sands with a silty matrix have been found, overlying coarse-grained bioclastic sands. Lithics are admixed with bioclastic fragments, represented by bivalves, gastropods, and echinids.
At water depths of 30 m, silty sands have been found, represented by coarse-grained sands with a silty matrix, including small bivalves and small lithic fragments. Proceeding towards increasing water depths (74 m) sandy silts occur, mostly represented by clays with a dark green sandy fraction, including volcaniclasts (volcanic lithics and pumices) and bioclasts (bivalves, rhizomes of Posidonia Oceanica, gastropods, and echinoids). On the parasitic vent of the Ischia Bank, sands and gravelly sands occur, represented by bioclastic fragments immersed in a silty matrix, with glass fragments, branched bryozoans and red algae.
Offshore of Procida, sands and gravelly sands have been found, in particular sands with a silty matrix, which overlie coarse-grained bioclastic sands. The bioclastic component is composed of fragments of bivalves, small lamellibranch, red algae, small gastropods, rhizomes of Posidonia Oceanica and branched bryozoans. In the same area, silty sands have been found, with a very fine-grained sandy fraction, which consists of bioclasts, pumiceous clasts, and lithics. In the first meters below the seafloor, remains of Posidonia Oceanica have been found, constituting masses of rhizomes. Moreover, coarse-grained bioclastic sands occur, whose bioclastic component is mainly composed of red algae and branched bryozoans.
The grain size of deposits cropping out at the sea bottom was determined by means of the sedimentological analysis carried out at the CNR ISMAR sedimentological laboratory for the Ischia samples ( Figure 5) and the Procida-Pozzuoli samples ( Figure 6). The location of samples is reported in Figure 4 (Ischia-Procida).
At Ischia, the recognized grain sizes are gravelly sands, sands, silty sands, muddy sands, sandy silts, and silts. Ternary plots were constructed ( Figure 5), respectively taking into account clay-sand-silt and gravel-sand-silt. The samples were split into the samples collected in the oceanographic cruise GMS02_01 (Figure 5a,b) and in the oceanographic cruise GMS06_03 (Figure 5c,d).
The samples were split into the samples collected in the oceanographic cruise GMS02_01 ( Figure 5; inset a; inset b) and in the oceanographic cruise GMS06_03 (Figure 5; inset c; inset d).
We can refer to the paper of Aiello, 2021 [37] for the description of the sedimentological data on Ischia, which have been here reported with the aim to be compared with the sedimentological data of Procida-Pozzuoli.
At Procida-Pozzuoli, the grain sizes recognized at the sea bottom based on grainsize-analysis include sands, silty sands, muddy sands, silts, sandy muds, and muds. The percentage of the gravelly grain size is lower than that observed in the Ischia samples [37].
Ternary plots were assembled (Figure 6), taking into account clay-sand-silt and gravelsand-silt, respectively. The samples were sub-divided into two groups, from sample 1A2 to sample 1054 (Figure 6a,c) and from sample BXMIS1 to sample PM27 (Figure 6b,d). The ternary plot of Figure 6a considers the variables: clay, sand, and silt. This plot has shown that the main grain sizes are the clayey silts, the sandy silts, the silty sands, and the sands. The ternary diagram of Figure 6b considers the variables: clay, sand and silt. This plot has shown that the main grain sizes consist of sands and subordinately of silty sands. The ternary plot of Figure 6c considers the variables: gravel, sand and silt. The plotted diagram suggests that the main grain sizes are the sands, the silts, and the silty sands and only subordinately, the gravelly sands. The ternary diagram of Figure 6d was constructed considering the variables: gravel, sand, and silt. This plot has shown that the prevalent grain sizes are represented by sands and silty sands, and, only subordinately, by gravelly sands.
Based on the description of the samples and on the sedimentological analysis, three facies associations have been distinguished.

Facies 1-Sands and Gravelly Sands
The first facies association is composed of sands and gravelly sands, with a significant bioclastic component. It has been found in the infralittoral zone and in the uppermost part of the circalittoral zone, at average water depths ranging between 20 and 40 m (southwestern Ischia, Ischia Bank, Ischia harbor, Casamicciola, Procida, Miseno Bank).

Facies 2-Silty Sands
The second facies association is represented by silty sands and is composed of finegrained sands with a bioclastic component. It has been detected in the infralittoral zone, at average water depths ranging between 40 and 70 m (Citara, Punta Imperatore, Ischia Bank, Ischia harbor, Casamicciola, Monte Vico).

Facies 3-Sandy Silts and Silts
The third facies association is represented by sandy silts and silts and is composed of silt with a fine-grained sandy fraction, both volcaniclastic and bioclastic. It has been detected in the infralittoral zone, at average water depths ranging between 70 m and 145 m (Ischia Bank at the base of the volcanic edifice, Citara, Punta del Chiarito, Procida, Forio Bank at the base of the volcanic edifice, Cava Grado, Punta Caruso).

Seismo-Stratigraphic Data
A new and detailed interpretation of the previously interpreted Sparker seismic profiles available around Ischia was carried out, focusing on the Ischia Channel and offshore Casamicciola (Figures 7 and 8) [36]. The aim of this interpretation is the identification of the bioclastic deposits and the definition of their stratigraphic relationships with the volcanic and the sedimentary units recognized offshore of Ischia [77].
Offshore Casamicciola, the rhodolith deposits are probably inter-layered in a wide seismic unit, characterized by parallel and continuous seismic reflectors and interpreted as highstand deposits (Figure 7) [97]. The Late Quaternary depositional sequence consists of both highstand deposits and forced regression prograding wedges (Figure 7). In this area, the transgressive system tract has not been recognized; perhaps there are no hints that the rhodolith deposits could have been deposited during the transgressive system tract. Previous seismo-stratigraphic studies have defined the stratigraphic architecture of offshore northern Ischia (Casamicciola), where progradational deposits, ascribed to forced regression, have been recognized [97].
Moreover, in the Gulf of Naples, the transgressive system tract is very reduced, or completely lacking, as shown by previous seismo-stratigraphic data [97,98]. For this reason, we cannot interpret the rhodolith deposits of the Gulf of Naples as deposited during transgressive phases, as highlighted for the Styrian Basin [99]. Friebe [99] considered that Middle Miocene rhodolith platforms of the Styrian Basin (Austria) grew during transgressive phases and were then onlapped by siliciclastic deposits during the early highstand phase. In the Styrian Basin, shallowing-upward siliciclastic parasequences have been recognized, terminating with a bank of rhodolith limestones, which formed during the late highstand system tract [99]. However, the age of deposits (Middle Miocene) and the stratigraphic context appear to be quite different when compared to the Gulf of Naples [100] and the rhodolith deposits described in this paper, which seem to be very recent, as they were found in sea bottom samples.
Several volcanic seismic units have been recognized (V1, V2, V3 in Figure 7). The V1 unit is characterized by acoustically transparent seismic facies and has been interpreted as an undetermined volcanic acoustic basement, topped by an erosional unconformity ( Figure 7).
The V2 unit is characterized by acoustically transparent seismic facies and by a mounded-shaped external geometry and corresponds to buried volcanic edifices, which appear to be partly tectonically controlled by normal faults, having a little vertical throw (Figure 7).
The V3 unit, interpreted as an undetermined pyroclastic unit, locally filling depressions, is characterized by discontinuous seismic reflectors and forms a mounded-shaped pyroclastic edifice (Figure 7). Both the M1 unit and the M2 unit are characterized by parallel and sub-parallel seismic reflectors (Figure 7).
Similarly to what was observed offshore northern Ischia, offshore eastern Ischia, the rhodolith deposits can be located in the upper part of the seismo-stratigraphic unit interpreted as the highstand deposits ( Figure 8) [97]. The geological interpretation of the seismic profile has shown a main volcanic feature ("Il Pertuso" volcano), extending from the Ischia Channel towards the Procida continental shelf. This feature, whose seismic facies is acoustically transparent, is probably composed of pyroclastic deposits and lavas, which are genetically related to the relict volcanic edifices of the Ischia Channel. Other significant seismo-stratigraphic units ( Figure 8) include a unit, probably pyroclastic in nature and a progradational unit, representing a relict prograding wedge, probably Late Pleistocene in age. The location offshore Procida of the pyroclastic unit suggests a possible correlation with the Solchiaro Tuff (18 ky B.P.).
The geologic interpretation of a Sparker line located in the Procida Channel was carried out in order to show the stratigraphic architecture of this area (Figure 9). Four main seismic units were recognized based on the geological interpretation. The lowest ones are represented by the pyroclastic sediments deposited before the Campanian Ignimbrite (39 ky B.P.; Figure 9). These units are overlain by the Campanian Ignimbrite (CI in Figure 9). The CI is overlain, in turn, by buried volcanic domes, whose emplacement seems to be younger than 39 ky B.P. The youngest seismic unit is represented by the pyroclastic sediments deposited after the Campanian Ignimbrite (CI; Figure 9).
Based on seismo-stratigraphic interpretation, there are no hints of the occurrence of bioclastic deposits in the Procida Channel since the pyroclastic units and the volcanic domes seem to prevail (Figure 9).
Quaternary 2021, 4, x FOR PEER REVIEW 17 pyroclastic edifice (Figure 7). Both the M1 unit and the M2 unit are characterized by allel and sub-parallel seismic reflectors (Figure 7). Similarly to what was observed offshore northern Ischia, offshore eastern Ischia rhodolith deposits can be located in the upper part of the seismo-stratigraphic unit in preted as the highstand deposits (Figure 8) [97]. The geological interpretation of the mic profile has shown a main volcanic feature ("Il Pertuso" volcano), extending from Ischia Channel towards the Procida continental shelf. This feature, whose seismic faci acoustically transparent, is probably composed of pyroclastic deposits and lavas, w are genetically related to the relict volcanic edifices of the Ischia Channel. Other signifi seismo-stratigraphic units (Figure 8) include a unit, probably pyroclastic in nature a progradational unit, representing a relict prograding wedge, probably Late Pleistocen  The geologic interpretation of a Sparker line located in the Procida Channel was c ried out in order to show the stratigraphic architecture of this area (Figure 9). Four m seismic units were recognized based on the geological interpretation. The lowest ones represented by the pyroclastic sediments deposited before the Campanian Ignimbrite ky B.P.; Figure 9). These units are overlain by the Campanian Ignimbrite (CI in Figure  The CI is overlain, in turn, by buried volcanic domes, whose emplacement seems to The Campanian Ignimbrite is an important pyroclastic deposit, genetically related to a collapsed caldera [101][102][103][104][105][106][107]. In this caldera, volcanic activity of little entity developed up to 1538 B.P. with the eruption of Monte Nuovo volcanic edifice. Based on seismo-stratigraphic interpretation, there are no hints of the occurrence of bioclastic deposits in the Procida Channel since the pyroclastic units and the volcanic domes seem to prevail (Figure 9).

Discussion
The sedimentological and seismo-stratigraphic results related to offshore Ischia and Procida-Pozzuoli were compared with previous studies on Mediterranean rhodolith layers at a regional scale [1][2][3][4]. At a local scale, these results have been compared with previous data on rhodolith deposits in the Southern Tyrrhenian sea (Pontine islands) [5], in the Gulf of Naples, and offshore of the Cilento Promontory.
The rhodolites were discovered in the Mediterranean at several sites and range in water depth between 30 m and 75 m [4]. They have been found in most of the marginal areas of the Mediterranean Sea, where they have a wide spatial distribution.

Discussion
The sedimentological and seismo-stratigraphic results related to offshore Ischia and Procida-Pozzuoli were compared with previous studies on Mediterranean rhodolith layers at a regional scale [1][2][3][4]. At a local scale, these results have been compared with previous data on rhodolith deposits in the Southern Tyrrhenian sea (Pontine islands) [5], in the Gulf of Naples, and offshore of the Cilento Promontory.
The rhodolites were discovered in the Mediterranean at several sites and range in water depth between 30 m and 75 m [4]. They have been found in most of the marginal areas of the Mediterranean Sea, where they have a wide spatial distribution.
The zonation of Mediterranean benthic assemblages established by Péres and Picard (1964) [108] has allowed for a better understanding of the lithologies and paleoenvironmental interpretations of the Miocene rhodolithic limestones [6]. Significant facies associations are represented by the "Détritique du large" (offshore deposits), by the rhodolith deposits (facies à pralines), by the "Détritique Cotier" (inshore deposits), and by the maërl deposits [6,108].
In the Pontine archipelago (Ponza and Palmarola), the carbonate sedimentation has been mainly controlled by red coralline algae and is concentrated in a channel located between the two islands [5]. The obtained results match the biocenosis of the "Détritique Cotier". The red coralline algae have controlled the carbonate deposition, whose deposits range in depth between 40 m and 70 m. This bathymetric interval is consistent with the bathymetric interval of the rhodoliths offshore Ischia. Another significant point of discussion is represented by the maërl facies and by its occurrence in Naples Bay, mainly if compared with offshore Cilento [109]. The maërl deposit displays the entire and/or fragmented thalli of calcareous algae, which amass in the depressions of the bedrock.
The term "maërl" was proposed by Crouan and Crouan [110] as a Breton word for unattached, branched corallines, alive or dead and loose lying, which often occurr in wide deposits or algal gravels found off the north-western coast of France (maërl beds or banks). The rhodolith layers frequently include the maërl facies.
In the Gulf of Naples, maërl deposits are rare and have only been detected in a few outcrops of the Gulf of Pozzuoli (Miseno Bank) [14] or offshore northwestern Ischia [111,112]; instead, they widely crop out at the sea bottom offshore of the Cilento promontory [109]. The maërl deposits predominate on the depositional terraces at water depths ranging from 42 m to 52 m [109].
The preferential distribution of the depositional terraces at these water depths has probably been controlled by background currents, generated by the local circulation, preventing the deposition of siliciclastic deposits. The sampling data have shown that the maërl deposits are composed of bioclastic gravels, which mainly consist of alive or dead thalli of red coralline algae [109]. In addition, there are sandy and sandy-muddy facies in transition to the maërl facies. This transition is similar to what is observed offshore Ischia, where sandy silts and silty sands are often associated with the bioclastic deposits, more properly gravelly sands ("Détritique Cotier") [108].
Offshore Procida-Pozzuoli, the sedimentological data of the analyzed sea-bottom samples show a lack of carbonate sedimentation, while mainly siliciclastic deposits prevail. In the Procida Channel, the sampling data have outlined that sandy silts and silty sands prevail. There exists a possibility that a hypothetical maërl facies, which is frequently associated with this kind of lithologies, may occur in this area, but it was not detected by the collected sea bottom samples. Further sampling data are perhaps required in order to support this working hypothesis.
Some bioclastic deposits have been mapped offshore of Procida island and have been outlined by the explanatory notes to the geological map [35], namely, the psgb facies [35]. The psgb facies is composed of bioclastic gravels and sands and crops out at the sea bottom at water depths above 30 m, in correspondence with the relict volcanic edifices of the Ischia Channel ( Figure 8) and with protected areas surrounding the island of Procida [35].
Offshore Ischia, towards offshore Procida-Pozzuoli, a transition from bioclastic to siliciclastic deposits has been observed. Moreover, the geological interpretation of seismic profiles located in the Procida Channel ( Figure 9) has shown that pyroclastic units and volcanic domes prevail in this marine sector. Based on the seismo-stratigraphic data, the most important geological event controlling the stratigraphic setting of the Procida Channel is the emplacement of the Campanian Ignimbrite (CI seismic units in Figure 9). However, important pyroclastic units were emplaced before and after the CI (pre-and post-CI seismic units in Figure 9).
Offshore Ischia-Procida, mixed carbonate-siliciclastic systems occur. These systems have been recently examined from a sequence stratigraphic point of view [110]. Zecchin and Catuneanu [110] have highlighted that in this kind of system, significant vertical shifts of the facies occur in the system tracts. A "reciprocal sedimentation pattern" was recognized, with transgressive/highstand carbonates and forced regression/lowstand siliciclastic deposits. A mixed depositional sequence is composed of a lower part, mainly consisting of carbonates, which can be interpreted as the Transgressive System Tract or the Transgressive System Tract plus the Highstand SystemTract, and of an upper part, mainly consisting of siliciclastic deposits, which can be interpreted as the Forced Regression System Tract and the Lowstand System Tract or as the Highstand System Tract plus the FSST/LST [110]. In this kind of depositional sequence, the carbonate deposition is typified by carbonate platform/ramp facies, which was deposited as controlled by a reduced siliciclastic supply during the transgression and the highstand, while the siliciclastic deposition is typified by regressive shelf and shoreface deposits, or by a deltaic wedge [110].
Mixed carbonate/siliciclastic depositional systems have been recognized and analyzed in the central Tyrrhenian sea [111,112]. Brandano and Civitelli [111] recognized six sedimentary facies on the continental shelf of the Pontine islands based on component analysis, grain size percentage, sorting, carbonate content, and authigenic mineralization rate. These facies include the planktonic silt, the skeletal sands (poorly sorted skeletal sand, silty skeletal sands), the foraminifer-crinoid sands, the maërl facies, the siliciclastic sands (terrigenous sands, siliciclastic skeletal sands) and the relict sands. The maërl facies represents the depositional environment with the highest sedimentation rates [111]. These facies may be compared with the facies associations distinguished in this paper and are represented by the facies 1 (sands and gravelly sands), facies 2 (silty sands) and facies 3 (sandy silts and silts). However, further sedimentological and paleontological studies, including the component analysis, are required in order to perform a detailed facies interpretation, which are planned in the future.
As a general rule, this kind of stratigraphic architecture for the mixed carbonatesiliciclastic depositional systems, as delineated by Zecchin and Catuneanu [110], fits well both with the seismo-stratigraphic data analyzed in this paper (Figures 7 and 8) and with previous data on the seismo-stratigraphic architecture in the Gulf of Naples [97,98,113]. This stratigraphic architecture is characterized by thick forced regressive and lowstand system tracts, very reduced or lacking transgressive system tracts, and well-developed highstand system tracts [113].

Conclusions
In the Gulf of Naples, mixed carbonate-siliciclastic systems were recognized based on the analyzed data. These systems were deeply studied, referring both to their sequence stratigraphic framework [110] and to the control factors triggering the siliciclastic deposition, or alternatively, the carbonate deposition. The siliciclastic systems are mainly controlled by the sediment supply, by the relative sea level fluctuations, by the climate, by the physiography and by the hydrodynamic setting [84,85,90,91]. The carbonate systems are also influenced by the trophic and chemical factors, in addition to the control factors mentioned above in the case of the siliciclastic systems [6,7,9,11,13].
In this paper, three main facies have been recognized, consisting of sands and gravelly sands (facies 1), silty sands (facies 2) and sandy silts and silts (facies 3). These facies appear to be concentrated in specific bathymetric intervals (respectively from 20 to 40 m, from 40 to 70 m, and from 70 to 145 m). The bathymetry appears to be one of the most important control factors acting on the distribution of these facies associations ( Figure 4).
Offshore Ischia island, the carbonate deposits composed of red coralline algae have a contrasting distribution. In particular, the rhodolite deposits have mainly been found on the relict volcanoes, topographic saddles, and roughly topographic areas of the sea bottom, which were established consequently to the emplacement of debris avalanches at the seafloor.
In the Procida Channel, background currents probably acted, so hindering the carbonate sedimentation. Moreover, in this area, pyroclastic units and lava domes prevailed, as a consequence of the main geologic and volcanologic event of the Campanian Ignimbrite eruption (39 ky B.P.). Diurnal alternating land and sea breezes occur in the coastal sectors of the Gulf of Naples (Figure 3). It is probable that this oceanographic condition acted in the Procida Channel, where along-shore currents could have been controlled by the sea breeze, so preventing the formation of bioclastic deposits and influencing mainly a siliciclastic deposition.
In the Gulf of Naples siliciclastic-carbonate mixing occurs in the infralittoral zone. Here, the erosional processes have been active in correspondence with the rocky shorelines, particularly abundant at the Ischia and Procida Islands, which have produced lithoclastic and volcaniclastic components. These components were probably then reworked by the action of the wave currents. The resulting detrital deposits are interlayered with the carbonate deposits, which were preferentially deposited on the relict volcanic edifices and in the bays alternating with the rocky promontories. On the relict volcanic edifices, in particular at the Ischia Bank, the facies relationships between facies 1, located at the top of the volcanic edifice and facies 2, located on the flanks and at the base of the volcanic edifices, can be supposed. While facies 1 deposited as layers located at the top of the relict volcano, facies 2 represents a drape located on the flanks of the volcanic edifice.
Author Contributions: G.A.: writing-original draft preparation, data interpretation, writingreview and editing; M.C.: methodology; software; data curation. All authors have read and agreed to the published version of the manuscript.