Reconstruction of last millennia Sea Surface Temperature (SST) evolution is challenging due to the difficulty retrieving good resolution marine records and to the several uncertainties in the available proxy tools. In this regard, the Roman Period (1 CE to 500 CE) was particularly relevant in the socio-cultural development of the Mediterranean region while its climatic characteristics remain uncertain. Here we present a new SST reconstruction from the Sicily Channel based in Mg/Ca ratios measured on the planktonic foraminifer Globigerinoides ruber. This new record is framed in the context of other previously published Mediterranean SST records from the Alboran Sea, Minorca Basin and Aegean Sea and also compared to a north Hemisphere temperature reconstruction. The most solid image that emerges of this trans-Mediterranean comparison is the persistent regional occurrence of a distinct warm phase during the Roman Period. This record comparison consistently shows the Roman as the warmest period of the last 2 kyr, about 2 °C warmer than average values for the late centuries for the Sicily and Western Mediterranean regions. After the Roman Period a general cooling trend developed in the region with several minor oscillations. We hypothesis the potential link between this Roman Climatic Optimum and the expansion and subsequent decline of the Roman Empire.
The semi-enclosed configuration of the Mediterranean Sea makes this an extremely vulnerable region to modern and also to past climate changes1.The strategic transitional zone that occupies the Mediterranean, between North Africa and European climates, from the arid zone of the subtropical high to the humid northwesterly air flows, provides to the Mediterranean a particular interest to unravel climate tele-connections during times of climate variability2,3,4. Several studies carried out in different marine sites, have focused the attention to the short-term climate variability over the last millennia but so far, any general reconstruction of the regional temperature evolution has not been attained5,6,7,8,9,10,11,12. Compared to other regions of the world, the Mediterranean is characterised by a wealth of archaeological studies and historic documents, as well as a paleoclimatic data that makes it a perfect case study to investigate the potential influence of climate on civilisations13. In fact, this time period is particularly challenging since coincided with important cultural changes (human civilizations) that developed around the Mediterranean area14. The study of the fossil archives remains the only valid tool to reconstruct past environmental and climatic changes during those times15. However, its application in the marine realm is compromised by the difficulty to obtain marine records with resolution enough, and clear proxy signals that can be reproduced among different records. Climate variability of this period is often close or within the proxy errors and uncertainties in its interpretation in terms of seasonality and/or local oceanographic processes, make difficult always its reading in terms of regional climate evolution11,12,16,17,18. Nevertheless, this is a critical information to identify past interactions between climate changes and evolution of human societies and their adaptive strategies19,20,21,22,23,24,25. In addition, the last report of the Intergovernmental Panel on Climate Change (IPCC 2018) underlines the requirement to assess climate feedbacks during past episodes of moderately warmer (1.5–2 °C) conditions26,27.
In this framework, we present a new generated SST record reconstruction from the central part of the Mediterranean Sea based on Mg/Ca ratios measured in the planktonic foraminifera Globigerinoides ruber covering the last 5 kyr BP. This new record is compared with other previously published SST records from the Alboran Sea, Menorca Basin and Aegean Sea based on different geochemical proxies and with a reconstruction of the NAO index28. We argue that this trans-Mediterranean compilation provides the basis to discuss main SST features during the last two millennia focussing the attention into the Roman period since it underscores as the warmest period with interesting implications on civilization development in the Mediterranean region.
The Mediterranean Sea is an anti-estuarine semi-enclosed sea that can be subdivided into two sub-basins, the western and eastern Mediterranean separated by the Sicily Strait sill29. Low salinity surface waters, called Modified Atlantic Water (MAW), enter in the Mediterranean Sea from the Strait of Gibraltar and occupy the first 0–100 m of the water column which overlies the outflowing saltier Levantine Intermediate Water (LIW) formed in the Levantine basin29. The MAW entrained by an intense jet coupled to a wave-like front30,31, cross the Alboran Sea and exits in its eastern end along the African coast32.
MAW flows along the Algerian coast as the coastal Algerian Current and separates into two branches at the entrance of the Sicily Strait33 (Fig. 1). The Strait of Sicily represents a physical barrier (about 500 m deep) of the eastern Mediterranean and implies considerable control over the biogeochemicals processes occurring within the eastern basin34. Most of the MAW passes through the Sicilian Channel35 dividing into two streams while the rest flows into the Tyrrhenian Sea36,37. The northern branch, called the Atlantic Ionian Stream (AIS), forms the MAW transport into the eastern Mediterranean off the southern coast of Sicily and baths the studied location38.
In detail, measurements of the main physical seawater parameters on the studied core location (SW104-ND11) in July 2014 (Fig. 1) show that surface temperatures (0–15 m) range from 23.0 °C to 22.0 °C. The base of the thermocline (17.5 °C) is located at a depth of 22 m and the salinity values are between 37.36 and 38.88 psu.
Results and discussion
The SST estimates from the Mg/CaG.ruber ratio range from 16.4 °C ± 1.5 °C to 22.7 °C ± 1.5 °C with a mean value of 19.5 °C ± 1.5 °C (Fig. 2); to document the main trends of variability [95% Confidence Interval (CI)] in SST reconstruction we adopted a Monte Carlo approach that use a non-parametric regression (LOESS function, see Material and Methods paragraph). The Mg/CaG.ruber SST reconstruction from core SW104-ND11 shows a progressive warming trend of 6.3 °C ± 2.0 °C from 3300 BCE (base of the sequence) to 330 CE, middle Roman Period when SST maxima are attained (Fig. 2). This long-term warming trend is punctuated by several short-term oscillations of different amplitude and duration (Fig. 2). From the Roman Period to 1700 CE, the SST shows a turnover into a cooling trend of 4.5 °C ± 2.1 °C (Fig. 2). The SST record ends from 1700 to 2014 CE with a short warming trend (Fig. 2).
This record is compared with other previously published SST reconstructions from the Mediterranean Sea (Fig. 2): an alkenone-SST record from the Alboran Sea39, a Mg/CaG.bulloides SST-stack record that integrates five Mg/Ca-SST records from the North Balearic Islands16 and a composite of two alkenone-SST records from the Aegean Sea40,41. This comparison is complemented with a record of north Hemisphere temperature anomalies42 and a reconstructed NAO record28. Thus, the compared SST records involve different proxies and calibrations with their own uncertainties. In order to facilitate their comparison, and eliminating potential biases in the absolute SST reconstruction associated with the applied calibrations or methods, we decided to compare the SST records transferring them into SST anomaly records in relation to a common reference period. Considering that the time interval covered by the considered marine records was different, the temperature anomalies (°C) were calculated respect to the only period shared by all the records, from 750 BCE to 1250 CE (Fig. 3).
Although some other Mediterranean SST records exist for the last millennia, previous compilation efforts have shown significant discrepancies in the main trends among different regions11,43. In this exercise, we have focussed on those records that are not directly influenced by large river outflows, human impact and local oceanographic circulation that could potentially alter the wide regional signal (See the supplementary material for detailed arguments on the criteria for record selection, paragraph S2). Different proxies reflect different seasons and depth habitats and thus SST signals are not expected to be identical among the chosen records. Nevertheless, this comparison exercise focuses on the main patterns and aims to evaluate their regional perseverance. In this regard, the most perseverant feature among the compared records is the well developed SST maximum recorded during the Roman Period (1–500 CE) (Fig. 3). The discussion on the SST evolution previous and posterior to this period is mostly focussed in the new Sicily SST record here presented.
The comparison of the Mediterranean SST records evidences very distinct regional patterns before the onset of the Greek period, with an overall cooling trend in the Aegean Sea, a warming trend in the Sicily Strait and more stable conditions in the Alboran Sea (Fig. 2). In the Sicily Channel, between the base of the record and the beginning of the Bronze Age period, the Mg/CaG.ruber SST signal shows two warm events, at ca. 2913 BCE and ca. 2040 BCE, respectively (Fig. 2). The first warm event, documented in the study record of Sicily Channel, although it does not fall within the 95% CI (Fig. 2), has been tentatively associated to the Copper Age warm phase (ca. 2913 BCE). This event chronologically corresponds to a warming in the Aegean Sea SST record (Fig. 2) and agrees to the beginning of a gradual aridification process described in northern Egypt44.
The second one (Early Bronze Age, ca. 2040 BCE, Fig. 2) is associated with a further aridity phase, as documented by the strong decrease in arboreal pollen in the central Mediterranean45,46. This latter event, chronologically corresponds to the fall of the ancient Egyptian Reign44 and to the end of the Mesopotamian civilization, associated with strong famines related to a dry climate and a strong aridity47. Later, Mg/Ca G.ruber SST warmed between 1800-1100 BCE (Fig. 2) indicating relative warm conditions during Late Bronze Age (ca. 1100 BCE) documented in the whole Mediterranean basin.
Significant cultural changes are documented during this period, in fact, historical records indicate the collapsed several civilizations48,49. Most of the Greek Bronze Age Palatial centres were destroyed and/or abandoned50. The Palatial centers were hit hard by the increase in aridity and the collapse of agriculture production48 making impossible for the population to sustain themselves50,51. The transition from Bronze Age to Iron Age chronologically approximate the short-term cooling event associated to the Homeric (ca. 800 BCE) grand solar minimum (Fig. 2). During this period, associated to negative NAO values, the climate condition were favourable to agriculture expansion in the eastern Mediterranean44. Compared to the subsequent Roman period, the Mediterranean was characterized by a colder phase from ca. 500 BCE to 200 BCE (Fig. 2), predating the rise of the Roman Empire. From a chronological point of view, this interval corresponds to the beginning of the so-called “sub-Atlantic phase”47,52,53, characterized by a cool climate and rainy winters which was propitious for the expansion of the Greek, Etruscan and Roman civilizations54. During this period, also global glacier advances are documented55 and a negative NAO phase is recorded (Fig. 2). The cool and humid climate of the sub-Atlantic phase lasted until ca. 100 BCE and covered the entire period of the Monarchy in Rome47. In addition, this period is characterised by a short-term cooling event associated with the Greek (ca. 350 BCE) solar minimum (Fig. 2).
However, at ca. 400 BCE, cultural changes were synchronised across the Mediterranean region14. The Greek and Phoenician colonies expanded and Rome and Carthage began their epic rise14, a situation coincident with the establishment of more homogeneous temperature conditions across the Mediterranean regions.
All the considered SST reconstructions, from the Alboran basin to the Aegean Sea, show homogeneous climate conditions at regional scale with the occurrence of a distinct warming phase from 1 CE to 500 CE (Fig. 3), coincident with the Roman Period and covering the whole Roman Empire archaeological period. This pronounced warming during the Roman Period is almost consistent with other marine records from Atlantic Ocean56,57,58 and with the continental anomaly reconstruction from Europe from PAGES 2 K, 201359. This climate phase corresponds to the so-called “Roman Climatic Optimum” characterized by prosperity and expansion of the empire.
The described warm climate conditions of the “Roman Climatic Optimum” appear independently of the considered SST proxies. In one side, Mg/CaG. ruber SST are expected to reflect temperatures from late summer to autumn60 while of Mg/CaG. bulloides SST in the Menorca area has been interpreted to reflect spring SST conditions16. On the other side, alkenone-SST signal is interpreted as mean annual averages12,16,39,61. The good agreement between the different SST proxies and Mediterranean regions reflecting warm conditions during the whole Roman Period suggests warm climate conditions during the whole annual cycle and at regional level. Nevertheless, some differences in the intensity of the warming arise between the different proxy records. This warming in SST reconstructions range from 19.6 ± 1.5 °C to 22.7 ± 1.5 °C in the Mg/CaG. ruber SST of the Sicily Channel (Fig. 2), from 16.6 to 18.5 °C in the Mg/CaG. bulloides SST from the Menorca basin record16 and from 14.4 to 16.1 °C in the Uk37 SST anomaly from both the Alboran and Aegean seas39,40. The overall larger Holocene variability in Mg/Ca-SST records in relation to alkenone-SST records have been previously recognized in the Alboran Sea and discussed to reflect the rather averaged and smoothed annual SST signal of the alkenone measurements61. Thus, our inter-basin comparison suggest that the warming was either more perseverant during the summer months or it was stronger in the central part of the Mediterranean in contrast to the eastern and western basins. The warmest part (ca. 250–420 CE) of the Mediterranean Roman Period is almost coincident with the warm phase in the north Hemisphere temperature anomalies42 (Fig. 2), suggesting a teleconnection between marine and continental system. Indeed, the “Roman Climatic Optimum” reveals itself as a phase of warm stable temperatures across much of the Mediterranean heartland of empire and covers the whole phase of origin, expansion and decline of the Roman Empire (Roman Climate Optimum, Roman Transitional Period and Late Antique Little Ice Age, according to Harper, 2019)62. Nevertheless, the expansion of the Roman Empire did not occur synchronously throughout the Mediterranean area because of the political situation, geographical pattern and, probably, also the specific climatic conditions in the various areas. However, starting from ca. 50 CE, the Roman Empire began its ascent towards north of Alps and the conquest of Gaul, probably assisted by the favourable climatic conditions that allowed this great feat. But, a shift from negative to positive NAO index values (Fig. 2) during the Roman period suggests a transfer of North Atlantic precipitation towards the central and northern Europe. In fact, positive NAO could indicate drier conditions by the second half of the Roman Period penalising southern Europe and North Africa economic development (i.e. agriculture)63,64. This feature agrees with the new finding by recent studies on vegetation patterns46,65 about the opening of vegetation in central/eastern Mediterranean areas since the ca. 100 CE. Thereafter the Rome expansion in central Europe and Great Britain was inexorable.
Therefore, it is possible to hypothesize that the auspicious climatic conditions of this period, the so-called “Roman Climatic Optimum”, may have contributed to the expansion of the Roman Empire in the various areas of the Mediterranean Sea probably linked to new needs and opportunities. By the end of this period, all the Mediterranean SST records show a generally consistent SST cooling trend (Fig. 3). During this period, both the Aegean Sea and Sicily Channel were laid siege by pirates and vandals. In fact, at the end of the IV century, the roman world was experiencing the beginning of the end that will culminate just few centuries later66. Sicily was a major supplier of grain throughout the Empire, and held the greatest part of the agrarian economy67. The island, deprived of military defences, was besieged by a band of pirates Illyrian (438 CE) which devastated a large stretch. In the same time, the Vandals of Genserico (440–477 CE) invaded the Eastern part of the Mediterranean Sea67.
Post roman period
The post-Roman Period was characterised by a progressive cooling trend of 4.5 °C ± 2.1 °C that ended by the Little Ice Age (LIA) interval (Fig. 2) and was punctuated by two short-term events well-known as the Late Antique Little Ice Age (LALIA) and the Medieval Warm Period (MWP). The LALIA cooling event23, at ca. 650–700 CE (Fig. 2), has been recently documented in marine records of the western Mediterranean basin9,10 and corresponds to a period characterised by a decrease in alboreal pollen index in the south-central Mediterranean, suggesting the occurrence of cold and dry climate condition45,46. The LALIA event had an additional environmental stress led by the establishment of the Justinian plague while the eastern Roman Empire transformed, the Sasanian Empire collapsed, and population movements occurred out of the Asian steppe and Arabian Peninsula, which spread Slavic-speaking peoples and political upheavals in China23. Instead, the MWP, dated at ca. 1300 CE (Fig. 2) and characterised by overall warmer SST than the LALIA event, is considered by several authors7,9,22,68,69,70 as a relatively stable and warm period. Recently, Margaritelli et al.10 documented the occurrence of this warm event at Mediterranean regional scale.
The cooling of 2.8 °C ± 2 °C associated to the LIA event occurring between ca. 1300 CE and ca. 1700 CE affects the whole Mediterranean area (Fig. 2). This event is also well documented in several marine records of the Mediterranean5,6,7,8,9,10,16,71,72,73. After 1700 CE, the SST from Sicily Channel, tend to warm (Fig. 2) although interrupted by a short cooling at 1980 CE74. This warming phase fits with the Aegean Sea SST record and with the north Hemisphere temperature anomalies5,40,75,76. This feature is probably associated to the onset of the Industrial Period/Modern Warm Period.
The new generated Mg/CaG. ruber SST record from Sicily channel allows to identify a series of climate events that can be associated to different remarkable socio-cultural developments of ancient Mediterranean civilizations for the last five millennia. Warm events are associated to historical periods such as the Cooper age, the Early and Late Bonze age while cold events are associated to the Homeric and Greek periods. The comparison of this new SST record with previous published Mediterranean SST records from Alboran Sea, Minorca Basin and Aegean Sea, highlight the overall perseverance of warm conditions during the Roman Period (1 CE to 500 CE). These warm conditions were particularly intense in the Sicily record, that reflect summer months, and correspond to the so called “Roman Climatic Optimum”. During this period, in fact, developed the greatest ancient civilization of all time, the Roman one. We hypothesise the relevance that these climate conditions may had in the expansion of the Roman Empire and its collapse with the general development of colder conditions. A cooling trend dominated after the Roman period reaching minimum values, in the whole Mediterranean, by the end of the LIA, minor oscillations punctuated this cold trend and were also associated with socio-cultural changes in central Mediterranean region.
Materials and methods
Core SW104-ND11 (37°01′57″N, 13°10′54″E – water depth of 475 m) is located in the northwestern part of the Sicily Channel (Fig. 1) and was recovered during NEXTDATA 2014 oceanographic expedition on board the R/V CNR-Urania. The sedimentary sequence has been recovered using the SW104 gravity corer system, which allows the recovery of undisturbed and very well preserved water-sediment interface. The core length is of 119 cm and consists of homogeneous grey hemipelagic sediments.
About 60 selected individuals of Globigerinoides ruber white variety sensu scricto (>125 µm) were used per sample. The tests were crushed under glass slides to open the chambers and carefully cleaned applying a sequence of clay removal, reductive, oxidative and weak acid cleaning steps77. A preliminary cleaning test was performed in order to adjust the cleaning protocol to the sample characteristics. In base to the Al/Ca and Mn/Ca ratios in the samples and also in the residuals of the cleaning steps it was decided that the reductive step was needed (Fig. S1a) but also additional clay removal steps were added in order to reduce the relatively high Al values. Additionally, the appropriate performance of the chosen cleaning protocol was further tested by visual inspection of the interior of the foraminifera walls by SEM-microscopy (Fig. S1b) in some chosen samples previously and after their cleaning. These SEM pictures (Fig. S1b) show an efficient removal of diagenetic and detritic carbonates attached to the foraminifera wall after the application of the described full cleaning protocol. Procedural blanks were routinely measured to discard any potential contamination problem during cleaning and dissolution. Instrumental analyses were performed in an inductively coupled plasma mass spectrometer (ICP-MS Perkin Elmer ELAN 6000) in the Scientific and Technological Centers of the University of Barcelona (CCiT-UB). A standard solution with known elemental ratios was used for sample standard bracketing (SSB) as a correction for instrumental drift. The Mn/Ca and Al/Ca ratios were always measured to identify potential contaminations due to the presence of manganese oxides and/or aluminosilicates77,78,79,80 (Fig. S1a). Measured samples always show Mn/Ca and Al/Ca values below 0.1 mmol/mol and 0.4 mmol/mol respectively (thresholds defined by some previous studies77,78,80) (Fig. S1a). However, for this study those Mn/Ca ratios above 2σ (0.09 mmol/mol; over standard deviations of the averages Mn/Ca values) were removed. Also those Al/Ca ratios with values below 2σ (0.21 mmol/mol). Those values bellow 2σ have been accepted as a good result since any clear covariance with Mg/Ca ratios exist (Fig. S1a). The Mg/Ca ratios were converted in Sea Surface Temperature (SST) values according to the G. ruber calibration from Elderfield and Ganssen (2000)81. The SST error has been estimated after propagating the errors associated to the Mg/Ca ratio measurements and to the used SST calibration according to the procedure described in the supplementary material (S1). Therefore, the resulting propagated errors from the Core SW104-ND11 covers an error range between ±1.35 °C to ±1.54 °C.
The estimated core top Mg/CaG.ruber SST of the study record is 20 ± 1.4 °C and this temperature is consistent with the measured temperature at ca. 20 metres in July 2014 during the sampling (Fig. 1c). This value is also in agreement with the seasonal temperature distribution at ca. 20 metres in the Sicily Channel evaluated with the World Ocean Atlas 201382 dataset that averages measurements from 1955 to 2012. Even the living G. ruber data from Pujol and Vergnaud Grazzini (1995)59 and Mallo et al.83 support the depth habitat of G. ruber in the Sicily Channel at ca. 20 metres and 20 °C.
Some previous studies have argued that Mg/Ca ratios in the Mediterranean were anomalously high due to the effect of high salinities in the region84. However, later studies indicate that such a high Mg/Ca ratios in warm oversaturated waters resulted from the chemical overprint of attached secondary calcite85,86. These observations are also in agreement with later culture experiments that reported relatively low sensitivity of the Mg/Ca ratio to high salinity conditions87. This study has taken particular care in ensuring the removal of any potential digenetic calcite by applying the so-called reductive cleaning step, which has demonstrated to promote the release of secondary calcite overgrowths77,80. The efficiency of the cleaning protocol has further been checked through the visual inspection of the cleaned foraminifera walls by SEM-microscopy (Fig. S1b). In addition, the obtained Mg/Ca ratios in the core top provide coherent SST for the late summer season when G. ruber growths in the region60. They are also coherent with the G. ruber-Mg/Ca values previously reported in the same region for a record covering the last 5 centuries76. Therefore, we can conclude that our Mg/CaG.ruber represents a reliable SST record with no major interferences of diagenetic artefacts in the signal.
The chronology for the uppermost 10 cmbsf (cm below seafloor) of core SW104-ND11 is estimated by radionuclides 210Pb and 137Cs. The 210Pb and 137Cs analysis were carried out at ISMAR-CNR Bologna, according to the procedures reported in Bellucci et al.88. The 137Cs activity-profile does not show a clear peak corresponding to 1963 CE (maximum 137Cs fallout from nuclear testing) in the first few cmbsf, suggesting the presence of a surficial mixed layer from 0 to 7 cmbsf. Considering that 137Cs from tests of nuclear weapons was first detected in the atmosphere in 1954 CE, we assigned this age at 7 cmbsf, where this radionuclide first appears (Fig. S2). Using the 210Pb profile below the mixed layer and applying the CF–CS model (Constant Flux–Constant Sedimentation), we calculate a mean sedimentation rate of 50 yr/cm.
The four AMS 14C radiocarbon analysis on mixed planktonic foraminifera (Tab. S1) were performed at iCONa Lab of the Dipartimento di Scienze e Tecnologie Ambientali dell’Università della Campania and at the laboratory LABEC (Laboratorio di tecniche nucleari per l’Ambiente e i Beni Culturali) of the INFN – Florence. Radiocarbon ages were calibrated using the MARINE13 calibration curves89. The calibrated age ranges are reported in years CE and refer to 2σ.
The age model for core SW104_ND11 was constructed using the Bayesian statistics software Bacon with the statistical package R90,91 for marine sediments using the four AMS 14C radiocarbon ages and the tie-point associated to the first detection of 137Cs in the study core. This age model shows a near-constant sedimentation rate with a mean value of 50 yr/cm from the top down to the base of the study core (Fig. S3).
Regression analysis of temperatures temporal profile
In order to statistically validate the observed trends in temporal profile, temperature values were smoothed by adopting a non-parametric regression approach. The adopted procedure was similar to the one applied in Martı́nez-Boti et al.92 and allow to validate the observed trend by taking into account all sources of uncertainty in temperature values while maintaining fixed the age model. Specifically, a LOESS (LOcally Estimated Scatterplot Smoothing93) fitting was carried out. LOESS belongs to the family of non-parametric fitting procedure and does not require the a priori specification of the relationship between dependent variable and predictors. In LOESS regression the user must specify the so called “span” (or “alpha”) parameter, determining the proportion of observation used in each local regression. The span parameter, ranging between 0 and 1, strongly influence the regression results; higher span values tends to produce too much smoothed curves, thus hiding the informative part of the regression, while lower span values will produce too much wiggling curves thus fitting the noise rather than the dominant signal. In order to identify the optimal span parameter, a v-fold cross validation and a generalized cross-validation (GCV) methods were used. Optimal span values were 0.17 and 0.19 as inferred by v-fold cross-validation and GCV respectively, thus a span value of 0.18 was used in LOESS regression. Finally, once the optimal span values was determined, a Monte Carlo simulation was adopted to take into account all the possible source of variability. In particular, 10000 possible realization of the temperature time series were generated by randomly sampling each Mg/Ca observation within its uncertainty bound (2σ), and computing the new temperature value according to Elderfield and Ganssen81 equation:_Mg/Ca(mmol·mol-1)=0.52exp0.1 T. The LOESS regression was then carried out for each realization and for each data-point the smoothed value was computed as the mean value of the 10000 smoothed values. Accordingly, the 95% CI for the smoothed curve were computed as the 2.5% and 97.5% quantile of the 10000 smoothed values. All the calculations were carried out in R statistical environment (R Core Team, 2019) by using packages fANCOVA94 and paleoMAS95 to perform generalized cross-validation and v-fold cross validation respectively. However, it is important to state that the discussions and interpretations of SST variability documented in the identified climate events (in the study core SW1044-ND11, Sicily Channel) refer to the measured data points in °C; the LOESS regression and its CI obtained by using Monte Carlo simulation in fact, has been used only to validate the main trends observed in the SST data.
The data from Alboran Sea (Rodrigo-Gamiz et al.39), Minorca Basin (Cisneros et al.16), from Aegean Sea (Kontakiotis, 2016; Gogou et al.40,41) and the north Hemisphere temperature reconstruction (Ljungqvist et al.42) were recovered from the original publications.
The Mg/CaG.ruber data from the Sicily Channel are available on Pangaea data publisher.
Giorgi, F. Climate change hot-spots. Geophysical Research Letters 33, L08707 (2006).
Corte-Real, J., Zhang, X. & Wang, X. Downscaling GCM information to regional scales: A non-parametric multivariate regression approach. Climate Dynamics 11, 413–42 (1995).
Xoplaki E. Climate Variability over the Mediterranean. Ph.D. Thesis, University of Bern, 193 pp. Available at, http://sinus.unibe.ch/klimet/docs/phd_xoplaki.pdf (2002).
Lionello, P., Malanotte-Rizzoli, P. & Boscolo, R. (Eds). Mediterranean Climate Variability. Elsevier, pp 438. ISBN: 0-444-52170-4, (2006).
Grauel, A. L., Goudeau, M. L. S., de Lange, G. J. & Bernasconi, S. M. Climate of the past 2500 years in the Gulf of Taranto, central Mediterranean Sea: a high-resolution climate reconstruction based on δ18O and δ13C of Globigerinoides ruber (white). The Holocene 23, 1440–1446 (2013).
Goudeau, M. L. S. et al. Seasonality variations in the Central Mediterranean during climate change events in the Late Holocene. Palaeogeography, Palaeoclimatology, Palaeoecology 418, 304–318 (2015).
Lirer, F. et al. Planktonic foraminifera as bio-indicators for monitoring the climatic changes that have occurred over the past 2000 years in the southeastern Tyrrhenian Sea. Integrative Zoology 9, 542–554 (2014).
Bonomo, S. et al. Reworked Coccoliths as runoff proxy for the last 400 years: The case of Gaeta Gulf (central Tyrrhenian Sea, Central Italy). Palaeogeography, Palaeoclimatology, Palaeoecology 459, 15–28 (2016).
Margaritelli, G. et al. Marine response to climate changes during the last five millennia in the central Mediterranean Sea. Global and Planetary Change 142, 53–72 (2016).
Margaritelli, G. et al. Climatic variability over the last 3000 years in the central – western Mediterranean Sea (Menorca Basin) detected by planktonic foraminifera and stable isotope records. Global and Planetary Change 169, 179–187 (2018).
Jalali, B., Sicre, M. A., Bassetti, M. A. & Kallel, N. Holocene climate variability in the North-western Mediterranean Sea (Gulf of Lions). Climate of the Past 12, 91–101 (2016).
Jalali, B. et al. High-resolution Holocene climate and hydrological variability from two major Mediterranean deltas (Nile and Rhone). The Holocene 1–11 (2017).
Labuhn, I., Finné, M., Izdebski, A., Roberts, N., Woodbridge, J. Climatic Changes and Their Impacts in the Mediterranean during the First Millennium AD. Environment and Society in the Long Late Antiquity, (Late Antique Archaeology 12), (Leiden 2018), Adam Izdebski and Michael Mulryan (eds), 65–88 (2018).
Roberts, N., Brayshaw, D., Kuzucuoglu, C., Perez, R. & Sadori, L. The mid-Holocene climatic transition in the Mediterranean: Causes and consequences. The Holocene 21, 3 (2011).
McCormick, M. et al. Climate Change during and after the Roman Empire: Reconstructing the Past from Scientific and Historical Evidence. Journal of Interdisciplinary History 2, 169–220 (2012).
Cisneros, M. et al. Sea surface temperature variability in the central-western Mediterranean Sea during the last 2700 years: a multi-proxy and multi-record approach. Climate of the Past 12, 849–869 (2016).
Neukom, R., Steiger, N., Gómez-Navarro, J.J., Wang, J., Werner, J.P. No evidence for globally coherent warm and cold periods over the preindustrial Common Era. Nature 571 (2019).
Neukom, R. et al. Consistent multidecadal variability in global temperature reconstructions and simulations over the Common Era. PAGES 2k Consortium. Nature Geosciences 12, 643–649 (2019).
Magny, M. & Combourieu Nebout, N. Holocene changes in environment and climate in the central Mediterranean as reflected by lake and marine records. Climate of the Past 9, 1447–1454 (2013).
Holmgren, K. et al. Mediterranean Holocene climate, environment and human societies. Quaternary Science Reviews 136, 1–4 (2016).
Sadori, L. et al. Climate, environment and society in southern Italy during the last 2000 years. A review of the environmental, historical and archaeological evidence. Quaternary Science Reviews 136, 173–188 (2016).
Büntgen, U. & Tegel, W. European tree-ring data and the Medieval Climate Anomaly. PAGES 19, 14–15 (2011).
Büntgen, U. et al. Cooling and societal change during the Late Antique Little Ice Age from 536 to around 660 AD. Nature Geoscience 9, 231–236 (2016).
Hodell, D., Brenner, M., Curtis, J. H. & Guilderson, T. Solar forcing of drought frequency in the Maya lowlands. Science 292, 1367–1370 (2001).
deMenocal, P. B. & Peter, B. Cultural Responses to Climate Change During the Late Holocene. Science 292, 667–673 (2001).
IPCC, 2018: Summary for Policymakers. In: Global warming of 1.5 °C. An IPCC Special Report on the impacts of global warming of 1.5 °C above pre-industrial levels and related global greenhouse gas emission pathways, in the context of strengthening the global response to the threat of climate change, sustainable development, and efforts to eradicate poverty. World Meteorological Organization, Geneva, Switzerland, 32, V. Masson-Delmotte, et al. (eds.) (2018).
Fischer, H. et al. Palaeoclimate constraints on the impact of 2 °C anthropogenic warming and beyond. Nature Geoscience 11, 474–485 (2018).
Faust, J. C., Fabian, K., Milzerc, G., Giraudeau, J. & Kniesa, J. Norwegian fjord sediments reveal NAO related winter temperature and precipitation changes of the past 2800 years. Earth and Planetary Science Letters 435, 84–93 (2016).
Pinardi, N. et al. Mediterranean Sea large-scale low-frequency ocean variability and water mass formation rates from 1987 to 2007: A retrospective analysis. Progress in Oceanography 132, 318–332 (2015).
La Violette, P. E. The western Mediterranean circulation experiment (WMCE): introduction. Journal of Geophysical Research 95, 1511–1514 (1990).
Viudez, A., Tintore, J. & Haney, R. L. Circulation in the Alboran Sea as determined by quasi-synoptic hydrographic observations. Part 1. Three-dimensional structures of the two anticyclonic gyres. Journal of Physical Oceanography 26, 684–705 (1996).
Viudez, A. & Tintore, J. Time and space variability in the eastern Alboran Sea from March to May 1990. Journal of Geophysical Research 100, 8571–8586 (1995).
Millot, C. Circulation in the western Mediterranean Sea. Oceanologica Acta 10, 143–149.4 (1987).
Krom M.D., Groom S. & Zohary T. The Eastern Mediterranean — Biogeochemistry of Marine Systems. Blackwell Publishing (2003).
Bethoux, J.-P. Mean water fluxes across sections in the Mediterranean Sea, evaluated on the basis of water and salt budgets and of observed salinities. Oceanologica Acta 3 l, 79–88 (1980).
Robinson, A. R. et al. The Atlantic Ionian Stream. Journal of Marine Systems 20, 129–156 (1999).
Bèranger, K. et al. The dynamics of the Sicily Strait: A comprehensive study from observations and models. Deep-Sea Research Part II 51, 411–440 (2004).
Incarbona, A. et al. Holocene millennial-scale productivity variations in the Sicily Channel (Mediterranean Sea). Paleoceanography 23, PA3204 (2008).
Rodrigo-Gámiz, M., Martínez-Ruiz, F., Rampen, S. W., Schouten, S. & Sinninghe Damsté, J. S. Sea surface Temperature variations in the western Mediterranean Sea over the last 20 kyr: A dual-organic proxy (UK′37 and LDI) approach. Paleoceanography 29, 87–98 (2014).
Kontakiotis, G. Late Quaternary paleoenvironmental reconstruction and paleoclimatic implications of the Aegean Sea (eastern Mediterranean) based on paleoceanographic indexes and stable isotopes. Quaternary International 401, 28–42 (2016).
Gogou, A. et al. Climate variability and socio-environmental changes in the northern Aegean (NE Mediterranean) during the last 1500 years. Quaternary Science Reviews 136, 209–228 (2016).
Ljungqvist, F. C. A new reconstruction of temperature variability in the extra-tropical Northern Hemisphere during the last two millennia. Geografiska Annaler Series A: Physical Geography 92A, 339–351 (2010).
Jalali, B. et al. Deltaic and coastal sediments as recorders of Mediterranean regional climate and human impact over the past three millennia. Paleoceanography and Paleoclimatology 33, 579–593 (2018).
Kaniewski, D. et al. Late Second-Early First Millennium BC abrupt climate changes in coastal Syria and their possible significance for the history of the Eastern Mediterranean. Quaternary Research 74, 207–21 (2010).
Di Rita, F. et al. Late Holocene forest dynamics in the Gulf of Gaeta (central Mediterranean) in relation to NAO variability and human impact. Quaternary Science Reviews 179, 137–152 (2018).
Di Rita, F. et al. Holocene forest dynamics in central and western Mediterranean: periodicity, spatio-temporal patterns and climate influence. Scientific Reports 8, 8929 (2018).
Behringer, W. A Cultural History of Climate. London, Polity Press, ISBN: 9780745645292, (2009).
Weiss, B. The decline of the Late Bronze Age civilization as a possible response to climate change. Climatic Change 4, 173–198 (1982).
Knapp, B. & Manning, S. W. Crisis in Context: The End of the Late Bronze Age in the Eastern Mediterranean. American Journal of Archaeology 120, 99–149 (2016).
Drake, B. L. The influence of climatic change on the Late Bronze Age Collapse and the Greek Dark Ages. Journal of Archaeological Science 39, 1862–1870 (2012).
Finnè, M. et al. Late Bronze Age climate change and the destruction of the Mycenaean Palace of Nestor at Pylos. PLoS ONE 12(12), e0189447 (2017).
Zolitschka, B., Behre, K. E. & Schneider, J. Human and climate impact on the environmental as derived from colluvial, fluvial and lacustrine archives-examples from the Bronze Age to the Migration period, Germany. Quaternary Science Reviews 22, 81–100 (2003).
Kotthoff, U. et al. Reconstructing Holocene temperature and salinity variations in the western Baltic Sea region: a multi-proxy comparison from the Little Belt (IODP Expedition 347, Site M0059). Biogeosciences 14, 5607–5632 (2017).
Shaw, B.D. Climate, Environment, and History: the Case of Roman North Africa. Chap. 16 [in] T. M. L. Wigley, M. Ingram, and G. Farmer eds., Climate and History: Studies in Past Climates and their Impact on Man. Cambridge University Press, pp. 379-403 (1981).
Mayewski, P. A. et al. Holocene climate variability. Quaternary Research 62, 243–255 (2004).
Bond, G. et al. Persistent solar influence on North Atlantic climate during the Holocene. Science 294, 2130–2136 (2001).
Sicre, M. A. et al. A 4500- year reconstruction of sea surface temperature variability at decadal time-scales off North Iceland. Quaternary Science Reviews 27, 2041–2047 (2008).
DeMenocal, P., Ortiz, J., Guilderson, T. & Sarnthein, M. Coherent high- and lowlatitude climate variability during the Holocene warm period. Science 288, 2198–2202 (2000).
PAGES 2K Consortium. Continental-scale temperature variability during the past two millennia. Nature 6, 339–346 (2013).
Pujol, C. & Vergnaud Grazzini, C. Distribution patterns of live planktic foraminifers as related to regional hydrography and productive systems of the Mediterranean Sea. Marine Micropaleontology 25, 187–217 (1995).
Català, A., Cacho, I., Frigola, J., Pena, L. D. & Lirer, F. Holocene hydrography evolution in the Alboran Sea: a multi-record and multi-proxy comparison. Climate of the Past 15, 927–942 (2019).
Harper, K. The Fate of Rome, Climate, Disease, and the End of an Empire. The Princeton History of the Ancient World Series. Princeton: Princeton University Press (2019).
Noti, R. et al. Mid- and late-holocene vegetation and fire history at Biviere di Gela, a coastal lake in southern Sicily, italy. Vegetation History and Archaeobotany 18, 371–387 (2009).
Bisculm, M. et al. Holocene vegetation and fire dynamics in the supra-mediterranean belt of the Nebrodi Mountains (Sicily, Italy). Journal of Quaternary Science 27, 687–698 (2012).
Di Rita, F. D. & Magri, D. An overview of the Holocene vegetation history from the central Mediterranean coasts. Journal of Mediterranean Earth Science 4, 35–52 (2012).
Gibbon, E. The History of the Decline and Fall of the Roman Empire. Strahan & Cadell, London (1776).
Momigliano, A. La caduta senza rumore. Ed. scientifiche italiane, Italy (1973).
Lamb, H.H. Climate: Present, Past and Future vol. 2. Methuen & Co, London (1977).
Jones, P., Mann, D. & Mann, M. E. Climate over past millennia. Reviews of Geophysics 42, RG2002 (2004).
Mann, M. E. et al. Global signatures and dynamical origins of the Little Ice Age and Medieval Climate Anomaly. Science 326, 1256–1260 (2009).
Piva, A., Asioli, A., Trincardi, F., Schneider, R. R. & Vigliotti, L. Late Holocene climate variability in the Adriatic Sea (Central Mediterranean). The Holocene 18, 153–167 (2008).
Incarbona, A. et al. The Impact of the Little Ice Age on Coccolithophores in the Central Mediterranea Sea. Climate of the Past 6, 795–805 (2010).
Sicre, M. A. et al. Sea surface temperature variability in the North Western Mediterranean Sea (Gulf of Lion) during the Common Era. Earth and Planetary Science Letters 456, 124–133 (2016).
Marullo, S., Artale, V. & Santoleri, R. The SST Multidecadal Variability in the Atlantic–Mediterranean Region and Its Relation to AMO. Journal of Climate 24, 4385–4401 (2011).
Versteegh, G. J. M., de Leeuw, J. W., Taricco, C. & Romero, A. Temperature and productivity influences on U37 K0 and their possible relation to solar forcing of the Mediterranean winter. Geochem. Geophys. Geosyst. 8, Q09005 (2007).
Incarbona, A. et al. Mediterranean circulation perturbations over the last five centuries: Relevance to past Eastern Mediterranean Transient-type events. Scientific Reports 6, 29623 (2016).
Pena, L. D., Calvo, E., Cacho, I., Eggins, S. & Pelejero, C. Identification and removal of Mn-Mg-rich contaminant phases on foraminiferal tests: implications for Mg/Ca past temperature reconstructions. Geochemistry Geophysics Geosystems 6, Q09P02 (2005).
Barker, S., Cacho, I., Benway, H. & Tachikawa, K. Planktonic foraminiferal Mg/Ca as a proxy for past oceanic temperatures: A methodological overview and data compilation for the Last Glacial Maximum. Quaternary Science Reviews 24, 821–834 (2005).
Lea, D. W., Pak, D. K. & Spero, H. J. Climate impact of late Quaternary equatorial Pacific sea surface temperature variations. Science 289(5485), 1719–1724 (2000).
Pena, L. D. et al. Characterization of contaminant phases in foraminifera carbonates by electron microprobe mapping. Geochemistry, Geophysics, Geosystems 9, Q07012 (2008).
Elderfield, H. & Ganssen, G. Past temperature and 18O of surface ocean waters inferred from foraminiferal Mg/Ca ratios. Nature 405, 442–445 (2000).
Boyer, T. & Mishonov, A. World Ocean Atlas 2013 Product Documentation (2013).
Mallo, M., Ziveri, P., Mortyn, P. G., Schiebel, R. & Grelaud, M. Low planktic foraminiferal diversity and abundance observed in a spring 2013 west–east Mediterranean Sea plankton tow transect. Biogeosciences 14, 2245–2266 (2017).
Ferguson, J. E., Henderson, G. M., Kucera, M. & Rickaby, R. E. M. Systematic change of foraminiferal Mg/Ca ratios across a strong salinity gradient. Earth and Planetary Science Letters 265, 153–166 (2008).
Hoogakker, B. A. A., Klinkhammer, G. P., Elderfield, H., Rohling, E. J. & Chris Hayward, C. Mg/Ca paleothermometry in high salinity environments. Earth and Planetary Science Letters 284, 583–589 (2009).
van Raden, U. J., Groeneveld, J., Raitzsch, M. & Kucera, M. Mg/Ca in the planktonic foraminifera Globorotalia inflata and Globigerinoides bulloides from Western Mediterranean plankton tow and core top samples. Marine Micropaleontology 78(3-4), 101–112 (2011).
Hönisch, B. et al. The influence of salinity on Mg/Ca in planktic foraminifers – Evidence from cultures, core-top sediments and complementary d18O. Geochimica et Cosmochimica Acta 121, 196–213 (2013).
Bellucci, L. G. et al. 210Pb and 137Cs as chronometers for salt marsh accretion in the Venice Lagoon — links to flooding frequency and climate change. Journal of Environmental Radioactivity 97, 85–102 (2007).
Reimer, P. J. et al. IntCal13 and Marine13 radiocarbon age calibration curves 0–50,000 years cal BP. Radiocarbon 55(4), 1869–1887 (2013).
Blaauw, M. & Christen, J. A. Flexible Paleoclimate Age-Depth Models Using an Autoregressive Gamma Process. Bayesian. Analysis 6(3), 457–474 (2011).
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria., https://www.R-project.org/ (2019).
Martínez-Boti, M. A. et al. Boron isotope evidence foroceanic carbon dioxide leakage during the last deglaciation. Nature 518, 49–54 (2015).
Cleveland, W.S., Grosse, E. & Shyu, W.M. Local regression models. Chapter 8 of Statistical Models in S eds J.M. Chambers and T.J. Hastie, Wadsworth & Brooks/Cole (1992).
Wang, X,F. fANCOVA: Nonparametric Analysis of Covariance. R package version 0.5-1, https://CRAN.R-project.org/package=fANCOVA (2010).
Correa-Metrio, A., Urrego, D.H., Kenneth, R. Cabrera and Mark B. Bush. paleoMAS: Paleoecological Analysis. R package version 2.0-1. (2012).
The core SW104-ND11 has been collected by IAMC-CNR (Napoli) aboard of the CNR-Urania vessel during the oceanographic cruise NEXTDATA-2014. This research has been financially supported by the Project of Strategic Interest NextData PNR 2011–2013 (http://www.nextdataproject.it/). This study has received funding from the European Research Council (ERC) through the ERC-Consolidator project TIMED (Grant Agreement No. 683237). I.C. and A.C. also thank Generalitat de Catalunya for grant 2017 SGR 315 to GRC Geociències Marines and for the ICREA-Academia programme support to I.C. Ester Garcia-Solsona is thanked for her assistance in the error propagation calculation.
The authors declare no competing interests.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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/.
About this article
Cite this article
Margaritelli, G., Cacho, I., Català, A. et al. Persistent warm Mediterranean surface waters during the Roman period. Sci Rep 10, 10431 (2020). https://doi.org/10.1038/s41598-020-67281-2