1 A composite annual-resolution stalagmite record of North Atlantic climate over the last three millenniaAndy Baker1, John C. Hellstrom2, Bryce F. J. Kelly1, Gregoire Mariethoz3 & Valerie Trouet4Annually laminated stalagmites can be used to construct a precise chronology, and variations in laminae thickness provide an annual growth-rate record that can be used as a proxy for past climate where precipitation is sensitive to North Atlantic climate variability and the winter North Atlantic NAO states, during the Medieval Climate Anomaly (MCA). Another persistently low growth period occurring at 290-550 CE coincides with the European Migration Period, and a subsequent period of sustained fast growth-rate (negative NAO) from 600-900 AD provides the climate context for the Viking Age in northern and western Europe.Stalagmites contain a variety of potential paloclimate proxies, and the sensitivity of any individual stalagmite proxy to climate is determined by many factors, including its location within the cave, the cave morphology, the hydrological ow path, and the overlying soil and vegetation1. In NW Scotland, Uamh an Tartair (Roaring Cave) is overlain by blanket peat, which accumulated over the last ~4000 years2, and has a rainfall regime that is highly correlated with the winter NAO3. Because of the caves shallow depth, an overlying soil carbon and water store, and a blanket peatland that experiences substantial seasonal water table variations, stalagmites within Uamh an Tartair contain annual uorescent laminae4 that are continuous for 103-104 yrs5,6. ese annual laminae provide both an annually resolved chronology and a climate proxy, with the dominant control on annual growth rate being the water table in the overlying peat, which aects soil CO2 production and concentration, limestone dissolution, and subsequent stalagmite calcite precipitation. As a result, high (low) lamina growth rates typically reect warm and/or dry (cool and moist) conditions5. Using a combination of lamina counting and Udating, stalagmite growth rate records from a region can be combined to produce long continuous paleoclimate records. One stalagmite from the site has been previously combined with tree-ring data in a 947-year annual resolution reconstruction of the NAO (NAOms), which suggested an unusual persistence of the positive phase of the NAO during the MCA (1050-1400 CE)7. Here, we produce a composite annual growth-rate record based on ve stalagmites collected from the Uamh an Tartair cave. Four of the stalagmites have previously been reported6,7,8 and we add a new addition (SU963) to the series.1Connected Waters Initiative Research Centre, UNSW Australia, Sydney, NSW 2052, Australia. 2School of Earth Sciences, The University of Melbourne, Melbourne, VIC 3010, Australia. 3University of Lausanne, Institute of Earth Surface Dynamics (IDYST), UNIL-Mouline, Geopolis, 1015 Lausanne, Switzerland. 4Laboratory of Tree-Ring Research, University of Arizona, Tucson AZ 85721, USA. Correspondence and requests for materials should be addressed to A.B. (email: email@example.com)Received: 22 January 2015 A ccepted: 08 April 2015 P ublished: 11 June 2015OPEN
2 ResultsBuilding a stalagmite growth-rate composite record. Stalagmite SU963 grew on the same owstone shelf as SU961 and SU962, and contained a continuous lamina sequence for 1802 years with no hiatuses. e growth-rate was determined from multiple manual analyses of images collected using UV microscopy (an example is provided in Supplemental Figure S1). e annual nature of the laminae was conrmed9 using Uanalyses (see Supplementary Table 1 and Methods). Five Uanalyses constrained the timing of deposition of SU963, between ~1600 and ~400 years before present. Using a combination of the Uages and cross correlation of annual lamina thickness with coeval stalagmites SU961, SU962, and SU967, the annual lamina sequence of 1802 laminae in SU963 was constrained to the period from 22 BCE to 1781 CE (see Methods). is produced a composite stalagmite record of annual growth-rate that has overlapping stalagmites for the period 1475 BCE-1996 CE. Limited replication of the record occurs in the early part of the record, and therefore we focus on the last 3000-yr composite growth-rate record (SUcomp). For the instrumental period, SUcomp has the strongest annual correlation with the winter NAO. A positive winter NAO is characterised by increased precipitation in NW Europe, including NW Scotland, and drier conditions in the western Mediterranean. is pattern is observed in the spatial correlation between SUcomp and winter (January-March) precipitation, with correlation p 5% using CRU TS3.2210 (Fig.1a). Spatial correlations are weak with winter (January-March) temperature (Fig.1b). Strong spatial correlations, with p 5%, are observed with winter (January-March) sea level pressure (Fig.1c). is observed pattern is typical of the strong correlations of opposite sign in the two NAO centers of action (the Icelandic Low and the Azores High). Final evidence of a NAO signal contained within the SUcomp series is the temporal correlation with the winter (January-March) NAO index11 (Fig.1d, Stykkishlmur-Azores, 1901-2004; r 0.46, p 0.01). e error on this correlation coecient is determined to be normal with a standard deviation of 0.08, based on a bootstrap analysis12, (see Methods). e winter season NAO signal therefore has the predominant correlation with SUcomp (climate correlations with other seasons are presented in Supplemental Figure S2). Despite this, an annual correlation of 0.46 indicates substantial noise will be present in this winter NAO signal. We suggest that this is due to several factors that include: (1) an observed weak, but signicant, annual correlation between SUcomp and spring temperature (See Supplemental Figure S2), which is unrelated to the NAO. (2) Mixing of water in the karst unsaturated zone which smooths the Figure 1. Pearsons correlation map between SUcomp and average elds of Jan-Mar precipitation (a), temperature (b), and sea level pressure (SLP); c). Precipitation and temperature time series were derived from the CRU TS3.22 data set, SLP from HadSLP2r. All correlations were calculated over the period 19012004 and climate variable distributions are signicant at the 95% condence level. e SUcomp (inversed) and Jan-Mar NAO (Stykkishomur-Azores index) time series (1901-2004; r 0.46, p 0.01) are plotted in (d). Correlation maps were generated using Climate Explorer32.
3 annual climate signal. (3) Variability of the sensitivity of individual stalagmite growth-rate series to the winter NAO, which could be better constrained with a larger stalagmite composite record. e 3000-year growth-rate record for the ve stalagmites is presented in Fig.2a. e availability of a composite stalagmite growth-rate record permits the screening of the dataset for climate complacency and non-climatic eects on stalagmite growth-rate. Many techniques have been developed for annual proxy series such as tree rings13, but methodologies appropriate to stalagmites including variogram analysis and growth acceleration14 are required here. Climatically sensitive stalagmites have growth-rates that icker around a mean state representative of the volume and chemistry of the karst water store supplying the stalagmite. An absence of ickering growth in SU967 at 1615-1645 CE demonstrated the dominance of non-linear karst hydrological responses (e.g. fracture storage being full and not climate sensitive) over this time period. Periods of no ickering were also identied in SU963 at 1428-1448 CE. ese periods have been removed from the composite series. Similar hydrological eects have been previously noted7 when stalagmites die due to a lack of stored water supply, with monotonous decreases or increases in growth-rate. With replicate samples available, these eects are clearly visible at the end of each deposition phase in stalagmites SU961 and SU962, and where observed in these specimens (the last 14 years of deposition) the data are removed from the lamina thickness record. Finally, variogram analysis of the information content of the stalagmite growth-rate series suggests that samples all had high information content (SU961, SU962, SU963 and SU967 have information contents of 75, 87, 93 and 70%, respectively, and variogram ranges values of 250, 200, 400 and 60 years, respectively). e new addition to SUcomp, SU963, has the highest information content and longest range of all samples in the composite, indicating that it is more likely to preserve multi-decadal variability in lamina thickness. e fully screened SUcomp is in Fig.2b (and available as Supplementary Data). In general, where multiple stalagmites are deposited at the same time, their growth-rates are similar. Despite statistical screening of the data, there are some time periods where one stalagmite has a diver gent growth rate. SU963 has faster growth-rate than other samples during periods of below average growth, and SU961 has a slower growth-rate during periods of above average growth. Between sample growth-rate variability should be anticipated due to known non-linear karst processes such as karst hydrology through changes in ow routing14 and the non-linear relationship between drip-rate and growth-rate16. Replicated, composite series such as SUcomp therefore provide the best approach to constraining speleothem growth-rate records with appropriate uncertainty bounds.Relating growth-rate persistence and periodicity to climate. Stalagmite SU967 has been previously utilised in the NAOms reconstruction, where a persistently positive NAO was identied in the MCA1. SUcomp (Fig.2b) replicates the SU967 record over the past 1000 years, and conrms the persistently low Figure 2. (a) Five standardised stalagmite growth-rate series. (b) Screened SUcomp. Data are shown as mean (black) and 1 variability (grey). Periods where there are low growth runs are shown by the vertical shading.
4 growth-rates reective of cool and moist conditions during the MCA. To investigate whether continuous low stalagmite growth during this period has occurred by chance, we analysed the annual growth-rate distribution of each stalagmite. Figure3a highlights periods when annual growth-rate was persistently below the 10th percentile. Our analysis conrms that persistently low growth is observed for multiple stalagmite samples during the MCA and it is also observed at 250-550 CE. We performed wavelet analysis on the growth-rate series for all stalagmites. e global wavelet spectra (Supplementary Figure 3) demonstrate signicant spectral power at periods greater than 40 years for all stalagmites. Wavelet analysis demonstrates a non-stationarity of spectral power over time, whose magnitude (Fig.3b) is consistent between replicate stalagmite records. Specically, multi-decadal to centennial-scale periodicity is present throughout all stalagmite records, but multi-decadal-scale periodicity decreases in magnitude from ~1100 to ~1450 CE (SU967 and SU963) and from ~0 to ~500 CE (SU961, SU962, and SU963). Other records from the North Atlantic region also show multi-decadal variability17,18, attributed to the Atlantic Multidecadal Oscillation (AMO). Comparison of the timing of persistent slow stalagmite growth (1080-1460 CE and 290-550 CE) and that of low magnitude of decadal to centennial-scale periodicity (~1100 to ~1450 CE and ~0 to ~500 CE) demonstrates that the two are coincident, with the opposite relationship in fast growth periods. Over the past 450 years, evidence has been presented that the AMO is externally forced and is probably related to the Atlantic Meridional Overturning Circulation (AMOC) strength, and also that this forcing is non-stationary over time19. If the decadal-scale periodicity observed in our composite stalagmite record is a response to the AMO, then a weakened AMO may have led to persistently positive NAO conditions not just aer the Little Ice Age19, but also at 1080-1430 CE (the Medieval Climate Anomaly) and at 290-550 CE (the Dark Ages).Comparison with NAO proxy records. Given the previously published evidence for the winter NAO sensitivity of annual growth-rates at this site and the strong correlations between SU and winter NAO over the instrumental period (Fig.1), we interpret SUcomp as one of North Atlantic climate and winter NAO. SUcomp shows good agreement with other published NAO-sensitive records (Fig.4), including winter drought-sensitive tree-ring20 (1049-1995 CE; r 0.40; n 926; p 0.001 at annual scale and r 0.56; n 31; p 0.01 at tri-decadal scale) and stalagmite records21 (748-1953 CE; r 0.47; n 197; p 0.001; 5-yr resolution) from Morocco and a West Greenland lake sediment record17 (985 BCE-1634 CE; r 0.34, n 167, p 0.001; 20-yr resolution). Comparison of the correlations between all records for a common period (1094-1634 CE) is presented in Supplemental Table 2. e SUcomp, Moroccan tree-ring, and Greenland lake sediment archives show a persistently positive winter NAO during the MCA (1080-1430 CE) that is unique within the last 1000 years (although we note that the Greenland lake sediment chronology was tuned to the NAOms series7, and would therefore also be expected to correlate with the SUcomp reconstruction from 1100 CE onwards). ese three Figure 3. (a) Distribution of stalagmite continuous low-growth periods. Low-growth periods are dened as a growth rate below the 10th percentile value of a given stalagmite. (b) Wavelet magnitude at periodicities of 64 (blue), 76 years (red) and 90 (black) years (red). Wavelet analysis was performed on all ve stalagmites using a Morlet transform. Shading shows the two periods of persistently slow stalagmite growth and positive NAO.
5 Figure 4. Comparison (991 BCE 2004 CE) of SUcomp (this study) with other proxy records of past NAO variability: a drought-sensitive tree-ring record from Morocco20, a drought-sensitive Sr stalagmite record from Morocco21, and a West Greenland ice core record17. All records were normalised over their full period of record and all records but the ice core record were inverted so that positive values reect positive NAO conditions. e Moroccan stalagmite and ice-core records are not continuous and these series were smoothed with a cubic spline that closely reects their average sampling interval (5 years for stalagmite data, 20 years for ice-core data). SUcomp was smoothed with a 20-year cubic spline (red and dark blue colors) to match the average resolution of the ice-core record. e most recent portion of SUcomp (orange and light blue colors) as well as the tree-ring record were smoothed with a 5-year cubic spline to match the average resolution of the Moroccan stalagmite record. e transition period (14301500) between a positive NAO phase during the MCA and a negative NAO phase during the LIA is indicated by a yellow bar.
6 records do not observe a more variable NAO state in the MCA as observed in the trace element record of a Moroccan stalagmite21. e latter NAO reconstruction was based on Mg/Ca and Sr/Ca ratios in a stalagmite, and we note that these trace element proxies are poorly understood in semi-arid regions where wind-blown dust can complicate signal calibration22. A distinct shi from this positive NAO phase into a negative phase during the Little Ice Age occurs simultaneously in all records at ~1430 CE (Fig.4). e timing of this climate dynamical MCA-LIA transition thus appears synchronous across the North Atlantic realm, although temporal osets have been recorded in dierent sub-continental regions23. Prior to the last 1000 years, we observe periods of good agreement in the reconstructed NAO between the two available records, the West Greenland lake sediment record and SUcomp, for example between 200 BCE 500 CE. However, there are also periods of dierence, especially before 500 BCE. Reasons for the dierences could include dating uncertainties in both series (the lake record is constrained by radiocarbon analyses and has a mean age model uncertainty of 54 years) or, more likely, limitations in the strength NAO proxy calibration in the lake sediment record. ere, the NAO-like pattern is found in the third principal component, which explains only 9.4% of the variability in the data. Only further long, high-resolution proxy records can help to develop improved NAO reconstructions.DiscussionOur 3000-year record conrms persistently low growth-rates, reective of positive NAO states, dur ing the MCA. Another persistently low growth period occurred at 290-550 CE, and wavelet analysis demonstrates that period of slow (or fast) growth are coincident with periods of relatively low (or high) magnitude of multi-decadal variability, indicating a potential relationship between the strength of the Atlantic Multidecadal Oscillation and the NAO at these times. e phase of persistent positive NAO that occurs at ~290-550 CE coincides with a period of increased climate variability and an unusually wet growing season (April-June) in central Europe as observed in tree rings24. We hypothesise that composite stalagmite and tree ring proxies are responsive to the same North Atlantic climate forcing, which leads to a persistently positive NAO. Assuming that we can utilise the SUcomp record as a proxy for the NAO over the last 3000 years, this allows us to start building hypotheses relating to the impact of climate variability on environmental and societal change in NW Europe. For example, although we recognise that climatic changes alone are unlikely to be the sole driver for societal change25, it is noted that the positive NAO phase at ~290550 CE is coincident with the decline of the Roman Empire and the European Migration period. is was followed by generally negative NAO conditions in the period 600-900 AD, which would provide suitable conditions for westward expansion in the north Atlantic region. In NW Europe, the timing coincides with the start of the Viking Age of expansion into northern and western Europe. e cause and precise timing of this is contested26, and our SUcomp record provides a useful climatic context of colder Scandinavian climate and weaker westerly circulation for this important phase of early European history.Lamina Counting. e uorescent laminae for stalagmite SU963 were counted and their thicknesses measured in duplicate using techniques detailed previously5. No hiatuses were observed and a total of 1704 continuous laminae were counted, with a mean lamina thickness of 27.0 m and a mean standard deviation between replicate measurements of 5.2 m (n 1704). e mean lamina thickness for all stalagmites reported in this paper is archived at the World Data Centre Paleoclimatology.Stalagmite SU963 was dated with Umass spectrometry using methods previously reported28. Uages and isotope data are presented in Supplementary Table 1. e Uages were used to constrain the oating lamina chronology to the Uages. Following established methods29, the relative lamina number was subtracted from each Uage to obtain ve determinations of the date of the end of deposition. is approach yielded a date for the end of deposition of 220 years before present (AD2010) (95% Condence Interval: 143-315 years before present). We compared this date with that calculated from a linear regression of the Uage vs the relative lamina number (Uage 288.5 0.9088 x lamina number; r 0.99, n 5). e slope of the regression (0.91, with a 95% condence range of 0.71 to 1.10) is indistinguishable from zero, provided independent conrmation of the continuous, annual nature of lamina deposition. e intercept provided a date for the end of deposition of 289 years before present (95% Condence Interval: 124-454 years before present). e Ubased and lamina-based age-depth relationships are presented in Supplemental Figure S4.Correlation of growth-rate series. To further constrain the timing of deposition of SU963, the annual lamina thickness series was cross-correlated with that of stalagmites SU961, SU962 and SU967. Cross-correlations were calculated for theoretical end-of-deposition dates between 120 and 320 years before present. e strongest cross-correlations were observed when the period of divergent lamina thickness in SU963 of 150 years at ~1550 years before present was removed from the analysis: for SU961 with an end-of-deposition date of 210 year before present (r 0.25); for SU962 with an end-of-deposition date of 229 years before present (r 0.46); and for SU967 with an end-of-deposition date of 223 years before present (r 0.34). e mean of these three approaches yields a date of growth cessation of 221 years before present. We have placed the SU963 lamina thickness series into the chronology with the end
7 of deposition at 229 years before present, which is identical to the mean of the Upredicted date of the end of deposition and with the strongest cross-correlating stalagmite (SU962).Geostatistical methods. Bootstrapping of instrumental and composite series, variogram analysis, growth acceleration and wavelet analyses of the annual growth-rate series were performed in MATLAB, Wolfram Mathematica and Microcal Origin, and wavelet analyses visualised using the IDL Wavelet Toolkit31. Bootstrapping was carried out by generating 100,000 replicates of the data by sampling with replacement years of the instrumental period, each set of dates yielding a new time series for both NAO and SUcomp. Correlation analyses were performed using Microcal Origin, SPSS, Excel and Climate Explorer32. To allow correlation analysis between proxy records of dierent temporal resolution, we smoothed the annual resolution SUcomp record using a cubic smoothing spline to the approximate resolution of the lower resolution proxies (5 years for the Morocco stalagmite21, 20 years for the West Greenland lake sediment record17, and 30 years for NAOms 1 and the Morocco tree-ring record20). To reduce the eects of autocorrelation, we sampled the smoothed time series every 30 years for NAOms and the Morocco tree-ring record20 or at the years when data points were recorded in the other proxies17,21. We then calculated Pearsons correlation coecients for time series consisting only of these sequential samples.References1. Fairchild, I. J. & Baer, A. Speleothem Science. Wiley, London (2012). 2. Charman, D. J., Caseldine, C. J., Baer, A., Gearey, B. & Hatton, J. Palaeohydrological records from peat proles and speleothems in Sutherland, NW Scotland. Quat. es 55, 223 (2001). 3. Hurrell, J. W. Decadal trends in the North Atlantic Oscillation: regional temperatures and precipitation. Science 269, 676 (1995). 4. Baer, A., Smart, P. L., Edwards, L. & ichards, D. A. Annual banding in a cave stalagmite. Nature 364, 518 (1993). 5. Proctor, C. J., Baer, A., Barnes, W. L. & Gilmour, M. A. A thousand year speleothem proxy record of North Atlantic climate from Scotland. Clim. Dyn. 16, 815820 (2000). 6. Proctor, C. J., Baer, A. & Barnes, W. L. A three thousand year record of north Atlantic climate. Clim. Dyn. 19, 449 (2002). 7. Trouet, V. et al. Persistent positive North Atlantic Oscillation mode dominated the Medieval Climate Anomaly. Science 324, 78 (2009). 8. Baer, A., et al. High resolution 18O and 13C records from an annually laminated Scottish stalagmite and relationship with last millennium climate. Glob. Planet. Chang 79, 303 (2011). 9. Shen, C-C. et al. Testing the annual nature of speleothem banding. Sci. ep. 3, 2633 (2013). 10. Harris, I., Jones, P. D., Osborn, T. J. & Lister, D. H. Updated high-resolution grids of monthly climatic observations the CU TS3.10 Dataset. Int. J. Climatol. 34, 623 (2013). 11. alnay, E. & coauthors, e NCEP/NCA eanalysis 40-year Project. B. Am. Meteorol. Soc. 77, 437 (1996). 12. Efron, B. Bootstrap Methods: Another Loo at the Jacnife. Ann. Stat. 7, 1 (1979). 13. Esper, J., Coo, E. ., rusic, P. J., Peters, & Schweingruber, F. H. Tests of the CS method for preserving low-frequency variability in long tree-ring chronologies. Tree-ing es. 59, 81 (2003). 14. Mariethoz, G., elly, B. & Baer, A. Quantifying the value of laminated stalagmites for paleoclimate reconstructions. Geophys. es. Lett. 39, L05407 (2012). 15. Baer, A., Bradley, C. & Phipps, S. J. Hydrological modelling of stalagmite 18O response to glacial-interglacial transitions. Geophys. es. Lett. 40, 1 (2013). 16. Dreybrodt, W. Chemical inetics, speleothem growth and climate. Boreas 28, 347(1999). 17. Olsen, J., Anderson, N. J. & nudsen, M. F. Variability of the North Atlantic Oscillation over the past 5,200 years. Nat. Geosci. 5, 808 (2012). 18. lafsdttir, ., Geirsdottir, A., Miller, G. H., & Larsen, D. J. Evolution of NAO and AMO strength and cyclicity from a 3-a varve-thicness record from Iceland. Quat. Sci. ev. 69, 142 (2013). 19. nudsen, M. F., Jacobsen, B. H., Seidenrantz, M-S. & Olsen, J. Evidence for external forcing of the Atlantic Multidecadal Oscillation since termination of the Little Ice Age. Nat. Commun. 5, 3323 (2014). 20. Esper, J., Fran, D., Buntgen, U., Verstege, A. & Luterbacher, J. Long-term drought severity variations in Morocco. Geophys. es. Lett. 34, L17702 (2007). 21. Wassenburg, J. A. et al. Moroccan Speleothem and tree ring records suggest a variable positive state of the North Atlantic Oscillation during the Medieval Warm Period. Earth Planet. Sci. Lett. 375, 291 (2013). 22. utlidge, H. et al. Dripwater organic matter and trace element geochemistry in a semi-arid arst environment: implications for speleothem paleoclimatology. Geochim. Cosmochim. Acta 135, 217 (2014). 23. Ahmed, M. et al. Continental-scale temperature variability during the past two millennia. Nat. Geosci. 6, 339 (2013). 24. Bntgen, U et al. 2500 Years of European climate variability and human susceptibility. Science 331, 578 (2011). 25. Diaz, H. & Trouet, V. Some perspectives on societal impacts of past climatic changes. History Compass 12, 160 (2014). 26. Barrett, J. H. What caused the Viing Age? Antiquity 82, 671 (2008). 27. von Holstein I. C. C. et al. Searching for Scandinadians in pre-Viing Scotland: molecular ngerprinting of Early Medieval combs. J. Archaeol. Sci. 41, 1 (2014) 28. Hellstrom, J. apid and accurate U/ dating using parallel ion-counting multi-collector ICP-MS. J. Anal. At. Spectrom. 18, 1346 (2003). 29. Asrat, A., et al. A high-resolution multi-proxy stalagmite record from Mechara, Southeastern Ethiopia: palaeohydrological implications for speleothem palaeoclimate reconstruction. J. Quat. Sci. 22, 53 (2007). 30. Dominguez-Villar, D., Baer, A., Fairchild, I. J. & Edwards, L. A method to anchor oating chronologies in annually laminated speleothems with Udates. Quat. Geochronol. 14, 57 (2012). 31. Torrence, C. & Compo, G. P. A practical guide to Wavelet Analysis. B. Am. Meteorol. Soc. 79, 61 (1998). 32. Trouet, V. & Van Oldenborgh, G. J. NMI Climate Explorer: A web-based research tool for high-resolution paleoclimatology. Tree-ing es. 69, 3 (2013).
8 Author ContributionsLamina measurements were carried out by A.B., and Uanalyses were undertaken by J.H. Geostatistical analyses were made by G.M., B.K. and A.B. Proxy and instrumental climate comparisons were undertaken by V.T. and A.B. wrote the manuscript with input from all authors.Additional InformationSupplementary information accompanies this paper at http://www.nature.com/srep Competing nancial interests: e authors declare no competing nancial interests. How to cite this article: Baker, A. et al A composite annual-resolution stalagmite record of North Atlantic climate over the last three millennia. Sci. Rep. 5, 10307; doi: 10.1038/srep10307 (2015). is work is licensed under a Creative Commons Attribution 4.0 International License. e images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/