Quaternary Evolution of the Lower Calore and Middle Volturno Valleys (Southern Italy)

: The lower Calore and middle Volturno valleys preserve stratigraphical and morphological evidence and tephrostratigraphic markers particularly suitable for reconstructing the long-term geomorphological evolution of the central-southern Apennines. Aim of our study is to identify the main steps of the Quaternary landscape evolution of these valley systems and to improve knowledge about the relationships between fluvial processes and tectonics, volcanic activity, climatic and human influences. To this purpose, we carried out an integrated geomorphological and chrono-stratigraphical analysis of identified fluvial landforms and related deposits, integrated by 230 Th/ 234 U datings on travertines from the Telese Plain area. The study highlighted in particular: (1) fluvial sedimentation started in the Middle Pleistocene (~650 ka) within valleys that originated in the lower Pleistocene under the control of high-angle faults; (2) extensional tectonics acted during the Middle and Upper Pleistocene, driving the formation of the oldest fluvial terraces and alluvial fans, and persisted beyond the emplacement of the Campanian Ignimbrite pyroclastic deposits (~39 ka); and (3) from the late Upper Pleistocene onwards (<15 ka), the role of tectonics appears negligible, while climatic changes played a key role in the formation of three orders of valley floor terraces and the youngest alluvial fans.


Introduction
Physical landscapes and related landforms result from complex sequences and overlaps of geological and geomorphological processes that have contributed to its long-term shaping, therefore providing valuable information on past morphodynamics and geological events.
Considering in particular the context of a young and active mountain chain, such as the Apennines [1,2], long-term landscape evolution appears complex and diverse in space, mainly responding to tectonic events that have implicated diverse regimes and spatio-temporal scales for different chain sectors (e.g., [2,3]), and environmental changes due to volcanic events and/or climate oscillations and related sea level changes (e.g., [4,5] and references therein).
Analyzing such a long-term landscape can provide valid data about the Quaternary geomorphological evolution of a certain area and the roles played by tectonics and climate changes (e.g., [6,7] and references therein). This is particularly true for terraced fluvial landscapes that can provide broad evidence about the relations of surface processes, tectonics, and climate influences [8][9][10][11]. Fluvial responses to phases of regional uplift (especially valley downcutting and widening) (e.g., [10]) and to climate fluctuations (especially warm-cold and cold-warm transitions) resulting in alternated phases of valley filling and incision [5,12,13] are scientific issues of particular interest and high complexity (e.g., [14] and references therein).
Starting from these considerations, to improve our knowledge on the Quaternary evolution of the central-southern Apennines, we searched for a study area that could offer the following advantages: (i) being located in a sector of the Apennine chain that is part of a wider areal context for which essential data about geological history and tectonics are already available; (ii) being characterized by a rich, diversified fluvial landscape; iii) being characterized by stratigraphical and geomorphological evidence providing morpho-chronological constraints suitable for landscape reconstruction. Once examined several regional contexts and relative data available, we focused our attention on the Volturno River system (Figure 1) that occupies a central position in the central-southern Apennines. As highlighted by several studies, the Volturno valley preserves rich evidence about the morpho-evolutive events that have affected this sector of the Apennine chain during the Quaternary [10,15,16]. It represents a key area that hosts stratigraphical and morphological evidence along with several tephrostratigraphic markers suitable for investigating the long-term evolution of this active chain sector.
The upper Volturno River valley has been subject of various studies [15,[17][18][19][20][21][22][23]] that have provided a basic framework of its tectonic and sedimentary Quaternary evolution. Studies on the middle Volturno River valley mainly focused on the southern Matese mountain front and the adjoining Alife basin [16,[24][25][26][27][28] (Figure 1b) and different reconstructions of the Quaternary tectonic framework are proposed [16,26], but no specific studies on fluvial terraces have been performed. Numerous studies concern especially the lower Calore River valley [15,[29][30][31][32][33][34], the main left tributary of Volturno River, mainly aimed at the reconstruction of the morphotectonic framework and the identification of seismogenic faults. Different models for the morphotectonic and paleoenvironmental evolution have been proposed highlighting a number of open questions regarding mostly the origin and ages of the fluvial terraces that characterize this valley. However, except for the paper of Massa and Zuppetta [35], which focuses on the recent tectonic activity in the confluence area of Calore and Volturno rivers, studies about the evolution of the Calore and Volturno rivers valleys and related chronological aspects are completely lacking. Especially the relations between morphogenetic events, climate and tectonic drivers and the largest Quaternary volcanic explosive event in the Mediterranean region [36] that has given origin to the Campanian Ignimbrite (CI) pyroclastic flow deposits (39 ka [37]), have not been fully addressed and clarified.
Based on already acquired data and knowledge gaps, briefly outlined above, we selected the middle Volturno River valley including the lower Calore valley (Figure 1a) as our study area.
Given the purpose to reconstruct the Quaternary evolution of the study area, major specific goals of this study are (i) to reconstruct and characterize in detail the stratigraphical sequences of fluvial terraces and alluvial fans; (ii) to add morpho-chronological constraints suitable for framing the geomorphologic evolution of the study area; and (iii) to better understand the possible roles of tectonics, climate, and volcanic events in the development of the fluvial landscape.

Study Area
The Volturno River, with a length of ~175 km and a catchment area of 5500 km 2 , is the largest river in southern Italy. The Calore River, with a length of 208 km and a catchment area of 3058 km 2 , is its main left tributary (Figure 1a).
The upper and middle stretches of Volturno River flow through a wide valley that crosses the axial zone of the central-southern Apennine chain. The lower segment of the Volturno River crosses the northern sector of the Campania coastal alluvial plain which is part of the broad peri-Tyrrhenian Piana Campania Graben [38,39]. The sector of the southern Apennine chain that includes the middle Volturno River valley is part of a Neogene fold and thrust system generated by complex tectonic deformations started during the Miocene [40,41]. The tectonic units consist of carbonate platform and basinal Meso-Cenozoic sedimentary successions [42,43] (Figure 1b), which were progressively incorporated in the Apennine accretionary wedge [44]. This sector of the Apennine chain emerged since the Pliocene [45]. Starting from the Lower Pleistocene, it was affected by NW-SE extensional tectonics that generated high-angle fault systems mainly composed of NW-SE trending strike-slip and NE-SW trending normal faults, responsible for the onset of Quaternary fluvial and slope sedimentation [46]. In the Middle Pleistocene, a change in the direction of extension, from NW-SE to NE-SW, caused the reactivation of previous structures, especially those oriented in the NW-SE direction [32,[47][48][49].
The pre-Quaternary bedrock mainly consists of Mesozoic-Cenozoic limestones and dolostones belonging to carbonate platform units, which are overlain by the terrigenous successions of the Upper Miocene Molise and Caiazzo Flysch [50] (Figure 1b). The latter are covered by basinal successions of the Sannio units that crop out in the eastern sector of the study area.
Quaternary deposits are widespread in both the middle Volturno and the lower Calore valleys (Figure 1b). They consist of fluvial, alluvial fan, eluvial-colluvial, slope, travertine, and volcanic deposits. The latter are prevailingly represented by the CI pyroclastic flow deposits (for details, see Section 2.1).
For this study, the following three valley sectors A, B, and C (Figure 1b), were considered.
Sector A coincides with the portion of the lower Calore River valley that extends between the villages of Ponte and Tuoro S. Antuono (Figure 1b). In this sector, the Calore River has a meandering planform morphology and crosses a structurally controlled, E-W elongated valley depression, limited to the north by the gentle marginal slopes of Matese Massif, mainly made of terrigenous and basinal successions of the Molise Flysch and Sannio Unit, and to the south by the steep Camposauro carbonate mountain front ( Figure  1b).
Sector B includes the portion of the lower Calore River valley located between Tuoro S. Antuono and the confluence with the Volturno River, and the stretch of the middle Volturno River valley located between the confluences of the Titerno Stream and the Calore River (Figure 1b). In this sector, the Calore River maintains a meandering planform morphology. The Volturno River has a length of ~10 km and a sinuous to meandering channel [51] and crosses a NW-SE elongated structurally controlled valley depression, mainly underlain by Miocene flysch deposits.
Finally, sector C consists in the stretch of the middle Volturno valley that extends between the Calore confluence and the Triflisco dam ( Figure 1b). Upstream the confluence of Isclero Stream, the Volturno has a length of 9 km, a sinuous morphology, and follows a prevailing NE-SW direction, crossing a valley mainly underlain by Miocene Flysch deposits. Downstream the Isclero confluence, the Volturno channel has a length of ~17 km and a prevailing meandering morphology; it maintains a prevailing NE-SW direction and is mainly underlain by Mesozoic carbonate rocks.

Literature Overview of Quaternary Deposits
To summarize pre-existing data and highlight related knowledge gaps, we performed a detailed literature review of Quaternary deposits.
The oldest outcropping Quaternary deposits are preserved in the lower Calore valley (sectors A and B in Figure 1b). These are mainly terraced alluvial deposits, widespread on the northern valley slope, and slope and alluvial fan deposits prevailing on the southern slope. Among the latter, the oldest deposits are the tectonized and well-cemented Lower Pleistocene Laiano Synthem slope breccias [45], which are widespread along the northern slope of Camposauro Mt (Figure 1b). Alluvial fan deposits were found widely in the Calore valley and several generations have been distinguished by different authors [15,[31][32][33][34][35], but without agreeing on the number of generations and related ages. In particular, Massa et al. [31], Di Bucci et al. [32,33] and Massa and Zuppetta [35] have distinguished two generations referred, respectively, to the Upper Pleistocene and the Early Holocene, while Magliulo et al. [34] highlighted the presence of three generations but without providing chronological constraints. Recently Amato et al. [15] highlighted the presence of four generations of superimposed or entrenched alluvial fan systems grown above the Laiano Synthem slope breccias. The first generation was chronologically constrained to the late Middle Pleistocene-early Upper Pleistocene, based on 158 ± 6−113 ± 7 ka 40 Ar/ 39 Ar dated tephra layers interbedded in its upper part. The second generation was referred to the Upper Pleistocene thanks to the presence of a 48 ± 7 ka 40 Ar/ 39 Ar dated tephra layer and a CI tephra layer (~39 ka) interbedded in the upper part. The third generation was constrained to the end of the Upper Pleistocene-Holocene based on the presence of a Neapolitan Yellow Tuff tephra layer (NYT hereinafter; ~15 ka [52]) in its upper part. Finally, the fourth generation was referred to the Holocene.
Ancient terraced gravelly-sandy alluvial deposits crop out extensively along the right Calore valley flank between the villages of Ponte and Solopaca [29][30][31][32][33][34] (sector A, Figure 1b) between ~65 m and ~220 m a.s.l. According to Magliulo [30] and Magliulo et al. [34], these deposits are the product of a single, prolonged phase of valley aggradation during the Middle Pleistocene (560 ± 2 ka 39 Ar/ 40 Ar age of an interlayered ignimbrite deposit [30]) followed by several phases of terrace formation due to valley incision during the Middle and Upper Pleistocene. Conversely, Massa et al. [31] and Di Bucci et al. [32,33] consider these terrace orders the result of separated depositional/erosional events. According to Di Bucci et al. [32], the sedimentation of these ancient alluvial deposits can be chronologically constrained between the Middle and the Upper Pleistocene, thanks to the presence of volcanic deposits (a 674 ± 10 ka 40 Ar/ 39 Ar dated pyroclastic layer, a tephra layer in the upper part of the succession with an 40 Ar/ 39 Ar age of 97 ± 25 ka, CI deposits covering the succession). A Middle Pleistocene age for the oldest alluvial terraced deposits is indicated also by Amato et al. [15] based on the 648 ± 6 ka 40 Ar/ 39 Ar age of an interbedded tephra layer. In sector B, ancient terraced alluvial deposits are represented by smaller outcrops in the piedmont area of the Camposauro Mt. (Solopaca area, Figure  1b), referred to the Middle-Upper Pleistocene S. Leonardo sub-Synthem [45], and by larger surfaces located on the right flank of the Calore valley (Telese Plain area, between Telese and Amorosi, Figure 1b) referred to the Upper Pleistocene [33,53].
In the Telese and Amorosi areas (sector B, Figure 1b), the ancient alluvial deposits are covered by travertines related to the springs of Telese and Grassano, located along the southern slope of Pugliano Mt. [54,55] and fed by the Matese Mts. Aquifer. The Amorosi travertines are mainly characterized by phytohermal facies and covered by CI pyroclastic flow deposits forming a wide terraced surface, the Amorosi terrace. These travertines are referred to the Amorosi sub-Synthem [45] and dated to the Upper Pleistocene [45,53,56,57]. The Telese travertines form a large, terraced surface (hereinafter Telese Terrace) on which the town of Telese is placed. They mainly consist of lacustrine facies passing laterally and upwards to swamp facies [56]. Near the Telese springs, travertines are still forming today and Aiello et al. [56] and Del Prete et al. [57] suggest the possibility that the formation of the Telese Terrace occurred even after the CI deposition.
CI pyroclastic flow deposits often form terraced surfaces and are widespread in the Calore and Volturno confluence area and in the middle Volturno River valley (sectors B and C) where they unconformably overlie the oldest Quaternary deposits or the pre-Quaternary bedrock. Reworked CI tephra layers embedded in alluvial deposits prevail, instead, in sector A [15]. Alluvial valley floor deposits in the lower Calore River valley (sectors A and B in Figure 1b) are mainly represented by silty sands and gravels. Two terrace orders are distinguished by [30][31][32][33][34]: Massa et al. [31] and Di Bucci et al. [32,33] consider them as depositional in origin, referring the older one to the late Upper Pleistocene and the younger one (containing reworked terracotta fragments with maximum ages < 7 ka) to the Holocene, while Magliulo [30] and Magliulo et al. [34] interpret the upper terrace as erosional in origin and the lower one as depositional, referring them respectively to the late Upper Pleistocene-Early Holocene and the Middle-Late Holocene. In sector A, the two valley floor terraces complete the flights of fluvial terraces that characterize the right flank of the lower Calore River valley, made of a total of five terrace orders referred to the Middle Pleistocene-Holocene time interval [30][31][32][33][34].
Alluvial valley floor deposits present downstream the Calore confluence (sector C in Figure 1b) are referred to the Upper Pleistocene-Holocene Limatola Synthem [45]. They consist of sandy silts and sands of prevailing pyroclastic nature, with carbonate or polygenic gravelly lenses. Near the Calore confluence, Massa and Zuppetta [35] distinguished within these deposits two orders of fluvial terraces that they referred, respectively, to the Upper Pleistocene and the Holocene.
However, the examined literature highlights several divergences between the various authors as for the correlation of terrace remnants and their inclusion in one order or another. Such differences in grouping terrace remnants, regarding above all the oldest ones, are probably mainly related to their high fragmentation due to valley downcutting by the right tributaries of Calore River, and to tectonics. In particular, as reported by Di Bucci et al. [32,33] and Magliulo et al. [34], high-angle faults affect the oldest terraces orders.

Materials and Methods
The general picture of the data that emerged from the critical revision of literature represents an essential, even if preliminary result of the present study. On its basis, methods and related research activities were defined to front existing knowledge gaps and improve the overall knowledge framework.

Geomorphological Analysis
The geomorphological analysis mainly focused on the identification of Quaternary fluvial landforms.
We performed a multi-temporal analysis based on large-scale topographic maps and Volturno valley between the Calore confluence and the Triflisco dam, sector C). All historical maps were scanned with a resolution of 300 dpi and georeferenced in the UTMWGS84 coordinate system using the Arc-GIS 10.7 software. The multitemporal analysis of topographic maps allowed dating the most recent valley floor terraces. According to this approach, the terraces that formed later than the year of edition of one of the historical maps (1834-1876), have been referred, for uniformity, altogether to the post-1870s.
Using the ArcGIS 10.7 software, we produced a geomorphological sketch focused on fluvial landforms of the study area. All mapped elements were checked with the help of field surveys and stratigraphic data derived from boreholes. To better correlate the mapped fluvial terrace remnants (especially the recent valley floor terraces), longitudinal terrace and thalweg profiles of Volturno and Calore rivers were reconstructed and plotted on a Cartesian plane reporting the distances on the x-axes and the elevation above sea level along the y-axes. For each terrace, the longitudinal average elevations found along the inner edges were plotted on a plane.

Stratigraphical Analysis
The stratigraphical analysis mainly focused on the definition of the morpho-stratigraphical framework of the study area and the characterization of both the valley fillings and the stratigraphical features of the terraced fluvial landforms. This analysis was essentially based on field surveys and the examination of 82 borehole logs provided by public institutions and private companies. Stratigraphical data extracted from boreholes were analyzed and compared with those derived from outcropping successions. Field surveys also allowed identifying and characterizing different tectonic features that affect the investigated Quaternary successions.

230 Th/ 234 U Datings on Travertines
Uranium-series disintegration method ( 230 Th/ 234 U) was performed on five travertine samples to obtain a chronological framework of the studied deposits. These samples were collected (for precise location see Section 4.2.2. and related figure) in the areas of Telese and Amorosi (sector B, Figure 1b), respectively, along the banks of Telese Lake in the uppermost part of the succession (Tr1 at 47 m a.s.l. and Tr2 at 51 m a.s.l.) and along the scarps of the Amorosi terrace in the middle part of succession (Tr3, at 54 m a.s.l., Tr4 at 53.5 m a.s.l. and Tr5 at 53 m a.s.l.).
The radiometric ages were obtained through alpha-spectrometry using an ORTEC OCTETE PLUS spectrometer equipped with eight BR-024-450-100 detectors from the Geosciences Barcelona (Geo3BCN-CSIC, Barcelona, Spain). The chemical separation of the radioisotopes and purification from travertine samples (~20 g) was conducted following the procedure described by Bischoff et al. [58] and isotope electrodeposition was done according to the method of Talvitie [59], modified by Hallstadius [60]. Absolute ages were obtained employing the software designed by Rosenbauer [61].

230 Th/ 234 U datings
Only the Tr2 sample contained > 15% wt of insoluble residue that impeded obtaining a reliable age. The other four samples yielded the radiometric ages reported in Table  1. Among these four travertine samples, Tr1, Tr4, and Tr5 samples contain a significant amount of siliciclastic material that may contribute with terrigenous 232 Th adsorbed in the octahedric layers of detrital clays. However, the relation between the 230 Th from the disintegration of the uranium and the detrital 232 Th from Amorosi travertines is higher than two, which means that contamination can be considered negligible. That is not the case of the Tr1 sample. Moreover, the ratio 230 Th/ 234 U of Amorosi samples is 0.55 ± 0.02, indicating that the radiometric clock has remained closed since the formation of these travertines whereas this ratio is 1.30 ± 0.05 for the Tr1 sample, suggesting that it might have suffered post-depositional alterations. Therefore, from an analytical point of view, the obtained radiometric ages for the Amorosi samples can be considered valid because the detrital contamination is relatively small and the radiometric clock of the U-Th system has been closed since the travertine formation, while the radiometric age obtained for Tr1 must be considered with caution.

Morpho-Stratigraphical Data
From a morpho-stratigraphical perspective, different landforms and relative deposits characterize individually the three study sectors as shown in Table 2 that also reports the associated codes hereinafter used in the text to indicate them. For each sector, the illustration of overall acquired data is supported by a geomorphological sketch focused on fluvial landforms, a related schematic geological cross profile and borehole logs considered in this study.  Some of the highest HFT remnants are cut into bedrock and thinly veneered by a few meters of polygenic fluvial gravels. Frequently, the HFT are reduced to narrow ridges, separated by up to 80 m deep valleys incised by the tributaries of Calore River. These HFT are quite difficult to correlate, given their poor state of preservation and their high fragmentation, which is due-in agreement with other authors (see Section 2.1)-also to tectonics. Therefore, a careful revision of HFT orders was carried out in this study. Thus, in the geomorphological sketch (Figure 2), we have distinguished two groups of HFT. The first group (HFT-I group) includes terrace remnants located between 100 and 220 m a.s.l. (at 60 to 140 m above the present channel; a.p.c. hereinafter). The second group (HFT-II group) is made of gently inclined terrace remnants located at 80/90 m a.s.l. (at 10 to 35 m a.p.c.).
The HTF-I group includes many of the I and II order and some of the III order terrace remnants of Massa et al. [31] and Di Bucci et al. [32,33]. It also includes the I, II, and III order terraces located between Ponte and Masseria Acquafredda, and the I and II order terraces between Masseria Acquafredda and Tuoro S. Antuono of Magliulo et al. [34].
The ancient alluvial (AL) deposits crop out from 65 to 220 m a.s.l. They mainly consist of moderately well sorted gravels in silty-sandy matrix covering the Meso-Cenozoic bedrock. Their maximum thickness in outcrop is ~80 m but exceeds even 100 m in borehole (see B36 in Figure 3). These could be distinguished in two different depositional cycles, named alluvial unit 1a and alluvial unit 1b (Figures 3 and 4) and associated respectively to the HFT-I and HFT-II group.  Alluvial unit 1a mainly consists of polygenic, heterometric, predominantly well-rounded pebbles, poorly cemented, with scarce to abundant silty-sandy and/or silty-clayey yellowish matrix, and lenses and layers of coarse to medium sands (few decimeters to some meters thick) (Figure 5a,b). Horizontal bedding imbrications are locally present in gravels levels and reddish erosion surfaces often mark the contact between sandy levels and gravels. In addition, the borehole successions highlight the presence of fine-grained intervals, containing beds of blackish organic substance and/or peaty layers. The overall textural characteristics of these deposits suggest a braided stream environment. Locally, paleosols and tephra layers are intercalated in unit 1a. In particular, at Masseria Acquafredda (S16 in Figure 2; Figure 5c) a 30 cm-thick yellow-brownish 648 ± 6 ka 40 Ar/ 39 Ar aged tephra layer [15], most likely corresponding to the pyroclastic layer dated to 674 ± 10 ka by Di Bucci et al. [33] (S17 in Figure 2), is interbedded in these deposits ( Figure 4a). Furthermore, unit 1a contains the 560 ± 2 ka aged ignimbrite deposit of Magliulo [30] and the 97 ± 25 ka aged tephra layer of Di Bucci et al. [33] (Figures 2 and  4a). Therefore, a Middle-Upper Pleistocene age between ~650 ka and ~97 ka can be assumed for this unit and the HFT-I group. This age is also supported by the presence of the CI (~39 ka) pyroclastic flow deposits (Figures 4a and 5d) filling the V.ne Ariola stream valley, which deeply cuts the HFT-I group, confirming that its formation took place long before the CI deposition.
Ancient alluvial fan bodies, made of massive, coarse, and medium, poorly rounded polygenic gravels in sandy-silty matrix (Figure 5e), close the fluvial successions that form the highest remnants of the HFT-I group, evidencing that the formation of the latter was influenced also by right tributary inputs of Calore River.
Our field surveys have evidenced, in agreement with other authors [15,[32][33][34], that high-angle faults have displaced the alluvial unit 1a and HFT-I group (Figure 5f,g,h). These faults also controlled the development of the rectilinear, about 25 m high fluvial scarps that terraced the HFT-I group. At the foot of these scarps, the HFT-II group and related alluvial unit 1b are found (Figure 4a).
The HFT-II group includes some of the terraces attributed to the III order by Massa et al. [31], Di Bucci et al. [32,33] and Magliulo et al. [34].
Collected stratigraphical data highlight that the alluvial unit 1b, associated with the HFT-II group, reaches at least 30 m in thickness, and is characterized by coarser deposits, mainly gravelly-sandy, with predominantly carbonate pebbles, which unconformably overlie both unit 1a and the Miocene siliciclastic deposits. The overall textural characteristics of unit 1b indicate braided stream environments.
In Tuoro S. Antuono locality, unit 1b is covered by CI deposits [15] (S2 in Figure 2) that crop out at the top of a terrace referred to the HFT-II group (Figure 4b). Therefore, in accordance with Di Bucci et al. [32,33], unit 1b and the HFT-II group can be constrained to the Upper Pleistocene, between ~97 ka and 39 ka (Figure 4a). The HFT-II group terraces are limited towards the Calore valley floor by degraded rectilinear scarps that are controlled by mainly NW-SE trending faults ( [15] and reference therein), highlighting that tectonics acted even after the CI deposition.
Three orders of valley floor terraces (VFT, Table 2; T1, T2, and T3, Figure 2) have been recognized in all sectors ( Figure 6). Based on previous and new morpho-stratigraphical data and considerations, we interpret these VFT, in agreement with other authors [31][32][33]35], as cut and fill terraces. However, collected borehole and field data did not always allow us to clearly distinguish the alluvial fill associated with each order of VFT; therefore we referred all alluvial deposits underlying the single terraces (T1-T3 alluvial deposits, Figure 4) to a generic alluvial unit 2 that can be clearly distinguished from the ancient alluvial deposits (alluvial unit 1, Figure 7).  In sector A, these terraces are steeper and convergent along the first 7.5 km of the longitudinal profile of the Calore River (Figure 6), then show gentler gradients and remain approximately parallel. Average slopes slightly decrease from T1 to T3 (Table 3). T1 is the terrace with the largest areal extension and is present on both valley sides. T2 is well developed especially on the right valley side and reaches a maximum width of 1000 m. Finally, T3 is rather continuous and present on both sides of the valley. It converges downstream towards the Calore River channel. Its average width ranges between about 50 m in the Ponte area and 700 m downstream. Topographic historical maps show that this terrace order formed after the 1870s.
Only the uppermost layers (~30 cm) of the deposits forming these three terraces are directly observable. They consist of loose sandy silts with gravels and include reworked terracotta pebbles in T3. Some more information was obtained for T1 by examining borehole data. The deposits forming T1 (alluvial unit 2 in B87, B33, B32, B88; left to right in Figure 7) are 5 to 9 m thick and rest above Middle-Upper Pleistocene, up to ~35 m thick alluvial deposits (unit 1 undifferentiated and unit 1b in Figure 7) and consist of both fine and coarse sediments. Fine sediments consist of silts, sandy-silts, silty-sands, and silty-clays with plane-parallel lamination (B87 and B33, left to right in Figure 7), sometimes including polygenic pebbles in the upper part, while coarse deposits essentially consist of loose clast-supported heterometric and predominantly rounded carbonate pebbles (B32 and B88, left to right in Figure 7).
In some places, colluvial deposits cover the inner edge of T1 terrace. In the Masseria Acquafredda area, these deposits contain a reworked NYT tephra layer [15] (S8 in Figure  2; Figures 4a,b and 8a).  Table 3; Figure 4a). In our geomorphological sketch (Figure 2), we have distinguished three generations of AF.
The first and second generation essentially correspond to the first and second generation of Amato et al. [15], the third one includes the third and fourth generation of these authors. The first generation of AF consists of the largest and steepest alluvial fans found in the study area, sometimes deeply incised and Middle Pleistocene-early Upper Pleistocene in age [15]. In agreement with Amato et al. [15], the gravel deposits of this fan generation either are in lateral heteropy with or overly the upper part of alluvial unit 1a (Figure 4a). Like unit 1a, this generation of AF hosts reworked or primary pyroclastic layers (Figures 2 and 4a) and paleosols. Furthermore, it seems to have registered the same tectonic phase of the HFT-I group. In fact, it is affected by primarily ~E to W and secondarily NE-SW oriented high-angle faults, which have also controlled the development of largely degraded, rectilinear scarps (15 to 30 m high) terracing in some places the alluvial fans.
The second generation of AF is entrenched in the first generation and, in some places, incised by channels and affected by E-W and NE-SW high-angle faults, as reported also by [15]. This second generation of AF consists of thinner deposits and is constrained to the late Upper Pleistocene [15] based on the occurrence of the 48 ± 7 ka 40 Ar/ 39 Ar dated tephra layer (Figure 4a) and the CI layer interbedded in its upper part. Moreover, it rests on the Middle-early Upper Pleistocene alluvial unit 1a, while the three orders of valley floor terraces (VFT) are entrenched within it. Therefore, the growth of this generation of AF is most likely contemporary up to subsequent (it contains reworked CI tephra deposits) to the formation of the Upper Pleistocene alluvial unit 1b and the HFT-II group.
Finally, the third generation of AF includes small and low gradient fans that are often connected downslope to the Calore River alluvial plain. On the left valley side, this generation of AF mostly coincides with the third generation of Amato et al. [15], chron-ologically constrained to the end of the Upper Pleistocene-Holocene based on the presence of reworked NYT tephra layers in their upper part. These fans are either superimposed or entrenched in the second AF generation and their deposits are either in lateral heteropy with or overlie the deposits of T1 (Figures 4a and 8b), while T2 terraces are mostly entrenched in them (Figure 4a).
The third generation of AF, unlike the other two, is present also on the right Calore valley side (Figure 4a) and corresponds to the fourth generation of Amato et al. [15]. It consists of alluvial fans made of gravelly-sandy deposits with heterometric, sub-angular to sub-rounded pebbles, containing paleosols, archaeological remains and recent anthropic artifacts (Figure 8c). The deposits of this generation of AF cover T1 and T2 terraces and, in some places, interfinger with the T2 deposits (Figure 4b). This generation of AF is chronologically constrained to the Holocene.

Sector B
This sector (Figure 9) is dominated by the Telese Plain that occupies a large morpho-structural depression (Figure 1b). In the Telese Plain area, the ancient fluvial terraces (HFT) are represented by large surfaces extending between 95 (inner rim) and 75/70 (outer rim) m a.s.l. (Figures 9 and 10), gently sloping toward and located about 50-30 m above the Volturno and Calore valley floors. NW-SE and NE-SW oriented rectilinear scarps (about 5 m high) bound downstream these terraces. Borehole data (B47, B53, B54, and B55 in Figure 11) highlight that their deposits (at least 30 m thick and referred to alluvial unit 1a), are made of gravelly-sandy successions. Gravel intervals mainly consist of polygenic, heterometric, predominantly well-rounded pebbles with silty-sandy and/or silty-clayey yellowish matrix and levels of coarse to medium sands. The fine-grained, often dominant intervals (up to 13 m thick), consist of alternations of clayey silts, clays, silty-sandy clays, and sandy silts.   The HFT are incised by deep, partly wide valleys filled with CI pyroclastic flow deposits. This suggests the correlation of these terraces with the Middle-Upper Pleistocene HFT-I group found in sector A, confirming the attribution of their deposits to alluvial unit 1a.
On the left side of the Calore River valley, the HFT are represented by a few small remnants located between 150 and 90 m a.s.l. at the toe of the Camposauro Mt. in the Solopaca area (Figures 9 and 10). Gravels made of polygenic, well-rounded, heterometric pebbles with abundant silty-sandy matrix and local lenses of sands are associated with these surfaces. These deposits either are in lateral heteropy with or covered by the Middle Pleistocene-early Upper Pleistocene deposits of the first generation of alluvial fans (AF, Figure 10, see later), while the late Upper Pleistocene-Holocene second and third generations of AF cover the associated HFT surfaces. Therefore, the chrono-stratigraphical data allow correlating also these terraces with the Middle-early Upper Pleistocene HFT-I group, and the related deposits with alluvial unit 1a. Rectilinear scarps, up to 40 m high and NW-SE and NE-SW oriented, border downstream these HFT and highlight high-angle faults [15] affecting unit 1a. Instead, no fluvial terraces referable to HFT-II group were found. CI pyroclastic flow deposits form terraces (CI-T, Table 2) which are represented by very large surfaces lying at 65-70 m a.s.l. both along the right side of the Calore valley and the left side of the Volturno valley, as well as by smaller strips located at 60-65 m a.s.l. on the left side of the Calore valley ( Figure 9).
The maximum thickness of exposed CI deposits is 25 m, while borehole data indicate a thickness of more than 25 m (B50, B56, upper line, left to right in Figure 11). The CI deposits unconformably overlie alluvial unit 1a (B48, upper line, B6, B7, B64, B72, second line, left to right in Figure 11; Figure 12a) and reach their greatest thicknesses where they filled valleys incised into the HFT-I group, such as the large and deep valley south of San Salvatore Telesino (B50 to B56, Figures 9, upper line in Figure 11). On the left side of the Calore River valley (Camposauro mountain front), CI-T develop at the foot of NE-SW oriented, rectilinear, and partly degraded scarps, and are covered by the second generation of alluvial fans (AF, Figure 9, see below).
The CI deposits that form the Amorosi terrace (70 m a.sl.) cover the Amorosi travertine deposits (Figures 10 and 12b, ATD in Table 2 and Figure 11). The latter mainly consist of phytohermal travertines and travertine sands reaching a maximum thickness of 5 m in boreholes and at least 15 m in outcrop. The ATD overlie alluvial deposits that we referred to the alluvial unit 1 undifferentiated (B42, B65, and B84, bottom line, left to right in Figure 11) since field and borehole data do not allow to establish which of the two depositional cycles (unit 1a or unit 1b) of the ancient alluvial deposits (AL) they represent.
The new 230 Th/ 234 U age of about 80 ka (samples Tr3-Tr5, Table 1) allows to chronologically constrain the ATD to the late Upper Pleistocene. Considering their chrono-stratigraphical position, the deposition of these travertines is most likely contemporaneous to that of the alluvial unit 1b found in sector A. Locally, alluvial-marshy deposits ( Figure 11) are in lateral heteropy with the ATD and CI deposits.
The Telese travertine terrace (T-T, Table 2, and Figure 9) extends between 65 and 50 m a.s.l. and marks the depositional surface of the Telese travertines which consist of phytohermal travertines and travertine sands (Figure 12c) reaching a maximum thickness of about 15 m (Figure 13). These travertines rest on alluvial-marshy deposits that in turn overlie fluvial deposits referable to alluvial unit 1; more rarely they directly overlie allu-vial unit 1 (B67 and B87, left to right in Figure 13). In the upper part of these travertines, locally alluvial-marshy deposits made of sandy silts containing volcaniclastic material and reworked CI clasts (B26, B25, B62, B76, and B75, left to right in Figure 13) are interbedded. They suggest an age younger than 39 ka for the travertines overlying them, confirming, in agreement with Aiello et al. [56] and Del Prete et al. [57], that the T-T is at least in part younger than the ATD. Therefore, considering that in this area the sedimentation of travertines is still active, a late Upper Pleistocene-Holocene age can be assumed for the T-T. Along the Calore River, the valley floor terraces (VFT) are essentially represented by T2 and T3 (Figures 9 and 6). T1 terrace is only marginally present to the east of Telese, at 52 m a.s.l. (Table 3), and disappears immediately upstream Telese travertine terrace (T-T). Its altitude, close to that of the lower part of the T-T, suggests that its alluvial deposits are at least partly in lateral heteropy with the Telese travertines. T2 reaches a maximum width of 750 m. T3 terrace is approximately continuous and has an average width of 800 m immediately downstream of the T-T.
Along the Volturno River, T1 is lacking and T2 and T3 are convergent ( Figure 6). T2 is present between 49 and 36 m a.s.l. and mostly developed on the right valley side, where it reaches a maximum width of 800 m. T3 is well represented on both valley sides at elevations between 48 and 36 m a.s.l.
In both river valleys, T2 and T3 terraces are entrenched in the CI terraces (CI-T), located at the foot of vertical scarps, not more than 15 m high in the Volturno valley and up to 25 m high in the Calore valley. In addition, in the Calore River valley the T2 terrace is also entrenched in the Telese travertine terrace.
The deposits of T2 (alluvial unit 2 in Figure 7) reach a maximum thickness of 6.5/7 m and mostly consist of medium to coarse sands and sandy-clayey silts predominantly of pyroclastic nature. These deposits rest on CI deposits (with a maximum thickness of 8.5 m in the Volturno valley) that in turn overlie Middle-Upper Pleistocene alluvial unit 1 (at least 23.5 m and 16.5 m thick respectively in the Calore and Volturno valleys).
At the foot of the Amorosi terrace, T2 deposits include, as reported by Calandini et al. [62], a layer containing the Pomici di Avellino deposits (3.95 ka [63,64]) and a silty clay layer rich in ceramic fragments and organic matter. The latter has been radiocarbonically dated to 3500 ± 60 BP [65], i.e., to 3770 ± 150 cal BP (age obtained by using CALIB 8.2 [66] for calibration). These layers, covered by at least two alluvial strata, unconformably overlie an accumulation of CI blocks [62] and highlight that the age of T2 deposits is at least partly referable to the Bronze Age. Finally, a post-1870 age can be attributed to T3 in both river valleys based on editing dates of I.G.M.I. historical topographic maps.
In this sector, all three generations of AF are present on the left side of Calore valley, while on the left Volturno valley side only the third generation is found. As mentioned above, the upper part of the Middle Pleistocene-early Upper Pleistocene first generation of AF covers the HFT-I group (Figures 9 and 10). The late Upper Pleistocene second generation of AF covers both the HFT-I group and the CI-T, highlighting that its growth started before the emplacement of the CI pyroclastic flow deposits and continued even after. The late Upper Pleistocene-Holocene third generation of AF rests on the HFT-I group and on T2 and T3 terraces and alternatively covers or is entrenched in the CI-T, as can be easily observed along the left Volturno side.

Sector C
In this sector, the ancient fluvial terraces (HFT) are lacking and the CI terraces (CI-T) are the highest and oldest surfaces ( Figure 14). However, several boreholes located in the valley floor area intercept ancient alluvial deposits (AL). Furthermore, the three orders of valley floor terraces (VFT) are well represented, while only the third generation of alluvial fans (AF) is found. Borehole data suggest that the Quaternary alluvial infill in this sector of the Volturno valley reaches at least 40 m in thickness (Figures 7 and 15). In addition, a borehole drilled on the T1 terrace in the Piana di Monte Verna area (borehole number 16952 [67]), crosses even 70 m of undifferentiated Quaternary sediments before reaching the carbonate bedrock. Part of this infill is represented by pre-CI alluvial deposits (B27 and B29 in Figure 7), which reach a thickness of at least 24 m and can be referred to the AL alluvial unit 1 undifferentiated. The latter consist of polygenic, heterometric, well-rounded gravels in silty-sandy matrix, containing layers of coarse to medium sands and fine-grained intervals made up of sandy silts, clays, and silty clays, frequently plane-parallel laminated, with beds of blackish organic substance and/or peaty layers. CI-T are represented by discontinuous, small remnants standing between 44 and 66 m a.s.l. and 16-18 m a.p.c. (Figure 14). These terraces are often limited towards the valley floor by up to 15 m high fluvial scarps.
CI deposits forming these surfaces have maximum thicknesses of about 20/25 m in outcrop and rest on fluvial deposits referable to the alluvial unit 1 undifferentiated (Figures 7 and 15).
All three VFT orders are entrenched in the CI-T (Figure 15). T1 is present exclusively between Limatola and Piana di Monte Verna (Figure 14), at an altitude of 37-36 m a.s.l. (Table 3; Figure 6). It is represented by small terrace remnants, often of limited width, except for those in the Piana di Monte Verna area that reach a width of 800 m. T2 ( Figures  6 and 14) is the most continuous and widest terrace (up to 1700 m wide in Piana di Monte Verna area) and ranges in elevation between 36 and 29 m a.s.l. T3 (Figures 6 and 16a) is also continuous and its elevation ranges between 33 and 27 m a.s.l.
T1 alluvial deposits are mainly of pyroclastic nature and made of fine to coarse sands, thinly layered sandy clayey silts, silty clays, and silty sands. Borehole logs evidence that these deposits (alluvial unit 2 in Figure 7) reach a thickness of about 20 m and rest on alluvial unit 1 (Figure 15).
Boreholes drilled on the T2 terrace ( Figure 7) evidence for alluvial unit 2 a thickness between 4.5 and 17.5 m, except for the area next to the Triflisco dam where it reaches even 35 m. Alluvial unit 2 rests on CI deposits near the Calore confluence (B29 and B27 in Figure 7), downstream the CI-T and in correspondence of the inner edge of the T2 terrace. CI deposits are 7.3 to 10 m thick and cover alluvial unit 1. Borehole logs show that alluvial unit 2 is made of fine-grained sediments with a strong volcaniclastic component, mainly consisting of clayey silts, sandy silts, sandy clays, and silty clays, but also of fine to coarse sands. Only in borehole B29 (Figures 7 and 14), probably because of its location near the Maltempo Stream confluence, unit 2 consists of clast-supported, sub-rounded carbonate heterometric gravels with pebbles even exceeding a size of 10 cm.
Medium-fine sands and greyish clay-silts containing a reworked block of NYT, present in the uppermost part of the B27 borehole sequence, confirm that the age of T2 is certainly younger than 15 ka. This is also confirmed by the fact that along the Biferchia Stream the T2 terrace is entrenched in alluvial deposits that can be chronologically constrained between the deposition of CI and NYT [45]. Furthermore, the presence of fragments of pottery and building material referred to Roman farms dated between the II century BCE and II century CE, found on T2 remnants located south of Caiazzo and Piana di Monte Verna villages [68], allow establishing that T2 terraces surely formed before the II century BCE.
T3 deposits (Figure 16b,c) are exposed near the Volturno riverbed along up to 4 m high fluvial scarps and made of sandy silts and fine sands, often thinly laminated, of predominantly pyroclastic origin. In some areas, alluvial fans cover both the T1 and T2 terraces. West of Piana di Monte Verna, the upper parts of these fans contain reworked NYT tephra deposits, confirming their attribution to the third AF generation.

Discussion
According to Amato et al. [15], the formation of the Volturno and Calore paleovalleys started during the Lower Pleistocene and was controlled by high-angle faults. In the Calore river valley (sectors A and B), extensional tectonics favored the progressive deepening and widening of the valley floors, enhancing in time the local relief that favored the deposition of the Laiano Synthem slope breccias.
Overall now available morpho-stratigraphical data and chronological constraints allow proposing a reconstruction of the evolution of the lower Calore and middle Volturno valleys from the Middle Pleistocene onwards. Several morpho-evolutive steps (Stages A-G, Figure 17) have been individuated, providing new clues on the role of Quaternary tectonic phases, climate changes and volcanic activity in the evolution of these valley systems. Below, the reconstructed stages are illustrated. Stage A. This stage represents the first phase of valley evolution during the Middle Pleistocene. In the Calore River valley (sectors A and B, Figure 17), ~650 ka ago a braided stream depositional environment had already started to develop, as highlighted by the oldest radiometric age (~650 ka) and the textural characteristics of alluvial unit 1a. During the Middle Pleistocene, the Calore valley was in prevailing aggradation. This aggrada-tion phase most likely encompasses several climate cycles of Middle Pleistocene and was favored by the activity of high-angle faults mainly affecting the Camposauro mountain front, as highlighted by the faulted slope breccias of the Laiano Synthem.
Pending the lack of stringent chronological constraints for the age of alluvial unit 1 undifferentiated, we hypothesize that sector C was occupied only by a narrow valley affected by prevailing erosion during this stage ( Figure 17).
Stage B. Between ~650 ka and ~97 ka (Figure 17), the formation of the most ancient fluvial terraces (HFT-I group) took place in the Calore River valley and, towards the end of the Middle Pleistocene (~150 ka), also the growth of the first alluvial fans (first AF generation) was fully underway. The deposition of the lower part of these alluvial fans most likely started already in an earlier period of the Middle Pleistocene [15]. The 100 ka climatic variations of the Middle Pleistocene certainly contributed to the development of the HFT-I group and the first AF generation. In particular, the penultimate glaciation (MIS 6) may have significantly concurred to the development of the upper part of this alluvial fan generation.
In sector A, the deposition of alluvial unit 1a and the formation of the related HFT-I group occurred in an asymmetrical valley, affected by a progressive southward migration of its axis, in sector B, instead, in a large morpho-structural depression (Telese Plain area). During this time interval, extensional tectonics mainly acted along NW-SE and NE-SW trending high-angle faults that affected both the HFT-I group and the first AF generation. Tectonics also controlled the development of the straight fluvial scarps, tens of meters high that border the HFT-I group and the first AF generation towards the modern valley floor. This tectonic activity is most likely responsible for the fragmentation of the remnants of HFT-I group, as well as for the different altitudes at which its remnants are located in sector B on opposite sides of the Calore valley.
In sector C, the thick (up to 70 m) pre-CI alluvial infill (alluvial unit 1 undifferentiated; Al-unit 1 in Figure 17) highlights that a watercourse crossed the Volturno valley long before 39 ka. Based on the data available to date, it is difficult to establish if fluvial sedimentation started in this sector already at the beginning of this stage; probably it started only during the late Middle Pleistocene. Moreover, also Valente et al. [16] highlight the existence of a valley already drained by the Volturno River during this time interval in the sector immediately upstream of the study area.
Stage C. Between ~97 ka and 39 ka, a new phase of valley floor aggradation occurred in the study area. In sector A, the deposition of alluvial unit 1b (Al-unit 1b, Figure 17) took place at the toe of structurally controlled fluvial escarpments bordering the HFT-I group. The textural characteristics of this alluvial unit suggest the persistence of braided stream environments. In sector B, the deposition of the ~80 ka aged Amorosi travertines occurred. These travertines are in lateral heteropy with alluvial-marshy deposits that highlight the presence of fluvio-palustrine environments in the Telese area. Along the right Calore side (Camposauro mountain front), the growth of the second generation of AF, covering the alluvial unit 1a and the HFT-I group, started. Finally, in sector C, valley floor aggradation is still recorded by the deposition of Al-unit 1. This phase of valley floor aggradation probably encompasses several episodes of alluviation that occurred during periods of cold/wet climate and reduced tree cover; in particular, significant river aggradation events in Mediterranean valleys are reported during MIS 5b (or at MIS 5b/5a boundary) and during the first half of MIS 3 [14,69].
Stage D. At ~39 ka, the CI explosive eruption caused the emplacement of large volumes of pyroclastic flow deposits. The latter covered alluvial unit 1b (sector A) and the Amorosi travertines (sector B), and filled the valleys incised in alluvial unit 1a. The CI pyroclastic flow severely modified the fluvial landscape and completely filled both the Calore (in sector B) and the Volturno valleys (sector C) as well as their tributary valleys. The lack of significant CI deposits in sector A might suggest that the Camposauro Mt. partially sheltered it.
After 39 ka, a deep valley incision occurred. The formation of CI terrace (CI-T) occurred most likely immediately after the emplacement of CI deposits and can represent the fluvial response to the valley obstruction caused by the eruption, consisting in the rapid remobilization and erosion of these deposits and consequent restoration of the drainage systems [70][71][72][73]. This phase of incision is well recorded in the lower Volturno alluvial plain. Here, a deep (up 30 m) incised paleovalley of the Volturno River, following the modern river axis, is cut in CI deposits, and sealed by a first generation of alluvial deposits dated to 37 ka [74]. However, in our hypothesis most of the valley incision occurred only later, during the LGM (~21 ka [75]), as also recorded in the lower Volturno Plain [74,[76][77][78][79]. We relate this incision to the upstream propagation of regressive erosion waves from the Tyrrhenian coast, triggered by the base level fall during the LGM low sea level stand. Tectonic subsidence affecting the Campania plain [76,77,79,80] must have further favored the development of retrogressive erosion waves as highlighted by some knickpoints along the longitudinal profile of the Volturno River. Inside the main valleys, instead, the role of intra-chain tectonics is completely marginal and negligible. Due to this phase of valley downcutting, aside from the CI-T formed in sectors B and C, probably most of the HFT-II group formed in sector A. Along the Camposauro mountain front alluvial fan prevailed: the deposition of the second AF generation, which contains reworked CI deposits and covers the CI-T in sector B, continued.
A moderate post-39 ka extensional tectonic phase is highlighted in sectors A and B by NW-SE and by E-W and NE-SW high-angle faults affecting respectively the HFT-II group and the second AF generation. Previous data show that tectonics also affected CI terraces and deposits in the Volturno valley, in sector C and immediately upstream [35,81]. In addition, post-CI fault activity is documented in the Campania Plain [82,83].
Stage E. During this stage, the evolution of the valley floor area and the formation of the valley floor terraces (VFT) started. Furthermore, the deposition of the Telese travertines, which are at least in part in lateral heteropy with the T1 deposits, had already started.
Since no tectonic dislocations are recorded in the three orders of VFT (T1-T3), and in the third alluvial fan generation, climate, and land use changes (natural and man-induced) are supposed to be the main drivers of fluvial dynamics from late Upper Pleistocene onwards.
Regarding the age of T1 deposits, no clear indications are actually available, and the relationships between them and the Neapolitan Yellow Tuff (NYT) remain unclear: in sector A, T1 terraces are covered by colluvium containing reworked NYT deposits, and the T1 deposits are either in lateral heteropy or covered by the third generation of AF that also contains reworked NYT tephra. In our hypothesis, the sedimentation of T1 deposits most likely occurred between the Late Glacial and the early Holocene and, therefore, can probably be related to Vita-Finzi's Older Fill [12,84]. Valley aggradation was likely the response to the deglaciation in the Apennines [85] and favored by increased fluvial discharges coupled with high availability of sediments, conditions that are well documented for this period [14,86].
A successive phase of valley incision caused the formation of T1 terrace. This phase of incision may probably be related to the cold-to-warm climate transition [87,88] that occurred between the Younger Dryas/Holocene boundary and the Holocene climate optimum and was characterized by rainfall increase and progressive expansion of temperate forest, well documented in the central-southern Mediterranean area [89,90].
Stage F. The sedimentation of T2 deposits is linked to a phase of valley floor aggradation that occurred from the Mid to Late Holocene (and certainly extended beyond 3.5 ka, as suggested by the 3770 ± 150 cal BP of organic matter and the Pomici di Avellino level included in their upper part; see Section 4.2.2). This phase of valley aggradation probably encompasses several episodes of alluviation favored by climate deterioration and land cover changes (especially forest reduction). Alluviation was probably favored by cooler conditions such as those recorded in central-southern Italy during the period 7700-5000 cal. BP ( [91] and references therein), and by prevailing wetter, partially also cooler conditions (highlighted by higher lake levels and/or a decrease in fire activity) as those reported for the period 7700-6400 cal. BP [92] and from 5000 to 3500 (3000) cal. BP in central Italy [93]. Furthermore, a link between this phase of valley aggradation and periods of increased flooding recorded in the Mediterranean area and southern Italy in particular ( [94] and references therein), seems likely. Especially, periods of major flooding occurred in southern Italy at 4200-4100 BP and 3450-3000 cal. BP, which overlap with lake highstands at 4200-4000 and 3500-3200 cal. BP in central Italy [94], may have played an important role during the terminal part of said phase of alluviation.
On the other hand, the instauration of drier conditions causing vegetation cover and, especially, forest decline may have favored soil erosion and the increase of sediments deliverable to the watercourses also in latter moments. Such periods of climate aridification characterized the central-southern Mediterranean during the Mid Holocene in various moments, and especially around 5.5 ka and at 4.2 ka [89,90,92,[95][96][97][98]. Aside from climate controls, also human activities such as agricultural expansion, cattle breeding, and timber extraction may have contributed to forest decline. In fact, human presence is well documented in the Calore River valley since Neolithic times [65,99], and human impact became significant at least since the late Holocene [65].
Successive fluvial incision led to the formation of T2 terrace that occurred well after 3.5 ka and before the II century BCE, as highlighted by Roman remains found on the terrace. This phase of incision can possibly be related to the transition towards the "Roman Climatic Optimum" (1 CE-500 CE, [100]).
Stage G. The deposition of T3 deposits may have already started in the late Roman period. It highlights phases of valley floor aggradation due to increased fluvial activity and sediment supply during colder periods, namely the Late Antique Little Ice Age (LALIA, 536-660 CE, [101]) and the Little Ice Age (LIA, 16th-mid 19th century, [102]), corresponding to the so called Vita-Finzi's Younger Fill [12].
The formation of T3 terraces, as already stated, occurred after the 1870s. In the Calore River valley, these terraces largely include the sector of the alluvial plain that, according to Magliulo and Cusano [103], formed during the <1870-1998 CE time span.
The late Upper Pleistocene-Holocene third generation of AF grew during the formation of VFT. In accordance with [104], in the Calore River valley many of the alluvial fans referred to this generation, are still active or have been recently reactivated during the destructive flash flood of 15th October 2015. During this flood event, most of the surfaces referred to T2 and T3 terraces in the Calore River valley were flooded [105].

Conclusions
The present study provides an overall revision of pre-existing and partially contrasting literature data (Section 2.1) on the Quaternary evolution of the Middle Volturno River valley system in a broad sense. This activity allowed getting an order of data, to merge data coming from separate or only partially overlapping area, and was the basis for addressing new research activities, to front existing knowledge gaps and improve the overall frame.
New research activities (field surveys, borehole analysis, 230 Th/ 234 U datings on travertines, etc.) added essential stratigraphical, morpho-structural, and geomorphological data, as well as precious chronological constraints. In particular, the 230 Th/ 234 U ages obtained for the Amorosi travertines allowed constraining them to about 80 ka and, along with morpho-stratigraphical data, separating their deposition clearly from that of Telese travertines.
The integrated analyses of pre-existing and new data allowed defining the major stages of the Quaternary evolution of the study area, highlighting:


Fluvial sedimentation started in the Middle Pleistocene (at least 650 ka ago) within valleys whose formation occurred already during the lower Pleistocene under the control of high-angle faults.  Extensional tectonics acted during the Middle and Upper Pleistocene, driving the progressive widening and deepening of the valleys, and controlling the formation of the HFT and the first and second generation of AF in both the lower Calore valley and the Calore confluence area (sectors A and B). The HFT and alluvial fans were affected by high-angle faults until the deposition of the CI in sectors A and B occurred, while in the Volturno River valley (sector C), probably due to the extensive deposition of the CI pyroclastic flow deposits, no clear evidence of these faults could be found.  At ~39 ka, the emplacement of large volumes of CI pyroclastic flow deposits severely modified the landscape, completely filling the Volturno valley and the Calore valley in sector B.  Extensional tectonics acted also after the emplacement of the CI deposits.  From the end of the Upper Pleistocene onwards (< 15 ka), the role of tectonics in the genesis of terraces and alluvial fans appears to be negligible. Climatic and land-use changes played a key role in the formation of the three orders of valley floor terraces and the third generation of alluvial fans.
In conclusion, our study provides an integral comprehensive picture of the Quaternary landscape evolution of the investigated Apennine valley sector and suitable to efficiently address future studies.
New investigations can in particular contribute to improve the data resolution of the valley floor areas, to better clarify the relationships between alluviation/incision cycles and environmental changes (natural or man-induced) during the late Quaternary. For this purpose, the execution of new boreholes on the valley floor terraces, connected paleo-environmental analyses and radiometric datings (extended also to the Telese travertines), can surely contribute to better characterizing the late Quaternary and, especially, Holocene fluvial dynamics and geomorphological evolution of the study area. Data obtainable can increase our knowledge also on the sensibility to environmental changes of the studied river systems in the frame of global climate change.
Funding: This research was financially supported by University of Naples Parthenope research found grant to P.P.C. Aucelli.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available in the article itself and references cited.

Acknowledgments:
The authors wish to thank the Special Issue Editor and the anonymous reviewers whose suggestions contributed to improving the manuscript.

Conflicts of Interest:
The authors declare no conflict of interest.