- Open Access
Potential carbon loss associated with post-settlement wetland conversion in southern Ontario, Canada
Carbon Balance and Management volume 13, Article number: 6 (2018)
Natural wetlands can mitigate ongoing increases in atmospheric carbon by storing any net balance of organic carbon (peat) between plant production (carbon uptake) and microbial decomposition (carbon release). Efforts are ongoing to quantify peat carbon stored in global wetlands, with considerable focus given to boreal/subarctic peatlands and tropical peat swamps. Many wetlands in temperate latitudes have been transformed to anthropogenic landscapes, making it difficult to investigate their natural/historic carbon balance. The remaining temperate swamps and marshes are often treated as mineral soil wetlands and assumed to not accumulate peat. Southern Ontario in the Laurentian Great Lakes drainage basin was formerly a wetland-rich region that has undergone significant land use change since European settlement.
This study uses southern Ontario as a case study to assess the degree to which temperate regions could have stored substantial carbon if it had not been for widespread anthropogenic land cover change. Here, we reconstruct the full extent and distribution of natural wetlands using two wetland maps, one for pre-settlement conditions (prior to 1850 CE) and the other for modern-day patterns of land use (2011 CE). We found that the pre-settlement wetland cover decreased by about 56% with the loss most significant for marshes as only 11% of predicted pre-settlement marshland area remains today. We estimate that pre-settlement wetlands held up to ~ 3.3 Pg of carbon relative to ~ 1.3 Pg for present-day (total across all wetland classes).
By not considering the recent carbon loss of temperate wetlands, we may be underestimating the wetland carbon sink in the pre-industrial carbon cycle. Future work is needed to better track the conversion of natural wetlands globally and the associated carbon stock change.
Natural wetlands can act as a long-term terrestrial carbon sink by storing peat (plant organic litter undergoing very slow decomposition) under waterlogged conditions . Also, as an active methane (CH4) source, the past expansion of global wetlands has been linked to Holocene carbon dioxide (CO2) and CH4 fluctuations [2, 3]. However, poor quantification of paleo-wetland extent and carbon balance limits robust coupling . Significant wetland conversion by human land use in the heavily populated temperate zone started earlier than the industrial era in many regions and must have reduced the natural wetland carbon storage since that time [5, 6]. The implications of anthropogenic wetland conversion on carbon stocks need to be more explicitly quantified at the global scale, but regional studies are a critical first step (e.g., [7,8,9,10]).
The majority of the global peatland carbon pool is contained within northern boreal and sub-arctic regions located poleward of 45°N. In the form of bogs and fens, these regional wetlands are usually characterized by thick and extensive peat deposits, either moss-covered or sparsely vegetated with shrubs and herbaceous plants. With increasing concerns about future climate change and the fate of such large organic carbon pools in climatically sensitive areas, efforts to improve knowledge about global peatlands have centered upon northern boreal/subarctic regions ( and ref. therein). The perceived importance of northern peatlands has resulted in less emphasis on to other wetland types of mid-latitude regions, such as swamps and treed-fens, which containing larger and more abundant trees compared to the northern open fen classification (e.g., ). While the environmental conditions associated with tree growth (e.g., improved soil aeration and nutrient-rich water supply) seem less optimal in terms of long-term preservation of organic matter, some recent studies suggest comparable peat accumulation and high organic carbon densities of the tree-covered temperate wetlands [12,13,14].
Since many wetlands were drained before inventories of carbon pools took place, studies are limited by uncertain pre-disturbance extent and lack of carbon density measurements for pristine conditions as indicated by Bridgham et al. . Southern Ontario (Canada), particularly the southwestern portion, is notable for its dramatic loss of natural wetlands since European settlement [15, 16]. Some wetlands survived and have been the focus of paleoenvironmental studies, mostly with peat cores containing Holocene-age basal sediments [17,18,19,20,21,22,23]. Analyses of these wetland cores indicate considerable peat thicknesses and organic matter densities (see Fig. 1 for locations), however there has never been an effort to systematically quantify net carbon uptake. In the 1980s, the still relatively abundant wetlands in southeastern Ontario (SEO) were investigated for potential peat extraction and resource development . The SEO peatland survey reported significant peat accumulation from many regional swamp and freshwater marsh sites, prompting a reconsideration of these systems, not as overall mineral-soil wetlands (e.g., [6, 7]), but as part of northern peatlands. Studies compiling data from soil cores suggest that various wetland sites other than bogs and fens are also important in wetland organic carbon stock accounting [9, 10].
In this study, we aim to reassess southern Ontario wetlands as an underestimated middle-latitude peat and organic carbon stock. Maps of pre-settlement large (> 10 ha) wetlands and current detailed wetland distributions were combined to reconstruct the extent and distribution of wetlands in five categories (tree swamp, shrub swamp, marsh, fen, and bog) before anthropogenic land use conversion since 1850 CE. The resulting natural wetland cover represents the available potential for Holocene carbon storage, as it integrates carbon density data for the different wetland classes from the SEO peat resource survey . Decreased natural wetland cover corresponds to anthropogenic land carbon emission over the settlement period. This study estimates wetland carbon loss by first categorizing current land uses for historic wetlands. By doing so, we show that the treed wetlands (swamps) and marshes which were so prevalent in the study region prior to 1850 CE, have the potential to store large amounts of organic carbon, wetland types that have been overlooked in analyses of Holocene peatland carbon stock.
We focus our study on the terrain falling south of the Precambrian shield in Ontario, Canada (Fig. 1). Based on the ecological land classification (ELC) of Ontario , our study area belongs to the Mixedwood Plains Ecozone. This Ecozone is composed of two Ecoregions (6E: Lake Simcoe-Rideau; 7E: Lake Erie-Lake Ontario), in turn sub-divided into Ecodistricts, of which we include 6E−1 to 6E−16 (excluding 6E−3 of Manitoulin Island) and 7E−1 to 7E−6 in our study. The climate is slightly different between 6E (annual mean temperature, 4.9–7.8 °C) and 7E (6.3–9.4 °C) but overall the mildest in Canada and humid due to proximity to the Great Lakes. Growing season length exceeds 200 days (up to ~ 240 days) and precipitation ranges from 720 to 1000 mm per year .
The Canadian Wetland Classification system classifies most part of southern Ontario as ‘Eastern Temperate Wetlands’ . Dominant wetland types are hardwood swamp with maple trees (Acer saccharinum) and freshwater marsh with cattails (Typha spp.). Those wetlands were naturally established as a part of the post-glacial landscape, but drainage and conversion to croplands have been widespread since European settlement . This study contrasts two periods of natural wetlands: one before European settlement (‘pre-settlement wetlands’ in 1800 CE) and the other with modern-day land use (surveyed in 2011).
Map overlay analysis
In the absence of historical wetland maps for pre-disturbance conditions, soil maps can be used as an alternative as wet soils are proxies for wetland locations. This was the approach used by Snell  (Additional file 1: Figure S1 for details). Ducks Unlimited Canada (DUC) , building upon the work of Snell , digitized the soil maps using GIS software, and added topographic and hydrological data to better approximate the potential distribution of wetlands in the pre-settlement period.
Wetland conversion was identified by areas where the pre-settlement wetlands overlap with non-wetland areas in the modern day (2011 CE) land use map (Additional file 1: Figure S1). In Snell’s work, maps for 1967 and 1982 land uses were superimposed on the soil maps, and later the digitized analyses of DUC  also included recent trends between 1982 and 2002 using the provincial land cover information system [Southern Ontario Land Resource Information System (SOLRIS) Version 1.0]. In our study, we utilize the most recent version of SOLRIS (V2.0) for the year 2011; the modern-day (2011) wetland map was obtained by extracting wetland polygons in five classes: tree swamp, shrub swamp, marsh, bog, and fen (Additional file 1: Table S1).
Because obtaining the maximum potential area of Holocene natural (pre-settlement) wetlands and carbon storage was our objective, the 2011 wetland map was used to adjust the DUC’s pre-settlement map. By overlaying the 2011 wetland layer (SOLRIS V2.0) to the DUC’s map layer, we assigned the five wetland classes of SOLRIS V2.0 on the pre-settlement wetland polygons for more explicit carbon stock estimation by wetland class. Each pre-settlement polygon had spatial proximity with type-classified polygons from the 2011 wetland map, and the spatially closest one was allocated to the same wetland class. Central to this methodology is the assumption that the current wetlands are fragments of the larger past wetland. Our pre-settlement layer also included small wetlands (< 10 ha), which were excluded from the previous report . These smaller wetlands were excluded from the DUC report due to the coarser map scale used in the data of Snell . The smallest wetland in our pre-settlement map was 0.5 ha for the polygons from the SOLRIS layer (the minimum mapping unit in the source data), and 1 ha from the DUC layer (determined to avoid slivered polygons; larger features assumed to be wetland features by the original map overlay analysis by DUC). Overall, wetlands found in the current land cover but not in the pre-settlement layer (summing up to ~ 4656 km2 in extent), mostly from those small wetlands, were added to our pre-settlement wetland map. It was unlikely that this additional area (~ 4656 km2) was caused by creation of new wetlands since the onset of the settlement period and land use expansion.
By deleting the areas recognized as modern-day wetlands from the full pre-settlement cover, the map layer for ‘converted wetlands’ was created. Then, the converted wetland layer was overlaid onto the SOLRIS V2.0 (2011) land cover layer to display the wetland conversion. Different land uses were summarized by each ecodistrict and by our five wetland classes.
Wetland carbon stock estimation
The pre-settlement wetland map represents the maximum capacity of natural carbon held by wetlands for the Holocene by assuming the cumulative increase since the deglaciation (e.g., [2, 3, 11]), and modern-day map is comparable to how much has been altered by human land use change (e.g., pre-settlement wetland converted into farmland). Wetland carbon stocks were estimated for both pre-settlement and post-settlement (modern-day) periods using the cumulative carbon mass data obtained from the SEO peat survey (using peat depth, bulk density, and carbon content; see Additional file 1: Appendix S1 for detailed procedures).
Pre-settlement, converted, and current wetlands
From the map overlay analysis, three wetland maps were created for southern Ontario (Fig. 2) with classification into five wetland classes (Additional file 1: Table S1 for detailed information). The total area of the study region is 83,810 km2 (Fig. 1). Of this total, the extent of pre-settlement wetlands was estimated at ~ 24,984 km2 (~ 30% of total land cover), reflecting the sum of the current wetlands (~ 11,032 km2) and converted wetlands (~ 13,953 km2) as in Table 1. Among the five wetland types, marsh is the largest portion of the pre-settlement cover (Table 1 and Fig. 2a), most notably the extensive marshes of southwestern Ontario (currently agricultural lands, Fig. 2c), and the extant marshes south of Lake St. Clair (still existing as marsh, Fig. 2b). Tree swamp is now the most frequent wetland class (Fig. 2b) and occupies 75.6% of the total current wetlands by areal extent, compared to 40.5% of the total historical wetland extent in Table 1.
Agricultural land use is the most common driver of the pre-settlement wetland conversion regardless of the wetland type (Table 1), but there are some local exceptions. In the densely populated areas along the coast of Lake Ontario (7E−3 and 7E−4 in Fig. 1), conversion of natural wetlands is due mostly to urban development rather than agriculture (Fig. 3c). Some ecodistricts with relatively small wetland loss (e.g., 6E−9 at the boundary of the Precambrian Shield; Fig. 1) have a higher proportion of transition to other natural land types (Fig. 3d) mostly attributable to afforestation (Table 1). The wetland-to-forest conversion may include the current swamp not recognized by SOLRIS mapping and possibly the past wetlands where lowering of water levels promoted increasing tree cover (see Additional file 1: Table S1 for tree swamp classification).
Southern Ontario wetlands carbon storage
The parameters for wetland carbon stock calculations were obtained from the SEO peatland survey and presented in Tables 2 and 3 (see Additional file 1: Appendix S1 for data acquisition and calculation). The average peat depths by different wetland types (Table 2) were applied to convert the organic carbon density (Table 3) to the cumulative carbon mass for each wetland type (Table 4).
Carbon stock changes were calculated for southern Ontario wetlands for the two periods, the pre-settlement 3.3 ± 1.2 Pg C (1 SD) and the current 1.3 ± 0.4 Pg C (1 SD) resulting in the difference of 1.9 ± 1.3 Pg C (1 SD) as the possible carbon loss due to wetland conversion (Table 4). The carbon loss contribution by each wetland class was proportional to the extent of conversion, and hence most significant change in carbon storage is in the conversion of marshes that explain 82% of the estimated total. However, wetland classes vary in terms of carbon storage potential. In Table 4, for example, organic carbon mass per unit area of coniferous swamp (180 ± 104 kg C m−2, 1 SD) is estimated almost twice that of deciduous swamp (94 ± 76 kg C m−2, 1 SD). Thus, a critical step in obtaining accurate estimates of potential wetland carbon storage is the assignment to correct wetland classes, as each class has different carbon storage potential. In this study, if the wetlands were not categorized and the single average carbon mass ‘(All)’ multiplied to ‘Total’ wetland extent, the resulting carbon stock values and uncertainties (1 SD) are larger than the sum of classified estimates for both pre-settlement and current wetlands (Table 4).
Reconstruction of pre-settlement wetlands
Our effort to include small wetlands (< 10 ha) in our pre-settlement reconstruction increased the total estimated extent of pre-settlement wetlands by 25% compared to previous work [15, 16]. Our methodology results in a lower estimate for percentage wetland loss; for example, 72% wetland loss was reported in DUC  versus our value of 55.8% in Table 1. Despite our attempt to include small wetlands (1–10 ha) that were excluded by previous studies [15, 16], we may still be underestimating the amount of ‘small’ wetlands present historically. Small wetlands are clearly important for estimating overall wetland extent . If those small wetlands (< 1 ha) could be fully mapped for pre-settlement wetlands, the loss ratio might be higher than the estimated 55.8%.
In the pre-settlement reconstruction, the marsh wetland class is overwhelmingly dominant and swamp wetlands are proportionately much smaller; both observations contrast to modern-day (Fig. 2). In the absence of actual historical maps of wetland extent, some larger-scale marshes from the past can be determined from the region’s settlement history, such as the Holland Marsh (now croplands), Long Point (conservation and recreational use), and Toronto Harbour (infilling and railway construction) . In our reconstruction (Fig. 2), the Holland Marsh (near Lake Simcoe) and Long Point (Lake Erie) areas remain as marshes. Also, the southeastern shoreline of Lake St. Clair was known for its extensive marsh cover prior to the intensive drainage and agricultural expansion since the late 1800s . Following drainage and conversion, some small fragmented areas of marshland may have continued to exist in the region [26, 32], with a relatively extensive marsh remained along the lakeshore of Lake St. Clair (Fig. 2b). From the workload perspective of early European settlers, marshes would have been easier to convert to cropland than forested swamp landscapes simply due to the reduction in heavy woody plants. Indeed, the surrounding treed swamps of Holland Marsh farmlands (Fig. 2 insets) have remained almost untouched since their discovery (e.g., the northwestern ‘tamarack swamp’ in an early surveyor’s note in 1804; ).
The transition to the current swamp-dominated landscape may also be associated with natural wetland succession that happened in this region [18, 26, 33] and not with selective anthropogenic disturbance alone. The abundance of swamps may be self-perpetuating because once trees become dominant within the wetland, they can act to stabilize the system; disturbances by beavers notwithstanding [33,34,35]. The tendency for trees to eventually dominate a wetland landscape can be observed in SEO peat core data, where successional shifts are documented from sedge to woody peat through time (Additional file 1: Appendix S2 for further discussion).
Therefore, current dominance of swamps likely resulted from both natural wetland succession and anthropogenic disturbance. This change implies that the wetland conditions identified at the time of modern-day investigations do not necessarily reflect past wetland distribution nor the associated ability to sequester carbon. Due to the lack of information on historical conditions and peat accumulation, potential peat deposits and cumulative organic carbon stores should not be calculated based only on sparse observations of the more typically studied peatland types (e.g., open bogs and fens).
Wetland conversion and carbon loss potential
Our results suggest that more than half of the natural wetlands have transformed to other land types, eventually forcing up to ~ 1.9 Pg C of organic carbon out of the original waterlogged systems (Table 4). This amount supports an underestimation of soil organic carbon (SOC) loss due to anthropogenic land use, especially from regions found to have had abundant wetlands in the past. Due to possibly much higher organic matter content in the form of peat and its vulnerability to degradation (induced by drainage, oxygenation, and compaction), cumulative carbon loss by direct use of wetlands (up to carbon mass values in Table 4 if the entire peat layer has been extracted) would far exceed the upper limit of worldwide SOC loss by human land use (54 Mg C ha−1 in 2-m depth; 95% confidence interval) in the prediction model of Sanderman et al. . The ~ 1.9 Pg of carbon stock, only from our study region’s wetland loss, is comparable to top ten listed SOC losses by country, in which Canada is not included (in supplementary information of ).
Modernization of the current landscape is estimated to have occurred between 1835 CE to 1967 CE in southern Ontario . Over this same period, fossil fuel emissions within Canada are estimated to be 2.15 Pg C . This industrial carbon emission is estimated without considering past drainage of wetland systems, therefore the value may be significantly underestimated . According to our results, the amount of land carbon released due to the settlement and wetland conversion in southern Ontario would almost double this national fossil carbon emission . Other anthropogenic disturbances related to land-use change such as deforestation may have also resulted in a release of terrestrial carbon to the atmosphere (e.g., ).
However, the rate of carbon loss associated with wetland conversion will ultimately depend on the type of anthropogenic disturbance (i.e., gradual or abrupt). Drainage of wetlands (gradual) is far more important than peat extraction (abrupt) in southern Ontario because the peat was found to be mostly unsuitable for use as fuel or horticultural resources . Instead, agricultural use of wetland surfaces (i.e., requiring drainage) was deemed a more productive and affordable land use alternative . A very small proportion of the pre-settlement wetland removal is in fact attributable to direct peat extraction (1.8% of the total ‘Developed’ in Table 2).
Wetland drainage and row crop harvesting began during the settlement period, and are generally represented as ‘Tilled’ land use classification (tillage or tilling for soil management; cf. ‘Tile drainage’ in Fig. 1 is a common method to lower the water table prior to tillage or other land uses). Wetlands managed for farming gradually lose the surface peat by oxidation and compression, resulting in overall subsidence of the land surface. Soil cores from a historic fen in Switzerland, with similar agricultural management over ~ 140 years, recorded 16–49% loss of ‘pre-drainage’ organic carbon (‘Table 2’ of ). Although Leifeld et al.  caution against applying a single default value all the southern Ontario wetland-to-cropland area is assumed here to have experienced the same level of subsidence and peat oxidation. This approach was taken due to lack of detailed information on rates of subsidence or oxidation for our region. Thus, we applied a ross rate of 49%, as was done for the Staatswald site  based on the similar pre-drainage carbon mass of southern Ontario marshes (142 kg C m−2 as in Table 4) and the high subsidence rate (~ 3.3 cm year−1) recorded from the Holland Marsh site in Ontario, which remains one of the only sites where detailed subsidence rates are available for drained wetlands in the study region .
Accordingly, the ~ 1.9 Pg of estimated carbon loss can be partitioned into ‘completed’ (c) and ‘potential’ (p) loss. Proportional to the 2011 land classes of the converted wetlands (Table 1), 57.0% of the total disturbed wetland carbon is under agricultural land use (1.1 Pg C), 7.4% developed (0.14 Pg C), 10.2% natural (0.20 Pg C), and 25.5% unspecified lands (0.48 Pg C). Among these, the developed wetland is likely accompanied by peat extraction and completed the carbon loss (Dc: 0.14). If all the drained peats under the croplands have now reached the 49% carbon loss by subsidence, 0.54 Pg C has been released (Agriculture completed, Ac: 0.54), leaving 0.56 Pg C in soils either translocated downwards or mixed with mineral substrates (Agriculture potential, Ap: 0.56). The former wetlands now in the “natural cover” class (10.2%) most likely retain the wetland accumulated carbon relatively intact (Np: 0.20). The 0.48 Pg C under the undifferentiated lands is unknown with respect to carbon release, so here it is considered as ‘potential’ loss (Up: 0.48). In summary, among the total ~ 1.9 Pg C, only 0.7 Pg C (0.54Ac + 0.14Dc) can be said to have been ‘completely’ lost, and 1.2 Pg C (0.56Ap + 0.20Np + 0.48Up) needs further investigation for possible on-going or future loss.
Underestimated Holocene peat accumulation
A comprehensive database for Canadian peatlands (Tarnocai database hereafter) has been available since its first release date in 2005, with revisions made in 2011 . The total peatland carbon stock in that database has been used in synthesis studies of northern peatlands [3, 11]. We compare our results to the information given for southern Ontario in the Tarnocai database. We argue that all wetland classes have the potential to store considerable carbon in the form of peat, so our carbon values represent the upper range of wetland peat carbon stock for the region. According to the Tarnocai database, southern Ontario has a total organic carbon content of 0.27 Pg C, a value much smaller than our estimate of ~ 1.3 Pg C (Table 4). The reason that our estimated value of carbon is significantly higher than the value indicated in the Tarnocai database has to do with differences in spatial extent of different wetland types. In other words, the Tarnocai database uses similar range of carbon mass per area (171.7 ± 170.9 kg m−2, mean ± 1 SD, not specified for peatland types) as we use in our study (Table 4) indicating the difference has to do with how much more land is covered by different wetland types. Marsh and swamp in our estimates were one tenth the spatial extent indicated in the Tarnocai database, while fen and bog in the Tarnocai database were 8.3 and 3.4 times greater than our estimates, respectively.
Future research refining the extent of swamps and marshes, as well as their peat accumulating abilities, will help to decrease discrepancies between different carbon accounting databases. For example, in the Tarnocai database  swamp-rich regions like southern Ontario might be underrepresented for post-glacial peat accumulation and the Holocene wetland carbon sink. Future research should focus on using detailed substrate profiles to identify ‘past peatland conditions’ from apparently non-peat forming wetlands. Because often peat is not floating at the surface in swamps and marshes, these wetland types are not fully accounted for in the carbon balance database. Sediment coring is needed to determine if peat is in fact held in below-surface sediments (e.g., [18, 22, 42]). Our study which uses the SEO wetland peat cores to distinguish different carbon storage potentials by wetland class highlights the likeliness that total Holocene peat carbon accumulation has been underestimated in previous studies across this region. We extend this statement to include all middle latitude (temperate) regions in the Northern Hemisphere as we believe swamp and marsh sites have not been fully considered in published global-scale syntheses (e.g., [3, 5, 6, 11]).
How can we increase our knowledge of the extent and peat accumulating potentials of swamps and marshes in the past? Paleoecological studies may provide some insight into long-term carbon accumulation in swamp and marsh peats (i.e., see Fig. 1 for locations). Some Great Lake coastal marshes exhibit 1–2 m of organic-rich sediment layers (or peat) found between strata of mineral sediments, indicative of high energy environments [22, 43]. Bunting et al.  investigated Oil Well Bog in Ontario with shallow peat cores (~ 30–50 cm) for late Holocene wetland history; the study site where the authors suggest “the official name is misleading” as according to the vegetation found here, the site would be better classified as a swamp than a bog. Yu et al.  report on a peat core from a hardwood swamp near Rice Lake, Ontario. The sediment profile exhibited a peat layer ~ 2 m in thickness, but the top 50 cm was a silt-dominant mineral layer; the swamp is not currently accumulating peat but did so in the past.
Ott and Chimner  examined some coniferous swamp sites in Michigan and Minnesota (USA) and discovered significant woody peat accumulation in situ. These systems have not been included in the data synthesis for northern high-latitude peatlands like Loisel et al. . In northern Ontario, the vast pristine forest was studied for its carbon stock behavior in the past and under future climate projections ; that study noted that belowground carbon was conservatively estimated as peaty forest soils have not been fully accounted. The treed-wetlands within the forest (~ 30% of the study area of ) need to be studied for peat deposits for Holocene long-term carbon sink potential like the extensive fen-bog peatlands dominant in more northern latitudes . A preliminary examination for peat-rich sediments may help locate overlooked peat-forming swamps from extensive forested areas. For example, Beamish  used penetration of airborne gamma radiation to distinguish forest peat deposits from mineral soils.
While the current study examined a small region, the approach presented here can be applied to other regions to better quantify the role of global wetlands in Holocene land carbon history. To fully consider the potential for natural wetlands to contribute to past land carbon sinks, our study shows that current surface conditions cannot be used to estimate past wetland carbon stocks. Also, many forested wetlands and marshes are likely to be underestimated in regional carbon inventories as they are often considered coeval with upland forest systems or mineral soil wetlands, with below ground carbon stocks considered much less significant than typical peatlands (bogs and fens).
Temperate wetlands have experienced high anthropogenic pressures, prior to efforts to document their important contribution to land carbon sequestration. Proper assessment of the converted wetlands can help recognize the overlooked portion of Holocene peatland carbon sink. Mid-latitude peat deposits have possibly acted as significant land carbon sink under natural conditions but have been recently transformed into carbon sources by wetland drainage and anthropogenic land use. There also remains the potential for on-going release of old wetland carbon as underlying peat continues to oxidize and wetlands continue to be degraded.
ecological land classification
Ducks Unlimited Canada
Southern Ontario Land Resource Information System
Gorham E. Northern peatlands: role in the carbon cycle and probable responses to climatic warming. Ecol Appl. 1991;1:182–95.
MacDonald GM, Beilman DW, Kremenetski KV, Sheng Y, Smith LC, Velichko AA. Rapid early development of circumarctic peatlands and atmospheric CH4 and CO2 variations. Science. 2006;314:285–8.
Yu Z, Loisel J, Brosseau DP, Beilman DW, Hunt SJ. Global peatland dynamics since the Last Glacial Maximum. Geophys Res Lett. 2010;37:L13402.
Finkelstein SA, Cowling SA. Wetlands, temperature, and atmospheric CO2 and CH4 coupling over the past two millennia. Glob Biogeochem Cycles. 2011;25:1002.
Armentano TV, Menges ES. Patterns of change in the carbon balance of organic soil-wetlands of the temperate zone. J Ecol. 1986;74:755–74.
Maltby E, Immirzi P. Carbon dynamics in peatlands and other wetland soils regional and global perspectives. Chemosphere. 1993;27:999–1023.
Bridgham SD, Megonigal JP, Keller JK, Bliss NB, Trettin C. The carbon balance of North American wetlands. Wetlands. 2006;26:889–916.
Euliss NH, Gleason RA, Olness A, McDougal RL, Murkin HR, Robarts RD, et al. North American prairie wetlands are important nonforested land-based carbon storage sites. Sci Total Environ. 2006;361:179–88.
Ausseil AGE, Jamali H, Clarkson BR, Golubiewski NE. Soil carbon stocks in wetlands of New Zealand and impact of land conversion since European settlement. Wetl Ecol Manag. 2015;23:947–61.
Nahlik AM, Fennessy MS. Carbon storage in US wetlands. Nat Commun. 2016;7:13835.
Loisel J, Yu Z, Beilman DW, Camill P, Alm J, Amesbury MJ, et al. A database and synthesis of northern peatland soil properties and Holocene carbon and nitrogen accumulation. Holocene. 2014;24:1028–42.
Cai S, Yu Z. Response of a warm temperate peatland to Holocene climate change in northeastern Pennsylvania. Quat Res. 2011;75:531–40.
Ott CA, Chimner RA. Long-term peat accumulation in temperate forested peatlands (Thuja occidentalis swamps) in the Great Lakes region of North America. Mires Peat. 2016;18:1–9.
Sleeter R, Sleeter BM, Williams B, Hogan D, Hawbaker T, Zhu Z. A carbon balance model for the great dismal swamp ecosystem. Carbon Balance Manag. 2017;12:2.
Snell EA. Wetland distribution and conversion in southern Ontario. Ottawa: Environemnt Canada; 1987.
Ducks Unlimited Canada (DUC). Southern Ontario. Wetland conversion analysis. Barrie: Ducks Unlimited Canada; 2010.
Dinel H, Richard PJH, Levésque PEM, Larouche A. Origine et évolution du marais tourbeux de Keswick, Ontario, par l’analyse pollinique et macrofossile. Can J Earth Sci. 1986;23:1145–55.
Yu Z, McAndrews JH, Siddiqi D. Influences of Holocene climate and water levels on vegetation dynamics of a lakeside wetland. Can J Bot. 1996;74:1602–15.
Campbell DR, Duthie HC, Warner BG. Post-glacial development of a kettle-hole peatland in southern Ontario. Écoscience. 1997;4:404–18.
Bunting MJ, Morgan CR, Van Bakel M, Warner BG. Pre-European settlement conditions and human disturbance of a coniferous swamp in southern Ontario. Can J Bot. 1998;76:1770–9.
Bunting MJ, Warner BG. Late Quaternary vegetation dynamics and hydroseral development in a shrub swamp in southern Ontario, Canada. Can J Earth Sci. 1999;36:1603–16.
Finkelstein SA, Davis AM. Paleoenvironmental records of water level and climatic changes from the middle to late holocene at a Lake Erie coastal wetland, Ontario, Canada. Quat Res. 2006;65:33–43.
Shiller JA, Finkelstein SA, Cowling SA. Relative importance of climatic and autogenic controls on Holocene carbon accumulation in a temperate bog in southern Ontario, Canada. Holocene. 2014;24:1105–16.
Riley JL. Peat and peatland resources of southeastern Ontario. Micellaneous Paper 154. Ontario Geological Survey; 1994.
Crins WJ, Gray PA, Uhlig PWC, Wester MC. The ecosystems of Ontario, part 1: ecozones and ecoregions. Technical Report SIB TER IMA TR-01. 2009.
National Wetlands Working Group. Wetlands of Canada. Ecological land classification series no. 24. Ottawa: Environment Canada; 1988.
SOLRIS Technical Team. Southern Ontario land resource information system (SOLRIS) Version 2.0: data specifications. 2015.
Riley JL. Laboratory methods for testing peat—Ontario peatland inventory project. Miscellaneous Paper 145. Ontario Geological Survey; 1989.
Serran JN, Creed IF. New mapping techniques to estimate the preferential loss of small wetlands on prairie landscapes. Hydrol Process. 2016;30:396–409.
Riley JL. The once and future Great Lakes country: an ecological history. Montreal: McGill-Queen’s University Press; 2013.
Weaver KHA, Petrie SA, Richman SE, Palumbo MD, Dyson ME, Brisco P, et al. Waterfowl and wetlands of the Lake St. Clair region: present conditions and future options for research and conservation. Unpublished report. 2015.
van Vuuren W, Roy P. Economic evaluation of wetland preservation. In: Bardecki MJ, Patterson N, editors. Wetlands: inertia or momentum. Don Mills: Federation of Ontario Naturalists; 1989.
Bunting MJ, Warner BG. Hydroseral development in southern Ontario: patterns and controls. J Biogeogr. 1998;25:3–18.
Mitsch WJ, Gosselink JG. Wetlands. 5th ed. New Jersey: Wiley; 2015.
Bunting MJ, Warner BG, Aravena R. Late Quaternary vegetation dynamics and hydroseral development in a Thuja occidentalis swamp in southern Ontario. Can J Earth Sci. 1996;33:1439–56.
Sanderman J, Hengl T, Fiske GJ. Soil carbon debt of 12,000 years of human land use. Proc Natl Acad Sci. 2017;114:9575–80.
Boden TA, Marland G, Andres RJ. Global, regional, and national fossil-fuel CO2 emissions. Oak Ridge: Oak Ridge National Laboratory; 2017. https://doi.org/10.3334/CDIAC/00001_V2017.
Gonsamo A, Chen JM, Colombo SJ, Ter-Mikaelian MT, Chen J. Global change induced biomass growth offsets carbon released via increased forest fire and respiration of the central Canadian boreal forest. J Geophys Res Biogeosci. 2017;122:1275–93.
Leifeld J, Müller M, Fuhrer J. Peatland subsidence and carbon loss from drained temperate fens. Soil Use Manag. 2011;27:170–6.
Mirza C, Irwin RW. Determination of subsidence of an organic soil in southern Ontario. Can J Soil Sci. 1964;44:248–53.
Tarnocai C, Kettles IM, Lacelle B. Peatlands of Canada. Ottawa: Geological Survey of Canada; 2011.
Beilman DW, Vitt DH, Bhatti JS, Forest S. Peat carbon stocks in the southern Mackenzie River Basin: uncertainties revealed in a high-resolution case study. Glob Change Biol. 2008;14:1221–32.
Finkelstein SA, Peros MC, Davis AM. Late Holocene paleoenvironmental change in a Great Lakes coastal wetland: integrating pollen and diatom datasets. J Paleolimnol. 2005;33:1–12.
Packalen MS, Finkelstein SA, McLaughlin JW. Climate and peat type in relation to spatial variation of the peatland carbon mass in the Hudson Bay Lowlands, Canada. J Geophys Res Biogeosci. 2016;121:1104–17.
Beamish D. Peat mapping associations of airborne radiometric survey data. Remote Sens. 2014;6:521–39.
EB, SF, and SC designed the study. PB contributed data from DUC. EB assembled datasets, performed analyses, and wrote the manuscript. All authors reviewed the results, discussed, and revised the manuscript. All authors read and approved the final manuscript.
We thank John Riley for collecting the data for the SEO peatland survey, and we acknowledge Peatland Research and Data housed at McGill University. We also thank Jim McLaughlin and Maara Packalen for helpful discussions.
The authors declare that they have no competing interests.
Availability of data and materials
Consent for publication
All co-authors have consented to publication.
Ethics approval and consent to participate
Funding for this research was provided by the Connaught International Scholarship for Doctoral Students (University of Toronto) to E. Byun and Grants from the Natural Sciences and Engineering Research Council (Canada) to S. Finkelstein and S. Cowling.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.