# Environmental evolution, faunal and human occupation since 2 Ma in the Anagni basin, central Italy

Fabio Florindo, Fabrizio Marra, Diego E Angelucci, Italo Biddittu, Luciano Bruni, Federico Florindo, Mario Gaeta, Hervé Guillou, Brian Jicha, Patrizia Macrì, et al.

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

## To cite this version:

Fabio Florindo, Fabrizio Marra, Diego E Angelucci, Italo Biddittu, Luciano Bruni, et al. Environmen- tal evolution, faunal and human occupation since 2 Ma in the Anagni basin, central Italy. Scientific Reports, 2021, 11 (1), 10.1038/s41598-021-85446-5. http-03207638

HAL Id: hal-03207638

https://hal.science/hal-03207638

Submitted on 26 Apr 2021

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_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

Check for updates

OPEN

# Environmental evolution, faunal and human occupation since 2 Ma in the Anagni basin, central Italy

Fabio Florindo<fnref n="1" />,²✉, Fabrizio Marra<fnref n="1" />, Diego E. Angelucci³, Italo Biddittu⁴, Luciano Bruni⁴, Federico Florindo⁵, Mario Gaeta⁶, Hervé Guillou⁷, Brian Jicha⁸, Patrizia Macri<fnref n="1" />, Caterina Morigi⁹, Sebastien Nomade⁷, Fabio Parenti⁴,¹⁰, Alison Pereira¹¹,¹² & Stefano Grimaldi³,⁴

We present the study of a composite, yet continuous sedimentary succession covering the time interval spanning 2.6–0.36 Ma in the intramontane basin of Anagni (central Italy) through a dedicated borecore, field surveys, and the review of previous data at the three palaeontological and archaeological sites of Colle Marino, Coste San Giacomo and Fontana Ranuccio. By combining the magneto- and chronostratigraphic data with sedimentologic and biostratigraphic analysis, we describe the palaeogeographic and tectonic evolution of this region during this entire interval. In this time frame, starting from 0.8 Ma, the progressive shallowing and temporary emersion of the large lacustrine basins and alluvial plains created favorable conditions for early hominin occupation of the area, as attested by abundant tool industry occurrences and fossils. This study provides new constraints to better interpret the hominin migratory dynamics and the factors that influenced the location and spatial distribution during the early occupation of this region.

Since the late 1960s, scholars of the Italian Institute of Human Palaeontology (hereafter, the IIPU) have carried out field surveys, archaeological excavations, and geomorphological studies in the Latin Valley, Central Italy (Fig. 1). This area is renowned for the abundant palaeontological, archaeological, and palaeoanthropological occurrences, including the finding in 1994 of a calvarium at the Ceprano locality, initially considered as the oldest evidence of human presence in the Italian peninsula<fnref n="1" />⁻³. Despite the large number of reports, few chronologic constraints have been achieved for the faunal and lithic assemblages recovered at several sites in the Latin Valley until recently. The first integrated, modern approach to the study of these sites was performed at Fontana Ranuccio (FR) and at Coste San Giacomo (CSG), where an outstanding vertebrate faunal record of Villafranchian age was obtained⁴⁻¹⁰.

Through combined sedimentological, palaeomagnetic, palynological and palaeontological investigations, an age of ca. 2 Ma for CSG was suggested and the site was proposed as the type-locality for the homonymous Faunal Unit of the Middle Villafranchian⁴,⁵. More recently, the application of the ⁴⁰Ar/³⁹Ar dating method on the volcanic deposits intercalated within the sedimentary successions of the Latin Valley has led to provide precise ages to several important Middle Pleistocene archaeological sites¹¹,¹². Among these, the abovementioned Ceprano site, as well as FR, where one of the richest bone-tool assemblages of Italy was recovered⁶,¹⁰,¹³,¹⁴, along with four human teeth¹⁵⁻¹⁷, lithic assemblage¹⁸, and a large vertebrate bone fossil record, which constitutes the local fauna for the homonymous, Late Galerian Faunal Unit of the Aurelian Mammal Age¹⁹,²⁰.

Other localities in this region contain hominin remains attributed to *Homo heidelbergensis* such as those recovered during excavations at the Cava Pompi site of Pofi in 1961 and 1976 along with lithic and bone industries²¹,²². Similar assemblages also occur at Lademagne and Isoletta⁶,²³.

<footnote n="1">¹Istituto Nazionale di Geofisica e Vulcanologia, Rome, Italy. ²Institute for Climate Change Solutions, Via Sorchio snc, 61040 Frontone, Italy. ³Department of Humanities, University of Trento, Trento, Italy. ⁴Istituto Italiano di Paleontologia Umana, Anagni, Italy. ⁵Sapienza Università di Roma, Piazzale Aldo Moro 5, 00185 Roma, Italy. ⁶Dipartimento di Scienze della Terra, Sapienza Università di Roma, Piazzale Aldo Moro 5, 00185 Rome, Italy. ⁷LSCE/IPSL, UMR CEA-CNRS-UVSQ 8212, Laboratoire des Sciences du Climat et de l'Environnement, CEA Saclay, Bat 714, Chemin de Saint Aubin - RD 128, 91191 Gif sur Yvette, France. ⁸Department of Geoscience, University of Wisconsin-Madison, Madison, USA. ⁹Department of Earth Sciences, University of Pisa, Via S. Maria 53, 56126 Pisa, Italy. ¹⁰Universidade Federal do Paraná, Curitiba, Brazil. ¹¹CNRS Laboratoire GEOPS, Université Paris-Saclay, Orsay, France. ¹²Département Hommes et Environnements, Muséum national d'Histoire naturelle, Paris, France. ✉email: fabio.florindo@ingv.it</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 1.** (a) Digital Elevation Model (DEM) image TINITALY/01 square WA 6570, of central Italy showing location of the Latin Valley and (b) of the investigated archaeological sites. Used with permission of the Istituto Nazionale di Geofisica e Vulcanologia, Rome.</figcaption></figure>

A much younger age than previously thought was determined for “Ceprano man” based on a 353 ± 4 ka date of an ash layer interbedded in the lacustrine deposits hosting the human calvarium¹¹. An age of 408 ± 10 ka was suggested for the FR Acheulian lithic assemblage, revising previous estimations of ca. 458 ka⁶, along with similar ages around 400 ka for the Lademagne, Isoletta, and Cava Pompi sites by¹².

In contrast to the great improvements achieved on the chronologic constraints of the archaeological sites mentioned above, no new investigations have been delivered at Colle Marino (CM) since its discovery⁶, despite the potential for groundbreaking implications of the findings reported at this site, which is located two kilometers east of CSG and FR. Indeed, a lithic assemblage—comprising uni- and bifacial choppers, polyhedrons, flakes, and pebble scrapers—was found on surface⁶, and its resemblance with the Oldowan (Mode I) industry, the oldest-known stone tool industry, was also suggested²⁴. The potential archaeological layer of CM was correlated with a sedimentologically similar horizon, occurring at the same elevation, and immediately overlying the Middle Villafranchian fossiliferous layer of CSM, suggesting quite an ancient age (i.e., more than 0.7 Ma) for the lithic assemblage, consistent with the alleged Oldowan features⁶. Despite the presence of faunal remains of the ancient hyena *Pachycrocuta brevirostris* in the continental deposits of CM¹⁰, which correspond to a wide
time interval between 1.8 and 0.8 Ma, the lithic assemblage recovered at CM was never re-investigated and its puzzling occurrence has been met with skepticism²⁵ probably due to the lack of suitable chronological markers.

To further understand the CM archaeological horizon and provide reliable age constraints for it, the IIPU and the Istituto Nazionale di Geofisica e Vulcanologia (hereafter, INGV) set up a research project in this area. In August 2020, a dedicated borehole was drilled at this locality in addition to field investigations and a systematic review of the log data from previous boreholes performed at the CSM and FR archaeological sites (Fig. 1).

The aim of the coring was to recover the sedimentary successions in the elevation range and correlate them with the one from CSG and with the archaeological layer outcropping at CM, to carry out litho-, bio-, and magneto-stratigraphic analyses, combined with $^{40}$Ar/$^{39}$Ar dating of possible tephra layers, as well as of detrital sanidine (see²⁶ for this methodology).

**Geological setting.** The foothills region of the Central Apennines is located between an arc-shaped fold-and-thrust belt and the Tyrrhenian Sea margin. The region hosts several palaeontological and archaeological sites which yielded a record of the faunal and early human occupation during the last two million years. The distribution of such occurrences is direct consequence of the palaeogeographic evolution of this region, a back-arc extensional domain (e.g.,²⁷,²⁸), which underwent to a progressive continentalization due to the uplift of the Apennine chain and to the southwest migration of the Tyrrhenian Sea basin and its coastline. Due to the temporal and spatial variations of the continental environments, the outcrops of sedimentary deposits featuring the oldest terrestrial faunal records, which are Middle to Late Villafranchian in age, are scattered within a narrow belt at the western foot of the Apennine range. This is an elongated area parallel to the coastline of the maximum marine ingression culminated during the Santernian (i.e., the local marine stage of the Lower Pleistocene spanning 1.8–1.5 Ma) (Fig. 2). The Latin Valley, a tectonic depression subjected to an extensional regime since the late Pliocene that led to the formation of lacustrine basins (e.g. Anagni and Ceprano), is located at the southern margin of this area²⁹⁻³¹.

**Colle Marino site location and lithology.** Drilling at CM was carried out in the same area (Fig. 1) where the Oldowan-like lithic assemblage⁶ tentatively dated to an age > 0.7 Ma was found. An 8 cm diameter core (from here on, CM1 core) was taken (41°45′32.30″N, 13°07′13.80″E, 227 m elevation; exact location of the drilling site is reported in Fig. S1 of Supplementary Material #1) by standard rotary drilling techniques to a depth of 35 m and the recovered sections lack an azimuthal orientation. The recovered log (Fig. 3; see Fig. S2 for detailed stratigraphy and sampling) includes: the topsoil, from 0 to 0.45 m below the ground surface (bgs); sand with carbonate concretions and travertine layers from 0.45 to 5.0 m bgs, which corresponds to the archaeological horizon⁶; loose silty sand from 5 to 6.98 m bgs, a buried soil from 6.98 to 7 m bgs and compact grey silty- to sandy-clays from 7 to 35 m bgs. Photographs of the cored sediment in the interval 0–9 m bgs are shown in Fig. S3.

In the area surrounding the drilling site⁶, reported the occurrence of a volcanic succession consisting of a lithified tuff at the base, overlain by a several meter-thick pozzolane rich pyroclastic-flow deposit, directly above the sand and travertine horizon. Field work performed for the present study allowed us to recognize the entire volcanic succession as part of the Pozzolane Rosse pyroclastic-flow deposit (456 ± 4 ka³²,³³).

**Field investigations.** The investigated site of CM is about 1.5 km from FR and 2 km from CSG (Fig. 1). Chronostratigraphy of the FR site is here reconstructed on the basis of the direct age determination by¹² and field investigations in the area surrounding the archaeological site. This allowed us to revise the stratigraphic schemes reported by⁶ and¹⁰, and to frame the stratigraphic setting at FR within the geochronological picture of the Alban Hills volcanic succession³². In addition, unpublished detailed stratigraphic logs by A.G. Segre of the two boreholes carried out by¹⁴ are also used to investigate the subsoil geology.

In order to re-interpret the stratigraphic setting above the palaeontological layer of CSG, we use here data from a large exposed section that we have surveyed at La Vignaccia quarry site, located ca. 0.7 km to the SW (Fig. 1). Here, a variably preserved succession of distal Alban Hills volcanic district pyroclastic-flow deposits unconformably overlies the upper portion of the ca. 10 m thick succession of calcareous silt deposits and travertine layers occurring above the CSG palaeontological horizon. The volcaniclastic succession comprises the primary to sub-primary deposits of the main eruption cycles of the Tuscolano-Artemisio phase, including the Tufo Pisolitico di Trigoria (561 ± 2 ka), Pozzolane Rosse (456 ± 4 ka), Pozzolane Nere (407 ± 3 ka), and Villa Senni (365 ± 4 ka) caldera-forming eruptions³².

**Sampling.** The clay succession occurring between 5 and 35 m bgs recovered from the CM1 core was sampled for palaeomagnetism using standard 8 cm³ plastic cubes, with the arrow on the sampling cube's bottom pointing up-core, at the IIPU (Anagni, Italy) where the core is stored. We collected 81 samples from the center of the split-core sections, as this is less affected by coring disturbance, with an average spacing of 34 cm. The samples were oriented only with respect to the vertical; the geocentric axial dipole (GAD) field at the latitude of the coring site has an inclination of ±61°, which makes it feasible to reconstruct palaeomagnetic polarity using only Characteristic Remanent Magnetization (ChRM) inclinations. To minimize sample dehydration, samples were packed in sealed bags and stored in a refrigerated room until they were processed at the INGV in Rome.

In order to provide geochronologic constraints to the sampled interval, the core was initially scrutinized in order to detect macroscopic tephra layers. One sub-primary fallout deposit was observed between 29.70 and 29.80 m bgs and sampled for $^{40}$Ar/$^{39}$Ar dating. Four additional bulk sedimentary samples were collected at 2.60, 3.90, 5.10 and 11.50 m bgs for detrital sanidine dating (see methods section). Finally, in the fine-grained section between 7 and 35 m bgs, ten samples for biostratigraphic and sedimentological analysis were extracted.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Palaeogeographic evolution of the Tyrrhenian Sea margin of central Italy according to⁵⁴ and data from this work. GPTS = Geomagnetic Polarity Time Scale.</figcaption></figure>

## Laboratory procedures and data analysis

**Palaeomagnetism.** Natural and artificial magnetizations were measured using a narrow-access pass-through cryogenic.

magnetometer (2-G Enterprises model 755R) with internal diameter of 4.2 cm, equipped with three DC SQUID sensors (noise level $3 \times 10^{-9}$ A m² kg⁻¹), housed in a Lodestar Magnetics shielded room (ambient field <200 nT).

After measurement of the Natural Remanent Magnetization (NRM), samples were AF-demagnetized at successive peak fields of 5, 10, 15, 20, 25, 30, 40, 50, 60, 80, and 100 mT. The stability of the NRM was assessed using vector component diagrams³⁴. ChRM directions were determined using principal component analysis with linear best fits calculated from 3 or more demagnetization steps using the PuffinPlot palaeomagnetic analysis software³⁵.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Simplified stratigraphic log of the Colle Marino borehole (CM1). For detailed lithological description and sampling, see Figure S2 in Supplementary Material #1.</figcaption></figure>

Following AF demagnetization of the NRM, mineral magnetic analyses were conducted on the same set of discrete samples. Low-field magnetic susceptibility (k) was measured using an AGICO KLY-5 Kappabridge magnetic susceptibility meter with a sensitivity of 2 × 10⁻⁸ SI.

An anhysteretic remanent magnetization (ARM), which is particularly effective in activating finer magnetic grains, was imparted by using a 0.05 mT direct current (DC) bias field superimposed on a 100 mT peak AF and by translating samples through the AF and DC coil system at 10 cm/s, which is the lowest speed allowed by the control software.

Temperature dependence of the magnetic susceptibility was measured on a few milligrams of powdered samples with a furnace-equipped KLY-5 Kappabridge, in order to determine the characteristic Curie or Néel temperatures of magnetic minerals and discriminate ferromagnetic mineralogy³⁶. Samples were measured in air up to a maximum temperature of 700 °C and the thermomagnetic curves were analyzed using the Cureval8 program.

**⁴⁰Ar/³⁹Ar dating methods.** Different crystal phases were extracted from five samples of sediment (CSG-13 2A, CM1-2.6, CM1-3.9, CM1-5.1, CM1-11.5), aimed to provide terminus *post-quem* (maximum age) to the time of deposition, and from one sample of a primary fallout deposit (CM1-29.75).

Samples CM1-2.6, CM1-3.9, CM1-5.1, CM1-11.5 were prepared (crushing, sieving, cleaning, and handpicking) at the Laboratoire des Sciences du Climat et de l'Environnement of Gif-sur-Yvette (France) facility using the same procedure as in¹¹. Unfortunately, samples CM1-2.6, CM1-5.1, CM1-11.5 did not yield K-bearing mineral phases suitable for ⁴⁰Ar/³⁹Ar dating. Pristine sanidines and few leucites ranging in size between 1 mm and 300 µm were otherwise found in the sample CM1-3.9.

The selected minerals were irradiated for 120 min in the Cd-lined, in core CLICIT facility of the Oregon State University TRIGA reactor (IRR CO007). Interference corrections were based on the nucleogenic production ratios given in³⁷. After irradiation, individual crystals were transferred into a copper 134 pits sample holder placed into a differential vacuum Teledyne Cetac window connected to the extraction line. Minerals were fused one by one using a 100 Watts Teledyne Cetac CO₂ laser during 15 s at 2.3 W. Before fusion, each crystal underwent a 10 s long sweeping at 0.2 W to remove unwanted gas trapped superficially. Extracted gas was purified by a SAES GP50 cold getter for 90 s and then for 210 s by two hot SAES GP50 getters. The five Argon isotopes (i.e., ⁴⁰Ar, ³⁹Ar, ³⁸Ar, ³⁷Ar and ³⁶Ar) were measured using a NGX 600 mass spectrometer equipped with 9 ATONA Faraday cups and an electron multiplier. Technical specifications and performances of the NGX 600 ATONA detector array are presented in detail by³⁸.
$^{40}$Ar, $^{39}$Ar, $^{38}$Ar, $^{37}$Ar and $^{36}$Ar isotopes were collected simultaneously while the $^{37}$Ar was measured in a second time. In the first run, $^{40}$Ar, $^{39}$Ar and $^{38}$Ar were measured simultaneously on 3 ATONA and $^{36}$Ar using the electron multiplier. Then, just following this first run $^{37}$Ar was measured using the electron multiplier after peak switching. Each isotope measurement corresponds to 15 cycles of 20-s integration time. Peak intensity data were reduced using ArArCALC V2.4$^{39}$.

Neutron fluence J for each sample was calculated using co-irradiated Alder Creek sanidine standard ACs-2 at 1.1848 Ma$^{40}$ according to the K total decay constant of$^{41}$ ($\lambda_{e,c.} = (0.580 \pm 0.007) \times 10^{-10}$ yr$^{-1}$ and $\lambda_{\beta.} = (4.884 \pm 0.049) \times 10^{-10}$ yr$^{-1}$). To determine the neutron flux of IRR CO007 flux monitors were placed in small pits framing the sample. Two standards from 3 pits around the unknown were measured to calculate CM1-3.9 J-value ($J = 0.0005519 \pm 0.00000055$).

Mass discrimination was monitored by analysis of 30 air pipettes of various beam sizes to verify the detectors linearity. These measurements are done automatically during the nights before and after the unknown measurements. Discrimination is calculated according to the $^{40}$Ar/$^{36}$Ar ratio of 298.56$^{42}$. Procedural blank measurements were achieved after every four unknowns. For typical 5 min time blank backgrounds are between $1.7 \times 10^{-4}$ V and $2.0 \times 10^{-4}$ V for $^{40}$Ar and 65 cps for 36 Ar ($9.5 \times 10^{-7}$ V equivalent).

Separates from sample CSG-13 2A and CM1-29.75 were co-irradiated with the 1.1864 Ma Alder Creek sanidine standard$^{43}$ at the Oregon State University TRIGA reactor in the Cadmium-Lined In-Core Irradiation Tube. Single crystal fusion analyses were performed at the WiscAr laboratory at the University of Wisconsin-Madison using a 60 W CO$_2$ laser and a Noblesse multi-collector mass spectrometer following$^{43}$.

Results for all the dated samples are reported with 2σ analytical uncertainties according to the respective laboratory standards (see above). Difference due to the two different calibrations affects the third decimal figure and are negligible within the aims of this work. Full analytical data are reported in Supplementary Material #2.

**Sedimentological and biostratigraphic analysis.** Samples were oven-dried at 50 °C, weighed, and subsequently washed and sieved through a 63 µm mesh and dried again at 50 °C. Separates were observed under an optical microscope LEICA M165 and selected fossils were analyzed with the Scanning Electron Microscope HITACHI TM3030 Plus. The semiquantitative analyses carried out on the microfossil content, including ostracods, and small bivalves, led to the rough indication of the occurrence of the most indicative taxa, in order to identify palaeoenvironmental changes.

**Petrographical and SEM analysis.** Mineral compositions were analysed at the CNR-Istituto di Geologia Ambientale e Geoingegneria (Rome), with a Cameca SX50 electron microprobe equipped with five wavelength dispersive spectrometers (WDS). Quantitative compositional analyses were performed using 15 kV accelerating voltage and 15 nA beam current. As standards we employed metals for Mn and Cr, jadeite for Na, wollastonite for Si and Ca, orthoclase for K, corundum for Al, magnetite for Fe, rutile for Ti and periclase for Mg. Counting times were 20 s for elements and 10 s for backgrounds. Light elements were counted first to prevent loss by volatilization. The PAP correction method was used.

## Results

**Magnetic properties.** Within CM1 core, NRM, κ and ARM do not show significant changes with a single peak in correspondence of a thin tephra layer at 29.75 m bgs (Fig. 4). NRM intensity values range between $3.89 \times 10^{-5}$ A/m and $4.36 \times 10^{-3}$ A/m with an average of $1.22 \times 10^{-3}$ A/m. The single sample at the tephra level has an intensity of $5.12 \times 10^{-2}$ A/m. Down-core variations of κ are associated with similar changes in ARM intensity which suggests that these fluctuations are mainly controlled by changes in magnetic mineral concentration. κ values range between $8.22 \times 10^{-5}$ SI and $3.69 \times 10^{-4}$ SI with an average of $1.78 \times 10^{-4}$ SI (data are reported in Supplementary Material #3).

The increase in κ and ARM in the interval between 18 and 22 m corresponds to a portion of the core characterized by decimetric alternations of light sandy clays and yellowish silty sands.

At the tephra level $κ = 1.17 \times 10^{-2}$ SI. ARM values range between $1.31 \times 10^{-3}$ A/m and $4.32 \times 10^{-2}$ A/m with an average of $7.44 \times 10^{-3}$ A/m. The single sample at the tephra level has an intensity of $1.70 \times 10^{-1}$ A/m.

Thermomagnetic curves on three selected samples (Fig. 4e) were produced via progressive heating. Samples CM29 (19.35 m bgs) and CM62 (33.50 m bgs) show a noisy susceptibility versus temperatures curves due to a weak κ but it has been still possible to observe a major drop occurring at ca. 580 °C (temperatures corresponding to the inflection point during the decrease in susceptibility) which indicates the ubiquitous presence of Fe spinels with a composition similar to magnetite. A more complex mineralogy occurs for tephra sample CM55 (29.5 m bgs) where a decrease at ca. 415 °C is followed by a second drop at ca. 530 °C which indicate the presence of unstable cation deficient titanomagnetite exsolved into Ti-rich phases upon heating in air$^{44}$.

**Demagnetization behavior and polarity zonation.** Stepwise AF demagnetization enabled isolation of the ChRM component for 63 (78%) of the samples analyzed. Most of the samples are characterized by noisy demagnetization paths in orthogonal vector diagrams but are sufficiently good to define a magnetic polarity zonation. The unstable behavior is reflected in a with maximum angular deviation (MAD) value comprised between 1.0° and 20.6°, with an average of 7.9°. Some samples, mainly from the sandy horizons, are characterized by aberrant behavior (Fig. 5) due to the fact that the deposition of large particles would be controlled by gravitational rather than magnetic forces. Thus, their orientation could not be expected to represent the geomagnetic field at or near the time of deposition and have been discarded.

In order to avoid the shallowing bias which can be introduced when taking the arithmetic mean of the inclination values, mean ChRM inclinations were computed separately for normal and reverse polarity populations
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Mineral magnetic data for CM1. (a) lithology, (b) stratigraphic variations of NRM intensity, (c) κ and (d) ARM. (e) Temperature dependence of κ (up to 700 °C) for three representative samples: CM29 (19.35 m), CM55 (29.5 m), CM62 (33.5 m). For samples CM29 and CM62 a three-point moving average is applied (red curves).</figcaption></figure>

using a maximum likelihood method⁴⁵ as implemented in the PuffinPlot application³⁵. The normal polarity samples and the reverse polarity samples have a mean inclination of 57.8° (N=13; a95 = 21.8°) and -55.6° (N=49; a95 = 5.6°), respectively. These values are not significantly different from antipodal at the 95% confidence level, affording a positive reversal test, and suggesting that the mean directions could indeed represent a reliable polarity record.

The ChRM inclinations enable delineation of eight magnetozones, which are defined using at least two consecutive samples with inclinations distinctly different from neighboring intervals (Fig. 6). Boundaries between magnetozones are placed at the mid-point between successive samples of opposite polarity. Starting from the bottom, the lower 5.73 m (from 34.55 to 28.82 m bgs) of the magnetic polarity record has reverse polarity (magnetozone R1); the lower boundary of this magnetozone was not sampled. Magnetozone N1, from 28.82 to 28.00 m bgs, has two samples with stable normal polarity. The boundary between magnetozones N1 and R2 lies within an 0.65 m interval (indicated in grey in Fig. 6). Above magnetozone N1, a reverse polarity interval (R2) extends from 28.00 to 20.12 m bgs and it is followed, from 20.12 to 19.15 m bgs, by magnetozone N2 defined by two samples with stable normal polarity. The boundary between magnetozones R2 and N2 lies within an 0.67 m interval with no samples. The reversed polarity magnetozone R3 is well defined and extends from 19.15 and 10.22 m bgs. The uppermost 10.22 m of the studied sequence is dominated by two normal polarity intervals (N3
<figure><img src="image_8.png" /><figcaption>**Figure 5.** AF (alternating field) demagnetization behavior for 9 representative samples from CM1. For the vector component diagrams, open (closed) symbols represent projections onto the vertical (horizontal) plane. The stereoplots are equal-area projections, with solid (open) symbols representing points projected onto the lower (upper) hemisphere. The cores were not azimuthally oriented; declinations are reported in the laboratory coordinate system with respect to the split face of the core. Plots were produced using PuffinPlot³⁵.</figcaption></figure>

and N4) separated by a reverse polarity interval (R4), from 9.47 to 8.82 m bgs. Samples at 9.52 (normal polarity) and 9.41 m bgs (reversed polarity) define the boundary of magnetozone N3 with R4. The boundary position between R4 and N4 has a range of uncertainty of about 1 m due to the lack of samples in this interval with acceptable demagnetization behaviors useful for defining the ChRMs. Data are reported in Supplementary Material #3.

⁴⁰Ar/³⁹Ar dating. CM1-3.9. Twenty-one sanidine crystals were dated and provided very scattered ages ranging between about 955 and 689 ka (see probability diagram in Fig. 7a and Supplementary Material #2). The
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Lithology (symbols as in Fig. 3), downcore variations of ChRM inclination, MAD, and magnetic polarity zonation (black=normal polarity and white=reverse polarity).</figcaption></figure>

youngest crystal is dated to 688.5 ± 2.7 ka (2σ analytical uncertainty). A population constituted by 9 on 21 dated crystals is dated to 843.3 ± 1.6 ka (MSWD = 1.20; P = 0.30, 2σ analytical uncertainty including J value).

CM1-29.75. Nineteen of 22 K-bearing plagioclase crystals give a weighted mean $^{40}$Ar/$^{39}$Ar age of 2.214 ± 0.034 Ma (Fig. 7b). Three dates are markedly younger and are excluded from the weighted mean calculations. These younger dates may be due to alteration or Ar loss.

CSG-13 2A. Eighteen of 22 single plagioclase fusions from sample CSG-13 2A give a weighted mean $^{40}$Ar/$^{39}$Ar age of 2.233 ± 0.032 Ma (Fig. 7b). Four older dates are excluded from the weighted mean calculations and are interpreted to be antecrysts or partially degassed xenocrysts.

**Age model.** We refer the geomagnetic polarity time scale (GPTS) of¹⁶ to develop an age model for the studied interval. Our magnetostratigraphic interpretation is presented in Fig. 8 along with $^{40}$Ar/$^{39}$Ar ages from both primary and reworked volcanic materials and the uncertainties associated with each datum. An $^{40}$Ar/$^{39}$Ar age of 2.21 ± 0.034 ka (2σ analytical uncertainty) for a tephra layer at 29.75 m bgs constrains correlation of magnetozone R1 to Chron C2r.2r (2.148–2.588 Ma). The youngest sanidine crystal extracted from a sediment sample at 3.9 m bgs, is dated to 688.5 ± 2.7 ka providing a maximum age for this horizon and constrains correlation of
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Relative probability diagrams showing the results of the 40Ar/39Ar dating experiments performed at LSCE (a) and at the Wiscar Laboratory (b). Full analytical data are reported in Supplementary Material #2, #3.</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>**Figure 8.** Correlation of the new magnetic polarity zonation from CM1 with the GPTS of$^{46}$.</figcaption></figure>

magnetozone N4 to Chron C1n (Brunhes; 0–0.781 Ma). This is in agreement with the absence of primary volcanic deposits ascribable to the Alban Hills and Monti Sabatini volcanic districts, confirming that the top of the sequence, from 0 to 3.9 m bgs, is older than ~ 0.6 Ma.

The two $^{40}$Ar/$^{39}$Ar ages and the absence of products of the paroxysmal volcanic activity, support a one-by-one correlation of observed magnetozones with each polarity of the GPTS from Chron C1n to C2r.2r with variable sedimentations rate up to 4.81 cm/kyr.

**Biostratigraphy.** The investigated clayey succession is remarkably barren, with a total absence of any bio-
genic fraction in all the samples collected between 12 and 35 m bgs, as well as in one sample collected at the very
top, 7.5 m bgs (Table 1).

A scanty ostracod assemblage has been found only in a sample at 9.5 m bgs, represented by *Ilyocypris gibba*
and *Candona* sp., along with a few undetermined bivalve shells (Fig. S4a in Supplementary Material #1). The
occurrence of *Ilyocypris gibba* with *Candona* sp. indicates fresh—oligohaline inland water body, with a mud
substrate$^{4,47}$. However, idiosyncratic calcified Charophytes thalli have been also observed in this sample (Fig. S4b).

Indeed, Charophytes (Characeae), are macroscopic green algae occurring mainly in the littoral zone of lakes.
Photosynthetic activity of charophytes leads to precipitation of autochthonous carbonates, which significantly
contribute to lacustrine sedimentation. The occurrence of most charophytes is limited between about 1 m to
10 m of depth in fresh water bodies with clear, alkaline waters, while they become rare along with increasing
trophy level and decreasing light availability (Garcia, 1994 in$^{48}$).

**Tephrostratigraphic and biochronologic data.** The tephra at 29.75 m bgs (sample CM1-29.75) con-
tains abundant amphibole and plagioclase crystals. Amphibole consists of submillimeter, euhedral to subeuhe-
<table><caption>**Table 1.** Composition of the mineralogical and biogenic fraction > 63 µm in the investigated samples.</caption><thead><tr><td>Sample-depth</td><td>Terrigenous fraction</td><td>Biogenic fraction</td></tr></thead><tbody><tr><td>CM 1 – 7.2</td><td>quartz+++, feldspar +, biotite-, oxidised carbonates concrections</td><td>Barren</td></tr><tr><td>CM 1 – 9.5</td><td>quartz+++ feldspar +, pyroxene-, biotite-,</td><td>Ostracods (Candona sp. and Cypridea torosa), rare undetermined Bivalves, calcified remains of charophyte cortex (Characee)</td></tr><tr><td>CM 1 – 12.7</td><td>quartz + + +, feldspar +, pyroxene-, biotite-, white pumice + +</td><td>Barren</td></tr><tr><td>CM 1 – 15.50</td><td>quartz + + +, feldspar + +, biotite +</td><td>Barren</td></tr><tr><td>CM 1 – 18.9</td><td>quartz + + +, feldspar +, pyroxene +, biotite +, yellowish pumice + +</td><td>Barren</td></tr><tr><td>CM 1 – 21.6</td><td>quartz + + +, feldspar +, biotite + +, pyroxene +, yellowish pumice + +</td><td>Barren</td></tr><tr><td>CM 1 – 26.60</td><td>quartz + + +, feldspar + +, biotite +</td><td>rare vegetal fragments</td></tr><tr><td>CM 1 – 29.3</td><td>quartz + + +, feldspar +, amphibole + +, pyroxene +, biotite +, pyrite +</td><td>Barren</td></tr><tr><td>CM 1 – 33.00</td><td>quartz + + +, feldspar +, biotite + +, pyroxene + +</td><td>Barren</td></tr><tr><td>CM 1 – 35</td><td>quartz + + +, feldspar +, pyroxene-, biotite-, pyrite-</td><td>Barren</td></tr></tbody></table>

dral, deep green crystals including titanomagnetite. According to the classification scheme of⁴⁹, the composition of amphibole is edenite with an Mg# restricted to a nearly constant value of 60–64. Plagioclase is found as submillimeter, euhedral to subeuhedral crystals, including apatite and is characterized by continuous normal zoning with anorthite (An₄₇–₅₇) content variable from core-to-rim.

The mineral chemistry of the tephra is typical of calcalkaline rocks. In particular, the Mg# of amphibole and the anorthite content of plagioclase indicate a crystallization from an andesite magma⁵⁰. The calcalkaline rocks underlying the Phlegrean Fields volcanic suites show similar K/Ar age⁵¹ of the CM1-29.75 tephra. Nevertheless, buried Campanian volcanics show porphyritic textures with phenocrysts of zoned plagioclase but without amphiboles. The latter phase occurs in the calcalkaline rocks of Capraia island but are much older (> 4.6 Ma;⁵²).

The CSG-13 2A sedimentary sample contained a significant amount of plagioclase crystals; 18 out of 22 yielded a homogeneous population age of 2.233 ± 0.032 Ma, virtually undistinguishable from that of 2.214 ± 0.034 Ma yielded by the primary tephra sample CM1-29.75, suggesting the sub-primary origin for the volcanic material included in the sediment. Indeed, photographs of the original excavation from which sample CSG-13 2A was collected (Fig. S5 in Supplementary Material #1) show the occurrence of cross-laminated, dark grey sands, which may be regarded as deriving from the reworking of a primary fallout deposit emplaced shortly before, as suggested by the abundant, homogeneous plagioclase crystals population recovered in the sedimentary sample. Consistent with this inference, the palaeontological layer of CSG occurs at similar, yet slightly higher elevation with respect to the tephra layer recovered from CM1 (Fig. 9). Therefore, the age of 2.233 ± 0.032 Ma can be considered the same of the faunal assemblage of Middle Villafranchian age recovered at CSG⁴.

The CSG faunal unit belong to the MNQ17b biozone the reference site of which is located in the French *Massif Central* (Senèze) and is dated by ⁴⁰Ar/³⁹Ar between 2.10 and 2.18 Ma⁵³. The age obtained for the fossil-bearing horizon at CSG suggests an earlier arrival of some species in Italy than in France and give a minimum age of 2.23 Ma for the base of the MNQ 17b biozone.

## Discussion

The chronostratigraphic constraints achieved for the sedimentary succession recovered from the borehole and documented in the field at the three palaeontological/archaeological sites show that these successions cover three distinct, partially overlapping, time intervals spanning ca. 2.6 through 0.36 Ma (Fig. 9). The tephra recovered at CM1 provides exact correlation with the fossil-bearing horizon of CSG, where a nearly homogeneous K-plagioclase crystal population extracted from a sand sample yielded an identical age, within the analytical errors, of 2.233 ± 0.032 Ma.

Remarkably, the clay section underlying the fossiliferous horizon at CSG1 borehole, between 198 and 193 m a.s.l., shows lithological and biostratigraphic features very similar to the clay interval underneath the tephra layer from the CM1 borehole. The clay section was described by⁴ as part of a floodplain facies association, characterized by dominant greenish/bluish, fine-grained, laminated sediments, with common CaCO₃ nodules, rootlets, and very scanty organic material. Few ostracods (*Pseudocandona rostrata* with subordinate *Candonopsis* juv.) have been recovered only in the upper, coarsening upward part of the facies association (201–198 m a.s.l.), which makes transition to the fluvial-channel facies association (206–201 m a.s.l.) in which a Middle Villafranchian vertebrate faunal assemblage occurs⁴⁻¹⁰.

In contrast, the lower portion of the clay section recovered from 193 to 163 m a.s.l. at CSG1 displays distinct features from those detected from 220 to 192 m a.s.l. (from 7 to 35 m bgs) at CM1, consistent with its different, older age range.

By combining the abundant palynological and ostracod records with facies analysis of the sedimentary succession⁴, recognized two near-shore facies associations at the top and at the bottom of the cored sediments, separated by a thick marsh facies association, with abundant organic-rich (lignite) layers.

No lower geochronologic boundary is achievable for the CSG1 clay section. A minimum age ≤ 2.58 Ma can be inferred based on the reversed polarity yielded by the whole analyzed set of samples after⁵, with the exception of two. Still, once the high sedimentation rate deriving from this assumption is taken into account, the possibility that the lowest normal polarity sample occurring at the base of the clay section may represent the top of the Gauss normal Chron C2An.1n (2.581 Ma) cannot be discarded.
<figure><img src="image_12.png" /><figcaption>**Figure 9.** (a) Reconstruction of the stratigraphic setting of the investigated area through the cross-correlation of the borehole data achieved in the present study and in previous work at Colle Marino, Fontana Ranuccio and Coste San Giacomo. (b) Sketch showing the environmental evolution of the Latin Valley in the interval 2.6–0.6 Ma, after restoration of the original elevation of borehole FR2 before tectonic dislocation (see text for explanations). Drawings by FM.</figcaption></figure>

The low elevation at which the Brunhes-Matuyama reversal is found in borehole FR2, combined with the large offset affecting the lithological boundary between the white calcareous silt and the grey lacustrine clay from cores FR1 and FR2<sup>14</sup>, suggests the occurrence of fault displacement (Fig. 9a). In addition, both the large thickness of the Pozzolane Rosse pyroclastic-flow deposit (456 ± 4 ka;<sup>32</sup>) and the occurrence of a deeply zeolitized basal por- tion point to deposition within an incised palaeovalley corresponding to a tectonic structure (“graben”) which would host the present Fosso delle Mole streambed (Fig. 1). Similarly, a thick succession of reworked volcanic material, older than the Pozzolane Rosse emplaced in fluvial-lacustrine setting at FR2, constrains the tectonic activity within the 600–450 ka time interval.

After the abovementioned tectonic offset is restored, the stratigraphic setting around 600 ka is reconstructed in Fig. 9a, while a sketch of the palaeogeographic conditions controlling the geologic evolution of this area since 2.6 Ma is shown in Fig. 9b.

Biostratigraphic features of the lower clay section recovered at CSG1 reflect the regional relative sea-level oscillations of the Tyrrhenian Sea basins during the Piacenzian (Fig. 10;⁵⁴ and references therein). The nearshore deposits found at the base and at the top of this succession testify approaching of the coastline during the Pia-enzian ingression 3.6 through 2.6 Ma, and during the early stages of the Gelasian/Santernian ingression, since 2.2 Ma. This evolution reflects the palaeogeographic conditions of the Latin Valley, which was an eastward stretching, deep gulf of the Tyrrhenian Sea during the two ingressive phases, while it became an estuarine, large alluvial plain during the Gelasian regressive phase, when the thick marshy facies association was deposited.

However, as suggested in⁵⁴, the Santernian ingression in the Latin Valley was aborted shortly after 2.2 Ma, in consequence of a tectonic "cutoff" of this engulfed region. Very likely, this process was due to the re-activation of a N-S strike-slip fault (Fig. 10c; Cisterna Fault Zone—CFZ⁵⁵), which created a structural threshold hindering the sea ingression, and forming a large lacustrine basin occupying the whole Sacco Valley (Figs. 1, 2, 10).

Previous literature⁵⁶,⁵⁷ interpreted the eastward "inversion" of the drainage in the Sacco Valley, from an original westward direction towards the Aniene and Tiber rivers catchment basins, has due to the emplacement of
<figure><img src="image_13.png" /><figcaption>**Figure 10.** Paleogeographic sketches showing the sea ingression and regression phases occurring in the investigated area 2.5 Ma through 1.8 Ma. See text for explanations. DEM images used with permission of the Istituto Nazionale di Geofisica e Vulcanologia, Rome.</figcaption></figure>

large thickness of Alban Hills pyroclastic-flow deposits since 450 ka. Palaeogeographic reconstruction based on the CSG1 and CM1 borehole data discussed here suggest that a significant change in the hydrographic catchment occurred much earlier, around 2.2 Ma, creating the conditions for the successive inversion of the river flow direction.

Indeed, rather than evolving towards a marine environment, like what happened along the northern stretch of the "Santernian coastline" (⁵⁴,⁵⁸ for a detailed discussion) (Fig. 10), a continental, lacustrine environment characterized the Sacco Valley from 2.2 through 0.6 Ma, at the least.

The lacustrine deposits spanning this age interval represent the entire, 35 m thick sedimentary succession recovered between 192 and 227 m a.s.l. at the CM1 borehole, and they constitute the uppermost 12 m, from 195 to 206 m a.s.l., at CSG1. However, the geologic section exposed at La Vignaccia allows a correlation of the thick succession of calcareous silt with frequent carbonate encrustations cropping out between 210 and 225 m a.s.l. with the coeval lacustrine succession in CM1. The top of the calcareous silt deposit reaches an elevation of 225 m a.s.l. in La Vignaccia, thus leading to estimate a thickness of ca. 20 m for this pack of homogeneous lacustrine sediments overlying the fossiliferous horizon of CSG. The "volcanic" age boundary at ~0.6 Ma shows that this sedimentary sequence accumulated during a time span of ca. 1.5 Myr, since 2.21 Ma. The thick suite of primary to deeply reworked volcanic materials exposed at La Vignaccia is correlated to the deposit described in the stratigraphic log of FR's boreholes 1 and 2, occurring between 232–227 and 217–202 m a.s.l., respectively. This sequence of reworked volcanics rests on the eroded sedimentary substrate, while its base deepens to the east, towards the depocenter of the lacustrine basin, that is, where FR and CM are located. In general, this correlation allows us to understand the palaeoenvironmental features of the lacustrine basin. With a deep depocenter near CM1, barren clay sediments accumulated in this area, and a peripheral, shallow coastal setting, characterized by clear waters allowing photosynthetic activity of charophytes, led to precipitation of autochthonous carbonates that substantially contributed to the lacustrine sedimentation.

These environmental conditions are framed within a significant climate event which possibly is the driving factor of the Gelasian/Santernian ingression and to the dispersal events that led into Europe African taxa such as, among others, *Homo* sp.⁵⁹. Moreover, the CSG faunal assemblage⁴ testifies the early-middle Villafranchian mammal turnover, characterized by the dispersal in Europe of species that adapted to open environments such as elephants and horses⁶⁰.

Since ~0.8 Ma, progressive shallowing of the lake resulted in an eastward "progradation" of the "Charophyte facies". Consistent with appearance of Charophytes and the abundant carbonate precipitation testified by the dif-fuse nodules occurring in the uppermost 3 m of the clay succession, a sharp transition to white, calcareous sandy silts occurs at 220 m a.s.l. at CM1 (Fig. 9). The transition is marked by a decimeter-thick interval of alternating thin layers of white carbonatic silt and dark-brown, organic-rich clay sediment. At 200 m a.s.l. the calcareous silt horizon passes upwards (the lithologic passage was not recovered in the core due to the partial loss of sampling caused by the unconsolidated features of the silt deposits) to yellow sands with abundant carbonate encrustations and travertine layers, in which two red soils occurs. The youngest sanidine crystal out of 21 extracted from a sediment sample collected within this horizon (CM1-3.9) provided a maximum age of 688.8 ± 2.7 ka for this layer.

We interpret the progressive shallowing and temporary emersion of the lacustrine basin as the combined effect of the sedimentary filling and regional uplift that occurred since 0.8 Ma and culminated around 0.6 Ma, in coincidence with the onset of the volcanic activity in the Roman Magmatic Region. In this context, the lacustrine basin of the Sacco Valley can be regarded as a fault-controlled, subsiding basin, within the overall uplifting regime of the Apennine chain during Middle Pleistocene times (e.g.,⁶¹⁻⁶³).

This location, which contains the earliest evidence of human occupation in central Italy (e.g.,⁶⁴,⁶⁵) is characterized by large lacustrine basins and surrounding alluvial plains with prairies and grassland and mild climatic conditions, and represent an ideal environment for game and their hunters.
## Conclusions

The present study provides accurate geochronologic constraints to the sedimentary succession recovered from a drilling at CM and unambiguously correlate it with those previously investigated at the nearby fossiliferous sites of CSG and the archeological site of FR. The main findings of this study are:

* Correlation of the sequence of magnetozones with the GPTS is constrained by $^{40}$Ar/$^{39}$Ar ages obtained from both in situ and reworked volcanic materials.

* The chronostratigraphic constraints show that the investigated successions cover three distinct time intervals, partially overlapping, spanning ca. 2.6–0.36 Ma.

* The lower clay succession recovered at CSG records two regional sea-level oscillations corresponding to the "Piacenzian" (i.e., 3.6–2.6 Ma) and to the early "Gelasian/Santernian" (i.e., 2.2–1.5 Ma) ingressions.

* The second ingressive phase is abruptly aborted around 2.2 Ma, as attested by the lack of marine deposits of Santernian age (1.8–1.5 Ma) in the CM1 borehole.

* A large lacustrine basin occupied the Latin Valley since 2.2 Ma providing an environment favorable to the dispersal of several "African" taxa which, among others, constitute the Middle Villafranchian Faunal Unit of CSG which is dated in the present work at 2.233 ± 0.032 Ma.

* The age of CSG fossiliferous horizon slightly pre-dates previous assessment at the French type-locality of the MNQ 17b biozone, suggesting an earlier occurrence of the African taxa in Italy than in France.

* Starting from 0.8 Ma, the progressive shallowing and temporary emersion of the large lacustrine basins created favorable conditions for early hominin occupation of the area, as it is possibly documented by the lithic assemblages from CM.

* The CM layer dated at about 0.7 Ma, tentatively suggested as being the origin of the lithic collection found in the site, could provide a strong reference for the pre-Acheulean technocomplexes in Mediterranean Europe.

* The human occupation of this region reaches a climax around 0.4 Ma as constrained by the chronostratigraphic framework reconstructed here for the FR site, consistent with geochronologic data from recent literature achieved at this site and several other sites throughout the Latin Valley, spanning 0.41–0.35 Ma.

These results contribute to better interpretation of the hominin migratory dynamics and the factors that influenced the location and spatial distribution during the early stages of occupation in these regions.

Received: 11 January 2021; Accepted: 24 February 2021

Published online: 29 March 2021

## References

1. Ascenzi, A., Biddittu, I., Cassoli, P. F., Segre, A. G. & Segre Naldini, E. A. calvarium of late Homo erectus from Ceprano, Italy. *J. Hum. Evol.* **31**, 409–423 (1996).

2. Ascenzi, A., Mallegni, F., Manzi, G., Segre, A. G. & Segre Naldini, E. A. a re-appraisal of Ceprano calvaria affinities with Homo erectus, after the new reconstruction. *J. Hum. Evol.* **39**, 443–450 (2000).

3. Manzi, G., Mallegni, F. & Ascenzi, A. A cranium for the earliest Europeans: phylogenetic position of the hominid from Ceprano, Italy. *Proc. Natl. Acad. Sci.* **98**, 10011–10016 (2001).

4. Bellucci, L. et al. The site of Coste San Giacomo (Early Pleistocene, central Italy): Palaeoenvironmental analysis and biochronological overview. *Quatern. Int.* **267**, 30–39 (2012).

5. Bellucci, L. et al. Evidence of late Gelasian dispersal of African fauna at Coste San Giacomo (Anagni Basin, central Italy): early Pleistocene environments and the background of early human occupation in Europe. *Quatern. Sci. Rev.* **96**, 72–85 (2014).

6. Biddittu, I. et al. Anagni, a K-Ar dated lower and middle pleistocene site. *Quat. Nova* **21**, 53–71 (1979).

7. Cassoli, P. F. & Segre Naldini, E. *Nuovo contributo alla conoscenza delle faune villafranchiane e del Pleistocene medio del bacino di Anagni (Frosinone)* 115–118 (Atti XXIV Riun. Sc. IIPP, 1984).

8. Segre, A. G. Guida ai giacimenti Pleistocenici e Villafranchiani del Lazio meridionale. *Quat. Nova* **7**, 205–232 (2004).

9. Segre Naldini, E. & Valli, A. M. F. Villafranchian Cervids from Central Italy. *Quaternaria Nova* **7**, 159–204 (2004).

10. Segre Naldini E, Muttoni G, Parenti F, Scardia G, Segre AG (2009) Nouvelles recherches dans le bassin Plio-Pléistocène d'Anagni (Latium méridional, Italie). *L'Anthropologie* **113**, 66–77.

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

12. 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).

13. Biddittu I., Segre A. G. & Segre Naldini, E. Fontana Ranuccio in *I primi abitanti d'Europa. 1.500.000 - 100.000 anni* (ed. Biddittu I., Segre A. G., Segre Naldini, E.) 132–134 (De Luca, 1985).

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

15. Ascenzi, A., Biddittu, I., Cassoli, P. F., Segre, A. G. & Segre Naldini, E. Dal Paleolitico inferiore al Bronzo nel bacino di Anagni in *Dives Anagnia dalla preistoria alla protostoria. Archeologia della Valle dei Sacco* (ed. Gatti, S.) 25–47 (L’Erma di Bretschneider: Roma, 1993).

16. Ascenzi, A. & Segre, A. G. 1996) Artefacts and human teeth at the Fontana Ranuccio Middle Pleistocene Site (Central Italy. *Anthropologie* **34**, 39–46 (1996).

17. Rubini, M., Cerroni, V., Festa, G., Sardella, R. & Zaio, P. A revision of hominin fossil teeth from Fontana Ranuccio (Middle Pleistocene, Anagni, Frosinone, Italy). *J. Hum. Evol.* **77**, 204–216 (2014).

18. Grimaldi, S., Santaniello, F., Angelucci, D. E., Bruni, L. & Parenti, F. A Techno-Functional interpretation of the lithic assemblage from Fontana Ranuccio (Anagni, Central Italy): an insight into a MIS 11 human behaviour. *J. Paleolithic Archaeol.* https://doi.org/10.1007/s41982-020-00064-3 (2020).
19. Gliozzi, E. et al. Biochronology of selected Mammals, Molluscs and Ostracods from the Middle Pliocene to the Late Pleistocene in Italy. The state of the art. *Rivista Italian. Paleontol. Stratig.* **103**, 369–388 (1997).

20. Petronio, C., Bellucci, L., Martinetto, E., Pandolfi, L. & Salari, L. Biochronology and Palaeoenvironmental Changes from the Middle Pliocene to the Late Pleistocene in Central Italy. *Geodiversitas* **33**, 485–517. https://doi.org/10.5252/g2011n3a4 (2011).

21. Biddittu, I. & Segre, A. G. Giacimenti preistorici e quaternario della Provincia di Frosinone. *Boll. Ist. St. e Arte del Lazio Merid.* **9**, 21–44 (1977).

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

23. 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 (ed. Ghini, G. & Mari, Z.) 437–443 (Edizioni Quasar, 2012).

24. Cauche, D., Celiberti, V., Barsky, D., Notter, O., Biddittu, I., Lumley, H. de. Les plus anciennes industries lithiques du Latium, Italie. *Actes du XIVème Congrès de l'U.I.S.P.P, Université de Liège, Belgique, 2–8 septembre 2001, BAR International Series* **1272**, 49–57 (2004).

25. Villa, P. Early Italy and the colonization of Western Europe. *Quat. Int.* **75**, 113–130 (2001).

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

27. 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).

28. Kinematic relations since late Tortonian times. Patacca, E., Sartori, R. & Scandone, P., 1990. Tyrrhenian Basin and Apenninic Arcs. *Mem. Soc. Geol. It.* **45**, 425–451 (1990).

29. Carrara, C., Frezzotti, M. & Giraudi, C. Stratigrafia plio-quaternaria, in *Lazio Meridionale, Sintesi delle ricerche geologiche multi-disciplinari* (ed. Carrara, C.) 62–85 (ENEA, Dip. Ambiente, Serie Studi e ricerche, 1995).

30. Segre, A. G. & Asceni, A. Fontana Ranuccio: Italy's Earliest Middle Pleistocene Hominid Site. *Curr. Anthropol.* **25**, 230–233 (1984).

31. Marra, F. et al. The Volsci Volcanic Field (central Italy): an open window on continental subduction processes. *International Journal of Earth Sciences*, in press (2020).

32. Marra, F., Karner, D. B., Freda, C., Gaeta, M. & Renne, P. R. Large mafic eruptions at the Alban Hills Volcanic District (Central Italy): chronostratigraphy, petrography and eruptive behavior. *J. Volc. Geoth. Res.* **179**, 217–232. https://doi.org/10.1016/j.jvolg.eores.2008.11.009 (2009).

33. Freda, C. et al. CO₂-driven large mafic explosive eruptions: the Pizzolane Rosse case study from the Colli Albani Volcanic District (Italy). *Bull. Volc.* **2010**, https://doi.org/10.1007/s00445-010-0406-3 (2010).

34. Zijderveld, J. D. A. A. C. demagnetization of rocks: Analysis of results, in *Methods in Palaeomagnetism*, Dev. In *Solid Earth Geophys* (eds Collinson, D. W. et al.) 254–286 (Elsevier, 1967).

35. Lurcock, P. C. & Florindo, F. New developments in the Puffin Plot paleomagnetic data analysis program. *Geochem. Geophys. Geosyst.* **20**, 5578–5587. https://doi.org/10.1029/2019GC008537 (2019).

36. Hrouda, F. A technique for the measurement of thermal changes of magnetic susceptibility of weakly magnetic rocks by the CS-2 apparatus and KLY-2 Kappabridge. *Geophys. J. Int.* **118**, 604–612 (1994).

37. Balbas, A., Koppers, A. A. P., Kent, D. V., Konrad, K. & Clark, P. U. Identification of the short-lived Santa Rosa geomagnetic excursion in lavas on Floreana Island (Galapagos) by ⁴⁰Ar/³⁹Ar geochronology. *Geology* **44**, 359–362. https://doi.org/10/f8kdpg (2016).

38. Cox, S. E., Hemming, S. R. & Tootell, D. The Isotopx NGX and ATONA Faraday amplifiers. *Geochronology* **2**, 231–243. https://doi.org/10.5194/gchron-2-231-2020 (2020).

39. Koppers, A. A. P. ArArCALC – software for ⁴⁰Ar/³⁹Ar age calculations. *Comput. Geosci.* **28**, 605–619 (2002).

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

41. 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).

42. Lee, J.-Y. et al. A redetermination of the isotopic abundances of atmospheric Ar. *Geochim. Cosmochim. Acta* **70**, 4507–4512 (2006).

43. 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).

44. Readman, P. W. & O'Reilly, W. Magnetic properties of oxidized (cation-deficient) titanomagnetites (Fe, Ti, Z)₃O₄. *J. Geomag. Geoelec.* [object Object]-[object Object]. (1972).

45. Arason, P. & Levi, S. Maximum likelihood solution for inclination-only data in paleomagnetism. *Geophys. J. Int.* **182**, 753–771. https://doi.org/10.1111/j.1365-246X.2010.04671.x (2010).

46. Gradstein, F. M., Ogg, J.G., Schmitz, M.D. & Ogg, G.M. The Geologic Time Scale 2012, pp. 1144 (Elsevier, Oxford 2012).

47. Frenzel, P., Keyser, D. & Viehberg, F. A. An illustrated key and (palaeo)ecological primer for Postglacial to recent Ostracoda (Crustacea) of the Baltic Sea. *Boreas* **39**, 567–575 (2010).

48. Apolinarskaa, K., Pelechatyb, M. & Pukacz, A. CaCO₃ sedimentation by modern charophytes (Characeae): can calcified remains and carbonate δ¹³C and δ¹⁸O record the ecological state of lakes? A review. *Stud. Linnol. Telmatol.* **5**, 55–66 (2011).

49. Leake, B. E. et al. Nomenclature of amphiboles: report of the subcommittee on amphiboles of the International Mineralogical Association, commission on new minerals and mineral names. *Can. Mineral.* **35**, 219–246 (1997).

50. Tecchiato, V. et al. Petrological constraints on the high-Mg basalts from Capo Marargiu (Sardinia, Italy): evidence of cryptic amphibole fractionation in polybaric environments. *J. Volcanol. Geotherm. Res.* **349**, 31–46 (2018).

51. Barbieri, M., Di Girolamo, P., Locardi, E., Lombardi, G. & Stanzione, D. Petrology of the calc-alkaline volcanics of the Parete-2 well (Campania, Italy). *Per. Mineral.* **48**, 53–74 (1979).

52. Peccerillo, A. Cenozoic Volcanism in the Tyrrhenian Sea Region in Advances in Volcanology, pp. 399; https://doi.org/10.1007/978-3-319-42491-0 (Springer International Publishing, 2017).

53. Nomade, S. et al. ⁴⁰Ar/³⁹Ar constraints on some French landmark Late Pliocene to Early Pleistocene large mammalian paleofaunas: Paleoenvironmental and paleoecological implications. *Quat. Geochronol.* **21**, 2–15 (2014).

54. Marra, F. et al. A review of the geologic sections and the faunal assemblages of Aurelian Mammal Age of Latium (Italy) in the light of a new chronostratigraphic framework. *Quat. Sci. Rev.* **181**, 173–199. https://doi.org/10.1016/j.quascirev.2017.12.007 (2018).

55. Frepoli, A. et al. Seismicity, seismogenic structures and crustal stress field in the greater area of Rome (Central Italy). *J. Geophys. Res.* **115**, B12303 (2010).

56. De Rita, D., Funiciello, R. & Rosa, C. L'influenza dell'attività del Complesso vulcanico dei Colli Albani sullo sviluppo idrografico della Campagna Romana: l'alta valle del Fiume Sacco. *Rend. Soc. Geol. It.* **12**, 50 (1989).

57. De Rita, D., Funiciello, R. & Rosa, C. Volcanic activity and drainage network evolution of the Colli Albani area (Rome, Italy). *Acta Volc.* **2**, 185–198 (1992).

58. Mancini, M., D'Anastasio, E., Barbieri, M. & De Martini, P. M. Geomorphological, paleontological and ⁸⁷Sr/⁸⁶Sr isotope analyses of early Pleistocene paleoshorelines to define the uplift of Central Apennines (Italy). *Quatern. Res.* **67**, 487–501 (2007).

59. O'Regan, H. J., Turner, A., Bishop, L. C., Elton, S. & Lamb, A. L. Hominins without fellow travellers? First appearances and inferred dispersals of Afro-Eurasian large mammals in the Plio-Pleistocene. *Quatern. Sci. Rev.* **30**, 1343–1352. https://doi.org/10.1016/j.quascirev.2009.11.028 (2011).

60. Rook, L. & Martinez-Navarro, B. Villafranchian: the long story of a Plio-Pleistocene European large mammal biochronologic unit. *Quatern. Int.* **219**, 134–144 (2010).
61. Cavinato, G. P., Cosentino, D., De Rita, D., Funiciello, R. & Parotto, M. Tectonic-sedimentary evolution of intrapenninic basins and correlation with the volcano-tectonic activity in Central Italy. *Memorie Descrittive della Carta Geologica d'Italia* **49**, 63–76 (1994).

62. Hyppolite, J. C., Angelier, J. & Roure, F. A major geodynamic change revealed by Quaternary stress patterns in the southern Apennines (Italy). *Tectonophysics* **230**, 199–210 (1994).

63. 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). *Geol Soc. Am. Bull.* **129**, 1314–1336 (2017).

64. Orain, R. et al. Hominin responses to environmental changes during the Middle Pleistocene in central and southern Italy. *Clim. Past* **9**, 687–697 (2013).

65. Nicoud, E. et al. Preliminary data from valle Giumentina pleistocene site (Abruzzo, Central Italy): a new approach to a clactonian and acheulian sequence. *Quat. Int.* **409**, 182–194. https://doi.org/10.1016/j.quaint.2015.08.080 (2016).

**Acknowledgements**

The fieldwork in the Anagni basin has been carried out under the formal permission provided to the IIPU by the MIBACT—*Ministero per i Beni e le Attività Culturali e per il Turismo, Direzione Generale Archeologia, Belle Arti e Paesaggio* and by the *Soprintendenza Archeologia, Belle Arti e Paesaggio per le province di Frosinone, Latina e Rieti* (DG-ABAP/07-05-2020/622). The new NGX 600 mass spectrometer installed in November 2019 at the LSCE Facility and dedicated to 40Ar/39Ar geochronology was funded by a SESAME Ile de France grant to S. Nomade as well as additional funding from the CNRS and CEA Saclay. This research has been financially supported by the DTD-MIUR NextData project. Drilling at Colle Marino was performed by C.N.G. S.r.l., Roma. We thank Massimo Lenoci for technical support. Authors are also grateful to Mr. Enrico Missori, owner of the Colle Marino locality, for his kindness.

**Author contributions**

Fa.F., F.M., designed the geological study, performed stratigraphic investigations, and wrote the paper; I.B., S.G., D.A., F.P., L.B., Fe.F. designed the archaeological study and performed stratigraphic investigations; A.P, S.N, B.R.J., H.G. performed the 40Ar/39Ar age determinations; Fa.F., Fe.F., P.M., performed the paleomagnetic analyses; M.G. performed the petrographic analyses; C.M. performed the sedimentological and biostratigraphic analysis. All authors contributed to the writing of the text, discussed the results of the study, and commented on the manuscript.

**Funding**

Funding was provided by PNR 2011-2013.

**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-021-85446-5.

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

**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_14.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) 2021