Climate induced human demographic and cultural change in northern Europe during the mid-Holocene

Warden, L.; Moros, M.; Neumann, T.; Shennan, S.; Timpson, A.; Manning, K.; Sollai, M.; Wacker, L.; Perner, K.; Häusler, K.; Leipe, T.; Zillén, L.; Kotilainen, A.; Jansen, E.; Schneider, R. R.; Oeberst, R.; Arz, H.; Sinninghe Damsté, J. S.

doi:10.1038/s41598-017-14353-5

Download PDF

Article
Open access
Published: 10 November 2017

Climate induced human demographic and cultural change in northern Europe during the mid-Holocene

L. Warden¹,
M. Moros²,
T. Neumann²,
S. Shennan³,
A. Timpson⁴,
K. Manning⁵,
M. Sollai¹,
L. Wacker⁶,
K. Perner²,
K. Häusler²,
T. Leipe²,
L. Zillén⁷,
A. Kotilainen⁸,
E. Jansen⁹,
R. R. Schneider ORCID: orcid.org/0000-0003-1453-9181¹⁰,
R. Oeberst¹¹,
H. Arz ORCID: orcid.org/0000-0002-1997-1718² &
…
J. S. Sinninghe Damsté ORCID: orcid.org/0000-0002-8683-1854^1,12

Scientific Reports volume 7, Article number: 15251 (2017) Cite this article

6392 Accesses
67 Citations
28 Altmetric
Metrics details

Subjects

Abstract

The transition from hunter-gatherer-fisher groups to agrarian societies is arguably the most significant change in human prehistory. In the European plain there is evidence for fully developed agrarian societies by 7,500 cal. yr BP, yet a well-established agrarian society does not appear in the north until 6,000 cal. yr BP for unknown reasons. Here we show a sudden increase in summer temperature at 6,000 cal. yr BP in northern Europe using a well-dated, high resolution record of sea surface temperature (SST) from the Baltic Sea. This temperature rise resulted in hypoxic conditions across the entire Baltic sea as revealed by multiple sedimentary records and supported by marine ecosystem modeling. Comparison with summed probability distributions of radiocarbon dates from archaeological sites indicate that this temperature rise coincided with both the introduction of farming, and a dramatic population increase. The evidence supports the hypothesis that the boundary of farming rapidly extended north at 6,000 cal. yr BP because terrestrial conditions in a previously marginal region improved.

Human forager response to abrupt climate change at 8.2 ka on the Atlantic coast of Europe

Article Open access 02 May 2022

Asier García-Escárzaga, Igor Gutiérrez-Zugasti, … Patrick Roberts

Mid-to-late Holocene climate variability in coastal East Asia and its impact on ancient Korean societies

Article Open access 16 September 2023

Jungjae Park, Junbeom Bahk, … Jieun Choi

Quantitative comparison of geological data and model simulations constrains early Cambrian geography and climate

Article Open access 23 June 2021

Thomas W. Wong Hearing, Alexandre Pohl, … Thijs R. A. Vandenbroucke

Introduction

The transition from foraging to farming provided the basis for a massive increase in food production and human population density, which in turn drove the accumulation of technology and specialization, ultimately enabling the establishment of modern day societies¹. The reasons and mechanisms for this Neolithic Revolution occurring at different times in different places continues to be debated. In Northern Europe cereal agriculture and animal husbandry did not begin until ca. 6,000 cal. yr BP, i.e. 1,500 years later than in neighboring populations to the south^2,3,4,5,6. It has been proposed that foraging societies were successful in relying largely on the marine resources from the Baltic Sea and had no incentive to adopt an agrarian lifestyle⁵, which required more intensive labor with no immediate increase in living standards⁷. This leaves the question of why the transformation did eventually occur at 6,000 cal. yr BP. It has been suggested that climate change may have been responsible for the late onset of farming in the northern European regions²; however, due to the lack of high resolution climate reconstructions the effects of environmental change on human societies in these regions have been unclear. In order to examine this at a fine temporal scale, we employed an approach combining sediment proxy studies with 3D ecosystem modeling to reconstruct climate conditions in northern Europe and environmental conditions in the Baltic Sea and compare this with human population and subsistence changes from archeological sites from the region.

Results and Discussion

We examined key-sedimentary records from the central Baltic Sea (Fig. 1) from 7,100-3,000 cal. yr BP (Table 1). Complementary information for the entire Baltic was obtained from a number of other sediment cores (Figs 1, S1 and S2; Table 1). Our detailed age-depth model (Fig. S3) was developed based on radiocarbon dating of benthic foraminifera (Table S1) that are present in sediments from the deepest part of the Gotland Basin when inflowing saline waters from the North Sea reach the basin. This age-depth model was projected onto shallow water sites using LOI and/or organic carbon downcore profiles (see methods). The deep sites are, however, influenced by lateral transport of fine-grained material and, therefore, not suitable for sediment proxy reconstruction of past changes in in-situ surface water conditions. Therefore, a core from site 303600 (175 m water depth), located outside the inflowing water pathway on the western side of the eastern Gotland Basin (Fig. S1), was selected for determining the palaeotemperature record. This core was analyzed in high resolution using the TEX₈₆-paleothermometer, a proxy for reconstructing past sea surface temperature (SST) based on the distribution of Thaumarchaeotal membrane lipids⁸. We applied a local calibration of TEX₈₆ to reconstruct summer SST in the Baltic Sea⁹. Our record shows that summer SST varied between 14.5-17.5 °C during the study period (Fig. 2B). A prolonged unstable but generally colder phase is terminated at 6,000 cal. yr BP by a distinct and rapid warming phase (Fig. 2B). Two SST maxima at 5,600 and 4,500 cal. yr BP are evident during this regional Holocene Thermal Optimum lasting from ca. 5,900 to 4,400 cal. yr BP. Throughout the mid- to late-Holocene the general summer SST trend declines, following summer solar insolation (Fig. 2B).

Table 1 Details on the core sites studied located in almost all basins of the Baltic Sea. The location (see Figure S1 for a map), present-day and estimated paleo-water depths (6,500 cal. yr BP) are indicated, along with cruises and years of sampling of sites. Key-sites for chronology and proxy studies are marked in bold font. Loss on ignition (LOI) data are available from all sites. LOI data from cores marked with (X) are shown in Figs S2–S4.

Full size table

The present-day Baltic Sea is persistently stratified with anoxic bottom waters in the deep basins and oxygen-depleted (hypoxic) waters in coastal areas (so-called dead-zones^10,11). Summer SST is thought to have controlled hypoxia in the Baltic Sea over the last 1000 yrs through various feedback loops^9,12. A comparison of SST, the total organic carbon (TOC) content and the degree of lamination of the sediments indicates that the critical 16 °C temperature threshold for the occurrence of severe hypoxia proposed for the last 1000 yrs⁹, also acted earlier in the Holocene (Fig. 2). The marked rapid warming phase at 6,000 cal. yr BP co-occurs with the onset of bottom water hypoxia (Fig. 2). This intensification of hypoxia at 6,000 cal. yr BP is evident in the sedimentary records from the entire Paleo-Baltic (Fig. S4). The stable hypoxic conditions of the Thermal Optimum phase, when summer SST is persistently >16 °C, are widespread in the entire Baltic Sea (Fig. S4) and terminated by a marked cooling after 4,500 cal. yr BP when summer SST dropped below 16 °C (Fig. 2).

Our analysis using a coupled hydrodynamic-biogeochemical model applying a reconstructed paleo-bathymetry (Fig. S5) and forced by a delta change approach (Fig. S6) suggests that the SST change as recorded in the central Baltic Sea was not local and affected environmental conditions throughout the Baltic. The model provides a causal mechanism between temperature change and environmental conditions as reflected in the sedimentary record (laminated vs. homogenous sequences), and simulations under this model show that during the more unstable and rather cold phase before 6,000 cal. yr BP the hypoxic area (laminated sediments) in the Paleo-Baltic was restricted to the formerly deep northern basin and covered an area of only 20,000 km². After the rapid warming at 6,000 cal. yr BP and during the regional Thermal Optimum the hypoxic area increased (resulting in deposition of laminated sediments) to more than 100,000 km² and spread from the northern basins to the entire Baltic (Fig. 3A; Fig. S4). Therefore the simulation results support a theoretical study¹² arguing for a positive feedback mechanism between temperature and hypoxia (laminated sediments). Two factors are crucial for this positive feedback loop: i) blooms of N₂-fixing cyanobacteria which can only develop at summer SSTs >16 °C (refs^13,14), and ii) phosphate availability. Phosphorus mineralized in sediments under oxic (cold) conditions becomes to a certain extent bound in phosphorus-iron complexes and is retained in the homogenous sediment. A change to hypoxic/anoxic conditions as a result of temperature induced cyanobacteria blooms liberates the iron bound phosphorus and it is released as phosphate to the water column^15,16. This surplus phosphate in the water column in turn stimulates further cyanobacterial blooms resulting in a substantial increase in primary productivity during warm climatic phases. This additional production of biomass (in the modern central Baltic approximately half of the sedimentary organic carbon is estimated to be derived from these blooms¹⁷) results in an increased carbon flux into the bottom water and in turn the oxygen demand for heterotrophic processes increases. Bottom-water oxygen concentrations decline and increasingly larger areas become covered by anoxic waters (laminated sediments), resulting in even more phosphate release. Nutrient concentrations increase and as a final result of this positive feedback loop a eutrophic state of the Baltic Sea is established. Our simulations suggest that this state is almost independent of external nutrient loading and can only be changed by reducing summer SSTs (<16 °C), preventing heavy cyanobacterial blooms.

The regional temperature rise causing widespread environmental changes in the entire Baltic region must have also influenced terrestrial settings. Indeed, our Baltic Sea Holocene Thermal Maximum (6,000 to 4,500 cal. yr BP) is also seen in pollen based reconstruction of summer temperature for NW Europe¹⁸ and Finland¹⁹. We reconstructed fluctuations in the local human population through time using a Summed Probability Distribution (SPD) of 1,960 terrestrial radiocarbon dates from 608 archaeological sites (Fig. 4A), taken as a proxy for variation in human population levels in the Baltic region through time. The rationale behind this approach is that to a first order of approximation, periods of increased population will leave more dateable material behind in the archaeological record^3,20 (see SI for further discussion). Evidently, the large climatic and environmental change at 6,000 cal. yr BP that took place in northern Europe is concurrent with the start of a massive increase in population levels. Between 6,000 and 5,600 cal. yr BP the summer climate substantially increased by 2 °C in parallel with a 3.3 fold increase in population levels (Fig. 4). This represents a population growth rate of 6.1% per 25 year generation compared to 1.9% previously between 7,000 and 6,000 cal. yr BP, though these may be slight overestimates given the possible effect of taphonomic loss over time. For comparison, if Surovell’s proposed taphonomic correction formula²¹ is applied, these numbers fall slightly to a 3.0 fold increase, with growth rates of 5.7% and 1.5% per generation respectively. Additionally, the strong relationship between temperature and population is evident throughout the entire overlap in both time-series from 7,000 to 4,000 cal. yr BP (Pearson’s R = 0.73, p-value < 0.00002; see methods). Furthermore, cross-correlation analysis shows the population lagged behind climate by 60 years, suggesting a rapid human response to climate changes.

The mechanism for how the summer temperature increase at 6,000 cal. yr BP could cause a rapid population increase in this particular region can be understood in terms of the improved environmental conditions extending the growing season for domestic crops in a previously marginal area. Nearby regions to the south were already farming, and would have been able to rapidly expand their agrarian economy north as conditions improved. This results in a substantial increase in the available terrestrial biomass per unit area for human consumption, compared to a closed-canopy temperate forest²². The availability of more food resources then supports a greater birth rate²³, which in turn creates a greater demand for and dependency on food production. Our analysis of the radiocarbon record (Fig. 3C) is consistent with this, showing both a substantial population increase and a cultural transition to farming both occurring at 6,000 cal. yr BP.

Evidence of this economic revolution is seen in archaeobotanical and faunal evidence showing the simultaneous introduction of wheat, barley and domestic animals throughout southern Scandinavia ca. 6,000-5,700 cal. yr BP²⁴. This transition reflects an important cultural shift, from the Ertebølle culture, representing hunter-gatherer-fisher societies, to the agricultural Funnel Beaker culture (Figs 3C and 4B), though some exploitation of marine resources seem to have continued (see SI). Modeling studies indicate that a warming climate in Scandinavia could induce shifts in the northern limit of cereal suitability by 100–150 km per °C (ref.²⁵). Ancient DNA^26,27 and archaeological evidence²⁸ suggests that the transition occurred as farming communities to the south responded to the warming climate by expanding north into southern Scandinavia. By around 5,600 cal. yr BP population levels and summer temperatures reach a maximum (Fig. 4B) and farming is very extensive, with large areas of cleared land²⁹.

However, this still leaves the question of whether the improved conditions for farming produced by increasing temperatures at 6,000 cal. yr BP (Fig. 3B) were solely responsible for the observed change from foraging to farming, or whether other factors were at work. In addition to temperature changes, the development of widespread hypoxia in the whole Baltic Sea including coastal areas (Fig. 3A) may also have played a role in the rapid transition from the Ertebølle to the Funnel Beaker culture in the Baltic region, by having a negative effect on the availability of major marine resources. At the beginning of the early Neolithic, ~6,000 cal. yr BP, summer SSTs started to increase rapidly and bottom water conditions in wide areas became persistently hypoxic (Fig. 3A), a state that continued for most of the next two millennia during a local Thermal Optimum (when SST >16 °C, Figs 4 and S4). Hypoxia in the present-day Baltic Sea has been shown to contribute to a reduction in fish populations^19,30 (see also SI), therefore it is plausible that a development of widespread hypoxia at 6,000 cal. yr BP would have had a similar detrimental impact on fish stocks, particularly the reproductive volume for cod (Fig. S7 and SI). To test this we carried out an analysis of the relative proportions of remains of aquatic and terrestrial fauna from 96 coastal sites, which reveals substantial variation in the aquatic proportion across different site-phases (Fig. S8). A permutation test cannot reject the null hypothesis of no change through time, with the average proportion of aquatic species in the period 6,500 to 6,000 cal. yr BP of 35.6% (95% CI = 26.6% to 44.5%) remaining very similar in the period 6,000 to 5,500 cal. yr BP at 36.6% 995% CI = 22.9% to 50.0%), n = 96, p = 0.941. This analysis, therefore, does not support the hypothesis of an overall reduction in aquatic resource. Further work will be needed in the future to explore this apparent contradiction with the evidently increased hypoxia.

After the peak at 5,500 cal. yr BP, the SPD population proxy shows a steady decline throughout the remainder of the record (Fig. 4A), synchronous with the declining summer insolation at 60°N and summer SST (Fig. 4C), indicating that climatic deterioration played a role in the decline in human population, probably due in part to poorer conditions for farming. At this time there is evidence of woodland regeneration, which would result from a decrease in human activity²⁹. This period also sees the southward expansion of the Pitted Ware Culture²⁶ (Fig. 4B), during a period when environmental conditions were worsening further north².

Our data reveals that climate change played a key role in demographic and cultural changes leading to the appearance of farming in mid-Holocene southern Scandinavia after a 1,500 yr delay via rapid warming ca. 6,000 cal. yr BP which improved terrestrial conditions for farming. This incited an expansion of the agrarian groups from the south²⁶ resulting in the introduction of farming and as a consequence a massive population increase in the Baltic region. Subsequently, the long term cooling from ca. 5,600 cal. yr BP (most pronounced after 4,500 cal. yr BP) drove a gradual decrease in farming productivity causing a reduction in population size and contributed to a later resurgence of hunter-gatherer-fisher communities. Our study suggests a rapid human reaction within only two to three generations to the new opportunities created by the improved conditions for farming, emphasizing the remarkable plasticity humans can demonstrate in response to climate changes.

Methods

Sampling and preparation

Sediment core 303600 was collected from the central Baltic Sea in the Gotland Basin just east of Gotland Island at 56°55.02′ N and 19°19.98′ W (Fig. 1, Table 1) using a gravity corer during a cruise of the research vessel (R/V) “Professor Albrecht Penck” in July 2009. The collected core was cut in 1-m sections and split lengthwise. The core halves were line-scanned (photographic documentation) and sampled at a 1 cm resolution. The samples were freeze-dried, ground, and homogenized by mortar and pestle for analysis. The sediment chronology was developed using gravity cores 370540-7 and P435-2-1 taken during cruises with R/V “Aranda” in 2009 and R/V “Poseidon” (P435) in 2012, respectively, in the eastern Gotland Basin. For comparative studies a number of cores were taken during several cruises between 2005 and 2014 in almost all Baltic Sea basins (Fig. S1, Table 1).

Radiocarbon dating

The sediment cores were radiocarbon dated using picked benthic forminiferal tests (i.e. Elphidium spp.), which were present in high numbers in cores P435-2-1 (220 m water depth) and 370540-7 (242 m water depth) from the deepest part of the eastern Gotland Basin as a result of strong saline water inflows³¹. Individually picked benthic foraminifer samples (69 in total) were radiocarbon dated using accelerator mass spectrometry (AMS) to obtain a highly resolved age-depth model. Eight samples of 3-5 mg foraminiferal test were leached to remove any potential surface contamination, before the remaining 70% was decomposed to CO₂ under vacuum and converted to graphite for AMS ¹⁴C measurement at the Poznan Radiocarbon Laboratory in Poland. The majority of the samples was too small for routine AMS measurements on graphite and was dated with a new technique at ETH Zürich employing an AMS system equipped with a sample introduction system based on a gas ion source³². Foraminiferal tests (0.3-1 mg) were leached with 100 µl of 0.02 M HCl to remove 100 µg carbonate (12 µg carbon), before the remaining carbonate was completely dissolved in phosphoric acid in 3 ml septa sealed vials and the formed CO₂ was introduced into the ion source of the AMS for ¹⁴C measurement³³. The foraminiferal abundance and quality in the lower core interval of 370540-7 is much higher than in the equivalent depth interval of P435-2-1. Therefore, high resolution dating of this older interval concentrated on core 370540-7. Nevertheless, in order to secure a sound and robust high-resolution composite (P435-2-1 and 370540-7) age-depth model, overlapping depths intervals in P435-2-1 were also dated (Figs S2 and S3). The raw AMS ¹⁴C dates of benthic foraminiferal samples were calibrated using the Marine13 database in the Calib7.01 (ref.³⁴) and adjusted with a standard reservoir age of 400 years commonly used in western Baltic Sea studies (Table S1). Since benthic foraminifera only occurred after strong saline water inflows in the saline, oxygen-rich deep water body the use of this reservoir age is deemed adequate.

Total organic carbon (TOC) content and Loss on ignition (LOI)

The TOC content of sediments of core 303600 and P435-2-1 was determined by measuring the total carbon (TC) using the EA 1110 CHN analyzer (CE Instruments) and the total inorganic carbon (TIC) using the Multi EA- 2000 Elemental Analyzer (Analytic Jena); the TOC content was calculated as the difference between TC and TIC. The LOI was determined by ashing freeze-dried samples at 550 °C for 3 h and calculating the resulting mass difference in wt.%. For Baltic Sea sediments, the LOI provides a proper estimate of the TOC content of the sediments³⁵ (see also Fig. S2).

Correlation of cores

Benthic foraminifera are absent in shallow water sediment cores (like 303600) since incursion of saline waters did not reach this shallower site. Central Baltic Sea cores can be correlated to by comparing loss on ignition (LOI) records, which reflect the burial flux of organic carbon to the sea floor, which is related to primary productivity and bottom water oxygenation. Hence, the age-depth relationship developed on the central Baltic Sea cores from the deeper cores, P435-2-1 and 370540-7 was transferred onto shallower cores from the central/southern Baltic by correlating LOI and/ or organic carbon profiles (Figs S2 and S3) as performed earlier^9,36. The 370540-7 depth scale was transferred onto P435-2-1 equivalent depth scale by correlating the LOI profiles using Analyseries Software. The resulting relationship for age versus P435-2-1 equivalent depth is shown in Fig. S3. A polynomial fit was applied for the age model and this model was transferred onto 303600 using high-resolution LOI data in Analyseries Software. AMS ¹⁴C dating of sediment cores from the northern Baltic Sea basins were performed in order to obtain information on the approximate duration of the early Littorina Sea Phase hypoxic stage (Table S1).

TEX₈₆ analysis

Typically 1-3 g of freeze dried sediment were extracted using the Dionex^TM accelerated solvent extraction with dichloromethane/methanol (9:1; v/v) as solvent. The total lipid extract was dried over a Na₂SO₄ column and then separated using Al₂O₃ column chromatography. The polar fraction was obtained using dichloromethane/methanol (1:1. v/v) as the eluent. The polar fraction was concentrated under N₂ gas and then dissolved in hexane/isopropanol (99:1, v/v) before being filtered using a 0.4 µm PTFE filter as described previously³⁷. The sample was then analyzed by high-performance liquid chromatography coupled with atmospheric pressure chemical ionization-mass spectrometry (HPLC/APCI-MS) as described previously³⁸. A local TEX₈₆ calibration for the Baltic Sea was previously established⁹ using nine surface sediment samples from the Baltic Sea and comparing the TEX₈₆ derived SSTs with satellite sea surface temperature (SST) data. The highest correlation was found using the TEX₈₆ ^L index³⁹ and with the SST of the months July through October. The resultant Baltic Sea summer SST calibration equation is:

$${\rm{SST}}={\rm{34.03}}\times {{{\rm{TEX}}}_{{\rm{86}}}}^{{\rm{L}}}+36.73({{\rm{R}}}^{{\rm{2}}}={\rm{0.89}},\,{\rm{n}}={\rm{9}})$$

Model Simulation

To simulate the climatic impact on biogeochemistry of the Baltic Sea, we used a coupled three dimensional hydrodynamic and biogeochemical model. The hydrodynamic part of the model is based on the circulation model MOM (version 5.1)⁴⁰ and has been adapted to the Baltic Sea with an open boundary condition to the North Sea, and freshwater riverine input. The MOM model is complemented with a sea ice model to estimate ice cover thickness and extent. Furthermore, a parametric surface wave model⁴¹ is coupled to MOM with the main purpose of controlling the resuspension of sedimentary material. The biogeochemistry was integrated using the ERGOM biogeochemical model (www.ergom.net), which has been developed for the Baltic Sea, and describes the dynamics of nitrogen, oxygen and phosphorus including the inorganic nutrients nitrate, ammonium and phosphate, particulate organic matter consisting of phytoplankton (autotrophs), dead organic matter (detritus) and zooplankton (heterotrophs). In the ERGOM model, in-situ organic matter is produced from the inorganic nutrients by three functional groups of phytoplankton: large cells, small cells and others, and cyanobacteria. Organic material sinks and enters the model sediment as benthic nitrogen and phosphorus. All model processes requiring electron acceptors (e.g. respiration) reduce the oxygen concentration in the water. When oxygen is depleted, nitrogen is first used as an electron acceptor (denitrification) and then subsequently sulfate becomes the electron acceptor resulting in the formation of hydrogen sulfide.

During the Holocene the northern Baltic Sea lifts up and becomes shallower due to the glacio-isostatic rebound, affecting its bathymetry. To reconstruct the Baltic Sea bathymetry about 6,500 cal. yr BP, we used published data⁴² for relative sea level change. These data were combined with a contemporary bathymetry of the Baltic region⁴³. The northern basins are sometimes >100 m deeper, while the southern Baltic remains nearly unchanged (Fig. S5). Due to the tremendous computational load (2,000 years have been simulated), we applied a relative coarse model grid of 6 nautical miles (Fig. S5). Vertically the model is resolved into 80 layers. Close to the surface the vertical resolution starts with 0.5 m extending to the bottom to a maximum of 6 m.

Atmospheric forcing was reconstructed using the coastDat data set⁴⁴, a dynamical downscaling re-analysis for the years 1948–2014. For the different climate states (warm vs. cold), we modified the atmospheric data by a delta change approach successfully applied earlier⁹. This approach is typically applied when meteorological data are not available in a sufficient temporal and spatial resolution. In order to simulate the biogeochemical consequences of the change in summer SST that occurred in the reconstructed SST record at 6,000 cal. yr BP, we applied the following conditions. For the warm climate state we increased the 10 m air temperature of the coastDat data by 1 K. To mimic a colder climate we reduced the 10 m air temperature by 2 K and decreased the downwelling shortwave radiation. The effect is a mean SST difference of 2 K between cold and warm climate periods. The temperature difference is higher in summer and amounts up to ca. 3 K (Fig. S6). This temperature change reflects the reconstructed temperature difference for the transition at 6,000 cal. yr BP (Fig. 2).

Pre-industrial nutrient loads to the Baltic Sea due to riverine discharge and atmospheric deposition are derived from previous estimates⁴⁵. We used their estimates for AD 1850 and reduced them by additional 10% to account for a smaller population density at 6,500 cal. yr BP.

We performed the simulations with our model in 50-year time slices, which were eventually concatenated into the complete simulated time-series. The first time slice was initialized using data from a pre-industrial reconstruction of the Baltic Sea’s environment⁴⁰, which is characterized by low nutrient and high oxygen concentrations in the waters of the Baltic. The output of the preceding time slice was used as the starting condition for the next time slice in the simulation. Each time slice was forced by the modified meteorological data producing respective warm or cold climatic conditions. The modeled cold period (350 yr) mimics the cold climate before 6,000 cal. yr BP and the modeled warm phase (400 yr) the warm climate after 6,000 cal. yr BP (Fig. S6).

Overall population and cultural composition proxies

Based on the premise that periods of increased human population will leave behind more dateable archaeological artifacts^1,13, local population fluctuations in the Baltic region were reconstructed using a Summed Probability Distribution (SPD) of 1,960 terrestrial radiocarbon dates from 608 archeological sites (see refs^46,47 for responses to criticisms of this approach; for comparing radiocarbon and other proxies see ref.⁴⁸ ¹⁴C dates (Fig. 4A black line, top) that fall within the Central and Western Baltic kml polygon (brown polygon Fig. 1) were obtained primarily from the EUROEVOL radiocarbon database (http://discovery.ucl.ac.uk/1469811/). Some additional dates, which became available after publication of the EUROEVOL database were also added. The polygon was specifically constructed to test the hypothesis of a demographic response to climatic change across the broad region of Central and Western Baltic where farming was adopted c.6000 BP. By utilizing the largest possible sample size (site n = 608, sample n = 1960) we can expect local variability (e.g. from endogenous causes) to destructively interfere, and any patterns common to the entire region to be maximized. All these ¹⁴C dates have an anthropogenic association and are of terrestrial origin from archaeological sites. The SPD was constructed using methods previously described^3,13. At each site, similar dates (within 200 ¹⁴C years) were binned, calibrated using Intcal13, summed and normalized to unity, such that multiple similar dates at a site had the same overall weighting in the SPD as a single date. This is a conservative approach to dealing with differential dating at sites due to, for example, a larger research budget permitting more samples to be radiocarbon dated. All binned SPDs were then summed and the final SPD was normalized to unity. The SPD provides a proxy for relative population levels through time within the polygon, but the density values are relative and therefore the units are arbitrary.

We also constructed SPDs from four subsets of ¹⁴C dates (Fig. 4B) from the Central and Western Baltic polygon, utilizing a total of 1144 of the 1960 ¹⁴C dates (58%) which had been assigned to one of the following cultures by the original excavator: Ertebølle (or Lihult), TRB (Funnel Beaker), Corded Ware, Pitted Ware. Dates that were considered transitional between two cultures were not included, i.e. Ertebølle/TRB (n = 79) and TRB/Bell Beaker (n = 1). Also not included were 24 dates associated with cultures beyond the temporal interest of this study (for example Early Bronze Age). The remaining 712 excluded dates had not been associated with a specific culture. SPD analysis was performed as described above for all dates in the polygon, with the exception of the final step of normalization to unity. This ensures a fair comparison of the magnitudes of each cultural SPD. These four probability distributions were then summed, and the percentage contribution of each of the four was calculated, thus rescaling each SPD to a percentage between zero and 100.

Correlation tests between SST and Central and Western SPD

Correlation between SST and Central and Western Baltic SPD was calculated using two methods. In both cases interpolation was required to simplify and match the SPD time series to exactly the same 99 unequal discrete SST time points that overlap with the SPD, from 4,036 to 7,099 yrs BP inclusive. The first method uses Pearson’s product moment correlation coefficient, r = 0.731 (95% CI = 0.624 to 0.812), p-value < 2.2e-16, t = 10.56, df = 97. However, there is the possibility of some spurious shared information between the SST and SPD, which may exaggerate the p-value. For example, SST core time points were estimated using a ¹⁴C time depth model from dates that were calibrated via the same Intcal13 curve. Similarly both time-series are subject to the same atmospheric variability in ¹⁴C. Therefore, a second method using a Monte Carlo simulation approach as previously employed³ generated 50,000 simulated radiocarbon datasets (using a uniform distribution between 4,000 and 8,000 cal. yr BP). An SPD was generated for each, which were then subjected to exactly the same correlation test with SST. This ensured each simulated SPD was processed using the same calibration and interpolation as the observed radiocarbon dataset and resulted in 50,000 Pearson’s product moment correlation coefficients ranging from −0.735 to + 0.721 (SD = 0.228), which were compared to the observed correlation coefficient. Not a single simulated correlation coefficient was equal or greater (1-tailed) and only one simulation was greater in magnitude but negative (2-tailed), providing a more conservative p-value estimate of 0.00002.

Both the SPD and SST time series were interpolated to annual points between 4,036 and 7,099 yrs BP, allowing a fine resolution cross-correlation analysis using Pearson’s product moment correlation coefficient in the lag range of ± 15% (460 years). A peak correlation was obtained when the SPD lagged SST by 61 years (r = 0.862 compared to r = 0.839 at lag = 0). Note only the relative values are informative as the magnitude of these values is inflated compared to the correlation test since both time series were heavily interpolated and trimmed to ensure a constant window length.

Analysis of faunal data

Zooarchaeological data from published Mesolithic/Neolithic archaeological sites in the Baltic region were recorded in a relational SQL database. The total Number of Identified Specimens (NISP) for all categories of anthropogenic faunal deposits were recorded, including wild and domestic mammals, birds, fish, molluscs, reptiles, amphibians and crustaceans. Intrusive faunal elements, where identified, were excluded from the database. Taxa were entered at the level of identification published by the original analysts, whether to species or genera, or at the more general level of family (or ‘type’), or according to body size (e.g. large/small mammal). Each taxonomic entry was then associated with a cultural phase e.g. Ertebølle, Pitted Ware etc, and assigned a unique phase identification code. Any published radiocarbon dates derived from the same archaeological contexts of that phase were recorded with the same phase identification code to ensure accurate site phase chronology (see SI). We combined multiple sources of chronological evidence in a Bayesian framework using OxCal, including radiocarbon dates at the target site-phase, date ranges estimated by the excavator using typological evidence, and radiocarbon dates from neighbouring site-phases from the same culture or period (see SI). We used a two-tailed permutation test to fully integrate the uncertainty in the date estimate of each site-phase, and sampling uncertainty in NISP counts (see SI).

Data availability

Upon publication the datasets generated in the current study will be available in the Pangea repository (www.pangaea.de).

References

Bellwood, P. First Farmers: The Origins of Agricultural Societies. (Blackwell, Oxford, 2005).
Bonsall, C., Macklin, M. G., Anderson, D. E. & Payton, R. W. Climate change and the adoption of agriculture in North-WestEurope. Eur. J. Archaeol. 5, 9–23 (2002).
Article Google Scholar
Shennan, S. et al. Regional population collapse followed initial agriculture booms in mid-Holocene Europe. Nat. Commun. 4, 2486 (2013).
Article PubMed PubMed Central CAS Google Scholar
Bogucki, P., Nalepka, D., Grygiel, R. & Nowaczyk, B. Multiproxy environmental archaeology of Neolithic settlements at Osłonki, Poland, 5500–4000 BC. Environ. Archaeol. 17, 45–65 (2012).
Article Google Scholar
Price, T. D. in Europe’s First Farmers. (ed. Price, T.D.) 260–300 (Cambridge University Press, Cambridge, 2000).
Krause-Kyora, B. et al. Use of domesticated pigs by Mesolithic hunter-gatherers in northwestern Europe. Nat. Commun. 4, 2348 (2013).
Article PubMed PubMed Central Google Scholar
Banks, W. E., Antunes, N., Rigaud, S. & d’Errico, F. Ecological constraints on the first prehistoric farmers in Europe. J. Archaeol. Sci. 40, 2746–2753 (2013).
Article Google Scholar
Schouten, S., Hopmans, E. C., Schefuß, E. & Sinninghe Damsté, J. S. Distributional variations in marine crenarchaeotal membrane lipids: A new tool for reconstructing ancient sea water temperatures? Earth Planet. Sci. Lett. 204, 265–274 (2002).
Article ADS CAS Google Scholar
Kabel, K. et al. Impact of climate change on the Baltic Sea ecosystem over the past 1,000 years. Nat. Clim. Change 2, 871–874 (2012).
Article ADS Google Scholar
Diaz, R. J. & Rosenberg, R. Spreading dead zones and consequences for marine ecosystems. Science 321, 926–929 (2008).
Article ADS CAS PubMed Google Scholar
Karlson, K., Rosenberg, R. & Bonsdorff, E. Temporal and spatial large-scale effects of eutrophication and oxygen deficiency on benthic fauna in Scandinavian and Baltic waters: A review. Oceanogr. Mar. Biol. 40, 427–489 (2002).
Google Scholar
Vahtera, E. et al. Internal ecosystem feedbacks enhance nitrogen-fixing cyanobacteria blooms and complicate management in the Baltic Sea. AMBIO J. Hum. Environ. 36, 186–194 (2007).
Article CAS Google Scholar
Wasmund, N. Occurrence of cyanobacterial blooms in the Baltic Sea in relation to environmental conditions. Int. Rev. Gesamten Hydrobiol. Hydrogr. 82, 169–184 (1997).
Article Google Scholar
Kononen, K. Dynamics of the toxic cyanobacterial blooms in the Baltic Sea. Finn. Mar. Res. 261, 3–36 (1992).
Google Scholar
Laima, M. J., Matthiesen, H., Christiansen, C., Lund-Hansen, L. C. & Emeis, K.-C. Dynamics of P, Fe and Mn along a depth gradient in the SW Baltic Sea. Boreal Environ. Res. 6, 317–333 (2001).
CAS Google Scholar
Neumann, T. & Schernewski, G. Eutrophication in the Baltic Sea and shifts in nitrogen fixation analyzed with a 3D ecosystem model. J. Mar. Syst. 74, 592–602 (2008).
Article ADS Google Scholar
Struck, U., Pollehne, F., Bauerfeind, E. & von Bodungen, B. Sources of nitrogen for the vertical particle flux in the Gotland Sea (Baltic Proper)—results from sediment trap studies. J. Mar. Syst. 45, 91–101 (2004).
Article ADS Google Scholar
Davis, B. A., Brewer, S., Stevenson, A. C. & Guiot, J. The temperature of Europe during the Holocene reconstructed from pollen data. Quat. Sci. Rev. 22, 1701–1716 (2003).
Article ADS Google Scholar
Tallavaara, M. & Seppä, H. Did the mid-Holocene environmental changes cause the boom and bust of hunter-gatherer population size in eastern Fennoscandia? The Holocene 22, 215–225 (2012).
Article ADS Google Scholar
Timpson, A. et al. Reconstructing regional population fluctuations in the European Neolithic using radiocarbon dates: a new case-study using an improved method. J. Archaeol. Sci. 52, 549–557 (2014).
Article Google Scholar
Surovell, T. A., Finley, J. B., Smith, G. M., Brantingham, P. J. & Kelly, R. Correcting temporal frequency distributions for taphonomic bias. J. Archaeol. Sci. 36, 1715–1724 (2009).
Article Google Scholar
Discamps, E. Ungulate biomass fluctuations endured by Middle and Early Upper Paleolithic societies (SW France, MIS 5-3): The contributions of modern analogs and cave hyena paleodemography. Quatern. Intern. 337, 64–79 (2014).
Article ADS Google Scholar
Bocquet-Appel, J. Paleoanthropological traces of a Neolithic demographic transition. Curr. Anthropol. 43, 637–650 (2002).
Article Google Scholar
Sørensen, L. & Karg, S. The expansion of agrarian societies towards the north – new evidence for agriculture during the Mesolithic/Neolithic transition in Southern Scandinavia. J. Archaeol. Sci. 51, 98–114 (2014).
Article Google Scholar
Carter, T. R., Saarikko, R. A. & Niemi, K. J. Assessing the risks and uncertainties of regional crop potential under a changing climate in Finland. Agric. Food Sci. 5, 329–350 (1996).
Google Scholar
Skoglund, P. et al. Genomic diversity and admixture differs for Stone-Age Scandinavian foragers and farmers. Science 344, 747–750 (2014).
Article ADS CAS PubMed Google Scholar
Mittnik, A. et al. The genetic history of northernEurope. bioRxiv113241, https://doi.org/10.1101/113241 (2017).
Rowley-Conwy, P. Westward Ho!: The spread of agriculture from Central Europe to the Atlantic. Curr. Anthropol. 52, S431–S451 (2011).
Article Google Scholar
Feeser, I., Dorfler, W., Averdieck, F.D. & Wiethold, J. New insight into regional and local land-use and vegetation patterns in eastern Schleswig-Holstein during the Neolithic. In: Siedlung, Grabenwerk, Grossteingrab. Studien zu Gesellschaft, Wirtschaft und Umwelt der Trichterbechergruppen im nordlichen Mitteleuropa. Fruhe Monumentalitat und Soziale Differenzierung (eds. M. Hinz and J. Muller) pp. 159–190 (2012).
Breitburg, D. L., Pihl, L. & Kolesar, S. E. Effects of low dissolved oxygen on the behavior, ecology and harvest of fishes: a comparison of the Chesapeake Bay and Baltic-Kattegat systems. In Coastal Hypoxia: Consequences for Living Resources and Ecosystems (eds. N. N. Rabalais and R. E. Turner), American Geophysical Union, Washington, D.C., pp. 241-267 (2001).
Lutze, G. F. Zur Foraminiferen-Fauna der Ostsee. Meyniana 14, 75–142 (1965).
Google Scholar
Wacker, L., Lippold, J., Molnár, M. & Schulz, H. Towards radiocarbon dating of single foraminifera with a gas ion source. Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. At. 294, 307–310 (2013).
Article ADS CAS Google Scholar
Bard, E. et al. AixMICADAS, the accelerator mass spectrometer dedicated to 1^4C recently installed in Aix-en-Provence, France. Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. At. 361, 80–86 (2015).
Article ADS CAS Google Scholar
Reimer, P. J. et al. IntCal13 and Marine13 radiocarbon age calibration curves 0-50,000 years cal BP. Radiocarbon 55, 1869–1887 (2013).
Article CAS Google Scholar
Leipe, T. et al. Particulate organic carbon (POC) in surface sediments of the Baltic Sea. Geo-Mar. Lett. 31, 175–188 (2011).
Article ADS CAS Google Scholar
Lougheed, B. C. et al. Using an independent geochronology based on palaeomagnetic secular variation (PSV) and atmospheric Pb deposition to date Baltic Sea sediments and infer ¹⁴C reservoir age. Quat. Sci. Rev. 42, 43–58 (2012).
Article ADS Google Scholar
Hopmans, E. C. et al. A novel proxy for terrestrial organic matter in sediments based on branched and isoprenoid tetraether lipids. Earth Planet. Sci. Lett. 224, 107–116 (2004).
Article ADS CAS Google Scholar
Schouten, S., Huguet, C., Hopmans, E. C., Kienhuis, M. V. & Sinninghe Damsté, J. S. Analytical methodology for TEX₈₆ paleothermometry by high-performance liquid chromatography/atmospheric pressure chemical ionization-mass spectrometry. Anal. Chem. 79, 2940–2944 (2007).
Article CAS PubMed Google Scholar
Kim, J.-H. et al. New indices and calibrations derived from the distribution of crenarchaeal isoprenoid tetraether lipids: Implications for past sea surface temperature reconstructions. Geochim. Cosmochim. Acta 74, 4639–4654 (2010).
Article ADS CAS Google Scholar
Griffies, S. M. Fundamentals of Ocean Climate Models. (Princeton University Press, Princeton, 2004). pp. 494.
Schwab, D. J., Bennett, J. R., Liu, P. C. & Donelan, M. A. Application of a simple numerical wave prediction model to Lake Erie. J. Geophys. Res. Oceans 1978–2012 89, 3586–3592 (1984).
Rosentau, A., Harff, J., Oja, T. & Meyer, M. Postglacial rebound and relative sea level changes in the Baltic Sea since the Litorina transgression. Baltica 25, 113–120 (2012).
Article Google Scholar
Seifert, T., B. Kayser and F. Tauber. “Bathymetry data of the Baltic Sea.” http://www.io-warnemuende.de/topography-of-the-baltic-sea.html.
Weisse, R. et al. Regional meteorological-marine reanalyses andclimate change projections: Results for Northern Europe and potential for coastal and offshore applications. Bull. Am. Meteorol. Soc. 90, 849–860 (2009).
Article ADS Google Scholar
Gustafsson, B. G. et al. Reconstructing the development of Baltic Sea eutrophication 1850–2006. Ambio 41, 534–548 (2012).
Article CAS PubMed PubMed Central Google Scholar
Timpson, A., Manning, K. & Shennan, S. Inferential mistakes in population proxies: A response to Torfing’s “Neolithic population and summed probability distribution of ¹⁴C-dates”. J. Archaeol. Sci. 63, 199–202 (2015).
Article Google Scholar
Tallavaara, M., Pesonen, P., Oinonen, M. & Heikki Seppä, H. The mere possibility of biases does not invalidate archaeological population proxies – response to Teemu Mökkönen. Fennosc. Archaeol. 31, 135–140 (2014).
Google Scholar
Lechterbeck, J. K. et al. Is Neolithic land use correlated with demography? An evaluation of pollen-derived land cover and radiocarbon-inferred demographic change from CentralEurope. Holocenev 24, 1297–1307 (2014).
Article ADS Google Scholar
Reinholdsson, M., Snowball, I., Zillén, L., Lenz, C. & Conley, D. J. Magnetic enhancement of Baltic Sea sapropels by greigite magnetofossils. Earth Planet. Sci. Lett. 366, 137–150 (2013).
Article ADS CAS Google Scholar
Häusler, K. et al. Mid‐to late Holocene environmental separation of the northern and central Baltic Sea basins in response to differential land uplift. Boreas 46, 111–128 (2017).
Article Google Scholar

Download references

Author information

Authors and Affiliations

NIOZ Royal Netherlands Institute for Sea Research, Department of Marine Microbiology and Biogeochemistry, and Utrecht University, 1790 AB Den Burg, Texel, The Netherlands
L. Warden, M. Sollai & J. S. Sinninghe Damsté
Leibniz Institute for Baltic Sea Research, Departments of Marine Geology and Physical Oceanography, Warnemünde, Germany
M. Moros, T. Neumann, K. Perner, K. Häusler, T. Leipe & H. Arz
Institute of Archaeology, University College London, 31-34 Gordon Square, London, WC1H 0PY, UK
S. Shennan
Department of Genetics, Evolution and Environment, University College London, Gower Street, London, WC1E 6BT, UK
A. Timpson
Department of Geography, King’s College London, London, WC2R 2LS, UK
K. Manning
Laboratory of Ion Beam Physics, ETH Zürich, Otto-Stern-Weg 5, 8093, Zürich, Switzerland
L. Wacker
Geological Survey of Sweden, Department of Marine Geology, Box 670, 751 28, Uppsala, Sweden
L. Zillén
Geological Survey of Finland, P.O. Box 96, 02151, Helsinki, Finland
A. Kotilainen
Bjerknes Centre for Climate Research/Department of Earth Science, University of Bergen Allégaten 70, 5007, Bergen, Norway
E. Jansen
Insitute of Geosciences, Kiel University, Ludewig-Meyn Strasse 10, 24118, Kiel, Germany
R. R. Schneider
Institute of Baltic Sea Fisheries, Thünen Institute, Alter Hafen Süd 2, 18069, Rostock, Germany
R. Oeberst
Utrecht University, Faculty of Geosciences, P.O. Box 80.021, 3508 TA, Utrecht, The Netherlands
J. S. Sinninghe Damsté

Authors

L. Warden
View author publications
You can also search for this author in PubMed Google Scholar
M. Moros
View author publications
You can also search for this author in PubMed Google Scholar
T. Neumann
View author publications
You can also search for this author in PubMed Google Scholar
S. Shennan
View author publications
You can also search for this author in PubMed Google Scholar
A. Timpson
View author publications
You can also search for this author in PubMed Google Scholar
K. Manning
View author publications
You can also search for this author in PubMed Google Scholar
M. Sollai
View author publications
You can also search for this author in PubMed Google Scholar
L. Wacker
View author publications
You can also search for this author in PubMed Google Scholar
K. Perner
View author publications
You can also search for this author in PubMed Google Scholar
K. Häusler
View author publications
You can also search for this author in PubMed Google Scholar
T. Leipe
View author publications
You can also search for this author in PubMed Google Scholar
L. Zillén
View author publications
You can also search for this author in PubMed Google Scholar
A. Kotilainen
View author publications
You can also search for this author in PubMed Google Scholar
E. Jansen
View author publications
You can also search for this author in PubMed Google Scholar
R. R. Schneider
View author publications
You can also search for this author in PubMed Google Scholar
R. Oeberst
View author publications
You can also search for this author in PubMed Google Scholar
H. Arz
View author publications
You can also search for this author in PubMed Google Scholar
J. S. Sinninghe Damsté
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

M.M. and J.S.S.D. designed the research. M.M. and H.A. conducted sampling campaigns. T.L. and K.H. obtained TOC/LOI data. L.W., M.S. and J.S.S.D. obtained the biomarker data. T.N. performed the numerical modelling and paleobathymetry reconstruction. S.S., K.M. and A.T. analysed the archaeological data and generated the population reconstruction curves. M.M., K.P., L.Z., A.K., R.R.S. and E.J. constructed the chronology. L.Wacker performed radiocarbon dating on small foraminifera samples. L.W., M.M., S.S., A.T., T.N. and J.S.S.D. interpreted the data. L.W., M.M., and J.S.S.D. wrote the paper with contributions from all other authors.

Corresponding authors

Correspondence to M. Moros or J. S. Sinninghe Damsté.

Ethics declarations

Competing Interests

The authors declare that they have no competing interests.

Additional information

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

Electronic supplementary material

Supplementary PDF File

Rights and permissions

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 license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Warden, L., Moros, M., Neumann, T. et al. Climate induced human demographic and cultural change in northern Europe during the mid-Holocene. Sci Rep 7, 15251 (2017). https://doi.org/10.1038/s41598-017-14353-5

Download citation

Received: 15 May 2017
Accepted: 09 October 2017
Published: 10 November 2017
DOI: https://doi.org/10.1038/s41598-017-14353-5

This article is cited by

Long-Term Urban and Population Trends in the Southern Mesopotamian Floodplains
- Nicolò Marchetti
- Eugenio Bortolini
- Federico Zaina
Journal of Archaeological Research (2024)
100 ancient genomes show repeated population turnovers in Neolithic Denmark
- Morten E. Allentoft
- Martin Sikora
- Eske Willerslev
Nature (2024)
A new perspective on Tegelbarg: character and chronology of a Late Neolithic shell midden in the western Baltic
- J.P. Kleijne
- D.J. Huisman
Archaeological and Anthropological Sciences (2023)
Approximate Bayesian Computation of radiocarbon and paleoenvironmental record shows population resilience on Rapa Nui (Easter Island)
- Robert J. DiNapoli
- Enrico R. Crema
- Terry L. Hunt
Nature Communications (2021)
Genome-wide study of a Neolithic Wartberg grave community reveals distinct HLA variation and hunter-gatherer ancestry
- Alexander Immel
- Federica Pierini
- Ben Krause-Kyora
Communications Biology (2021)

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.