Abstract
In the coming century, climate variability is projected to increase in Northeast Pacific coastal areas, increasing the need for land managers to understand how ecosystems are expected to change in response to new or enhanced disturbances. Previous research indicates that the Pacific coast of Canada, southern British Columbia (BC) experienced warmer and drier climate conditions than present, with higher than modern fire activity during the early Holocene xerothermic interval (9.5 – 7.0 kcal BP). In this study, we reconstructed past climate-fire-vegetation changes from a 13,000-year sediment record from Lost Lake in Vancouver’s Lower Seymour Conservation Reserve, BC. Contrary to other sites, the coastal western hemlock forest at this site remained cool and moist with low fire activity throughout the xerothermic period. Instead, peak fire frequencies were observed during the cool and moist Neoglacial period (4.5 kcal BP – present), when human activity became prevalent. These results have implications for the managed watershed’s resilience to fire and response to future warming conditions.
1 Introduction
Understanding past changes in climate and ecology is a useful component of managing the environmental shifts associated climate, land management, and anthropogenic climate change. By the year 2050, temperatures in British Columbia are estimated to increase by 1.3 – 2.4°C under the Shared Socioeconomic Pathway 5-8.5 (SSP 5-8.5) future climate scenario from the sixth phase of the Coupled Model Intercomparison Project (CMIP6) (), and regional wildfire activity and severity are expected to intensify (). While impacts of wildfires on contemporary forest and water resources along the northwest coast of North America have been well documented (; ; Zwolinksi, 2000; ; Sankey et al., 2017), the effects of changing climate-fire-vegetation interactions in South Coastal BC in response to increasing temperatures and wildfire severity are less clear.
Paleoecology is a well-regarded method of understanding past interactions between climate and vegetation. Quantifying past disturbances and assessing how ecosystems respond to change is a means of identifying and evaluating disturbance regimes, helps with our understanding of the timing and significance of ecosystem shifts, and supports planning for future land management and restoration projects (; Swetnam et al., 1999; ; Whitlock et al., 2015; ; Pellatt et al., 2015; ). The Holocene epoch (11.7 kcal BP – present; ) is of interest to the understanding of long-term changes in climate and wildfire activity throughout western North America because the scales of temperature and precipitation variability may have been larger than those seen in the past 200 years (; ). Two commonly used proxies that have been applied to better understand changes in vegetation and fire behavior are palynomorphs (; ; ; ; ; ; ; ; ; ; ; , ; ; ) and charcoal particles, respectively (; ; ; ; , ; Sugimura et al., 2008; ; ; ; Pellatt et al., 2015; ; ; Prichard et al., 2009 (Supplementary Table S1; Figure 1).
Figure 1
Regional Holocene records of the past ~12,000 years indicate that the forested ecosystems in and around the Fraser Valley of British Columbia (BC) experienced large disturbances over relatively short (centennial to several thousand year) timescales (Walker and Pellatt, 2003). Previous paleoclimatic studies of British Columbia have divided the Holocene into major periods of climatic change. The Younger Dryas, occurring approximately between 12.9 – 11.6 kcal BP, marked a period of cooling in much of the Northern Hemisphere around the North Atlantic Ocean (
Figure 2

Synthesis of regional climatic data in the Northern Hemisphere from 14,000 kcal BP to present. Grey bars represent timing of Lost Lake Pollen Zones LL-1 to LL-6b. Colored bars represent timing of major climatic events: Bølling-Allerød (B-A, purple,
From ca. 2.4 kcal BP to present, multiple shorter, discrete climate intervals have been observed in coastal BC. An increase in fire activity and drought frequency between 2.4 and 1.3 kcal BP has been referred to as the Fraser Valley Fire Period (FVFP) (
In addition to climate driven change in ecosystem structure throughout the Holocene, Indigenous people have managed their traditional lands to improve foraging capacity, and for gardening, hunting, and habitation for millennia (
In this study, we use palynological and charcoal analyses to understand climate-fire-vegetation changes within the Lower Seymour Conservation Reserve (LSCR) in North Vancouver, BC, from the very end of the Pleistocene period to the present day (approximately 13,900 years). We place this record in the context of regional climatic changes throughout the Holocene to interpret long-term ecosystem changes in the LSCR. Investigating the mechanisms behind ecosystem changes during the Holocene provides insight into how local watersheds may respond to climate change in the coming decades and can be used to inform adaptation strategies within the Metro Vancouver water supply area (WSA).
2 Materials and methods
2.1 Study area
Lost Lake is a small (3.8 ha) lake located within the Metro Vancouver-governed Lower Seymour Conservation Reserve (LSCR), situated in the lowlands of the southernmost range of the Coast Mountains in North Vancouver, British Columbia (Figure 1). Late fall and winter (October to March) in this region are characterized by northeast Pacific storms which deposit a cumulative average of 112 cm of precipitation between October and March (
Although no archaeological evidence of pre-Colonial habitation has been directly observed at Lost Lake, the Lower Mainland and surrounding areas lie on the unceded territories of several Coast Salish First Nations, including those of the Musqueam, Squamish, and Tsleil-Waututh. Archaeological evidence dating as far back as 9.5 kcal BP has been found in the Fraser River Delta, which lies directly to the south of Lost Lake (
2.2 Field methods
In September 2020, we retrieved a 3.64 m piston core from the deepest part of Lost Lake (12 m) in six separate drives using a modified Livingston piston corer (Wright et al., 1984). We also collected a 0.41m surface core from the same location in November 2020 using a Glew gravity corer (
2.3 Chronology
To estimate the sediment mass accumulation rate and age of the core, we created an age-depth model using tie points based on 210Pb and AMS-14C dates and regional stratigraphic markers (i.e., the Mazama tephra). Eleven samples were selected from the surface core and sent to Flett Research Ltd. in Winnipeg, Manitoba, for 210Pb analysis (Supplementary Table S1). Additional 137Cs analyses were conducted on four samples, and 226Ra analyses were conducted on three samples (Supplementary Table S1). The age model for the Lost Lake surface core was created based on a constant rate of supply (CRS) model (
Five macrofossil samples from the Livingston core (Supplementary Table S2) were sent to Beta Analytic in Marathon, Florida, for AMS-14C dating. In preparation for dating, the macrofossil samples were washed with distilled water and dried in an oven overnight at 30°C. As per Beta Analytic’s standard pretreatment protocol for plant material, the samples underwent a hot acid (HCl) wash to remove carbonates, then an alkali wash (NaOH) to eliminate secondary organic acids, followed by a final acid rinse (
The age model for the composite core (surface core + piston core) was constructed using the Bacon modelling program in R (
Figure 3

Bacon age model of the Lost Lake Core. Grey dotted lines represent error envelope at 95% confidence interval. Grey fill represents all likely age-depth models calculated by Bacon. Red line represents mean age (the “best” model selected by Bacon). Blue points represent AMS-14C samples (2σ probability distributions of calibrated 14C ages). Green points represent locked in 210Pb dates and Mazama ash date. Gray vertical line represents slump point at Mazama tephra to indicate instantaneous deposition.
The mean sample resolution for the composite core was 48.8 yr/cm, with a much higher resolution (5.3 yr/cm) in the first 41 cm of the core. The mean sedimentation rate for the composite core was 0.086 cm/yr, with a maximum of 1.2 cm/yr between 12 and 15 cm of the composite core.
2.4 Pollen analysis
The Lost Lake composite core was sub-sectioned for pollen analysis at 10 cm intervals along the entire length of the core, with a 1-cm3 sample being removed at each sample point. Pollen preparation methods were adapted from standard recovery techniques (
Figure 4

Pollen and spore grains from the Lost Lake core. Pollen grains are identified as: (A)Tsuga heterophylla (western hemlock); (B)Pseudotsuga menziesii (Douglas fir); (C)Abies (true fir);(D)Picea (spruce) in foreground and T. heterophylla behind; (E)Pinus (diploxylon) (vesicular grains) and Polypodium (polypody fern); (F); Alnus (alder);(G) Cupressaceae (cedar); (H)Acer (Maple, top center) and folded T. heterophylla. All grains identified using keys from
Figure 5

Pollen percent diagram of Lost Lake. Radiocarbon dates and lithology are shown on the left, fire frequency and total terrestrial pollen concentration are shown on the right. Tenfold exaggeration curves (light grey) are shown to highlight abundance of infrequent pollen and spores. Zones were delineated using constrained cluster analysis (CONISS). Solid horizontal lines represent zonal breaks. Dotted horizontal lines represent breaks in the core.
2.5 Charcoal analysis
Analysis of macroscopic charcoal followed methods commonly used in western North America modified from
Charcoal analysis was conducted using the free software package CharAnalysis (
Fire episodes were represented by the high-frequency CHAR (Cpeak) component of the record that exceeded a threshold value which isolated fire-related peaks from non-fire related peaks. To separate fire-related peaks from non-fire related peaks, a Gaussian mixture model was used to define the noise distribution, with threshold values limited to the locally defined 99th percentile of the noise distribution. The cut-off probability for minimum counts was set to 0.05, indicating that the minimum charcoal count within 75 years before a given peak was required to have a less than 5% chance of coming from the same Poisson distribution as the maximum count associated with said peak, otherwise the peak was removed (
3 Results
3.1 Lithology
We identified four lithological units along the 294 cm length of the composite Lost Lake core. The deepest seven centimeters of the core (294 – 287 cm) was composed of gyttja. A fine-grained, light-gray clay from 287 to 281 cm was inferred to be a result of glacial scouring, and indicated that the 294 – 287 cm section was likely re-cored younger sediment caused by shifting during extraction of the Livingston corer. An inverted radiocarbon date from this section supports this hypothesis, and the 7-cm section was removed from the final lithology. Undifferentiated, dark brown gyttja was present from 281 to 269 cm, and was overlain by a 5-cm layer of lighter brown gyttja intermingled with small (1 – 5 mm) angular clay inclusions from 269 to 260 cm. From 260 to 206 cm, the sediment was composed dominantly of dark brown gyttja, with a mottled, inconsistent layer of presumably displaced Mazama tephra observed from 233 to 228 cm. Sediment for the depth interval of 206 – 204 cm consists of a consolidated light brown tephra from the Mount Mazama eruption. From 204 – 0 cm, the sediment was an undifferentiated dark brown gyttja containing few macrofossils.
3.2 Pollen analysis
Six pollen assemblage zones were identified in the Lost Lake composite core (LL-1 to LL-6, Figure 5) based on a constrained cluster analysis using CONISS total sum of squares (
3.2.1 Zone LL-1 (287 - 251 cm; >12.2 kcal BP)
Pinus reached its highest pollen percentage (>50%) at the beginning of Zone LL-1 (13.6 kcal BP) and then decreased substantially over the subsequent ca. 1400 years, reaching 9% by the end of the zone. Pseudotsuga/Larix pollen abundance remained at <1% from the base of the core until 12.7 kcal BP when it rapidly increased to its maximum value (~28%) during the transition between LL-1 and LL-2 at ca. 12.2 kcal BP. Picea (inferred to be Picea sitchensis) pollen gradually increased to its maximum percentage of 16% at 13.1 kcal BP and then steadily decreased to 2% by the end of Zone LL-1. Tsuga heterophylla abundance increased from <1% during the oldest part of the core to 12% at 12.7 kcal BP.
3.2.2 Zone LL-2 (251 – 216 cm; 12.2 – 9.1 kcal BP)
Tsuga heterophylla increased from 12% at the beginning of Zone LL-2 to ~29% by 9.1 kcal BP. Alnus decreased from 37% to 30%. Pinus continued to decline to ~1.5% of the assemblage by 9.1 kcal BP. Abies and Picea pollen abundances remained steady but low between 0.5 and 2%, respectively.
3.2.3 Zone LL-3 (216 – 175 cm; 9.1 – 6.7 kcal BP)
Pinus percentages remained below 3.5% for the duration of Zone LL-3. Tsuga heterophylla and Cupressaceae increased to 40% and 10%, respectively. Pseutotsuga/Larix decreased to 7% by the end of Zone LL-3. Shrub and herb pollen species experienced their highest percentages, with Rosaceae, Salix, and Artemisia reaching 1.3%, 1% and 0.55%, respectively.
3.2.4 Zone LL-4 (175 – 110 cm; 6.7 – 3.9 kcal BP)
Tsuga heterophylla abundances decreased from ~40% to 28 – 34% during Zone LL-4. Cupressaceae increased from 10% to nearly 40% by 3.9 kcal BP. Alnus pollen fluctuated between 15 and 20%. Pseudotsuga/Larix declined to <5% by the conclusion of Zone LL-4. Pinus, still dominated by P. contorta, remained below 3.5% of the pollen assemblage.
3.2.5 Zone LL-5 (110 – 75 cm; 3.9 – 2.4 kcal BP)
Tsuga heterophylla increased slightly to ~39% by ca. 2.4 kcal BP. Cupressaceae declined briefly to 24% before again rising to >30%, while Alnus fluctuated between 13 and 19%. Pinus and Pseudotsuga/Larix abundances fluctuated between 0.5-4% and 2-5%, respectively.
3.2.6 Zone LL-6a (75 – 54 cm; 2.4 – 1.2 kcal BP)
Tsuga heterophylla reached its maximum contribution of 50% of the assemblage at 1.7 kcal BP, at the expense of Cupressaceae, which decreased from >30% to 9% during Zone LL-6a. Alnus abundances remained between 17 and 22% of the assemblage. Pinus percentage steadily increased from 2.5 to 4% by 1.2 kcal BP, while Pseudotsuga/Larix contributed 4-6% of the pollen assemblage during Zone LL-6a.
3.2.7 Zone LL-6b (54 – 0 cm; 1,200 – -71 cal yr BP)
Alnus abundance increased from 17-22% to 36% at ca. -35 cal yr BP (~1985 AD). Cupressaceae sharply increased to 32% at ca. 140 cal yr BP at the expense of Tsuga heterophylla, which decreased to 32% before climbing back to a peak of 42% at approximately 20 cal yr BP (1930 AD). Pinus and Pseudotsuga/Larix abundance remained between 3- 6% and 2.5-5%, respectively. In all arboreal species and the majority of herb and shrub species pollen accumulation rates (Supplementary Figure S1) appear to emphasize rapid increases in influx in the final 200 years of the record. This pattern was not observed in the pollen percentage records of most species, but the percentages of Alnus rubra and Poaceae display moderate increases that follow the trend of their PARs.
3.3 Charcoal analysis
A total of 23 significant fire episodes were detected during the 13,900 years of the Lost Lake charcoal record (Figure 6; Supplementary Table S4). The mean CHAR for the Lost Lake composite record, interpolated to 44-year sample intervals, is 14 pieces/cm2/year. The period of highest CHAR occurred from -71 cal yr BP to 61 cal yr BP during Zone LL-6b (Figure 6). The current time since last fire (TSLF) is 144 years, and the most recent detected fire event occurred in 1845 AD. The longest interval of no recorded fire activity was approximately 1230 years between 12.5 and 11.2 kcal BP in Zone LL-2, while the mean fire-free interval was 555 years. The mean fire return interval (mFRI) for the composite core was 598 years with a natural range of variability of 466 - 735 years.
Figure 6

Synthesis of Lost Lake charcoal-inferred fire activity and major anthropogenic and climatic changes in southwestern British Columbia since the Late Glacial. Red fill represents interpolated CHAR; black solid line represents fire frequency; dotted black line represents interpolated background CHAR; X symbols represent fire events; Vegetation Zones LL-1 – LL-6b and major climatic intervals are shown at the top of the figure. Vertical colored bars represent the FVFP (orange), MWP (pink) and LIA (blue). Dotted grey vertical lines represent breaks in the Lost Lake core.
3.3.1 Zonal charcoal record
The charcoal record of Lost Lake was assessed in greater detail by using the CharAnalysis output data to calculate the fire characteristics for the individual pollen assemblage Zones LL-1 to LL-6 (Supplementary Table S4).
Fire frequency in the deepest section of the Lost Lake core was moderate (approximately 1.3 fires/1000 years) and remained low during the post-glacial, pine-dominated assemblage of Zone LL-1. Fire Frequency began increasing ca. 11.6 kcal BP and displayed a generally increasing trend throughout Zones LL-2 to LL-4. Fire frequency peaked at ca. 6.0 kcal BP when T. heterophylla and Cupressaceae dominated the pollen assemblage. A rapid decrease in fire frequency occurred between ca. 6.0 and 3.4 kcal BP during Zone LL-4 and was followed by a similarly rapid increase. Fire frequencies reached maximum values at 2.3 kcal BP. Fire frequency subsequently declined after 2.3 kcal BP with a minor increase from 550 cal yr BP to the present.
CHAR values have fluctuated substantially throughout the record but generally remained low-to-moderate in the earliest half of the core, with intermittent peaks throughout. CHAR began noticeably increasing in Zone LL-3 at approximately 7.3 kcal BP and remained elevated until ca. 4.0 kcal BP. Following 4.0 kcal BP, CHAR declined until its final small increase between approximately 550 cal yr BP and the present. Background CHAR generally aligns with the CHAR trend: it remained moderate-to-low during the earliest part of the Lost Lake core until approximately 7.5 kcal BP, at which point it increased substantially. Background CHAR levels peaked at approximately 3.6 kcal BP and then declined and remained low for the remainder of the core.
4 Discussion
4.1 Zonal interpretations
4.1.1 Zone LL-1 (287 - 251 cm; >12.2 kcal BP)
With an estimated basal age of 13.9 kcal BP, Zone LL-1 records the transition from the last glacial period into the Holocene period (11.7 kcal BP– present) (
As the climate moistened ca. ~13.0 kcal BP (Walker and Pellatt, 2003;
The Younger Dryas period (YD; 12.8 – 11.5 kcal BP) was a period of cooling in the Northern Hemisphere which was first observed in northwestern Europe (Rind et al., 1986), and later linked to coeval ecosystem-scale shifts in vegetation in the Rocky Mountains (Reasoner et al., 1994; Reasoner and Jodry, 2000), western North America (
4.1.2 Zone LL-2 (251 – 216 cm; 12.2 – 9.1 kcal BP)
At the end of the Younger Dryas chronozone, temperature began to increase across much of the Northern Hemisphere (
The inferred shift to warmer, drier conditions around Lost Lake after 12.0 kcal BP aligns with rising solar insolation values, which peaked at approximately 12.0 kcal BP (Figure 2), and increasing SST in the northeast Pacific between ca. 12.0 and 10.0 kcal BP (Figure 2;
The fire-sensitive species Tsuga heterophylla gradually increased to near-modern levels beginning ca. 10.0 kcal BP (Figure 5). This increase aligns temporally with a decline in summer insolation values and increase in winter insolation values at 60°N (Figure 2), and is supported by similarly timed T. heterophylla increases further north along the coast of BC (
4.1.3 Zone LL-3 (216 – 175 cm; 9.1 – 6.7 kcal BP)
The early Holocene xerothermic period (ca. 9.5 – 7.0 kcal BP,
These results indicate that the xerothermic interval was not felt as strongly at Lost Lake in comparison to other sites in the Fraser Valley (
4.1.4 Zone LL-4 (175 – 110 cm; 6.7 – 3.9 kcal BP)
During Zone LL-4, the inferred vegetation assemblage was dominated by T. heterophylla and Cupressaceae, with lesser amounts of Betula (swamp birch) and Salix (willow), indicating a moist climate continued to dominate from 6.7 to 3.9 kcal BP. Pseudotsuga declined to near modern levels by ca. 6.0 kcal BP, indicating that the climate was no longer dry enough to support its growth.
A local peak in fire frequency of ~3 fires/1000yrs, with an average mFRI of 308 years, occurred ca. 6.0 kcal BP (Figure 6; Supplementary Table S4), and was accompanied by elevated values of background CHAR between 6.5 and 4.0 kcal BP (Figure 6). These changes in fire frequency, mFRI, and background CHAR indicate that biomass burning increased locally around Lost Lake. Evidence of human activity has been abundant along the coast during the early-to-mid Holocene (
While fire frequency decreased during the last half of Zone LL-4, CHAR remained high despite the inferred cool, moist climate. A potential cause for this elevated CHAR is burning that occurred either upwind of the lake or outside the watershed entirely (
4.1.5 Zone LL-5 (110 – 75 cm; 3.9 – 2.4 kcal BP)
Zone LL-5 occurred shortly after the beginning of the Neoglacial period (ca. 4.5 kcal BP), during which several glacial advances occurred in the Coast Mountains likely in response to declining summer insolation in the Northern Hemisphere (Figure 2; Ryder and Thomson, 1986;
Fire activity began to increase at Lost Lake at ca. 3.5 kcal BP, which coincided with the intensification of the positive (El Niño) phase of El Niño Southern Oscillation (ENSO) that is believed to have developed between 3.5 and 2.5 kcal BP (
4.1.6 Zone LL-6a (75 – 54 cm; 2.4 – 1.2 kcal BP)
The increase in fire frequency at Lost Lake also characterizes the beginning of Zone LL-6a. These increases in fire activity between ca. 3.5 and 2.3 kcal BP has similarly been observed in other paleoclimate studies in southwestern BC (Wainman and Mathewes, 1987;
During Zone LL-6a, pollen abundances at Lost Lake show little change, but a decline in T. plicata percentage at ca. 1.8 kcal BP is consistent with a slightly offset vegetational response to the increased wildfire and drought conditions inferred between ca. 3.0 and 2.0 kcal BP. From ca. 2.0 kcal BP until the end of Zone LL-6a at ca. 1.2 kcal BP, the pollen assemblage remained indicative of cool, wet conditions while fire frequencies declined. Moisture indicators Cupressaceae, Betula and Salix increased at approximately ca. 1.3 kcal BP, while T. heterophylla and P. menziesii remained stable, signaling only minor change to the overall forest composition.
4.1.7 Zone LL-6b (54 – 0 cm; 1,200 – -71 cal yr BP)
Lost Lake’s vegetation composition did not experience any significant changes between 1200 and 100 cal yr BP, indicating a relatively consistent environmental controls. In the past ~170 years (i.e., ca. 100 cal yr BP to present), fire activity at Lost Lake increased substantially, as evidenced by the large increases in CHAR, rising fire frequency, and high levels of disturbance in the pollen record.
The high CHAR beginning at 105 cal yr BP (AD 1845) was likely due to the increased disturbance caused by settler logging, mining and slash burning in the Lower Mainland. Increases in Pinus, Alnus, Pteridium, and shrub growth are likely due to higher levels of disturbance and stand gaps allowing increased understory growth. Cupressaceae, P. menziesii and T. heterophylla abundances decreased between ca. AD 1920 and AD 1940, most likely a direct result of logging activity.
4.2 Management implications
Looking ahead to the next century, watershed management in the LSCR has the potential to change substantially as vegetation, fire activity, and natural ranges of variability are altered by climate change (
An important implication of this warming and drying is a potential vegetational shift away from the moisture-loving T. plicata and T. heterophylla and an increase in fire-adapted species such as P. menziesii, perhaps moving towards a composition more similar to the drier subvariants of the CWH zone or the coastal Douglas fir (CDF) zones that currently exist in the rain shadow of Vancouver Island. Sites in the cool, moist variants of the CWH zone, such as the current conditions at our site, became more Douglas-fir dominated during the xerothermic period (
The current closed canopy forests surrounding Lost Lake indicate a higher fuel load and potential for high-severity, stand-replacing crown fires. Decreasing precipitation in the coming decades will likely dry fuels out, resulting in higher burn likelihood in the event of an ignition. Alternatively, a transition to intermediate precipitation and moderate to high fuel loads could result in mixed-severity fire regimes, causing a patchwork distribution of ground and crown fires that results in variable tree mortality, similar to what is seen in montane forests today (
We consider multiple hypotheses for the cause of Lost Lake’s perceived resilience to fire during its warmest period. The first is that the coastal setting of the LSCR created a climatic buffering effect as the cool, moist air of the Pacific Ocean dampened the insolation-driven summer temperature variations, somewhat protecting the site from extreme temperature and drought (
Similar to the inferred anthropogenic burning during the FVFP, the potentially intensive effects of human activity on fire frequency may have implications for watershed management as population density around the watersheds continues to increase. Housing developments in North Vancouver already border the southern boundary of the LSCR, and as the population of the Lower Mainland is projected to increase to 4.1 million by the 2040s (
5 Conclusion
This study presents new sedimentary pollen and charcoal records from Lost Lake that fill an information gap in the paleoclimate history of CWHvm1 forests in the coastal Lower Mainland of British Columbia. Our pollen record indicates that the period of highest temperature at Lost Lake occurred when P. menziesii was at its highest abundance (18-26% between ca. 12.2 and 9.2 kcal BP; Figure 5), which was several thousand years before the commonly described xerothermic interval (9.5 – 7.0 kcal BP). P. menziesii forest likely propagated at this time due to a combination of the inferred dry climate in the early Holocene (Walker and Pellatt, 2003) and the strong summer insolation (Figure 2) causing very warm summer conditions suitable to growth of the species. The low fire frequency and CHAR values during this time (Figure 6) can perhaps be explained by (a) lower overall biomass following deglaciation, or by (b) Lost Lake’s coastal proximity giving it a “moisture buffer” (
Climate at Lost Lake during the mid-Holocene was cool and moist, producing western hemlock-western redcedar closed forests. Higher fire frequency was likely due to a combination of increased anthropogenic landscape modification and elevated charcoal influx due to increased precipitation. In the late Holocene (3.0 kcal BP onward), the vegetation assemblage at Lost Lake signaled a continuous temperate and moist climate, but fire frequency peaked during the Fraser Valley Fire Period at ca. 2.4 – 1.3 kcal BP. Plant assemblages were likely similar to modern day, but prolonged and frequent droughts may have occurred in summer due to a stronger summer Pacific High and weakened winter Aleutian Low (
The interpreted climatic changes of this study indicate that natural vegetational succession, disturbance, broad-scale climatic changes, and human impact all contributed to the development of the current conditions found at Lost Lake. Based on our pollen and charcoal data and the province’s projected shift towards drier conditions within the century, the future vegetation and fire regime at Lost Lake may change in the following ways: (1) the fire season could lengthen, potentially with less available moisture during summers; (2) fire-adapted, drought-tolerant species such as P. menziesii and Pteridium may increase in abundance, and the ecosystem could shift towards a drier variant of the CWH zone; (3) fuel availability may rise due to increased mortality of arboreal species, potentially increasing fire risk if not managed; (4) continued population growth in the Greater Vancouver Area may increase the risk of ignitions in the WSA.
In conclusion, understanding the role of past climatic changes and ranges of variability provides a valuable perspective when considering future changes to the watershed that may occur as climate change progresses. Local-scale controls (i.e., topography/elevation, weather, human influence, fuel availability) undoubtedly play a part in Lost Lake’s paleoclimate history and should be considered when applying these results broadly to CWH forests on British Columbia’s southern coast. The potential effects of coastal moderation, combined with the apparent resilience of the forests around Lost Lake indicates that they are able to withstand large amounts of disturbance without major changes in assemblage, which bodes well for the future of the watershed.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: PANGAEA Data Repository (pollen data: https://doi.pangaea.de/10.1594/PANGAEA.968499; Charcoal data: https://doi.pangaea.de/10.1594/PANGAEA.968500).
Author contributions
MD: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Project administration, Resources, Software, Validation, Visualization, Writing – original draft, Writing – review & editing, Methodology. MP: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – review & editing. KK: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This project was supported by joint funding from Mitacs Accelerate and the Metro Vancouver Regional District to MD, and a Natural Sciences and Engineering Research Council (Discovery Grant RGPIN342251) to KK. Parks Canada provided lab space and materials. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Acknowledgments
We are thankful for the support provided by the Metro Vancouver Regional District, Parks Canada, and the Mitacs Accelerate program. Thank you to Dr. Rolf Mathewes for providing materials and assistance with pollen identification. Thanks to Hasini Basnayake, Jacqui Levy, Cass Perreira, and other members of the Climate, Oceans and Paleo-Environments Lab for assistance with laboratory and field work. We deeply appreciate the field assistance, site knowledge, discussions, and editing provided by D. Dunkley, which greatly improved the quality of the manuscript. Thank you to Dr. Scott Mooney and a reviewer for constructive comments that improved the manuscript. This manuscript is based on MD’s Masters of Resource and Environmental Management thesis (
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fevo.2025.1504983/full#supplementary-material
References
1
AbatzoglouJ. T.WilliamsA. P. (2016). Impact of anthropogenic climate change on wildfire across western US forests. Proc. Natl. Acad. Sci. – PNAS113, 11770–11775. doi: 10.1073/pnas.1607171113
2
AgeeJ. K. (1993). Fire ecology of Pacific Northwest forests (Washington, D.C: Island Press).
3
AllenG. B. (1995). Vegetation and climate history of southeast vancouver island, british columbia (Victoria (BC): University of Victoria).
4
AlleyR. (2000). The younger dryas cold interval as viewed from central greenland. Quat. Sci. Rev.19, 213–226. doi: 10.1016/S0277-3791(99)00062-1
5
ApplebyP. G.OldfieldF. (1977). The calculation of lead-210 dates assuming a constant rate of supply of unsupported 210Pb to the sediment. Catena5, 1–8. doi: 10.1016/S0341-8162(78)80002-2
6
BarnoskyC. W. (1981). A record of late Quaternary vegetation from Davis Lake, southern Puget Lowland, Washington. Quat. Res. 16, 221–239. doi: 10.1016/0033-5894(81)90046-6
7
BarronJ. A.AndersonL. (2011). Enhanced Late Holocene ENSO/PDO expression along the margins of the eastern North Pacific. Quat. Int.235, 3–12. doi: 10.1016/j.quaint.2010.02.026
8
BergerA.LoutreM. F. (1991). Insolation values for the climate of the last 10 million years. Quat. Sci. Rev.10, 297–317. doi: 10.1016/0277-3791(91)90033-Q
9
Beta Analytic (2022). Beta Analytic Standard Pretreatment Protocols. Available online at: https://www.radiocarbon.com/pretreatment-carbon-dating.htm (Accessed December 15, 2024).
10
BlakeM. (2004). “Fraser Valley Trade and Prestige as Seen from Scowlitz,” in Complex Hunter-Gatherers: Evolution and Organization of Prehistoric Communities on the Plateau of Northwestern North America. Eds. PrentissW. C.KuijtI. (Salt Lake City: University of Utah Press), 103–112.
11
BlaauwM.ChristenJ. A. (2011). Flexible paleoclimate age-depth models using anautoregressive gamma process. Bayesian Anal.6, 457–474. doi: 10.1214/11-BA618
12
BoydR. (1999). Indians, fire, and the land in the Pacific Northwest (Corvallis, OR: Oregon State University Press).
13
BrauerA.HaugG.DulskiP.SigmanD.NegendankJ. (2008). An abrupt wind shift in western europe at the onset of the younger dryas cold period. Nat. Geosci.1, 520–523. doi: 10.1038/ngeo263
14
BrownK. J.HebdaR. J. (2002a). Ancient Fires on Southern Vancouver Island, British Columbia, Canada: A Change in Causal Mechanisms at about 2,000 ybp. Env. Archaeol: J. Hum. Palaeoecol.7, 1–12. doi: 10.1179/env.2002.7.1.1
15
BrownK. J.HebdaR. J. (2002b). Origin, development, and dynamics of coastal temperate conifer rainforests of southern Vancouver Island, Canada. Can. Jour. For. Res.32, 353–372. doi: 10.1139/x01-197
16
BrownK.HebdaN.SchoupsG.ConderN.SmithK.TrofymowJ. (2019). Long-term climate, vegetation and fire regime change in a managed municipal water supply area, British Columbia, Canada. Holocene29, 1411–1424. doi: 10.1177/0959683619854523
17
BrownK. J.HebdaR. J. (2003). Coastal rainforest connections disclosed through a late quaternary vegetation, climate, and fire history investigation from the mountain hemlock zone on southern vancouver island, british colombia. Canada. Rev. Palaeobot. Palynol.123, 247–269. doi: 10.1016/S0034-6667(02)00195-1
18
CarignanR.ArcyP. D.LamontagneS. (2000). Comparative impacts of fire and forest harvesting on water quality in Boreal Shield lakes. Can. J. Fish. Aquat. Sci.57, 105–117. doi: 10.1139/cjfas-57-S2-105
19
CarlsonR. L. (1996). “Early Namu,” in Early Human Occupation in British Columbia. Eds. CarlsonR. L.Della BonaL. (University of British Columbia Press, Vancouver).
20
CarlsonR. L. (1994). “Trade and exchange in prehistoric british columbia,” in: Prehistoric exchange systems in North America. Interdisciplinary contributions to archaeology. Eds. BaughT. G.EricsonJ. E. (Boston, MA: Springer).
21
ChaseM.BleskieC.WalkerI. R.GavinD. G.HuF. S. (2008). Midge-inferred holocene summer temperatures in southeastern british columbia, canada. Palaeogeogr. Palaeoclimatol. Palaeoecol.257, 244–259. doi: 10.1016/j.palaeo.2007.10.020
22
ChengH.ZhangH.SpötlC.BakerJ.SinhaA.DongX.et al. (2020). Timing and structure of the Younger Dryas event and its underlying climate dynamics. PNAS117, 23408–23417. doi: 10.1073/pnas.2007869117
23
ClagueJ. J.MathewesR. W.GuilbaultJ. P.HutchinsonI.RickettsB. D. (1997). Pre-Younger Dryas resurgence of the southwestern margin of the Cordilleran ice sheet, British Columbia, Canada. Boreas26, 261–278. doi: 10.1111/j.1502-3885.1997.tb00855.x
24
ClagueJ. J.MenounosB.OsbornG.LuckmanB. H.KochJ. (2009). Nomenclature and resolution in Holocene glacial chronologies. Quat. Sci. Rev.28, 2231–2238. doi: 10.1016/j.quascirev.2008.11.016
25
ClarkJ. S.RoyallP. D. (1995). Particle-Size evidence for source areas of charcoal accumulation in late holocene sediments of Eastern North American lakes. Quat. Res.43, 80–89. doi: 10.1006/qres.1995.1008
26
ClarkJ.RoyallP.ChumbleyC. (1996). The role of fire during climate change in an eastern deciduous forest at devil's bathtub, new york. Ecology.77, 2148–2166. doi: 10.2307/2265709
27
ConstantineM.MooneyS. (2022). Widely used charcoal analysis method in paleo studies involving NaOCl results in loss of charcoal formed below 400°C. Holocene32, 1358–1362. doi: 10.1177/09596836211041740
28
ConstantineM.ZhuX.CaddH.MooneyS. (2023). Investigating the effect of oxidants on the quantification and characterization of charcoal in two southeast Australian sedimentary records. Fire (Basel Switzerland)6, 54. doi: 10.3390/fire6020054
29
CwynarL. C. (1987). Fire and the forest history of the North Cascade Range. Ecology68, 791–802. doi: 10.2307/1938350
30
DaviesA.BuntingM. (2010). Applications of palaeoecology in conservation. Open Ecol. J.3, 54–67. doi: 10.2174/1874213001003020054
31
DerrK. (2014). Anthropogenic fire and landscape management on valdes island, southwestern BC. Can. J. Archaeol.38, 250–279.
32
DohertyC. T. (2005). The Late-Glacial and Holocene relative sea-level history of the Seymour Inlet Complex, British Columbia, Canada. Queen’s University, Belfast, Belfast (Ireland).
33
DuncanM. E. (2022). Holocene fire and paleoclimate history of a small lake in the Lower Seymour Valley, British Columbia. Burnaby (BC): Simon Fraser University.
34
EyquemJ. L.FeltmateB. (2022). Irreversible Extreme Heat: Protecting Canadians and Communities from a Lethal Future (University of Waterloo: Intact Centre on Climate Adaptation).
35
FaegriK.IversenJ. (1989). Textbook of Pollen Analysis. 4th ed (Caldwell, NJ: Blackburn Press).
36
FastieC. L. (1995). Causes and ecosystem consequences of multiple pathways of primary succession at Glacier Bay, Alaska. Ecology76, 1899–1916. doi: 10.2307/1940722
37
FrieleP. A.ClagueJ. J. (2002). Younger dryas readvance in squamish river valley, southern coast mountains, british columbia. Quat. Sci. Rev.21 (18), 1925–1933. doi: 10.1016/S0277-3791(02)00081-1
38
FritzS. C. (1996). Paleolimnological records of climatic change in north america. Limnol. Oceanogr.41, 882–889. doi: 10.4319/lo.1996.41.5.0882
39
GallowayJ. M.DohertyC. T.PattersonR. T.RoeH. M. (2009). Postglacial vegetation and climate dynamics in the Seymour-Belize Inlet Complex, central coastal British Columbia, Canada: palynological evidence from Tiny Lake. J. Quat. Sci.24, 322–335. doi: 10.1002/jqs.1232
40
GallowayJ. M.PattersonR. T.DohertyC. T.RoeH. M. (2007). Multi-proxy evidence of postglacial climate and environmental change at Two Frog Lake, central mainland coast of British Columbia, Canada. J. Paleolim.38, 569–588. doi: 10.1007/s10933-007-9091-4
41
GardnerJ. J.WhitlockC. (2001). Charcoal accumulation following a recent fire in the cascade range, northwestern USA, and its relevance for fire-history studies. Holocene11 (5), 541–549. doi: 10.1191/095968301680223495
42
GauvreauA.FedjeD.DyckA.MackieQ.HebdaC. F. G.HolmesK.et al. (2023). Geo-archaeology and Haíɫzaqv oral history: Long-term human investment and resource use at EkTb-9, Triquet Island, N̓úláw̓itxˇv Tribal Area, Central Coast, British Columbia, Canada. J. Archaeol. Sci. Rep.49, 103884. doi: 10.1016/j.jasrep.2023.103884
43
GavinD. G.BrubakerL. B. (1999). A 6000-year soil pollen record of subalpine meadow vegetation in the Olympic mountains, Washington, USA. J. Ecol.87, 106–122. doi: 10.1046/j.1365-2745.1999.00335.x
44
GavinD.BrubakerL.LertzmanK. (2003). An 1800-year record of the spatial and temporal distribution of fire from the west coast of vancouver island, canada. Can. J. For. Res.33, 573–586. doi: 10.1139/x02-196
45
GavinD. G.MclachlanJ. S.BrubakerL. B.YoungK. A. (2001). Postglacial history of subalpine forests, olympic peninsula, washington, USA. Holocene11, 177–188. doi: 10.1191/095968301670879949
46
GavinD. G.HallettD. J.FengS. H.LertzmanK. P.PrichardS. J.BrownK. J.et al. (2007). Forest fire and climate change in western North America: Insights from sediment charcoal records. Front. Ecol. Env.5, 499–506. doi: 10.1890/060161
47
GavinD. G.HendersonA. C. G.WestoverK. S.FritzS. C.WalkerI. R.LengM. J.et al. (2011). Abrupt Holocene climate change and potential response to solar forcing in western Canada. Quat. Sci. Rev.30, 1243–1255. doi: 10.1016/j.quascirev.2011.03.003
48
GervaisP. (1867–1869). Zoologie et paléontologie générales: nouvelles recherches sur les animaux vértetebrés vivants et fossils (Paris: Arthus Bertrand), 263.
49
GiffordR.BrownC.BaronC.ClementD.MelnychukN.NelsonH.et al. (2022). British Columbia Chapter in Canada in a Changing Climate: Regional Perspectives Report. Eds. WarrenF. J.LulhamN.LemmenD. S. (Government of Canada, Ottawa, Ontario).
50
GlewJ. (1988). A portable extruding device for close interval sectioning of unconsolidated core samples. J. Paleolimnol.1, 235–239. doi: 10.1007/BF00177769
51
Government of Canada (2022). Canadian Climate Normals 1981-2010 Station Data for Haney UBC RF Admin. Available online at: https://climate.weather.gc.ca/climate_normals/results_1981_2010_e.html?stnID=776&autofwd=1 (Accessed April 10, 2022).
52
Government of Canada (2024). Canadian Climate Normals 1991-2020 Data for Vancouver (Harbour) British Columbia. Available online at: https://climate.weather.gc.ca/climate_normals/results_1991_2020_e.html?searchType=stnProv&lstProvince=BC&txtCentralLatMin=0&txtCentralLatSec=0&txtCentralLongMin=0&txtCentralLongSec=0&stnID=359000000&dispBack=0 (Accessed December 30, 2024).
53
GreenR. N.BlackwellB. A.KlinkaK.DobryJ. (1999). Partial reconstruction of fire history in the Capilano watershed (Burnaby, B.C: Report to the Greater Vancouver Water District).
54
GreenR. N.KlinkaK. (1994). A field guide for site identification and interpretation for the vancouver forest region (Vancouver, BC: British Columbia Ministry of Forests).
55
GrimmE. C. (1987). CONISS: a FORTRAN 77 program for stratigraphically constrained cluster analysis by the method of incremental sum of squares. Comput. Geosci.13, 13–35. doi: 10.1016/0098-3004(87)90022-7
56
GrimmE. C. (1990). Tilia and tiliagraph, PC spreadsheet and graphics software for pollen data, INQUA work. Gr. Data Handl. Methods Newsl.4, 5–7.
57
GroveA. T. (2001). The “Little ice age” and its geomorphological consequences in mediterranean europe. Clim. Change48, 121–136. doi: 10.1023/A:1005610804390
58
HahmW. J.RempeD. M.DralleD. N.DawsonT. E.LovillS. M.BrykA. B.et al. (2019). Lithologically controlled subsurface critical zone thickness and water storage capacity determine regional plant community composition. Water Resour. Res.55. doi: 10.1029/2018WR023760
59
HallettD. J. (2001). Holocene fire history and climate change in southern british columbia, based on high-resolution analyses of sedimentary charcoal (Burnaby (BC): Simon Fraser University).
60
HallettD. J.HillsL. V. (2006). Holocene vegetation dynamics, fire history, lake level and climate change in the Kootenay Valley, Southeastern British Columbia, Canada. J. Paleolim.35, 351–371. doi: 10.1007/s10933-005-1335-6
61
HallettD. J.LepofskyD. S.MathewesR. W.LertzmanK. P. (2003). 11,000 Years of fire history and climate in the mountain hemlock rain forests of southwestern British Columbia based on sedimentary charcoal. Can. J. For. Res.33, 292–312. doi: 10.1139/x02-177
62
HallettD.WalkerR. (2000). Paleoecology and its application to fire and vegetation management in Kootenay National Park, British Columbia. J. Paleolim.24, 401–414. doi: 10.1023/A:1008110804909
63
HebdaR. J. (1983). Late-glacial and postglacial vegetation history at Bear Cove Bog, northeast Vancouver Island, British Columbia. Can. J. Bot.61, 3172–3192. doi: 10.1139/b83-355
64
HebdaR. J. (1995). British Columbia vegetation and climate history with focus on 6 ka BP. Géogr. Phys. Quater.49, 55–79. doi: 10.7202/033030ar
65
HebdaR. J. (1997). “Impact of Climate Change on Biogeoclimatic Zones of British Columbia And Yukon,” in Responding to Global Climate Change in British Columbia and Yukon, vol. 1-15 . Eds. TaylorE.TaylorB. (Ottawa, Ontario: Environment Canada and British Columbia Ministry of Environment, Lands and Parks).
66
HebdaR. J.MathewesR. W. (1984). Holocene history of cedar and native indian cultures of the north american pacific coast. Science225 (4663), 711–713. doi: 10.1126/science.225.4663.711
67
HeusserL. E. (1983). Palynology and paleoecology of post-glacial sediment in an anoxic basin, Saanich Inlet, British Columbia. Can. J. Earth Sci.20, 873–885. doi: 10.1139/e83-077
68
HeusserC. J.HeusserL. E.StreeterS. S. (1980). Quaternary temperatures and precipitation for the north-west coast of North America. Nature286, 702–704. doi: 10.1038/286702a0
69
HigueraP. (2009). CharAnalysis 0.9: Diagnostic and analytical tools for sediment-charcoal analysis. Available online at: https://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.513.64&rep=rep1&type=pdf (Accessed December 1, 2021).
70
HigueraP. E.PetersM. E.BrubakerL. B.GavinD. G. (2007). Understanding the origin and analysis of sediment-charcoal records with a simulation model. Quat. Sci. Rev.26, 1790–1809. doi: 10.1016/j.quascirev.2007.03.010
71
HigueraP. E.GavinD. G.BartleinP. J.HallettD. J. (2010). Peak detection in sediment–charcoal records: impacts of alternative data analysis methods on fire-history interpretations. Int. J. Wild. Fire19, 996–1014. doi: 10.1071/WF09134
72
IpF.LavoieS. (2019). PEOPLE 2020: BC Sub-Provincial Population Projections. Available online at: https://www2.gov.bc.ca/assets/gov/data/statistics/people-population-community/population/people_population_projections_highlights.pdf (Accessed April 2022).
73
IvanovicR.GregoireL.KageyamaM.RocheD. M.ValdesP. J.BurkeA.et al. (2016). Transient climate simulations of the deglaciation 21–9 thousand years before present (version 1) – PMIP4 Core experiment design and boundary conditions. Geosci. Mod. Dev.9, 2563–2587. doi: 10.5194/gmd-9-2563-2016
74
JenningsA. E.HaldM.SmithM.AndrewsJ. T. (2006). Freshwater forcing from the Greenland ice sheet during the Younger Dryas: evidence from southeastern Greenland shelf cores. Quat. Sci. Rev.25, 282–298. doi: 10.1016/j.quascirev.2005.04.006
75
JosenhansH.FedjeD.PienetzR.SouthonJ. (1997). Early humans and rapidly changing Holocene sea levels in the Queen Charlotte Islands-Hecate Strait, British Columbia, Canada. Science277, (71–74). doi: 10.1126/science.277.5322.71
76
KenadyS. M.WilsonM. C.SchalkR. F.MierendorfR. R. (2011). Late Pleistocene butchered Bison antiquus from Ayer Pond, Orcas Island, Pacific Northwest: Age confirmation and taphonomy. Quat. Int.233, 130–141. doi: 10.1016/j.quaint.2010.04.013
77
KennettD. J.KennettJ. P.WestG. J.ErlandsonJ. M.JohnsonJ. R.HendyI. L.et al. (2008). Wildfire and abrupt ecosystem disruption on California’s Northern Channel Islands at the Ållerød–Younger Dryas boundary (13.0–12.9 ka). Quat. Sci. Rev.27, 2530–2545. doi: 10.1016/j.quascirev.2008.09.006
78
KidwellS. M. (2015). Biology in the Anthropocene: Challenges and insights from young fossil records. PNAS112, 4922–4929. doi: 10.1073/pnas.1403660112
79
KienastS. S.McKayJ. L. (2001). Sea surface temperatures in the subarctic northeast Pacific reflect millennial-scale climate oscillations during the last 16 kyrs. Geophys. Res. Let.28, 1563–1566. doi: 10.1029/2000GL012543
80
KrajinaV. J. (1969). “Ecology of forest trees in British Columbia,” in Ecology of western North America. Vol. 1. Ed. V. J. Krajina (University of British Columbia, Department of Botany, Vancouver).
81
KrawchukM. A.MoritzM. A. (2011). Constraints on global fire activity vary across a resource gradient. Ecology92 (1), 121–132. doi: 10.1890/09-1843.1
82
LacourseT. (2005). Late Quaternary dynamics of forest vegetation on northern Vancouver Island, British Columbia. Quat. Sci. Rev.24, 105–121. doi: 10.1016/j.quascirev.2004.05.008
83
LandresP. B.MorganP.SwansonF. J. (1999). Overview of the use of natural variability concepts in managing ecological systems. Ecol. Appl.9 (4), 1179–1188. doi: 10.2307/2641389
84
LepofskyD.BlakeM.BrownD.MorrisonS.OakesN.LyonsN. (2000). The archaeology of the Scowlitz site, SW British Columbia. J. Field Archaeol.27, 391–416. doi: 10.1179/jfa.2000.27.4.391
85
LepofskyD.LertzmanK.HallettD.MathewesR. (2005). Climate change and culture change on the Southern coast of British Columbia 2400-1200 cal. B.P.: Hypothesis. Am. Antiq.70, 267–293. doi: 10.2307/40035704
86
LepofskyD.SchaepeD. M.GraeschA. P.LenertM.OrmerodP.CarlsonK. T.et al. (2009). Exploring Stó:Lō-coast salish interaction and identity in ancient houses and settlements in the Fraser Valley, British Columbia. Am. Antiq.74, 595–626. doi: 10.1017/S0002731600048988
87
LepofskyD.SmithN. F.CardinalN.HarperJ.MorrisM.GitlaE. W.et al. (2015). Ancient shellfish mariculture on the northwest coast of North America. Am. Antiq.80, 236–259. doi: 10.7183/0002-7316.80.2.236
88
LongC. J.WhitlockC.BartleinP. J.MillspaughS. H. (1998). A 9000-year fire history from the Oregon Coast Range, based on a high-resolution charcoal study. Can. J. For. Res.28, 774–787. doi: 10.1139/cjfr-28-5-774
89
LucasJ. D.LacourseT. (2013). Holocene vegetation history and fire regimes of pseudotsuga menziesii forests in the gulf islands national park reserve, southwestern british columbia, canada. Q. Res.79, 366–376. doi: 10.1016/j.yqres.2013.03.001
90
Lynch-StieglitzJ.SchmidtM. W.CurryW. B. (2011). Evidence from the Florida Straits for Younger Dryas ocean circulation changes. Paleoceanography26 (1), PA1205. doi: 10.1029/2010PA002032
91
MannM. E.ZhangZ.RutherfordS.BradleyR. S.HughesM. K.ShindellD.et al. (2009). Global signatures and dynamical origins of the little ice age and medieval climate anomaly. Science326, 1256–1260. doi: 10.1126/science.1177303
92
MathewesR. W.KingM. (1989). Holocene vegetation, climate, and lake-level changes in the interior douglas-fir biogeoclimatic zone, british columbia. Can. J. Earth Sci.26 (9), 1811–1825. doi: 10.1139/e89-154
93
MathewesR. W. (1973). A palynological study of postglacial vegetation changes in the University Research Forest, southwestern British Columbia. Can. J. Bot.51, 2085–2103. doi: 10.1139/b73-271
94
MathewesR. W. (1993). Evidence for Younger Dryas-age cooling on the North Pacific coast of America. Quat. Sci. Rev.12, 321–331. doi: 10.1016/0277-3791(93)90040-S
95
MathewesR. W.HeusserL. E. (1981). A 12,000 year palynological record of temperature and precipitation trends in southwestern British Columbia. Can. J. Bot.59, 707–710. doi: 10.1139/b81-100
96
MatsonR. G.CouplandG. (1995). The Prehistory of the Northwest Coast (San Diego, CA: Academic Press), 364.
97
McCormickE. L.DralleD. N.HahmW. J.TuneA. K.SchmidtL. M.ChadwickK. D.et al. (2021). Widespread woody plant use of water stored in bedrock. Nature597, 225–229. doi: 10.1038/s41586-021-03761-3
98
McCuneJ. L.PellattM. G.VellendM. (2013). Multidisciplinary synthesis of long-term human–ecosystem interactions: a perspective from the garry oak ecosystem of british columbia. Biol. Conserv.166, 293–300. doi: 10.1016/j.biocon.2013.08.004
99
McLachlanJ. S.BrubakerL. B. (1995). Local and regional vegetation change on the northeastern olympic peninsula during the holocene. Can. J. Bot.73 (10), 1618–1627. doi: 10.1139/b95-175
100
McLarenD. (2017). “The occupational history of the stave watershed,” inArchaeology of the lower fraser river region. Ed. RousseauM. K. (Simon Fraser University, Burnaby, BC: Archaeology Press), 149–158.
101
McLarenD.FedjeD.DyckA.MackieQ.GauvreauA.CohenJ. (2018). Terminal Pleistocene epoch human footprints from the Pacific coast of Canada. PloS One13, e0193522. doi: 10.1371/journal.pone.0193522
102
McManusJ. F.FrancoisR.GherardiJ. M.KeigwinL. D.Brown-LegerS. (2004). Collapse and rapid resumption of atlantic meridional circulation linked to deglacial climate changes. Nature428, 834–837. doi: 10.1038/nature02494
103
MeidingerD.PojarJ. (Eds.) (1991). Ecosystems of British Columbia (Victoria, BC: Ministry of Forests and Range Research Branch).
104
MenounosB.OsbornG.ClagueJ. J.LuckmanB. H. (2009). Latest Pleistocene and Holocene glacier fluctuations in western Canada. Quat. Sci. Rev.28, 2049–2074. doi: 10.1016/j.quascirev.2008.10.018
105
Metro Vancouver (2016). Climate Projections for Metro Vancouver. Available online at: https://metrovancouver.org/services/air-quality-climate-action/Documents/climate-projections-for-metro-vancouver-2016.pdf (Accessed April 1, 2022).
106
Metro Vancouver (2018). Regional parks land acquisition 2050. (Metro Vancouver Regional District). Available at: http://www.metrovancouver.org/services/parks/ParksPublications/RegionalParksLandAcquisition2050.pdf (Accessed May 15, 2022).
107
MooreP. D.WebbJ. A.CollinsonM. E. (1991). Pollen Analysis. 2nd ed (Oxford: Blackwell Science).
108
MoteP. W.SalathéE. P. (2010). Future climate in the Pacific Northwest. Climatic Change102, 29–50. doi: 10.1007/s10584-010-9848-z
109
MoyC. M.SeltzerG. O.RodbellD. T.AndersonD. M. (2002). Variability of El Niño/Southern Oscillation activity at millennial timescales during the Holocene epoch. Nature420, 162–165. doi: 10.1038/nature01194
110
MurphyS. F.PellattM. G.KohfeldK. E. (2019). A 5,000-year fire history in the strait of Georgia Lowlands, British Columbia, Canada. Front. Ecol. Evol. 7. doi: 10.3389/fevo.2019.00090
111
NASA (2024). World of Change: Global Temperatures. Available online at: https://earthobservatory.nasa.gov/world-of-change/global-temperatures:~:text=According%20to%20an%20ongoing%20temperature,1.9%C2%B0%20Fahrenheit)%20since%201880 (Accessed December 30, 2024).
112
North Greenland Ice Core Project Members (2007). 50 year means of oxygen isotope data from ice core NGRIP. Supplement to: North Greenland Ice Core Project Members, (2004): High-resolution record of Northern Hemisphere climate extending into the last interglacial period. Nature431, 147–151. doi: 10.1038/nature02805
113
Pacific Climate Impacts Consortium (2024). Available online at: https://www.pacificclimate.org/analysis-tools/plan2adapt (Accessed December 20, 2024).
114
PalmerS.WalkerI.HeinrichsM.HebdaR.ScudderG. (2002). Postglacial midge community change and holocene palaeotemperature reconstructions near treeline, southern british columbia (Canada). J. Paleolimnol. 28, 469–490. doi: 10.1023/A:1021644122727
115
PellattM. G.HebdaR. J.MathewesR. W. (2001a). High-resolution Holocene vegetation history and climate from Hole 1034B, ODP leg 169s, Saanich Inlet, Canada. Mar. Geol.174, 211–226. doi: 10.1016/S0025-3227(00)00151-1
116
PellattM. G.GedalofZ. (2014). Environmental change in garry oak (Quercus garryana) ecosystems: the evolution of an eco-cultural landscape. Biodivers. Conserv.23, 2053–2067. doi: 10.1007/s10531-014-0703-9
117
PellattM. G.SmithM. J.MathewesR. W.WalkerI. R.PalmerS. L. (2000). Holocene treeline and climate change in the subalpine zone near stoyoma mountain, cascade mountains, southwestern british columbia, canada. Arctic Antarctic Alpine Res.32 (1), 73–83. doi: 10.2307/1552412
118
PellattM. G.MathewesR. W. (1994). Paleoecology of postglacial tree line fluctuations on the Queen Charlotte Islands, Canada. Écoscience1, 71–81. doi: 10.1080/11956860.1994.11682230
119
PellattM. G.MathewesR. W. (1997). Holocene tree line and climate change on the Queen Charlotte Islands, Canada. Quat. Res.48, 88–99. doi: 10.1006/qres.1997.1903
120
PellattM. G.MathewesR. W.ClagueJ. J. (2001b). Implications of a late-glacial pollen record for the glacial and climatic history of the Fraser Lowland, British Columbia. Paleogeog. Paleoclim. Paleoecol.180, 147–157. doi: 10.1016/S0031-0182(01)00426-6
121
PellattM. G.McCoyM. M.MathewesR. W. (2015). Paleoecology and fire history of Garry oak ecosystems in Canada: implications for conservation and environmental management. Biodiv. Cons.24, 1621–1639. doi: 10.1007/s10531-015-0880-1
122
PrichardS. J.GedalofZ.OswaldW. W.PetersonD. L. (2009). Holocene fire and vegetation dynamics in a montane forest, north cascade range, washington, USA. Q. Res.72, 57–67. doi: 10.1016/j.yqres.2009.03.008
123
ReasonerM. A.JodryM. A. (2000). Rapid response of alpine timberline vegetation to the younger dryas climate oscillation in the colorado rocky mountains, USA. Geology28 (1), 51–54. doi: 10.1130/0091-7613(2000)028<0051:RROATV>2.3.CO;2
124
ReasonerM. A.RutterN. W.OsbornG. (1994). Age of the crowfoot advance in the canadian rocky mountains. a glacial event coeval with the younger dryas oscillation. Geology.22 (5), 439–442. doi: 10.1130/0091-7613(1994)022<0439:AOTCAI>2.3.CO;2
125
ReimerP. J.AustinW. E. N.BardE.BaylissA.BlackwellP. G.RamseyC. B.et al. (2020). The Intcal20 northern hemisphere radiocarbon age calibration curve (0-55 cal kBP). Radiocarbon62 (4), 725–757. doi: 10.1017/RDC.2020.41
126
RenssenH.MairesseA.GoosseH.MathiotP.HeiriO.RocheD.et al. (2015). Multiple causes of the younger dryas cold period. Nat. Geosci.8, 946–949. doi: 10.1038/ngeo2557
127
RenssenH.GoosseH.FichefetT.BrovkinV.EsschaertE.WolkF. (2005). Simulating the Holocene climate evolution at northern high latitudes using a coupled atmosphere-sea ice-ocean-vegetation model. Clim. Dy.24, 23–43. doi: 10.1007/s00382-004-0485-y
128
RindD.PeteetD.BroeckerW.McIntyreA.RuddimanW. (1986). The impact of cold North Atlantic sea surface temperatures on climate: implications for the Younger Dryas cooling (11-10 k). Clim. Dy.1, 3–33. doi: 10.1007/BF01277044
129
RitchieM.LepofskyD. (2020). From local to regional and back again: Social transformation in a Coast Salish settlement 1500–1000 BP. J. Anth. Arch.60, 101210. doi: 10.1016/j.jaa.2020.101210
130
RorigM. L.FergusonS. A. (1999). Characteristics of lightning and wildland fire ignition in the Pacific Northwest. J. App. Meteorol. (1988)38, 1565–1575. doi: 10.1175/1520-0450(1999)038<1565:COLAWF>2.0.CO;2
131
RosenbergC. M.WalkerI. R.MathewesR. W.HallettD. J. (2004). Midge-inferred Holocene climate history of two subalpine lakes in southern British Columbia, Canada. Holocene14 (2), 258–271. doi: 10.1191/0959683604hl703rp
132
RyderJ. M.ThomsonB. (1986). Neoglaciation in the southern Coast Mountains of British Columbia: chronology prior to the late Neoglacial maximum. Can. J. Ear. Sci.23, 273–287. doi: 10.1139/e86-031
133
SankeyJ. B.KreitlerJ.HawbakerT. J.McVayJ. L.MillerM. E.MuellerE. R.et al. (2017). Climate, wildfire, and erosion ensemble foretells more sediment in western USA watersheds. Geophys. Res. Let.44, 8884–8892. doi: 10.1002/2017GL073979
134
SchaepeD. M. (1998). Recycling Archaeology: Analysis of material from the 1973 excavation of an ancient house at the Mauer Site. Simon Fraser University, Burnaby (BC).
135
SchaepeD. M. (2009). Pre-colonial stó:lō Community organization: an archaeological perspective. University of British Columbia, Vancouver (BC).
136
SheaC. J.SteinmanB. A.BrownE. T.SchreinerK. M. (2022). A multi-proxy lake-sediment record of middle through late Holocene hydroclimate change in southern British Columbia, Canada. J. Paleolim.67, 163–182. doi: 10.1007/s10933-021-00231-8
137
StuiverM.ReimerP. (1993). Extended 14C database and revised CALIB radiocarbon calibration program. Radiocarbon35, 215–230. doi: 10.1017/S0033822200013904
138
SugimuraW. Y.SprugelD. G.BrubakerL. B.HigueraP. E. (2008). Millennial-scale changes in local vegetation and fire regimes on mount constitution, orcas island, washington, USA, using small hollow sediments. Can. J. For. Res.38, 539–552. doi: 10.1139/X07-186
139
SwetnamT. W.AllenC. D.BetancourtJ. L. (1999). Applied historical ecology: using the past to manage for the future. Ecol. App.9, 1189–1206. doi: 10.1890/1051-0761(1999)009[1189:AHEUTP]2.0.CO;2
140
TurnerN. (2017). “Re Styecwmenúlecws-kucw: How We Look(ed) after Our Land,” in Secwépemc People, Land, and Laws: Yerí7 re Stsq’ey’s-kucw. Eds. IgnaceM.IgnaceR. E. (McGill-Queen’s University Press, Montreal & Kingston), 145–219.
141
TurnerN.PeacockS. (2005). “Solving the Perennial Paradox: Ethnobotanical Evidence for Plant Resource Management on the Northwest Coast,” in Keeping it Living: Traditions of Plant Use and Cultivation on the Northwest Coast of North America. Eds. DeurD.TurnerN. (UBC Press, Vancouver), 101–150.
142
WainmanN.MathewesR. W. (1987). Forest history of the last 12000 years based on plant macrofossil analysis of sediment from Marion Lake, southwestern British Columbia. Can. J. Bot.65, 2179–2187. doi: 10.1139/b87-300
143
WalkerI. R.PellattM. G. (2003). Climate change in coastal British Columbia - A paleoenvironmental perspective. Can. Wat. Res. J.28, 531–566. doi: 10.4296/cwrj2804531
144
WalshM. K.MarlonJ. R.GoringS. J.BrownK. J.GavinD. G. (2015). A regional perspective on holocene fire–climate–human interactions in the Pacific Northwest of North America. Ann. Ass. Am. Geog.105, 1135–1157. doi: 10.1080/00045608.2015.1064457
145
WangT.HamannA.SpittlehouseD.CarrollC. (2016). Locally downscaled and spatially customizable climate data for historical and future periods for North America. PloS One11, e0156720. doi: 10.1371/journal.pone.0156720
146
WhitlockC.McWethyD. B.TepleyA. J.VeblenT. T.HolzA.McGloneM. S.et al. (2015). Past and present vulnerability of closed-canopy temperate forests to altered fire regimes. Bioscience65, 151–163. doi: 10.1093/biosci/biu194
147
WilliamsJ. (2006). Clam gardens: Aboriginal mariculture on Canada’s west coast (Vancouver: New Star Books).
148
WrightH. E.MannD. H.GlaserP. H. (1984). Piston corers for peat and lake sediments. Ecology65, 657–659. doi: 10.2307/1941430
149
ZdanowiczC. M.ZielinskiG. A.GermaniM. S. (1999). Mount Mazama eruption; calendrical age verified and atmospheric impact assessed. Geology27, 621–624. doi: 10.1130/0091-7613(1999)027<0621:MMECAV>2.3.CO;2
150
ZwolinksiM. J. (2000). “Role of Fire in Management of Watershed Responses,” in Land Stewardship in the 21st Century: The Contributions of Watershed Management (U.S. Dept. of Agriculture, Forest Service, Rocky Mountain Research Station, Fort Collins, CO), 367–370.
Summary
Keywords
paleoclimate, paleoecology, wildfire, disturbance, novel ecosystems
Citation
Duncan ME, Pellatt MG and Kohfeld KE (2025) Coastal moderation of Holocene fire and vegetation change on the Pacific coast of Canada. Front. Ecol. Evol. 13:1504983. doi: 10.3389/fevo.2025.1504983
Received
01 October 2024
Accepted
16 January 2025
Published
10 February 2025
Volume
13 - 2025
Edited by
Manel Leira, University of Santiago de Compostela, Spain
Reviewed by
Scott David Mooney, University of New South Wales, Australia
Maarten Blaauw, Queen’s University Belfast, United Kingdom
Updates

Check for updates
Copyright
© 2025 Duncan, Pellatt and Kohfeld.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Maggie E. Duncan, maggie.duncan@dri.edu
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.