Site Quality Models and Fuel Load Dynamic Equation Systems Disaggregated by Size Fractions and Vegetative States in Gorse and High Heath Shrublands in Galicia (NW Spain)

: Compatible model systems were developed for estimating fuel load dynamics in Ulex europaeus (gorse) and in Erica australis (Spanish heath) dominated shrub communities at stand level. The models were based on intensive, detailed destructive field sampling and were fitted simultaneously to fulfill the additivity principle. The models enable, for the first time, estimation of the biomass dynamics of the total shrub layer, size fractions and vegetative stage, with reasonably good accuracy. The approach used addresses the high variability in shrub biomass estimates by using a site index (SI) based on biomass levels at a reference age of 10 years. Analysis of the effect of climatic variables on site index confirmed the preference of gorse for mild temperatures and the ability of high heath communities to tolerate a wider range of temperatures. In the gorse communities, SI tended to increase as summer rainfall and the mean temperature of the coldest month increased. However, in the heath communities, no relationships were observed between SI and any of the climatic variables analyzed. The study findings may be useful for assessing and monitoring fuel hazards, updating fuel mapping, planning and implementing fuel reduction treatments and predicting fire behavior, among other important ecological and biomass use-related applications.


Introduction
Shrublands play an important role in ecological processes that are essential for life on earth, such as carbon sequestration, nutrient recycling, hydrological regime and soil protection, and they also provide biodiversity, wildlife habitats and bioenergy for human use [1][2][3].On a global scale, the area occupied by shrubland has increased greatly as a result of various factors, including the reduction in livestock and agricultural activities (due to rural abandonment) and more extreme temperatures and prolonged droughts associated with climate change [4][5][6][7].The increased frequency of forest fires and the consequent transformation of forests into shrubland has also become a widespread phenomenon [8,9], of particular importance in Mediterranean ecosystems [10,11].Many Mediterranean-type shrublands are considered fire-prone vegetation [12,13], in which fire is both a disturbance and a first-order ecological factor in the evolution and dynamics of vegetation and, therefore, of shrubs [14][15][16].On the one hand, shrub fuels provide the energy source for wildfires while fuel attributes are determinants for fire behavior and its effects, modulating fire intensity and severity [17][18][19].In fact, shrub productivity regulates fuel availability, driving, together with climate and ignition sources, and fire activity in shrublands [20,21], with fuel load being a determinant of fire intensity [22].The age of the shrubland can also modify the flammability of the community through changes in fuel structure [12,23].On the other hand, the fire regime, associated with productivity gradients, affects vegetation evolution and drives the development of adaptive traits in plants, favoring certain compositions of shrub communities [24,25].Thus, a complex feedback process can be established between shrub fuel dynamics and fire regime [26,27].Fire behavior is driven by the relationships between meteorological conditions, topography, human intervention and, most importantly, fuel characteristics [28].Studying the characteristics of fuels that affect fire behavior and the impact on the ecosystem is of particular interest as these are the only landscape elements that can be managed [29,30].Considerable efforts have thus been made in recent years to characterize fuel structure, in particular, to determine the distribution of fuel loads in different size fractions and in relation to physiological state [7,31,32].
Since fuel attributes change with age, knowledge of the dynamics of fuel load fractions of shrublands at the community level is essential to evaluate the effectiveness of fuel management treatments, prioritize fuel risk reduction actions and optimize treatment scheduling [33].Fuel dynamics is also important to assess the influence of climate and other environmental factors on that vegetation type, for a proper estimation of climate change impacts on vegetation and fire regimes and to design adaptation strategies [34][35][36].
However, modeling fuel dynamics in shrub ecosystems at the stand level is considered a challenge [34], since most of the available fuel consists of living vegetation, comprising species of different characteristics that are often more influenced by the environmental conditions of the site than forest ecosystems.Changes in community composition with age are also common, affecting fuel properties and complicating modeling.The stand-level approach in studies of fuel load dynamics in shrublands is critical to operational objectives; since in shrubland fires, the fuel is the whole vegetation pool rather than individual plants.
Remote sensing methods provide invaluable information on biomass accumulation and growth in shrublands and are the only viable alternative for fuel mapping [32,[37][38][39][40][41].However, detailed field inventories remain critical for knowledge of fine-scale shrub structural variables and for testing and refining remote sensing-based approaches [42,43].
By adapting the classic definition of site quality in forest science [44] to the study of shrublands, this index can be defined as the net biomass productive potential of a site for a particular species or shrubland community.Any index used to define site quality must fulfill a series of basic requirements: (1) it must be quantitative so that it provides a reference point for management diagnosis and prescription; (2) it must be objective and (3) it must be easy to determine [45].
In dominant height-age curves, which are commonly used to assess site quality, especially in pure and even-aged stands, the stand site index is defined as the value of the dominant height in the specific curve for that stand at a given reference age [46].However, the concept of dominant height is not easy to adapt to shrubland communities, especially in dense formations, in which individual measurements are difficult to carry out.Thus, the use of fuel load-age models seems a more appropriate choice in this case.
Efforts to model shrubland fuel load dynamics at the stand level have been limited, in terms of the types of shrublands examined, and have focused primarily on total fuel or some specific fraction [39,[47][48][49][50][51], including in some cases live or dead fuel [35,52,53].Dynamic models that consider total fuel load and various fuel fractions, differentiating by size and physiological state, at the stand level are even scarcer [54][55][56].
Although shrublands dominated by Ulex europaeus L. (gorse) cover an increasing area worldwide, there is a large gap in knowledge about the fuel dynamics of this species.The same also applies to Erica australis L. (Spanish heath), which is one of the main components, together with gorse, of the dry heathlands in the Iberian Peninsula (Western Europe).Both of these species play a key ecological role in this region and are involved in severe and recurrent wildfires.Gorse is considered a hazardous fuel due to its high biomass growth [1,57,58] and high aerial fuel flammability [59][60][61], favored by the production of very flammable volatile organic compounds [62].Overall, gorse-dominated stands are qualified as very fire-prone shrublands [63], which are capable of generating intense fires [64][65][66].Both, gorse-and high heath-dominated communities have high energy potential [67,68], and continuous stands dominated by high heath are also considered very flammable and prone to fires [69,70].
In Galicia (northwestern Spain), shrublands cover approximately one-third of the forest area [71], with gorse-and high heath-dominated communities being the most abundant with 77% of the total shrubland area.These formations are involved in approximately two-thirds of the forest fires occurring in the region [72]; so, knowledge of fuel dynamics in these communities is a priority in fire management.
In summary, the overall aim of this work was to develop site quality curves (fuel load-age relationships) for two of the main shrubland communities in NW Spain by using an extensive database derived from destructive sampling.The specific aims were as follows: to determine the site index of each sampling plot by using the previously developed curves; to evaluate the effect of topographic and climatic variables on the site quality; and, finally, to fit a compatible system of equations to estimate fuel loads disaggregated by size fractions and vegetative state (live or dead) for each shrubland community using site index and age as independent variables.

Study Area, Shrub Communities and Inventory Plots
The study was carried out in Galicia (NW Spain), where shrublands represent about 20% of the total area of land and 30% of the forest land area in the region.Most shrublands in Galicia are dry heaths [73] and are included in habitat 4030 of the European habitat classification (Council Directive 92/43 CEE), particularly in the Atlantic heathlands and Ibero-Atlantic heathlands categories.Gorse-dominated communities (Ulex sp.) cover about 55% and heath-dominated communities (Erica spp.) about 22% of the shrubland areas.
Two shrub communities were considered in the present study: a gorse-dominated community (Ue) and a high heath-dominated community (Ea).Both are primarily composed of perennial, multi-stemmed (or highly branched) evergreen woody species (shrubs and sub-bushes), which typically form dense (a large number of plants per area) and closed (high coverage) stands ranging in height from moderate to moderately high (0.5-3 m).
The dominant species in the gorse-dominated community is Ulex europaeus L., and other main secondary species are often present, such as Ulex gallii Planch., U. Minor Roth, Erica umbellata Loefle ex L., E. Cinerea L., Pteridium aquilinum (L.) Kuhn in Kersten and Pterospartum tridentatum (L.) Willk.The high heath-dominated community includes three dominant species (Erica australis L., Erica arborea L. and Erica scoparia L.) with some common structural characteristics allowing them to be grouped together.Other main secondary species frequently present in the high heath-dominated communities include Pterospartum tridentatum, Halimium alyssoides (Lam.)Greuter, Ulex europaeus and Pteridium aquilinum.
Sites occupied by these two shrub communities were randomly selected on the basis of the information on the treeless shrub-covered polygons provided in the Spanish Forest Map 1:25,000 [74].The number of sampling sites for each shrub community was approximately proportional to the area covered by each in Galicia [74].In each site selected, a circular sampling plot was established, and a total of 316 plots were inventoried (191 of gorsedominated communities and 125 of high heath-dominated communities).

Biomass Sampling
From the center of each circular sample plot, a random azimuth was generated to establish one diameter (20-30 m length, depending on shrub height), and another diameter was drawn perpendicular to the first.Four destructive sampling subplots (quadrats) were located in the center of the four plot radii corresponding to the diameters (Figure 1, right).The area of each quadrat ranged from 4 to 36 m 2 , depending on shrub height: For shrubs

Biomass Sampling
From the center of each circular sample plot, a random azimuth was generated to establish one diameter (20-30 m length, depending on shrub height), and another diameter was drawn perpendicular to the first.Four destructive sampling subplots (quadrats) were located in the center of the four plot radii corresponding to the diameters (Figure 1, right).The area of each quadrat ranged from 4 to 36 m 2 , depending on shrub height: For shrubs smaller than 1.0 m in height, 4 m 2 quadrats were destructively sampled; for shrubs taller than 1.0 m, the quadrat size varied from 3 × 3 m to 6 × 6 m.According to the approach proposed by Pearce et al. [48] for biomass sampling in continuous dense shrub communities, each quadrat was delimited by four wooden poles, and a linear transect was laid out (with a tape) along the perimeter and one diagonal of the quadrat.The vegetation growing inside the quadrat was carefully clipped along the lateral boundaries, taking care to exclude portions of plants growing within the quadrat but hanging outside the boundary.The horizontal lengths of the intercepted crown of the standing shrub species were measured (cm) with a graduated tape, along the transect.The values were used to determine the linear cover by the shrubs in terms of the percentage of the transect length intercepted (maximum 100%).Shrub height was determined as the vertical distance (cm) between the surface of leaf litter and the top of the plant canopy and was measured with a graduated tape every 50 cm along the transect.All vegetation portions in the vertical projection of the sampling quadrat area were carefully harvested at ground level and placed in bags, which were labeled appropriately and transported to the laboratory.
The shrub and litter biometric measurements in each of the four sampling quadrats and frames, respectively, in each circular plot were averaged to obtain a value per inventory plot of shrub cover (CovShr) and mean shrub height (ℎ ).According to the approach proposed by Pearce et al. [48] for biomass sampling in continuous dense shrub communities, each quadrat was delimited by four wooden poles, and a linear transect was laid out (with a tape) along the perimeter and one diagonal of the quadrat.The vegetation growing inside the quadrat was carefully clipped along the lateral boundaries, taking care to exclude portions of plants growing within the quadrat but hanging outside the boundary.The horizontal lengths of the intercepted crown of the standing shrub species were measured (cm) with a graduated tape, along the transect.The values were used to determine the linear cover by the shrubs in terms of the percentage of the transect length intercepted (maximum 100%).Shrub height was determined as the vertical distance (cm) between the surface of leaf litter and the top of the plant canopy and was measured with a graduated tape every 50 cm along the transect.All vegetation portions in the vertical projection of the sampling quadrat area were carefully harvested at ground level and placed in bags, which were labeled appropriately and transported to the laboratory.
The shrub and litter biometric measurements in each of the four sampling quadrats and frames, respectively, in each circular plot were averaged to obtain a value per inventory plot of shrub cover (Cov Shr ) and mean shrub height (h Shr ).

Laboratory Work
In the laboratory, the material of the standing shrub stratum was physically separated by size class into fine fuels (diameter < 0.6 cm, hereafter G1), medium fuels (0.6 cm ≤ diameter < 2.5 cm, hereafter G2) and coarse fuels (2.5 cm ≤ diameter < 7.5 cm, hereafter G3), with a go-no go gauge [75].The material was further subdivided on the basis of condition (live or dead), determined by visual inspection.The size categories were selected according to their surface-to-volume ratios and, for dead fuels, because they coincide with the size ranges defined by Fosberg et al. [76], due to the contrasting water absorption and desorption rates (time lags of 1 h, 10 h and 100 h, respectively).These three size ranges have also been used to construct custom fuel models for predicting fire behavior [77][78][79][80].Once classified, the material was weighed and dried in forced air-drying chambers (105 • C for 24 h for fine fuels and 48 h for coarse fuels) for determination of the dry biomass of each fraction.The fuel load of each fraction was obtained by dividing the dry biomass by its corresponding sampling area.Thus, five different loads of the various biomass fractions were computed: W Shr_G1_dead = dead fine shrub load, W Shr_G1_live = live fine shrub load, W Shr_G1 = fine shrub load (dead + live), W Shr_G23 = coarse shrub load and W Shr = total shrub load = W Shr_G1 + W Shr_G23 = AGB at stand level.Fractions G2 and G3 were grouped to prevent loss of data, as fraction G3 is infrequent in the juvenile stages of these shrub communities.
The age of the shrub communities was estimated from the date of the last natural or anthropogenic disturbance that resulted in complete regeneration, corroborated by counting annual growth rings on the basal section of randomly selected stems of individual plants in each sampling plot.The basal sections were polished with fine-grade sandpaper.Those plots in which it was not possible to reliably assign an age were excluded from the analysis.Finally, a total of 267 sample plots (165 of gorse-dominated communities and 102 of high heath-dominated communities) were used in this study.The basic descriptive statistics of shrub strata for the main structural characteristics of each shrub community are given in Table 1.
Table 1.Mean values of the standing shrub fuel strata characteristics.Std.dev.= standard deviation, n = number of plots, h Shr = shrub height, Cov Shr = shrub cover, W Shr = total shrub fuel loadW Shr_G23 = coarse shrub fuel load, W Shr_G1 = fine shrub fuel load, W Shr_G1_dead = dead fine shrub fuel load, W Shr_G1_live = live fine shrub fuel load.See definitions in the text.Ea = high heath-dominated communities and Ue = gorse-dominated communities.

Development of Total Fuel Load Growth Curves
W Shr -t curves were constructed on the basis of the information obtained in the sampling plots.This type of relationship has long been used in forestry as a proxy for site productivity, mainly in forest stands, giving rise to site quality or site index curves [44,81,82].According to Clutter et al. [44], most techniques used to construct site quality curves can be viewed as special cases of three general methods: (1) the guide-curve method, (2) the difference-equation method and (3) the parameter-prediction method.The latter two methods require remeasurement data from the sample plots; however, as we only have one measurement of W Shr and t for each sample plot, the site quality curves were developed using the guide-curve method with a growth function as the base model.
Growth functions describe variations in the global size of an organism or a population with age; they can also describe the changes in a particular variable of a tree, stand or shrubland with age, in this case, total shrub load (W Shr ).Numerous growth functions can be used in forestry, such as the 74 documented by Kiviste et al. [83].The following are the most important desirable attributes for site quality equations [84][85][86]: (1) polymorphism, (2) sigmoid growth pattern with an inflexion point, (3) horizontal asymptote at old ages, (4) logical behavior (e.g., W Shr should be zero at age zero), ( 5) theoretical basis and ( 6) baseage invariance.Fulfillment of these attributes depends on both the construction method and the base model used to develop the curves, and it cannot always be achieved.
Multiple asymptotes are also a desirable attribute [87,88], although some of the most frequently used functions have a common asymptote.This does not appear to be of great importance, as the curves are generally suitable for the range of ages that would be used in practice, and the common asymptote is usually achieved at very old ages.
where W shr is the total shrub load; t is the age considered as the time elapsed since the last natural or anthropic disturbance and a iH , a iK and a iB are parameters to be estimated.
In the first step, the parameters a i of each base model are estimated by non-linear regression using the entire W Shr -t database, obtaining the so-called guide curve, which is the "average" line representing the data used.The same number of families of curves as there are parameters in the base model (two in the case of the Hossfeld model and three in the other two cases) can be obtained from the guide curve.The curves of a given sample plot for each family are obtained from the guide curve equation by varying one site-dependent parameter (W Shr -t values of the sample plot) and holding the others constant.
Selection of the best base model and the best family of curves of this model was based on the analysis of the values of the model efficiency (ME, Equation ( 4)) and root mean squared error (RMSE, Equation ( 5)) goodness-of-fit statistics, on the fulfillment of the previously mentioned desirable attributes and on the visual inspection of the families of curves when plotted on the W Shr -t data of the sampling plots.
where Y i , Ŷi and Y are the observed, predicted and mean values of the dependent variable, and n is the number of observations used to fit the equation.

Compatible System for the Shrub Fuel Complex
Once the equation has been fitted to obtain the total fuel load growth curve as a function of age, the next step is to obtain the equations for the remaining fractions.For this purpose, the site index (SI) of each sample plot is determined as the value of the total fuel load (kg m −2 ) at a reference age of 10 years.
Equations were developed for estimating the load of each of the four fractions of the shrub fuel complex (W Shr_G1_dead , W Shr_G1_live , W Shr_G1 and W Shr_G23 ) for each community.Allometric models (y = b 0• X b i i ) for estimating fuel loads were tested for all the biomass fractions by considering the site index (SI) and age (t) as independent variables to be tested.
Allometric equations must fulfill the property of additivity, i.e., the sum of biomass predictions from separate fuel fractions must equal the biomass prediction from the total biomass model (the sum of W Shr_G1_dead and W Shr_G1_live estimates must equal W Shr_G1 estimates or the sum of W Shr_G1 and W Shr_G23 must equal W Shr estimates).Therefore, in the first step, the equation of each fuel fraction of each shrub community was fitted separately, and the complete system of four equations (one for fraction) was then fitted simultaneously for each shrub community to guarantee additivity: Two equations that discriminated between fine (W Shr_G1 ) and coarse fuel loads (W Shr_G23 ) by disaggregating W Shr equation were fitted as follows: The equation for estimating the fine fuel load was then obtained as follows: with b i = b ig23 − b ig1 ; and the equation for estimating the coarse fuel load was as follows: The two equations discriminating between dead fine (W Shr_G1_dead ) and live fine fuel loads (W Shr_G1_live ) loads by disaggregating Equation (7) were then fitted as follows: The equation for estimating the fine dead fuel load was then obtained as follows: with c i = c ig1_live − c ig1_dead ; and the following equation was fitted to estimate the live fine fuel load: Fire 2024, 7, 126 8 of 25 The condition number was used to evaluate the presence of multicollinearity among variables in the fitted equations.According to Myers [93], values higher than √ 1000 indicate problems associated with multicollinearity.The presence of heteroscedasticity was analyzed by the White test [94] and by visual inspection of studentized residuals plotted against fitted values.When heteroscedasticity was detected, each observation was weighted by the inverse of its estimated variance ( σ2 i ), assuming that this variance can be modeled as a power function of the independent variables [95], i.e., σ2 i = (X i ) k .The value of the exponential term k was optimized to provide the most homogeneous studentized residual plot by using the method proposed by Harvey [96].
The systems were fitted using the nonlinear seemingly unrelated regression (NLSUR) method, which considers the cross-equation correlations, and the cross-equation error covariance matrix obtained by ordinary least squares was used to initiate the iterative procedure.When necessary, the weighting factor for heteroscedasticity was programmed in the MODEL procedure of SAS/ETS ® [97].
The accuracy of estimates was evaluated using the values of the model efficiency (ME) and root mean square error (RMSE) goodness-of-fit statistics.

Assessing the Effect of Topographic and Climatic Variables on Site Quality
An explicit relationship between site index values and site environmental conditions, defined by a set of topographic and climatic variables, was evaluated.For this purpose, the elevation and slope of each sample plot were obtained from the field sampling.Mean values of annual temperature (T), temperature of the warmest month (Twm), temperature of the coldest month (Tcm), annual precipitation (P) and summer precipitation (Ps) of each of the 267 sample plots were obtained from a climatic map (spatial resolution, 1 km) constructed by Gonzalo [98] by ordinary cokriging of data from 5426 weather stations (series 1951-1999).Moreover, mean rainfall days per month (Rd) and per summer months (Rds) were obtained for each sample plot from the climatological series database of the nearest meteorological station in Galicia (https://www.meteogalicia.gal,accessed on 14 December 2023).
Additionally, the annual aridity index (P/T) and the summer aridity index (Ps/Twm) were derived from the original data for each of the 267 sample plots (Table 2).Elevation (m.a.s.l.); Slope (%); P annual rainfall (mm); Ps summer precipitation (mm); T annual temperature ( • C), Twm mean temperature of the warmest month ( • C); Tcm mean temperature of the coldest month ( • C); Ps mean monthly summer rainfall (mm); Rd number of rainy > 1 mm per month; Rds number of rainy days > 1 mm in summer per month.Different letters indicate significant differences in mean values (α = 0.05) between shrubland communities.The modeling approaches used should enable the development of parsimonious models, preventing overfitting and resulting in relationships, where ecological coherence can be easily assessed.The multivariate adaptive regression splines (MARS) approach was used to model the relationship.This approach has been successfully used for fitting this type of relationship to establish forest stand site quality [99,100].
MARS is a nonparametric technique based on fitting piecewise linear regressions by intervals of the independent variable space [101].The generalized cross-validation Fire 2024, 7, 126 9 of 25 (GCV) criterion was used to select the best set of independent variables using the earth package [102] implemented in R software 4.2.2 [103].
The climatic data correspond to mean values obtained from a long series of annual observations.Only those sample plots of age 5 years or more were therefore considered in the model development to prevent including growth data that could be affected by anomalous climatic conditions occurring during short periods of time.Finally, data from 132 sample plots of gorse-dominated communities and from 62 sample plots of hearthdominated communities were used to develop the MARS models.

Development of Total Fuel Load Growth Curves
The parameter estimates for each base model and the goodness-of-fit statistics are given in Table 3.All parameters were found to be significant at the 5% level.The values of the goodness-of-fit statistics are very similar for the three base models in the case of both gorse-dominated and high heath-dominated communities; although in the latter case, the values indicate more accurate estimates.In both cases, the base model with the highest percentage of variability explained and the lowest RMSE was the Korf model, which was therefore selected for generating the W Shr -t growth curves in both communities.Three different families of curves can be obtained from the guide curve derived from the Korf model, depending on which parameter is considered site-dependent (a 0 , a 1 or a 2 , with the other two held constant).The three families of curves for each community overlaid on the observed W Shr -t data are shown in Figure 2. The curves represented correspond to W Shr values of 1.5, 3, 4.5, 6 and 7.5 kg m −2 at a reference age (t ref ) of 10 years for gorsedominated communities and 1, 2.5, 4, 5.5 and 7 kg m −2 at a reference age (t ref ) of 10 years for high heath-dominated communities.
The best family of curves was selected according to the performance of the site form curves over the data by considering the polymorphism of the resulting curves as a fundamental criterion.The family in which parameter a 0 varies when considered site-dependent results in anamorphic curves with different asymptotes (Figure 2, left); when parameter a 1 is considered site-dependent, the family of curves is polymorphic with a common asymptote (Figure 2, middle), and the same occurs when parameter a 2 is considered site-dependent (Figure 2, right).a1 is considered site-dependent, the family of curves is polymorphic with a common asymptote (Figure 2, middle), and the same occurs when parameter a2 is considered sitedependent (Figure 2, right).

Gorse-dominated communities
High heath-dominated communities Considering the trends in the different families of curves in Figure 2, the most representative for both shrubland communities is that in which parameter a1 depends on the site, resulting in a polymorphic family.The mathematical expression of the curve for a sample plot, in which the total shrub fuel load (WShr_t1) at a given age (t1) is known, would thus be obtained as follows: In addition, the site index (SI) of that same sample plot, defined as the value of the total shrub fuel load at the reference age (tref) of 10 years, would be obtained as follows: Considering the trends in the different families of curves in Figure 2, the most representative for both shrubland communities is that in which parameter a 1 depends on the site, resulting in a polymorphic family.The mathematical expression of the curve for a sample plot, in which the total shrub fuel load (W Shr_t1 ) at a given age (t 1 ) is known, would thus be obtained as follows: In addition, the site index (SI) of that same sample plot, defined as the value of the total shrub fuel load at the reference age (t ref ) of 10 years, would be obtained as follows: A reference age (t ref ) of 10 years was selected to improve the accuracy of predictions, reducing the prediction bias associated with sample plots, where t differs greatly from the reference t ref [104].
Site index values for the gorse-dominated formations ranged from 0.84 to 8.23 kg m −2 at a reference age of 10 years, with a mean value of 3.79 kg m −2 (sd = 1.15 kg m −2 ), while for the high heath-dominated formations, the site index ranged from 1.01 to 8.25 kg m −2 with a mean of 3.68 kg m −2 (sd = 1.61 kg m −2 ).

Effect of Topographic and Climatic Variables on Site Quality
Mean values of topographic and climatic variables considered for the sampled shrubland communities dominated, respectively, by U. europaeus and E. australis in the study sites (Table 2) were significantly different, except for monthly number of rainy days < 1 mm for summer.In general, the ranges of variables were similar, except for elevation, with high heath-dominated communities having a higher elevational limit than gorse-dominated communities, while the opposite was found for Tcm.Additionally, the aridity index P/T was larger in heath than in gorse.
Preliminary analysis indicated no significant linear correlations between the physiographic and climatic variables and the site index for both communities under study.However, the results of the MARS model for gorse-dominated communities showed a significant effect of both mean summer precipitation (Ps) and the mean temperature of the coldest month (Tcm): where SI is expressed in kg m −2 , Ps in mm and Tcm in • C; I Ps is a dummy variable with a value of 1 for Ps greater than 167 mm and 0 otherwise; I Tcm1 is a dummy variable with a value of 1 for Tcm greater than 8.5 • C and 0 otherwise, and I Tcm2 is a dummy variable with value 1 for Tcm less than 8.5 • C and 0 otherwise.For equal values of Tcm, the threshold for summer precipitation was 167 mm, below which SI values are similar and above which SI increases at an average rate of 0.84 kg m −2 for every 10 mm increase in precipitation.In the case of Tcm, for equal Ps values, an increase in its value implies a moderate increase in SI until a threshold of 8.5 • C, beyond which SI growth increases considerably.These thresholds and growth rates should be considered with caution as the predictive capacity of the model is limited (ME = 0.11; RMSE% = 25.58%),although the model is useful from an explanatory point of view.

Compatible System for the Shrub Fuel Complex
Once the SI of each sample plot was estimated using Equation ( 11), the compatible system of four equations relating age and site index to the biomass of the remaining fractions of the shrub fuel complex (W Shr_G1_dead , W Shr_G1_live , W Shr_G1 and W Shr_G23 ) was fitted for each shrub community.The mathematical expression of the system (Equations ( 6) to ( 9)) fitted to gorse-dominated and high heath-dominated shrub communities and the goodnessof-fit statistics are given in Table 4.All parameters were found to be significant at the 5% level.
The results of the White test and the values of the condition number did not indicate any problems of heteroscedasticity or multicollinearity, respectively, for any fuel fraction in either of the shrub communities.
The best results were obtained for the high heath-dominated community, with percentages of observed variability explained varying between 73% for dead fine shrub load (W Shr_G1_dead ) and 94% for coarse shrub load (W Shr_G23 ); for the gorse-dominated community, the respective percentages varied between 48% of live fine shrub load (W Shr_G1_live ) and 86% of coarse shrub load (W Shr_G23 ).
Vega et al. [7] also obtained better goodness-of-fit statistics for the high heath-dominated community versus the gorse-dominated community in the same study area by fitting a system of compatible equations to estimate shrub loads by fractions using the mean height of the shrub, shrub cover and the depth of the litter and duff layer as independent variables.Table 4. Mathematical expressions and goodness-of-fit statistics of the system of four equations fitted to estimate the biomass of the shrub fuel complex in gorse-and high heath-dominated communities using site index and age as independent variables.0.9377 0.2792

Gorse-Dominated
The complete set of fitted equations (Equations ( 6) to ( 11)) can also be used to estimate other variables of interest such as the fine ratio (W Shr_G1 /W Shr ) and the dead ratio (W Shr_G1_dead /W Shr_G1 ) with the age and the site index of the shrub community or the mean and current fuel load annual increment of each of the five fractions (W Shr_G1_dead , W Shr_G1_live , W Shr_G1 , W Shr_G23 and W Shr ).The mathematical expression of the fine and dead ratios would be as follows: As parameters b 1 and b 2 are both positive for both shrubland communities under study (Table 3), the fine ratio decreases with age and also with increasing site quality (increasing SI).This behavior is biologically expected, as the fine shrub load increases in absolute terms with age and site quality, but does so at a slower rate than the coarse shrub load (shrub fraction with no upper limit of thickness).On the other hand, as c 2 is negative and c 1 is zero for both shrub communities analyzed (Table 3), the dead ratio varies exclusively with the age of the community, increasing as the latter increases.In our study, the fine ratio varied in gorse-dominated communities between 0.292 and 0.989 (mean value = 0.673), and in high heath-dominated communities, the range was similar (0.266-0.988) with a mean value slightly higher (0.727).
The mean annual increment (MAI, kg m −2 year −1 ) in the shrub load of any of the fractions analyzed can be obtained by dividing the expression of the corresponding fitted equation by age; likewise, the current annual increment (CAI, kg m −2 year −1 ) can be obtained by calculating the derivative of the same equation relative to age.For example, the MAI and CAI of the total shrub load (W Shr ) can be obtained as follows: Fire 2024, 7, 126 13 of 25 Total fuel loading is an important component in operational fire hazard systems [105,106]; however, fine fuel loading, being the amount of fuel consumed in the flaming phase of combustion, has a primary influence on fireline intensity [23,33,107].Within this fraction, the load and proportion of fine dead fuel, combined with its low moisture content, favor flammability in heaths [108] and gorse [61], fire initiation [66,109], fire spread and duration [23,109,110].For these reasons, Figure 3 shows the W Shr_G1 -age curves and the W Shr_G1_dead -age curves for five different site indices for both shrubland communities analyzed (1.5, 3, 4.5, 6 and 7.5 kg m −2 of W Shr at a reference age of 10 years for gorse-dominated communities and 1, 2.5, 4, 5.5 and 7 kg m −2 of W Shr at a reference age of 10 years for high heath-dominated communities).
equation by age; likewise, the current annual increment (CAI, kg m −2 year −1 ) can be obtained by calculating the derivative of the same equation relative to age.For example, the MAI and CAI of the total shrub load (WShr) can be obtained as follows: Total fuel loading is an important component in operational fire hazard systems [105,106]; however, fine fuel loading, being the amount of fuel consumed in the flaming phase of combustion, has a primary influence on fireline intensity [23,33,107].Within this fraction, the load and proportion of fine dead fuel, combined with its low moisture content, favor flammability in heaths [108] and gorse [61], fire initiation [66,109], fire spread and duration [23,109,110].For these reasons, Figure 3 shows the WShr_G1-age curves and the WShr_G1_dead-age curves for five different site indices for both shrubland communities analyzed (1.5, 3, 4.5, 6 and 7.5 kg m −2 of WShr at a reference age of 10 years for gorse-dominated communities and 1, 2.5, 4, 5.5 and 7 kg m −2 of WShr at a reference age of 10 years for high heath-dominated communities).The pattern of variation of fine fuel load with age is very similar in both shrubland communities, with a first phase of increase until reaching a maximum and a subsequent decline phase less pronounced.The maximum value of WShr_G1 is reached earlier in the best site qualities for both communities but is more premature in high heath-dominated communities (4 years for a SI of 7 kg m −2 ) than in gorse-dominated communities (6 years for a The pattern of variation of fine fuel load with age is very similar in both shrubland communities, with a first phase of increase until reaching a maximum and a subsequent decline phase less pronounced.The maximum value of W Shr_G1 is reached earlier in the best site qualities for both communities but is more premature in high heath-dominated communities (4 years for a SI of 7 kg m −2 ) than in gorse-dominated communities (6 years for a SI of 7.5 kg m −2 ).Thereafter, the decrease is more pronounced in the best site qualities, and the curves are even cut at advanced ages (more than 20 and 16 years for gorse-and high heath-dominated communities, respectively).
Regarding the W Shr_G1_dead curves, the pattern is also similar in both shrubland communities, with an initial growth to a higher maximum in the better site qualities and a posterior decrease.The maximum is reached earlier in the best site qualities and, again, is more premature in high heath-dominated communities (5 years for an SI of 7 kg m −2 in these communities vs. 9 years for an SI of 7.5 kg m −2 in gorse-dominated communities).It is important to emphasize that these W Shr_G1_dead maxima are delayed with respect to the W Shr_G1 maxima, and these delays are longer as site quality becomes worse (from 3 years for an SI of 7.5 kg m −2 to 8 years for an SI of 3 kg m −2 in the gorse-dominated communities and from 1 year for an SI of 7 kg m −2 to 3 years for an SI of 2.5 kg m −2 for the high heath-dominated communities).

Discussion
Although remote sensing methods provide invaluable information on biomass accumulation and growth in shrublands, detailed field inventories continue to be critical to obtain information on shrub structural variables at stand level and at a fine scale.Each method provides complementary information, and while field data are necessary to obtain accurate fuel data and test and refine remote sensing models, remote sensing methods represent the only feasible approach for geospatial assessment of these shrub attributes.
To our knowledge, the models developed in this study are the first to enable detailed estimation of the dynamics of total and different fuel load fractions in gorse and high heath communities at stand level.The few scientific examples of dynamic models that disaggregate shrubland fuel loads by fractions refer to the individual plant level [111,112].However, in multi-stemmed, dense shrub stands in temperate ecosystems, individual-based biomass equations are impractical because measurements of plant density and individual plant attributes are not feasible [7].Furthermore, our approach, using a site index based on biomass levels at a reference age of 10 years, addresses the issue of the typical wide variability in shrubland biomass, with a reasonably good level of accuracy.

Effect of Topographic and Climatic Variables on Site Quality
Elevation and climatic data for gorse (Table 2) are consistent with the biogeographic information classifying U. europaeus as a maritime climate-influenced [113][114][115] and mesotemperate species, preferably thriving in the colline and montane elevational levels [116].Compared to its current global climatic niche [117], the gorse stands included in this study receive abundant annual rainfall, while the annual mean temperature is close to the global mean value.Gorse plants in NW Spain also have an advantageous climatic position for growth, relative to those growing in the Western European habitat [118].
The significant differences between gorse-and high heath-dominated communities in almost all climatic parameters considered (Table 2) are consistent with the distinct habitats.The high heath-dominated communities grow preferentially in the eastern part of the region, further from the oceanic effect, more under the influence of the Mediterranean climate, and at a higher elevation than the gorse-dominated communities.Although the data were partly derived from mixed stands with the presence of other typical species of dry heathlands, they confirm the known biogeographic differences between high heath and gorse [115,119].
The detection of a summer rainfall threshold affecting SI suggests, however, that gorse growth is responsive to an increase in water input in summer, irrespective of the high average annual precipitation rate.This also appears to be consistent with the water shortage for the gorse area in summer in NW Spain (Table 2).Summer rainfall is important for gorse biomass dynamics, as part of the growth occurs in summer [57], coinciding with increased evapotranspiration.Gorse has a tap root that can effectively explore subsurface soil water [113] and can withstand some summer droughts (if not very intense and widespread) [114,117,120].The resistance, however, will probably lead to a reduction in metabolism, as found in Ulex airensis [121], with growth ultimately affected.Furthermore, gorse mainly grows in mountainous areas in NW Spain, where shallow, coarse-textured soils are most abundant, generally resulting in low water storage and availability in summer [122,123].On the other hand, the rainfall pattern may be more important in explaining the impact of drought on growth than the amount of rainfall [124], and the mean monthly number of summer rainy days in the study area is comparatively low (4.8).Reduced growth under drought has been reported for U. europaeus [125] and for U. parviflorus [126], and high drought sensitivity of U. europaeus has been noted at the seedling stage [127].
The Tcm threshold seems compatible with the gorse sensitivity to low temperature [113,117,118,128].Overall, the plasticity of gorse [117,118,129] may also partly explain the modest effect of climatic variables on SI observed in the present study.
The lack of any apparent relationships between SI and climatic variables in heathdominated communities is not surprising, given the inconsistent information regarding the relationship between water availability and resprouting growth after disturbance in these species [130][131][132][133] as well as in other Mediterranean high heath species such as E. multiflora [134].The inconsistent response to moderate drought and its well-established strong resilience to various types of disturbance [135][136][137][138][139][140], although sensitive to very frequent clipping [138,141], suggest lower sensitivity to climate than in gorse.The lack of growth response to rainfall and temperature detected in the high heath-dominated communities may also be partly due to the favorable mean annual rainfall-mean annual temperature ratio (150 ± 4.4 mm) for the high heath habitat in NW Spain.Meanwhile, high heath-dominated communities may be sensitive to soil characteristics [142], which also affect resprouting [130], as found for different levels of soil burn severity [143].
Despite the reduced response of biomass growth to climatic variables in gorse-and high heath-dominated communities, the significant increase in mean annual temperature observed in the last few decades in NW Spain and the absence of any changes in precipitation [144] indicate reduced soil water availability.This situation, together with the observed increase in the duff moisture code (strongly and negatively related to the moisture content of the soil organic layer and topsoil) [145], suggests probable adverse changes in plant growth conditions, also leading to increased fire severity in shrubland wildfire.
This study did not address the influence of other factors on the biomass dynamics of gorse and high heath communities.Such factors, along with the study limitations, may also have partly constrained the detection of a stronger influence of climate on growth.One possible limitation of the study is the representativeness of climatic data in relation to the sampled sites and the temporal range of growth considered.Although gorse and high heaths are the dominant species in their respective communities, the spatial composition varied widely and may have affected the biomass growth patterns, increasing the variability.Finally, many other factors, such as topographic characteristics, soil properties, plant competition level and the plant physiological status, can strongly affect temporal patterns of plant growth.The negative consequences of the NW Spain disturbance regime, with an active millennia-long human presence in the region [146,147] and a long history of intensive use of shrub biomass, have probably affected biomass dynamics.Disturbances due to frequent burning, wildfire, grazing and rotational cultivation in prehistoric [148] and historic times [149,150] have probably also influenced the dynamics.Although different types of disturbances, particularly fire, are considered key factors in the perpetuation of heathlands [151,152], some of the disturbances have been recurrent, frequent and severe, dramatically reducing heathland area in the region, and may have had long-lasting effects on plant growth patterns [153].All of the above may have had an overwhelming influence Fire 2024, 7, 126 on both gorse-and high heath-dominated communities, which may have altered the climatic effect on the growth of both shrublands.Further research is required to shed light on this point in order to refine current models.

Biomass Accumulation and Mean Annual Increment (MAI)
The similarity of the values of biomass mean annual increment (MAI) in gorsedominated communities in the sample plots of this study with those reported for native Western Europe stands [154] seems consistent with similar climatic conditions and history of land use and disturbance in those habitats [155].However, data reported by Carswell et al. [49] for this species in treeless stands in New Zealand show high growth maintained even beyond 15 years, which could be related to a niche shift, or simply due to the short time interval explored.Since Oceanic climatic conditions are similar to the native gorse region, the response observed in New Zealand does not seem caused by a climatic niche shift [117] and may have been due to a more complex environmental niche shift.This can entail less competition and fewer natural enemies than in the native niche [129,156], enhancing resource assignation to reproduction and growth (EICA hypothesis: [157,158]).However, the effect of differences in soil properties between gorse from native areas and from sites measured in New Zealand cannot be ruled out in regard to explaining differences in the respective MAIs.
The MAI values of the Mediterranean gorse (Ulex parviflorus)-dominated communities at ages 9 and 17 years [12] were only slightly lower than those from our guide curve for the same age.This was not expected, since the two species have different strategies of regeneration; while Mediterranean gorse is an obligate seeder [12,31,159], U. europaeus is a facultative resprouter [160], which a priori provides an advantage for faster growth.However, at the individual plant level, U. parviflorus biomass between 5 and 26 years showed a pattern of continued growth [161].
The guide curve for the high heath-dominated communities in our study shows a lower MAI at juvenile stage (5 years of age) than that reported by Fernandes et al. [47] and similar to those found by Fernandez-Abascal et al. [162] for the same species and by Montero et al. [51] in E. arborea in sites under Mediterranean climate; although, they surpass the value of these other Mediterranean communities from 10 years of age onwards.
The pattern of biomass accumulation with age in the two shrubland communities considered in this study seems compatible with the traditional summer agricultural burning by a smoldering fire that produces ash, which is used as fertilizer for a subsequent rye crop.The rotation cycles of 27, 36 and 41 years generally applied to this practice, depending on the site quality [163], ensure the accumulation of sufficient fuel to cause severe burns and the recovery of favorable edaphic properties for cultivation.

Fine Ratio
Information on the dynamics of fine fuel loading in U. europaeus was lacking, until now.Our data suggest that, in addition to age, the structure of oceanic gorse biomass is quite sensitive to local conditions, as reflected in the influence of SI on fine fuel dynamics.Through SI, the findings also demonstrate a weak but appreciable influence of climate on gorse fine fuel ratio, with an increase in Tcm and Ps, both leading to an increase in total standing shrub biomass but a decrease in the fine fuel ratio.The findings also indicate that with a single fuel inventory and determination of stand age, it is possible to substantially improve the accuracy in fine fuel ratio estimation over time.Baeza et al. [12,31] examined the variation in the fine fuel load ratio with age in U. parviflorus in E. Spain, reporting mean values of 0.95, 0.71 and 0.48 for 3-, 9-and 12-year-old plants, respectively.These values are generally higher than those observed for U. europaeus in our study.
The dynamics of the fine fuel ratio in the high heath-dominated community is very similar to that in gorse-dominated communities, with slightly lower values for the same SI and age.For dry low heaths, with Pterospartum tridentatum and E. umbellata as dominant species, Fernandes and Rego [55] reported a fine fuel (<2.5 mm) ratio decreasing from 0.9 at 2 years to 0.70 at 22 years in N. Portugal, the latter value being much higher than that observed for the guide curve in our case at the same age (0.38).

Dead Ratio
The apparent independence of the dynamic of this parameter from SI in both gorseand heath-dominated communities was unexpected and may indirectly indicate that this plant trait has an evolutionary character with a genetic basis [164] that increases plant flammability [25,165].The dead ratio increases rapidly in gorse-dominated stands in the early years of age, with these communities becoming highly flammable [59,60].This is consistent with the widespread historical use of burning in cycles as short as six years in the region [149,150], given the importance of the amount of fine dead fuel and moisture for fire initiation and spread in gorse [66,109] and also in other shrublands [166].
Overall, the increase in the dead ratio with age was quite low and with a similar age-related variation in gorse-and high heath-dominated communities.In gorse, the mean dead ratio for the guide curve increased from 0.32 to 0.45 at between 2 and 22 years of age, whereas in high heath-dominated communities, it increased from 0.22 to 0.32.The dead ratio varied at a greater rate at the juvenile stage, up to about five years of age, and then scarcely increased up to the mature stage (≈10 years).In contrast, Baeza et al. [12,167] observed greater age-related variation in the dead ratio in Ulex parviflorus, reaching a maximum value (0.56) at 9 years of age, followed by a decrease to 0.40 at 26 years of age, similar to the value obtained for the guide curve at the same age in this study (0.46).However, the dead ratio in the Mediterranean tall heath Erica multiflora increased continuously from 0 to 0.15 [167], in contrast to the present observations for E australis.Fernandes et al. [56] also observed a linear increase in the dead ratio (<2.5 mm thickness material) with age in a dry low heathland of age between 2 and 22 years in Portugal, leading to a significant change in the estimated fire behavior.

Implications for Fuel Hazard Management
The site-specific equations presented in this study for estimating the future development of different shrubland formations at stand level can be useful in different stages of operational fuel hazard management.Only the value of W Shr and the age of a sample in the target community are required to improve substantially the accuracy in the estimation of different parameters related to shrubland fuel load and predict their change with age.Our results also offer guidance for fuel treatment planning, based on the age at which different fuel fractions reach the maximum fuel load.They suggest conducting fuel reduction treatments in relatively young stands and, due to the fast recovery to high fuel loads, repeating the treatments frequently.However, fuel managers should balance information on the rate of fuel accumulation with the local risk of ignitions to optimize the interval between interventions.

Conclusions
The models developed in this study provide new information on biomass and fuel load dynamics in the two main shrubland communities in NW Spain, dominated by gorse and high heaths.These communities cover a large area in the region, and they are involved in frequent wildfires and play an important role in terms of biodiversity, C sequestration and potential energetic use, among other functions.Despite the extensive area of gorse at the global level, paradoxically, a review of the relevant literature revealed that knowledge about fuel load dynamics in this species is extremely scarce.Such knowledge is even more limited for E. australis.
Native gorse stands in NW Spain share a range of total fuel load dynamics that are similar to those of stands in geographically very distant areas, in some of which gorse is an introduced species.This study contributes, for the first time, to analyzing the influence of climatic factors on fuel load dynamics in gorse habitats.The results suggest that although gorse-dominated communities receive a fairly high amount of rainfall in the study area, growth showed slight sensitivity to increased summer rainfall, while the preference of gorse for mild temperatures was confirmed.However, the growth of high heath-dominated communities was not sensitive to any of the climate variables explored, in line with previous studies giving rise to inconclusive results regarding the effects of drought on growth.However, in the framework of global change, future studies should further explore the influence of climatic factors on biomass and fuel dynamics in the two communities under consideration.
Our model enables, for the first time, the estimation of biomass dynamics in two different fractions (<6 mm and >6 mm) in gorse and high heath-dominated shrublands and for two vegetative stages (dead and alive), with reasonably good accuracy.Given the generally high variability in fuel loads in these shrub communities, the available models so far provided estimates with a large range of uncertainty.The approach used addresses this by modeling the dynamics of fuel biomass and fractions as a function of site index (SI).In this approach, only the value of shrub fuel load and the age of a sample in the target community are required.Furthermore, the equations developed ensure the fulfillment of the additivity principle, which has not been addressed in other studies of dynamic fuel loading in shrublands.
Although this study focused on fuel load dynamics in the shrub layer, the litter layer can also actively contribute to fire behavior in shrublands.Furthermore, other fuel structural properties, such as fuel height and bulk density of shrub and litter layers, are important for fire behavior, and the dynamics were not addressed here.As in other research on biomass dynamics, our study was limited by the lack of quantification of the disturbance history influence on fuel accumulation, which has probably affected the vegetation development in the study area and may contribute to the high variability in the fuel loads in both gorse-and heath-dominated communities.The influence of soil on growth was also not considered here.The effects of all of these factors remain as challenges for future research.
On the other hand, our findings indicate that the dead fuel ratio in both gorse-and high heath-dominated communities depends on age alone, with no apparent influence of site quality, which facilitates the estimation of this important variable in the two communities examined.We did not find a clear steady state for W Shr ; therefore, future studies may increase the effort in older stands to better clarify the longer-term fuel load dynamics.
The information obtained in our study may be useful for fire management planning.Thus, fuel hazard assessment and monitoring, updating of fuel maps, fuel reduction treatment planning and implementation as well as fire behavior prediction in wildfire suppression can benefit from our findings.Furthermore, the study findings can improve current knowledge on C cycling and C sequestration dynamics in shrublands, assessment of their bioenergetic potential and the dynamics of competition between species after a disturbance that may favor the expansion of invasive species, among other aspects.Finally, the approach used, which combines intensive destructive field sampling, detailed laboratory work to classify the fuel load into components and the construction of accurate empirical models, can be applied to other dense thickets with high canopies and may be useful in combination with remote sensing approaches.

Figure 1 .
Figure 1.Geographical location of the 316 inventory plots in Galicia encompassing two shrub communities (left).Layout of inventory plot showing transects and location of four destructive sampling subplots (quadrats in green) in each inventory plot (right).

Figure 1 .
Figure 1.Geographical location of the 316 inventory plots in Galicia encompassing two shrub communities (left).Layout of inventory plot showing transects and location of four destructive sampling subplots (quadrats in green) in each inventory plot (right).

Figure 2 .
Figure 2. Three different families of WShr-t curves generated from the guide curve (red line) of the base model proposed by Korf (cited in Ludqvist [90]) for five different site indexes (1.5, 3, 4.5, 6 and 7.5 kg m −2 at a reference age of 10 years for gorse-dominated communities and 1, 2.5, 4, 5.5 and 7 kg m −2 at a reference age of 10 years for high heath-dominated communities).

Figure 2 .
Figure 2. Three different families of W Shr -t curves generated from the guide curve (red line) of the base model proposed by Korf (cited in Ludqvist [90]) for five different site indexes (1.5, 3, 4.5, 6 and 7.5 kg m −2 at a reference age of 10 years for gorse-dominated communities and 1, 2.5, 4, 5.5 and 7 kg m −2 at a reference age of 10 years for high heath-dominated communities).

Fire 2024, 7 , 25 Figure 3 .
Figure 3. WShr_G1-age curves (left column) and WShr_G1_dead-age curves (right column) for the SI corresponding to the guide curve (red) and five different site indexes: 1.5, 3, 4.5, 6 and 7.5 kg m −2 of WShr at a reference age of 10 years for gorse-dominated communities (upper row) and 1, 2.5, 4, 5.5 and 7 kg m −2 of WShr at a reference age of 10 years for high heath-dominated communities (lower row).

Figure 3 .
Figure 3. W Shr_G1 -age curves (left column) and W Shr_G1_dead -age curves (right column) for the SI corresponding to the guide curve (red) and five different site indexes: 1.5, 3, 4.5, 6 and 7.5 kg m −2 of W Shr at a reference age of 10 years for gorse-dominated communities (upper row) and 1, 2.5, 4, 5.5 and 7 kg m −2 of W Shr at a reference age of 10 years for high heath-dominated communities (lower row).

Table 2 .
Summary statistics of the physiographical and climate variables extracted from the raster datasets corresponding to the 267 sample plot locations.Standard deviations are shown in parentheses.Different letters indicate significant differences between mean values (α = 5%).

Table 3 .
Parameter estimates and goodness-of-fit statistics for the three base models fitted to W Shr -t data from gorse-and high heath-dominated communities.