Abstract
Small areas of the wetter parts of southeast Australia including Tasmania support high-biomass “wet” eucalypt forests, including “mixed” forests consisting of mature eucalypts up to 100 m high with a rainforest understorey. In Tasmania, mixed forests transition to lower biomass rainforests over time. In the scientific and public debate on ways to mitigate climate change, these forests have received attention for their ability to store large amounts of carbon (C), but the contribution of soil C stocks to the total C in these two ecosystems has not been systematically researched, and consequently, the potential of wet eucalypt forests to serve as long-term C sinks is uncertain. This study compared soil C stocks to 1 m depth at paired sites under rainforest and mixed forests and found that there was no detectable difference of mean total soil C between the two forest types, and on average, both contained about 200 Mg·ha−1 of C. Some C in subsoil under rainforests is 3000 years old and retains a chemical signature of pyrogenic C, detectable in NMR spectra, indicating that soil C stocks are buffered against the effects of forest succession. The mean loss of C in biomass as mixed forests transition to rainforests is estimated to be about 260 Mg·ha−1 over a c. 400-year period, so the mature mixed forest ecosystem emits about 0.65 Mg·ha−1·yr−1 of C during its transition to rainforest. For this reason and because of the risk of forest fires, setting aside large areas of wet eucalypt forests as reserves in order to increase landscape C storage is not a sound strategy for long-term climate change mitigation. Maintaining a mosaic of managed native forests, including regenerating eucalypts, mixed forests, rainforests, and reserves, is likely to be the best strategy for maintaining landscape C stocks.
1. Introduction
Natural forest ecosystems contain about 1500 Pg of carbon (C) [1], and on a worldwide basis, forest soils contain about twice the C of the vegetation they support [2]. The proportions of C in vegetation and soil varies greatly depending on latitude, altitude, vegetation, climate, and soil [2–4], but in boreal and temperate forests, the average C content of soils and biomass is about equal (120 Mg·ha−1) [5]. Australian forests as a whole contain the lowest forest biomass and soil C for midlatitudes (mean C of 45 Mg·ha−1 for forests and 83 Mg·ha−1 for soils to 1 m depth) [2] because of climatic limitations and the relative infertility of large areas of forest-supporting soils. However, the natural forests of Australia’s wet southeast include the tallest hardwoods in the world (up to 100 m high) (Figure 1), and according to some authors (e.g., Keith et al. [6]), the aboveground C content of these ecosystems is over 2500 Mg·ha−1, which is over 50 times the mean Australian forest aboveground value (45 Mg·ha−1) and also greatly exceeds mean values for other midlatitude forests (range 32–114 Mg·ha−1; [2]), the IPCC figure (96 Mg·ha−1; [8]) for biomass C in world temperate forests, and the mean aboveground C for mountain ash (Eucalyptus regnans) forests in Victoria (246 Mg·ha−1; [7]). However, Sillett et al. [9] found that aboveground biomass C in forests similar to those described by Keith et al. [6] was at most 706 Mg·ha−1, but probably lower (438 Mg·ha−1) due to the loss of mass in decayed hollow trunks and limbs, and questioned the assertion that the tallest southeast Australian eucalypt forests are the most C-dense forests in the world, quoting measurements of North American redwood (Sequoia sempervirens) forests which contain more than 2000 Mg·ha−1 of biomass C.

(a)

(b)
A typical value for biomass C (including roots) for Tasmanian mature eucalypt forests >55 m high, with >40% crown cover (Class 1 forests in the Tasmanian classification), was obtained by Moroni et al. [10]: these contain 470 Mg·ha−1 of C. If all forests supporting all mature trees >55 m are included (i.e., all Class 1 and 2 forests), the figure is lower (387 Mg·ha−1) which is close to the 378 Mg·ha−1 figure obtained by Moroni et al. [11] in a 14-site paired comparison of aboveground biomass (trees and coarse woody debris, but not litter) under mixed forests and rainforest but four times the IPCC figure for biomass C in world temperate forests (96 Mg·ha−1 [8]). The mean figure for C in all mature Tasmanian “wet” eucalypt forests (those with a dense understorey of shrubs and/or rainforest trees) is 232 Mg·ha−1 [10].
In the northern hemisphere, it has been found that temperate old-growth forests continue to accumulate C as they mature [12–14] and protection of old-growth temperate forests has been promoted as a C sequestration policy. A similar policy has been argued for wet eucalypt forests of southeast Australia, in order to moderate climate change induced by the effects of increased concentrations of CO2 in the atmosphere. For example, Mackey et al. ([15], p. 39), in reference to eucalypt forests, proposed that “the remaining intact natural forests constitute a significant standing stock of C that should be protected from carbon-emitting land-use activities.” However, although C stocks will undoubtedly increase as young trees mature, such a policy ignores the fact that old-growth eucalypt forests in Tasmania are not the end point of forest succession [11] and that all eucalypt forests are highly susceptible to fire.
Provided there is no stand-destroying fire and there is a seed source nearby (for example, in a rainforest-dominated gully, on a shady slope, or in a nearby unburned patch of forest), wet eucalypt forests will transition to rainforests, a process described in the classic paper by Gilbert [16]. After fire, a eucalypt forest will normally establish rapidly. By the time the eucalypt forests are about 100 years old, shade-tolerant rainforest species may already be established in the understorey. As the eucalypts reach maturity, the rainforest will form a continuous understorey and a “mixed” tall eucalypt/rainforest results ([16], Figure 1.10). The eucalypts cannot regenerate under the rainforest canopy; hence, once the eucalypts reach their maximum age of 300–500 years and progressively die, they are replaced by the rainforest understorey, although eucalypt woody debris may persist for up to two centuries on the forest floor. Peak aboveground wood volumes occur when the eucalypts are tallest and contain an understorey of rainforest species (Figure 1). Rainforests, on average, contain about 45% of the aboveground wood volume and biomass C of mixed forest [11]. To halt the progression to rainforests and maintain eucalypt cover, fire must return before the eucalypt overstorey (the eucalypt seed source) dies [16, 17].
When considering C trends in Tasmanian forests, an important question to ask is whether the lower biomass C in rainforest than in mixed mature eucalypt forest also indicates lower C in the total ecosystem under rainforest. The answer to this question depends on soil C content, does it change with vegetation cover or does it remain the same? This unanswered question prompted Norris et al. [18] to comment that “the uncertainty around soil carbon represents a priority for further work.” One might assume that as mixed forests contain about twice the biomass of rainforests [11], they should return more biomass to the forest floor than rainforests and thus accumulate more soil C than soils under rainforest. Alternatively, one could argue that frequent fires in eucalypt forests deplete soil C stocks, negating any soil C increase caused by greater cycling of C to the forest floor by the maturing eucalypts of greater biomass. A complicating issue can be the accumulation of pyrogenic C in soils under fire-susceptible forests [19, 20].
Previous work on soil C under eucalypts and rainforest in Tasmania has been limited to soil profile analysis undertaken in soil surveys covering a range of soils and parent material and forest types [21–23], supplemented by two additional profiles ([24]; P. D. McIntosh, unpublished data). The surveys were not specifically designed to compare amounts of soil C under different vegetation types but contain sufficient information to enable soil C to 1 m depth under different forest types to be calculated (Table 1). Wet eucalypt forest soils contain about 50% more C to 1 m depth than dry forests, and rainforest soils contain about 50% more C than wet forests (Table 1). However, the mean rainforest soil C figure (226 Mg·ha−1) and the C/N ratio of 19 for soil to 1 m depth under this forest type are based on analysis of only three profiles in northeast Tasmania and cannot be regarded as typical. The highest mean C/N ratio (24) is found in soils under dry forests, which experience frequent understorey fires [25], and in which, charcoal is generally visible on the soil surface and in topsoils. Rainforest sites, which have probably not been burnt for at least 500 years, have topsoil C/N ratios of 15, probably a result of higher proportion of microbial residues in relation to plant remains than in topsoils under eucalypts. About half the total soil C to 1 m depth is held below 30 cm depth.
Interestingly, the mean figure of 123 Mg·ha−1 for soil C to 1 m depth under all eucalypt forests in Tasmania over a range of soil types (Table 1) is indistinguishable from the IPCC default figure of 122 Mg·ha−1 [8] for all temperate forests, higher than the “typical” temperate forest soil figure of 100 Mg·ha−1 by Lal and Lorenz [4] but lower than the c. 150 Mg·ha−1 figure for temperate broadleaf forests studied by Duarte-Guardia et al. ([26], Figure 5(d)). Previous estimates of the soil C pool in forest soils of Tasmania and mainland southeast Australia under different vegetation types by Mackey et al. ([15], p. 28), Dean and Wardell-Johnson [27], and May et al. [28]) did not take into account published data and exceeded the figures summarised in Table 1 by factors of up to 3.5 [29].
Fedrigo et al. [30] were the first to compare measured soil C in nearby rainforest, mixed rainforest-wet sclerophyll stands (“ecotone forest”), and wet sclerophyll forest (“eucalypt forest”) in the Yarra Ranges of Victoria, but only to 30 cm depth. They found no significant difference of soil C between these forest types: values were 163 Mg·ha−1 under wet sclerophyll forest and 149 Mg·ha−1 under rainforest. Unfortunately, no information was provided on whether soils were closely matched across different vegetation plots. Consequently, effects of vegetation or fire on soil C storage may have been conflated with those resulting from inherent soil differences such as depth, drainage, and parent material.
As the information summarised in Table 1 and in the literature is insufficient to answer the question of whether soil C increases, decreases, or remains the same as wet eucalypt forests transition to rainforest over time, a study was undertaken to measure soil C in paired sites under mixed forest and rainforest in Tasmania. This paper presents the new soil C data obtained in this study. It also discusses the ecological processes affecting total ecosystem C and assesses the value of wet eucalypt forests in Tasmania for mitigating the effects of CO2 emissions on climate.
2. Methods
2.1. Site Selection
At 14 site pairs, Moroni et al. [11] compared aboveground C stocks in mixed forests and rainforests representative of mature Tasmanian forest ecosystems. (Mixed forest is defined as a forest of dominant mature eucalypts (e.g., Eucalyptus obliqua, E. regnans, E. delegatensis, E. globulus, and E. viminalis) generally 40–100 m high with a rainforest understorey; rainforest lacks standing eucalypts and the canopy is dominated by species such as Nothofagus cunninghamii, Atherosperma moschatum, and Phyllocladus aspleniifolius.) These forest pairs provided a framework for a soil C comparison. In a desktop study, all 14 pairs were examined to establish which were best matched with respect to landform, slope, aspect, elevation, and geology. For a short list of pairs, these attributes, together with soil depth and stoniness, were field checked. Six pairs showing least difference between these basic soil physical attributes and which were nonstony or contained only a few stones, and which covered representative areas of mixed eucalypt forest in Tasmania and widespread soil parent rocks, were selected. A seventh mixed forest/rainforest pair, not studied by Moroni et al. [11], was selected using the same site criteria, to provide information on a soil developed on Precambrian siltstone, a rock type extensive in northwest Tasmania. No site had evidence of human-induced soil disturbance or vegetation disturbance in the form of recent fires or clearfelling. Site locations are shown in Figure 2, and brief site details are shown in Table 2.

2.2. Field Sampling
Mäkipää et al. [34] and Muukkonen [35] recognised that in boreal forests, 10–20 samples per plot or more are required in order to reasonably estimate amounts of topsoil C on an area basis, and Cunningham et al. [36] recommended at least 30 samples for accurate determination of topsoil C (0–30 cm depth) in plantations on ex-farmland in Victoria, Australia. In practice, ideal sampling strategies such as the mechanised site sampling designed for flat land in New South Wales [37] are usually modified for logistical reasons (e.g., sampling in steep country, manual transport of samples, and time and resources available). For example, despite their earlier studies, Mäkipää et al. [34] and Muukkonen et al. [35] recommended only four samples per plot in their design for soil C assessment in Tanzania. McKenzie et al. [38] recommended at least four topsoil replicates per site for soil C assessment in Australia but gave no guidance for subsoil sampling intensity. In a national survey of topsoil C in Australia [39] designed for C accounting purposes, 10 topsoil replicates were taken at each site, at each of 3 depths (0–10 cm, 10–20 cm, and 20–30 cm), and samples were bulked for each depth increment. An identical sampling approach was used in a national survey of soil C in rainforests of Papua New Guinea [40], and in addition, at each site, one sample was taken from the subsoil at 30–60 cm depth and another at 60–100 cm depth. McIntosh et al. [41] found that 10 topsoil samples per plot, separately analysed, were sufficient to demonstrate significant () differences between soil C concentrations under different vegetation types and land uses in loessial soils of southern New Zealand.
For reasons of logistics mentioned above, including the limitations imposed by the necessity for manual sampling by a small team in dense forests, sampling followed procedures developed for quantitative (per hectare) soil C measurement in forests of Papua New Guinea [40], except that samples were individually analysed rather than bulked. For topsoil sampling, an Eijkelkamp split-tube sampler taking a 48 mm diameter core was used. Ten cores to 30 cm depth in mineral soil (excluding any litter or O horizon) were sampled within a typical and undisturbed area under each forest type (mixed forest or rainforest). Each core was split into 3 depth increments (0–10 cm, 10–20 cm, and 20–30 cm) and separately bagged for later analysis. A soil pit was dug to 1 m depth or a rock contact and photographed, described and classified. Subsoils were sampled at 30–60 cm and 60–100 cm depths, or to a rock contact (sites 1 and 7). Subsoil bulk density samples were taken with an aluminium tube (75 mm internal diameter and length 100 mm), at 40–50 cm and 75–85 cm depths, where the rock contact allowed.
Samples were air-dried at about 25°C and weighed. For each depth increment, subsamples were oven-dried at 105°C and passed through a 2 mm sieve to remove stones and determine the effective bulk density (EBD) of respective <2 mm diameter fractions.
2.3. Organic Matter Analysis
Acid soil pH values and tests with HCl indicated that no free carbonates were present in any profile. At sites 1–5, total C and N were analysed by combusting finely ground 0.1 g subsamples in a Vario Macro Cube elemental analyser, and the labile C fraction was measured using a hot water extraction method as outlined by Ghani et al. [42]. At sites 6 and 7, a similar combustion method was used: analysis was performed with a Euro EA Elemental Analyser (HEKAtech GmbH, Wegberg, Germany).
Physical soil organic matter (SOM) fractionation was undertaken on sieved air-dried soil samples for two contrasting sites: Sumac Road and Eleven Road (Table 2). A composite sample (30 g) derived from the right, middle, and left soil profile wall was prepared for the uppermost (B1 or A1) and lowermost (C or Bhs) horizons at these two sites, respectively, to identify differences between forest types as well as between topsoils and subsoils (the Bhs horizon at Eleven Road is a B horizon rich in humus and sesquioxides). The fractionation procedure combined density and particle-size fractionation (Kögel-Knabner [43]). The bulk sample was homogenised in a glass beaker and sprayed with water to reduce hydrophobicity and remove entrapped air. 150 ml of sodium polytungstate (SPT) solution (1.8 g·cm−3) was added gradually to saturate the sample and finally submerge it. The suspension was left standing overnight to allow separation of lighter and heavier fractions. Free-floating particles, representing free particulate organic matter (fPOM), were removed with a vacuum pump and washed using pressure filtration and 0.22 μm membrane filters (Berrytec GmbH) until the outflow had an electrical conductivity (EC) of <2 μS cm−1. Subsequently, the fPOM was freeze-dried. The heavier fraction was subjected to ultrasonication (Sonoplus HD 2200, Bandelin, Berlin, Germany) applying 600 J·ml−1 at 70% intensity to break up soil aggregates and release occluded particulate organic matter (oPOM). To avoid strong heating and consequent mutation of released fragments, the suspension was cooled with water during this step. Samples were centrifuged (3500 rpm, 30 min), and the free-floating oPOM was removed by suction. The extracted oPOM was partitioned into >20 μm (oPOM>20) and <20 μm (oPOM<20) fractions by sieving. Both fractions were washed until outflows had an EC of 2 μS cm−1, and then freeze-dried. The mineral residue was repeatedly purged with deionized water and centrifuged (3500–5000 rpm, 30 min) until the EC of the supernatant was below 50 μS cm−1. The residue was wet sieved to obtain coarse sand (2000–630 μm), medium sand (630–200 μm), fine sand (200–63 μm), and coarse silt (63–20 μm) fractions. The outflow was collected to further separate medium silt (20–6.3 μm) from fine silt and clay (<6.3 μm) by gravity sedimentation in Atterberg-type cylinders, making use of Stokes’ law. The two size fractions were freeze-dried, and all remaining mineral fractions were dried at 65°C. Thereupon, mineral fractions >20 μm were ground with a swing mill (3 min). The fPOM and oPOM>20 fractions were ground by hand using a mortar. Processed samples were put in glass vials and stored at room temperature until elemental analysis.
The chemical composition of fPOM, oPOM>20, oPOM<20, and clay fractions was analysed using 13C-CPMAS NMR spectroscopy. A Bruker Advance III 200 Spectrometer was used, coupled with a Bruker MAS II control unit, and operated with Bruker BioSpin software. The NMR rotor was spun at a MAS spin rate of 6.8 kHz with varying recycle delay time and scan counts. (The recycle delay time refers to the time between individual scans, while scan counts refer to the total number of scans required until an appropriate spectrum could be identified.) Fractions with a relatively high C concentration, i.e., mainly fPOM, oPOM>20, and oPOM<20, for which sufficient material was available from the fractionation procedure, were scanned with a recycle delay time of 1.0 s while scan counts ranged between 3,400 and 91,400. For fractions showing a low C concentration, i.e., mainly clay fractions, or fractions containing only small amounts of material from the fractionation procedure, the recycle delay time was set to 0.4 s and up to 370,000 scans were required to obtain quality spectra. The clay fraction of the samples of the Ferralsol at Sumac Road required most scan counts: this fraction combined low C contents with high contents of pedogenic Fe. Pedogenic Fe has the potential to reduce the 13C signal and thus lead to less distinct peaks [44], which increases the number of scans required. After the samples were scanned, the spectra were processed with line broadening ranging from 0 to 75 Hz, depending on the distinctness of their peaks. Furthermore, the spectra were phase adjusted and baseline corrected (in an automatic mode). The interpretation of spectra followed principles outlined in detail by Kögel-Knabner et al. [45].
2.4. Statistical Analysis
Standard deviations were calculated for all means (Table 3). The effect of forest type on bulk soil C/N ratios in both top- and subsoil was tested based on estimation graphs as demonstrated by Gardner and Altmann [46]. The graphs show the 95% confidence interval, which is determined via bootstrap resamples (5000 repetitions) of the effect size. Bootstrap was bias-corrected to account for skew. This was accomplished using the Python (version 3.7) module DABEST [47].
3. Results
3.1. Carbon Content of Site Pairs and C/N Ratios
The soils under rainforest and mixed forest contain 110 to 326 Mg·ha−1 (Table 3; Figure 3) of total C. No statistical difference between soil C values under rainforest and mixed forest (t-test, ) was detected: mean values were 203 Mg·ha−1 in soils under rainforest and 199 Mg·ha−1 in soils under mixed forest. The latter figure is close to the mean value measured by Cotching ([50], p. 86) to 1 m depth for six soil orders under unspecified native forest in Tasmania (193 Mg·ha−1; range 94–273 Mg·ha−1) but higher than the mean soil figure obtained for wet eucalypt forests in earlier surveys (147 Mg·ha−1; Table 1) probably because the earlier surveys were not confined to mixed forest containing mature eucalypts and included very stony profiles (the stoniest profile contained c. 70% stones by volume and 37 Mg·ha−1 of C). Such profiles were not sampled in this study.

In deep soils sampled to 100 cm depth (Table 3, sites 2–6), 39–76% (mean 50%) of the soil C is held below 30 cm depth (Table 3; Figure 4). The clay-rich Styx Road Cambisols developed in Quaternary colluvium contained most C (300–326 Mg·ha−1).

There was no significant difference between C/N ratios in uppermost (A1 and B1) horizons sampled under rainforest and mixed forest; mean values were 20–22 (Figure 5). However, in lowermost horizons sampled (generally, B3, C, BCg, or Bh horizons), C/N ratios were significantly higher in soils under mixed forests (17 ± 5) than under rainforests (13 ± 6) (Figure 5).

(a)

(b)
During sample preparation, 50% of the samples (n = 102) taken from mixed forest profiles were observed to contain identifiable pieces of pyrogenic C (charcoal) in the >2 mm fraction, but only 16% of the samples (n = 92) taken from rainforest profiles contained identifiable pyrogenic C.
Mean total C pools in mixed forests and rainforests were estimated (Table 4) using the soil data of this paper and the biomass data of Moroni et al. [11] and Moroni et al. [10], who assumed that roots contain 25% of measured aboveground biomass [51, 52].
3.2. Soil Organic Matter Characterisation
The distribution of C among SOM fractions varied between sites and with soil depth. In the uppermost (B1) horizons at Sumac Road, most C was found in the clay fraction, with slightly more under mixed forest (68%) than rainforest (60%) (in Ferralsol profiles under both rainforest and mixed forest, an A1 horizon as defined by the National Committee on Soil and Terrain [53] was not present and the B1 horizon was the uppermost mineral horizon). The remaining C was largely in the oPOM form with a negligible contribution of fPOM. In contrast at Sumac Road under rainforest, 75% of total C in the C horizon was held in the oPOM<20 fraction compared to 6% in this fraction below mixed forest. This subsoil difference may reflect different SOM stabilisation mechanisms in subsoils under the two forest types, as observed in forests elsewhere [54] and in agricultural systems [55]. In the uppermost (A1) horizon at Eleven Road under mixed forest, most C was held in the oPOM<20 fraction, but under rainforest, most was held in the clay fraction. The Bhs horizons under both forest types had highest C contents in the clay fraction.
In the Ferralsols at Sumac Road (site 5), the fPOM extracted from the <2 mm B1 horizon under mixed forest had a distinct NMR peak in the aromatic aryl-C region (130 ppm), indicating the presence of pyrogenic C [56, 57] in this horizon under mixed forest; this peak was not detected in the B1 horizon under rainforest (Figure 6). In the fPOM fraction of the B1 horizon under mixed forest 28% of the total contribution was derived from aryl-C, whereas under rainforest, the respective aryl-C contribution was only 15%. In Ferralsol subsoils (Figure 7), no clear differences of aryl-C under rainforest and mixed forest were apparent.


(a)

(b)
Eluviation of pyrogenic C through profiles is confirmed by the weak aryl-C peak in NMR analyses of organic/clay coatings in the Bhs horizon of the rainforest Podzol at Eleven Road (Figure 7) which contains 3.4% C. Unfortunately the amount of fPOM and oPOM material in Podzol subsoil samples was insufficient for all organic fractions to be analysed, and usable spectra for samples under mixed forest were not obtained. The NaOH-insoluble fraction of the organic matter in the rainforest Bhs horizon (Figure 7(a)) was radiocarbon dated by accelerator mass spectrometry (AMS) to 2882 ± 24 yr BP (Wk49496, uncalibrated). In another study (P. McIntosh, unpublished data), organic matter from a podzol pan in northeast Tasmania was dated 3133 ± 38 BP (Wk17421, uncalibrated), after similar NaOH pretreatment. If it is assumed that C illuviation is a continuous process under the present climate, then these dates are the average age of C illuviation in these soils, and the subsoil horizons analysed may well be around 6000 years old or older. However, some profiles are younger: buried charcoal at 1 m depth in the Umbrisol under mixed forest in granitic colluvium under mixed forest at site 6 (South Esk Road) (Table 1 and Figure 2) was dated 1838 ± 14 yr BP (Wk49497, uncalibrated), and probably dates localised erosion, for example, fall of a large eucalypt tree with its intact root-ball of soil, following a fire.
4. Discussion
4.1. Soil Carbon Stocks
There is no evidence that soil C either increases or decreases during the transition of mixed forests to rainforests (Table 3). Soil C stocks are within the range of previous profile values (Table 1; [50]). It is most unlikely that, even if significant soil C differences were detected by more detailed and statistically robust studies, they would be sufficiently large to balance the loss of total biomass C (mean value 259 Mg·ha−1; Table 4) resulting from the transition of mixed forests to rainforests. In addition, as rainforest sites had little litter cover and litter in mixed forests was mainly in the form of bark accumulations at the base of large eucalypts, taking account of litter (estimated by Mackey et al. ([15], p. 22) to be 2% of total biomass and soil C in mature eucalypt forest) would probably accentuate the biomass C difference between rainforest and mixed forest, rather than lessen it.
The mean amount of soil C held below 30 cm depth for all profiles is 50%, demonstrating the importance of sampling full profiles to 1 m depth or a rock contact when calculating ecosystem C storage. The mean C contained at 0–30 cm depth under mixed eucalypt forests in this study (Table 3) is 95 Mg·ha−1 (range 52–144 Mg·ha−1); previously published estimates of forest soil C under mixed forests and rainforest are too high and not based on information available at the time of publication. For example, the unreferenced 271 Mg·ha−1 value quoted in 2010 by Dean and Wardell-Johnson [27] for soil C at 0–30 cm depth under tall old-growth forest in Tasmania is almost three times the figure measured in this study and over twice the highest measured forest C value at 0–30 cm soil depth in Tasmania available in 2010: 121 Mg·ha−1 for the Stronach soil (Umbrisol) formed in granite colluvium in northeast Tasmania [21]. Likewise, the unreferenced 369 Mg·ha−1 of C at 0–30 cm depth quoted by Dean and Wardell-Johnson ([27], Table 1) for rainforest soils in Tasmania is over three times too high: in 2010, the mean measured C value for 0–30 cm soils at three rainforests sites in Tasmania (Table 1) was 102 Mg·ha−1, almost identical to the mean value of 105 Mg·ha−1 measured in this study (Table 3). In addition, Dean and Wardell-Johnson [27] provided no evidence to support their assertion that C levels in soils at 0–30 cm depth soils under rainforest exceed those under old-growth (eucalypt) forest by 98 Mg·ha−1; no significant difference of 0–30 cm soil C between these two forest types was found in this study (Table 3) or in the study conducted by Fedrigo et al. [30].
The unreferenced c. 700 Mg·ha−1 “baseline” soil C figure for the “full soil profile” of undefined depth under Tasmanian E. regnans forest used by Dean and Wardell-Johnson ([27], Figure 8) is more than twice the maximum measured soil C value under E. regnans available to these authors in 2010 which was 273 Mg·ha−1 in the previously mentioned Stronach profile [21]. The 700·Mg·ha−1 value also exceeds the highest published forest soil C figure known to the authors for a forested nonswamp site in the Australasian region (600 Mg·ha−1) measured in a Papua New Guinea soil formed in gabbroic alluvium [58] and is more than three times the mean mixed forest soil C figure in this study (Table 3).
There are numerous references to total soil C in Tasmanian forests by May et al. [28] but the soil depth to which their figures refer is not specified, although one estimate ([28], p. 17) refers to “<30 cm” soils and another refers to “A and B horizons” (p. 274). Whether their average soil C figure of 280 Mg·ha−1 “assumed to represent all eucalypt forests in southeast Australia” ([28], p. 274) refers to 0–30 cm or 0–1 m depth, it is at least two times too high for Tasmanian eucalypt forests, in which the mean soil C in information available to May and co-authors in 2012 for all eucalypt forests in Tasmania was 123 Mg·ha−1 for 0–1 m depth and 63 Mg·ha−1 for 0–30 cm depth (Table 1). Consequently modelled potential C sequestration figures published by May et al. ([28], p. 54, Figure 4 and Table 5), unfortunately promoted with the aim to “develop a clear accounting framework for carbon stocks and flows in Tasmania’s forests” ([28], p. 10), are overestimated and incorrect. Similarly, the modelled losses in soil C from “old-growth” (mixed) eucalypt forests after land-use change (e.g., the 300 Mg·ha−1· C loss modelled by Dean and Wardell-Johnson [27]) based on the unreferenced 700 Mg·ha−1 baseline is, in turn, overestimated; the estimated loss exceeds the total C held in most forest soils as determined in this and previous studies (Table 1).
In their baseline map of soil C at 0–30 cm depth, designed to support national carbon accounting in Australia, Viscarra Rossell et al. [39] estimated that the largely forested lands of western Tasmania contained 161–220 Mg·ha−1 of C in soils at 0–30 cm depth. When compared to the range of C values found in this study for soils with few or no stones (52–144 Mg·ha−1; Table 3), the modelled range appears to be an overestimate, especially as large areas of forested western Tasmania are steeplands with stony soils. Consequently, the mean C figure for all Tasmanian soils at 0–30 cm depth (134 Mg·ha−1, [39]) may also have been overestimated.
The similarity of mean measured soil C stocks under mixed forest and rainforest, when contrasted with the large difference in the biomass C of the two forest types (Table 4) and likely greater biomass and C inputs to the soil under mixed forest than under rainforest, raises three questions: (1) Is the similarity of mean soil C values under mixed forest and rainforest a result of soils under rainforest not yet having attained equilibrium with their “new” lower-biomass forest cover? (2) Has soil C under mixed forest been enhanced (above levels expected in forests subject to frequent fire) by an intermittent supply of relatively inert charcoal? (3) Does the soil C measured in these studies represent the balance of C input and breakdown over millennia rather than centuries and is what we measure now the result of inherited soil C accumulation independent of present vegetation cover?
Although NMR results for a greater range of sites would be advantageous for determining organic matter processes in these soils, the NMR results obtained, together with the observation that charcoal is present in both eucalypt and rainforest soils, and the radiocarbon ages on buried charcoal, allow tentative answers to be provided for the 3 questions posed above. All three propositions appear to be correct: (1) rainforest soils have not yet attained equilibrium with their “new” lower-biomass forest cover; they still contain a chemical signature of past fires; (2) NMR results show that soil C levels under mixed forest (and to a lesser extent under rainforest) have been enhanced by intermittent supply of relatively inert (pyrogenic) carbon, as well as other stable carbon fractions; and (3) measured soil C under rainforest retains characteristics of mixed forest cover, and both the amount and character of the C in the soil under these vegetation types is probably largely inherited and has been determined by processes acting over thousands of years.
However, soil organic matter characteristics are also influenced by present-day vegetation. The greater relative amounts of pyrogenic C (and aryl-C) in Ferralsol surface B1 horizons under mixed forest and the higher mean C/N ratio of subsoils under mixed forest than under rainforest is likely to be a consequence of greater fire frequency in the eucalypt-dominated mixed forests than in the adjacent rainforests. These effects are discussed below.
4.2. Changes of Soil C after Fire
Measures of short-term change of soil C stocks after fires have yielded ambiguous results. In a global meta-analysis, Johnson and Curtis [59] noted an increase of soil C after wildfire (attributed to the sequestration of pyrogenic C and contributions from postfire nitrogen-fixing plants) but a decrease after prescribed fuel-reduction burns. Pennington et al. [60] noted significant declines of topsoil C concentration at 0–5 cm and 5–10 cm depth (from 8.04 to 5.40% and from 4.12 to 3.41%, respectively) after harvest and a regeneration burn in tall wet eucalypt forests at a southern Tasmania site. However, the difference between unburnt and burnt topsoils at 0–30 cm depth (5 Mg·ha−1) was not statistically significant once bulk density had been taken into account. No correction was made for stones, but stone content “was shown to be very low” [60]. In contrast, Slijepcevic [61] noted an average 3.5% increase of topsoil C (from 145 to 150 Mg·ha−1 at 0–10 cm depth) after harvest and burning in three wet eucalypt coupes (harvest areas). The results were not statistically analysed, but Slijepcevic ([61], p. 285) considered that they “did not provide any evidence of carbon loss or gain from the upper layers after burning.” Ellis and Graley [62] found that immediately after “hot” (intense) regeneration burns following eucalypt harvest at two locations in Tasmania, C in topsoil (0–10 cm depth) decreased by 6.4–7.4 Mg·ha−1, but long-term trends were not investigated. In Victoria, Australia, Polglase et al. [63] found that in E. regnans forests regenerating after fire, soil C reached 86% of equilibrium (steady state) values when the forest was aged 30 years, with “true equilibrium” not being reached until about 150 years. The effect of regular fires on soil properties in wet eucalypt forests was addressed by Guinto et al. [64] who measured a significant 1.9% decline in C in soils at 0–10 cm depth in wet sclerophyll forest soils burnt every two years for 20 years in southeast Queensland. Assuming a topsoil EBD of 0.5 Mg·ha−1 for these 0–10 cm soils, this 20-year decline equates to a C loss of 9 Mg·ha−1. In contrast, Krishnaraj [19] found that topsoil C under eucalypt forests of the Otway Ranges of southern Victoria increased by 2% after a forest fire due to accumulation of charcoal.
The different results obtained by these authors probably reflect the problem of adequately covering postdisturbance soil variation by the sampling strategies used. Until more comprehensive work is done, results of these case studies should not be used to generalise about the short-term effect of harvest and burning on soil C. The eucalypt forests investigated in this study must have experienced stand-destroying fires more recently than the rainforests, but the effects of fires is not evident in mean total C values.
4.3. Sequestration of Carbon by Wet Eucalypt Forest
Many authors (e.g., Keith et al. [6], Mackey et al. [15], Dean and Wardell-Johnson [27], May et al. [28], and Dean et al. [65]) have promoted the protection of eucalypt forests to mitigate the effects of climate change. The weaknesses in the arguments presented by these authors have been highlighted by Moroni et al. [67] and McIntosh and Moroni [29] who pointed out that (1) mature eucalypt forests are not the steady-state end point of forest succession and (2) that a concept originally designed to assess C sequestration potential in ecosystems (Gupta and Rao [68]) cannot be applied to landscapes as a management tool. The relevance of this study to these two issues is discussed below.
4.3.1. Forest Succession
The soil studies summarised in Table 3 and Figures 3 and 4 indicate that during the succession from mixed forest to rainforest, no change in mean soil C stocks is detectable. Given the long-term stability of the soils as indicated by NMR analyses and the lower biomass (and probably lower biomass inputs into the soil) under rainforest, and previous measurements of total C in Tasmanian profiles (Table 1), it appears unlikely that more detailed and statistically robust studies would demonstrate that the loss of biomass C (259 Mg·ha−1) during the mixed forest/rainforest transition is compensated by an increase in soil C. (The mean value of soil C measured under rainforests in this study would need to be 2.5 times higher to compensate for the biomass C loss.) Consequently, it can be concluded that the C loss resulting from the mixed forest/rainforest transition is around 260 Mg·ha−1. If we assume that wet eucalypt forests are tallest (and contain most C) when they are 400–500 years old, that all eucalypt trees have fallen after another 150 years, and that fallen trees have rotted away after another 250 years, then we can calculate that during the mixed forest/rainforest transition, 260 Mg·ha−1 of C will be released to the atmosphere over an approximately 400-year period, i.e., C is lost at a rate of about 0.65 Mg·ha−1 year−1. Thus, during natural succession, mature wet eucalypt forests are significant C emitters, not C accumulators.
4.3.2. Forest Management
Despite acknowledging that “the role of fire in natural forests is complex and must be considered on a landscape-wide basis” (Mackey et al. [15], p. 20), these authors underestimated the effect of fires on C accumulation at the landscape scale, even though eucalypt ecosystems are fire-induced. For example, Mackey et al. [15] (p. 23 and 26) excluded bushfire-affected areas when calculating gross primary production, and living biomass C was calculated only for “undisturbed” forests. The limitation of fires on C accumulation is evident in areas of Tasmania and the southeast Australian mainland affected by the 2019 and 2020 bushfires (discussed below): while not all burnt eucalypts are killed by fire, the overall effect of crown fires is to halt landscape-scale C accumulation or reverse it. Because eucalypt forests in their natural state will either burn and regenerate or transition to rainforest in the absence of fire, it is not correct to assume that mature eucalypt forests are a “baseline” against which to judge land-use change, as done, for example, by May et al. ([28], p. 274). Neither the biomass present before disturbance (whether the result of burning, harvesting, or storm damage) nor that present after disturbance represents the potential biomass in the landscape as a whole, and the true biomass potential at a landscape scale will always lie between the disturbed and undisturbed extremes. The biomass C in the landscape is always in a state of flux and cannot be defined.
Even the concepts of “undisturbed” (Mackey et al. [15]) and “natural” (Keith et al. [66]) must be used cautiously, since Aboriginal management has profoundly affected forest structure in Tasmania for many thousands of years. Early European settlers remarked on fire use by the Aboriginal population [69–72]. Frequent fires probably lit by the Aboriginal population maintained eucalypts and grasslands in upland northeast Tasmania and prevented expansion of rainforests [73]. The pollen and charcoal record in deep peats on Surrey Hills in northwest Tasmania and on the Nicholas Range in northeast Tasmania provides evidence of fires and fire-induced open eucalypt woodlands being present for at least the last 12 000 years [74–76].
We can conclude that many eucalypt forests had a more open character under Aboriginal management than they do at present, and that regular low-intensity fires (mostly burning understorey only) periodically released only modest amounts of C to the atmosphere. We can be reasonably certain that at least during the Holocene, the Tasmanian lowland and midaltitude vegetation pattern before European arrival was a mosaic of fire-induced forest, rainforest, open woodland, grassland, and sedge-dominated moorlands, a pattern noted by early surveyors such as Lawrence [77] in the Florentine Valley near Maydena and Hellyer [71] in northwest Tasmania, which in part was anthropogenically induced and which, incidentally, makes description of the Tasmanian wet forests as a “wilderness,” e.g., “Tasmanian Wilderness World Heritage Area,” a misnomer. For a brief period from approximately 1825 (by which time, the Aboriginal population was severely depleted) to 1920 (when mechanised harvest began in earnest), forest cover was neither controlled by frequent Aboriginal land management nor disturbed by clearfelling operations, although timber was extracted from more accessible sites [78]. In this period, for the first time in 12000 years, many forests previously managed by Aboriginal burning grew more densely and regenerated without large-scale human intervention, greatly increasing aboveground biomass and fuel for fires.
Possibly partly because of this post-1825 fuel increase and cessation of low-intensity burning, there have been a number of large fires since European settlement in Tasmania. The largest was in 1898 when over 1 million ha and possibly up to 2 million ha of forests in the southwest burned [28]. Other very large fires were in 1934 when 800 000–900 000 ha burned [79] and in 1967 when approximately 250 000 ha burned [80]. Recently, over 2400 recorded lightning strikes in Tasmania in December 2018 and January 2019 caused 72 vegetation fires extending over 205 000 ha including c. 100 000 ha of wet eucalypt forests and 14% (296 ha) of Tasmania’s tall (>70 m high) eucalypt forests and 7000 ha of rainforest [81].
5. Conclusions
Previous studies showed that soil C to 1 m depth (or a rock contact) under wet eucalypt cover ranged from 35 to 325 Mg·ha−1, depending on stoniness, soil depth, and parent material, but provided little data on soil C in rainforest soils, so the soil C trend (if any) during the ecological transition of mixed eucalypt forests (mature tall wet eucalypts with a rainforest understorey) to rainforests could not be established.
This study of 7 paired sites found no consistent difference between soil C stocks under mixed forests or rainforest in Tasmania, and the mean soil C value under both forest types is about 200 Mg·ha−1. Consequently, at a landscape scale, the loss of biomass C (259 Mg·ha−1) during the transition of mixed eucalypt forests (mean biomass C, 473 Mg·ha−1) to rainforest (mean biomass C, 214 Mg·ha−1) (Table 4) approximates to the total ecosystem loss (269 Mg·ha−1) during this ecological transition (Table 4).
If it is assumed that the transition to rainforest takes 400 years, then the mixed forests are emitting about 0.65 Mg·ha−1 yr−1 until they reach pure rainforest. During this ecological transition, the soils remain (on average) relatively stable: deep soil horizons under rainforest retain an imprint of pyrogenic C from past fires, and the soils appear to be buffered against the effects of vegetation change and variable C inputs.
The observation that Tasmanian tall wet eucalypt forests become net C emitters as they mature and transition to rainforests is contrary to trends monitored in northern hemisphere temperate forests in which C continues to accumulate as forests age [12–14]. It has significant negative implications for land managers expecting C credits as a result of reserving wet eucalypt forests. Although at a local scale, a young forest will accumulate C, over the long term, and at the landscape scale, there is unlikely to be a net C benefit from setting aside these forests as reserves: the unmanaged forests will either burn, or they will transition to rainforests containing less C than the eucalypt forests they replace.
Wetter Tasmanian forests have always consisted of a mosaic of eucalypt age classes and rainforests and more open areas, as a result of natural disturbance (fires and storms) and Aboriginal landscape management [70–76], and we can be confident that they have never consisted of uniform old-growth eucalypt forest. It is unlikely that the theoretical maximum percentage cover of 10–15% tall wet forests (Class 1–3 forests of Moroni et al. [10]) could ever be achieved. As previously argued by Harmon [82], Norris et al. [18], and Poynter and Ryan [83], maintaining a mosaic of managed native production forests, including regenerating eucalypts, mixed forests, rainforests, and reserves for biodiversity, soil and water, and cultural values, and also for recreational use, while at the same time ensuring good access to facilitate wildfire control, is likely to be the best strategy for maintaining landscape C stocks in the long term.
Data Availability
The survey data on which this study is based are available in the theses by Hardcastle [47] and Klöffel [48]. Calculations of C stocks from the soil data of Grant et al. [21], Hill et al. [22], and Laffan et al. [23] are available as an excel spreadsheet from the corresponding author.
Disclosure
This paper is partly based on the unpublished research theses of J. Hardcastle and T. Klöffel.
Conflicts of Interest
None of the authors have conflicts of interest regarding the conduct of this study or the interpretation of its results.
Acknowledgments
The authors are grateful to Elizabeth Brewer, Adrian Slee, and Carsten Müller who assisted with fieldwork, Peter Schad for checking FAO soil classifications, and Carsten Müller, Mark Poynter, and two anonymous reviewers for their valuable comments on drafts of the text. The Board of the Forest Practices Authority provided funds to support the research and travel costs of James Hardcastle, Tobias Klöffel, and Carsten Müller. The authors thank Sustainable Timber Tasmania and the Department of Primary Industries, Parks, Water, and the Environment for allowing access to sites.