Forest Transformation Following European Settlement in the Saguenay-Lac-St-Jean Valley in Eastern Québec, Canada

Introduction

Forests of northeastern North America have undergone major transformations since European settlement. During the 19th and 20th centuries, deforestation for agriculture, fire escapes from settlements and subsequent fire suppression, industrial logging and urbanization have modified the natural disturbance regimes (Whitney, 1994; Lorimer, 2001; Nowacki and Abrams, 2008). As a result, modern forest landscapes are younger, more fragmented and homogenized, and composed of an increased proportion of early-successional species compare to pre-settlement times (Foster et al., 1998; Schulte et al., 2007; Boucher et al., 2009; Thompson et al., 2013; Terrail et al., 2019). Although climate change has also influenced forest composition during that period (Fisichelli et al., 2014; Pederson et al., 2015; Boisvert-Marsh et al., 2019), the impact of anthropogenic activities was stronger or more prevalent in several regions (Nowacki and Abrams, 2015; Danneyrolles et al., 2019), such as the temperate-boreal ecotone (Brice et al., 2019, 2020).

Knowledge of forest composition and natural disturbances prior to European settlement can provide targets for ecosystem-based management and restoration of forest ecosystems (Landres et al., 1999; Gauthier et al., 2009). This is particularly important in heavily transformed regions where unmanaged forests no longer exist. While our ability to reduce gaps between pre-settlement and modern forest composition could become unrealistic under severe climate forcing (Boulanger et al., 2019), understanding how forests have been transformed by human activities can help make forests more resilient to future disturbances and climate change (Millar and Stephenson, 2015; Johnstone et al., 2016) as well as to restore functional targets. For example, growth rate, longevity and shade tolerance of pre-settlement species could be considered, regardless of strict species identity, in order to maintain or restore ecosystem function.

Land survey records are a reliable source of data to evaluate changes in forest composition since pre-settlement (Schulte and Mladenoff, 2001; Terrail et al., 2014). Several studies have documented changes in frequency and dominance of individual taxa at different spatial scales by comparisons with modern inventories (e.g., Friedman and Reich, 2005; Dupuis et al., 2011; Thompson et al., 2013), and have shown that taxa assemblages also responded to environmental gradients (e.g., elevation, surface deposits) and disturbances (White and Mladenoff, 1994; Schulte et al., 2007; Danneyrolles et al., 2016, 2018; Terrail et al., 2019; Thomas-Van Gundy et al., 2020).

In the northern part of the temperate zone, forests were historically dominated by late-successional and long-lived conifers, particularly spruces (Picea glauca and Picea mariana), balsam fir (Abies balsamea), white cedar (Thuja occidentalis), and pines (Pinus strobus and Pinus resinosa), along with deciduous companion species like maples (Acer spp.) and yellow birch (Betula alleghaniensis) (Lorimer, 1977; Jackson et al., 2000; Dupuis et al., 2011; Terrail et al., 2019). Sectors where fires were more frequent included a larger component of pines (notably jack pine, Pinus banksiana), poplars (Populus tremuloïdes and Populus balsamifera) and white birch (Betula papyrifera) (Pinto et al., 2008; Hanberry et al., 2013; Danneyrolles et al., 2016; Marchais et al., 2020). The most noticeable vegetation change since pre-settlement is the general expansion of early-successional and fast-growing species such as poplars, maples, and white birch to the detriment of conifers (Friedman and Reich, 2005; Pinto et al., 2008; Danneyrolles et al., 2016; Terrail et al., 2019; Marchais et al., 2020). However, only a few studies have used land survey records to evaluate forest composition changes in the temperate-boreal ecotone in Ontario (Pinto et al., 2008) and western Québec (Danneyrolles et al., 2016; Marchais et al., 2020), and none in eastern Québec.

Located at the ecotone between temperate and boreal forests in eastern Québec, the Saguenay-Lac-St-Jean valley (hereafter SLSTJ) has experienced rapid colonization since the mid-19th century, resulting in the clearing of large areas of forest for agriculture in parallel with numerous human-induced fires (Girard and Perron, 1989; Boucher et al., 2014). The region has also been extensively exploited by the forest industry, such that unmanaged forests are now small and scarce. While fire had a relatively minor role in disturbance regime in temperate forests of eastern Québec and northeastern United States just before European settlement (Lorimer, 1977; Fortin, 2018; Terrail et al., 2019), the SLSTJ probably experienced higher fire activity due to its location at the very southern fringe of the boreal biome, which is characterized by higher fire activity than the temperate forest (Erni et al., 2020).

Although land survey archives are available throughout SLSTJ in the form of tree taxon lists, which are ranked tree taxa enumerations along survey lines (Dupuis et al., 2011; Terrail et al., 2014), this data source has not been exploited yet. Accordingly, the objectives of this study were to (1) document the pre-settlement forest composition of SLSTJ at both regional and landscape scales, (2) assess changes through comparison with modern forest inventories, while paying special attention to settlement fires and environmental gradients, and (3) discuss how this new knowledge can be used to plan forest management in the temperate-boreal ecotone.

Study Area

The study area (19,400 km²) is located in eastern Québec, between 69°18′ and 72°59′W and 47°47′ and 49°15′N (Figure 1). It is situated within the Saguenay-Lac-St-Jean valley, a depression in the eastern rim of the Canadian Shield (Robitaille and Saucier, 1998). According to the provincial vegetation classification, the Saguenay Lac-St-Jean plains and the low hills surrounding both Saguenay and St. Lawrence Rivers correspond to mixed temperate forests, whereas boreal mixedwood forests dominate the more elevated outer edge of the study area (Robitaille and Saucier, 1998). The natural disturbance regime is characterized by low intensity outbreaks of spruce budworm (Choristoneura fumiferana) (Blais, 1983; Simard et al., 2006), and fires (Boucher et al., 2014; Pilon and Payette, 2015). Fertile agricultural lands are present in the lowlands where clay deposits left by the postglacial Laflamme Sea (Elson, 1969) mix with sandy and organic deposits, while glacial tills and rocky outcrops are the main surface deposits on the uplands (Robitaille and Saucier, 1998). Annual mean temperature and precipitations (1981–2010) range from 0.8 to 3.4°C and from 833 to 1 180 mm (with 24–35% as snow), respectively (Environment Climate Change Canada, 2019).

Figure 1. Location of the study area (inset) in eastern Canada, along with maps of townships and vegetation zones (A), and of the elevation bands (metres above sea level) showing the Saguenay-Lac-St-Jean Plains (B).

From the 17th to the early 19th centuries, human activities in the Saguenay Lac-St-Jean valley were related to the fur trade between the hunter-gatherer Innu First Nation and Europeans (Girard and Perron, 1989). Nomadic people have frequented the valley for more than 6,000 years (Langevin and Plourde, 2018) and used its hydrographic network to travel and make exchanges with other First Nations (Girard and Perron, 1989). Even though they may have influenced the fire regime, the land-use practices of these populations and their impact on the regional forest are poorly documented (Oberndorfer, 2020). It was only in 1842 that the region was officially opened to European settlement with the aim of clearing new agricultural land and exploiting forests to supply timber to England (Girard and Perron, 1989). Selective cuttings during the 19th century were directed toward large white pine, red pine and white spruce trees along rivers. Extensive logging began at the turn of the 20th century with the rise of pulp and paper industry targeting conifers (Girard and Perron, 1989). The important river network allowed the construction of several hydroelectric dams and the region experienced an industrial golden age during the 20th century (Girard and Perron, 1989), with the population increasing about tenfold from about 28,800 to 262,500 inhabitants between 1891 and 1961 (Pouyez et al., 1983).

Materials and Methods

Historical and Modern Data

The reconstruction of pre-settlement forest composition is based on 146 logbooks recording the first subdivision of 74 townships from 1842 to 1935, of which 80% were surveyed before 1900. The typical SLSTJ township is a square of 10 × 10 miles (16 km × 16 km) subdivided into parallel ranges of 1 mile (1.6 km) deep, that are in turn subdivided into lots of 13 chains (1 chain = 20.1 m) wide. Observations of forest composition along range lines are mostly taxon lists (e.g., spruce, fir, and white birch) and specific cover types (e.g., spruce stand and cedar stand). General cover types such as “mixed wood” or “hardwood” are infrequent and were not considered in the analyses. We assigned a rank to each taxon listed, assuming that taxa were enumerated according to their relative basal area in the forest stands (Terrail et al., 2014). Specific cover types were considered equivalent to pure stands of the corresponding taxon in order to include them as taxon lists.

In the surveyors’ logbooks, each observation comes with a distance in chains from the start of the survey line, or from the last planted post. We classified these observations into two geometric types, points or lines, according to how the surveyors recorded them. Most observations are points, written next to a distance but with no clear beginning or end along the survey line. Points are either regularly or irregularly spaced. Regularly spaced point observations are often made at lot corners or at every 5, 10, or 20 chains (101, 201, and 402 m). Linear observations are also frequent and contain both a start and an end position, generally delineated by lot limits (260 m). In order to combine points and lines into the same database, both types were weighted according to the size of the corresponding stands along surveyed lines. Each line observation was weighted by its length and each point was weighted by the mean of the distances to the previous and next observations (Dupuis et al., 2011). Weighted observations were used in all calculations involving historical data. Observations were georeferenced as points or lines using modern cadastral maps built from the same land survey records.

Because the occurrence of fires was important in the logbooks reporting the original surveys of townships (see section “Results”), fires mentioned in resurvey mission logbooks were also tallied in order to document the prevalence of forest fires over the entire settlement period. Resurvey missions were required to inspect or correct previous works, and to restore original lines and lot corners that had disappeared over time. To complete the picture, recent fires (1950–2015) were compiled from the Québec government fire database (MFFP).

The present-day forest composition was assessed using the last three decadal forest inventories (1980–2010) conducted by the Québec government. These inventories are based on circular plots of 0.04 ha, stratified according to the different types of productive forests (Ministère des Forêts, de la Faune et des Parcs secteur des forêts [MFFP], 2016). Within the plots, all stems >9 cm DBH (diameter at breast height: 1.3 m) are tallied by species and 2 cm DBH classes, thus allowing calculation of the basal area (m² ha^–1) by species. To match the taxa mentioned by surveyors, spruce, maple, poplar, and pine species other than jack pine (listed as “cypress” by surveyors) were grouped at the genus level (i.e., Picea spp., Acer spp., Populus spp., and other pines). In modern plots, all taxa are recorded in all plots, whereas in historical observations less abundant taxa may have been overlooked by surveyors at some sites. Thus, taxa representing less than 5% of the total basal area of a plot were excluded from the modern dataset so that the average number of taxa per observation is equivalent between the two periods (Supplementary Figure S1). A rank was then assigned to each taxon according to its relative basal area within modern plots. Only modern plots that were located within a maximum distance of 3 km from the nearest historical observation were retained. In sectors with low density of historical observations, this maximum distance was set to 1 km. Finally, taxa that occurred in less than 3% of both data sets were grouped as “other.”

Data Analyses

We used three different metrics to describe the regional forest composition for the pre-settlement and modern periods. First, we calculated the frequency of each taxon as the proportion of all observations (weighted length or modern plots) mentioning the corresponding taxon, regardless of its rank in those observations. Second, we calculated a rank frequency (RF_ir) for each taxon occurring in the first four ranking positions (i.e., r = 1, 2, 3, 4; Scull and Richardson, 2007) using the formula:

where N_i is the number of observations mentioning the taxon i and M_i is the total number of observations. Third, the dominance represents the frequency of occurrence of each taxon i at the first ranking position (i.e., r = 1).

To compare pre-settlement and present-day forest composition at the landscape scale, we divided the study area into a grid of 5 km × 5 km cells. We first retained cells containing at least four historical observations (n = 579), of which we discarded cells containing less than four modern plots (99 cells or 17%, Supplementary Figure S2), resulting in 480 cells with a mean of 42 historical observations and 24 modern plots. Most discarded cells on the plains did not contain enough modern plots due to settlement and low present-day forest cover (Supplementary Figure S2), while the low density of historical observations at the periphery of the study area was the main factor restricting the retention of cells with enough modern plots.

In each cell, and for both periods, we first calculated the frequency of each taxon as the proportion of all observations (weighted length or modern plots) mentioning the corresponding taxon in the grid cell. Then, to facilitate comparisons between taxa and cells, we calculated a relative frequency for each taxon into each cell by scaling its frequency relative to the sum of the frequencies of all taxa. We mapped the relative frequency for each of the 11 retained taxa.

Two subsequent analyses were performed using the grid cells of taxa relative frequency. We first identified landscape tree assemblages for both time periods using a Ward-like hierarchical clustering based on Bray–Curtis dissimilarity index. This clustering method minimizes the number of groups and the intra-group variation (Legendre and Legendre, 2012). The 480 cells from the pre-settlement and modern periods were combined into one dataset and Bray–Curtis dissimilarity was calculated between every pair of cells based on their taxa relative frequency. The clustering was performed in a single run, so that the outcome was one set of clusters for the 960 cells, which allowed a direct comparison of tree assemblages between periods. An optimal number of five clusters was retained in order to minimize the number of clusters while maximizing the proportion of variance explained by the clustering (Supplementary Figure S3).

Second, to visualize compositional changes along environmental gradients, a principal component analysis (PCA) was performed on environmental variables (surface deposits including clay, organic, till, rocks, sand; elevation; slope) for the 480 cells. We also considered a “settlement” variable as a combination of agricultural and urban (including roads) areas. All variables are expressed as percentages of land cover per cell as determined from modern forestry maps (photo-interpretation 2006–2011; 1:20,000 scale; Ministère des Forêts, de la Faune et des Parcs secteur des forêts [MFFP], 2015), except elevation and slope which are means from a TIN (Triangular Irregular Network) built from governmental hypsometric maps (scale 1:20 000 with 10-m contour intervals; Ministère de l’Énergie et des Ressources Naturelles du Québec [MERN], 2007). The PCA was performed on a correlation matrix of the standardized environmental variables. No axes rotation was used to facilitate the interpretation of the data. In order to evaluate the transition of landscape tree assemblages across environmental gradients, each cluster for each period was plotted as the centroid of cell scores of the corresponding cluster on the two principal axes of the PCA. PCA biplots showing environmental variables and cell colors for the clusters of each period are presented in Supplementary Figure S4. Ward’s hierarchical clustering and PCA analysis were performed using R software (R Core Team, 2019 version 3.6.1) using the stats (hclust function) and FactoMineR (PCA function) packages, respectively.

Results

Fires occurred throughout the study area during the settlement period (Figure 2). Surveyors reported at least one fire in 72 of the 74 townships under study between 1842 and 1971. Fires from resurveys represent 58% of all fires (2,164 km), and most of them happened during the 1920 and 1930 decades (Figure 2A). Sectors with no fire reported by surveyors could still have burned because fire compilation maps prior to 1950 (Québec government, not shown), which are also incomplete, show several burned areas not mentioned by surveyors. Concerning recent fires, they mostly occurred during the 1950 decade, which account for 79% of the total area burned between 1950 and 2015 (Figure 2B). Recent fires were mostly concentrated in the northern sector where they overlapped several fires previously recorded by surveyors. The cause of the fires during the 1950–1970 period is unknown, while 90% of the fires reported between 1970 and 2015 were of human origin.

Figure 2. Maps showing (A) surveyor’s fire observations (points and lines) tallied from logbooks of original surveys and resurveys, (B) distribution and size of recent fires from mapped fire polygons, and (C) present-day land use in the study area. Histograms represent kilometers of fire along range lines (surveyor’s observations) and burned areas (recent fires) by decade.

During the pre-settlement period, spruce, white birch, and balsam fir were the most frequent (58.6–80.4%) and dominant (12.7–41.5%) taxa (Table 1), across the study area (Figure 3). Yellow birch, poplar and jack pine were also common with frequency ranging from 21.7 to 17.0%. Poplar displayed a patchy distribution across the grid cells, while jack pine and yellow birch were more prevalent in the northwest and southeast sections, respectively. All remaining taxa were comparatively less frequent and dominant and more locally distributed. The spruce taxon was more frequently mentioned at rank 1 than ranks 2–4, indicating that it tended to be dominant when present. Conversely, balsam fir and white birch tended to be co-dominant as they were mostly listed at positions 2 and 3, respectively (Table 1).

Table 1. Taxa frequencies and ranked frequency index (RF_ir) for pre-settlement (1842–1935) and modern (1980–2010) periods.

Figure 3. Maps of taxa relative frequency aggregated to 25-km² grid cells for pre-settlement (1842–1935) and modern (1980–2010) periods. The difference is calculated by subtracting the relative frequency of the two periods (modern minus pre-settlement).

Comparison with modern plots indicated substantial changes in forest composition since pre-settlement time (Table 1 and Figure 3). Poplar and maple showed the most important frequency increases (+33.3 and +14.2%, respectively), while spruce (−23.1%) and yellow birch (−13.8%) experienced the largest decreases. Concerning dominance, poplar showed by far the largest increase (+20.7%), alongside jack pine (+3.6%), maple (+3.2%) and white birch (+2.3%), while spruce and yellow birch declined by 22.1 and 2.9%, respectively. The remaining conifers decreased in frequency and dominance, especially red and white pines, two species that almost disappeared from the region (Table 1). At the cell scale (Figure 3), pattern of changes in relative frequency were rather uniform for most increasing (jack pine, maple, and poplar) or decreasing taxa (spruce, yellow birch, white cedar, larch, and other pines), except white birch and balsam fir which displayed patchy patterns of increases and decreases.

Landscape Tree Assemblages

Cluster analysis revealed considerable reorganization of tree assemblages across the study area during the last century (Figure 4). Clusters 1, 2, and 4 were the most widespread before settlement (83% of all cells) and were dominated by spruce, balsam fir and white birch (Figure 4A). Although similar in overall composition, these three clusters differ from each other on some points. In cluster 2, spruce, balsam fir and white birch dominate largely over the remaining taxa. Cluster 1 shows a more diverse composition and is the cluster with the largest mean relative frequency of yellow birch (16%), white cedar (7%), and other pines (5%) (Figure 4B). Cluster 4 contains noticeably more poplar and jack pine than clusters 1 and 2. In contrast to pre-settlement times, the present-day landscape is dominated by cluster 3 (36% of all cells) with high relative frequency of white birch (24%), poplar (19%), and maple (8%), and the lowest relative frequency of spruce (16%) among clusters. This cluster was almost absent during pre-settlement times (14 cells) and became widespread throughout the study area. The jack pine cluster (5) also increased in frequency, mostly on the plains, while clusters 1 and 2 have been completely excluded from the plains and confined to uplands at the periphery of the study area (Figure 4A).

Figure 4. Landscape tree assemblages identified by hierarchical clustering of the 25-km² grid cells. (A) Maps of the clusters for the two time periods with histograms showing the number of cells in each cluster; (B) Table of cluster composition showing taxon mean relative frequency and standard deviation (bars); (C) Transition among clusters from pre-settlement (top) to modern (bottom) periods. The width of each box is proportional to the corresponding number of cells.

The transition diagram (Figure 4C) and the PCA (Figure 5) provide insights into trajectories of tree assemblages between time periods. Most cells belonging to clusters 1 and 2 in pre-settlement times shifted to clusters 3 and 4 in the present-day, which are characterized by high poplar relative frequency. Cluster 3 has remained ubiquitous across environmental gradients, although slightly more associated with the uplands today. Conversely, cells from clusters 1 and 2 are now much more associated with till deposits, steep slopes and high elevations than they were during pre-settlement times. The jack pine cluster gained cells almost exclusively from cluster 4, where jack pine was already present (Figure 4B), and strengthened its relationship with sandy deposits. In contrast, cluster 4 exchanged many cells with other clusters and became less strongly associated with the surface deposits of the plains.

Figure 5. Principal component analysis (PCA) where axes 1 and 2 explain 46.8% and 19.6% of the environmental variability, respectively. Biplot shows cell scores (grey dots), correlations of environmental variables with axes 1 and 2, and clusters (circles) for both time periods as the centroids of the corresponding cell scores. Circles size is proportional to the corresponding number of cells in each cluster. Arrows between clusters show the transition from pre-settlement to modern periods. PCA1 mainly represents the plains on the left, where clay and sandy deposits are predominant, and the uplands on the right, where glacial till is abundant and settlement is absent. PCA2 is more associated with soil texture and reflects more subtle differences among cells.

Discussion

Pre-settlement Tree Assemblages

The landscape tree assemblage dominated by spruce, balsam fir and white birch (cluster 2) was the most common before settlement and was mostly associated with the uplands of the boreal zone. Since the mean relative frequency was high for balsam fir, and very low for poplar and jack pine, the disturbance regime governing this assemblage was probably characterized by forest openings created by windthrow and spruce budworm outbreaks (Kneeshaw and Bergeron, 1999; Bergeron, 2000; Lesieur et al., 2002), along with occasional fires. Recent research in the region supports our results and shows low fire activity on uplands during the preindustrial era (Couillard et al., 2013; Boucher et al., 2014). Southeast of Lake St-Jean and along the Saguenay River, the spruce-fir-white birch assemblage was often accompanied by more temperate species, such as yellow birch, white cedar and red and white pines (cluster 1). This was probably due to the more temperate conditions associated with the Saguenay River valley, where low elevations create milder winters and warmer summers compared to uplands (Girard and Perron, 1989; Audet et al., 2012).

Two fire-adapted tree assemblages (clusters 4 and 5) prevailed north and west of lake St-Jean. The presence of trembling aspen (Populous tremuloïdes, Supplementary Figure S5) and jack pine in these clusters suggests that fires were frequent before settlement. Fire was probably more influential in the jack pine cluster (5) considering the very low relative frequency of balsam fir, a fire-sensitive species (Frank, 1990), along with the high relative frequency of jack pine, often considered as fire-dependant (Burns and Honkala, 1990). Moreover, the fire-adapted black spruce (Viereck and Johnston, 1990) was probably an important species in cluster 5, as modern plots show a very low present-day frequency of white spruce compared to black spruce (Supplementary Figure S5 and Supplementary Table S1). This sector has few firebreaks as its topography is flat and lake abundance is low. This may have facilitated the spread of fires (Lynch et al., 2014), thus promoting the presence of jack pine, especially on poor and well-drained sandy deposits where fires are more likely to occur (Mansuy et al., 2010).

Innu people may also have influenced fire activity in the SLSTJ region, as did other First Nations in southern Québec (Blarquez et al., 2018), Labrador (Oberndorfer, 2020), Ontario (Munoz and Gajewski, 2010), and United States (Roos et al., 2018; Abrams and Nowacki, 2019). For example, Innu people possibly used fire to manage wild animal populations, to clear small patches to set up camps, to stimulate productivity of berry patches or the recruitment of white birch, which was essential for the construction of dwellings and canoes (Day, 1953; Lewis, 1980; Lewis and Ferguson, 1988; Girard and Perron, 1989; Williams, 2000; Oberndorfer, 2020). Although mostly of local importance (Russel, 1983; Oswald et al., 2020), these fires may have occasionally escaped and affected forest composition on large areas on the SLSTJ’s plains, where fire-adapted species were frequent before European settlement.

Changes in Forest Composition

Our results show an increase of fast-growing and early-successional species in poplar, maple and jack pine since pre-settlement times, to the detriment of long-lived, late-successional species such as spruce, yellow birch, cedar, and other pines (red and white pines). This trend has been widely documented across temperate and boreal zones of eastern North America as a consequence of human activities (Friedman and Reich, 2005; Pinto et al., 2008; Hanberry et al., 2013; Thompson et al., 2013; Danneyrolles et al., 2016; Boucher et al., 2017; Terrail et al., 2019).

Although the effects of climate change cannot be discarded (Pederson et al., 2015; Boisvert-Marsh et al., 2019), anthropogenic fire disturbance most likely played an important role in these changes. During the settlement period, settlers cleared the forest and usually burnt the logging waste to facilitate the establishment of agriculture (Blanchet, 2003). This operation was often sloppy, and settlers did not take the necessary precautions to prevent fire from spreading to nearby forests. Aided by drought and high winds, these fires could have spread over large areas (Boucher et al., 2014), as it was the case in 1870 when a fire ignited by settlers, west of Lake St-Jean, escalated into a major conflagration that ran over 120 km on the plains (Blanchet, 2003). Despite the enactment of a fire exclusion policy enforced by law in 1870 (Blanchet, 2003), subsequent governmental reports indicated that more than 67% of the forest fires were anthropogenic between 1906 and 1941, with 6% caused by lightning, and 27% from unknown origin (Québec, 1906–1941). High fire activity during the settlement period has been documented in several regions of North America (Weir and Johnson, 1998; Lefort et al., 2003; Hessl et al., 2011; Stambaugh et al., 2018; Terrail et al., 2020) and Scandinavia (Lehtonen and Huttunen, 1997; Groven and Niklasson, 2005). While severe fire weather conditions may have facilitated fire spread (Lefort et al., 2003; Grenier et al., 2005), anthropogenic ignition appears to be the main factor explaining the high fire activity during the settlement period (Niklasson and Granström, 2000; Blanchet, 2003; Boucher et al., 2014; Terrail et al., 2020).

Consequently, pyrophile tree species benefited the most from the European settlement in the SLSTJ region. Trembling aspen (Supplementary Figure S5) increased throughout the study area and is now the most dominant tree. This species is known to quickly invade burned sites, both through wind-dispersed seeds and suckers (Perala, 1990; Bergeron and Charron, 1994; Bergeron, 2000; Fortin, 2008). As for jack pine, with its serotinous cones and early sexual maturity, it probably took advantage of successive fires to outcompete fire-sensitive species such as white spruce, balsam fir and larch (Larix laricina) (Burns and Honkala, 1990), especially on sandy deposits where it was already present and highly competitive (Greene et al., 2004; Mansuy et al., 2010). In fact, the increase of jack pine is a particularity of our study area because this taxon has generally decreased at the temperate-boreal ecotone since pre-settlement times (Radeloff et al., 1999; Friedman and Reich, 2005; Pinto et al., 2008; Hanberry et al., 2013). In these regions, fire suppression and logging have been identified as the main drivers of forest composition changes, whereas in SLSTJ the increase in jack pine is most likely related to high fire activity during the settlement period. Red maple is another taxon that may have benefited from fire disturbance. Red maple is by far the most abundant maple species in the study area today (Supplementary Figure S5) and is known to be favored by fire at its northern range limit (Zhang et al., 2015).

Logging has been identified as the main cause of vegetation change since pre-settlement in many regions at the temperate-boreal transition (Jackson et al., 2000; Friedman and Reich, 2005; Boucher et al., 2006; Dupuis et al., 2011; Hanberry et al., 2013) and was probably an additional causal factor of forest changes in SLSTJ. In a first phase, selective cutting that started with settlement quickly depleted red and white pines and then white spruce, which became the principal species harvested during the 19th century (Girard and Perron, 1989). During the second phase, extensive clear-cut logging during the 20th century targeted all conifers. In addition to directly decreasing conifer dominance, large openings in the canopy may have favored the invasion of fast-growing species such as poplars, white birch and red maple (Harvey and Bergeron, 1989; Fei and Steiner, 2009; Laquerre et al., 2009). Logging and fire may also have acted synergistically to favor disturbance-adapted species (Boucher et al., 2017).

The variable trend of balsam fir abundance in our study area concur with previous studies as its decline (Jackson et al., 2000; Danneyrolles et al., 2016), its stability (Pinto et al., 2008; Terrail et al., 2019), or its increase (Friedman and Reich, 2005; Hanberry et al., 2013) have all been reported. The increase of balsam fir’s relative frequency in the eastern sector of our study area can be explained by the maintenance of a large seedlings bank, which allows the species to rapidly re-establish following clear-cuts (Archambault et al., 2006; Prévost, 2008) or spruce budworm outbreaks (Baskerville, 1975; Blais, 1983). However, the sensitivity of balsam fir to fire probably led to its decline in the western part of our study area, where the fire-adapted jack pine, aspen and white birch (Burns and Honkala, 1990) greatly increased. Similarly, we observed increasing and decreasing trends of white birch abundance, paralleling other studies in boreal and temperate-boreal transition forests where the species was abundant prior to settlement [Friedman and Reich, 2005; Danneyrolles et al., 2019 (see Supplementary Figure S2), Marchais et al., 2020]. Although white birch in our study area may have benefited from fire, and from openings created by spruce budworm outbreaks and logging, it may have been outcompeted by aspen on burnt sites (Terrail et al., 2020), and likely affected by episodes of dieback during the 20th century (Pomerleau, 1991; Auclair et al., 1997).

Although logging and repeated fires have most likely played a role in the exclusion of the spruce-fir-white birch assemblages (clusters 1 and 2) from the plains, land cover change was probably a major contributing factor as urban and agricultural areas now occupy 32% of the plains’ surface (Figure 2C). Indeed, most of the plains’ 25-km² cells discarded from the analyses due to the lack of modern plots have had on average 62% of their surface transformed to agricultural or urban land uses and were dominated by clusters 1 and 2 before settlement (Supplementary Figure S2). Residual cells of these clusters now occur chiefly on the uplands or close to the Saguenay River, where human population density is low and agriculture not developed. Thus, several cells of clusters 1 and 2 on the plains have been either lost to land-use change (agriculture, urbanization) or replaced by new tree assemblages (clusters 3 and 4) better adapted to human disturbances (fire, logging, and land clearing).

Forest Management Implications

In light of our results, strict ecosystem-based forest management in SLSTJ should aim to increase spruce frequency and dominance, and decrease disturbance-adapted (fast-growing, shade intolerant) species such as aspen and red maple. Since spruce is still present throughout the territory, conversion from aspen to coniferous stands through natural succession or partial cutting (Brais et al., 2019; Prévost and Charette, 2019) are potential management solutions, although landscape fragmentation and continuing forest harvesting are likely to favor the retention of shade intolerant species such as white birch and aspen (Fortin, 2008). Thus, tree planting programs guided by past substrate preferences may be needed for spruce as well as for restoring the less abundant yellow birch, pines (red and white) and cedar to higher frequency and dominance levels.

However, in light of anticipated climate change, along with current uncertainties on future disturbance rates (Bergeron et al., 2006, 2010; Boulanger et al., 2019) and fire suppression effectiveness (Wotton et al., 2017), more flexible management options may be needed. Pre-settlement forest composition can provides targets to maintain or restore functional diversity and resilience (Boulanger et al., 2019). For example, species sensitivity or adaptations to disturbances can be considered in order to maintain functionally diverse and spatially heterogeneous stands and landscapes similar to pre-settlement conditions (Danneyrolles et al., 2020). Yet, it should be noted that our study only provides reference conditions for the period just before the European settlement, and that longer paleoecological records should also be considered to identify a wider range of vegetation targets for future climate trends and disturbance regimes.

Data Availability Statement

The datasets generated for this study are available on request to the corresponding author.

Author Contributions

SD, DA, and VD designed the study and methodology. SD extracted and compiled the data. SD analyzed the data with the help of VD. SD wrote the manuscript with substantial inputs from DA, YB, VD, and JL. All authors contributed to the article and approved the submitted version.

Funding

This study was financed by the Québec’s Ministère des Forêts, de la Faune et des Parcs (project no. 142332152-B).

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.

Acknowledgments

We thank two reviewers for their constructive comments, which greatly improved the first version of the manuscript. We also thank Rachel Bouserhal for English revision.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fevo.2020.00257/full#supplementary-material

References