Analyzing the Joint Effect of Forest Management and Wildfires on Living Biomass and Carbon Stocks in Spanish Forests

Research Highlights: This is the first study that has considered forest management and wildfires in the balance of living biomass and carbon stored in Mediterranean forests. Background and Objectives: The Kyoto Protocol and Paris Agreement request countries to estimate and report carbon emissions and removals from the forest in a transparent and reliable way. The aim of this study is to forecast the carbon stored in the living biomass of Spanish forests for the period 2000–2050 under two forest management alternatives and three forest wildfires scenarios. Materials and Methods: To produce these estimates, we rely on data from the Spanish National Forest Inventory (SNFI) and we use the European Forestry Dynamics Model (EFDM). SNFI plots were classified according to five static (forest type, known land-use restrictions, ownership, stand structure and bioclimatic region) and two dynamic factors (quadratic mean diameter and total volume). The results were validated using data from the latest SNFI cycle (20-year simulation). Results: The increase in wildfire occurrence will lead to a decrease in biomass/carbon between 2000 and 2050 of up to 22.7% in the medium–low greenhouse gas emissions scenario (B2 scenario) and of up to 32.8% in the medium–high greenhouse gas emissions scenario (A2 scenario). Schoolbook allocation management could buffer up to 3% of wildfire carbon loss. The most stable forest type under both wildfire scenarios are Dehesas. As regards bioregions, the Macaronesian area is the most affected and the Alpine region, the least affected. Our validation test revealed a total volume underestimation of 2.2% in 20 years. Conclusions: Forest wildfire scenarios provide more realistic simulations in Mediterranean forests. The results show the potential benefit of forest management, with slightly better results in schoolbook forest management compared to business-as-usual forest management. The EFDM harmonized approach simulates the capacity of forests to store carbon under different scenarios at national scale in Spain, providing important information for optimal decision-making on forest-related policies.


Introduction
Forests store large amounts of carbon [1] and are a critical component of the global carbon cycle as they store over 80% of global terrestrial above-ground carbon [2]. According to the State of European Forests [3], in Europe, an average of 35.6% of carbon is stored in living biomass, comprised Inventory (SNFI). We hypothesized that there will be a combined effect of, on the one hand, the increase in living biomass/carbon stored in the Spanish forests following the abandonment of traditional land uses and forest expansion, and on the other, the living biomass/carbon losses associated with both a more intensive management and an increase in forest wildfires based on future climatic scenarios. These forecasts provide key evidence for policy makers to disentangle the potential trade-offs between the effects of forest management and global change on Mediterranean forests. This is the first time that natural disturbances such as forest wildfires have been considered in the balance of living biomass and carbon stored in Mediterranean forests.

SNFI Data and Auxiliary Information
The SNFI is the primary source of forest data for national and large-area assessments because of the high number of monitored plots in the forests [34]. Since the Second SNFI, denoted SNFI2 (1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996), a continuous inventory has been conducted with permanent plots established in forested areas with an intensity of approximately of one plot per km 2 . Each plot is re-measured approximately every 10 years. A total of 81,024 sample plots measured in the Third Spanish National Forest Inventory (SNFI3, 1997(SNFI3, -2007, were used as the input data for the EFDM (Figure 1). For validation purposes, the 13,667 plots measured to date in the Fourth Spanish National Forest Inventory (SNFI4, 2008-to date) were used (Table 1).
future projections of the above-ground carbon and that associated with fellings in 23 European countries.
The objective of this study is to use the EFDM to forecast the living biomass and carbon of Spanish forests for the period  under three forest wildfires scenarios [33] and two different forest management regimes (business-as-usual allocation and schoolbook allocation). To evaluate the results, a validation process was carried out using available data from the latest Spanish National Forest Inventory (SNFI). We hypothesized that there will be a combined effect of, on the one hand, the increase in living biomass/carbon stored in the Spanish forests following the abandonment of traditional land uses and forest expansion, and on the other, the living biomass/carbon losses associated with both a more intensive management and an increase in forest wildfires based on future climatic scenarios. These forecasts provide key evidence for policy makers to disentangle the potential trade-offs between the effects of forest management and global change on Mediterranean forests. This is the first time that natural disturbances such as forest wildfires have been considered in the balance of living biomass and carbon stored in Mediterranean forests.

SNFI Data and Auxiliary Information
The SNFI is the primary source of forest data for national and large-area assessments because of the high number of monitored plots in the forests [34]. Since the Second SNFI, denoted SNFI2 (1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996), a continuous inventory has been conducted with permanent plots established in forested areas with an intensity of approximately of one plot per km 2 . Each plot is re-measured approximately every 10 years. A total of 81,024 sample plots measured in the Third Spanish National Forest Inventory (SNFI3, 1997(SNFI3, -2007, were used as the input data for the EFDM (Figure 1). For validation purposes, the 13,667 plots measured to date in the Fourth Spanish National Forest Inventory (SNFI4, 2008-to date) were used (Table 1).  In each plot, trees are measured in four concentric circular subplots (5 m, 10 m, 15 m and 25 m radius) according to the diameter at breast height of the trees (dbh ≥ 7.5 cm, dbh ≥ 12.5 cm, dbh ≥ 22.5 cm and dbh > 42.5 cm, respectively). The following data are recorded in each plot: tree species, diameter at breast height, total height, and distance and azimuth from the plot center. SNFIs also record information on silvicultural practices, such as felled trees, soil preparation, pruning, and forest disturbances, such as forest fires or windthrows [35].
The EFDM is parameterized by matrices with dynamic and static dimensions or factors [36]. To analyze the dynamics of carbon storage in Spanish forests, SNFI plots were classified according to five static (forest type, known land-use restrictions, ownership, stand structure and bioclimatic region) and two dynamic factors (quadratic mean diameter and total volume) from SNFI field data recorded in each plot or using auxiliary cartographical information (Table 1). One of the factors that most affects growth and management is forest species composition. The SNFI plots have been classified into seven different forest types: Broadleaf forests, Conifer forests, Mixed forests, Dehesas (open woodlands), Other Conifer forests, Conifer plantations, and Broadleaf plantations.
With regards to known land-use restrictions, forested areas can be subdivided into forest available for wood supply (FAWS) and forest not available for wood supply (FNAWS) [37]. The concept "Forests Available for Wood Supply" (FAWS) is defined as "forest where any legal, economic, or specific environmental restrictions do not have a significant impact on the supply of wood"; i.e., all forests except those with administrative restrictions. Therefore, "Forests Not Available for Wood Supply" (FNAWS) refers to "forest where legal, economic or specific environmental restrictions prevent any significant supply of wood" [38].
Ownership was classified as public or private according to data from the Spanish Nature Data Bank. Stand structure was classified as even or uneven-aged according to the Spanish Forest Inventory Service. The map of biogeographical regions of Europe was used as a reference for classification [39]. Spanish forests belong to four biogeographical regions: Alpine, Atlantic, Mediterranean and Macaronesian ( Figure 1).
For forest wildfire scenarios, Vázquez de la Cueva et al. [33] developed estimates for future fire activity in peninsular Spain using regression models between monthly meteorological variables and the recorded fire activity in the period 1974-2005. With regards to results for the frequency of wildfires, a variable used in our modeling approach, the B2 scenario (characterized by medium-low greenhouse gas emissions, IPCC [40]) was found to cause a 2-fold increase and the A2 scenario (characterized by medium-high greenhouse gas emissions, IPCC [40]) a 2.5-fold increase in fire occurrence in peninsular Spain (see Vázquez de la Cueva et al. [33] for more detail).

EFDM Model
The EFDM [28] is a Markov chain and area-based matrix model, meaning that forest areas or strata (not trees or stands) are transiting between elements of a set of fixed states. The matrices represent forest areas classified according to ecological and socio-economic factors. The model simulates the development of the forest area distribution as a product of its initial state, proportions of areas expected to be managed according to different silvicultural practices, and the corresponding transition probabilities. The transition probabilities are conditioned by the activities, both of which can differ between factors. The parameterization of the Markov model is described in detail in Vauhkonen and Packalen [30].
Specific materials used in our simulations were: • Spanish Forest Map 1:50,000 [41] providing the cartographic base to assign the equivalent area to each plot.

•
The initial state (defined by the five static factors and two dynamics factors for each plot) was determined based on the SNFI3. The data processing involved the initial estimation and formatting of the data for the EFDM.

•
The dynamics of any-aged forest management were simulated by using the forest area classified in the quadratic mean diameter and volume matrix as an input for the EFDM, plus the five static factors. Starting from the statements provided by experts [28] and then following by an iterative process, static factors were weighted selecting those that minimized the error. The selected weights were as follows: forest type (1), bioclimatic region (0.4), known land-use restrictions (0.2), ownership (0.2) and stand structure (0.2).

•
Modeling the transitions due to natural processes was performed using pairwise observations from plots measured in both the SNFI2 and SNFI3 (51,676 plots out of 81,024 plots). These two data sets were also used to derive the transition probabilities matrix. Some plots, especially in the Atlantic bioregion, could not be included in the model because they were not re-measured in two consecutives inventories.

•
The activity probabilities were defined in two steps with two assumptions. First, the initial allocation of the harvests to the different types of forests was assumed to follow either the proportion of harvests carried out during the SNFI2-SNFI3 ten-year period ("business-as-usual allocation", ABAU) or the application of future harvests in accordance with the specific silvicultural recommendations ("schoolbook allocation", ASB) which are given in Serrada et al. [42] per species at national level. In a second step, the final values for the activity probabilities under both the allocations were obtained by iteratively adjusting the initial probabilities to produce the harvest levels in future, large-scale scenarios.
After running the EFDM, forest area for each time step in total volume, scaled-up to a hectare, is obtained for each stratum defined by the five static factors and two dynamic factors. From these results, standing volume, and fellings can be calculated.
Above-and below-ground biomass was calculated for all SNFI3 plots using species-specific allometric equations and diameter given by Montero et al. [43]. Biomass/volume ratios were then computed specifically for each stratum. Above-and below-ground biomass is calculated as follows: where W: above-or below-ground tree biomass (dry weight, Mg), R 1 : biomass/volume ratio Mg m −3 , and V: volume (m 3 ). Above-and below-ground carbon were also estimated through a similar process, using a species-specific carbon conversion factor (given by Montero et al. [43]) for the biomass previously calculated through allometric equations. Having calculated carbon/volume ratios for each stratum, above-and below-ground carbon is obtained: where C: carbon (Mg), R 2 : carbon/volume ratio Mg m −3 , and V: volume (m 3 ).
Having derived projections for the future development of forests under both management allocations according to present scenario wildfire frequency, the simulations are then repeated, varying wildfire frequency proportions by 2 and 2.5 times the present frequency for B2 and A2 scenarios, respectively, given by Vázquez de la Cueva et al. [33]. It was assumed that wildfire severity was always severe, so the forests affected were forced to transition to the beginning of the rotation.
Therefore, a total of six combinations of management regimes and wildfire scenarios were considered in the study ( Table 2). Business-as-usual allocation (ABAU) B2 3 Business-as-usual allocation (ABAU) A2 4 Schoolbook allocation (ASB) Present 5 Schoolbook allocation (ASB) B2 6 Schoolbook allocation (ASB) A2 To determine the accuracy of the model predictions, the relative bias (rBias, %) of the model was calculated [44] using the statistics: where est i : ith estimated value; obs i : ith observed value.

Business-as-Usual Allocation (ABAU) Management
A comparison of the changes in living biomass and carbon stock in Spanish forests as a result of the effects of wildfire is shown in Figure 2 and Table 3. If present climatic conditions (Combination 1) are considered, the simulation of total carbon stock reveals an increase from 629.0 M Mg in 2000 to 1115.4 M Mg in 2050 (+77%). In the same period, total biomass goes up from 1325.8 M Mg to 2366.7 M Mg (+78%), the growing stock increases from 962.0 M m 3 to 1844.0 M m 3 (+91%), and harvesting from around 75 to 183 M m 3 of roundwood over a ten-year period (+143%). Because biomass and carbon behavior are so similar, we limit the discussion below to the results for carbon.
Carbon sequestration varies significantly among static factors (Table 4) under the present scenario between 2000 and 2050. The greatest increases in carbon stock values (in percentages) are those for Other Conifer (+116%) and Conifer (+113%) forests, and the lowest values correspond to plantations (Broadleaf +44% and Conifer +41%) and Dehesas (+38%). We found similar carbon stock increments for FNAWS and FAWS (+78 and +77%), and ownership, with +79% for private forests and +75% for public forests. Even-aged forests (+97%) display a larger carbon stock increment than uneven-aged forests (+69%). Mediterranean areas present the highest increment, +95%, followed by Alpine areas (+59%). The increase in carbon stocks in Atlantic and Macaronesian areas is far lower, at +38% and +23% respectively.    The higher fire frequency under the B2 and A2 scenarios (2 and 2.5 times the present frequency, respectively) results in a declining trend in the carbon stocks. The B2 scenario (Combination 2) shows a total carbon stock of 972.8 M Mg (+55%) at the end of the simulation, which is 22.7% lower than was simulated in the present scenario. Under the A2 scenario (Combination 3), the final value is 908.9 M Mg (+44%), which is 32.8% lower than under the present conditions scenario.
The most affected forest type under the two climatic change scenarios which consider increased wildfire frequency (B2 and A2) is Mixed forests (−27% and −39%) followed closely by Broadleaf forests, Conifer forests and Broadleaf plantations (−25% and −36%). In contrast, Dehesas (−1.6% and −2.4%) display greater stability in the face of changing fire occurrence (more details in Appendix A, Figures A1-A7). As regards wood supply, there is little variation in availability, with carbon stocks decreasing in 2050 by around −23% for scenario B2 and −33% for scenario A2 for FAWS, whereas the decrease is −24% and −35% in the case of FNAWS. Changes in carbon stocks are very similar for private and public forests, with around −23% under the B2 scenario and −33% for the A2 scenario. Even-aged forests seem more affected by forest fires (−25% and −37%) than uneven-aged forests (−21% and −31%). The Macaronesian region is the most affected, the decrease being −32% under the B2 scenario, and −44% under the A2 scenario. The Alpine region, with −5% and −7% under the B2 and A2 scenarios, appears to be the least affected by changes due to forest wildfires (Table 4).  Table 5 show carbon stock for time series over the present and future 50-year period for the three different wildfire scenarios under ASB management. Over recent decades, ABAU in Spain has been lower than that associated with the silvicultural recommendations, so ASB management means an increase in harvests and consequently, a decrease in growing stock.

Schoolbook Allocation (ASB) Management
Simulation of total carbon stock for Spanish forests under ASB management (Combination 4) indicates an increase from 642.1 M Mg in 2000 to 1103.7 M Mg in 2050 (+72%). Compared to ABAU, the ASB resulted in a −3.0% drop in growing carbon stock (above-ground and below-ground carbon) in 2050 (from 1044.3 M Mg C to 1013.1 M Mg C). The accumulated harvest in a ten-year period increases from around 75-183 M m 3 to 102-230 M m 3 , an increase of 36% in the first ten-year period and 26% in 2050.
The simulations of climate scenarios reveal lower carbon stocks under both scenarios, although the decrease is a little less than under the business-as-usual allocation. Under the B2 scenario (Combination 5), the total carbon stock is predicted to be 968.6 M Mg (+51%) by the end of the simulation, this being 20.7% lower than for the simulation of the present condition scenario. Under the A2 scenario (Combination 6), carbon stocks reach 908.5 M Mg (+42%), although this is 29.9% lower than under the present scenario.
Considering static factors, the results are similar to the business-as-usual allocation, but there are slight differences. As for forest types, the B2 and A2 scenarios have a greater effect on Broadleaf, Conifer, Mixed forests and Broadleaf plantations (around −24% and −34% in all cases). The two scenarios have little if any influence on Dehesas (−1% and −2%). FNAWS are more sensitive (−26% and −38% carbon stock in 2050 for B2 and A2 scenarios) than FAWS stands (−20% and −29%). As regards ownership, the results point to a similar decreasing trend in carbon stock under both scenarios, around −20% under the B2 scenario and −30% under the A2 scenario. Under the B2 scenario, even-aged and uneven-age forest carbon decreases by around −21%. This similarity is maintained under the A2 scenario, with around −30% in both cases. As occurred under ABAU, the Macaronesian area is the most affected under both the B2 and A2 scenarios (−33% and −45%), followed by the Mediterranean (−23% and −33%) and Atlantic (−19% and −27%) areas. The Alpine area presents a much smaller decrease, with −4% under the B2 scenario and −6% under the A2 scenario (Table 6). Table 7 summarizes all the differences between management allocation in total carbon stock in 2050 considering static factors and the different scenarios. With the exceptions of FNAWS (as mentioned, with a strict no felling policy) and the Macaronesian region, all cases show slightly better results with the schoolbook allocation than with the business-as-usual allocation. Schoolbook allocation shows a total improvement of 2% in Present-B2 scenarios and 2.9% in Present-A2 scenarios against business-as-usual allocation. The greatest differences were found in Mixed forests (3.4% and 5%), Conifer plantations (3.7% and 6%) and even-aged stands (4.2% and 6.3%), respectively.     Table 7. Total carbon stock differences (%) between management allocations in 2050 per static factor and wildfire scenarios in Spain.  Table 8 presents the relative bias rBias (%) when EFDM is validated using the 13,667 plots currently available belonging to the SNFI4 using values from the SNFI2 as the initial state (20-year prediction). If we consider the total forest area, EFDM underestimates the total volume (M m 3 ) by 2.2%. Analyzing forest types, Mixed forest stands out with a high relative bias (11.5%). For the rest of the forest types, except for Broadleaf plantations (7.7%), the relative bias is less than ±5%. Wood supply and ownership present a similar error (1.7-2.2%) in all cases, while for stand structure there are differences between even-aged (5.5%) and uneven-aged forests (0.2%). The Macaronesian area shows the greatest underestimation (13.4%), followed by the Atlantic area (8.4%). The relative bias for the Alpine (0.4%) and Mediterranean (−1.9%) areas are much lower.

Discussion
The Kyoto Protocol [5] requires every industrialized country to have a transparent and verifiable method for estimating the size and evolution of the carbon stored in forest ecosystems. More recently, the Paris Agreement [6] encourages action to be taken to implement and support both policy approaches and positive incentives for activities relating to reducing emissions. Alternative policy approaches are also proposed, such as joint mitigation and adaptation approaches for the integral and sustainable management of forests. This is the first study to simulate the future living biomass and carbon stocks of the Spanish forests under any-aged management, using an area-based matrix model that considers different/alternative management and wildfires scenarios. The study focuses on harmonized definitions, assumptions, and methodology to account for the administrative restrictions affecting forest use.
A significant interaction was found between the disturbance effects of higher wildfire frequency (B2 and A2 scenarios) and decreasing forest management, which supports the hypothesis that warming and abandonment of forest management have a synergistic effect on the forest carbon balance, as reported in previous studies conducted in Spain [45,46]. The thinning response hypothesis [47] asserts that thinning does have a significant influence on stand volume stock at the end of the rotation period for a wide range of thinning intensities or stocking densities, whereas heavier thinning beyond the considered range reduces volume stock. Hence, although there are no large differences in total carbon stock, our results reveal the potential benefit of forest management regarding biomass and carbon stocks. Moreover, Ruiz-Peinado et al. [46] conclude that in Mediterranean areas, forest management is an effective way to maintain and enhance high carbon sequestration rates. Vauhkonen & Packalen [36] states that the carbon stored by Finnish forests can be higher by applying less intensive management systems, but the extensive shifts from conventional even-aged management to alternative management systems exhibited a strong impact on harvesting costs. In eastern North America (State of New York), a multicriteria decision analysis approach identified and quantified trade-offs among three important ecosystem services-timber, carbon and habitat quality-under commonly used management scenarios. The most intensive management yielded greater timber volumes but resulted in the weakest carbon and habitat quality scores. For carbon storage, No Management resulted in higher utilities than the two other scenarios [48].
One of the objectives of sustainable forest management is to reduce the competition among trees through silvicultural operations and release resources to the remaining trees [49]. In Spain, there has been a progressive abandonment of forest management operations over recent decades, which could have increased competition for resources (mainly water) and increased the vulnerability of forests to climate change [50]. With regards to the findings for the ASB management scenario, several studies have concluded that regeneration methods and thinning treatments that maintain a large proportion of mature trees are better, in terms of maintaining carbon stores, than those involving more intensive removals [51]. However, such treatments can provide the opportunity for long-term off-site carbon storage in forest products made from the felled wood [23]. In contrast, extended rotations with high stocking levels to maximize on-site stores lead to the lowest levels of live carbon increment [11] and higher carbon loss in the case of wildfire disturbance.
Regarding forest types, the results show Mixed forests to be one of the most affected by higher wildfire occurrence, although its complexity in terms of forest structure may foster self-regulation and provide greater adaptability to cope with increasing uncertainty due to global change [52]. Nevertheless, the stability of Dehesas is noteworthy under all scenarios, with a −2% variation under both management approaches. The open distribution of the mature oak tree population in the Dehesa woodlands protects them from the risk and severity of wildfires. The Spanish Dehesa is a noteworthy example of a managed, anthropogenic, Mediterranean ecosystem. Dehesas are agrosilvicultural open woodlands comprised mainly of Quercus sp., with holm oak (Quercus ilex L.) being the most common species, followed by Quercus suber L. These oak stands are not suited to traditional, intensive forestry because of the poor sandy soils and variable, dry Mediterranean climate [53,54].
Our findings show that the Macaronesian area (located in the Canary Islands) is the most affected by wildfires, with almost four times more occurrence than in the other three bioclimatic areas. According to the SNFI3, almost 12% of the plots located in the Macaronesian area experiences a wildfire between 1990 and 2000, as opposed to less than 3.5% in the other three areas. Our model shows that the effect of increasing forest wildfire occurrence on the Macaronesian area is particularly high, with carbon stock decreasing by −32% and −44% in 2050 under B2 and A2 scenarios, respectively. The frequency of wildfires in the Canary Islands has risen since the arrival of the human population to the islands [55] and the number of anthropogenic related fires quadrupled between 1970 and 2010 [56]. Burned subtropical cloud forests are more prone to increased wildfire events due to a greater presence of pioneer species, a higher density of trees, climatic variables which tend to range more widely throughout the day, and the long time span needed for recovery to a pre-fire state [57]. However, the Canary Island pine, which accounts for more than 30% of Macaronesian forest area, is characterized by its fire adaptation strategies [58].
The simulation results underestimate the values for all forest types in the validation data set. Simulation from SNFI2 data for the total area had a mean volume of 2.2% less than the SNFI4 real data (20-year simulation). Higher relative biases correspond to forests with broadleaf species (Broadleaf plantations and Mixed forests). Broadleaf plantations are particularly difficult to project because the rotation is usually around 12 years, close to the ten-year measurement carried out in the SNFI. Moreover, not all forest types considered, such as Dehesas or Other conifer forest, are well represented in the validation database. Other studies using the SNFI to analyze forest dynamics and ecological processes covering a wide regional or national scale, including high variability, also found low or moderate predictive ability [53,59,60].
The production of forest biomass as a carbon-neutral alternative for renewable energy production [61] can be increased by using specific biomass crops or by changing forest management in plantations of exotic or native species [62]. Broadleaf and Conifer plantations together make up just over 7% of the total area and accounted for almost 12% of total carbon in 2000. They are not greatly affected by forest wildfire scenarios under ABAU, which suggests that the production of biomass from Eucalyptus species as well as products from other fast-growing species used in plantations such as Pinus pinaster Ait. and Pinus radiata D. Don, will be maintained.
Our simulations considered carbon stocks in above-and below-ground forest biomass, focusing on standing living trees included in the Spanish NFIs. Tree biomass represents one of the main carbon sinks in forests [63], but shrub biomass could be also important in terms of carbon sequestration, particularly in the context of the Mediterranean basin [64]. In Spain, the amount of area covered by shrubs accounts for more than 35% of the forested land [65] and will likely increase in the future due to the abandonment of crops and disappearance of forests, caused by increased temperatures and drought and the probable increase in the occurrence of forest fires [64]. Other important carbon reserves not included in the current study are soils, litter, and dead wood, which should be included in the future to correctly model the total capacity of forests as carbon sinks.
In this research, the fire hazard under different climatic scenarios was accounted for using the forecasts reported by Vázquez de la Cueva et al. [33]. These authors only considered the direct effects of climate change on wildfire activity for different IPCC climate change scenarios, although any conclusions will be limited by not including the possible effect of other factors, besides that of climate change, on fire occurrence probability. Weather and climate are among the main factors influencing wildfire potential, but human activities such as land-use change [66] and socio-economic factors [67] greatly influence fire regimes. Additionally, fire has other direct effects on the vegetation that we have not considered, such as species distribution, migration, substitution or extinction [68]. Furthermore, apart from wildfire risk and severity, changes in climate are expected to affect other ecological processes in the Mediterranean basin which are not considered in this research. Some examples of these processes include regeneration, growth and mortality rates [69].
The EFDM is a well-developed research tool to support the assessment of forest policy effects through management and the generation of scenarios for sustainable management of forests in Europe at sub-national, national and international level [28]. Nevertheless, it should be stressed that EFDM simulation results are assumed to be deterministic, so uncertainties associated with the state of the space, activities, transition probabilities, and output coefficients are not considered [28,30]. Additionally, the EFDM assumes fixed cycles (10 years in our study), but the measurements can often go on for longer or shorter periods than that.
According to Regulation 2018/841 of the European Union, the forest reference level (FRL) must be estimated by each EU Member State considering similar forest management practices to those of the period 2000-2009. Forsell et al. [70] state that there is an interval of the FRL in which predictions may vary depending on several factors such as the ecological bioregions. They also show that assumptions related to climate change, such as forest fires and the allocation of forest management practices, have the greatest impacts on the FRL estimates, while the start year of the projection or the stratification of the forest have lower impacts on the FRL estimates. Our results show that national projections can differ considerably due to the effects of disturbances depending on the different scenarios considered, especially for the Macaronesian and Mediterranean bioregions.
Because of the harmonized approach, the results can be easily compared with other projections at European Member State level.  presented future projections of the forest above-ground carbon in 23 European countries using EFDM. Their results varied between 105.6 Mg C ha −1 (Switzerland) and 31.8 Mg C ha −1 (Finland) in 2040, with a growth ratio up to 59% (Ireland) between 2015-2040. Our results show 45.4 Mg ha −1 of above-ground carbon in 2040 and a growth ratio of 29.2% between 2015-2040. Furthermore, they can also be compared with those computed for the same area using less harmonized approaches such as national projection models. Previous studies have focused on quantifying biomass and carbon stocks in Spanish forests. Montero et al. [43] developed the first species-specific allometric equations for above-ground and below-ground fractions to calculate tree biomass and carbon content. Their calculation was 1806 M Mg CO 2 (492 M Mg C) in 1990 (from SNFI2 data) and an estimation of 2858 M Mg CO 2 (779 M Mg C) in 2004 based on annual carbon increments. A study by González-Diaz et al. [71] quantifies the mean value for total carbon stored as 43.35 Mg C ha −1 for the SNFI3, dated in 2000, using the same species-specific allometric equations developed by Montero et al. (2005). Both estimations are slightly higher than those of our simulations, which were 629.0 M Mg C (38.6 Mg C ha −1 ) in 2000 and 764.9 M Mg C (46.9 Mg C ha −1 ) in 2010, where biomass/volume and carbon/volume ratios specific to static and dynamic factors were applied to calculate biomass and carbon instead. Nevertheless, the national living biomass carbon value reported by Spain for FOREST EUROPE [3] for 2010 was 569.56 (423.6 M Mg C for above-ground and 140.1 for below-ground), which is lower than our estimations. This comparison may provide useful insights into the influence that harmonizing may have on both national and international policy and decision-making.

Conclusions
The EFDM harmonized approach allows us to simulate the capacity of forests to store carbon under different management and wildfire scenarios at the national scale in Spain. The results can be easily compared with other projections either at individual European Member State level, groups of countries, or at European level. However, the EFDM underestimated the total volume (m 3 ) by 2.2% in the tested validation for a 20-year simulation.
Since this is the first time that this type of model has been applied to Mediterranean forests, the consideration of forest wildfires has been key to achieving more realistic simulations. The most stable forest type under both wildfire scenarios would be Dehesas. Considering bioregions, the Macaronesian area is the most affected and the Alpine region, the least affected.
Regarding forest management regimes, the results show the potential benefit of forest management as regards biomass and carbon stocks, with slightly better results in ASB management compared to ABAU management.
The results for projected biomass and carbon stock variability under different forest wildfire and management scenarios provide valuable information, permitting optimal decision-making for the forest-based bioeconomy and ecosystem services and consequently, for forest-related policies.