the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
High-Resolution Paleo-Storm Reconstruction from Eastern Canada Aligns with Late-Holocene Northwestern Atlantic Hurricane Records
Abstract. Atlantic Canada experiences frequent major storms, particularly tropical cyclones transitioning into post-tropical storms. Events such as Hurricane Fiona (2022), Dorian (2019), and Juan (2003) have caused significant damage, loss of life, and coastal erosion, exacerbated by sea level rise and warming waters. Despite this, centennial- to millennial-scale storm records in the region remain scarce. Existing studies in North America focus primarily on marine and coastal overwash records, with limited use of aeolian mineral inputs in ombrotrophic peatlands as storm proxies. Here, we address these gaps by analysing grain-size and geochemical data from two peatlands in Quebec, Canada’s Magdalen Islands.
Our two peat records reveal consistent storm signals over the past 4000 years, with three key periods of heightened activity: 800–550 BCE, 600–800 BCE, and 1300–1700 CE. These signals align with marine and overwash records spanning the past 2000 years across eastern Canada, the US, and the Bahamas, indicating low storm activity during the Medieval Climate Anomaly, followed by increased activity during the Little Ice Age. Our findings suggest that storm records in these regions are influenced by local climatic factors. Negative phases of the Atlantic Multidecadal Variability, which typically suppress hurricane activity in the North Atlantic, are associated with conducive hurricane formation and intensification north of the Bahamas. Additionally, the position of the Bermuda High seems to play a more significant role in directing storm tracks during different climatic phases. Our findings highlight the potential antiphase relationship in storm activity between regions north of the Bahamas and those in the Gulf of Mexico, suggesting broader climatic mechanisms that warrant further investigations.
Despite the similarities between our two sites, discrepancies in geochemistry and mineralogical profiles highlight the importance of site-specific conditions in interpreting the storm record from peatlands, namely the distance of the sites to the coast and source of aeolian sediment, as well as peatland size. Challenges also remain in calibrating peat-based proxies with historical storm records, as identifying specific events from the past 150 years remains difficult.
Competing interests: At least one of the (co-)authors is a member of the editorial board of Climate of the Past. Dr. Pierre Francus is a member of the editorial board of Climate of the Past.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this preprint. The responsibility to include appropriate place names lies with the authors.- Preprint
(8600 KB) - Metadata XML
- BibTeX
- EndNote
Status: open (until 16 Apr 2025)
-
RC1: 'Comment on egusphere-2025-400', Anonymous Referee #1, 10 Mar 2025
reply
General Comments
The manuscript by Lachance and coauthors presents two new records of past storm activity from northeastern Canada. Methods of Aeolian Sediment Influx (using a sieving method) and ITRAX XRF analysis are used to measure the past aeolian sand deposits to the two bogs, spanning the late Holocene. The manuscript is very well written, with clear language, well presented figures and appropriate numbers of citations in most parts. The statistical analysis of the data appears to me (a non expert) to be appropriate, although as detailed below I have questions about the use of unnormalized itrax count data in parts of the paper and the chosen method of detecting storm events. The discussion is particularly interesting and well written, and includes comparisons between the storm records and individual recent storms, as well as a very good comparison of records from a latitudinal transect. The conclusion includes interesting take home messages, but could be strengthened by detailing more of the results and findings.
Major Corrections
1) Unnormalized ITRAX count data: the authors do a good job with the statistical analysis of the ITRAX data, converting it to Centered Log Ratios before doing the statistical analysis, because ITRAX data is compositional data. Therefore there is an acknowledgement that changes in one element can affect the counts of another element, and hence the need to normalize it. It was therefore surprising that following the statistical analysis the Ti records that are shown and on which the storm records are based were raw counts per second (kcps) results, rather than the CLR results. Previous studies (e.g. Löwemark et al. 2011; Weltje and Tjallingii, 2008; Croudace et al. 2019) have highlighted that ITRAX results should be normalized because otherwise variations in the amount of undetected organic matter, water content, surface roughness and counts of other elements can all influence the results. Appendix D shows photos of the cores, which have a very rough surface – was the surface smoothed at all prior to ITRAX scan as this would surely influence the ITRAX results greatly? This makes it important that the results are normalized as the variability in the surface would have likely influenced the measured Ti kcps. The CLR results (or natural log ratios with a carefully selected denominator) should provide the basis for the storm records.
2) Potential misinterpretation and overreliance on the Ti results
- Why was Ti selected as the element most indicative of mineral content over other elements? I think more justification of this decision is needed. The Principle Component Analysis in figure 4 shows that for TLM Ti is not strongly correlated with PC1, when the ASI results and other elements are. Figure 7 which shows the ASI overlaying the Ti record does support the use of Ti and perhaps this comparison should be done earlier in the paper (and using the full records), as currently much of the ITRAX results are discounted with little explanation in the results section.
- The event frequency was calculated using the number of times that Ti exceeded a threshold. The data was detrended using a 10/30 year moving average window and peaks above a threshold were marked as storm events. The authors state ‘To avoid overestimating the number of events due to consecutive Ti measurements associated with a single Ti peak, consecutive measurements were grouped together, and only the maximum value within each group was retained’. I wonder however whether some of the identified storm events are represented by single data points or small groups, which could be outliers or erroneous data points rather than sand layers. As noted above, in data that has not been normalized other factors could cause an increase or decrease in the kcps that is unrelated to mineral content in the peat. The authors should check the Mean Squared Error and total counts per second measurements through the records to see if any of the peaks in Ti are associated with anomalies and should be excluded as artefacts (see Lowemark et al. 2019).
- A point related to this one about the detection of storm events, is that the moving window used may be too narrow and single large events may be identified as several small events. If you look at Figure 5b and the record of TAC the red line of the moving average threshold goes up over some of the increases in Ti as a result of the concentrations increasing over a >10 year period (depending on the accumulation rate 10 years could be just 1-2cm within the peat). My concern is that if there was a large deposit of sand on the bog by a single large magnitude storm, then sand may settle on the bog at different depths due to the uneven surface, with some falling into pools and sinking down, some at the roots, others on the top of plants etc.. Further mixing could occur with root bioturbation potentially. You acknowledge in the discussion ‘downward movement of sediments in the spongey peat matrix may have dated sediments associated to Hurricane Dorien to an older interval’ to explain the timing of an event, so a similar mechanism could spread the sand from a single large event over a few cm’s within the consolidated peat. Therefore, this could lead to overestimation of the number of events if these larger increases in sand are removed by detrending and then small variations in Ti on these peaks are counted as events. Either the method should be adjusted to address this or these limitations should be discussed, with perhaps more cautious statements about the frequency of events.
3) Exclusion of part of the TAC record
You should be consistent about whether or not the period between ~600 BCE and ~1000 CE in TAC is useful as a storm record, or not given you interpret it as being not ombrotrophic. In Figure 5 you shaded out the section but in Figure 8 it is included. In the discussion you also observe that the Ti and storm events records of the two cores are similar over this section. While you have concluded that TAC was not ombrotrophic at this time, that is not to say that the peaks in Ti and sand shown by the TAC record were not deposited by the same storms that caused the deposits in TLM. Perhaps better not to discount this section of the TLM record but to rather use it with the caveat that deposits may come from other sources in addition to storms.
Minor and Moderate Corrections
Title – I would put ‘storm reconstruction’ as plural, and not sure Late Holocene should be hyphenated
Line 36 – ‘in Atlantic Canada ever’, I would suggest that this is changed to ‘on record’ or ‘to date’ rather than ‘ever’
Line 46 – ‘Answering these questions is aided by a long-term perspective’, the wording of this is a little unusual, perhaps ‘A long term perspective can help to answer these questions,’
Line 65 – the sentence ending ‘cyclone strikes’ needs citations. Currently these are put after the following sentence, but they should go after the first sentence, or the two sentences could be merged.
Line 68 – sentence ending ‘centennial data’. As with the above point, the citations should go after the first sentence rather than the following one.
Line 82 – NAO needs writing in full the first time it is mentioned
Line 97 – ‘composed of quartz and being the main sand source’, this is phased in an unusual way. I would suggest ‘ …quartz, which is the main sand source..’
Line 100 – ‘steady winds from all directions’. Remove this, as in a following sentence you discuss wind directions being seasonal and show this with the results in the appendix.
Line 111 – it depends on the journal guidelines perhaps, but shouldn't the ages be 'ka before present' and ka BP rather than just ka?
Figure 1d and 1e - TLM is in the center of the island so could have sand from a few directions. Perhaps also include a couple more arrows with distances to beaches. Also the distance of TAC to the cliffs may not be as relevant to highlight as the distance to sand sources.
Line 182 – I don’t think Bjorck has an ‘l’ at the end
Line 202 – in equation 2, were the denominator measurements all the elements measured, or just those terrigenous mineral elements listed above.
Figure 3 – some of the y axis labels are too close together and are overlapping
Line 345 – The statement that the storm records are based solely on the Ti results needs more justification
Figure 5 - 1) Could the colored dots on the storm events be made smaller? It is hard to see the results. 2) While the results on this figure are convincing, as they show similar patterns, I would like to also see the ASI results presented, as these results are surely also as relevant as Ti to showing past storminess. 3) The AMV results should be included on 8, rather than 5, as this is where they are discussed.
Line 402 – this sentence says that sand layers in TAC at 600BCE and 810 CE were not visible in TLM. But there was an increase in identified storm events coinciding with these times in TLM, so could it be that there was an enhanced sand deposition related to storms at these times, but just more sand reached TAC?
Line 403 – The sentence says that there was no abrupt contacts between the sand and overlying and underlying peat, supporting a gradual accumulation of sand over time. I am not sure about this interpretation as peat bog environments don’t seem to often have sharp boundaries in the same way as lakes. The surface of bogs are uneven and I could see that when sand blows in a single event over the surface it could fall into pools, sit on top of plants but also land or wash down to the base of the plant and so be incorporated within different depths giving a gradual boundary even for a single large event.
Line 406 – for the sentence about the likely explanation is the in-situ concentration of mineral matter I think you should say that this interpretation is based on the hiatus shown by the age model, and any other evidence you have. I was confused at first why this would be a more likely explanation than the proximity of the sand sources.
Line 413 - I would not include the Netherlands example, as the Netherlands are so far away and not relevant to the regional climate. Lots of other sites around the world probably show no change at this time.
Line 422 – the sentence suggests the proximity of TAC to the cliffs is the reason for the larger particles at TAC. But are the cliffs the sediment source during storms? Are sand sized particles being eroded from the cliffs and transported in land during storms? It looks like TAC is slightly closer to the beach to the west than TLM is, but the satellite photos seem to show the environment around TAC is dryer, so would exposed soil potentially be providing a source of minerals to TAC?
Line 426 – in this sentence again the assumption is that the windblown sands came from the cliffs and/or beaches. I expect the soil on the island would have a similar elemental composition to the bedrock, so could wind blown soil also be a contributor and spatial variations in vegetation a factor?
Line 431 – remove typo (sand written twice)
Line 432 - ‘Ti emerged as the common aeolian sand indicator in both TLM and TAC.’ I think this needs more justification earlier in the paper.
Line 446 – remove the link to Fig 6c because it doesn’t show evidence that supports the use of ASI and Ti
Line 458 – ‘winter westerlies storms’, maybe 'storms associated with the winter westerlies' would sound better here
Line 461 – recent calendar years need CE throughout this section
Line 467 – ‘despite it coming ashore’, its not clear what this means - did waves come ashore?
Line 526 – ‘active hurricane period’ perhaps
Section 5.3 – this is a really nice comparison of the different records along the transect
Figure 8 – 1) you could just have CE on the x axis label. 2) in the caption you write that there are similarities between the Ti record from TAC and the storm activity in TLM. The fact that you have two records close to each other showing similar patterns increases confidence in the records and I think this point should be made more prominent in the text. It also gives some support that the non-ombrotrophic section of TAC is still capturing a storm signal.
Line 560 - a citation is needed here
Conclusion – this part highlights some interesting challenges and lessons but doesn’t go into much detail about the results of the paper, in the same way as the abstract does for example. I would like to see the results summarized in the conclusions either at the start or incorporated in. An example is the paragraph at line 631, where there could be another sentence covering the point you make about regional versus tropical Atlantic SST's being important for hurricane impacts in different regions.
Appendix A - 30-40, 40-50 with hyphens might make this clearer
Citation: https://doi.org/10.5194/egusphere-2025-400-RC1 -
RC2: 'Comment on egusphere-2025-400', Anonymous Referee #2, 24 Mar 2025
reply
SUMMARY: The authors present a detailed multi-proxy reconstruction of climate and environmental change over late Holocene timescales based on two cores from two different peat bogs on the Magdalen Islands, eastern Canada, and compare these to recorded events of storms in the past in order to facilitate a storminess reconstruction for the last 4000 years. The analyses involve a description of lithology, 14C- and 210Pb-based chronologies, XRF-core scanning, as well as measurements of LOI and grain size. Single elements identified in the XRF data, the inferred mineral content and the relative contribution of “sand” (based on the “Aeolian Sand Index”, Björck & Clemmensen, 2004) for the fractions <63µm and <125µm are used as proxies for storminess. Overall, the manuscript is well-written and thoroughly prepared and builds on sound and timely methods. There is a lot of positive things about this manuscript, but my main concern is the (uncontested) validity of the data as a storm record. This is for three main reasons: (1) The authors do not discuss local landscape evolution and potential changes regarding the sources of the allochthonous material blown into the bogs by storms or other dynamics that could have influenced the amounts deposited (such as exposure or storm direction). (2) The used grain sizes for ASI differ from other applications (which of course is not an issue), but are neither explained nor argued for, in terms of source properties or the processes determining transport (i.e., storminess). (3) There are a couple of issues with the chronologies (which is not uncommon for a natural system), but my critique is that the inferred age models are used a little too uncritical in the reconstruction of storminess and the discussion of climate in the big picture. That said, a huge strength of the manuscript is innovative combination of different parameters and methods (to my surprise very little GIS/geodata), which is extremely useful for depositional environments at the land-ocean interface and has a lot of potential. An ideal correspondence between all parameters cannot be expected, but I would regardless of that strongly encourage the authors to revise the manuscript thoroughly and argue for their data set more carefully and convincingly. I have added a couple of comments, both major and minor, below (ordered by line number).
L79: Small typo here and again further below. The author’s name is “Björck”.
L131ff: A little later on you mention a peat dome. Is there anything for the field site that is worth mentioning regarding to local geomorphological context (do you have access to LiDAR data or another type of digital elevation model for the island)? A quick look at Google Earth reveals that there is a prograded beach-ridge system located to the NE of your site TLM. Is the age of the system known? Its presence suggests that beaches have been as close to the sampling site as ~1200 m at some point in the past. Is there anything known about local coastal evolution at TLM and TAC (or for the island as a whole)? (There are a couple of really interesting studies by Billy et al. for the nearby island of Miquelon). The youthern coast around TAC appears to be a rock coast mostly, but there are smaller (pocket) beaches nearby. Have beaches ever been wider here? In analogy to that: Is there anything known about the past vegetation or land use history of the Magdalene Island (e.g., have the forests ever been cleared and has the peat ever been harvested for fuel)?
L144ff: How long are the cores overall and what’s their surface elevation? Did you reach the bottom of the bog? If a moraine, a marine sequence or something similar was reached, this would be very useful as supporting information of the chronology.
L149: A bit of unnecessary information here. Important is that they were stored cold.
L150: What is the “correction factor”? A percentage? This is a bit unclear.
L161ff: I am not sure to understand. This means that there is no valid Pb-based chronology for the TLM monolith, right?
L173ff: I would suggest to stick with a simpler nomenclature for the LOI processing stages and to use dry weight (DW or DW105) and LOI550 for combusted samples. Also, I feel that drying peat samples at 105°C overnight appears a little short.
L179: Unlike LOI, the ASI approach is not a very well established as a method and does require a more detailed description or more thorough and precise referencing. Did you consider the whole grain-size distribution of the samples or was there a significant amount of larger minerogenic particles present? It would be worthwhile to reference a bit more here and argue for the approach. Also, it reads here as if the fraction <125 µm is used as a proxy, but in the discussion, it is clearly stated that >125 µm is used. What did you use? Also, I would suggest to separate LOI and ASI, or change the title to “Variations in mineral content and ASI”. LOI is often used to identify the content of organic material in a sample, and I was a little confused that you went straight to the residual here.
L192: Often (but not always) LOI correlates well with measured values of TOC based on element analysis, so arguably LOI is an absolute measure (with a bit of uncertainty and the occasional outlier). Would the results look much different would you skip the data processing as described?
L228: You mention the identification of (non-)aeolian processes here. Did you by any chance measure full grain-size distributions on the samples from the core or on a few of the samples from the bottom of the cores at least? This would be very useful to pinpoint the differences in sources and processes.
L229f: The last sentence is very interesting, but does not really fit into the methods and would be more suitable for the discussion.
L243ff: In a similar fashion, the interpretation of the Ti outliers should be moved to the discussion.
L249ff: This sounds very promising and is a huge effort on your part! I look forward to read, how that integrates into indications from the lead-dated monoliths.
L261f: The sand layers at the bottom of the core are relatively thin. How can you be certain that they are a basal sand layer that marks the bottom of the peat sequence (and not just an interbed, e.g. from overwash)?
L265: Arguably, the marked drop in PAR could be interpreted as a hiatus rather than a stage of very slow accumulation (which slightly overstretches the ability of the 14C chronology). Also, there is a lot happening in the core at about the same time (or slightly thereafter). The peat looks much brighter in the core images and the ASI picks up up-core from here. Wouldn’t it be advisable to add an uncertainty to the rates and to subdivide a little more, especially in the case of TAC?
L272: Figure 2. I do not find the top panels above the age-depth models particularly informative (also, the last bit of text to the right is cut off) and you should consider to remove them. It would be helpful for the readers understanding to have a table of the radiocarbon ages with uncalibrated and calibrated ages (ideally as ranges), the 2-sigma probabilities, and info on the material dated and possibly the lithology of the unit dated. (NB: Sorry! I just found the table in the appendix – it could be useful to have this integrated into the main text). While the age model for TLM looks fairly straight forward, there is a tricky hiatus (or continued deposition followed by erosion?) at 150 cm. Would there be a good way (visually) to scale the y-axes the same way? The fact that the oldest peat deposits have a surprisingly similar age is easily overlooked!
L280ff: I feel it could be sensible to exclude the ASI (as a separate proxy) from the description of the lithology, which is (mostly) descriptive otherwise.
L300: Figure 3. I find it difficult to compare the sites visually. Would it be an option to display the data scaled for age and displayed with a linear axis (even distance of the tick-marks between e.g. centuries)? In a similar way, the values for ASI and mineral content should be scaled the same way and you could consider to use a log scale to make the small variations more visible. Some of the labels are not entirely visible in the current version (e.g., Si, ASI).
L347: I am a little surprised, that zone 3 of the fore from TAC was not excluded. The hiatus in the record (and the lack of information to fully explain it) is a bit puzzling, but if your storm proxy is the input of allochthonous material into the peat bog, a hiatus can not be ignored, as the entire zone 3 is potentially the time-integrated results of almost 2kyr of all kinds of processes (limited plant growth, a higher ground water table, flooding with marine water, etc.) potentially paired with a near constant influx of sand.
L365: I would suggest to remove the “and hence storminess” but here in the results section. The result is that you detected periods of enhanced mineral inputs, the rest is an interpretation and should ideally be discussed.
L383: There is a very nice match between the two reconstructed event densities for the 880-550 BCE window, that unfortunately falls into the hiatus at TAC (and hence in a time window, where the age model may be highly unreliable. It is very interesting regardless, but I really feel that you should discuss this more. For the 600-800 CE window, I don’t see a way the age model can be trusted for TAC. The phase is nicely visible in the Ti record for TLM, but bare visible in the ASI record.
L430: Not sure, I can follow the argument. Earlier you argued that the XRF values are relative and (unless I missed it) no absolute values were measured. In the data displayed in Figure 3, I do not see an elevated Fe content in the TAC core. The highest relative values are in fact found in the bottom section of core TLM.
L430: A depletion in minerals other than quartz is not a necessity and highly dependent on e.g., local properties (lithology of the source rock), the sediment budget of a site (and hence the time of exposure in the system prior to deposition), and transport distance. In order to make this understandable for all readers, the local properties need to be described more.
L433: It would be useful to show the correlation between the two different (!) ASI values and elements in the manuscript or the supplement.
L440ff: My feeling is, that this section is a little too sloppily written for the strong conclusion made in the end. You should lay out more of your considerations and explain the acting processes.
L441: “relatively unconsolidated“ is difficult to understand as a description. If the sources allow for transport of sediment to a distal site during storms, it would be interesting to know what happens at a more proximal location. Did you observe the formation of primary dunes on the cliff top or the beaches and what are the grain-size properties here? In other words, what happens in the source areas during other conditions than storms?
L449ff: It would be useful to know more about the typical wind directions during these different storm events. This does have a huge a huge impact on how the sources and your storminess archives connect.
L456: I like to use dashes in more informal contexts, like emails, but think that commas would be more appropriate here.
L462: This is news! My understanding up to this point was, that the beaches NE of TLM would be the major source of sediment to the site or those immediately to the south.
L448ff: This is a very interesting section and you are really working at the limits of the ASI approach here. You are considering a time window of ca. 170 years with an averaged accumulation rate on the order of 3 mm year (much less for TLM based on the age model shown in Fig. C1). It would be really interesting to hear your opinion here: Is it at all possible to use the mineral content of peat to get to a storminess reconstruction (ideally with the attribution of single events) for these time scales? Have you considered to lump the data (e.g., to decades)?
L507-608: This section is very important, but I struggle a bit with three things here: (1) The storminess reconstruction stands a bit on shaky grounds for a couple of reasons (understandable reasons, not necessarily flaws!), but the identified storm periods are here taken as fact and I miss the careful arguing a bit, (2) Your results are placed in general climate context here, but while reading through this I repeatedly wondered, if the purpose is to argue for the validity of ASI (and your age model for TLM) or to show the added value of your records, and (3) Several different scales are naturally mixed here and your hypothesis is stuck between the coast of West Africa and distance relationships between source and sink in sediment transport. To summarize, I like the section overall, but I am not sure of its purpose.
Citation: https://doi.org/10.5194/egusphere-2025-400-RC2
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
121 | 155 | 5 | 281 | 4 | 3 |
- HTML: 121
- PDF: 155
- XML: 5
- Total: 281
- BibTeX: 4
- EndNote: 3
Viewed (geographical distribution)
Country | # | Views | % |
---|---|---|---|
United States of America | 1 | 129 | 44 |
Canada | 2 | 39 | 13 |
China | 3 | 18 | 6 |
Ireland | 4 | 18 | 6 |
France | 5 | 17 | 5 |
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1
- 129