ANALYZING CARBON STOCKS IN A MEDITERRANEAN FOREST ENTERPRISE: A CASE STUDY FROM KIZILDAG, TURKEY

Carbon dioxide is one of the most important greenhouse gases responsible for global climate change. Terrestrial ecosystems are widely recognized as useful mechanisms, as they store large quantities of carbon in their bodies. In this context, the present study aims at analyzing carbon stocks in forest (both productive and degraded), grassland, and agricultural ecosystems of Kizildag Planning Unit. To this end, forest inventory data, published statistics, and GIS tools were used for modeling carbon densities of these three ecosystems in a spatially explicit manner. Moreover, carbon contents were estimated for different pools including above- and below-ground biomass, litter layer, deadwood, and soil components. Results showed that a total of nearly 3.5 million tons of carbon stored in the entire study area. When standardized to per unit area, productive and degraded forests stocked the most significant amounts followed by grassland and cropland land use/land cover (LULC) classes. The average density was 120.8 tons ha-1 for productive forests, while it was 32.3 tons ha-1


INTRODUCTION
The greenhouse effect has become harmful for the global climate as a consequence of the increased anthropogenic activities in recent centuries. It is considered that the essential part of this effect stemming from carbon dioxide (CO 2 ) gases. Thus, international scientific research has been focused on the dynamics of CO 2 in the last decades. They indicate that CO 2 concentration in the atmosphere has increased from 280 ppm to 412 ppm after the industrial revolution (NOAA, 2013;Sharma et al., 2011;URL-1). In this context, many global mechanisms such as REDD + , Kyoto and LULUCF have been introduced for widely investigate the carbon deposits of terrestrial ecosystems as well as their trading potential to the countries (Aryal et al., 2013).
Among terrestrial ecosystems; forests, grasslands, and agricultural areas are the primary deposits of carbon. These ecosystems play an essential role in the global CO 2 cycle since they can sequester carbon and store it in their bodies effectively. Biomass expansion factors are mostly used for determining carbon amounts in these ecosystems. Forests, for instance, have an enormous storage capacity of aboveground biomass compared to other terrestrial ecosystems. They store almost 80 % of global aboveground carbon on their own (Waring and Running, 1998;Asan, 1999). Accordingly, a total of 650 billion tons of carbon is stored within forests all over the world (FAO, 2010). In a forest ecosystem, carbon is stored in five different pools including (i) aboveground biomass, (ii) belowground biomass, (iii) dead wood, (iv) litter layer and (v) soil. Within these pools, approximately 75 % of the net primary productivity is realized in forest ecosystems worldwide (Melillo et al., 1993).
Grassland ecosystems, on the other hand, are another vital pool, as they contain about 12% of terrestrial carbon stocks globally (Adams et al., 1990;Ojima et al., 1993). Unlike forests, they store carbon mostly in their belowground components. That is why grasslands are dominated by herbaceous (non-woody) vegetation and they are relatively short-lived due to their perennial lifespan as well as grazing activities. However, their extensive root systems may store almost 80% of the total carbon in the soil (Glatzle, 2012). Many studies have been carried out to estimate gross grassland carbon using different methods at various scales (Foody et al., 2003;Kumar et al., 2015). Biomass harvesting is common in these studies (Dinc et al., 2018) although it has significant limitations, as discussed by Jobbagy et al. (2002). Another limitation relates to the high spatial variation of aboveground biomass on grasslands. To overcome such limitations, researchers attribute the aboveground biomass with the data of precipitation, temperature, and soil (Sala et al., 1988;Epstein et al., 1997;Dinc et al., 2018). Along with such datasets, remote sensing products are useful for effectively mapping the complex spatial patterns in large areas (Lauenroth, 1979;Sala et al., 1988;Burke et al., 1997;Jobbagy et al., 2002;Ni, 2004;Fang et al., 2005;Dinc et al., 2018).
Regarding agricultural areas, carbon levels stored both in above-and below-ground biomass strictly depends on the crop characteristics. Therefore, they may markedly differ based on both plant species and its developmental stage (Lichaikul, 2004;Jenkins et al., 2003). More or less, crops are capable of carbon sequestration in the form of accumulated biomass in their stems and soils (IPCC, 2003;IPCC, 2006;FAO, 2012). Thus, the agriculture sector has also the potential to mitigate global climate change. The critical point here is that decisionmakers should always consider not to have any adverse effect on food production while they are planning to carbon management activities. Otherwise, global food security may be at risk in the future.
The countries that have signed the Kyoto Protocol are obliged to declare their forest resources and their contribution to global carbon cycle. In this context, Turkey has periodically reported both its carbon storage and oxygen production amounts based on landscape-level forest management plans (Seki et al., 2017). However, they are being calculated only for the forested lands and identified by total amounts for forest enterprises. Therefore, carbon storage capacities for neither nonforest lands (e.g., openings without tree) nor other ecosystem types (e.g., grassland, agriculture, etc.) can be determined. Moreover, spatially explicit presentation of these capacities are lacking in the management plans. Yet, it is of crucial importance to know the ratio of storage capacities per unit area belonging to different land use/ land cover (LULC) classes.
The objective of this study is to analyze the carbon storage capacities of forest, grassland, and agricultural ecosystems in the Kizildag Forest Planning Unit (FPU) located at the southern Turkey. To this end, we calculated and mapped the carbon stocks within these ecosystems based on field data and published statistics using spatial analysis methods in a Geographical Information Systems (GIS) software. Furthermore, total carbon stocks were subdivided into its above-and below-ground components such as vegetation, deadwood, and soil. Obtained results are expected to be useful to natural resource managers in their decision-making processes. Another importance of this study lies in the study area, as it is in the most vulnerable sub-region (i.e., the western Mediterranean) to future climate change in Turkey Semerci et al., 2018;Türkeş et al., 2018).

CERNE
DINC et al. Kizildag FPU stretches from 4125732-4153251 N to 353172-4138184 E (UTM Zone36N, WGS84 Datum) and is located between Antalya and Konya provinces in the Mediterranean region of Turkey ( Figure 1). Its mountainous lands cover 53290 ha with an average slope rate of 31 % ( Figure 2). LULC is very heterogeneous with highly diverse tree species (Table 1). Moreover, land covers within forest and grassland are mostly degraded. Human-induced degradation is attributed to illegal logging, improper land use and overgrazing (OGM, 2016). The region shows a transition climate type between Mediterranean and continental. Accordingly, summer seasons are hot and dry, while the winters are cold and rainy. The monthly average temperature ranged from 1.8 0 C in January to 23.4 0 C in July with a total annual average precipitation of 564 mm (SMS, 2018). Based on Erinc's aridity index (Erinc, 1984), the study area is in the arid region due to rain mostly drop in December and January ( Figure 3). Summer drought affects a significant part of the planning unit during the vegetation period. It is an essential factor influencing plant growth and its yield when coupled with the sunny aspects ( Figure 2). The dominant soil type in the study area is brown forest soil (Ozkan et al., 2007).

Field sampling and modeling
Forest inventory data were used to determine the biomass amounts in the forest ecosystem. For this purpose, 264 sampling plots were distributed all over the forested lands with a systematical sampling scheme. Circular plot sizes were 400 m 2 , 600 m 2 or 800 m 2 based on the canopy cover classes. Ground measurements were conducted by  professional forest engineers in the summer of 2015. Then, stand parameters including species mixture, diameter at breast height (dbh) and stem quality were recorded into inventory sheets. No measurement was performed in other land use classes except for forests. Accordingly, information on agricultural areas was retrieved from annual statistics published by TUIK (2016). Grassland data, on the CERNE DINC et al.

ANALYZING CARBON STOCKS IN A MEDITERRANEAN FOREST ENTERPRISE: A CASE STUDY FROM KIZILDAG, TURKEY
(BEF). WD and BEF coefficients developed by Turkish researchers for native tree species were given in Table  2. Then, carbon amounts were calculated by multiplying aboveground biomass (AGB) with a coefficient of 0.51 for coniferous and 0.48 for deciduous forests (IPCC, 2006). Finally, these values were aggregated to plot level and converted to per unit area (i.e., 1 ha).

Calculation of carbon stocks
The amount of carbon stored in Kizildag FPU was determined using the inventory-based carbon stock change method (IPCC, 2003;2006;Serengil, 2018). Accordingly, biomass amounts were first calculated for 5 carbon pools (i.e., aboveground, belowground, litter, deadwood, and soil) in the ecosystem. Then, they were converted to carbon content as explained in the next subsections.

Carbon content in aboveground biomass
In forested lands, aboveground biomass for each tree in a sample plot was first calculated by multiplying wood density (WD) and biomass expansion factor In grasslands, total biomass was first calculated by Equation1 (Dinc et al., 2018) based on the precipitation data obtained from Şarkikaraağaç (MGM, 2019). Then, Lales et al. (2001)'s coefficient (i.e., 0.5) was used for conversion to carbon. Finally, 27 % of this amount was considered as aboveground, while the rest of it stands for belowground carbon content, as suggested by Dinc et al. (2018) for grasslands of Turkey.
In agricultural lands, both cultivated areas and crop types on them were obtained from the Turkish Statistical Institute's reports (TUIK, 2016). Based on these statistics, aboveground biomass was calculated by the related coefficients in Table 3 for per unit area. Then, carbon amounts were calculated by multiplying AGB with 0.45, which was a coefficient suggested by Surur et al. (2013).

Carbon content in belowground biomass
In forested lands, belowground biomass (BGB) and its carbon equivalent were calculated by the help of IPCC guidelines (i.e., Agriculture, forestry, and other land use guideline) (IPCC, 2003;2006). Accordingly, root to shoot ratios and carbon factors for temperate forests were utilized depending on forests' structural characteristics such as species and growing stock. Relevant coefficients used in this guideline were given in Table 4.   Forest inventory data were used for calculating deadwood biomass of each tree at the plot level. Then, the carbon content was calculated by multiplying the deadwood biomass with 0.47 as reported by the Turkish Forest Management Guideline (FRA, 2010;OGM, 2017). No calculations were made for non-forested lands.

Carbon content in soil
In forested lands, soil carbon was determined using the average values given in Table 6. These values were developed by tree species for both degraded and productive forests based on extensive field studies conducted in Turkey between 1974 and 2008. In grasslands, on the other hand, the values reported by Dinc (2017) were used. For agricultural lands, finally, the values suggested by Parlak et al. (2008) were directly employed.
As for grasslands, total biomass was first calculated by Equation 1 given in the previous section. Then, 73 % of this amount was considered as BGB as suggested by Dinc et al. (2018). Finally, the carbon content was multiplied with 0.5 which was a conversion coefficient reported by Lales et al. (2001).
In agricultural lands, BGB was calculated by multiplying AGB with 0.12, which was a conversion coefficient (Surur et al., 2013). Subsequently, the carbon content was calculated by multiplying BGB with 0.45, which was a conversion coefficient for agricultural lands (Surur et al., 2013).

Carbon content in litter
In forested lands, the weights of the litter layer were first determined by each stand depending on the dominant tree species for per unit area. Organic matter in these weights was then calculated. Finally, carbon content in the organic matter was found by the help of the average values given in Table 5 , and Equation4 were used for forest, grassland and agriculture ecosystems, respectively, where TC is the total carbon stock, AGC is the aboveground carbon stock, BGC is the belowground carbon stock, DWC is the carbon stock in deadwood, LC is the carbon stock in litter layer, SC is the carbon stock in the soil.
Although there are some marginal values which are lower than 37.9 tons ha -1 (Dixon et al., 1994;Evrendilek, 2014), our finding is slightly smaller than the essential literature. The main reason for that may be the harsh environmental conditions within the study area. As seen in Figure 2, almost half of Kizildag FPU has a slope rate of more than 40 %. This causes the soil to remain shallow. As a result, the soil cannot store enough water which is essential for plant growth. Indeed, large quantities of tree mortalities are widespread in the Mediterranean region. In a recent study conducted in this region (Semerci et al., 2018), almost 30 % of the growing stock was found to be dieback, especially in Taurus Fir stands. They attributed this phenomenon to severe drought stress in the last decades. Taurus Fir was dominant tree species in Kizildag FPU, as well.
Specific emphasis should also be given to other carbon pools including soil, litter, and deadwood in forest ecosystems. The soil carbon, in particular, is generally overlooked in total carbon calculations due to lack of relevant data. However, the amount of carbon stored in forest soils is much higher than those of biomass. It was determined by the dynamic interaction between organic matter entry and decomposition of this matter (Tolunay and Çömez, 2007). Fallen leaves from trees are the first entry for abundant organic matter in the soil (Pausas, 1997). However, decomposition of fine roots also contributes to this amount (Berg et al., 2003). Aside from the soil, the litter layer is another carbon pool in a high forest. A mean of 8.7 tons . ha -1 carbon was stored within the litter layer in Kizildag FPU. In another study by Mısır et al. (2013), it ranged from 2 tons ha -1 to 7 tons ha -1 in Giresun, Turkey. The differences can be attributed to both stand types and tree species. The researchers studied in a pure Beech forest which was intensively managed by relatively short rotations. The stand types in Kizildag, on the other hand, is dominated by mixed-or Fir-species managed primarily for protective functions.
Unlike forests, the lowest carbon density was observed in croplands. It was between 30.2 and 45.7 tons ha -1 with a mean of 32.3 tons . ha -1 (Table 7). More than [3] [4]

Carbon densities by LULC classes
The carbon density in productive forest class ranged between 73 and 246.5 tons . ha -1 with a mean of 120.8 tons . ha -1 (Table 7). It had the highest density among all LULC classes. Kauppi et al. (1992), Gifford et al. (1992) and Lee et al. (2015) also found mean carbon densities of 122 tons . ha -1 in Europe, 128 tons ha -1 in China, and 130 tons . ha -1 in Korea's forests, respectively. They related the higher carbon densities to the areas dominated by mixed forest ecosystems. In another study conducted near the Kizildag FPU, higher carbon storages were estimated in mixed forest stands, as well (Mirici et al., 2015). Similarly, the higher carbon stocks were observed in Fir-Cedar mixed stands in this study, as discussed in the next section. Thus, it can be said that our results -in general -are in good agreement with the relevant literature.
The biomass carbon density (i.e., the summed of above-and below-ground carbon values) was 37.9 tons ha -1 for productive forests in the present study. Tolunay (2011), in another study, found a mean value of 41.7 tons . ha -1 for Turkey's forests while UN-ECE/FAO (2006) reported it as 43.9 tons . ha -1 for European forests. DINC et al. 90 % of this amount (i.e., 29.7 tons . ha -1 ) stored only in the soil. Similarly, Shah et al. (2015) reported that 32.71 tons . ha -1 carbon stored in the soil of Pakistan's croplands. Grassland LULC class, on the other hand, showed a carbon density of 40.1 tons . ha -1 in the present study. As seen in Table 7, there was no range within grasslands due to the calculation method we used. However, it can be said that they have nominal storage capacity in their different carbon pools. There are many grassland studies conducted in Turkey reporting similar results to ours (Dinc, 2018;Kendir, 1999;Yavuz et al., 2013;Bilgin, 2010). Laswai (2011), on the other hand, found a carbon density value of 65.8 ton . ha -1 for Tanzania's grasslands. These kind of higher values should be considered ordinary since water is an essential part of plant growth in grasslands and Kizildag suffers from summer drought as can be seen in Figure 3. Furthermore, the rate of AGB is expected to decrease in the summer season, while it increases for BGB. Therefore, both temporal and seasonal changes should be taken into account in such studies, as already reported by Kahmen et al. (2005).

Density maps for different carbon pools
The finding of the present study revealed that both aboveground and belowground carbon densities were the highest in productive forests. The spatial distribution of these two essential carbon pools can be seen in Figure 4a and Figure 4b. The similarity between two maps was evident in the figures. However, carbon densities of them were totally different as seen in their legends. The green hotspots showing high storage capacities located at both southern and northern parts of the Kizildag FPU. According to forest stand type map, these areas generally consisted of mixed (Cedar-Fir) forest stands that were at higher age-classes from 80 to 140 years old. Thus, they had a lot of large-girth trees whose dbh classes were d (36-52 cm) and e (>52 cm) in the Turkish forest management system (OGM, 2017).
The litter carbon densities were mapped for degraded and productive forest LULC classes ( Figure  4c). Accordingly, litter carbon hotspots were generally aggregated in the northwestern parts of Kizildag. Pure fir or fir-dominated mixed stands were common in these areas. Another distinctive feature of these stands was that they were almost fully-covered at all. The lowcapacity regions, which were shown in red color on the map, located at both the center and the eastern parts of the study area. These areas generally consisted of degraded (i.e., canopy cover < 10 %) forest stands.
The deadwood carbon density was mapped only for productive forests. Regarding deadwood carbon, Kizildag showed a relatively homogenous pattern as seen in Figure 4d. Red colors were dominant in the available areas indicating very low-capacity of carbon storage in deadwood. Few hotspots existed (e.g., in orange color) in the regions with mature forests. These forest stands overlapped with the hotspots of above-and below-ground carbon density maps showed in Figure  4a and Figure 4b. It may be attributed to the extensive management practices on these sites. Indeed, these sites are mainly allocated to soil protection function in the forest management plan (OGM, 2016) due to the high slope rates higher than 60 %. Therefore, no or minimal intervention is performed to these forest sites. Higher deadwood amounts can be explained in this way. The deadwood carbon densities in Figure 4d ranged between 0 and 41 tons ha -1 with an average of 1.2 tons ha -1 . It was about 4 % of the average aboveground carbon density as seen in Table 7. Karahalil et al. (2017) conducted a comprehensive study on deadwood near Kizildag. They obtained similar findings to ours. According to this study, the mean deadwood volume constituted only 2.4 % of the total standing volume in Koprulu Canyon National Park, Antalya. The slight difference may be attributed to different tree species in two study areas. Unlike us, they studied in a Calabrian pine-dominated forest.
As for soil carbon density, a non-homogenous spatial pattern was observed in the study area. Different patches of both high-and low-capacity sites existed together. However, the darkest green hotspots dominated the northwestern parts of the Kizildag FPU (Figure 4e). These sites spatially correlated with litter carbon density map (Figure 4c). There were also many hotspots in the eastern parts. In general, they consisted of mixed stands of Juniper and Fir with various developmental stages. This finding is in agreement with Özkan (2003)'s study. He has investigated different site factors in forest ecosystem near the Beysehir Lake which is very close to our study area. He found that soil carbon amount in Fir-Cedar mixed forest was the highest among other forest types. Interestingly, many degraded stands showed a high-capacity of carbon storage in the soil in our case.
Total carbon density at the landscape level In order to map the total carbon stock, all maps in Figure 4 were overlapped and aggregated into one GIS layer. It was seen that red and yellow colors were dominant in the final map ( Figure 5). This was indicating very low-and moderate-capacities of carbon storages were typical for the entire FPU. In general, the CERNE DINC et al.

ANALYZING CARBON STOCKS IN A MEDITERRANEAN FOREST ENTERPRISE: A CASE STUDY FROM KIZILDAG, TURKEY
areas presenting moderate-capacities consisted of the productive forest with different characteristics in terms of tree species and developmental stages. Grasslands and croplands, on the other hand, located at the very low-capacity sites colored in red. Degraded forests dominated the orange sections showing low-capacities. There were some exceptions, as expected. Some of the green hotspots, for example, consisted of old-growth forest sites which were designated for nature protection in the forest management plan (OGM, 2016). The mixed forest stands dominated other hotspots. These hotspots had large-girth Fir, Cedar, and Black Pine trees that were common in the region. Mirici et al. (2015), Kauppi et al. (1992), Gifford et al. (1992) and Lee et al. (2015)'s studies were also consistent with the present study. They found the highest carbon values in mixed forest ecosystems in their study areas. Finally, the grey areas on the map had no carbon stocks. They were irrelevant LULC classes; settlement, mine sites, rocky areas etc. Among these, the power line located in the east-west direction can be seen in the upper part of the study area in Figure 5.
The total carbon stock was calculated as 3404828 tons for the entire study area (Table 8). Forest was the best LULC class in Kizildag FPU in terms of maximum carbon storage capacity. Nearly 90 % of the total carbon stored only in forest ecosystems both in their degraded and productive types. Grassland, on the other hand, was the worst LULC class as it could save only 2 % of the total carbon amount. The remaining carbon (ca. 8 %) was retained by croplands covering 16 % of the whole study area (Table 8). As expected, vast differences were observed amongst the different LULC classes. Nevertheless, these differences were affected by the area coverage of LULC classes across the planning unit, as well.
There were also significant differences among carbon contents in the different components (C pools) of LULC classes at the landscape level. In productive forests, for example, approximately 60 % of the total carbon  (Table 8). Taking all LULC classes together, the distributions into carbon components were app. 73.0 % for soil, 15.4 % for aboveground, 6.5 % for litter, 4.5 % for belowground, and 0.6 % for deadwood. As seen in Table 8, almost ¾ of the total amount stored in the soil in Kizildag FPU.

CONCLUSION
In the present study, the amounts and spatial distribution of carbon stocks for forest, grassland and agriculture ecosystems were analyzed within the boundaries of the Kizildag FPU. The total stock was found to be as almost 3.5 million tons at the landscape level. Carbon hotspots were generally located at the southwestern parts of the study area which were adjacent to Kimyos Plain. These parts were the old-growth and mixed forest sites designated for protective functions. When standardized and allocated into ecosystem types, average stocks were 120.8 tons ha -1 , 51.3 tons ha -1 , 40.1 tons ha -1 , and 32.3 tons ha -1 for productive forest, degraded forest, grassland, and cropland, respectively. Regarding carbon pools, the soil was seen as an essential component for stocking. It was obviously the best pool followed by aboveground biomass in terms of carbon storage in its body.
It is concluded that forest ecosystems with their productive and degraded forms are crucial for carbon management. They have a considerable potential for stocking significant amounts of carbon in their different components. By utilizing the density maps presented in this study, forest managers could improve the capacity of degraded ecosystems covered more than one-third of the Kizildag's lands by proper forestry activities such as rehabilitation, afforestation, and protection. Thus, alleviating the harmful effects of climate change might be a viable option at least at the landscape level. In time, we believe that more effort will be spent and successful results will be achieved in combatting climate change both at national and global levels.   Series A, v.29, n.1, p.14-53, (in Turkish), 1979a.