Global change impacts on forest and fire dynamics using paleoecology and tree census data for eastern North America

Abstract

Key message

The tree census, paleopollen, fossil charcoal, human population, and climate data presented here provide unique support for important anthropogenic influences on fire over the last 2000 years in the eastern USA. This includes multiple instances of climate fire anomalies that may be best explained by the role of human-caused burning.

Context

The coupling of paleoecological and tree census data to address larger global change questions is a novel research approach to describe and ascribe recent vegetation dynamics vis-à-vis the climate versus disturbance debate.

Aims

The aims of the study are to (1) compile and compare pre-European settlement versus modern upland arboreal pollen and tree survey data from a large number of studies in various forest regions in the eastern USA, (2) analyze fossil charcoal dating back 2000 years for the northern versus central/southern tiers of the eastern USA, and (3) compare and contrast compositional and ecophysiological attributes for both datasets and temporal changes to known climate or disturbance phenomena to elucidate global change impacts and the drivers of forest change.

Methods

We analyzed paleoecological (pollen and charcoal) and tree census studies to compare protohistoric and modern vegetation assemblage for eastern North America, including the drivers of forest change. A total of seven forest types in the north and central regions of the eastern USA were used to co-analyze fossil pollen, fossil charcoal, and tree survey data.

Results

Disparities and consistencies existed when independently assessing witness tree and pollen records. Although forests north of the tension zone line (TZL) contained mostly Fagus, Pinus, Tsuga, and Acer witness trees, pollen records were dominated, as expected, by high-pollen-producing Pinus, Quercus, Tsuga, and Betula. Here, present-day pollen and tree survey data revealed significant declines in Fagus, Pinus, Tsuga, and Larix and increases in Acer, Populus, Fraxinus, Quercus, and Abies. South of the TZL, both witness tree and pollen records pointed to Quercus and Pinus domination, with declines in Quercus and Castanea and increases in Acer and Betula based on present-day data. Modern assemblages comprise tree genera that are increasingly cool-adapted, shade-tolerant, drought-intolerant pyrophobes. Paleocharcoal data from 1 to 1750 AD indicate a slight increase in burning in southern forests and stable levels in the north, despite the increasing cold associated with the Little Ice Age. The most significant increase in burning followed the dramatic increase in human population associated with European settlement prior to the early twentieth century.

Conclusion

Post-1940, fire suppression was an ecologically transformative event in all datasets. Our analysis identifies multiple instances in which fire and vegetation changes were likely driven by shifts in human population and land use beyond those expected from climate alone.

1 Introduction

During the last five centuries, vegetation in eastern North America has been impacted by a suite of global change phenomena (Abrams and Nowacki 2008, 2015; Munoz and Gajewski 2010; McEwan et al. 2011; Woodall et al. 2013). This includes the initiation and rapid expansion of Euro-American settlement, Native American depopulation, and dramatic changes in the magnitude, extent, and type of land use practices, abrupt shifts and/or reversals in fire regimes, outbreaks of insect and disease, and significant climate change (Crosby 1976; McAndrews 1988; Denevan 1992; Whitney 1994; Parshall and Foster 2002; Foster 2004; Munoz et al. 2010). At the beginning of this five-century period, the northern hemisphere was in the midst of the Little Ice Age (LIA; ca. 1400 to 1850; Fig. 1; Mann et al. 2009) and the initial stages of European settlement took place within this climatic regime. The LIA was followed by abrupt warming associated primarily with the end of a natural cooling period coupled with increased anthropogenic modifications to atmospheric chemistry (Ruddiman 2005; Mann et al. 2009; IPCC 2013). Human populations and their impacts on vegetation through land use have also changed appreciably during the late Holocene (Denevan 1992; Abrams and Nowacki 2008; Munoz and Gajewski 2010; Nowacki et al. 2012; Munoz et al. 2014). Nevertheless, we still have only a marginal understanding of the role of climate and disturbance and their interactions with vegetation dynamics, past and present, for most regions (Rhemtulla et al. 2009). The importance of climate versus human impacts on protohistoric ecosystems in eastern North America is a highly debated issue (Munoz et al. 2010; Pinter et al. 2011; Marlon et al. 2013). One argument emphasizes the role of climate driving fire and vegetation dynamics (Parshall and Foster 2002; Shuman et al. 2004; Pederson et al. 2015); another argues that human-caused disturbance, including intentional burning, has been the primary driver, particularly during the second half of the Holocene (Guyette et al. 2006; Steyaert and Knox 2008; Nowacki and Abrams 2008, 2015). A more complete understanding of past human-fire-climate-vegetation dynamics and their anomalies requires additional research (Munoz et al. 2014; Abrams and Nowacki 2015).

Fig. 1

Climatic data (from year 500 to 2000, left to right) for the eastern USA including temperature anomaly for the northeastern USA and Palmer drought severity index (PDSI) from three sites north (top panel) and three sites south (middle panel) of the tension line. The PDSI data from the three northern locations were in Vermont, Michigan, and Wisconsin, while the three southern sites were in Virginia, south central Pennsylvania, and Kentucky. The annual PDSI values were smoothed using a 20-year running average. Climate data are from Mann et al. (2009) and the North American Drought Atlas (http://iridl.ldeo.columbia.edu/SOURCES?.LDEO/.TRL/.NADA2004/.pdsi-atlas.html)

The study of fossil charcoal as an indicator of fire has helped elucidate disturbance regimes and their impacts in pollen interpretation (Patterson and Sassaman 1988; Parshall and Foster 2002). Charcoal data from sediment records can be used to provide information about past fire activity from local to global scales (Marlon et al. 2008). Nevertheless, the role of fire, including its origin, drivers, and extent in various forest types and regions, remains a contentious idea among historical ecologists (Pinter et al. 2011; Marlon et al. 2013; Abrams and Nowacki 2008, 2015). Opinions differ about the roles and relative strength of climate versus human (anthropogenic) drivers of fire, including the Early Anthropocene Burning Hypothesis. It deals with the possibility that the early human use of fire was profound, frequent, and prevalent, resulting in it becoming an early ecological driver (Ruddiman 2005; Marlon et al. 2013; Abrams and Nowacki 2015). This presumes that small populations of humans were able to burn expansive areas. Other issues that need to be resolved are the extent of lightning as an ignition source and whether burning was localized or ubiquitous in various vegetation types in the eastern USA (Ruffner and Abrams 1998; Nowacki et al. 2012).

The impact of present and future climate change and other global change phenomena on world ecosystems is one of the premier research topics for ecologists and environmental scientists of this era (IPCC 2013). Research techniques in the field of historical ecology have particular relevance for assessing the impacts of climate and humans on ecosystems because what has happened in the past may provide important clues of what will happen in the future (Foster 1998; Egan and Howell 2001). For example, it is possible to study how various ecosystems responded to abrupt climate change (cooling and warming) and shifts in disturbance regimes in the past (Jackson 2006; Booth et al. 2012; Nowacki and Abrams 2015). Studies of this type have often involved the use of early land survey (witness trees and associated line notes) or paleoecology data (sediment pollen and charcoal in lakes, bogs, and caves). Comparative tree censuses using witness and modern tree data can provide fairly reliable information about forest change spanning several centuries, whereas paleoecology chronologies often span many millennia (Whitney 1994; Schulte and Mladenoff 2001; Foster 2004; Munoz et al. 2010).

Paleoecology and tree census studies have vastly contributed to our understanding of vegetation dynamics and the impacts of climate and disturbance as vegetation drivers. These include broad-scale and long-term climate change over centuries or millennia, as well as changes and interactions from both anthropogenic and natural changes in climate and perturbations. The reaction of ecosystems to global change factors can be assessed by examining the impacts of changing environmental conditions in the past, including the rise and fall of various vegetation. Indeed, recent paleoecological modeling studies have simulated responses to future environmental change (Jackson et al. 2009; Hannah et al. 2014). This type of information is critical to ecologists, land managers, and policy makers trying to make informed decisions about conservation, environmental change, and ecosystem resilience.

Many comparative tree studies have been amassed to assess post-European settlement changes to investigate the impacts of recent climate change on vegetation dynamics (see Appendix S1 of Nowacki and Abrams 2015). The same can be said for paleoecology studies (Davis 1963; Webb III et al. 1981; Schwartz 1989; Paciorek and McLachlan 2009). The coupling paleoecological and tree census data to address larger global change questions seems to us an important, but underutilized, research approach. The potential exists for studies to more accurately describe and ascribe recent vegetation dynamics to its causal factors, particularly in relation to the climate versus disturbance debate. A recent paper by Nowacki and Abrams (2015) presented a unique approach to untangle the role of climate versus disturbance as forest change drivers over multiple centuries. This was accomplished by categorizing major tree species/genera into temperature, shade tolerance (succession), and pyrogenicity classes based on ecophysiological characteristics, and then applying those classes to comparative studies of past and present tree censuses. In contrast, paleoecological studies have typically concentrated on documenting vegetation dynamics and their relationship to environmental changes over the Late Quaternary (Davis 1969; Webb 1983; Russell et al. 1993; Shuman et al. 2004; Williams et al. 2004). Relatively few paleoecological studies investigated the impacts of late Holocene climatic variability on vegetation dynamics.

For understanding the relative roles of climatic variability and disturbance regimes on vegetation communities, we use an approach that integrates historical and modern tree surveys with fossil pollen and charcoal data (cf. Dawson et al. 2016). Witness trees (trees used to demarcate property corners in early historical land surveys; also referred to as corner trees and bearing trees) provide a “snapshot” of forest composition during the time of land survey, spanning from the 1600s (eastern colonies) to the early 1900s (Minnesota) (Whitney 1994; Schulte and Mladenoff 2001). Although there are some shortcomings (Schulte and Mladenoff 2001; Whitney and DeCant 2001), witness tree data probably provide the best depictions of arboreal vegetation during the onset of European settlement, especially when compiled at the landscape scale (Delcourt and Delcourt 1996). Witness tree data are extensively used in historical ecology studies throughout the eastern USA (Whitney 1994). When paired with current vegetation surveys (e.g., Forest Inventory and Analysis data), changes in arboreal composition and structure can be deciphered and associated ecological drivers assigned (Nowacki and Abrams 2008, 2015).

Paleoecological data (biological and geochemical remains preserved in sedimentary archives) record environmental changes over centuries to millennia through the use of proxies. In this study, we use fossil pollen as a proxy of past vegetation composition and change (Davis 1969; Russell et al. 1993; Shuman et al. 2004; Williams et al. 2004) and fossil charcoal as a proxy of past biomass burning (Power et al. 2010; Marlon et al. 2013). The combination of these different data types should give a broader perspective of how past and projected climatic change will impact vegetation communities. More specifically, the objectives of this study are to (1) compile and compare pre-European settlement versus modern upland arboreal pollen and tree survey data from large number of studies in various forest regions in the eastern USA, (2) analyze fossil charcoal dating back 2000 years for the northern versus central/southern tiers of the eastern USA, and (3) compare and contrast compositional and ecophysiological attributes for both datasets and temporal changes to known climate or disturbance phenomena to elucidate global change impacts and the drivers of forest change.

2 Methods

A total of seven forest types in the north and central regions of the eastern USA were used to co-analyze fossil pollen, fossil charcoal, and tree survey data (Tables 1 and 2; Fig. 2). These forest types were assembled, and the tension zone line (TZL) struck using ecological subsections as base units (Cleland et al. 2007). Four of the forest types are located north of the TZL, and three forests types are located to its south (Fig. 2). The TZL is a boundary between two distinct floristic zones, in this case conifer-northern hardwood and sub-boreal to the north and oak-pine to the south (Curtis 1959). Most of the pollen chronologies used for study were extracted from Neotoma V1.0, a relational multiproxy paleoenvironmental database (Grimm 2008), with additional data extracted from published literature (Table 1). Using a criteria similar to Munoz et al. (2010), we excluded a pollen record if it (a) did not contain dated radiocarbon years within the study periods (pre-European and Little Ice Age versus modern), (b) contained a hiatus of > 500 years of accurately dated samples within the study period, (c) had no chronological control (other than the top of the core) within the study period, (d) had samples from large bodies of water (> 115 ha) or very small ponds (< 1 ha) (leaving only larger ponds and small lakes for our analysis), and (e) had samples located in present-day urban areas that might have impacted modern vegetation assemblages. By applying these criteria, we filtered our initial population of 129 pollen chronologies to 76 for this study.

Table 1 List of lake basins within each forest type included in the study with the two periods of comparison (pre-European settlement versus present day), the % of tree pollen versus total pollen (all upland vegetation) in the sample and citation. NA = not available in the Neotoma database. GL_BM = Great Lakes beech-maple; GL_CNH = Great Lakes conifer-northern hardwood; NE_CNH = Northeast conifer-northern hardwood; SC = sub-boreal conifer; Cent_OP = Central oak-pine; GL_OP = Great Lakes oak-pine; NE_OP = Northeast oak-pine

Table 2 Sites used for paleocharcoal analysis by region. The northern and southern zones were divided at 43° N latitude across the eastern USA (see map in Fig. 2)

Fig. 2

Pre-European settlement forest biomes of North America assembled from ecological subsections of USA, pollen study site, and tension zone locations (Cleland et al. 2007; Nowacki and Abrams 2015). Black dots are the location of the paleopollen study sites (Table 1), and red dots are the location of paleocharcoal study sites (Table 2)

We used the Neotoma age-depth information for our chronologies, which typically use the European settlement horizon as a chronological control. Whenever possible, we used age-depth models based on calibrated 14C dates, but for this time period, there is not a major difference between radiocarbon and calibrated ages. Arboreal pollen percentages for selected taxa (Table 3) were calculated using the terrestrial pollen sum (Table 1). First, we eliminated all herb and shrub pollen; then, the remaining tree pollen were relativized for each species with the total tree pollen count equal to 100%. Taxonomic resolution of pollen data were reduced to the genus level in arboreal taxa where species-level identification were made (e.g., Acer saccharum and Acer rubrum were aggregated to Acer) because most records are not identified beyond the genus level. The protohistoric period encompasses Native American occupation during the Little Ice Age (ca. 1400 to 1850 AD), while the contemporary period reflects post-Euro-American land use and climate change after ca. 1850.

Table 3 Common eastern North American tree genera classified by temperature, shade tolerance, and pyrogenicity (adapted from Nowacki and Abrams (2015)). Average annual range temperatures (in parentheses) derived from the Tree Atlas database (Prasad et al. 2007–ongoing) were used to define temperature classes for species (cold ≤ 6.6 °C; cool = 6.7–10.7 °C; warm = 10.8–14.0 °C). Range temperatures for genera were averaged from component species. Pollen representation in sediment samples is based on Prentice and Webb III (1986). Pyrophilic genera are also classified as drought tolerant (including moderately drought tolerant, indicated with DT

Charcoal data are from the Global Charcoal Database (GCD, v. 2) (Daniau et al. 2012), including 13 records from north of 42° N latitude (2430 samples) and 11 records south of 42° N latitude (428 samples), and were used to reconstruct trends in biomass burning for the eastern USA during the past 2000 years (Table 2). Data from the GCD are typically standardized to produce anomalies from a shared base period for multiple records. Standardizing the charcoal values makes comparisons possible across many records, despite differences in specific particle size classes measured, lab, and analytical methodologies. For standardization, each record was “presampled” to a common temporal window (10 years in this case). A weighted mean was calculated if multiple values existed within a given 10-year period for each record. The presampled records were subsequently rescaled, transformed, and converted to z scores following the protocol detailed in Power et al. (2010) and implemented via the PfTransform (paleofire) function in the paleofire R package (Blarquez et al. 2014). A base period of 2000–200 cal year BP was used to calculate the z scores. The standardized and normalized charcoal records were then smoothed using a 200-year window width with the pfCompositeLF (paleofire) function to create a composite for the northern and southern zones, divided at 42° N latitude, and for the entire eastern USA, along with 5% and 95% bootstrap confidence intervals.

Climatic data was compiled and synthesized for the eastern USA, including temperature anomaly for the northeastern USA and Palmer drought severity index (PDSI) from three sites north and three sites south of the tension line. The east–west gradient of PDSI data are from three northern locations that were in Vermont, Michigan, and Wisconsin while the three southern sites were in Virginia, south central Pennsylvania, and Kentucky. The annual PDSI values were smoothed using a 20-year running average. The mean annual temperature and PDSI data were obtained from Mann et al. (2009) and the North American Drought Atlas (http://iridl.ldeo.columbia.edu/SOURCES?.LDEO/.TRL/.NADA2004/.pdsi-atlas.html), respectively. Temperature data were obtained by extracting data from the paleoclimate dataset reconstructed by Mann et al. (2009) covering grid cells in the eastern USA.

For tree census data, we identified 52 studies in the eastern USA that produced a total of 192 datasets available for analysis (Nowacki and Abrams 2015). Early tree surveys chronicle the westward progression of European land acquisition, with some dating back to the seventeenth century along the Atlantic Coast (Whitney and Davis 1986; Foster et al. 1998; Thompson et al. 2013). These are divided into two time periods, the early historic (1620 to 1850 AD) and contemporary (1940 to 2000 AD), making these as similar as possible with the paleopollen data.

Population estimates for the eastern USA were obtained from the HYDE 3.1 global population dataset (Klein Goldewijk et al. 2011) by aerially averaging data from grid cells between 25° N and 49° N latitude and 95° E and 65° E longitude. Temperature data were obtained by extracting data from the paleoclimate dataset reconstructed by Mann et al. (2009) covering grid cells in the eastern USA.

2.1 Data analysis

Major tree genera in eastern North American (16 genera distributed across the seven forest types in our study region) were classified by temperature, shade tolerance, and pyrogenicity (including drought tolerance) based on available data, published literature, and authors’ knowledge (Table 3). Temperature classes were established using actual temperature data from the Climate Change Tree Atlas (Prasad et al. 2007–ongoing). Each tree genus was classified into one of four temperature classes (cold, cool, warm, and hot) based on the average annual temperature within its ecological range (US distribution; Nowacki and Abrams 2015). Genera were also categorized by shade tolerance (intolerant, intermediate, and tolerant) and pyrogenicity (pyrophilic and pyrophobic) based on their known life history and physiological characteristics (Table 3). Individual tree percentages from pollen and survey data were tallied by temperature (cold, cool, warm), shade tolerance (intolerant, intermediate, tolerant), and pyrogenicity (pyrophilic, pyrophobic) classes based on the lowest taxonomic level reported for each dataset. All of the genera classified as pyrophilic were also classified as drought tolerant (Burns and Honkala 1990); thus, changes in pyrogenicity are indicative of changes in climate conditions and drought tolerance. Absolute percentage changes were then calculated for each class by subtracting protohistoric from contemporary abundances, and two-tailed t tests were performed between the two time periods to test for statistical differences in both the pollen and tree survey data.

3 Results

PDSI was consistently higher (less droughty) in northern forests than southern, where it was typically above and below zero, respectively (Fig. 1a, b). However, PDSI generally increased after about 1900 in southern regions, but within the historic range of variability. Temperatures across the eastern USA declined from 500 to ~ 1700 AD and abruptly increased thereafter (Fig. 1c).

Pollen assemblages from the Great Lakes beech-maple (Fagus-Acer) forest are dominated by oak (Quercus) and pine (Pinus) in the protohistoric and contemporary periods because these genera are high-pollen producers (Tables 3 and 4). Inconsistent with the principal components by which this forest type was named (sensu Braun 1950), beech and maple, with intermediate and low pollen production, were vastly underrepresented in the pollen record. However, protohistoric tree surveys supported the forest type designation, with elevated occurrences of beech and maple followed by oak, although this is based on only two sites. Arboreal pollen assemblages exhibited little change over the study period, but tree survey data document a large decline in beech (− 29%). Pollen assemblages in the protohistoric Great Lakes conifer-northern hardwood biome were dominated by high-pollen producers of pine and birch (Betula) and intermediate pollen producer hemlock (Tsuga), whereas the tree survey data indicate that these forests were dominated by pine, hemlock, and maple (Table 4). In this region, tree survey data document significant changes in 11 species with the largest increases for maple, aspen (Populus), and oak and largest decreases for pine and hemlock. Although no statistically significant changes in the relative abundances were found in any of the pollen taxa, the direction of change for the important genera (maple, pine, and hemlock) correspond to that found in the tree survey data.

Table 4 Pre-European settlement (protohistoric) and present-day percentages and temporal change of pollen and tree surveys for the major tree genera in forest types north of the tension zone line. Number of sites for pollen and survey trees, respectively, is in parentheses. Forest type abbreviations are found in Table 1

Protohistoric pollen assemblages in the Northeast conifer-northern hardwoods exhibited high abundances of birch, pine, hemlock, oak, and beech (Table 4), whereas tree survey data indicated that these forests were disproportionally beech dominated, followed hemlock and maple. Over the period of study, the pollen data document statistically significant declines in the abundances of hemlock and beech and increases in birch and maple, all of which correspond to trends observed in the tree survey record. Protohistoric pollen assemblages in sub-boreal conifer forests were again dominated by high-pollen producers pine and birch, while tree survey data document a diverse forest community dominated by pine, larch (Larix), spruce (Picea), and birch (Table 4). Over the period of record, sub-boreal forest pollen assemblages document statistically significant increases in oak, aspen, and ash (Fraxinus), whereas tree survey data statistically document a more robust array of compositional shifts, with significant increases in aspen, ash, maple, and fir (Abies) and decreases in pine, larch, and birch.

Pollen assemblages in the oak-pine forests south of the TZL (including the Central, Great Lakes, and Northeast regions) are dominated by oak and pine over the period of study, with few significant shifts in the abundances of any genera documented (Table 5). Correspondingly, tree survey data recorded the same dominants over the period of record. However, tree survey data did record significant declines in oak and chestnut and significant increases in maple. In most cases, the direction of change recorded in tree surveys paralleled those in the pollen record.

Table 5 Pre-European settlement (protohistoric) and present-day percentages and temporal change of pollen and tree surveys for the major tree genera in forest types south of the tension zone line. Number of sites for pollen and survey trees, respectively, is in parentheses. Forest type abbreviations are found in Table 1

Compositional differences between pollen and tree survey data were generally consistent with differences in pollen productivity (Table 3), with high-pollen-producers pine, birch, and oak overrepresented and low-pollen-producers maple and aspen underrepresented relative to tree survey data. Overall, forests north of the TZL were dominated by genera typical of conifer-northern hardwoods, including pine, birch, beech, hemlock, and maple. Forests south of the TZL were dominated by oak followed by pine in both data types. Here, high-pollen-producing pine and birch were overrepresented and low-pollen-producing maple and aspen were underrepresented in the pollen record relative to tree surveys.

Tree genera were partitioned into temperature, shade tolerance, pyrogenicity (including drought classes) and applied to pre-European fossil pollen and witness tree data (Fig. 3). The Great Lakes beech-maple pollen data were dominated by warm-adapted trees that were evenly split between two shade-tolerant classes (intolerant and intermediate) and both fire/drought classes. Witness trees from the same forest type had few cold-adapted, intolerant pyrophiles. The Great Lakes conifer-northern hardwood pollen data were dominated by cold-adapted, shade-intolerant (early successional) trees that were nearly evenly split between pyrophiles and pyrophobes. Witness trees differed from this by having had a higher proportion of cool-adapted, shade-tolerant, drought-intolerant/pyrophobic trees. In protohistoric Northeast conifer-northern hardwoods, pollen data was dominated by cold-adapted, drought-intolerant pyrophobes that were evenly split between shade intolerant and tolerant trees. Witness trees had a higher proportion of cool- and warm-adapted, shade-tolerant pyrophobes. The sub-boreal forest was greatly dominated by cold-adapted, shade-intolerant trees with more drought-tolerant pyrophiles in the pollen data and more drought-intolerant pyrophobes in the witness tree data.

Fig. 3

Tree genera percentage in each of the ecophysiological classes for pollen versus witness tree data from forest biomes north of the tension zone line. Pyrophilic genera are also rated as drought tolerant; pyrophobic are drought intolerant (low to moderate)

South of the TZL, pre-European Central oak-pine forests were composed mostly of warm-adapted, shade-intermediate, drought-tolerant pyrophiles, with a slightly higher proportion of cold and shade-intolerant trees in the pollen versus witness tree data (Fig. 4). The Great Lakes oak-pine forest too was dominated by warm-adapted, shade-intermediate, drought-tolerant pyrophiles in both datasets. The Northeast oak-pine forest had a fairly even distribution of tree species across all ecophysiological classes, but tended to be dominated by warm-adapted, shade-intermediate pyrophiles, more so in the witness tree record (Fig. 4). Pollen data had somewhat more cold-adapted, intolerant and drought-intolerant pyrophobic trees. Overall, forests south of the TZL contained a higher proportion of warm-adapted, shade-intermediate, drought-tolerant and pyrophilic trees prior to European settlement, whereas northern forests were dominated by cold-adapted, shade-intolerant species that were fairly evenly split between pyrophobes and pyrophiles. However, in the north, witness tree data consistently had higher values for shade-tolerant, drought-intolerant, and pyrophobic genera, which can be attributed to differential pollen production among genera.

Fig. 4

Tree genera percentage in each of the ecophysiological classes for pollen versus witness tree data from forest biomes south of the tension zone line. Pyrophilic genera are also rated as drought tolerant; pyrophobic are drought intolerant (low to moderate)

Comparing mean pre- and post-European settlement changes in the transformed ecophysiology data, the Great Lakes beech-maple forests exhibited a decrease in warm-adapted, shade-tolerant, drought-intolerant pyrophobes in both data types (Table 6). Temporal changes in the Great Lakes conifer-northern hardwood forest included an increase in warm-adapted, shade- and drought-tolerant pyrophilic survey trees, whereas pollen data showed a decline in drought-tolerant pyrophilic trees and neutral shade tolerance changes. Northeast conifer-northern hardwoods had a decrease in warm-adapted, shade-tolerant trees and an increase in drought-tolerant pyrophiles in both datasets. The sub-boreal conifers exhibited an increase in cool and warm-adapted and shade-tolerant trees and a decline in drought-tolerant pyrophiles, the latter only in pollen data. South of the TZL, Central oak-pine forest type had a decline in warm, intolerant, drought-tolerant/pyrophilic trees in both datasets (Table 6). In the Great Lakes oak-pine forest type, tree composition exhibited an overall decline in temperature and drought/fire adaptation and an increase in shade tolerance in the tree survey data. Ecophysiological changes in Northeast oak-pine forests varied between pollen and witness tree data, while both data types showed a general cooling of temperature adaptation they differed in the direction or lack of changes in shade and drought tolerance and pyrogenicity. Across all three ecophysiological classes, the magnitude of temporal change (as an absolute value) was significantly greater (P < 0.001) for tree surveys than for pollen data. The overall temporal change was also somewhat higher in southern versus northern forests (Table 6). In most cases, the direction of change was similar between pollen and witness tree data. Northern forests were highly variable in terms of temperature, tolerance, and pyrogenicity/drought changes, going up in two and down in two forest types. In contrast, southern forests exhibited an overall decline in temperature adaptation and drought/ pyrogenicity but were inconsistent in shade-tolerant changes. In five of seven forest types, aggregated tree genera registered temperature declines from past to present and moved in the opposite direction of the prevailing post-1820 warming (Table 6; Fig. 1).

Table 6 The mean percent change (as an absolute value) and the overall direction of change in the three ecophysiological classes (pyrogenicity and drought tolerance combined) from pre-European settlement to modern day in pollen and tree survey data for each forest type north (n = 4) and south (n = 3) of the tension zone line, including a grand mean for each biome. Forest type cells without an arrow had neutral changes for that parameter. See Table 1 for forest type abbreviations

Trends in standardized charcoal accumulation rates (the charcoal index) show different trajectories in biomass burning in the northern and southern zones of the eastern USA (Fig. 5). In the north, variations in standardized charcoal accumulation rates show relatively minor variations from 1 to 1750 AD, with the maximum occurring between 1 and 200 AD and the minimum occurring around 1200 AD. After European settlement, around 1750 AD, charcoal levels increased rapidly towards present. Confidence intervals in the most recent decades expand substantially from their previous range, indicating increased variability in the time series. In the southern zone, the charcoal records are much lower resolution (i.e., each sediment core was sampled less frequently than what is typical in the north), and as a result, there are less than one fifth of the samples than what is available in the north. Nevertheless, the combined data from the southern zone show a gradual upward trend from 1 to 1750 AD, with the minimum around 1 AD and the pre-European maximum around 700 AD. After 1750 AD, charcoal levels in the southern zone form a parabola, increasing rapidly, cresting around 1900, then declining thereafter. Taken as a whole (24 records; 2858 samples), the eastern US standardized charcoal index shows a flat trend, until the large increase in burning during the European settlement era. A decline in recent centuries is also evident in the south and the eastern USA as a whole. Population from the HYDE 3.1 database shows a gradual increase from 1 to 1500 AD. Population declines from 1500 to 1600 AD (consistent with the Native American pandemic), increases slightly from 1600 to 1700 AD, and shows a rapid increase subsequently towards present after European colonization. The initial increase in charcoal predated abrupt warming that started about 1820 (Figs. 1 and 5).

Fig. 5

Eastern North America biomass burning (charcoal index) trends for regions a north of tension zone line, b south of tension zone line, and c entire eastern USA, all smoothed with a 250-year (black) line and shown with 95% bootstrap confidence interval (gray lines) and a linear regression line for data between year 0 and 1700. Panel d is the estimated human population for the eastern USA from year 0 to 2000 (from the HYDE 3.1 global population dataset (Klein Goldewijk et al. 2011))

4 Discussion

The disparities of using pollen-based versus tree census data, and their inherent differences, in vegetation reconstruction have been previously recognized by paleoecologists (Dawson et al. 2016). Pollen data are parsed from accumulations of particulate matter cast from large catchments that collect incrementally over long-time periods. When considering species differences in pollen production, dissemination and preservation, differential rates of deposition and mixing of sediments over time, radiocarbon-dating uncertainties, and the segment sizes required to gain adequate pollen counts, it is not surprising that the resulting data are rather low spatiotemporal resolution (within-core samples often spanning multiple years representing hundreds of hectares of depositional input). Exceptions occur where varves exist in sediment samples. In contrast, survey or census data are obtained from the direct measure of living trees surrounding a specific point at a specific time, thus resulting in high spatiotemporal resolution. However, this method has its limitations as well. For instance, witness trees have their own set of species selection biases contingent on tree size, wood durability and commercial value, surveyor instructions and partialities, and ease of reach and scribing survey trees (Bourdo 1956; Whitney 1994; Black and Abrams 2001; Schulte and Mladenoff 2001). Despite their respective weaknesses, both data types are invaluable to the understanding of historical ecology, particularly when used in the proper context. This compilation adds to this discussion by quantifying the disparities among these two data types in forest reconstruction and tracking compositional change over time, with cautions and remedies.

Our reconstruction of forest composition from fossil pollen versus tree surveys revealed important similarities and differences between the two data types. The disparities can be largely explained by genera variations in pollen production and dispersal (Table 3). Almost universally, high-pollen producers (birch, pine, and oak) were over represented, intermediate pollen producers (beech, hemlock, and elm) were adequately represented, and low-pollen producers (fir, maple, and aspen) were underrepresented in pollen records relative to tree surveys. Discrepancies were most stark and numerous north of the TZL where low-pollen-producing genera (fir, maple, cedar, ash, larch, and aspen) were more abundant, and thus largely underrepresented in pollen records relative to other genera (Table 4). Indeed, the biggest disparity between the two data types existed in the Great Lakes beech-maple protohistoric forest, which was dominated by oak and pine as pollen versus beech and maple as witness trees (the later consistent with the forest type’s namesake). Discrepancies were less apparent in oak-pine forests south of the TZL, having reasonably similar pollen and tree survey composition and dominance except for the high-pollen-producing birch. Paleoecologists have long recognized the incongruence between tree pollen counts and tree census data (Fagerlind 1952; Potzger et al. 1956; Davis and Goodlett 1960), and have made corresponding adjustments in their interpretation of forest composition and dominance, including pollen correction factors (Davis 1963; Andersen 1970; Prentice and Webb III 1986; Jackson and Williams 2004). We used correction factors to categorize pollen-to-tree census relationships (Table 7). Although correction factors seem like an attractive method to recalibrate pollen to better equate/correspond to tree census data, problems exist with that approach as (1) the relationship between vegetation and pollen percentage is not linear, (2) a linear correction factor assumes an intercept of zero (e.g., zero pollen percent = tree absence), and (3) the relationship between vegetation and pollen percentage is not constant across space or time (pers. comm, Dr. Samuel Munoz, November 26, 2014). However, correction factors are valuable in that they capture the magnitude of difference (thus the scope of the problem) between pollen records and tree census data at the species/genera level. Palynologists are developing more robust and elaborate methods to close the gap between pollen and actual tree census data (Paciorek and McLachlan 2009; Sugita et al. 2006, 2010).

Table 7 Correction factors by tree genera by individual study and summarized (averaged) by continent and for all studies. Pollen representation relative to tree census data were based on the following correction factor ranges: underrepresented > 1.25 (blue), representative = 0.75–1.25, and overrepresented < 0.75 (orange). Final genera designations of underrepresented, representative, or overrepresented were based on correction factors from all studies (last column)

When evaluating temporal variation in forest composition and dominance, we found that the finer resolution of tree census data captured many important changes in genera abundance not seen in the pollen data. This included the well-documented declines in oak and chestnut and increases in maple census trees south of the TZL (Nowacki and Abrams 2008). North of the TZL, census trees routinely showed significant declines in beech and pine and increases in maple and aspen, changes that were only marginally evident in pollen data. While these compositional shifts make sense ecologically (see below), it is important to point out that witness tree surveys are very different in methodology and point location from current tree surveys. The greatest methodological departures occur in the far eastern portions of the USA (original 13 colonies) where metes-and-bounds surveys were conducted, and property corners were oftentimes marked by widely and irregularly spaced, single trees (Thomas-Van Gundy and Nowacki 2013). The General Land Office Survey (GLO), conducted from Ohio westward, generally recorded two to four witness trees at section and quarter-section corners and at river crossings (“meander posts”) (Bourdo 1956; Whitney 1994; Black and Abrams 2001). The point-based tree data recorded by GLO surveys differ substantially from modern-day surveys (e.g., Forest Inventory and Analysis) which often use replicated fixed-area plots within forest stands. In contrast, tree variation in fossil pollen is recorded using a uniform analytical technique often from a single sediment core from a lake or bog, which minimizes spatial-temporal errors. Witness trees often only make up a small portion of any one forest or location, whereas paleoecological data are derived from a larger geographic catchment representing both local and regional vegetation (Jackson 1994; Fuller et al. 1998). However, individual witness tree warrants were often combined into connect drafts to increase number of survey trees and geographic area relative to a single tract of land. Nevertheless, when interpreting these types of data, one should be aware of the inherent differences among pollen, witness tree, and modern forest surveys.

Many palynological studies have reported substantial changes in forest composition from the pre-European settlement to present day. After European settlement, forests in the northeast had considerably less beech, chestnut, hickory, and hemlock and more birch, spruce, and fir (Russell and Davis 2001; Munoz et al. 2010). In this study, we also found that hemlock and beech pollen declined and that of birch increased in the Northeast conifer-northern hardwood forest. Former chestnut-hemlock forests in north-central Massachusetts are now dominated by oak as a result of human activity, including the introduction of chestnut blight (Fuller et al. 1998; Foster et al. 2002b). Following European settlement in New England, conifer-hardwood forests had declines in hemlock and beech and increases in white pine or pitch pine that invaded the open and abandoned forests (Parshall and Foster 2002). In our study, pine pollen dominated most forest types, but did not change significantly, likely reflecting its status as a high-pollen-producing species (Davis and Goodlett 1960). In southeastern New York, forests exhibited a decline in pine (mostly Pinus strobus) and an increase in birch and oak (Pederson et al. 2005). A substantial increase in yellow birch (Betula alleghaniensis) pollen during the late Holocene in the western Great Lakes region is attributed to a wetter climate (increased lake levels; Jackson and Booth 2002). In this study, birch pollen did not increase significantly following European settlement in the Great Lakes forests. A broad-scale study by Gajewski (1987) reported that hemlock and beech declined in northern New England while expanding in the Midwest over the past 500 years. This is consistent with the significant pollen declines in our Northeast conifer northern hardwoods, but we found no such increase of these genera in the Great Lakes beech-maple and conifer-northern hardwoods (spanning the Midwest). Davis et al. (1998) reported that former pine forests in northern Wisconsin and Michigan became dominated by hemlock and northern hardwoods in more recent years. Our Great Lakes-conifer-northern hardwood forests were dominated by these species (including birch), both pre- and post-European settlement, with pine and hemlock exhibiting non-significant declines.

Pollen records indicate that strong climate-vegetation relations existed prior to European arrival in certain segments of North America (Watts 1979; Shuman et al. 2004; Gill et al. 2009; Munoz et al. 2010). However, climatic controls of vegetation seemingly dissipated during the onset of European settlement, with humans (disturbance) becoming the primary driver of compositional change and leading to regional vegetation homogenization (Fuller et al. 1998; Foster et al. 1998). This post-European disruption of climate-vegetation relationships was recently corroborated in comparative tree survey data, whereby compositional changes have been attributed to changes in land use history rather than climate (Nowacki and Abrams 2015). The synchrony of European population expansion and Native American depopulation differed on the North American continent. Native American populations started declining from the very instant of European arrival driven largely by introduced diseases (Crosby 1976; Dobyns 1993). Starting along travel route of DeSoto’s 1539–1542 expedition and from early settlements along the East Coast, an ensuing pandemic was unleashed, quickly sweeping across America along major travel corridors and decimating Native American populations throughout the sixteenth and seventeenth centuries (Richter 2001; Mann 2005). Coupled with intertribal conflicts and social upheavals, population reductions of 80–90% likely happened among Native Americans. In many locations, this “wave of death” preceded European settlement by a hundred years or more, especially west of the Appalachians.

Extensive settlement of the eastern USA was not complete until about 1850, effectively expanding westward from the East Coast (Hart and Buchanan 2012). Concurrent with European settlement, the exploitation of forests occurred at an ever-increasing pace until the early twentieth century (Whitney 1994; Abrams 2003; Nowacki and Abrams 2008; Rhemtulla et al. 2009). Most of the eastern seaboard and the Ohio Valley were already logged at least once by the mid-nineteenth century. By 1920, approximately 99% of the original forest was gone. Not only were the original forests cut, but there was a large loss of forest area to land clearing for agriculture during the nation-building period. The “Great Cutover” logged billions of board feet of timber in the eastern and western USA and produced vast areas covered in “slash” (logging debris; Whitney 1994; Abrams 2003; Nowacki and Abrams 2008). As the slash dried, huge wildfires followed, which burned with an intensity not experienced in the original forest (Pyne 1982; Fig. 2). These fires ushered in the fire suppression (Smokey Bear) era in the USA starting in the 1930s. The most noticeable feature of the 2000-year fire history here is a large increase in burning during the past ~ 300 years associated with European land clearance for agriculture. Biomass burning peaked between 1800 and 1900 AD based on historical records, followed by substantial declines after ca. 1940.

The charcoal index between 1 and 1750 AD in this study suggests stable levels of burning in the north and a slight increase in the south. It is interesting to note that southern forests exhibited a slight increase in burning associated with the Medieval Warm Period (ca. 950 to 1250 AD) but did not decline with the LIA (ca. 1350–1850; Mann et al. 2009). In northern forests, the charcoal index dipped between 1000 and 1200 AD and then increased at the start of the LIA, the opposite of what would be expected from climate control. However, the greater abundance of samples in the north makes north/south comparisons difficult. Munoz et al. (2010) reported a steep rise in burning after 1000 AD that continued through the LIA in the northeastern USA. A synthesis of pre-European fire return intervals for the eastern USA report a north to south trend, with generally longer fire frequency in the cool north (> 16 years) and shorter in the hot south (< 2 years; Guyette et al. 2012). A slight increase in fires in the southern zone prior to 1750 AD corresponds with a gradual increase in Native American populations estimated during that time, so the increased burning may have been due to a growing number of human-set fires.

Based on fire-scar data for the eastern USA, the major factors controlling late Holocene differences in fire regimes are human population density, culture, and annual drought (Guyette et al. 2006, 2012). For the period between 1650 and 1930, these authors report stable levels of burning in the Northeast and Midwest and increased burning in the Upper Lake States and Central Plains, with no decline as a result of the Little Ice Age. The fossil charcoal trend here also indicates that there was no significant decline in burning in either the northern or southern zone, which would be expected if changing Native American populations controlled charcoal levels. The absence of a decline is also somewhat surprising given that the general ubiquity of reduced biomass burning observed elsewhere in the Americas has a result of simultaneous depopulation and cooling climate (Power et al. 2012; Stambaugh et al. 2013). In southeast New York, for example, transitioning to the LIA resulted in pine forests converting to spruce and hemlock and a corresponding decline in paleocharcoal relative to during the Medieval Warm Period (Pederson et al. 2005). Variations in temperature or drought, however, are not a good explanation for the pre-European settlement increase in biomass burning in the southern zone because the paleoclimate data suggest that conditions became progressively cooler after 1100 AD (Fig. 1), which would not explain increased fires. This suggests a disconnect whereby human burning often trumped climate, even in the face of depopulation (Abrams and Nowacki 2015). Moreover, it has been argued that few people can be responsible for burning large areas (Kay 2007; Bond and Keeley 2005; Pinter et al. 2011). Climate warmed after 1600 before an abrupt decline in the early 1800s without much noticeable impact on charcoal abundance. This was followed by another period of accelerated warming after ca. 1840 that persists to the present day. The abrupt and sustained increase in charcoal starting in the early 1700s predates the latest warming phase but is highly consistent with the increase in human population and associated activities (cf. Parshall and Foster 2002; Munoz et al. 2010; Marlon et al. 2013). The large decrease in fire after 1940, mainly attributed to active fire suppression via the Smokey Bear campaign (Abrams 2010), occurred during a significant warming period and may represent another important fire-climate disconnect. However, the decline of burning may have been facilitated by lessening frequency and intensity of drought after 1930 (Fig. 1; McEwan et al. 2011). The impact of severe drought and drought lessening is apparently quite important in the ecological history of the eastern USA. In the Big Woods of Minnesota, a transition from savanna woodlands to closed deciduous forest has been attributed to early (ca. AD 1300) drought reducing fuels and the impact of fire (Shuman et al. 2009). Pulse of tree recruitment at the subcontinental scale during droughts of the late 1600s was severe enough to open large canopy gaps in the eastern broadleaf forests and impact forest composition (Pederson et al. 2014).

The resilience of eastern forests that followed the clearcut and catastrophic fire era, as indicated by the results of this study, reinforces the disturbance attributes for most of its species (Nowacki and Abrams 2008). Nevertheless, there was the near eradication of chestnut, which was once an associate of oak along the Appalachian Mountains from the chestnut blight (Cryphonectria parasitica; Whitney 1994). However, generally low levels of chestnut were recorded in the pollen and witness tree records in this study, and its relative pollen production is unknown (Table 2; Steve Jackson personal communication October 2014). It has been suggested that the abundance of chestnut might have been overemphasized in the past, especially at regional scales (Hanberry and Nowacki 2016).

We recorded a large decline in pine, hemlock, and larch (in the sub-boreal forest), called “deconiferization,” which is attributed to intensive logging and their inability to reproduce sprout vegetatively (Abrams 2001; Nowacki and Abrams 2008). The Great Lakes region experienced a large increase in aspen, an aggressive pioneer tree species, as a direct result of clearcutting and burning the original conifer-northern hardwood forests (Nowacki and Abrams 2008). Large increases also occurred in maples, with cool-based sugar maple being the principal maple species regenerating within mesic conifer-northern hardwoods and warm-based red maple within more xeric pine-northern hardwoods (Nowacki and Abrams 2015). While evident in our witness tree data, the apparent lack of this response in pollen data is likely due to maple having low-pollen presence in sediment records (Prentice and Webb III 1986). Moreover, the paleochronologies used here had an average ending date of 1966 and thus often do not capture the last half-century of forest change.

The more aggressive of the maple species, particularly over the last 50 years, is red maple which is now a dominant tree over much of the eastern USA (Abrams 1998). The dramatic rise of red maple has been attributed to a suite of factors, including the extensive and repeated logging during the late nineteenth and early twentieth centuries, the loss of chestnut, the suppression of understory burning, and increased precipitation (Abrams 1992; Abrams 1998; Fei and Steiner 2007; McEwan et al. 2011). We can add increase in temperature as a driver for this warm-adapted maple as compared with sugar maple. The increase in red maple is believed to be one of several important factors responsible for oak decline in the eastern USA, which is mostly evident south of the TZL in this study. The coincidence of active fire suppression, a warming climate, and lessening of drought during the twentieth century make moisture attributions difficult. For example, there is little evidence in the paleoecology literature that maple (including red maple) significantly increased to become a forest dominant during past wet periods prior to twentieth century (Gajewski 1987; Russell et al. 1993; Foster et al. 2002; Pederson et al. 2005). In contrast, tree survey studies indicate that red maple responded quite quickly to fire suppression after 1940 (Larsen 1959; Lorimer 1984; Abrams 1992; Shumway et al. 2001). Following European settlement, the pollen and tree survey data used here revealed changes in the ecophysiological environment indicative of an overall cooling and loss of pyrogenicity and drought and increase in shade tolerance. This is consistent with a mesophication of forests in the eastern USA, which has been primarily attributed to the decline of burning after 1940 (Nowacki and Abrams 2008, 2015). The overall cooling is particularly noteworthy because forests are becoming more mesic despite general warming over the period of study.

In conclusion, paleoecological (pollen and charcoal), tree census, land use history, and climate data used in concert in this study provided a more robust interpretation of historical ecology than possible with only one data type by compensating for inherent weaknesses in each data type. The charcoal, human population, and climate data provided important support for the anthropogenic fire hypothesis over the last 2000 years in the eastern USA, including multiple instances of climate-fire disconnects that may be best explained by the role of human-caused burning.