Abstract
The forcing mechanisms responsible for centennial to millennial variability of mid-latitude storminess are still poorly understood. On decadal scales, the present-day geographic variability of North-Atlantic storminess responds to latitudinal shifts of the North-Atlantic westerly wind-belt under the prime control of the North-Atlantic Oscillation (NAO). An equivalent mechanism operating at centennial to millennial time scales during the Holocene is still to be ascertained, especially owing to the lack of high-resolution and continuous records of past-storminess extending far enough in time. Here we present a reconstruction of past storminess activity based on a high-resolution record of wind-blown sand retrieved from a near-coastal wetland. Our record extends back to ca. 10,000 B.P. and allows to continuously document fluctuations in the frequency of Holocene storm-force winds at our study-site at a mean high temporal resolution of 40 years. Large similarities between our record and palaeo-oceanographic records of Holocene climate changes in the North-Atlantic suggest that our past-storminess record reproduces a signal of significance for the North-Eastern Atlantic realm. We find that Holocene North-Atlantic storminess is dominated by robust millennial (≈2,200-year) to centennial (≈450, 300 and 200-year) periodicities. These changes in storminess were accompanied by changes in the precipitation regimes over northern Europe, evidencing large-scale shifts in the latitudinal positions of the Atlantic westerlies akin to present-day NAO patterns. We propose that these shifts originate from changes in the position and extent of the Azores high-pressure system and Polar vortex, as supported by climate model simulations. Finally, we demonstrate that enhanced zonal storminess activity over the North-Atlantic was the driver of millennial and centennial-scale changes in North-Atlantic oceanic circulation, while ocean dynamics most likely influenced back the atmospheric circulation at millennial time-scales. This may vouch for the instrumental role played by North-Atlantic storminess in triggering abrupt climate change at centennial scales during the Holocene.
Similar content being viewed by others
Introduction
Westerly wind-belts and associated storminess are prime agents of North-Atlantic mid-latitude climate modulation. During recent years, efforts have been made to predict changes in storminess activity related to present and future climate-change notably through the use of coupled ocean-atmosphere models. However, there is currently no clear consensus on the future latitudinal position of extra-tropical storm tracks. Indeed, climate projections alternatively suggest that a northward or a southward shift of extra-tropical storm-tracks over the North-Atlantic may be experienced1, bearing large consequences for the builing of reliable strategies of adaptation to future climate change in North-western Europe.
At present, a large proportion of the climate variability in North-western Europe is controlled by the North Atlantic Oscillation (NAO)2. The NAO especially influences the latitudinal position of the westerly wind-belt and associated storminess over the North-Atlantic. A common paradigm is that positive phases of the NAO (NAO+) induce a zonal northerly located Atlantic Westerly Jet that drives storm-tracks and humid air masses over northern Europe, while the westerly circulation over southern Europe is blocked3. On short (monthly to decadal) time-scales, tight connections have been evidenced between NAO variability, activity of westerly winds and oceanic circulation in the North-Atlantic. NAO+ situations were for instance shown responsible for meridian displacements of the subpolar (SPG) and subtropical (STG) gyres in connection with the onset of a northward relocated and strengthened westerly wind stress (e.g.4). Modelling studies5 suggested that prolonged NAO+ situations and associated wind anomalies over the North Atlantic could have been instrumental in triggering of the Holocene abrupt cold climatic anomalies through oceanic feedbacks6, thus promoting storminess to a central role in the modulation of the Earth’s climate. However, the control exerted by the NAO on the position of North-Atlantic westerly storm tracks as well as the mechanisms controlling the NAO variability have been repeatedly discussed7 and are currently challenged. An emerging view is that the leading role of the NAO may have been overestimated, and that other modes of atmospheric variability such as the East-Atlantic pattern (EA), responsible for atmospheric blocking events over the North-eastern Atlantic, and the Scandinavian (SCAND) pattern may play important roles in controlling climate variability and storminess over North-Western Europe8,9,10. A three-poles linking blocking events (as part of the EA pattern), Total Solar Irradiance and Subpolar Gyre (SPG) dynamics has notably been suggested to have dominated atmospheric and climate variability at decadal to centennial scale over the last millennium8.
Better understanding the links between NAO variability and storminess over long time-scales therefore remains a central question, as it would allow to refine the projections produced by climate-models and to better address the challenges raised by ongoing climate-change. Up to now, the simulations of climate-models have been mostly validated against records of recent and historical meteorological data ranging from years to a few centuries in length11. Validation by robust storminess data extending over pluri-millennial time-periods is still lacking to ascertain the long-term spatio-temporal variability of storminess predicted by the models. How and under which control North-Atlantic storminess has varied in the past at shorter centennial-scales is still poorly understood. A reason for this may be that traditional wave-induced palaeo storminess indicators can hardly produce storm reconstructions of sufficiently wide spatial significance and of sufficiently high temporal resolution over the Holocene12. During the last decade however, promising reconstructions of past storminess have been derived from proxies of past aeolian activity (aeolian sand influx in near-coastal and perched peat-bogs), notably in Northern Europe and Scandinavia13,14,15 but also worldwide16. The use of aeolian sand influxes preserved within low-energy organic sedimentary sinks offers numerous advantages over the proxies of storminess derived from wave-influenced coastal stratigraphies. On the one hand, peatlands are very much stable sedimentary environments and offer the opportunity (i) to derive continuous records that can be precisely dated and (ii) to obtain perfectly synchronous palaeo-environmental data on the same core17. On the other hand, the process of aeolian sand transport is controlled to a high degree by the relationships between wind-energy and the size of the material available for transport, so that it has become possible to evaluate the occurrence of winds of different intensities over time13. Despite promising results, the use of aeolian proxies of past-storminess activity has for now remained rather limited.
Here, we investigate the links between storminess in North-western Europe and ocean-atmosphere variability in the North Atlantic during the Holocene on the basis of past influxes of aeolian sand recorded in the sedimentary sequence of Filsø, a coastal wetland of Western Denmark (Fig. 1, Supp. Info. 1). The record extends back to the early Holocene (ca. 10,200 B.P.), making it the longest continuous and high-resolution record of past storminess for North-western Europe. Combined with palaeo-oceanographic records available for the North Atlantic (Fig. 1), this unique record allows us to display and discuss the forcing mechanisms behind the variability of North-eastern Atlantic Holocene wind climate.
Regional settings
The study area forms an ideal candidate for storminess reconstruction on the basis of wind-transported material. At an early stage of evolution, Filsø formed upon a glacio-fluvial topography and is considered to have experienced marine conditions in its deepest parts during the transgressive maximum of the regional mid-Holocene relative sea-level (around 7,000 B.P.18,19,20). From 7,000 B.P., the connection to the open North Sea closed rapidly by barrier-spit development driven by an intense southward-orientated longshore wave transport19,21 fed with the products of the erosion of local moraines, proglacial glacio-fluvial deposits and smoother periglacial landscapes18,21. A low-energy organic depositional environment settled fast behind this barrier22. At the same time, stacked beach ridge systems formed and prograded towards the west (seaward) and south under a continued and high sediment supply, completely isolating Filsø from direct marine influences19,21. Thick aeolian sand dunes subsequently developed upon the wave-built deposits and covered most of the area. First local occurrences of aeolian sand deposition were recorded south of our study site ca. 15,000 B.P., most probably originating from the reworking of glacio-fluvial sands of Saalian and Weichselian ages (as part of the European Sand Belt22,23,24). Periods of intense large-scale inland aeolian sand movement and coastal dune formation were then observed in western Denmark at around 8,000 B.P., 6,000 B.P. and ca. 4,300, 3,450, 2,800 and 2,500 B.P.19,22,23,25,26. At present, prevailing winds blow from the south-western to western sectors along the west coast of Denmark27. Storm activity is rather high, with winds of gale-force (8 on the Beaufort scale) and above occurring at annual frequencies of 5–10%. The drift potential reaches 1510 vector units (wind data in knots), thus making the wind-climate a climate of high-energy28.
Stratigraphic record
The stratigraphy of the site was reconstructed based on two cores (F06 and F02) retrieved from the western part of Filsø (55°42′N, 8°13′E, Fig. 1 and Supp. Info. 1 and 2) which were linked using µ-XRF data (Supp. Info. 2). The chronology was built using ten AMS 14C dates and two optically-stimulated luminescence (OSL) ages (Supp. Info. 2 and 3).
The reconstruction of past-storminess shown in the present study was derived from core F06 only. The inferred age of the sediment in the core F06 spans from 10,200 B.P. to about 2,500 B.P. It was sampled with a mean resolution of 40 years (i.e. 1-cm sampling interval; Fig. 2, Supp. Info. 2). Core F06 can be subdivided into four main sedimentary units (Fig. 2): (i) Bottom units I and II (10,200-6,800 and 6,800-4,300 B.P., respectively) are composed of well-decomposed gyttja deposits progressively gaining organic content towards the top. Based on the paleo-geographical context of the site, the plant macrofossil assemblages for this unit (Supp. Info. 4) and with reference to other reconstructions from the area22, we believe that unit I deposited in a small freshwater pond which surely developed within a depression in the undulated antecedent low-gradient glacio-fluvial topography. This pond evolved into a back-barrier freshwater lacustrine environment (unit II), that the plant macrofossils assemblages show to have remained very shallow with water-depths most likely lower than 1 m29 (Supp. Info. 4). At a depth of 1.72 m below core top (−1.43 m DVR 90), there is a sharp transition to a 15 cm-thick bed of dune-sand (unit III, Fig. 2, Supp. Info. 2) which was dated to 4100 ± 200 B.P. and undoubtedly corresponds to the period of enhanced aeolian activity and intense dune movement identified for the same period along the entire western coast of Denmark26. Above this sand layer, a 1 m-thick layer of a well-humified Sphagnum-peat (unit IV, Fig. 2) is encountered which is primarily composed of Sphagnum teres and Sphagnum palustre (Supp. Info. 4). The topmost 0.55 m of the core were disturbed by agricultural practices as evidenced by a clear ploughing surface (Supp. Info. 2) and were not considered in the investigation (water level in Filsø was artificially lowered for land reclamation around the mid-1800s until 2012 when the site was turned back into a lake).
Source of sand
There are several means by which sand-sized grains can be deposited within terrestrial organic sedimentary sinks. These include aeolian transport, remobilization of lake-shore material, fluvial input or mass wasting17. Our study site being concerned with two small rivers entering the site in its south-eastern termination, the source of the sand material found within core F06 organic sediments had to be tested for its aeolian origin so that the sand influx could be subsequently used as a trustful proxy of past aeolian activity and westerly storminess at the site. To do so, we used a combination of geochemical (µ-XRF) and paleo-ecological indicators (see the method section).
Modern coastal, dunes and river surface sands were sampled around the study site (Supp. Info. 5) and analysed for their composition in major and trace elements. The results revealed significant geochemical differences between coastal (beach and dunes) and river sand populations. The former are characterized by high Si and K contents, while the latter are marked by high Fe and Mn contents (Supp. Info. 5). Thus, the very good correlation observed between the ASI at our site and peaks in K/Fe ratio (Spearman correlation = 0.46, p < 0.0001 between K/Fe and Total ASI across unit I and II) suggests that the sand found in the core is primarily composed of coastal sand (Fig. 2). This statement is reinforced by the paleo-geographical information obtained from the plant macro-fossils assemblages which show that the site evolved into a very shallow lake (water-depth <1 m, 22, Supp. Info. 4) between ca. 7,000 and 4,300 B.P. It is thus unlikely that wave energy in the lake could have been high enough to provoke an erosion of the shores and large sand redistributions. As for the upper part of the sequence, the Sphagnum peat-bog character of the depositional environment (Supp. Info. 4) indicates that the occurrence of high-velocity water circulation can be ruled out at the coring site. The watershed of our study site is characterized by very gentle slopes (see elevations on Supp. Info. 1) and the two rivers of have very low transport competences. This makes it very unlikely that fluvial transport could have produced current velocities high enough to carry grain sizes >125 µm in suspension to the coring location. Based on these considerations, we infer that the sedimentary sequence of core F06 continuously recorded influx of aeolian sand blown downwind by storm-winds from NW to SW sectors between 10,200 and 2,500 B.P.
Storminess record
The sedimentary record of core F06 shows repetitive fluctuations between periods of higher and lower input of aeolian sand to the site between 10,200 and 2,500 B.P. The bottom part of the sequence (unit I, ca. 10,200-6,800 B.P.), is first characterized by a marked succession of two wide peaks in both total and >125 µm ASI separated by a clear low. In this early stage of lake basin development, the sand may have mostly originated from the reworking of local glacio-fluvial sand sources (observed at the base of core F02, see Supp. Info. 2). From ca. 6,800 B.P. and onwards (unit II), the content of both total and >125 µm ASI progressively and steadily increased until the occurrence of the massive 4,300 B.P. sand invasion (unit III). The well sorted and rounded sand-grain population shows that the source of aeolian sand may have become dominated by coastal dunes during this period, following the establishment of coastal barriers directly west of the site after the mid-Holocene deceleration of the relative sea-level rise19,21. The gradual increase of aeolian sand influx across this unit II is most likely the expression of the transgressive coastline and coastal dunes systems moving eastward (and upwards), thus getting closer to the coring-site. The most prominent features of unit II are the well-defined and cyclic alternations between periods of very high (above ca. 75%) and low (25%) proportion of sand >125 µm in the total ASI (Fig. 2). Following the deposition of the sandy unit III (4,300-4,100 B.P.) and throughout most of unit IV, we observe a stable decrease of ASI, and notably of the ASI >125 µm before it rises again sharply around 3,200 B.P. An abrupt drop in the proportion of coarse sand in the total ASI to the site occurs at the transition between unit III and unit IV, with mean percentages of sand >125 µm dropping at around 15-20% and peaks hardly reaching 50% of the total ASI. Clear and well-defined peaks in the percentage of sand >125 µm nonetheless appear, still showing very good consistency with the peaks in total sand influx and suggesting periods of intensified storminess activity. In opposition to the steady increase of ASI observed throughout unit II, which we link to the mid-Holocene transgression, we argue that the gradual decrease in ASI shown for unit IV may be the expression of an intense westward (seaward) progradation undergone by the exposed coastline after 4,300 cal. B.P.19,21.
The sedimentary record from Filsø can be divided into eight main periods of increased storminess (“Filsø Storm Periods” FSP) centred on 9,700-9,100 (FSP1), 8,800-8,600 (FSP2), 8,300-7,100 (FSP3), 6,900-6,700 (FSP4), 6,400-5,500 (FSP5), 5,100-4,700 (FSP6), 4,400-3,800 (FSP7) and 3,300-2,800 B.P. (FSP7) (Fig. 5A). These results are in general agreement with the Holocene Storm Periods previously proposed for NW Europe30, but show two additional periods of marked storminess for the early Holocene (FSP1 and FSP2, Fig. 5-I). Furthermore, the continuous high-resolution storm record produced in the present study allow us to go further by depicting high-frequency centennial-scale changes in storm activity in the north-Atlantic during the Holocene which were never described with such accuracy.
A multi-scale variability of Holocene North-Atlantic storminess activity
To evaluate the modulations of the north-eastern Atlantic storm activity over long- (millennial) and shorter-term (pluri-decadal to centennial) time periods, we performed time-series analyses on our records (Fig. 3 and methods). Spectral analyses were made to evidence the dominating wavelengths comprised in our signal (Fig. 3A-III,B-III), while Evolutive Harmonic Analyses (EHA) allowed us to follow changes in these periodicities through time (Fig. 3A-IV,B-IV). Spectral analysis of the lower frequencies shows storminess to be composed of five periodicities of ≈2,200, ≈1,350, ≈590, ≈450 and ≈200 years at the 90% significance level (Fig. 3A-III), the ≈2,200-yr period showing the strongest power. As higher frequencies are concerned, strong wavelengths of ≈600-, ≈500-, ≈450-, ≈280-, ≈200-, ≈150- and ≈80-yr periods are evidenced (Fig. 3B-III), demonstrating the persistence of pluri-decadal to centennial scale modulations of North Atlantic storminess during the Holocene. EHA show variable stationarity of these millennial and centennial periodicities. A relatively wide band of quasi-periodic variability containing the dominant ≈2,200-yr period clearly dominates the signal over most of the Holocene (Fig. 3A-IV). A notable non-stationarity of the spectral signature nonetheless appears between ≈6,000 and ≈4,500 B.P., when the signal observes a marked excursion toward shorter periods of ≈1,300 years (Fig. 3A-IV). This temporary excursion towards a ≈1,300-yr mode of periodicity at the mid-Holocene seems accompanied by a divergence in the high-frequencies wavelengths characterizing the storminess signal during the early-Holocene (Fig. 3B-IV). Indeed, from ca. 10,200 to ca. 6,800 B.P, EHA show our storm record to be dominated by robust quasi-stationary ≈400 to 500-yr and ≈300-yr wavelengths (Fig. 3B-IV). A clear shift operates around 6,800-6,500 B.P. when the ≈280-yr periodicity is observed to progressively vanish while the ≈500-yr period evolves towards a strong ≈650-yr and 450 ones.
Pluri-millennial periodicities have been repeatedly documented in Holocene palaeo-environmental studies of the North Atlantic domain. The origins of these latter have been intensively discussed with regards to external (solar) and internal (oceanic/atmospheric) variability30,31,32,33,34,35. Periodicities of ca. 2000-2500 years have been proposed to be strongly correlated with variations in the total solar irradiance34,36. The persistence of a clear and stationary ≈2,200-yr periodicity in our storm-record over most of the Holocene (Fig. 3A-IV) thus bring in contention that solar forcing may have been the main driver of the functioning of the climate system and storminess activity at a ca. 2,500-yr pace over the Holocene. Shorter periodicities of ≈1,300–1,500 yrs were suggested to be either sub-harmonics of the 2,500- and 1,000-yr cycles32 or the expression of a massive reorganization of North Atlantic thermohaline circulation at that time37. The non-stationarity of the low-frequencies we observe within our storm record between ca. 6,000-4,500 B.P., characterized by an excursion of the signal towards a shorter period of ≈1,300 yrs may thus be the expression of a temporary overtake of the influence of solar forcing on climate by the mid-Holocene stabilization ocean-atmosphere internal modes of variability as suggested by previous studies30,37,38,39,40.
Interpreting the origin of the pluri-centennial periodicities evidenced by spectral analyses in our storminess record is more difficult, owing to consequent trade-offs between the several internal forcing factors susceptible to have been active at these time-scales41,42,43. Periodicities of around ≈500, ≈400, ≈300 and ≈200 yrs have all been suggested to reflect fluctuations in the total solar irradiance42,43,44,45,46. The EHA especially depicts a great stationarity of the ≈200-yr period in our storm record throughout most of the Holocene, suggesting that it could have been the main mode of variability of North-Atlantic storminess at centennial scale. The ≈200-yr period is traditionally attributed to be the expression of the Suess-DeVries solar cycle (208-yr period44,47). How variations in solar irradiance are transmitted to atmospheric and surface wind circulation remain for now rather unclear. Solar irradiance may be transmitted to the westerly jet activity throughout a “top-down” coupling between the stratosphere and the troposphere48. Models of the effect of solar ultraviolet irradiance on North-hemisphere winds showed that this may be instrumental in fostering baroclinic gradients responsible for NAO variability at decadal scales, and demonstrated that solar minima may translate into strengthened surface mid-latitude westerlies activity during the winter season48,49,50. The amplification of solar variability through ocean-feedback patterns have also been suggested to possibly occur by changes in the thermohaline circulation51 and by solar heating of the sea-surface being transmitted to the atmosphere by “bottom-up” sea-air coupling52. Interestingly, the loss of power of the ≈200-yr period between ca. 6,500 and 4,500 yrs B.P. (Fig. 3B) occurs simultaneously with large reorganizations of the oceanic circulation reconstructed from proxy-records for the mid-Holocene, such as the onset of warm subsurface water transport from the tropics to the subpolar regions53 and the weakening of subpolar gyre circulation37, suggesting that these mechanisms may be linked to each other either by control links or by being the product of a common driver.
Zonality of Holocene North-Atlantic storminess controlled by atmospheric variability
Synchronous variations between our record of storminess and proxy-records of precipitations available from Scandinavia and southern Europe confirm that variations in storminess activity may have originated from persistent north to south fluctuations in the latitudinal position of North Atlantic westerly winds and storm-tracks during the Holocene (“seesaw” patterns). Periods of increased storminess documented by our record show for instance good agreement with the periods of increased humidity identified in south-western Norway for the Holocene54, implying that both regions recorded concurrent large-scale changes in the delivery of humid air masses during the Holocene (Fig. 5A). On the contrary, the abrupt drop in storminess shown in our record around 6,500 B.P. was accompanied by a rapid onset of more humid conditions over Portuguese wetlands as well as by a shift from an extra-tropical to a tropical provenance of the humid air masses affecting the Iberian Peninsula55. Thus, we infer that the marked drop of storminess observed at this period was controlled by a southward relocation of the Northern Atlantic westerly wind belt. In reference to the mechanisms behind present-day NAO variability, it is reasonable to hypothesize that such changes in the position of Holocene North-Atlantic westerly winds may have been produced by large-scale variations in the temperatures and pressure gradients across the North Atlantic under the primary control of changes in the position and extent of the Azores High.
To test this hypothesis, we compared our reconstruction with a transient simulation of the Holocene performed with the Earth system model LOVECLIM56. The simulated 800 hPa wind speed averaged over the North Sea (Fig. 4C and Supp. Info. 6) and over the North-Atlantic (Supp. Info. 6) are in overall agreement with our data. Periods of increased storminess over the North Atlantic also correspond to periods of relative increase of simulated precipitation over Northern Europe (Supp. Info. 6). In addition, the periods of enhanced storminess identified in our proxy record are associated with marked differences between the 800 hPa and 500 hPa geopotential heights around Iceland and the Azores (Supp. Info. 6), pointing towards periods of more positive NAO situations. In comparison to periods of low storminess activity, the periods of increased storminess are predicted by LOVECLIM model simulations to be characterized by prominent positive geopotential anomalies centred north of the Azores and extending over most of the North Atlantic and negative anomalies over southern Greenland (Fig. 4A,B). In other terms, LOVECLIM simulations support our hypothesis that Holocene periods of high storminess activity over Northern Europe were associated with a stronger and more northward-located Azores high-pressure centre and with a strengthened polar vortex, in a configuration resembling present-day positive NAO-situations (Fig. 4D). Between these two strong centres of action, a narrow corridor most likely funnelled the Atlantic Westerly Jet and associated cyclonic perturbations, thus driving storm-tracks towards northern Europe (Fig. 4D).
North-Atlantic storminess and ocean dynamics linkages
Comparisons of our storminess record and palaeo-oceanographic records of the North Atlantic demonstrate the close links that existed between North-Atlantic variability and storminess activity at multi-temporal scales during the Holocene. A striking similarity stands out between our record and changes in the abundance of the planktonic foraminifera Globigerinoides ruber (white) (G.ruber w.) in a deep-marine core taken south of the Azores53 (Fig. 5A-IV,B). G.ruber w. being most abundant in the North-Atlantic Sub Tropical Gyre (STG), its abundance in marine shelf sediments is considered as a proxy for the extent of the STG and for the position of the Azores front (AF), high and low abundances revealing northward and southward position of the STG, respectively53. Periods of increased and decreased storminess evidenced by our storm record are concurrent with respective northward and southward shifts of the location of AF and thus of the STG throughout the Holocene on both long (millennial) and shorter (centennial) time-scales (Fig. 5A-I and A-IV). Multi-taper coherence analyses (CMTM, see methods) performed on our storm record and on the record of the abundance of G.ruber w. show the two series to share common periodicities of ≈ 1,300, ≈ 700, ≈ 400, ≈ 240 yrs, the highest coherence estimates being obtained for the ≈ 1,300, 700 and ≈ 400-yr periods (Fig. 5B). These results echo the same periodicities identified within our storminess record (Fig. 3), thus showing that both storminess and the position of the AF/STG most likely responded synchronously to one or several drivers varying at these periods. Which of atmospheric or oceanic variability was the first driver remains a challenging question to which our data cannot provide conclusive answers. The results of the coherence analyses between storminess and the proxy record of the extent of the STG53 nonetheless show clear positive phasings (of 13, 95 and 28-degrees, respectively) at both millennial and centennial time-scales (Fig. 5B). Such phasing suggests that changes in storminess led STG variability, in line with what is observed nowadays at inter-annual to decadal scales4. Smaller phases are obtained for ≈1,300 and 700-yr periods (13 and 0 degrees, respectively, Fig. 5B) being in favour of a bottom-up forcing of oceanic variability on atmospheric processes at millennial time-scales. A strong pacing at millennial to centennial time-scales is also evident for most of the Holocene between our storminess record and the activity of the sub-polar gyre (SPG) as expressed by differences in the density of the water masses6 (Fig. 5A-VI). During periods of increased storminess, density difference increases, indicating a stronger and more east-west elongated SPG, thus contributing more to the inflow of water to the north-eastern Atlantic through the NAC (Fig. 1). This is mirrored in the fair synchronism that appears at the ca. 2,500-year pace between our data and the strength of the North Atlantic Water inflow into the Norwegian Sea (Fig. 5A-V). CMTM analyses further show these two records to be coherent and in phase at the ≈300-yr period. We propose that strengthened storminess over the North Atlantic would have enhanced the activity of the North-Atlantic Current (NAC, Fig. 1) and brought warmer water to the Nordic seas57. This would have fostered large SST differences between the eastern and western part of this area, in line with the hypothesis that SST differences in the Nordic Seas relates to NAO variability58 and, maybe more importantly, to North-eastern Atlantic storminess activity. All in all, these results support the proposition recently made that persistent positive NAO-like situations and associated enhanced North Atlantic storminess during the Holocene may have been instrumental in amplifying or even triggering abrupt pluri-centennial cold events by strengthening the SPG and causing a freshening of the North Atlantic inflow and, in turn, the slowdown of the Atlantic Meridional Overturning Circulation5,6.
A nearly perfect peak to peak synchronism is also apparent between our record and the stacked record of Ice-Rafted Debris (IRD) reconstructed in the North Atlantic over the 10,200-2,500 B.P. period31 (“Bond events” 6, 5, 4 and 3, Noted “BE” on Fig. 5A-II). The supply of IRD to the North Atlantic was proposed to be due to an extensive cooling at polar and sup-polar latitudes of the Northern Hemisphere that could have strengthened the polar vortex and triggered the drifting of sea-ice, as a result of strong northerly winds31. An alternative explanation associates peaks of IRD in the North Atlantic to a dual strengthening of the East Greenland Current (EGC, Fig. 1) and of the North Atlantic Current (NAC, Fig. 1), under the joint influence of strong northerly winds caused by the presence of a strong Icelandic low-pressure system and of strengthened westerlies59. Multi-taper coherence analyses of our storm record and the input of IRD to the North-eastern Atlantic show the two series to be highly coherent at the 1,500-yr period with a phase lag of 47 degrees (Fig. 5C). Despite the uncertainties in the phasing produced by the CMTM analysis, this positive phasing suggests a leading of storminess over the IRD supply to the North-Atlantic, which is incidentally evident when looking at the two time-series superimposed on each other (Fig. 5C). Our results are therefore in favour of the IRD-events in the North Atlantic to be (at least partially) the result of strengthened and northward-located westerlies and storm-tracks in the North Atlantic during NAO-like situations that would have enhanced the SPG and helped spreading Ice-Rafted Debris over the North- western Atlantic.
Conclusions
Our high-resolution record of aeolian sand influx to a coastal wetland in western Denmark allows the identification of multi-scale periodicities in North-Atlantic storminess over a large part of the Holocene. As such, it brings a significant contribution to the understanding of climate variability in the North-Atlantic realm during the past ten millennia. Our results further corroborate the hypothesis that North Atlantic westerlies and storm-tracks underwent millennial to centennial shifts in their latitudinal position during the Holocene. These changes appear to have been modulated by persistent long- and short-term NAO-like patterns, probably induced by changes in the position and strength of the Azores High pressure system and of the polar vortex. Furthermore, our findings stress that the close synergy observed today between North Atlantic oceanic circulation and the position of the westerlies wind belt in the North Atlantic also operated at millennial to centennial timescales during most of the Holocene. Storminess most likely drove the activity of North-Atlantic oceanic gyres at both millennial and centennial time-scales. This nexus highlights the major role of westerly wind belts not only as the products, but also as forcing-agents of climate variability. Further work would nonetheless be needed to disentangle the respective roles played by internal and external forcing agents in controlling the changes in the configuration of pressure centres and associated wind-shifts over the North-Atlantic.
Methods
Core collection and preparation
Two cores (F02, F06) were taken from a raft in the western part of Filsø (Western Denmark, see map in Supp. Info. 1) in September 2016 using a vibracoring system and aluminium tubes with an outer diameter of 80 mm. Cores were sealed and kept in a cold room at 4 °C until opening. Cores were then cut lengthwise and split in halves under subdued red-light conditions. Sand layers were sampled for OSL dating. The lithology/stratigraphy of the cores was described and pictures were taken.
Dates and age-model
Two OSL dates were established on sand layers present at the bottom and in the middle of core F02. Ten radiocarbon dates were established on plant macrofossils on both F02 and F06 cores (Supp. Info. 3). All dates presented in the manuscript were calibrated using Intcal13 calibration program60 and are reported within a two-sigma confidence interval. Age-models were constructed using Clam age-modelling package in R software61. The age-model used in this study is a compound age-model constructed using ages from both cores F06 and F02 (Supp. Info. 2 and 3). Correlations between the two-cores were drawn on the basis of concurrent stratigraphic markers and peak-to-peak correspondence in K content (Supp. Info. 2). Sedimentation rates were calculated as the difference between two ages at a 1 cm-scale and consequently used for the calculation of the Aeolian Sand Influx index (see below).
Storminess reconstruction
Storm activity was derived from the quantity of aeolian sand deposited within the organic sediment of the wetland. The first half of the core was cut into exact 1-cm thick slices. The dimensions of each slice were measured to obtain volumetric information. Slices were weighted, dried at 105 °C, cooled in a dessicator and re-weighted to obtain water content. Slices were then burnt at 550°, cooled in a dessicator and re-weighted to obtain loss-on-ignition (LOI). Post-LOI remains were washed with 10% hydrochloric acid and wet-sieved with 125 µm and 63 µm meshes to extract minerogenic content. The ≥ 63–< 125 µm (“fine” sand) and the ≥125 µm “coarse” sand fractions were then weighted. Building up on the Aeolian Sand Influx (ASI) method13, here we calculated the ASI by dividing the total sand weight comprised in each core slice by the sedimentation rate and the volume of the slice to obtain a sand influx to the site in g.cm−3.yr−1. Total ASI was calculated for all sand grains >63 µm, while ASI 125 was derived from the influx of sand grain of size >125 µm only. Following available laws of aeolian sand transport62, the influx of sand of grain-size ≥125 µm was considered to sign for storm events of wind speeds ≥22.5 m.s−1 at 10-m high (10 on the Beaufort scale and above). An additional indicator was used that is the percentage of sand of grain-size ≥125 µm comprised in the Total ASI, which was considered to sign for the frequency of storm events of wind speeds ≥22.5 m.s−1 at 10-m high (10 on the Beaufort scale and above). A storm index was then derived from the detrended the percentage of sand of grain-size ≥125 µm.
Plant macro-remains
The paleo-ecological evolution of our site was investigated using plants macro-remains. 1-cm sediment slices were dispersed and sieved at 200 µm. Sieving residuals were investigated wet under the microscope and macro-remains were picked and identified. For each species, macro-remains are reported in term of occurrence at a 1-cm resolution (Supp. Info. 4).
Tracing the sand sources
Core-F06 and F02 as well as modern beach, dune and watershed river domains surrounding our study site were analysed for their geochemical characteristics (Supp. Info. 5) using an ITRAXTM core-scanner (Natural History Museum of Copenhagen) mounted with an Rh tube. As for downcore µ-XRF profiles, these were obtained at 1-mm increments with a current of 30 kV/50 mA and an exposure time of 35 seconds. Modern surface sand samples were analysed as discrete samples within round cups using the same settings. For both core and modern samples, raw XRF results in counts per second were divided by the sum of the Raleigh and Compton components (Rh coh + inc) to account for down-core variations in water and organic content63. The K/Fe ratio was chosen as the sand source indictor following the results given by the PCA investigation of the geochemical characteristics of modern surface sand samples (Supp. Info. 5). K was preferred over Si due to Si being potentially influenced by biogenic silica. Very good correlation between K and Ti ensured us that K was of minerogenic origin.
Statistical analyses
Spectral analyses and wavelet transform analyses were performed on the coarse-sand fraction of the record using the Astrochron packages in R64. The signal was smoothed prior to spectral analysis by detrending the mean and the variance of the time-serie using a LOWESS smoothing. The variance was detrended by performing a Hilbert transform of the signal to get the instantaneous amplitude, calculating a LOWESS smoothing and by dividing the signal by the LOWESS smoothing. A Long-term linear detrending was performed using a LOWESS smoothing span of 1/1 and was applied to the record in order to investigate the low frequencies (Fig. 3). A LOWESS smoothing span of 1/15 was used on the record to investigate the high frequencies (Fig. 3). MTM spectrum analyses65 and Evolutive Harmonic Analyses (EHA66) were then ran on the interpolated series using a padding factor of 6 and a window size of 3,000 years. Confidence levels were estimated using a red noise fit (ML9667). Coherence Multi-Taper estimates analyses were performed between our data and literature palaeo-climatic/paleo-oceanographic time-series using cross-MTM method in Matlab (e.g.68). All time-series were linearly re-interpolated at 40 yrs from 2,488 to 10,200 yrs B.P. prior to the CMTM. The CMTM was estimated with 3 windows for the ASI-AF and ASI-NAC coherence estimations. Two windows were used for the ASI-Bond coherence estimate to better express the coherence at low frequencies (long time periods). All CMTM were run with 5000 Monte Carlo iterations for estimation of phase uncertainty.
Change history
06 August 2020
An amendment to this paper has been published and can be accessed via a link at the top of the paper.
References
Collins, M. et al. Chapter 12 - Long-termclimate change: Projections, commitments and irreversibility. In Climate Change 2013: The Physical Science Basis. IPCC Working Group I Contribution to AR5 (ed. IPCC, Cambridge University Press) (2013).
Hurrell, J. W. Decadal Trends in the North Atlantic Oscillation: Regional Temperatures and Precipitation. Science 269, 676–679 (1995).
Hurrell, J. W. & Deser, C. North Atlantic climate variability: The role of the North Atlantic Oscillation. Journal of Marine Systems 78, 28–41 (2009).
Barrier, N., Cassou, C., Deshayes, J. & Treguier, A.-M. Response of North Atlantic Ocean Circulation to Atmospheric Weather Regimes. J. Phys. Oceanogr. 44, 179–201 (2013).
Klus, A., Prange, M., Varma, V., Tremblay, L. B. & Schulz, M. Abrupt cold events in the North Atlantic in a transient Holocene simulation. Clim. Past Discuss, in review, https://doi.org/10.5194/cp-2017-106 (2017).
Thornalley, D. J. R., Elderfield, H. & McCave, I. N. Holocene oscillations in temperature and salinity of the surface subpolar North Atlantic. Nature 457, 711–714 (2009).
Bader, J. et al. A review on Northern. Hemisphere sea-ice, storminess and the North Atlantic Oscillation: Observations and projected changes, Atmospheric Research 101, 809–834 (2011).
Moffa-Sanchez, P., Born, A., Hall, I. R., Thornhalley, D. J. R. & Barker, S. Solar forcing of North Atlantic surface temperature and salinity over the past millennium. Nature geoscience 7, 275–278 (2014).
Comas‐Bru, L. & McDermott, F. Impacts of the EA and SCA patterns on the European twentieth century NAO–winter climate relationship. Q.J.R. Meteorol. Soc. 140, 354–363 (2014).
Sánchez-López, G. et al. Climate reconstruction for the last two millennia in central Iberia: The role of East Atlantic (EA), North Atlantic Oscillation (NAO) and their interplay over the Iberian Peninsula. Quaternary Science Reviews 149, 135–150 (2016).
Feser, F. et al. Storminess over the North Atlantic and North-western Europe—A review. Q.J.R. Meteorol. Soc 141, 350–382 (2015).
Goslin, J. & Clemmensen, L. B. Proxy records of Holocene storm events in coastal barrier systems: Storm-wave induced markers. Quaternary Science Reviews 174, 80–119 (2017).
Björck, S. & Clemmensen, L. B. Aeolian sediment in raised bog deposits, Halland, SW Sweden: a new proxy record of Holocene winter storminess variation in southern Scandinavia? The Holocene 14, 677–688 (2004).
Nielsen, P. R., Dahl, S. O. & Jansen, H. L. Mid- to late Holocene aeolian activity recorded in a coastal dunefield and lacustrine sediments on Andøya, northern Norway. The Holocene 26, 1486–1501 (2016).
Orme, L. C. et al. Past changes in the North Atlantic storm track driven by insolation and sea-ice forcing. Geology 45, 335–338 (2017).
Björck, S., Rundgren, M., Ljung, K., Unkel, I. & Wallin, Å. Multi-proxy analyses of a peat bog on Isla de los Estados, easternmost Tierra del Fuego: a unique record of the variable Southern Hemisphere Westerlies since the last deglaciation. Quaternary Science Reviews 42, 1–14 (2012).
DeVries-Zimmermann, S., Fisher, T. G., Hansen, E. C., Dean, S. & Björck, S. Sand in lakes and bogs in Allegan County, Michigan, as a proxy for eolian sand transport, Geological Society of America. Special Papers 508 (2014).
Jonassen, H. Bidrag til Filsøegnens naturhistorie. Meddelser fra Dansk Geologisk Forening 13, 192–205 (1957).
Clemmensen, L. B., Pedersen, K., Murray, A. & Heinemeier, J. A 7000-year record of coastal evolution, Vejers, SW Jutland, Denmark. Bulletin of the Geological Society of Denmark 53, 1–22 (2006).
Vink, A., Steffen, H., Reinhardt, L. & Kaufmann, G. Holocene relative sea-level change, isostatic subsidence and the radial viscosity structure of the mantle of northwest Europe (Belgium, the netherlands, Germany, southern North Sea). Quat. Sci. Rev. 26, 3249e3275 (2007).
Nielsen, S. T., Clemmensen, L. B. & Andreasen, F. The middle and late Holocene barrier spit system at Vejers, Denmark: Structure and development. Bull. Geol. Soc. Denmark 42, 105–119 (1995).
Fruergaard, M. et al. High-resolution reconstruction of a coastal barrier system: Impact of Holocene sea-level change. Sedimentology 62(3), 928–969 (2015).
Dalsgaard, K. & Odgaard, B. V. Dating sequences of buried horizons of podzols developed in wind-blown sand at Ulfborg, western Jutland. Quaternary 18(1), 53–60 (2001).
Zeeberg, J. The European sand-belt in Eastern Europe and comparison of Late Glacial dune orientation with GCM simulation results. Boreas. 27, 127–139 (1998).
Fruergaard, M. et al. Stratigraphy, evolution, and controls of a Holocene transgressive-regressive barrier island under changing sea-level: Danish North Sea coast. Journal of Sedimentary Research 85(7), 820–844 (2015).
Clemmensen, L. B., Murray, A., Heinmeier, J. & de Jong, R. The evolution of Holocene coastal dunefields, Jutland, Denmark: A record of climate change over the past 5000 years. Geomorphology 105, 303–313 (2009).
Clemmensen, L. B., Hansen, K. W. T. & Kroon, A. Storminess variation at Skagen, northern Denmark since AD 1860: Relations to climate change and implications for coastal dunes. Aeolian Research 15, 101–112 (2014).
Clemmensen, L. B., Andreasen, f, Nielsen, S. T. & Sten, E. The late Holocene coastal dunefield at Vejers, Denmark: characteristics, sand budget and depositional dynamics. Geomorphology 17, 63–77 (1996).
Hannon, G. E. & Giallard, M. J. The plant-macrofossil record of past lake-level changes. Journal of Paleolimnology 18, 15–28 (1997).
Sorrel, P. et al. Persistent non-solar forcing of Holocene storm dynamics in coastal sedimentary archives. Nature Geoscience 5, 892–896 (2012).
Bond, G. et al. Persistent Solar Influence on North Atlantic Climate During the Holocene. Science 294, 2130–2136 (2001).
Debret, M. et al. Evidence from wavelet analysis for a mid-Holocene transition in global climate forcing. Quaternary Science Reviews 28, 2675–2688 (2009).
Wanner, H., Solomina, O., Grosjean, M., Ritz, S. P. & Jetel, M. Structure and origin of Holocene cold events. Quaternary Science Reviews 30, 3109–3123 (2011).
Darby, D. A., Ortiz, J. D., Grosch, C. E. & Lund, S. P. 1,500-year cycle in the Arctic Oscillation identified in Holocene Arctic sea-ice drift. Nature Geoscience 5, 897–900 (2012).
Ojala, A. E. K., Launonen, I., Holmström, L. & Tiljander, M. Effects of solar forcing and North Atlantic oscillation on the climate of continental Scandinavia during the Holocene. Quaternary Science Reviews 112, 153–171 (2015).
Mayewski, P. A. et al. Holocene climate variability. Quaternary Research 62, 243–255 (2004).
Colin, C., Frank, N., Copard, K. & Douville, E. Neodymium isotopic composition of deep-sea corals from the NE Atlantic: implications for past hydrological changes during the Holocene. Quaternary Science Reviews 29, 2509–2517 (2010).
Renssen, H. et al. Simulating the Holocene climate evolution at northern high latitudes using a coupled atmosphere-sea ice-ocean-vegetation model. Clim Dyn 24, 23–43 (2005).
Debret, M. et al. The origin of the 1500-year climate cycles in Holocene North-Atlantic records. Clim. Past 3, 569–575 (2007).
Fletcher, W. J., Debret, M. & Goñi, M. F. S. Mid-Holocene emergence of a low-frequency millennial oscillation in western Mediterranean climate: Implications for past dynamics of the North Atlantic atmospheric westerlies. The Holocene 23, 153–166 (2013).
Wunsch, C. Abrupt climate change: An alternative view. Quaternary Research 65, 191–203 (2006).
Wanner, H. et al. Mid- to Late Holocene climate change: an overview. Quaternary Science Reviews 27, 1791–1828 (2008).
Sejrup, H. P., Haflidason, H. & Andrews, J. T. A Holocene North Atlantic SST record and regional climate variability. Quaternary Science Reviews 30, 3181–3195 (2011).
Knudsen, M. F., Jacobsen, B. H., Riisager, P., Olsen, J. & Seidenkrantz, M.-S. Evidence of Suess solar-cycle bursts in subtropical Holocene speleothem? 18O records. The Holocene 22, 597–602 (2012).
Knudsen, M. F., Seidenkrantz, M.-S., Jacobsen, B. H. & Kuijpers, A. Tracking the Atlantic Multidecadal Oscillation through the last 8,000 years. Nature Communications 2, 178 (2011).
Turney, C. S. M. et al. A 250-year periodicity in Southern Hemisphere westerly winds over the last 2600 years. Clim. Past 12, 189–200 (2016).
Knudsen, M. F. et al. Taking the pulse of the Sun during the Holocene by joint analysis of 14C and 10Be. Geophys. Res. Lett. 36, L16701 (2009).
Lockwood, M. et al. Top-down solar modulation of climate: evidence for centennial-scale change. Environ. Res. Lett. 5, 34008 (2010).
Ineson, S. et al. Solar forcing of winter climate variability in the Northern Hemisphere. Nature Geoscience 4, 753–757 (2011).
Thiéblemont, R., Matthes, K., Omrani, N.-E., Kodera, K. & Hansen, F. Solar forcing synchronizes decadal North Atlantic climate variability. Nature Communications 6, 8268 (2015).
Stuiver M. & Brazuinas T. Modeling atmospheric 14C influences and 14C ages of marine samples to 10,000 bc. In Stuiver, M., Long, A. & Kra, R. S. eds, Calibration, Radiocarbon 35(1), 137–189 (1993)
Meehl, G. A., Arblaster, J. M., Matthes, K., Sassi, F. & Loon, Hvan Amplifying the Pacific Climate System Response to a Small 11-Year Solar Cycle Forcing. Science 325, 1114–1118 (2009).
Repschläger, J., Garbe-Schönberg, D., Weinelt, M. & Schneider, R. Holocene evolution of the North Atlantic subsurface transport. Clim. Past 13, 333–344 (2017).
Bjune, A. E., Bakke, J., Nesje, A. & Birks, H. J. B. Holocene mean July temperature and winter precipitation in western Norvay inferred from palynological and glaciological lake-sediment proxies. The Holocene 15, 177–189 (2005).
Taylor, A. K., Benedetti, M. M., Haws, J. A. & Lane, C. S. Mid-Holocene Iberian hydroclimate variability and paleoenvironmental change: molecular and isotopic insights from Praia Rei Cortiço, Portugal. J. Quaternary Sci. 33, 79–92 (2018).
Kobashi, T. et al. Volcanic influence on centennial to millennial Holocene Greenland temperature change. Scientific Reports 7, 1441 (2017).
Giraudeau, J. et al. Millennial-scale variability in Atlantic water advection to the Nordic Seas derived from Holocene coccolith concentration records. Quaternary Science Reviews 29, 1276–1287 (2010).
Andersen, C., Koç, N., Jennings, A. & Andrews, J. T. Nonuniform response of the major surface currents in the Nordic Seas to insolation forcing: Implications for the Holocene climate variability. Paleoceanography 19, PA2003 (2004).
Moros, M., Andrews, J. T., Eberl, D. D. & Jansen, E. Holocene history of drift ice in the northern North Atlantic: Evidence for different spatial and temporal modes. Paleoceanography 21, PA2017 (2006).
Reimer, P. J. et al. IntCal13 and Marine13 Radiocarbon Age Calibration Curves 0–50,000 Years cal BP. Radiocarbon 55, 1869–1887 (2013).
Blaauw, M. Methods and code for ‘classical’ age-modelling of radiocarbon sequences. Quaternary Geochronology 5, 512–518 (2010).
Pye, K. & Tsoar, H. Aeolian sand and sand dunes. (Springer Netherlands, 1990).
Kylander, M. E., Ampel, L., Wohlfarth, B. & Veres, D. High-resolution X-ray fluorescence core scanning analysis of Les Echets (France) sedimentary sequence: new insights from chemical proxies. Journal of Quaternary Science 26, 109–117 (2011).
Meyers, S. R. Astrochron: An R Package for Astrochronology. https://cran.r-project.org/package=astrochron (2014).
Thomson, D. J. Quadratic-inverse spectrum estimates: applications to palaeoclimatology. Phil. Trans. R. Soc. Lond. A 332, 539–597 (1990).
Kodama, K. P. & Hinnov, L. A. Time Series Analysis for Cyclostratigraphy. in Rock Magnetic Cyclostratigraphy 52–89 (John Wiley & Sons, Ltd, 2014).
Mann, M. E. & Lees, J. M. Robust estimation of background noise and signal detection in climatic time series. Climatic Change 33, 409–445 (1996).
Huybers, P. & Denton, G. Antarctic temperature at orbital timescales controlled by local summer duration. Nature Geoscience 1, 787–792 (2008).
Acknowledgements
This paper benefited from a financial support from two Carlsberg Foundation post-doctoral fellowships granted to J. Goslin (CF-16-0349) and M.Fruergaard (CF15-0254). J.Monkenbusch is supported by Carlsberg project CF-16-0457. L.Menviel acknowledges funding from the Australian Research Council grants DE150100107 and DP180100048. We express our gratitude to the Carlsberg Foundation for their continuous support throughout this study. We thank A. Murray at the Nordic Laboratory for Luminescence Dating (Risø, Denmark) for conducting the OSL dating. J.Goslin would like to express his gratitude to P. Friberg and P. Christiansen (both Institute from the department of Natural Resources Management IGN, University of Copenhagen, CPH) for their invaluable help and advice before and during field campaigns as well as A. Thoisen (IGN, Uni CPH) for her help and advices during laboratory analyses. M.L. Siggaard-Andersen (Natural History Museum of Copenhagen) is thanked for her help in running the XRF core-scanner and for post-processing the data. P. Bernatchez (Université du Québec à Rimouski) is thanked for having facilitated the revision of the manuscript through a post-doctoral fellowship granted to J.Goslin. The Geological Survey of Denmark and Greenland (Copenhagen, Denmark) is acknowledged for providing access to their cold-storage and core photography facilities. N.Marriner and two anonymous reviewers are thanked for their suggestions that helped to improve the manuscript. Finally, we are grateful to M. Weinelt (Univ. Kiel, Germany), T. Kobashi (Univ. of Bern, Switzerland) and J. Giraudeau (UMR EPOC, Univ. Bordeaux, France) for showing interest in our study and providing data and contacts.
Author information
Authors and Affiliations
Contributions
J.G. designed and led the project, collected and analyzed the sedimentary cores, interpreted the data, drew the figures and wrote the manuscript. M.F. and L.S. assisted collecting the cores, helped in the description and in the sampling of the cores, co-interpreted the data and co-wrote the manuscript. M.G. analysed plant macrofossils for radiocarbon dating and plant macro-remains analyses and provided corrections to the manuscript. L.M. provided LOVECLIM model simulations and provided corrections to the manuscript. J.M. provided his expertise in using the Astrochron package for spectral and wavelet analyses and helped tuning the analyses and provided corrections to the manuscript. N.T. performed coherence analyses and provided corrections to the manuscript. L.B.C. initiated and co-designed the project, provided overall scientific guidance, co-interpreted the data and co-wrote the manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare 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
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
Goslin, J., Fruergaard, M., Sander, L. et al. Holocene centennial to millennial shifts in North-Atlantic storminess and ocean dynamics. Sci Rep 8, 12778 (2018). https://doi.org/10.1038/s41598-018-29949-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-018-29949-8
This article is cited by
-
Drought as a possible contributor to the Visigothic Kingdom crisis and Islamic expansion in the Iberian Peninsula
Nature Communications (2023)
-
Past 200 kyr hydroclimate variability in the western Mediterranean and its connection to the African Humid Periods
Scientific Reports (2022)
-
Striking forest revival at the end of the Roman Period in north-western Europe
Scientific Reports (2020)
-
Palaeo-productivity record from Norwegian Sea enables North Atlantic Oscillation (NAO) reconstruction for the last 8000 years
npj Climate and Atmospheric Science (2020)
-
Historical predictability of rainfall erosivity: a reconstruction for monitoring extremes over Northern Italy (1500–2019)
npj Climate and Atmospheric Science (2020)
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.