# Terrestrial records of deglaciation events during terminations V and IV in the central Apennines (Italy) and insights on deglacial mechanisms

F. Marra, A. Pereira, B. Jicha, S. Nomade, I. Biddittu, F. Florindo, G.

Muttoni, E M Niespolo, P R Renne, V. Scao

► To cite this version:

F. Marra, A. Pereira, B. Jicha, S. Nomade, I. Biddittu, et al. Terrestrial records of deglaciation events during terminations V and IV in the central Apennines (Italy) and insights on deglacial mechanisms. Scientific Reports, 2022, 12 (1), pp.18770. 10.1038/s41598-022-23391-7. hal-03866886

HAL Id: hal-03866886

https://hal.science/hal-03866886

Submitted on 23 Nov 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN

# Terrestrial records of deglaciation events during terminations V and IV in the central Apennines (Italy) and insights on deglacial mechanisms

F. Marra<sup>1,∗</sup>, A. Pereira<sup>2,3</sup>, B. Jicha<sup>4</sup>, S. Nomade<sup>5</sup>, I. Biddittu<sup>6</sup>, F. Florindo<sup>1</sup>, G. Muttoni<sup>7</sup>, E. M. Niespolo<sup>8,9,10</sup>, P. R. Renne<sup>8,9</sup> & V. Scao<sup>5</sup>

<sup>40</sup>Ar/<sup>39</sup>Ar geochronology constraints to aggradational phases and grain size variations show that the two large gravel beds occurring in the sedimentary filling of the Liri fluvial-lacustrine basin (central Italy) recorded the occurrence of deglaciation events synchronous within uncertainties with global meltwater pulses at ca. 450 and 350 ka. In particular, we find a precise match between the ages of gravel deposition and the occurrence of moderate sea-level rise events which anticipate those more marked during the glacial termination V and IV in the Red Sea relative sea level curve, as already verified by data from the Tiber River catchment basin. Such correspondence suggests that gravel deposition is facilitated by melting of Apennine mountain range glaciers, which provide the water transport energy and a surplus of clastic input to the rivers draining the mountain regions and flowing into the Tyrrhenian Sea. Therefore, the thick gravel beds intercalated in the sedimentary filling of the catchment basins of the major rivers in central Italy may be regarded as an equivalent proxy of large deglaciation events, similar to the ice-rafted debris in northern Atlantic. Consistent with this hypothesis, we also show the close correspondence between the occurrence of particularly mild (warmer) minima of the mean summer insolation at 65° N and these early aggradational phases, as well as with other anomalous early sea-level rises occurring c. 750 ka and 540 ka at the onset of glacial termination VIII and VI, and 40 ka at the onset of the so-called Heinrich events.

Geochronologically constrained records of glacial-interglacial variations in ice volume and sea level represent a fundamental tool to decipher Pleistocene global climate evolution. However, records characterized by global rather than local significance are very rare (e.g., ocean cores, ice cores, coral reefs, speleothems) and those with direct, precise radioisotopic age constraints, especially in the time interval > 500 ka, are even less common.

Over 30 years of dedicated studies have shown that fluvial-lacustrine sedimentation within the catchment basin of the Tiber River (Fig. 1) responded synchronously with changes in base-level induced by glacio-eustatic fluctuations during Middle-Upper Pleistocene<fnref n="1" />⁻¹⁵.

The sedimentary record of the Paleo-Tiber River consists of a series of fining-upward sequences of clastic sediments deposited above an erosional surface in response to sea-level rise during the last eight glacial-interglacial transitions ('aggradational successions')¹⁷ (and references therein). Each sequence is characterized by a basal interval of distinct gravel deposition that abruptly switches to clay deposits. In particular, ⁴⁰Ar/³⁹Ar age constraints on tephra layers interbedded with the sedimentary deposits indicate that coarse gravel beds at the base of each aggradational succession are deposited as a consequence of ice-melting during the glacial termination.

<footnote n="1">¹Istituto Nazionale di Geofisica e Vulcanologia, Rome, Italy. ²Université Paris-Saclay, CNRS UMR 8148, GEOPS, Gif-sur-Yvette, France. ³Département Hommes et Environnements, Muséum National d'Histoire Naturelle, Paris, France. ⁴Department of Geoscience, University of Wisconsin-Madison, Madison, USA. ⁵CEA Saclay, LSCE, UMR-8212, UVSQ-IPSL et Université Paris Saclay, Gif-Sur-Yvette Cedex, France. ⁶Istituto Italiano di Paleontologia Umana, Anagni, Italy. ⁷Department of Earth Sciences, University of Milan, Milan, Italy. ⁸Department of Earth and Planetary Science, University of California, Berkeley, USA. ⁹Berkeley Geochronology Center, Berkeley, USA. ¹⁰Department of Geosciences, Princeton University, Princeton, USA. ∗∗email: fabrizio.marra@ingv.it</footnote>
<figure><img src="image_3.png" /></figure>

**Figure 1.** 1) Catchment basins of the Tiber River (black border) and of the Sacco-Liri-Garigliano Rivers (green border); 2) areas with evidence of glacial landforms associated with MIS 14 through MIS 6 glaciations and (3) areas above 1500 m a.s.l. of potential expansion of the glaciers according to Giraudi and Giaccio¹⁶; 4) meltwater and clastic inputs reaching the Liri Basin during deglacial periods; 5) main border faults of the Sacco-Liri basins. The location of the investigated area represented in Fig. 2 is shown. DEM images: TINITALY/01 square WA 6570, used with permission of the Istituto Nazionale di Geofisica e Vulcanologia, Rome. Graphics handrawn by authors.
These results have highlighted that conditions for the accumulation of the coarse gravel deposits only coexisted at the onset of glacial terminations due to several concurrent factors:

1. low sea level at glacial maxima, which steepens the gradients and, in turn, enhances the river competence through the more deeply incised valley;

2. melting of Apennine mountain chain glaciers that releases large amounts of clastic material, increasing the sediment supply to the river drainage basin;

3. overall increase in regional precipitation.

Gravel starts accumulating at the end of the glacial period, when sea-level continues to fall, which caused re-incision of the valley floor and removal of the gravel transported during the regressive phase up to the glacial maximum. The abruptness of the "sedimentary switch" that marks the transition from the gravel bed (2–8 m in thickness), through a thin (<1 m in thickness) sand bed, into the thick (20–40 m in thickness) silt and clay section calls for the sudden establishment of a low gradient, consistent with fast sea-level rise (meltwater pulse) and subsequent development of a sea-level highstand.

Previous work focused on the aggradational successions in the coastal plain and the terminal tract of the Tiber River. More recent studies have highlighted the synchronicity, within uncertainties of age models, between the gravel/clay switch and the peaks of sea-level rise during Marine Isotopic Stage (MIS) 11 and MIS 9 in the higher portion of the Tiber catchment basin, as far inland as 50 km¹⁴ and 100 km⁹ from the coast (Sant'Abbondio and Civitella di Agliano locations in Fig. 1).

Giaccio et al.¹⁴ have remarked on the coincidence among melt-water pulse events, peaks in the ice-rifted debris (IRD) curve¹⁸, and deposition of the gravel beds of the MIS 11 aggradational succession (San Paolo Formation²), suggesting that gravel deposition in the catchment basins of the Tiber River can be regarded as an equivalent proxy of deglaciation events. This composite record of radioisotopically dated (¹⁴C and ⁴⁰Ar/³⁹Ar) morpho-sedimentary units can thus provide key geochronological constraints that are generally lacking in the Middle Pleistocene sea-level records, and they can be used to better evaluate the relationship between insolation changes and sea-level oscillations.

In order to verify the global significance of this kind of proxy and the possibility to use similar morpho-stratigraphic units from other regions to constrain the timing of the deglacial phases, we have investigated the possible response to the glacio-eustatic signal of sediment supply grainsize within the Sacco-Liri-Garigliano Rivers catchment basin, located in central-southern Italy (Fig. 1).

Here, a more than 50 m thick fluvial-lacustrine succession filled the tectonic depression of the Liri basin during the Middle Pleistocene¹⁹⁻²¹.

In this study, we identify two 3–4 m thick beds of coarse gravel deposits intercalated within clayey lacustrine sediment of the Middle-Pleistocene Sacco-Liri basin (Fig. 2). We correlate these gravel beds cropping out at three geologic sections (Colle Avarone, San Giorgio al Liri, Pignataro Interamna) and re-analyze the stratigraphy of six sites described in the literature (Campo del Conte, Cava Pompi, Campogrande, Isoletta, Lademagne, Pontecorvo) located in three different structural portions of the Sacco-Liri basin (Fig. 2). Detailed stratigraphic sketches of these sections are reported in Suppl. Mat. #1.

In order to provide age constraints to the time of deposition of these gravel beds, we produced 8 new ⁴⁰Ar/³⁹Ar ages both on primary volcanic samples and reworked deposits containing potassium (K)-feldspar crystals, which we integrated with 9 ⁴⁰Ar/³⁹Ar and 2 K/Ar dates previously produced from these sections.

A detailed description of the sedimentary features of the Sacco-Liri basin can be found in Devoto¹⁹. This author described the sedimentary succession filling the large lacustrine basin in the "Lower Liri Valley", between Ceprano and the Garigliano River confluence (Fig. 2), as composed of three lacustrine facies with vertical and lateral transition from one to another:

1. *Lower lacustrine mud.* Bedded white calcareous muds with frequent intercalated black tephra.

2. *Typical lacustrine facies.* White calcareous muds with alternating cross-bedded yellow sand layers, black and brown "tuffite" with slumpings, conglomerate.

3. *Late lacustrine facies.* Calcareous muds varying in color with occasional lignite and peat layers, transitioning laterally and vertically to travertine.

The occurrence of conglomerate deposits is only briefly mentioned by Devoto¹⁹ in the description of the "Typical lacustrine facies". Such high energy sediments within a large lacustrine basin may have a differing origin, possibly related to e.g., a local alluvial fan or sedimentary traps. However, our study evidenced that coarse gravel layers are concentrated within two horizons occurring at different elevation. The upper one, 2–3 m thick, can be laterally correlated for several kilometers through the investigated sections. The lower gravel bed, ca. 4 m thick, has been found in a borehole at Ceprano²², where it occurs at the base of a ca. 20 m thick package of lacustrine muds of the "Typical lacustrine succession", and is reported by Devoto¹⁹ to occur in the south-eastern portion of the basin within this sedimentary unit.

Such gravel occurrences are not repeated in the upper, younger levels of the Liri basin succession, suggesting that its emplacement is linked with significant variation in water capacity of transport that affected in two distinct periods of the past the Sacco-Liri catchment basin. We interpret these to correlate with melting of the glaciers occurring in the high ranges of the central Apennines. The presence of glaciers in the central Apennines during past glacial periods has been shown by a review of the existing data for the Middle Pleistocene glacial remnants and new stratigraphic and tephrochronological data which enabled Giraudi and Giaccio¹⁶ to recognize at least five glacial stages during MIS 14, 12, 10, 8 and 6. In particular, one of these areas is located at the northern
<figure><img src="image_4.png" /><figcaption>**Figure 2.** (a) Morpho-structural sketch of the Latin Valley with the Sacco-Liri catchment basin showing location of the investigated geologic sections (white dots). (b) Regional geologic map (created by authors) straddling the investigated area. All graphics handdrawn by authors.</figcaption></figure>

margin of the Sacco-Liri-Garigliano catchment basin (Aremogna, Fig. 1). Glacial landforms associated with the Apennine's glaciers have been reported to occur at elevations as low as 1500 m a. s. l.<sup>16</sup> (and references therein).

**Geological-structural setting.** The Latin Valley (Fig. 2) is part of the central Apennine fold-and-thrust belt, which originated from Late Tortonian-Early Messinian compressional phases and has been affected by extensional tectonics since the Pliocene<sup>23,24</sup>. The outcropping terrains belong to the Latium-Abruzzi neritic carbonate domain (upper Triassic-middle Miocene) and are covered by middle Miocene to lower Pliocene synorogenic siliciclastic deposits<sup>25</sup>.

The study area is organized in several NW-SE striking imbricate thrust sheets that overthrust onto Tortonian-lower Messinian terrigenous deposits, cross-cut by a system of conjugated synthetic and antithetic Quaternary normal faults which controlled the formation and growth of intramountain basins during the extensional phase<sup>26</sup>. The articulated catchment basin of the Sacco, Liri and Garigliano Rivers (Figs. 1, 2) developed in a graben structure within the Latin Valley<sup>26</sup>.

Locally, the study area has also been affected by N- to NNE- and E- to ENE-striking high angle faults with strike-slip kinematics up to the Middle Pleistocene<sup>21,26,27</sup>. The strike-slip tectonics were associated with the erup-tive centers of the Volsci Volcanic Field<sup>28</sup>, the activity of which occurred in three main phases<sup>21</sup>. An early eruptive phase spanning 761.5 ± 9.5 to 541.0 ± 14.0 ka (all errors 2σ) was characterized by long quiescent periods between isolated eruptive events; major eruptive activity occurred from 424.0 ± 13.0 to 349.5 ± 5.0 ka; finally, a late eruptive phase has less precise geochronologic constraints between 300.0 ± 28.0 and 231.0 ± 19.0 ka<sup>21</sup>.

**Results**

<sup>40</sup>Ar/<sup>39</sup>Ar data. Results for the eight samples dated in the present study are reported, along with those of 9 previous samples, as age probability diagrams in Fig. 3 and summarized in Table 1. Full analytical data are reported in Supplementary Material #2.

**Chronostratigraphic analysis of the Sacco-Liri basin.** *Campo del Conte.* A sedimentary sequence about 7 m in thickness, including four fluvial depositional cycles, crops out in the western sector of the Latina Valley in Campo del Conte<sup>31</sup> (Fig. 2). This sequence is characterized by the occurrence of *Mammuthus meridionalis* fossils along with those of a cervid belonging to the Pseudo-dama group<sup>31</sup> indicating an Early Pleistocene
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Results of $^{40}$Ar/$^{39}$Ar analyses on single crystals from this study presented as probability diagrams.</figcaption></figure>

age. More specifically, *M. meridionalis* lived in central Italy during the period chronologically framed between 2.6 and 1.6 Ma³².

Campo Grande (Ceprano boreholes). Two boreholes drilled in the western sector of the Liri basin in Campo Grande (Fig. 2) recovered 48 m of lacustrine succession between 108 and 60 m a.s.l., without reaching its bottom²². Three main successions, the lowest two separated by a ca. 5 m thick gravel layer, were tentatively correlated with the three lacustrine facies described by Devoto¹⁹. Although not strictly applicable everywhere, we
<table><caption>**Table 1.** $^{40}$Ar/$^{39}$Ar sample ages. All ages according to ACs at 1.1848 Ma$^{29,30}$.</caption><thead><tr><th>Sample/site</th><th>Age (ka) ± 2σ</th><th>n/N</th><th>MSWD</th><th>Latitude<br>Longitude</th><th>Elev<br>(m asl)</th></tr></thead><tbody><tr><td>BL-5</td><td>300 ± 12</td><td>2/17</td><td>–</td><td>41° 22' 22.16" N<br>13° 45' 58.36" E</td><td>48.5</td></tr><tr><td>PO-C6</td><td>337.5 ± 6.4</td><td>7/23</td><td>0.90</td><td>41° 33' 31.24" N<br>13° 26' 18.58" E</td><td>129.5</td></tr><tr><td>CE-2</td><td>389.6 ± 3.7</td><td>1/33</td><td>–</td><td>41° 31' 43.35" N<br>13° 28' 45.92" E</td><td>95</td></tr><tr><td>PI-2</td><td>400.9 ± 3.4</td><td>9/11</td><td>0.90</td><td>41° 26' 05.37" N<br>13° 47' 30.12" E</td><td>50</td></tr><tr><td>PI-1</td><td>513.2 ± 7.4</td><td>9/10</td><td>0.70</td><td>41° 26' 05.37" N<br>13° 47' 30.12" E</td><td>48</td></tr><tr><td>BL-4</td><td>516.8 ± 2.1</td><td>28/29</td><td>1.04</td><td>41° 22' 22.16" N<br>13° 45' 58.36" E</td><td>47</td></tr><tr><td>BL-1A</td><td>538.3 ± 7.0</td><td>7/19</td><td>0.56</td><td>41° 22' 22.16" N<br>13° 45' 58.36" E</td><td>40.5</td></tr><tr><td>CE-1</td><td>452.4 ± 1.8</td><td>1/30</td><td>–</td><td>41° 31' 43.35" N<br>13° 28' 45.92" E</td><td>72</td></tr><tr><td colspan="6"><strong>From previous literature</strong></td></tr><tr><td>CA-C1<sup>21</sup></td><td>345.9 ± 4.3</td><td>8/28</td><td>0.82</td><td>41° 32' 11.80" N<br>13° 28' 41.40"E</td><td>119.5</td></tr><tr><td>Ceprano<sup>33</sup></td><td>350.6 ± 8.0</td><td>14/17</td><td>1.50</td><td></td><td></td></tr><tr><td>CA-CGT<sup>21</sup></td><td>359.6 ± 6.5</td><td>3/24</td><td>0.04</td><td>41° 32' 11.80" N<br>13° 28' 41.40" E</td><td>112.5</td></tr><tr><td>Isoletta 3<sup>35</sup></td><td>362.6 ± 3.8</td><td>2/8</td><td>0.34</td><td>41° 31' 44.86" N<br>13° 34' 0.16" E</td><td>112</td></tr><tr><td>Lademagne 2<sup>35</sup></td><td>386.2 ± 4.6</td><td>7/17</td><td>1.60</td><td>41° 31' 16.68" N<br>13° 35' 0.72" E</td><td>105</td></tr><tr><td>Isoletta 2<sup>35</sup></td><td>373.2 ± 2.8</td><td>10/13</td><td>1.40</td><td>41° 31' 44.86" N<br>13° 34' 0.16" E</td><td>104</td></tr><tr><td>Cava Pompi<sup>35</sup></td><td>394.3 ± 8.2</td><td>8/8</td><td>0.25</td><td>41° 33' 31.24" N<br>13°26' 18.58" E</td><td>127</td></tr><tr><td>Isoletta 1<sup>35</sup></td><td>400.3 ± 3.0</td><td>11/12</td><td>0.85</td><td>41° 31' 44.86" N<br>13° 34' 0.16" E</td><td>92</td></tr><tr><td>Lademagne 1<sup>31</sup></td><td>402.4 ± 4.8</td><td>7/15</td><td>1.10</td><td>41° 31' 16.68" N<br>13° 35' 0.72" E</td><td>101.8</td></tr></tbody></table>

will use these three successions, here termed "lower lacustrine", "middle lacustrine-fluvial", and "upper fluvial-lacustrine", as reference chronostratigraphic units in this paper.

We have collected two samples from the Campo Grande borehole cores stored at Università degli Studi di Milano. Sample CE-1 was collected in core Ceprano 1 at 39.3 m depth within a coarse gravel layer with abundant sand matrix. The youngest crystal out of a population of 30 extracted from this sediment yielded a $^{40}$Ar/$^{39}$Ar age of 452.4 ± 1.8 ka. Sample CE-2 was collected at 15.1 m depth in borehole Ceprano 2, at the base of a coarse sand layer and yielded a youngest crystal date of 389.6 ± 3.7 ka. As discussed in the previous section, these ages can be considered good approximations of the time of deposition of these sediments (see also "Age interpretations").

Consistent with this hypothesis, a constant sedimentation rate of ~38 cm/ky is calculated from our sample ages combined with that of 350.6 ± 8.0 ka$^{33}$ on the primary volcanic layer cored at the top of the lacustrine succession (Fig. 4), in good agreement with previous estimation (30–40 cm/ky$^{22}$). While linear long-term sedimentation rates are likely to obscure high frequency changes in sediment accumulation as revealed by sedimentological observations, they provide useful constraints on the accumulation history of the basin during the overall lifetime of its existence. Indeed, gravel and cross-bedded sands represent sudden, high-energy inputs within the lacustrine basin. Also in consideration of their small thickness with respect to the silty-clayey sediments, the ages of the samples collected in these two horizons may be regarded as good spot constraints to the overall lacustrine sedimentation, which is characterized by relatively constant sedimentation rate.

These observations indicate that the ages of 452.4 ± 1.8 ka and 389.6 ± 3.7 ka are excellent time constraints on the deposition of the gravel horizon at the base of the "Middle lacustrine-fluvial succession", and to the start of coarse-sized sedimentary input (cross-bedded sands with sparse gravel) at the base of the "upper fluvial-lacustrine succession".

**Cava Pompi.** Temporary archaeological trench excavations at this location$^{34}$ exposed a ca. 3 m thick sedimentary succession overlying a primary pyroclastic-flow deposit (Suppl. Fig. S1). It was constituted, from bottom to top, by a basal, coarse sand-and-gravel volcaniclastic horizon, mostly deriving from reworking of the underlying volcanic deposit. It abruptly passed upwards to bedded, white lacustrine muds and massive, yellowish sandy silt with clay intercalation. A travertine horizon closes the succession.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Synthetic stratigraphic log of the two boreholes drilled in Ceprano<sup>22</sup> showing the samples dated in the present work (red dots) and in Nomade et al.<sup>33</sup> (black dot), which allow to assess the sedimentation rate of the lacustrine succession.</figcaption></figure>

Lower and upper age constraints to this sedimentary succession are provided by two <sup>40</sup>Ar/<sup>39</sup>Ar ages performed on the pedogenically weathered top horizon of the pyroclastic flow deposit occurring at the base of the succession<sup>35</sup>, and on a second, primary volcanic deposit (this work) unconformably overlying it (Suppl. Fig. S1). The lowest sample provides a *terminus post- quem* of 394.3 ± 8.2 ka to the beginning of sedimentation, while the upper sample (PO-C6) gives a *terminus ante- quem* of 337.5 ± 6.4 ka for it.

*Colle Avarone.* A ca. 8 m-thick succession is exposed at several sections at this locality (Suppl. Fig. S2). A coarse gravel in abundant sand-matrix horizon, > 2 m in thickness, occurs at the base and is overlain by a lacustrine succession in which three primary volcanic deposits are intercalated. Two <sup>40</sup>Ar/<sup>39</sup>Ar ages on one sample (CA-CGT) of sand matrix collected in the basal gravel horizon and on the uppermost pyroclastic-flow deposit (CA-C1) were performed<sup>21</sup>, bracketing the deposition of the sedimentary succession between 359.6 ± 6.5 and 345.9 ± 4.3 ka.

*Isoletta.* A more than 30 m-thick section was temporarily exposed during construction of the high-velocity railway in the 90's at this location<sup>36</sup>. A more than 10 m-thick grey clay-and-silt lacustrine to fluvial deposit was exposed at the base of the succession, followed by another ca. 10 m-thick package of cross-bedded to planar coarse sand (Suppl. Fig. S3). A ca. 3 m-thick layer of coarse gravel follows in the succession, transitioning upwards to sand and silt. A travartinaceous silt horizon closes the succession. Two sub-primary (i.e., reworked, in-place ashfall deposit) volcanic layers intercalated at the base of the clay and in the middle of the coarse sand horizons provided ages of 400.3 ± 3.0 and 372.2 ± 2.8 ka, respectively<sup>35</sup>. A reworked volcaniclastic layer collected in the coarse gravel horizon yielded a youngest population of two crystals providing a *terminus post- quem* age of 362.6 ± 3.8 ka for its deposition<sup>35</sup>.

*Lademagne.* The sedimentary succession is made up of a 2 m-thick horizon of sand with abundant medium-to coarse gravel, overlying a clay layer and topped by 1 m of sandy silt deposits and was described at this location<sup>37</sup>. Two sub-primary volcaniclastic layers collected in the vicinity of the original archaeological site, which stratigraphically constrain the sedimentary succession at the bottom and at the top (Suppl. Fig. S4), yielded ages of 402.4 ± 4.8 ka and 386.2 ± 4.6 ka, respectively<sup>35</sup>.
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Stratigraphic sketches of San Giorgio al Liri (a) and Pignataro Interamna (b) sections showing the sampled tephra layers and the average sedimentation rate.</figcaption></figure>

*Pontecorvo.* Two tephra layers interbedded in the carbonate-rich muds of the "Lower Lacustrine Succession" of the Liri basin cropping out in the surroundings of Pontecorvo village were dated by the K/Ar method at 583 ± 11 ka and 570 ± 11 ka¹⁹,³⁸. We have re-investigated this sector and detected a geologic section in which several tephra layers occur at elevation ranging 55–57 m a.s.l. (Suppl. Fig. S5).

*San Giorgio al Liri.* A more than 10 m-thick succession of white carbonate-rich muds, silts and travertine layers cropping out in the neighborhoods of San Giorgio al Liri village was discovered (Fig. 5a). Eight cm- to dm-thick tephra layers are intercalated in the lowest 8 m of the succession which is constituted by lacustrine white carbonate-rich muds. We have dated the lowermost (BL-1A) and uppermost (BL-4) of these tephra layers, which yielded ⁴⁰Ar/³⁹Ar ages of 538.3 ± 7.0 ka and 516.8 ± 2.1 ka, respectively, allowing correlation of the sedimentary deposits with the "Lower lacustrine succession". In contrast, a sedimentary sample collected in the travertine horizon at the top of the succession provided a *terminus post- quem* age of 300 ± 12 ka, evidencing a large sedimentary hiatus at the top of the lacustrine succession.

The ages of the two tephra layers allow an estimate of the sedimentation rate of ca. 31.5 cm/ky for the lacustrine basin during the corresponding time span.

*Pignataro Interamna.* A more than 20 m-thick lacustrine succession crops out in the surroundings of Pignataro Interamna village, which was attributed to the "Typical lacustrine facies"¹⁹. We have investigated the lowest, 10 m-thick portion of this succession cropping out few km southeast of Pignataro Interamna, constituted by white carbonate-rich muds and yellow silts, in which two clastic horizons are intercalated (Fig. 5b). We sampled and dated one primary tephra occurring in the lower part of the lacustrine deposits and the volcaniclastic sand matrix of a discontinuous, up to 20 cm thick gravel layer occurring in the upper portion of the succession. An anomalous, 50 cm-thick conglomeratic horizon constituted by poorly rounded, ≤1 cm sized carbonate fragments, chert and limestone pebbles within a silty matrix occurs ca. 1.5 m above the lowest tephra layer dated at 513.2 ± 7.4 ka (sample PI-1). The maximum age of 400.9 ± 3.4 ka provided by the sample PI-2 collected in the upper gravel layer, 50 cm above the conglomerate, implies that this clastic horizon marks a significant sediment-
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Results of $^{40}$Ar/$^{39}$Ar analyses on single crystals presented as probability diagrams along with plots of $^{40}$Ar$^{\circ}$ and K/Ca ratios. (d) Age populations of sample BL-5 are further distinguished by the amount of radiogenic Ar (%$40$Ar$^{\circ}$); see text for comments and explanation.</figcaption></figure>

tary hiatus, similar to that occurring at the top of the "Lower lacustrine succession" in San Giorgio al Liri. However, more lacustrine deposits occur above this hiatus in Pignataro Interamna, which are a lateral transitional facies of the "upper fluvial-lacustrine succession", as evidenced by the geochronologic constraints provided here.

**Discussion**

**Age interpretations.** Hereby we discuss the interpretation of the maximum ages derived from dating of four sedimentary samples: CA-CGT (dated in Marra et al.$^{21}$), CE-1, CE-2, BL-5 (dated in the present study).

*Sample CA-CGT.* 24 crystals extracted from the sand matrix of a ca. 2 m thick gravel layer cropping out at Colle Avorone locality (Sample CA-CGT) were dated by Marra et al.$^{21}$. 3 youngest crystals provided a weighted mean age of 359.5 ± 6.5 ka (2σ uncertainty) (Fig. 6a). Other 21 crystals provided two groups of ages, the largest one (14 grains) ranging 359–450 ka, with other six crystals spanning 700–850 ka (Fig. 6a).

The younger group of ages matches the duration of the main eruptive phase occurred at the Volsci Volcanic Field (VVF) 424 ± 13–349.5 ± 5.0 ka$^{21}$.

Regarding the older crystal ages, it should be noted that five samples from phreatomagmatic deposits analyzed in Marra et al.$^{21}$ yielding weighted mean ages in the interval 761.5 ± 9.5 to 349.5 ± 5.0 ka evidenced significant dispersion towards old ages. These products were mainly sourced from isolated, monogenetic eruptive centers. Therefore, it is unlikely that the older age values might have derived from earlier buried volcanic edifices or erupted deposits, and yet may provide evidence of older magma batches that cooled in sub-surface conditions.

*Samples CE-1, CE-2.* Sample CE-1 was collected in borehole core Ceprano 1 at 39.3 m depth within a coarse gravel layer with abundant sand matrix. The youngest crystal out of a population of 30 extracted from this sediment yielded a $^{40}$Ar/$^{39}$Ar age of 452.4 ± 1.8 ka (Fig. 6b). Sample CE-2 was collected at 15.1 m depth in borehole Ceprano 2, at the base of a coarse sand layer and yielded a youngest crystal date of 389.6 ± 3.7 ka (Fig. 6c).

These two maximum ages can be regarded as statistically significant even if based on one single crystal. The Ceprano boreholes were drilled in Campogrande which is located on the left hydrographic side of the Sacco catchment basin, a sector draining the most active and densely vent-populated volcanic area of the Volsci Volcanic Field. A climactic eruptive phase occurred at the VVF in the interval 420–350 ka$^{21}$, so the lack of crystals younger than 453 ka strongly suggests that the emplacement of the sand deposit occurred before the start of this volcanic phase. Consistent with this hypothesis, there is one crystal of 428 ± 10 ka along with one youngest crystal
of 390 ± 3.6 ka in the sample stratigraphically above. Moreover, these two ages along with that of 350.8 ± 8 ka on the primary layer occurring at the top of the sedimentary succession recovered in the Ceprano boreholes³³ provide a constant sedimentation rate of 38 cm/ky (see Fig. 4).

Therefore, it is reasonable to assume that likewise, the maximum ages derived from reworked sanidine crystals can be regarded as providing precise time constraints to sediment deposition, as the one on the primary volcanic layer.

Indeed, an age close to 453 ka for the gravel deposition during MIS 11 is in good agreement with the constraints provided from the Paleo-Tiber aggradational successions, which bracket it between 451 ± 2 and 445 ± 3 ka¹⁰,¹³.

**Sample BL-5.** Sample BL-5 was collected in a sandy-clayey travertine layer embedding several sub-cm sized, very altered volcanic scoriae. Crystals extracted from this sample yielded two distinct groups of ages ("a" and "b" in Fig. 6d).

One group of older ages ranging 300–600 ka ("a") is consistent with expected distribution for a deposit reworking the volcanic deposit of the Volsci Volcanic Field (VVF), the activity of which broadly spanned the interval 750–250 ka, as also observed in the sedimentary samples CE-1 and CE-2 (see Fig. 6b,c). The occurrence of one second group of ages ranging 120–30 ka ("b") is problematic to explain, given the much unlikely circumstance that it may reflect the real age of the deposit. Indeed, the travertine layer occurs on top of the "lower lacustrine succession", closely constrained at this location by tephra ages of 538 and 517 ka, and is part of the "upper fluvial-lacustrine succession" which previous¹⁹ and the present study constrain between 390 and ca. 300 ka.

Therefore, two hypotheses can be made to explain the second group of crystal ages. The anomalously young ages most probably represent much younger ash fall and/or eolian material that was worked into the deposit through fractures; indeed, the sampled deposit was exposed on a slightly inclined surface in the middle of a hillside, which is affected by continuous sliding of reworked sediments from the uphill section. Alternatively, these anomalous young ages may result from contamination in the field or laboratory. In any case, we consider the weighted mean age of 300 ± 12 ka yielded by the two youngest crystals in group "a" as a reliable maximum age for the deposit.

**Gravel origin and mechanisms of transport and deposition.** The Liri basin is a lacustrine basin of tectonic origin²⁴,²⁵, similar to other intra-mountain basins in the central Apennines³⁹, like the Fucino⁴⁰, L'Aquila⁴¹, Sulmona⁴², and Anagni⁴³ basins, characterized by deposition of thick successions of carbonate muds⁴³. The occurrence of discrete gravel beds in the Sacco-Liri basin were interpreted in the literature (e.g., Devoto¹⁹) as the evidence of intervening fluvial conditions. Our study shows that conglomerates represent episodes of high-energy water transport superposing the lacustrine conditions. The correlation of the investigated cross-sections shows that there are only two large gravel beds in the basin, which can be correlated laterally for several ten kilometers, intercalated with lacustrine muds. On the contrary, local mechanisms for gravel deposition, like occasional damming, may cause deposition of very limited fans. Moreover, in order to have gravel deposition as a consequence of damming, there must be strong capacity of transport within the catchment basin, which is not the case for a lacustrine basin such as the Sacco-Liri basin.

In contrast, during deglacial periods there is a strong release of gravel from the fans at the foot of the carbonate reliefs that border the upper portion of the catchment basin. This gravel is mobilized by the strongly increased water capacity of transport due to the increase in rainfall during the glacial termination, at the end of the cold-arid period leading into the warm-humid period.

Such processes have been documented and dated through the aggradational successions of the Tiber River catchment basin, which responded to variation in sea level. In contrast, an intra-mountain lacustrine basin is not expected to respond to the variation in base level, directly. However, during sea-level low stands an increased gradient in the lower portion of the basin (between the sea and the Liri basin) may cause erosion of the lacustrine deposits, with a shift into fluvial conditions which reverberates in the higher portion of the catchment basin, also causing an increase in the gradient. This adds increase in water capacity of transport allowing gravel eroded from the mountain slopes to be deposited in the Liri basin during the glacial termination at the beginning of the sea-level rise, when the baselevel is still low.

⁴⁰Ar/³⁹Ar ages provided here to constrain the gravel deposits demonstrate that their emplacement is coeval with global deglacial episodes, as reported in the δ¹⁸O curve⁴⁴.

**Morphostructural analysis of the Liri basin.** Figure 7A shows a NNW-SSE cross-section along the Sacco-Liri valley providing the chrono-stratigraphic correlation of the investigated sections, while a tentative reconstruction of the structural and sedimentary evolution of the Liri lacustrine basin is provided in Fig. 8.

Correlation of the sedimentary phases with the stack of globally distributed δ¹⁸O records⁴⁴ and the Red Sea relative sea-level (RSL) curve⁴⁵ is shown in Fig. 7b.

The occurrence of a Lower Pleistocene succession in Campo del Conte at higher elevation with respect to the Middle Pleistocene lacustrine succession in the Liri Basin provides a paleogeographic constraint to the tectonically subsiding sector (inferred fault A in Figs. 2 and 8a). Consistent with the present morphology of the Sacco Valley, confined within a narrow incision between the Meso-Cenozoic sedimentary and Middle Pleistocene volcanic terrains (dark green bordered sector in Fig. 2), no significant lacustrine sediments occur in this area, showing the absence of tectonic subsidence. In contrast, a more than 50 m-thick sedimentary succession spanning 450–337 ka (middle and upper lacustrine and fluvial successions) occurs in the sector between Cava Pompi and Lademagne, broadly corresponding to the Liri lacustrine basin (light green bordered sector in Fig. 2), and
<figure><img src="image_9.png" /><figcaption>**Figure 7.** (a) Schematic cross-section along the Liri basin showing the chronostratigraphic correlation of the investigated sections. See text for comments and explanation. (b) The age constraints to the sedimentary filling of the Liri basin provided by volcanic layers dated in this and in previous works (vertical red bars; shaded boxes are the $2\sigma$ uncertainties) allow comparison of sediment aggradation with the stack of $\delta^{18}\text{O}$ records$^{44}$ and the relative sea level (RSL) curve$^{45}$.</figcaption></figure>

overlies more lacustrine sediments (lower lacustrine succession), with ages spanning at least the 583–515 ka interval (Fig. 7a).

Extensive stratigraphic investigations conducted for the present work, summarized in the cross-section of Fig. 7, showed that the presence of coarse gravel beds cannot be related to local factors but concentrate in two distinct stratigraphic horizons. The uppermost of this horizon has average thickness of ca. 3 m over a large areal extent, covering at least 10 km, between Campo Grande and Lademagne (see Fig. 2), where it crops out at constant elevation between 110 and 115 m a.s.l. The lower gravel bed occurs at lower elevation, between 70 and 75 m a.s.l. in the Ceprano boreholes and is intercalated within the clayey sediments of the "Typical lacustrine succession" as far as Pontecorvo$^{19}$.

Taking into account all the above-mentioned chrono-stratigraphic constraints, we reconstructed a simplified structural evolution for the Liri Basin (summarized in Fig. 8).

A lower age constraint to the initiation of the tectonic subsidence of the Liri Basin is not available because the total thickness of the lower lacustrine succession is not known.

In Fig. 8a,b, we have represented only the portion of the lower lacustrine succession which has chronostratigraphic constraints, between 583 and 515 ka. Given the peculiar sedimentologic features of the lacustrine sediments, represented by carbonate-rich muds which are characteristic of a shallow water environment$^{19}$, this succession should have deposited at relatively constant elevation throughout the Liri basin (Fig. 8a). There is no reason to exclude that deposition continued until 500 ka (MIS 13.3 highstand, Fig. 7b) and that the upper portion of the succession was successively eroded during the MIS 12 regressive phase (Fig. 8b), as evidenced by the stratigraphy of the Pignataro Interamna section and by the age constraint of $452.4 \pm 1.8$ ka (sample CE-1) from the lower gravel horizon in the Ceprano boreholes (GH1 in Fig. 7a).

A close relation between the occurrence of MIS 12 lowstand and the deposition of the gravel layer GH1 at the base of the middle lacustrine-fluvial succession is evidenced by the age of sample CE-1 (Fig. 7b). $^{40}\text{Ar}/^{39}\text{Ar}$ constraints to the upper portion of this succession show that rapid aggradation of lacustrine to fluvial clay sediments occurred since the Glacial Termination V, which has an established age of 424 ka$^{44}$, through 402 ka in response to MIS 11 sea-level rise. The increased sedimentation rate based on the $^{40}\text{Ar}/^{39}\text{Ar}$ age constraints at the Ceprano boreholes is suggestive of enhanced tectonic subsidence in coincidence with the climactic eruptive
<figure><img src="image_10.png" /><figcaption>**Figure 8.** Proposed stepwise reconstruction of the sedimentary and tectonic evolution of the Liri basin through time. See text for comments and explanations.</figcaption></figure>

phase occurred at the Volsci Volcanic Field from 424 through 350 ka²¹. Starting from 390 ka, an increase in water transport energy within this sector of the Liri catchment basin resulted in an abrupt sedimentary switch, leading to the deposition of coarse sand with intercalated fine gravel sediments (Fig. 8d; upper fluvial-lacustrine succession). It is followed by the deposition of a 2–3 m thick horizon of well-sorted, coarse gravel, narrowly constrained in the interval 350.6 ± 8.0–345.9 ± 4.3 ka. This time interval encompasses the entire regressive phase of MIS 10 (Fig. 8e). Another abrupt sedimentary shift into clayey sediments is noted at Colle Avarone and Cava
Pompi sections, where lacustrine deposits overlying the coarse gravel horizons (Fig. 8f). However, thickness of the fine-grained portion of the upper fluvial-lacustrine succession is very limited and it passes upwards to wide travertine plateaus characterized by very thick deposits to the southeast¹⁹ (Fig. 8g).

We correlate this sedimentary change to the decrease/cessation of the tectonic subsidence in the Liri Basin, which may be due to the vanishing of the volcanic activity at the Volsci Volcanic Field (VVF) since 350 ka. Published ages of the VVF spanning 300–200 ka are poorly constrained, with the youngest reliable eruption age that of the Colle Borrello center occurring at 331.6 ± 3.0 ka²¹. The lack of a fine-grained aggradational succession deposited during MIS 9 highstand (337–325 ka, Fig. 7b) is likely related to the tectonic inversion within the Liri Basin.

Evidence for deep erosion of the upper fluvial-lacustrine succession occurs at San Giorgio al Liri section, where travertine deposits directly overlying the lower lacustrine succession yielded a *terminus post- quem* age of 300 ± 12 ka, consistent with occurrence of the regressive phase of MIS 8 since 320 ka (Fig. 7b). Contextual erosion and travertine progradational deposition is likely to have occurred throughout MIS 8, 320–270 ka (Fig. 7b), while the origin of a series of terraced paleo-surfaces throughout the Liri Basin is reflecting the subsequent regional uplift phase. Indeed, according to several authors¹⁶⁻⁵⁰, the Tyrrhenian Sea Margin of central Italy underwent variable uplift of several tens of meters in the last 250 ka.

**Relationships between sediment grainsize and sea-level fluctuations.** The ⁴⁰Ar/³⁹Ar constraints provided here to the deposits of the middle lacustrine succession bracket their deposition in the interval 452.4 ± 1.8 to 400.3 ± 3.0 ka, demonstrating that they correspond to an aggradational succession sensu Marra et al.⁵,⁷. In other words, the dated morpho-sedimentary unit formed in response to sea-level rise during the glacial-interglacial transition.

The complete chronostratigraphic record of the MIS 11 aggradational succession occurring in the coastal setting (San Paolo Formation²,⁷,¹⁰,¹³,¹⁴) reported in Fig. 8a shows that the equivalent deposit of the Liri Basin provides indistinguishable age constraints to the deposition of the basal gravel layer and to the completion of aggradation of the fine-grained sedimentary package, between 450 and 400 ka.

Remarkably, Giaccio et al.¹⁴ demonstrated the early aggradation of a first gravel layer at 443.1 ± 3.2 ka before the onset of Glacial Termination V and coinciding with a minor sea-level rise on the RSL curve⁴⁵ (Fig. 9a). Following this event, another gravel layer passing upwards to a thick clay succession was deposited after 437.1 ± 1.2 ka, in good agreement with the timing of the glacial termination and the sea-level rise at the onset of MIS 11 highstand (Fig. 9a). Giaccio et al.¹⁴ interpreted the early aggradation phase during MIS 12 as a first meltwater pulse (MWP1) event, preceding the larger amplitude meltwater pulse during Glacial Termination V (MPW2). These authors remarked that both MWPs coincide with episodes of ice-rafted debris deposition in the North Atlantic (Heinrich-like events⁵¹) and with attendant Southern Hemisphere warming, plausibly associated with the bipolar seesaw. Indeed, the occurrence of prominent ice-rafted debris (IRD¹⁸) peaks associated with these meltwater pulses alludes to episodes of extensive iceberg calving in the North Atlantic, consistent with sustained melting of the circum-North Atlantic ice sheets⁵² during these events.

The age of 452.4 ± 1.8 ka, which constrains the start of gravel deposition in the Liri Basin, provides strong indication of coarse clastic input to the river catchments of the Tyrrhenian Sea margin deriving from early deglaciation in the Apennine mountain range. It also provides further evidence for the validity of the sedimentary model of the "aggradational successions"⁵,⁷.

The reasons for such a far-field response also rely on the lesser elevation gain within the catchment basin, which has been successively increased by regional uplift over the last 250 ka⁴⁶⁻⁵⁰. In contrast, a comparison of the average uplift rate in the last 250 ka (0.24 mm/year) with the average sedimentation rate during the aggradational phases (e.g., 2.3 mm/year) suggests that glacio-eustasy overrides the tectonic effects, which only impact the accommodation space, and, in turn, the total thickness (rather than the timing of deposition) of each aggradational succession⁹.

The early aggradation of the MIS 9 sedimentary record (Aurelia Formation²) has been presented⁷ and has been constrained further by successive work⁹. In particular, a primary volcanic deposit dated at 345 ± 3 ka at the top of the basal gravel bed of the MIS 9 aggradational succession in the upper sector of the Tiber River catchment basin⁹ (Fig. 9b) provided a *terminus ante- quem* for its deposition, preceding the canonical age of 337 ka⁴⁴ for Glacial Termination IV. The narrow age constraints to the second gravel layer occurring in the upper lacustrine succession in the Liri Basin provide a striking match with those provided in the coastal and more inner sectors of the Tiber River basin²,⁷,⁹, evidencing an early deglaciation occurring during MIS 10 (Fig. 9b). Remarkably, this early aggradational phase also corresponds with a minor sea-level rise in the RSL curve preceding the glacial termination (Fig. 9b), suggesting the same triggering mechanism (i.e., early ice melting) as that hypothesized for the analogous eustatic event during MIS 12.

**Possible triggering mechanisms to deglaciation.** The occurrence of meltwater pulses/events during the glacial maxima that preceded both T-V and T-IV is challenging since the two glacial-interglacial cycles have very different glacial histories and insolation forcing (e.g., Spratt and Lisiecki⁵³). Moreover, T-V was characterized by slower rates of ice-sheet melting/sea-level rise than for most of the last five terminations⁴⁵.

While an in-depth analysis of the possible causes is beyond the scope of this paper, we remark on an intriguing aspect of the insolation curve during these glacial termination which was already pointed out as a possible contributing factor to trigger the deglacial process.

Associated with the introduction of the aggradational successions model, Marra et al.⁷ proposed a possible forcing mechanism based on the occurrence of particularly mild insolation minima which may be regarded as the pre-conditioning factor to trigger a glacial termination. We re-propose this notion (Fig. 10) based on the
<figure><img src="image_11.png" /><figcaption>**Figure 9.** $^{40}$Ar/$^{39}$Ar age constraint to sediment aggradation in the Liri basin (red vertical bars) for (a) MISs 11–12 across Glacial Termination V, and (b) MISs 9–10 across Glacial Termination IV are compared to those provided in previous works to the aggradational successions of the Paleo-Tiber River reported in the stratigraphic sections to the right. Red arrows pointing to the left indicate terminus post-quem (maximum) ages and those pointing to the right indicate terminus ante-quem (minimum) ages. See text for comments and explanation.</figcaption></figure>

comparison of the two geochronologically constrained sedimentary records of glacial terminations V and IV provided in this paper with the mean summer insolation curve at 65° N⁵⁴.

The interval 450–359 ka encompassing the two glacial terminations investigated in this paper is characterized by the occurrence of four such "mild" minima. Indeed, if we consider the insolation value of 362 W/m² for the insolation minimum occurring at 360 ka as a lower threshold (blue horizontal line in Fig. 10), we see that only six other minima were characterized by higher insolation, whereas all other minima in the last 800 ky were "colder". These six mild (warmer) minima are all associated with evidence of anomalous meltwater pulse events. The oldest one (369 W/m²), at ca. 740 ka, is associated with the early emplacement at ca. 750 ka of a thick gravel layer in the Paleo-Tiber delta, which was interpreted as a previously unrecognized glacial termination (i.e., MIS 18.3, TVIII-A⁴,⁵). Remarkably, the second oldest mild minimum at 543 ka is associated with an identical feature of the $\delta^{18}$O curve: a double isotopic peak. However, the only difference in this case is that the earliest peak (13.3) is associated in the literature with the eustatic event considered to be the true glacial termination VI, as opposed to MIS 17, for which the glacial termination has been associated with the younger 17.3 peak⁴⁴ (Fig. 10). Unfortunately, the lack of a detailed record of RSL in the interval > 500 ka hinders the possibility to verify if both these mild minima actually preceded meltwater pulses.
<figure><img src="image_12.png" /><figcaption>**Figure 10.** $^{40}$Ar/$^{39}$Ar age constraints (vertical red bars, shaded boxes are the 2σ uncertainties) to the coarse gravel beds occurring in the sedimentary filling of the Liri basin provided in this work, as well as at the base of the aggradational successions of the Paleo-Tiber River provided in previous studies$^{7,8,11}$. Red arrows as in Fig. 8. (a) Comparison with the glacial termination (GT) ages established through the calibration of the $\delta^{18}$O curve$^{44}$; (b) comparison with a set of "mild" minima (see text for explanation) of the mean summer insolation curve at 65°N$^{54}$; (c) comparison with the ice-rafted debris (IRD) record of core ODP 983$^{51}$ (re-drawn here for comparative purpose).</figcaption></figure>

In contrast, a striking coincidence between the mild minima at 452 and 360 ka and the early aggradational phases during MIS 12 and MIS 10, which are associated with the emplacement of the coarse gravel beds at ~453 and ~351 ka in the catchment basins of the Tiber and of the Sacco-Liri rivers, is evident by inspection of Fig. 10. A lesser of such events is also clearly associated with the mild minimum of 371 W/m$^2$ occurring at 380 ka, which triggered the grainsize inversion in the Liri Basin causing the emplacement of the coarse sand-and fine gravel horizon which has a *terminus post- quem* age of 390 ± 1 ka.

Moreover, the most outstanding observation is that the only other mild minima of 366 W/m$^2$ in the last 350 ka occurred at 40 ka, in close connection with the onset of the Heinrich events$^{51}$.

This is strong supporting evidence for relating the early aggradational phases responsible for the emplacement of gravel beds in the Tiber and Liri basins with deglaciation events (meltwater pulses) that triggered relatively minor sea-level rises, and which represent the equivalent of past Heinrich-like events in the Apennines of central Italy.

## Conclusions

$^{40}$Ar/$^{39}$Ar geochronology used to constrain the aggradational phases and grainsize variations of the sedimentary filling of the Liri fluvial-lacustrine basin provide strong evidence to support the "aggradational successions" model$^{5,7,14}$ as a powerful tool to detect the occurrence of deglaciation events that triggered global meltwater pulses.

We demonstrate a substantial synchronicity between the ages of gravel deposition in both the Tiber and Liri rivers catchment basins and the occurrence of moderate sea-level rise events, which anticipate those more marked during the glacial terminations V and IV in the Red Sea relative sea level curve$^{45}$.
We also show a striking correspondence among the occurrence of particularly mild (warmer) minima of the mean summer insolation at 65° N⁵⁴ and these early aggradational phases, as well as with other anomalous early sea-level rises occurring at the onset of glacial terminations VIII and VI, and at 40 ka at the onset of the so-called Heinrich events⁵¹.

Such patterns suggest that gravel deposition is triggered by melting of the Apennines mountain range glaciers, providing the water transport energy and a surplus of elastic input in the catchment basins of the rivers draining the mountain regions and flowing into the Tyrrhenian Sea. Such hydrologic/sedimentary processes, which occur in any mountain region of the globe, if correctly depicted and dated, can provide a large dataset of deglaciation proxies to unravel the chronology of glacio-eustatic events occurring in the last 1 Ma.

## Methods

⁴⁰Ar/³⁹Ar analysis. Samples for ⁴⁰Ar/³⁹Ar analyses were prepared at the Laboratoire des Sciences du Climat et de l'Environnement facility (CNRS-CEA, Gif-sur-Yvette), France, and at the University of Wisconsin-Madison.

Three distinct irradiations have been performed and the samples were dated in three facilities (Berkeley Geochronology Center, USA), Laboratoire des Sciences du Climat et de l'Environnement (CEA, Gif-sur-Yvette), and WiscAr Laboratory of Wisconsin University (USA). Samples PO-C6, BL-1A, BL-5, PI-1 and PI-2 were irradiated in the Cd-lined, in-core CLICIT facility of the Oregon State University TRIGA reactor.

Samples PO-C6, BL-1A and BL-5 were analyzed at the Berkeley Geochronology Center (BGC; California, USA), using a MAP 215-C mass spectrometer (MAP 1), following procedures described in Giaccio et al.⁵⁵. Samples BL-4, CE-1, and CE-2 were analyzed at the University of Wisconsin-Madison (USA), using a Noblesse 5-collector mass spectrometer, following procedures described in Jicha et al.⁵⁶. Samples PI-1 and PI-2 were analyzed at LSCE using a VG 5400 mass spectrometer (LSCE; Gif-sur-Yvette, France), following procedures described in Pereira et al.³⁵.

All ages are calculated according to the fluence monitor age of Alder Creek sanidine (40 K total decay constant of Min et al.⁵⁷; $ACs = 1.1848 \pm 0.0012 Ma^{29,30}$) and are reported to the precision level of 2σ standard deviation. Full ⁴⁰Ar/³⁹Ar data are reported in Supplementary File #2.

## Detrital sanidine dating approach

The implemented sensitivity of the modern mass-spectrometers permits to date with great precision Pleistocene grains smaller than 400–300 µm. Combined with the continuous eruptive activity that characterized the volcanic region of central Italy during the last 800 ka⁵⁸, and ref. therein, dating of sedimentary samples has become an extremely useful tool to assess the ages of aggradational successions deposited in response to sea-level rise during glacial terminations in the absence of intercalated, primary tephra layers (see Marra et al.⁵⁰ for an in-depth discussion). In fact, when a statistically significant number of crystals is dated (i.e., 30–40 grains), it is reasonable to assume that the age of the youngest crystal population, besides providing a maximum age for the sedimentary deposit, should also be regarded as documenting the lack of crystals from younger eruptive products, implying that no younger eruptions occurred before the time of deposition. Such assumption allows to consider the youngest crystal age also an approximate minimum age (*terminus ante- quem*) (within the recurrence time of the volcanic activity) to the time of deposition of the sediment. As discussed in Marra et al.⁵⁰ the youngest eruptions should be better represented in reworked, sedimentary deposits because their products crop out in wider areas than the older ones, which are buried under a longer sequence of strata. This consideration supports the principle that the age of a layer is bracketed between the ages of its youngest crystal population and of the next younger eruption, whose crystals do not occur in the layer but is documented in the area.

## Data availability

All data generated or analysed during this study are included in this published article [and its supplementary information files].

Received: 14 May 2022; Accepted: 31 October 2022

Published online: 05 November 2022

## References

1. Alvarez, W. et al. Quaternary fluvial-volcanic stratigraphy and geochronology of the Capitoline hill in Rome. *Geology* **24**, 751–754 (1996).

2. Karner, D. B. & Marra, F. Correlation of fluviodeltaic aggradational sections with glacial climate history: A revision of the Pleistocene stratigraphy of Rome. *Geol. Soc. Am. Bull.* **110**, 748–758 (1998).

3. Karner, D. B. & Renne, P. R. 40Ar/39Ar geochronology of Roman volcanic province tephra in the Tiber River valley: Age calibration of middle Pleistocene sea-level changes. *Geol. Soc. Am. Bull.* **110**, 740–747 (1998).

4. Marra, F., Florindo, F. & Karner, D. B. Paleomagnetism and geochronology of early Middle Pleistocene depositional sequences near Rome: Comparison with the deep sea $\delta^{18}O$ climate record. *Earth Planet. Sci. Lett.* **159**, 147–164 (1998).

5. Marra, F., Florindo, F. & Boschi, E. History of glacial terminations from the Tiber River, Rome: Insights into glacial forcing mechanisms. *Paleoceanography* **23**, 1–17. https://doi.org/10.1029/2007PA001543 (2008).

6. Marra, F., Bozzano, F. & Cinti, F. R. Chronostratigraphic and lithologic features of the Tiber River sediments (Rome, Italy): Implications on the Post-glacial sea-level rise and Holocene climate. *Glob. Planet. Change* https://doi.org/10.1016/j.gloplacha.2013.05.002 (2013).

7. Marra, F. et al. Independent 40Ar/39Ar and 14C age constraints on the last five glacial terminations from the aggradational successions of the Tiber River, Rome (Italy). *Earth Planet Sci. Lett.* **449**, 105–117. https://doi.org/10.1016/j.epsl.2016.05.037 (2016).

8. Marra, F., Jicha, B. & Florindo, F. 40Ar/39Ar dating of Glacial Termination VI: Constraints to the duration of Marine Isotopic Stage 13. *Sci. Rep.* **7**, 8908. https://doi.org/10.1038/s41598-017-08614-6 (2017).
9. Marra, F. et al. Combined glacio-eustatic forcing and volcano-tectonic uplift: Geomorphological and geochronological constraints on the Tiber River terraces in the eastern Vulsini Volcanic District (central Italy). *Global Planet. Change* **182**, 103009. https://doi.org/10.1016/j.gloplacha.2019.103009 (2019).

10. Marra, F., Pereira, A., Boschian, G., & Nomade, S., 2021. MIS 13 and MIS 11 aggradational successions of the Paleo-Tiber delta: Geochronological constraints to sea-level fluctuations and to the Acheulean sites of Castel di Guido and Malagrotta (Rome, Italy), Quaternary International.

11. Marra, F. & Florindo, F. The subsurface geology of Rome: Sedimentary processes, sea-level changes and astronomical forcing. *Earth Sci. Rev.* **136**, 1–20 (2014).

12. Florindo, F. et al. Radioisotopic age constraints for Glacial Terminations IX and VII from aggradational sections of the Tiber River delta in Rome, Italy. *Earth Planet. Sci. Lett.* **256**, 61–80. https://doi.org/10.1016/j.epsl.2007.01.014 (2007).

13. Pereira, A. et al. Tephrochronology of the central Mediterranean MIS 11c interglacial (~425 to 395 ka): New constraints from Vico volcano and Tiber delta, Central Italy. *Quatern. Sci. Rev.* **243**, 106470 (2020).

14. Giaccio, B. et al. Tephrochronological constraints on the timing and nature of sea-level change prior to and during glacial termination V. *Quatern. Sci. Rev.* https://doi.org/10.1016/j.quascirev.2021.106976 (2021).

15. Marra, F., Pereira, A., Boschian, G. & Nomade, S. MIS 13 and MIS 11 aggradational successions of the Paleo-Tiber delta: Geochronological constraints to sea-level fluctuations and to the Acheulean sites of Castel di Guido and Malagrotta (Rome, Italy). *Quatern. Int.* https://doi.org/10.1016/j.quaint.2021.12.016 (2021).

16. Giraudi, C. & Giaccio, B. Middle Pleistocene glaciations in the Apennines, Italy: New chronological data and preservation of the glacial record. *Geol. Soc. Lond. Spec. Publ.* **433**, 161–178. https://doi.org/10.1144/SP433.1 (2017).

17. Luberti, G. M., Marra, F. & Florindo, F. A review of the stratigraphy of Rome (Italy) according to geochronologically and paleomagnetically constrained aggradational successions, glacio-eustatic forcing and volcano-tectonic processes. *Quatern. Int.* **438**, 40–67. https://doi.org/10.1016/j.quaint.2017.01.044 (2017).

18. Barker, S. et al. Early interglacial legacy of deglacial climate instability. *Paleoceanogr. Paleoclimatol. 34*, 1455–1475. https://doi.org/10.1029/2019PA003661 (2019).

19. Devoto, G. Lacustrine Pleistocene in the lower Liri Valley (southern Latium). *Geol. Romana* **4**, 291–368 (1965).

20. Centamore, E. et al. Note illustrative del Foglio 402 Ceccano. *Carta Geologica d’Italia* 1:50,000 (Servizio Geologico d’Italia, 2010).

21. Marra, F. et al. The Volsci Volcanic Field (central Italy): An open window on continental subduction processes. *Int. J. Earth Sci.* **110**, 689–718. https://doi.org/10.1007/s00531-021-01981-6 (2021).

22. Muttoni, G., Scardia, G., Kent, D. V. D., Swisher, C. C. & Swanzi, G. Pleistocene magnetochronology of early hominin sites at Ceprano and Fontana Ranuccio, Italy. *Earth Planet. Sci. Lett.* **286**, 255–268 (2009).

23. Malinverno, A. & Ryan, W. B. F. Extension in the Tyrrhenian sea and shortening in the Apennines as results of arc migration driven by sinking of the lithosphere. *Tectonics* **5**, 227–245 (1986).

24. Patacca, E., Sartori, R. & Scandone, P. Tyrrhenian basin and apenninic arcs: Kinematic relations since late Tortonian times. *Mem. Soc. Geol. Ital.* **45**, 425–451 (1990).

25. Centamore, E., Di Manna, P. & Rossi, D. Kinematic evolution of the Volsci Range: A new overview. *Ital. J. Geosci.* **126**, 159–172 (2007).

26. Cardello, G. L. et al. Tectonically controlled carbonate-seated maar-diatreme volcanoes: The case of the Volsci Volcanic Field, central Italy. *J. Geodyn.* **139**, 101763 (2020).

27. Sani, F. et al. Tectonic evolution of the internal sector of the Central Apennines, Italy. *Mar. Pet. Geol.* **21**(10), 1235–1254 (2004).

28. Acocella, V., Facenna, C. & Funiciello, R. Elementi strutturali della media Valle Latina. *Boll. Soc. Geol. Ital.* **115**, 501–518 (1996).

29. Niespolo, E. M., Rutte, D., Deino, A. L. & Renne, P. R. Intercalibration and age of the Alder Creek sanidine 40Ar/39Ar standard. *Quat. Geochronol.* **39**, 205–213 (2017).

30. Nomade, S., Gauthier, A., Guillou, H. & Pastre, J. F. 40Ar/39Ar temporal framework for the Alleret maar lacustrine sequence (French Massif-Central): Volcanological and paleoclimatic implications. *Quatern. Geochronol.* **5**, 20–27 (2010).

31. Palombo, M. R., Magri, D., Molinari, D., & Pisano, V. 2000–2002. The Pleistocene sequence of Campo del Conte (Lower Sacco Valley, southern Lazio). *Geologica Romana*, 36, 289–309.

32. Gliozzi, E. et al. Biochronology of selected mammals, molluscs and ostracods from the middle Pliocene to the late Pleistocene in Italy. *Riv. Ital. Paleontol. Stratigr.* **103**, 369–388 (1997).

33. Nomade, S., Muttoni, G., Guillou, H., Robin, E. & Scardia, G. First $^{40}$Ar/$^{39}$Ar age of the Ceprano man (central Italy). *Quat. Geo-chronol.* **6**, 453–457 (2011).

34. Biddittu, I. & Segre, A. G. Paleolitico inferiore a Cava Pompi presso Pofi, Frosinone. *Quaderni Centro Studi Arch. Etrusco-Ital.* **1**, 78–79 (1978).

35. Pereira, A. et al. Geochronological evidences of a MIS 11 to MIS 10 age for several crucial Acheulian sites from the Frosinone province (Latium, Italy): Archaeological implications. *Quatern. Sci. Rev.* **187**, 112–129 (2018).

36. Biddittu, I., 2004. Guida del Museo Preistorico di Pofi. Quaderni di Argyl 1. Regione Lazio, Amministrazione provinciale di Frosinone, Comune di Pofi.

37. Biddittu, I., Canetri, E., Comerci, V., Germani, M. & Picchi, G. 2012. Nuove ricerche nel giacimento del Paleolitico inferiore di Lademagne, S. Giovanni Incarico (Frosinone). In *Lazio e Sabina 9* (eds Ghini, G. & Mari, Z.) 437–443 (Edizioni Quasar, 2012).

38. Narcisi, B. Ricerche di tefracronologia nella media e bassa Valle Latina. *Mem. Soc. Geol. Ital.* **35**, 909–912 (1986).

39. Centamore, E. & Rossi, D. Neogene-Quaternary tectonics and sedimentation in the Central Apennines. *Itali. J. Geosci.* **128**, 73–88 (2009).

40. Cavinato, G. P., Carusi, C., Dall'Asta, M., Miccadei, E., & Piacentini, T. Sedimentary and tectonic evolution of Plio-Pleistocene alluvial and lacustrine deposits of Fucino Basin (central Italy). *Sed. Geol.* **148**, 29–59. https://doi.org/10.1016/S0037-0738(01)00209-3 (2002).

41. Cosentino, D. et al. New insights into the onset and evolution of the central Apennine extensional intermontane basins based on the tectonically active L'Aquila Basin (central Italy). *Bull. Geol. Soc. Am.* **129**, 1314–1336. https://doi.org/10.1130/B31679.1 (2017).

42. Galli, P., Giaccio, B., Peronace, E. & Messina, P. Holocene Pale-oearthquakes and Early-Late Pleistocene slip rate on the Sulmona Fault (Central Apennines, Italy). *Bull. Seismol. Soc. Am.* **105**(1), 1–13 (2015).

43. Florindo, F. et al. Chronostratigraphy of the Anagni basin, central Italy. Environmental evolution, faunal and human occupation since 2 Ma. *Sci. Rep.* **11**, 7056. https://doi.org/10.1038/s41598-021-85446-5 (2021).

44. Lisiecki, L. E. & Raymo, M. E. A Pliocene–Pleistocene stack of 57 globally distributed benthic $\delta^{18}$O records. *Paleoceanography* **20**, PA1003. https://doi.org/10.1029/2004PA001071 (2005).

45. Grant, K. M. et al. Sea-level variability over five glacial cycles. *Nat. Commun.* **5**, 5076. https://doi.org/10.1038/ncomms6076 (2014).

46. Karner, D. B., Marra, F., Florindo, F. & Boschì, E. Pulsed uplift estimated from terrace elevations in the coast of Rome: Evidence for a new phase of volcanic activity?. *Earth Planet. Sci. Lett.* **188**, 135–148 (2001).

47. Ferranti, L. et al. Markers of the last interglacial sea-level high stand along the coast of Italy: Tectonic implications. *Quatern. Int.* **145**(146), 30–54 (2006).

48. Marra, F., Florindo, F., Anzidei, M. & Sepe, V. Paleo-surfaces of glacio-eustatically forced aggradational successions in the coastal area of Rome: Assessing interplay between tectonics and sea-level during the last ten interglacials. *Quatern. Sci. Rev.* **148**, 85–100. https://doi.org/10.1016/j.quascirev.2016.07.003 (2016).
49. Marra, F., Florindo, F. & Petronio, C. Quaternary fluvial terraces of the Tiber Valley: Geochronologic and geometric constraints on the back-arc magmatism-related uplift in central Italy. *J. Sci. Rep.* 7, 2517. https://doi.org/10.1038/s41598-017-02437-1 (2017).

50. Marra, F. et al. MIS 9 to MIS 5 terraces along the Tyrrenian Sea coast of Latium (central Italy): Assessing interplay between sea-level oscillations and tectonic movements. *Geomorphology* **346**, 106843. https://doi.org/10.1016/j geomorph.2019.106843 (2019).

51. Hemming, S. R. Heinrich events: Massive late Pleistocene detritus layers of the North Atlantic and their global climate imprint. *Rev. Geophys.* **42**, 1005. https://doi.org/10.1029/2003RG000128 (2004).

52. McManus, J. F., Oppo, D. W. & Cullen, J. L. A 0.5-million-year record of millennial-scale climate variability in the North Atlantic. *Science* **283**(5404), 971–975 (1999).

53. Spratt, R. M. & Lisiecki, L. E. A Late Pleistocene sea level stack. *Clim. Past* **12**, 1079–1092. https://doi.org/10.5194/cp-12-1079-2016 (2016).

54. Laskar, J. et al. A long-term numerical solution for the insolation quantities of the Earth. *Astron. Astrophys.* **428**, 261–285 (2004).

55. Giaccio, B. et al. First integrated tephrochronological record for the last ~190 kyr from the Fucino Quaternary lacustrine succession, central Italy. *Quatern. Sci. Rev.* **158**, 211–234. https://doi.org/10.1016/j.quascirev.2017.01.004 (2017).

56. Jicha, B. R., Singer, B. S. & Sobol, P. Re-evaluation of the ages of 40Ar/39Ar sanidine standards and supereruptions in the western US using a Noblesse multi-collector mass spectrometer. *Chem. Geol.* **431**, 54–66 (2016).

57. Min, K., Mundil, R., Renne, P. R. & Ludwig, K. R. A test for systematic errors in 40Ar/39Ar geochronology through comparison with U/Pb analysis of a 1.1-Ga rhyolite. *Geochim. Cosmochim. Acta* **64**, 73–98 (2000).

58. Marra, F. et al. Monti Sabatini and Colli Albani: The dormant twin volcanoes at the gates of Rome. *Sci. Rep.* **10**(1), 8666 (2020).

## Author contributions

F.M. designed the study, performed the field investigations, and wrote the paper; A.P. performed the ⁴⁰Ar/³⁹Ar analyses and field investigations; B.J.; performed the ⁴⁰Ar/³⁹Ar analyses; S.N. performed the ⁴⁰Ar/³⁹Ar analyses and field investigations; I.B. performed the field investigations; F.F. performed the field investigations; G.M. performed the borehole investigations in Campo Grande; E.M.N. performed the ⁴⁰Ar/³⁹Ar analyses; P.R.R. performed the ⁴⁰Ar/³⁹Ar analyses; V.S. performed the ⁴⁰Ar/³⁹Ar analyses. All authors contributed to the writing of the text, discussed the results of the study and commented on the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-23391-7.

**Correspondence** and requests for materials should be addressed to F.M.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_13.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022