Soil carbon modelling as a tool for carbon balance studies in forestry

Soils represent a remarkable stock of carbon, and forest soils are estimated to hold half of the global stock of soil carbon. Topical concern about the effects of climate change and forest management on soil carbon as well as practical reporting requirements set by climate conventions have created a need to assess soil carbon stock changes reliably and transparently. The large spatial variability of soil carbon commensurate with relatively slow changes in stocks hinders the assessment of soil carbon stocks and their changes by direct measurements. Models therefore widely serve to estimate carbon stocks and stock changes in soils. This dissertation aimed to develop the soil carbon model YASSO for upland forest soils. The model was aimed to take into account the most important processes controlling the decomposition in soils, yet remain simple enough to ensure its practical applicability in different applications. The model structure and assumptions were presented and the model parameters were defined with empirical measurements. The model was evaluated by studying the sensitivities of the model results to parameter values, by estimating the precision of the results with an uncertainty analysis, and by assessing the accuracy of the model by comparing the predictions against measured data and to the results of an alternative model. The model was applied to study the effects of intensified biomass extraction on the forest carbon balance and to estimate the effects of soil carbon deficit on net greenhouse gas emissions of energy use of forest residues. The model was also applied in an inventory based method to assess the national scale forest carbon balance for Finland’s forests from 1922 to 2004. YASSO managed to describe sufficiently the effects of both the variable litter and climatic conditions on decomposition. When combined with the stand models or other systems providing litter information, the dynamic approach of the model proved to be powerful for estimating changes in soil carbon stocks on different scales. The climate dependency of the model, the effects of nitrogen on decomposition and forest growth as well as the effects of soil texture on soil carbon stock dynamics are areas for development when considering the applicability of the model to different research questions, different land use types and wider geographic regions. Intensified biomass extraction affects soil carbon stocks, and these changes in stocks should be taken into account when considering the net effects of forest residue utilisation as energy. On a national scale, soil carbon stocks play an important role in forest carbon balances.

Fen ecosystem carbon gas dynamics in changing hydrological conditions

Northern peatlands are thought to store one third of all soil carbon (C). Besides the C sink function, peatlands are one of the largest natural sources of methane (CH4) to the atmosphere. Climate change may affect the C gas dynamics as well as the labile C pool. Because the peatland C sequestration and CH4 emissions are governed by high water levels, changes in hydrology are seen as the driving factor in peatland ecosystem change. This study aimed to quantify the carbon dioxide (CO2) and CH4 dynamics of a fen ecosystem at different spatial scales: plant community components scale, plant community scale and ecosystem scale, under hydrologically normal and water level drawdown conditions. C gas exchange was measured in two fens in southern Finland applying static chamber and eddy covariance techniques. During hydrologically normal conditions, the ecosystem was a CO2 sink and CH4 source to the atmosphere. Sphagnum mosses and sedges were the most important contributors to the community photosynthesis. The presence of sedges had a major positive impact on CH4 emissions while dwarf shrubs had a slightly attenuating impact. C fluxes varied considerably between the plant communities. Therefore, their proportions determined the ecosystem scale fluxes. An experimental water level drawdown markedly reduced the photosynthesis and respiration of sedges and Sphagnum mosses and benefited shrubs. Consequently, changes were smaller at the ecosystem scale than at the plant group scale. The decrease in photosynthesis and the increase in respiration, mostly peat respiration, made the fen a smaller CO2 sink. CH4 fluxes were significantly lowered, close to zero. The impact of natural droughts was similar to, although more modest than, the impact of the experimental water level drawdown. The results are applicable to the short term impacts of the water level drawdown and to climatic conditions in which droughts become more frequent.

Plant secondary compounds and soil microbial processes in carbon and nitrogen cycling in relation to tree species

The aim of this study was to explore soil microbial activities related to C and N cycling and the occurrence and concentrations of two important groups of plant secondary compounds, terpenes and phenolic compounds, under silver birch (Betula pendula Roth), Norway spruce (Picea abies (L.) Karst) and Scots pine (Pinus sylvestris L.) as well as to study the effects of volatile monoterpenes and tannins on soil microbial activities. The study site, located in Kivalo, northern Finland, included ca. 70-year-old adjacent stands dominated by silver birch, Norway spruce and Scots pine. Originally the soil was very probably similar in all three stands. All forest floor layers (litter (L), fermentation layer (F) and humified layer (H)) under birch and spruce showed higher rates of CO2 production, greater net mineralisation of nitrogen and higher amounts of carbon and nitrogen in microbial biomass than did the forest floor layers under pine. Concentrations of mono-, sesqui-, di- and triterpenes were higher under both conifers than under birch, while the concentration of total water-soluble phenolic compounds as well as the concentration of condensed tannins tended to be higher or at least as high under spruce as under birch or pine. In general, differences between tree species in soil microbial activities and in concentrations of secondary compounds were smaller in the H layer than in the upper layers. The rate of CO2 production and the amount of carbon in the microbial biomass correlated highly positively with the concentration of total water-soluble phenolic compounds and positively with the concentration of condensed tannins. Exposure of soil to volatile monoterpenes and tannins extracted and fractionated from spruce and pine needles affected carbon and nitrogen transformations in soil, but the effects were dependent on the compound and its molecular structure. Monoterpenes decreased net mineralisation of nitrogen and probably had a toxic effect on part of the microbial population in soil, while another part of the microbes seemed to be able to use monoterpenes as a carbon source. With tannins, low-molecular-weight compounds (also compounds other than tannins) increased soil CO2 production and nitrogen immobilisation by soil microbes while the higher-molecular-weight condensed tannins had inhibitory effects. In conclusion, plant secondary compounds may have a great potential in regulation of C and N transformations in forest soils, but the real magnitude of their significance in soil processes is impossible to estimate.

Chemical and source characterisation of size-segregated urban air particulate matter

Recent epidemiological studies have shown a consistent association of the mass concentration of urban air thoracic (PM10) and fine (PM2.5) particles with mortality and morbidity among cardiorespiratory patients. However, the chemical characteristics of different particulate size ranges and the biological mechanisms responsible for these adverse health effects are not well known. The principal aims of this thesis were to validate a high volume cascade impactor (HVCI) for the collection of particulate matter for physicochemical and toxicological studies, and to make an in-depth chemical and source characterisation of samples collected during different pollution situations. The particulate samples were collected with the HVCI, virtual impactors and a Berner low pressure impactor in six European cities: Helsinki, Duisburg, Prague, Amsterdam, Barcelona and Athens. The samples were analysed for particle mass, common ions, total and water-soluble elements as well as elemental and organic carbon. Laboratory calibration and field comparisons indicated that the HVCI can provide a unique large capacity, high efficiency sampling of size-segregated aerosol particles. The cutoff sizes of the recommended HVCI configuration were 2.4, 0.9 and 0.2 μm. The HVCI mass concentrations were in a good agreement with the reference methods, but the chemical composition of especially the fine particulate samples showed some differences. This implies that the chemical characterization of the exposure variable in toxicological studies needs to be done from the same HVCI samples as used in cell and animal studies. The data from parallel, low volume reference samplers provide valuable additional information for chemical mass closure and source assessment. The major components of PM2.5 in the virtual impactor samples were carbonaceous compounds, secondary inorganic ions and sea salt, whereas those of coarse particles (PM2.5-10) were soil-derived compounds, carbonaceous compounds, sea salt and nitrate. The major and minor components together accounted for 77-106% and 77-96% of the gravimetrically-measured masses of fine and coarse particles, respectively. Relatively large differences between sampling campaigns were observed in the organic carbon content of the PM2.5 samples as well as the mineral composition of the PM2.5-10 samples. A source assessment based on chemical tracers suggested clear differences in the dominant sources (e.g. traffic, residential heating with solid fuels, metal industry plants, regional or long-range transport) between the sampling campaigns. In summary, the field campaigns exhibited different profiles with regard to particulate sources, size distribution and chemical composition, thus, providing a highly useful setup for toxicological studies on the size-segregated HVCI samples.

Forest inventory-based large-scale forest biomass and carbon budget assessment: new enhanced methods and use of remote sensing for verification

In recent years, concern has arisen over the effects of increasing carbon dioxide (CO2) in the earth's atmosphere due to the burning of fossil fuels. One way to mitigate increase in atmospheric CO2 concentration and climate change is carbon sequestration to forest vegeta-tion through photosynthesis. Comparable regional scale estimates for the carbon balance of forests are therefore needed for scientific and political purposes. The aim of the present dissertation was to improve methods for quantifying and verifying inventory-based carbon pool estimates of the boreal forests in the mineral soils. Ongoing forest inventories provide a data based on statistically sounded sampling for estimating the level of carbon stocks and stock changes, but improved modelling tools and comparison of methods are still needed. In this dissertation, the entire inventory-based large-scale forest carbon stock assessment method was presented together with some separate methods for enhancing and comparing it. The enhancement methods presented here include ways to quantify the biomass of understorey vegetation as well as to estimate the litter production of needles and branches. In addition, the optical remote sensing method illustrated in this dis-sertation can be used to compare with independent data. The forest inventory-based large-scale carbon stock assessment method demonstrated here provided reliable carbon estimates when compared with independent data. Future ac-tivity to improve the accuracy of this method could consist of reducing the uncertainties regarding belowground biomass and litter production as well as the soil compartment. The methods developed will serve the needs for UNFCCC reporting and the reporting under the Kyoto Protocol. This method is principally intended for analysts or planners interested in quantifying carbon over extensive forest areas.

Geological factors affecting on after-use of Finnish cut-over peatlands : with implications on the carbon accumulation

In Finland, peat harvesting sites are utilized down almost to the mineral soil. In this situation the properties of mineral subsoil are likely to have considerable influence on the suitability for the various after-use forms. The aims of this study were to recognize the chemical and physical properties of mineral subsoils possibly limiting the after-use of cut-over peatlands, to define a minimum practice for mineral subsoil studies and to describe the role of different geological areas. The future percentages of the different after-use forms were predicted, which made it possible to predict also carbon accumulation in this future situation. Mineral subsoils of 54 different peat production areas were studied. Their general features and grain size distribution was analysed. Other general items studied were pH, electrical conductivity, organic matter, water soluble nutrients (P, NO3-N, NH4-N, S and Fe) and exchangeable nutrients (Ca, Mg and K). In some cases also other elements were analysed. In an additional case study carbon accumulation effectiveness before the intervention was evaluated on three sites in Oulu area (representing sites typically considered for peat production). Areas with relatively sulphur rich mineral subsoil and pool-forming areas with very fine and compact mineral subsoil together covered approximately 1/5 of all areas. These areas were unsuitable for commercial use. They were recommended for example for mire regeneration. Another approximate 1/5 of the areas included very coarse or very fine sediments. Commercial use of these areas would demand special techniques - like using the remaining peat layer for compensating properties missing from the mineral subsoil. One after-use form was seldom suitable for one whole released peat production area. Three typical distribution patterns (models) of different mineral subsoils within individual peatlands were found. 57 % of studied cut-over peatlands were well suited for forestry. In a conservative calculation 26% of the areas were clearly suitable for agriculture, horticulture or energy crop production. If till without large boulders was included, the percentage of areas suitable to field crop production would be 42 %. 9-14 % of all areas were well suitable for mire regeneration or bird sanctuaries, but all areas were considered possible for mire regeneration with correct techniques. Also another 11 % was recommended for mire regeneration to avoid disturbing the mineral subsoil, so total 20-25 % of the areas would be used for rewetting. High sulphur concentrations and acidity were typical to the areas below the highest shoreline of the ancient Litorina sea and Lake Ladoga Bothnian Bay zone. Also differences related to nutrition were detected. In coarse sediments natural nutrient concentration was clearly higher in Lake Ladoga Bothnian Bay zone and in the areas of Svecokarelian schists and gneisses, than in Granitoid area of central Finland and in Archaean gneiss areas. Based on this study the recommended minimum analysis for after-use planning was for pH, sulphur content and fine material (<0.06 mm) percentage. Nutrition capacity could be analysed using the natural concentrations of calcium, magnesium and potassium. Carbon accumulation scenarios were developed based on the land-use predictions. These scenarios were calculated for areas in peat production and the areas released from peat production (59300 ha + 15 671 ha). Carbon accumulation of the scenarios varied between 0.074 and 0.152 million t C a-1. In the three peatlands considered for peat production the long term carbon accumulation rates varied between 13 and 24 g C m-2 a-1. The natural annual carbon accumulation had been decreasing towards the time of possible intervention.

Below-ground processes in meadow soil under elevated ozone and carbon dioxide : - Greenhouse gas fluxes, N cycling and microbial communities

This thesis focuses on how elevated CO2 and/or O3 affect the below-ground processes in semi-natural vegetation, with an emphasis on greenhouse gases, N cycling and microbial communities. Meadow mesocosms mimicking lowland hay meadows in Jokioinen, SW Finland, were enclosed in open-top chambers and exposed to ambient and elevated levels of O3 (40-50 ppb) and/or CO2 (+100 ppm) for three consecutive growing season, while chamberless plots were used as chamber controls. Chemical and microbiological analyses as well as laboratory incubations of the mesocosm soils under different treatments were used to study the effects of O3 and/or CO2. Artificially constructed mesocosms were also compared with natural meadows with regards to GHG fluxes and soil characteristics. In addition to research conducted at the ecosystem level (i.e. the mesocosm study), soil microbial communities were also examined in a pot experiment with monocultures of individual species. By comparing mesocosms with similar natural plant assemblage, it was possible to demonstrate that artificial mesocosms simulated natural habitats, even though some differences were found in the CH4 oxidation rate, soil mineral N, and total C and N concentrations in the soil. After three growing seasons of fumigations, the fluxes of N2O, CH4, and CO2 were decreased in the NF+O3 treatment, and the soil NH4+-N and mineral N concentrations were lower in the NF+O3 treatment than in the NF control treatment. The mesocosm soil microbial communities were affected negatively by the NF+O3 treatment, as the total, bacterial, actinobacterial, and fungal PLFA biomasses as well as the fungal:bacterial biomass ratio decreased under elevated O3. In the pot survey, O3 decreased the total, bacterial, actinobacterial, and mycorrhizal PLFA biomasses in the bulk soil and affected the microbial community structure in the rhizosphere of L. pratensis, whereas the bulk soil and rhizosphere of the other monoculture, A. capillaris, remained unaffected by O3. Elevated CO2 caused only minor and insignificant changes in the GHG fluxes, N cycling, and the microbial community structure. In the present study, the below-ground processes were modified after three years of moderate O3 enhancement. A tentative conclusion is that a decrease in N availability may have feedback effects on plant growth and competition and affect the N cycling of the whole meadow ecosystem. Ecosystem level changes occur slowly, and multiplication of the responses might be expected in the long run.

Aerobic carbon-cycle related microbial communities in boreal peatlands: responses to water-level drawdown

Boreal peatlands represent a considerable portion of the global carbon (C) pool. Water-level drawdown (WLD) causes peatland drying and induces a vegetation change, which affects the decomposition of soil organic matter and the release of greenhouse gases (CO2 and CH4). The objective of this thesis was to study the microbial communities related to the C cycle and their response to WLD in two boreal peatlands. Both sampling depth and site type had a strong impact on all microbial communities. In general, bacteria dominated the deeper layers of the nutrient-rich fen and the wettest surfaces of the nutrient-poor bog sites, whereas fungi seemed more abundant in the drier surfaces of the bog. WLD clearly affected the microbial communities but the effect was dependent on site type. The fungal and methane-oxidizing bacteria (MOB) community composition changed at all sites but the actinobacterial community response was apparent only in the fen after WLD. Microbial communities became more similar among sites after long-term WLD. Litter quality had a large impact on community composition, whereas the effects of site type and WLD were relatively minor. The decomposition rate of fresh organic matter was influenced slightly by actinobacteria, but not at all by fungi. Field respiration measurements in the northern fen indicated that WLD accelerates the decomposition of soil organic matter. In addition, a correlation between activity and certain fungal sequences indicated that community composition affects the decomposition of older organic matter in deeper peat layers. WLD had a negative impact on CH4 oxidation, especially in the oligotrophic fen. Fungal sequences were matched to taxa capable of utilizing a broad range of substrates. Most of the actinobacterial sequences could not be matched to characterized taxa in reference databases. This thesis represents the first investigation of microbial communities and their response to WLD among a variety of boreal peatland habitats. The results indicate that microbial community responses to WLD are complex but dependent on peatland type, litter quality, depth, and variable among microbes.

The role of lakes in carbon cycling in boreal catchments

Lakes are an important component of ecosystem carbon cycle through both organic carbon sequestration and carbon dioxide and methane emissions, although they cover only a small fraction of the Earth's surface area. Lake sediments are considered to be one of rather perma-nent sinks of carbon in boreal regions and furthermore, freshwater ecosystems process large amounts of carbon originating from terrestrial sources. These carbon fluxes are highly uncer-tain especially in the changing climate. -- The present study provides a large-scale view on carbon sources and fluxes in boreal lakes situated in different landscapes. We present carbon concentrations in water, pools in lake se-diments, and carbon gas (CO2 and CH4) fluxes from lakes. The study is based on spatially extensive and randomly selected Nordic Lake Survey (NLS) database with 874 lakes. The large database allows the identification of the various factors (lake size, climate, and catchment land use) determining lake water carbon concentrations, pools and gas fluxes in different types of lakes along a latitudinal gradient from 60oN to 69oN. Lakes in different landscapes vary in their carbon quantity and quality. Carbon (C) content (total organic and inorganic carbon) in lakes is highest in agriculture and peatland dominated areas. In peatland rich areas organic carbon dominated in lakes but in agricultural areas both organic and inorganic C concentrations were high. Total inorganic carbon in the lake water was strongly dependent on the bedrock and soil quality in the catchment, especially in areas where human influence in the catchment is low. In inhabited areas both agriculture and habitation in the catchment increase lake TIC concentrations, since in the disturbed soils both weathering and leaching are presumably more efficient than in pristine areas. TOC concentrations in lakes were related to either catchment sources, mainly peatlands, or to retention in the upper watercourses. Retention as a regulator of the TOC concentrations dominated in southern Finland, whereas the peatland sources were important in northern Finland. The homogeneous land use in the north and the restricted catchment sources of TOC contribute to the close relationship between peatlands and the TOC concentrations in the northern lakes. In southern Finland the more favorable climate for degradation and the multiple sources of TOC in the mixed land use highlight the importance of retention. Carbon processing was intensive in the small lakes. Both CO2 emission and the Holocene C pool in sediments per square meter of the lake area were highest in the smallest lakes. How-ever, because the total area of the small lakes on the areal level is limited, the large lakes are important units in C processing in the landscape. Both CO2 and CH4 concentrations and emissions were high in eutrophic lakes. High availability of nutrients and the fresh organic matter enhance degradation in these lakes. Eutrophic lakes are often small and shallow, enabling high contact between the water column and the sediment. At the landscape level, the lakes in agricultural areas are often eutrophic due to fertile soils and fertilization of the catchments, and therefore they also showed the highest CO2 and CH4 concentrations. Export from the catchments and in-lake degradation were suggested to be equally important sources of CO2 and CH4 in fall when the lake water column was intensively mixed and the transport of sub-stances from the catchment was high due to the rainy season. In the stagnant periods, especially in the winter, in-lake degradation as a gas source was highlighted due to minimal mixing and limited transport of C from the catchment. The strong relationship between the annual CO2 level of lakes and the annual precipitation suggests that climate change can have a major impact on C cycling in the catchments. Increase in precipitation enhances DOC export from the catchments and leads to increasing greenhouse gas emissions from lakes. The total annual CO2 emission from Finnish lakes was estimated to be 1400 Gg C a-1. The total lake sediment C pool in Finland was estimated to be 0.62 Pg, giving an annual sink in Finnish lakes of 65 Gg C a-1.

Country-scale carbon accounting of the vegetation and mineral soils of Finland

The increase in global temperature has been attributed to increased atmospheric concentrations of greenhouse gases (GHG), mainly that of CO2. The threat of severe and complex socio-economic and ecological implications of climate change have initiated an international process that aims to reduce emissions, to increase C sinks, and to protect existing C reservoirs. The famous Kyoto protocol is an offspring of this process. The Kyoto protocol and its accords state that signatory countries need to monitor their forest C pools, and to follow the guidelines set by the IPCC in the preparation, reporting and quality assessment of the C pool change estimates. The aims of this thesis were i) to estimate the changes in carbon stocks vegetation and soil in the forests in Finnish forests from 1922 to 2004, ii) to evaluate the applied methodology by using empirical data, iii) to assess the reliability of the estimates by means of uncertainty analysis, iv) to assess the effect of forest C sinks on the reliability of the entire national GHG inventory, and finally, v) to present an application of model-based stratification to a large-scale sampling design of soil C stock changes. The applied methodology builds on the forest inventory measured data (or modelled stand data), and uses statistical modelling to predict biomasses and litter productions, as well as a dynamic soil C model to predict the decomposition of litter. The mean vegetation C sink of Finnish forests from 1922 to 2004 was 3.3 Tg C a-1, and in soil was 0.7 Tg C a-1. Soil is slowly accumulating C as a consequence of increased growing stock and unsaturated soil C stocks in relation to current detritus input to soil that is higher than in the beginning of the period. Annual estimates of vegetation and soil C stock changes fluctuated considerably during the period, were frequently opposite (e.g. vegetation was a sink but soil was a source). The inclusion of vegetation sinks into the national GHG inventory of 2003 increased its uncertainty from between -4% and 9% to ± 19% (95% CI), and further inclusion of upland mineral soils increased it to ± 24%. The uncertainties of annual sinks can be reduced most efficiently by concentrating on the quality of the model input data. Despite the decreased precision of the national GHG inventory, the inclusion of uncertain sinks improves its accuracy due to the larger sectoral coverage of the inventory. If the national soil sink estimates were prepared by repeated soil sampling of model-stratified sample plots, the uncertainties would be accounted for in the stratum formation and sample allocation. Otherwise, the increases of sampling efficiency by stratification remain smaller. The highly variable and frequently opposite annual changes in ecosystem C pools imply the importance of full ecosystem C accounting. If forest C sink estimates will be used in practice average sink estimates seem a more reasonable basis than the annual estimates. This is due to the fact that annual forest sinks vary considerably and annual estimates are uncertain, and they have severe consequences for the reliability of the total national GHG balance. The estimation of average sinks should still be based on annual or even more frequent data due to the non-linear decomposition process that is influenced by the annual climate. The methodology used in this study to predict forest C sinks can be transferred to other countries with some modifications. The ultimate verification of sink estimates should be based on comparison to empirical data, in which case the model-based stratification presented in this study can serve to improve the efficiency of the sampling design.