# Striking forest revival at the end of the Roman Period in north-western Europe

C. Lambert, A. Penaud, M. Vidal, C. Gandini, L. Labeyrie, Laurent

Chauvaud, A. Ehrhold

► To cite this version:

C. Lambert, A. Penaud, M. Vidal, C. Gandini, L. Labeyrie, et al. Striking forest revival at the end of the Roman Period in north-western Europe. Scientific Reports, 2020, 10 (1), 10.1038/s41598-020-77253-1. hal-03126799

HAL Id: hal-03126799

https://hal.univ-brest.fr/hal-03126799

Submitted on 1 Feb 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_1.png" /></figure>

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

Check for updates

OPEN

# Striking forest revival at the end of the Roman Period in north-western Europe

C. Lambert<sup>1,2✉</sup>, A. Penaud<sup>1</sup>, M. Vidal<sup>1</sup>, C. Gandini<sup>3,4</sup>, L. Labeyrie<sup>2</sup>, L. Chauvaud<sup>5</sup> & A. Ehrhold<sup>6</sup>

The Holocene period (last 11,700 years BP) has been marked by significant climate variability over decadal to millennial timescales. The underlying mechanisms are still being debated, despite ocean–atmosphere–land connections put forward in many paleo-studies. Among the main drivers, involving a cluster of spectral signatures and shaping the climate of north-western Europe, are solar activity, the North Atlantic Oscillation (NAO) varying atmospheric regimes and North Atlantic oceanic gyre dynamics. Over the last 2500 years BP, paleo-environmental signals have been strongly affected by anthropogenic activities through deforestation and land use for crops, grazing, habitations, or access to resources. Palynological proxies (especially pollen grains and marine or freshwater microalgae) help to highlight such anthropogenic imprints over natural variability. Palynological analyses conducted in a macro-estuarine sedimentary environment of north-western France over the last 2500 years BP reveal a huge and atypical 300 year-long arboreal increase between 1700 and 1400 years BP (around 250 and 550 years AD) that we refer to as the ‘1.7–1.4 ka Arboreal Pollen rise event’ or ‘1.7–1.4 ka AP event’. Interestingly, the climatic 1700–1200 years BP interval coincides with evidence for the withdrawal of coastal societies in Brittany (NW France), in an unfavourable socio-economic context. We suggest that subpolar North Atlantic gyre strengthening and related increasing recurrence of storminess extremes may have affected long-term coastal anthropogenic trajectories resulting in a local collapse of coastal agrarian societies, partly forced by climatic degradation at the end of the Roman Period.

The Holocene is a complex puzzle of climate variability<sup>1–11</sup>, characterized at a millennial-scale by a 2500-year solar-related cyclicity over the whole period, and the appearance of a 1500-year ocean-related cyclicity over the last 5000 years, i.e. during the Mid- to Late-Holocene<sup>12</sup>. Gyre dynamics in the Atlantic Ocean are considered as a primary mechanism controlling rapid climate changes, especially subpolar gyre (SPG) dynamics which modulate northward heat transport through the North Atlantic Current (NAC) as well as deep-water formation and sea-surface temperatures (SSTs) in the North Atlantic Basin<sup>8–10</sup> (Fig. 1). Ocean surface density differences reconstructed south of Iceland and foraminifera shell isotopic analyses conducted in a core from the northern Bay of Biscay (Sites 2 and 5, Fig. 1) provide evidence of SPG strengthening intervals<sup>9,11</sup> (Fig. 2C,D), which led to a greater northward penetration of the warmer NAC in the northern Atlantic Ocean to the Nordic Seas<sup>13–17</sup>. Moreover, the North Atlantic Oscillation (NAO) pattern, defined as the difference between sea-level pressure over the Azores and Iceland, exerts strong control on the latitudinal position and vigour of the westerly wind belt<sup>18</sup>. An ‘NAO-like’ mechanism has been commonly evoked to explain paleoenvironmental evidence of climate variability over the Holocene<sup>11,16,19,20</sup>. Drift-ice discharge events<sup>1</sup> corresponding to SPG strengthening phases<sup>9,21</sup> would then correspond to recurrent positive modes of the NAO, as deduced from increasing storminess and moisture over northern Europe<sup>16,19,22,23</sup>.

North-western European coastal environments are interesting case studies because of their connection to the main Atlantic atmospheric and oceanic patterns that govern regional and global climates over different timescales. In this context, palynological studies can provide essential information about past environmen- tal conditions through combined analyses of phytoplanktonic microalgae and pollen grains. High-resolution

<footnote n="1"><sup>1</sup>CNRS, UMR 6538 Laboratoire Géosciences Océan (LGO), Univ. Brest (UBO), 29280 Plouzané, France. <sup>2</sup>UMR 6538 Laboratoire Géosciences Océan (LGO), Univ. Vannes (UBS), 56000 Vannes, France. <sup>3</sup>Centre de Recherche Bretonne et Celtique (CRBC), Univ. Brest (UBO), Brest, France. <sup>4</sup>AOROC-UMR 8546, CNRS, École Normale Supérieure, Paris, France. <sup>5</sup>CNRS, IRD, UMR 6539 Laboratoire des Sciences de l'Environnement Marin (LEMAR), Univ. Brest (UBO), 29280 Plouzané, France. <sup>6</sup>Géosciences Marines, Centre de Brest, IFREMER, Plouzané, France. ✉email: clement.lambert@univ-ubs.fr</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Map of Sea Surface Temperatures showing the schematic surface circulation of the North Atlantic Ocean. Red arrows characterize the main branches of the North Atlantic Current (NAC), namely the Irminger Current (IC) and the Norwegian Current (NC), bringing Subpolar Gyre (SPG) and Subtropical Gyre (STG) waters into the northern latitudes of the North Atlantic Basin. Blue arrows characterize colder surface waters, the East Greenland Current (EGC), the West Greenland Current (WGC) and the Labrador Current (LC). Sediment cores discussed in this study are represented by red dots: (1) study cores “G” and “KS-02” in the Bay of Brest, (2) core CBT-CS1 in the northern Bay of Biscay²¹, (3) core KV14bis at the Loire estuary mouth³³, (4) Cueva de Asiul speleothem isotopic record in Northern Spain³⁴, (5) core RAPID-12-1K south of Iceland⁹. The map was created using QGIS software version 3.4.5 (https://www.qgis.org).</figcaption></figure>

coastal records available close to continental sources such as watershed outlets are suitable archives for exploring climate-society relationships through time. Indeed, human-environment interactions remain poorly understood for pre-industrial periods, although the study of these times would allow us to put contemporary changes into perspective and thus improve understanding of underlying natural and/or anthropogenic forcings. The purpose of the present study is thus to identify past environmental and climate changes in a macro-estuarine sedimentary environment of north-western France (north-eastern Atlantic Ocean) over the last 2500 years BP and to discuss their potential implications for the dynamics of coastal societies. This study is relevant for north-western Europe but also for the wider understanding of long-term socio-ecosystemic trajectories.

## Results and discussion

Two highly temporally resolved sediment cores retrieved from the Bay of Brest (BB; Site 1, Fig. 1; Supplementary information Fig. 1) are used to infer crossed land-sea changes regarding vegetation dynamics on watersheds and sea-surface changes over the last 2500 years BP. We analysed pollen and dinoflagellate cyst (dinocyst) assemblages on a composite recording, including the 'G' core (from 2530 to 1660 years BP; 48°19'14" N; 4°23'5" W; 7.4 m depth; 'Défis Golfe de Gascogne' cruise) and 'KS-02' core (from 1660 to 500 years BP; 48°18'46" N; 4°24'27" W; 8 m depth; 'EssCalico' cruise), with a temporal sampling resolution of about 35 years (Supplementary information Fig. 2). Over the last 2500 years BP, pollen results from the BB record show a great variability, with tree pollen percentages ranging between 45 and 90% (Fig. 2G). While tree percentages of 45% characterize the average values recorded in modern BB sediments²⁴ (dotted line on Fig. 2G), the maximal values reached between 1700 and 1400 years BP led us to refer to this interval as the '1.7–1.4 Arboreal Pollen rise Event' or '1.7–1.4 ka AP event' (green band on Fig. 2G). Tree percentages exceed even those values recorded close to the studied area during the Mesolithic at around 9000 years BP²⁵ (~ 80%; dotted line on Fig. 2G), which is a period characterized by few or no major human-forced environmental changes. The '1.7–1.4 ka AP event' is all the more striking as it takes place in a major increasing phase of landscape opening (detectable in pollen sequences since around 4000 years BP¹¹,²⁶). This atypical arboreal pollen increase, only noticed in a few records¹¹,²⁶,²⁷, thus raises the question of its climate and/or anthropogenic origin (i.e. agricultural abandonment, probably leading to a forest revival).

Between 1700 and 1400 years BP, the mixed oak forest peak, mainly supported by *Corylus* in the BB (Fig. 2H), could be partly attributed to moister and milder conditions in north-western Europe. During the 1.7–1.4 ka AP event, humid conditions are supported by increases in riparian pollen percentages as *Alnus* and *Salix* (Fig. 2I),
the latter being recently discussed as powerful fluvial discharge tracers in the studied area<sup>11,21,24,28</sup>. Temporary increases of total pollen fluxes during the climatic interval from 1700- to 1200 years BP also testify to the increase in fluvial terrigenous inputs, which further enhanced the transport of pollen grains (Fig. 2K). Continental evidence of moisture increase is also corroborated by dinocyst observations (Fig. 2J). *Lingulodinium machaerophorum*, an estuarine-proliferating species<sup>29,30</sup> that tolerates strong drops in salinity<sup>29,31</sup>, testifies to coastal stratified waters subject to strong continental influence<sup>11,21</sup> and, thus, to stronger fluvial discharges at the studied site. Humid conditions together with an increase of fluvial inputs may explain the *Corylus* pollen rise during the 1.7–1.4 ka AP event in the BB. Moreover, this interval is also characterized by wetter conditions over northern central Europe then reflecting an eastward extension of the climatic conditions reconstructed in the BB<sup>32</sup>.

Recent studies conducted over south-western Europe allow us to enlarge our observations. In the Bay of Biscay (Site 3, Fig. 1), a large fall in Titanium-XRF counts (Ti-XRF; Fig. 2F) was observed between 1700 and 1400 years BP in a sediment core retrieved at the outlet of the Loire River<sup>33</sup>. This signal was interpreted as a result of a smaller contribution of Loire River discharges and thus of lower runoff over Loire watersheds<sup>11,33</sup>. The decrease in markers of fluvial influence recorded in a marine core from the northern Bay of Biscay (Site 2, Fig. 1) also testify to the decrease of the Loire contribution<sup>11,21</sup> during the 1.7–1.4 ka AP event in the BB. Furthermore, between 2700 and 800 years BP, precipitation quantifications carried out in speleothems of northern Spain<sup>34</sup> (Site 4, Figs. 1; 2E; Cueva de Asiul) show comparable variations to those of Loire paleo-fluvial discharges with a decrease of moisture between 1700 and 1300 years BP<sup>35</sup>. An 'aridification event' was also detected in the southwestern Mediterranean area, characterized by forest retraction and a decrease in fluvial inputs<sup>35-37</sup>, suggesting similar hydro-climatic influences on southern European sites. We suggest that the obvious correlation between higher humidity indexes in the BB (tree pollen percentages, especially *Corylus* and *Alnus*, as well as *L. machaerophorum* dinocysts; Fig. 2G,I,J) and higher aridity indexes to the south (lower Loire terrigenous inputs to the sea, Fig. 2F; Mediterranean aridity markers<sup>35-37</sup>), could attest to the presence of a differential north-south behaviour similar to the one described by the NAO pattern today, responsible for contrasting weather and precipitation regimes over Europe<sup>18,38</sup>. After 800 years BP (black dotted line on Fig. 2), the Ti-XRF signature from the Loire marine core (Fig. 2F) and precipitations reconstructed in northern Spain (Fig. 2E) show the reverse pattern compared with that described above. From this threshold, the Ti-XRF terrigenous detrital proxy for Loire watersheds seems to become a marker of soil erosion caused by massive human deforestation following the Medieval Climatic Optimum<sup>11</sup>, rather than a paleo-fluvial proxy in a context of enhanced precipitations. Tree percentages in BB watersheds also sharply decline at that time and reach values below the average of the 2500 year-long dataset, pointing to a pronounced anthropogenic limit for the Late Holocene in the studied area (green dotted line on Fig. 2G).

These data therefore allow us to consider a climatic hypothesis for the 1.7–1.4 ka AP event. We assume that a reinforced SPG and the subsequent intensification of the NAC, promoting northward oceanic heat transport and SPG-related internal climate feedback (e.g., reduced winter sea-ice coverage, near-surface atmosphere warming, negative sea-level pressure leading to a northward positioning of the westerly belt; Fig. 3), are the main drivers for the warm, moist winters recorded at the end of the western European Roman Period, thus explaining the temperate forest growth.

In the climatic context of increasing fluvial discharges and pollen fluxes (Fig. 2K,L), the major peak of trees (exceeding the Mesolithic rate) is, moreover, observed synchronously with a sharp decrease in cultivated plants (*Cerealia*-type pollen grain concentrations, Fig. 2N, as well as anthropogenic pollen indicators, API, Fig. 2M; Supplementary Fig. S3). We thus hypothesize that the 1.7–1.4 ka AP event (Fig. 2, green band) may also reflect a local collapse of the agrarian system allowing reforestation of wetlands and abandoned agricultural plots especially around 1500 years BP. The end of the Roman Period was also perturbed by political, economic, social and military unrests<sup>39,40</sup>, which could have been locally exacerbated by the increased storminess known on a broader coastal European-scale, during the 1700–1200 years BP climatic interval<sup>23</sup> (Fig. 2A) or the larger 1900–1050 years BP interval<sup>19</sup> (Fig. 2B). From an archaeological point of view, fine ceramics from this period are rare and little known coupled with poor conservation of late levels. This bias can indeed hamper identification of the occupation phases of the late Roman Period<sup>41</sup>. For the best documented regions in north-western Europe (e.g., the centre of France<sup>42,43</sup>), the fifth century (i.e., after 1600 years BP) shows a sharp reduction in the number of occupied sites, while the creation of settlements seems to represent a marginal phenomenon. Nevertheless, the interpretation of this trend is still widely discussed as the rate and extent of withdrawals can vary considerably by region. Moreover, studies have shown changes since the beginning of the second century (i.e. around 1850 years BP), well before the main upheavals of the 3rd–4th (i.e. between 1750 and 1500 years BP) centuries, urging caution regarding the explanation of military disturbances and suggesting rather a gradual reorganisation of land use and the production system<sup>42,44</sup>. Thus, we suggest that the climatic degradation, which could for example have caused significant floods, could have been an amplifying factor in the collapse of the otherwise struggling agrarian society in the region.

A combination of oceanic (i.e., SPG and NAC strengthening) and atmospheric configurations (i.e., northward position of the westerly belt), strongly influence European climate over the Late Holocene (Fig. 3). The 1.7–1.4 ka AP event may result from both greater humidity and major phases of human withdrawals in a general context of regional climatic deterioration. Our pollen records thus agree with a 300-year-long collapse of coastal agrarian societies recorded in NW France and probably partly forced by climatic variability. At 800 years BP, a resumption of deforestation activities contributed to a landscape opening rate similar to that reconstructed before the 1.7–1.4 ka AP event.
◄**Figure 2.** Comparison of palynological and sedimentological signals from the composite "G-KS02" Bay of Brest sequence (G–N) with paleoclimatic records from the North Atlantic: (A) European Atlantic storm events (EASEs) over European Atlantic coasts²³, (B) Holocene storm periods (HSPs) over northern Europe¹⁹, (C) water density differences for SPG strength⁹, (D) δ18O on benthic foraminifera shells (A. falsobeccarii) from the core CBT-CS11 (Northern Bay of Biscay) for SPG strength²¹, (E) detrended and combined isotopic record from two Cueva de Asiul speleothems, Northern Spain³⁴, (F) Ti-XRF record from core KV14bis at the Loire estuary mouth³³, and main results from the composite "G-KS02" Bay of Brest sequence (this study): (G) sum of tree pollen percentages, (H) *Corylus* percentages, (I) sum of riparian tree (*Alnus*, *Salix* and *Betula*) percentages, (J) percentages of the dinocyst *Lingulodinium machaerophorum*, (K) total pollen grain fluxes, (L) sedimentation rates, (M) percentages of Anthropogenic Pollen Indicators (API: Asteraceae, Brassicaceae, *Fagopyrum*, *Plantago lanceolata*, Poaceae, *Polygonum aviculare*, *Rumex*) with (black) and without (grey) *Cerealia*-type, (N) *Cerealia*-type pollen concentrations. Light grey bands highlight intervals likely characterized by "NAO + like" atmospheric configuration (i.e., enhanced storminess and northward positioning of the westerly belt). The green band highlights the atypical forest cover increase in the Bay of Brest watersheds between 1700 and 1400 years BP. Dotted black line and grey arrows underline the divergence observed since 800 years BP between higher Loire runoff and lower precipitations reconstructed from northern Spain, likely pointing to an anthropogenic-related runoff forcing at that time.

<figure><img src="image_4.png" /><figcaption>**Figure 3.** Schematic conceptual models showing storm tracks positioning, North Atlantic Current (NAC) and Subpolar Gyre (SPG) vigor during (A) 'NAO+like' and (B) 'NAO- like' pluri-secular intervals and the associated influence on storms and precipitation patterns over North West Europe. Watersheds in blue represent increased precipitation and runoff. The map was created using QGIS software version 3.4.5 (https://www.qgis.org).</figcaption></figure>

## Methods

**Stratigraphy of the composite G-KS-02 study core.** Both studied cores (Supplementary Table S1) were taken in the Bay of Brest macro-estuary (NW France; Supplementary Fig. S1) onboard the RV 'Côtes de la Manche' in the framework of (i) the Ifremer program 'Défis Golfe de Gascogne' (2003, vibrocorer) for G core and (ii) the 'EssCALICO' cruise (2010, gravity corer) for the Ks-02 core. Descriptions, photographs and X-ray radiographies are available (Supplementary Fig. S1). A total of 14 radiocarbon dates were acquired on gastropods (*Turritella* spp., *Caliostoma zizyphinum*, *Bittium reticulatum*) for the G core (six ¹⁴C-AMS dates) and the KS-02 core (eight 14C-AMS dates). All 14C-AMS dates (Supplementary Table S2) were calibrated to calendar years with the CALIB 7.1 program using the Marine13 calibration curve⁴⁵ which considers a 400-year correction for the mean ocean surface reservoir age. Reservoir age can, however, vary spatially and temporally. The local deviation from the oceanic mean (ΔR) is estimated at –40 ± 23 years off the Brittany peninsula⁴⁶, as confirmed by a reservoir age test in the BB²⁵. We applied this additional correction to all dating obtained on marine carbonate material. Thus, all the dates were calibrated and are expressed hereafter in calendar age (Cal.) BP ('before present', before 1950).

The age model produced on the KS-02 core confirms that this sedimentary sequence provides information on the continuity of the very closed G core (whose upper part was missing), thus making it possible to constitute a long composite sequence, as illustrated by the continuity in the pollen data (Supplementary Fig. S3). The composite age model was established with the 'Bacon' age-depth modelling package⁴⁷ in R software (Supplementary Fig. S2) and makes it possible to document a period of around 2000 years (from 2528 to 488 years Cal. BP). In the manuscript, when referring to a precise age, 'Cal.' will not be systematically written, and ages will often be
specified as 'year BP' or 'year BC' to facilitate the reading of results by paleoenvironmental or archaeological
communities, respectively.

**Palynological analyses.** In this study, we analysed 58 samples (27 for core G, 31 for core KS-02). Palynological preparations on the <150 µm sediment fraction were carried out at the EPOC laboratory (University of Bordeaux, France). The procedure is based on chemical (cold 10–25–50% HCl and cold 40–70% HF) treatments to remove siliceous and carbonate fractions, as well as physical sieving (10-µm nylon mesh screen) to remove the mineral clayey fractions and concentrate organic micro-fossils (i.e., palynomorphs)⁴⁸. Final residues were mounted with glycerine between a slide and a coverslip. Pollen and dinocysts were identified using an optical microscope Leica DM2500 at X630 magnification⁴⁹⁻⁵¹.

For each sample, a minimum of 300 pollen grains and 150 dinocysts were systematically counted (Supplementary Table S3) in order to provide robust assemblages from a statistical point of view⁵². Percentages were calculated on a total pollen or total dinocyst sum with no taxon exclusion, and concentrations were expressed in number of palynomorphs/cm³ using the 'marker grain method', which consists in adding a known number of exotic spores (*Lycopodium clavatum*) to dried sediments before any palynological treatments⁵³.

**Data availability**

The datasets generated during this study are partly included in this published article and its supplementary
information and are available from the corresponding author upon reasonable request.

Received: 25 August 2020; Accepted: 3 November 2020

Published online: 15 December 2020

**References**

1. Bond, G. et al. Persistent solar influence on North Atlantic climate during the Holocene. *Science* **294**(5549), 2130–2136 (2001).

2. Mayewski, P. A. et al. Holocene climate variability. *Quat. Res.* **62**(3), 243–255 (2004).

3. Wanner, H., Solomina, O., Grosjean, M., Ritz, S. P. & Jetel, M. Structure and origin of Holocene cold events. *Quat. Sci. Rev.* **30**(21–22), 3109–3123 (2011).

4. Debret, M. et al. Evidence from wavelet analysis for a mid-Holocene transition in global climate forcing. *Quat. Sci. Rev.* **28**(25–26), 2675–2688 (2009).

5. Steinhilber, F., Beer, J., & Fröhlich, C. Total solar irradiance during the Holocene. *Geophys. Res. Lett.* **36**, L19704. https://doi.org/10.1029/2009GL040142 (2009).

6. Olsen, J., Anderson, N. J. & Knudsen, M. F. Variability of the North Atlantic Oscillation over the past 5,200 years. *Nat. Geosci* **5**(11), 808–812 (2012).

7. Morley, A., Rosenthal, Y. & DeMenocal, P. Ocean-atmosphere climate shift during the mid-to-late Holocene transition. *Earth Planet. Sci. Lett.* **388**, 18–26 (2014).

8. Hátűn, H., Sandø, A. B., Drange, H., Hansen, B. & Valdimarsson, H. Influence of the Atlantic subpolar gyre on the thermohaline circulation. *Science* **309**(5742), 1841–1844 (2005).

9. Thornalley, D. J., Elderfield, H. & McCave, I. N. Holocene oscillations in temperature and salinity of the surface subpolar North Atlantic. *Nature* **457**(7230), 711–714 (2009).

10. Moffa-Sánchez, P. & Hall, I. R. North Atlantic variability and its links to European climate over the last 3000 years. *Nat. Commun.* **8**(1), 1–9 (2017).

11. Penaud, A. et al. Oceanic versus continental influences over the last 7 kyrs from a mid-shelf record in the northern Bay of Biscay (NE Atlantic). *Quat. Sci. Rev.* **229**, 106135 (2020).

12. Debret, M. et al. The origin of the 1500-year climate cycles in Holocene North-Atlantic 239 records. *Clim. Past Discuss.* **3**, 679–692 (2007).

13. Giraudeau, J. et al. Millennial-scale variability in Atlantic water advection to the Nordic Seas derived from Holocene coccolith concentration records. *Quat. Sci. Rev.* **29**(9–10), 1276–1287 (2010).

14. Miettinen, A., Divine, D., Koç, N., Godtliebsen, F. & Hall, I. R. Multicentennial variability of the sea surface temperature gradient across the subpolar North Atlantic over the last 2.8 kyr. *J. clim.* **25**(12), 4205–4219 (2012).

15. Moreno-Chamarro, E., Zanchettin, D., Lohmann, K. & Jungclaus, J. H. An abrupt weakening of the subpolar gyre as trigger of Little Ice Age-type episodes. *Clim. Dyn.* **48**(3–4), 727–744 (2017).

16. Goslin, J. et al. Holocene centennial to millennial shifts in North-Atlantic storminess and ocean dynamics. *Sci. rep.* **8**(1), 1–12 (2018).

17. Dubois-Dauphin, Q., Colin, C., Elliot, M., Dapoigny, A. & Douville, E. Holocene shifts in sub-surface water circulation of the North-East Atlantic inferred from Nd isotopic composition in cold-water corals. *Mar. Geol.* **410**, 135–145 (2019).

18. Hurrell, J. W. & Deser, C. North Atlantic climate variability: The role of the North Atlantic Oscillation. *J. Mar. Syst.* **79**(3–4), 231–244 (2010).

19. Sorrel, P. et al. Persistent non-solar forcing of Holocene storm dynamics in coastal sedimentary archives. *Nat. Geosci.* **5**(12), 892–896 (2012).

20. Mojtahid, M. et al. High resolution Holocene record in the southeastern Bay of Biscay: Global versus regional climate signals. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **377**, 28–44 (2013).

21. Penaud, A. et al. A new midshelf record in the northern Bay of Biscay (NE Atlantic, CBT-CS11 core): Sedimentological, geochemical and palynological data over the last 7 kyrs. *Data Brief* **29**, 105323–105323 (2020).

22. Sorrel, P., Tessier, B., Demory, F., Delsinne, N. & Mouazé, D. Evidence for millennial-scale climatic events in the sedimentary infilling of a macrotidal estuarine system, the Seine estuary (NW France). *Quat. Sci. Rev.* **28**(5–6), 499–516 (2009).

23. Pouzet, P. et al. Chronology of Holocene storm events along the European Atlantic coast: New data from the Island of Yeu, France. *Prog. Phys. Geogr. Eart Environ.* **42**(4), 431–450 (2018).

24. Lambert, C. et al. Modern palynological record in the Bay of Brest (NW France): Signal calibration for palaeo-reconstructions. *Rev. Palaeobot. Palynol.* **244**, 13–25 (2017).

25. Lambert, C. et al. Palaeoenvironmental reconstructions during the Meso-to Neolithic transition (9.2–5.3 cal. ka BP) in North-western France: Palynological evidences. *Holocene* **29**(3), 380–402 (2019).

26. Fernane, A. et al. Coastal palaeoenvironmental record of the last 7 kyr BP in NW France: Sub-millennial climatic and anthropic Holocene signals. *Holocene* **24**(12), 1785–1797 (2014).

27. Barbier, D. & Visset, L. Logné, a peat bog of European ecological interest in the Massif Armorican, western France: Bog development, vegetation and land-use history. *Veg. Hist. Archaeobot.* **6**(2), 69–77 (1997).
28. Lambert, C. et al. Human-induced river runoff overlapping natural climate variability over the last 150 years: Palynological evidence (Bay of Brest, NW France). *Glob. Planet. Chang.* **160**, 109–122 (2018).

29. Morzadec-Kerfourn, M. T. Les kystes de dinoflagellés dans les sédiments récents le long des côtes Bretonnes. *Rev. Micropaleontol.* **20**(3), 157–166 (1977).

30. Ganne, A., Leroyer, C., Penaud, A. & Moj tahid, M. Present-day palynomorph deposits in an estuarine context: The case of the Loire Estuary. *J. Sea Res.* **118**, 35–51 (2016).

31. Reid, P. C. A regional sub-division of dinoflagellate cysts around the British Isles. *New Phytol.* **75**(3), 589–603 (1975).

32. Büntgen, U. et al. 2500 years of European climate variability and human susceptibility. *Science* **331**(6017), 578–582 (2011).

33. Durand, M. et al. Late Holocene record from a Loire River incised paleovalley (French inner continental shelf): Insights into regional and global forcing factors. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **511**, 12–28 (2018).

34. Smith, A. C. et al. North Atlantic forcing of moisture delivery to Europe throughout the Holocene. *Sci. Rep.* **6**, 24745 (2016).

35. Fletcher, W. J., Debret, M. & Goñi, M. F. S. Mid-Holocene emergence of a low-frequency millennial oscillation in western Mediterranean climate: Implications for past dynamics of the North Atlantic atmospheric westerlies. *Holocene* **23**(2), 153–166 (2013).

36. Zielhofer, C. et al. Atlantic forcing of Western Mediterranean winter rain minima during the last 12,000 years. *Quat. Sci. Rev.* **157**, 29–51 (2017).

37. Di Rita, F. et al. Holocene forest dynamics in central and western Mediterranean: Periodicity, spatio-temporal patterns and climate influence. *Sci. Rep.* **8**(1), 1–13 (2018).

38. Tréguer, P. et al. Large and local-scale influences on physical and chemical characteristics of coastal waters of Western Europe during winter. *J. Mar. Syst.* **139**, 79–90 (2014).

39. France, J., & Delaplace, C. *Histoire des Gaules* (VIe siècle av. J.-C.-VIIe siècle ap. J.-C.), Coll. Cursus. Histoire, Armand Colin, 5e edition, Paris (2016).

40. *Mondes ruraux en Orient et en Occident- I et II, 20-21*, Brepols Publishers, Turnhout, (2012–2013).

41. Ouzoulias, P., Pellecuer, C., Raynaud, C., Van Ossel, P., & Garmy, P. *Les campagnes de la Gaule à la fin de l'Antiquité*, Actes du IVe colloque Ager, Montpellier, 11–14 mars 1998, Ed. APDCA (2001).

42. Gandini, C. *Des campagnes gauloises aux campagnes de l'Antiquité tardive: la dynamique de l'habitat rural dans la cité des Bituriges Cubi* (IIe s. av. J.-C.-VIIe s. ap. J.-C.). 33e suppl. à la RACF, FERACF, Tours. (2008).

43. Gandini, C., Favory, F., & Nuninger, L. dir. *Seven millenia of territorial dynamics: settlement pattern, production and trades from Neolithic to Middle Ages, ArchaeDyn*, colloque de Dijon, 23–25 juin 2008, British Archaeological Reports, BAR 2370, 90 (2012).

44. Gandini, C. Mettre à l'épreuve les typologies basées sur des données de prospection: la caractérisation de l'habitat rural de la fin de l'Antiquité dans le centre de la Gaule, in *Tobalina Pulido, L., Campo, A., & Cabes, L. Traiter les données archéologiques tardo-antiques. Approches, methodes et traitement de l'information*, actes du séminaire organisé à l'université de Pau et des Pays de l'Adour, 12 avril 2019, *Latomus*, (in Press).

45. Reimer, P. J. et al. IntCal13 and Marine13 radiocarbon age calibration curves 0–50,000 years cal BP. *Radiocarbon* **55**(4), 1869–1887 (2013).

46. Mangerud, J., Bondevik, S., Gulliksen, S., Hufthammer, A. K. & Høisæter, T. Marine 14C reservoir ages for 19th century whales and molluscs from the North Atlantic. *Quat. Sci. Rev.* **25**(23–24), 3228–3245 (2006).

47. Blaauw, M. & Christen, J. A. Flexible paleoclimate age-depth models using an autoregressive gamma process. *Bayesian Anal.* **6**(3), 457–474 (2011).

48. de Vernal, A., Henry, M. & Bilodeau, G. Techniques de préparation et d'analyse en micropaléontologie. *Les cahiers du GEOTOP* **3**, 41 (1999).

49. Faegri, K. & Iversen, J. *Textbook of Pollen analysis* (Wiley, New-York, 1989).

50. Reille, M. *Pollen et spores d'Europe et d'Afrique du Nord* (Laboratoire de Botanique historique et Palynologie, Marseille, 1995).

51. Zonneveld, K. A. F. & Pospelova, V. A determination key for modern dinoflagellate cysts. *Palynology* **39**(3), 387–409 (2015).

52. Fatela, F. & Taborda, R. Confidence limits of species proportions in microfossil assemblages. *Mar. Micropaleontol.* **45**(2), 169–174 (2002).

53. Mertens, K. N. et al. Determining the absolute abundance of dinoflagellate cysts in recent marine sediments: The Lycopodium marker-grain method put to the test. *Rev. Palaeobot. Palynol.* **157**(3–4), 238–252 (2009).

## Acknowledgements

This study was supported by the French CNRS and is a contribution to the 2015-2016 INSU project EC2CO-LEFE “CAMOMI: Convergences/Approches croisées des signaux MOléculaires et Ml cripaléontologiques pour décrypter les forçages anthropiques et climatiques en milieu côtier (Rade de Brest)” and the 2014-2015 UBO-BQR project: “PARADE: Signature PAléoenvironnementale des séquences Holocène en RADE de Brest : forçages climatiques et anthropiques”. This work was supported by the Laboratoire d'Excellence LabexMER (ANR-10-LABX-19) and co-funded by a grant from the French government under the program Investissements d'Avenir. We thank the UMR CNRS 5805 EPOC (Talence) for palynological laboratory procedures (Muriel Georget), the UMS 2572 LMC14 (Saclay) for carbon dating via Artemis facilities, the laboratory IFREMER-Marine Geosciences (Plouzané) for the KS-02 core availability (Esscalico cruise; 2010), and the LEMAR (*Laboratoire des Sciences de l'Environnement Marin*, IUEM, Plouzané) for the “G core” availability (Ifremer program ‘Défis Golfe de Gascogne’, 2003). Main issues of this project are integrated within the theme ‘Dynamics of Human Settlement and Paleoenvironments’ of the Zone Atelier Brest Iroise (ZABRI, INEE-CNRS). Finally, we would like to thank the *Bureau de traduction* of the University of Brest for the improvement of English.

## Author contributions

We confirm that all authors have approved the final version of the manuscript. All authors made substantial contributions to the submission. Details of individual contributions: C.L. (sampling, core analyzes, palynological analyzes, interpretation of data, creation of figures, writing the original draft); A.P. (Palynologist: funding acquisition, coring and core analyzes, interpretation of data, reading and revision); M.V.: (Palynologist: interpretation of data, reading and revision); C.G. (Historian: writing the archaeological context, reading and revision); L.L. (Climatologist: reading and revision); L.C.: (Biologist: coring, reading and revision); A.E. (Sedimentologist: head of coring campaign, reading and revision).

## Competing interests

The authors declare no competing interests.
**Additional information**

**Supplementary information** is available for this paper at https://doi.org/10.1038/s41598-020-77253-1.

**Correspondence** and requests for materials should be addressed to C.L.

**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_5.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) 2020