Estimating Energy Concentrations in Wooded Pastures of NW Spain Using Empirical Models That Relate Observed Metabolizable Energy to Measured Nutritional Attributes

: Wooded pastures serve as a traditional source of forage in Europe, where forest grazing is valued as an efﬁcient tool for maintaining the diversity of semi-natural habitats. In a forest grazing setting with diverse diet composition, assessing the energy content of animal diets can be a difﬁcult task because of its dependency on digestibility measures. In the present study, prediction equations of metabolizable energy (ME) were obtained performing stepwise regression with data ( n = 297; 44 plant species) on nutritional attributes (Acid Detergent Fiber, lignin, silica, dry matter, crude protein, in vitro organic matter digestibility) from 20 representative stands of Atlantic dry heathlands and pedunculate oak woodlands. The results showed that the prediction accuracy of ME is reduced when the general model (R 2 = 0.64) is applied, as opposed to the use of the speciﬁc prediction equations for each vegetation type (R 2 = 0.61, 0.66, 0.71 for oak woodlands; R 2 = 0.70 heather-gorse dominated heathlands, R 2 = 0.41 continental heathlands). The general model tends to overestimate the ME concentrations in heaths with respect to the observed ME values obtained from IVOMD as a sole predictor, and this divergence could be corrected by applying the speciﬁc prediction equations obtained for each vegetation type. Although the use of prediction equations by season would improve accuracy in the case of a Winter scenario, using the general model as opposed to the prediction equations for Spring, Summer or Fall would represent a much smaller loss of accuracy. results showed that the prediction accuracy of the general model applied, opposed to the speciﬁc prediction the case of heathlands. The general model


Introduction
Grazing is one of the influential and sustainable management interventions for shaping and maintaining semi-natural habitats and promoting greater biodiversity and multifunctionality in forest understory and heathlands [1][2][3][4]. Horse grazing, for example, can decrease gorse dominance and promote heath composition in Atlantic heathlands [4][5][6], which are prized for their biodiversity, aesthetic and cultural value, and supported by the European Union (EU) within the High Nature Value farming framework. Grazing oak woodlands is included as a type of silvopasture that has been practiced throughout the ages across Europe and has gained growing recognition in the last decades [7][8][9][10][11][12]. Several studies have also pointed out the importance of livestock grazing as a management strategy that can improve rural sustainability while effectively controlling the accumulation of flammable woody vegetation derived from rural depopulation and the abandonment of management practices in the heathlands of NW Spain [5,6,13,14], thus contributing to the mitigation of greenhouse gas emissions by reducing the risk of fire [15]. Furthermore, in projected climate change scenarios, a stable forage supply in hot, dry years would be The aim of the current study was to obtain prediction equations of ME using a wide range of data on nutritional attributes (IVOMD, ADF, cellulose, lignin, silica, dry matter and crude protein) from various stands that represent prototypes of European Atlantic heathlands and pedunculate oak woodlands. We hypothesized that if the existing correlations between the measured nutritional characteristics and observed ME could result in obtaining significant empirical linear models for such vegetation types, then such prediction equations could be broadly applicable in stands with similar vegetation depending on their available data on chemical composition, thus facilitating grazing management decisions that are based on the ME requirements of the animals.

Nutritional Characterization of Vegetation
For the nutritional analysis, 44 plant species were selected according to their abundance in each forest site and their palatability for ruminants (e.g., roe deer, goat, beef cattle, sheep and horse). Their occurrence by vegetation type was as follows: 20 such species were present in O-Bl, 17 in O-My, 12 in O-Ru, 16 in Atlantic heathlands and 5 in continental heathlands. Plant samples for each species were composed of apical portions no more than 15 cm in length and 0.5 cm in diameter, and were clipped randomly from selected individuals before being combined into a single sample for analysis. Samples were collected quarterly during 2 years for most of the species.
For nutritional characterization, acid detergent fiber (ADF), permanganate lignin, cellulose and silica contents were obtained with sequential detergent analysis [38]. Crude protein was determined by analyzing the nitrogen content using the micro-Kjeldahl method, with the results being multiplied by 6.25. In vitro organic matter digestibility (IVOMD) was assessed using the technique of Tilley and Terry (1963) as modified by Alexander (1969) [39,40]. More details about the methods of sampling and analysis can be found in González-Hernández and Silva-Pando (1999) [17]. Additionally, the IVOMD values were used to estimate observed ME using the equation reported by Barber et al. (1984) and recommended by AFRC (1993) [28]: ME (MJ/kg DM) = 0.0157 × IVOMD (g/kg DM). ME concentrations obtained from applying this equation will be referred to in the text as observed ME concentrations or ME equations based on IVOMD as a sole predictor.
With the purpose of being a reference for the nature and range of the data from which the ME equations were obtained, the nutritional characteristics by plant species are summarized in Appendix A. We used the coefficient of variation (CV) to refer to the relationship between the size of the mean and the variability of each nutritional parameter. Values of observed ME for the same plant species were grouped into those cases where such species occurred in more than one type of oak woodlands, or was in the two types of heathlands, and when certain homogeneity in their ME values was observed (i.e., low CV). In the field of agricultural sciences, CV is broadly accepted when it is ≤20%. A scale of reference is often as follows: ≤10 is very good, 10-20 is good, 20-30 is acceptable and CV ≥30 is not acceptable. Because of the limited records of some forbs and grasses in our data set, we decided to pool them together.

Data Analysis to Obtain Prediction Equations of ME
To analyze the possible linear relationships between ME and nutritional characteristics, Pearson correlation coefficients were calculated and the trends of pairs of ME data with each possible dependent variable were visualized. In view of these trends, it was decided to also include the logarithmic transformation of nutritional characteristics as possible dependent variables. Cellulose was excluded as an independent variable because it is a parameter rarely consider in the nutritional reports from laboratory analysis. It is also included as a fraction of the parameter ADF, and its occurrence as a good predictor in attempts of ME prediction equations is not common.
Empirical linear models that relate the observed ME to the measured nutritional characteristics were fitted for each one of the 5 vegetation types for each season, and for the pooled data using the SAS/STAT ® software [41]. Construction of the models involved the following two steps: (i) selection of the best set of independent variables and (ii) evaluation and cross-validation of the fitted models by graphical and numerical analysis of the residuals.
To select the best set of independent variables, the stepwise variable selection method [42], combined with the values of the Pearson correlation coefficients and an understanding of the fitting process, was used. Multicollinearity is a problem that could be associated with the proposed models that violate the fundamental least squares assumption of independence and equal distribution of errors with zero mean and constant variance. Multicollinearity refers to the existence of high intercorrelations among the independent variables in multiple linear or nonlinear regression analysis, because some of the variables represent or measure similar phenomena. Although the least squares estimates of regression coefficients remain unbiased and consistent under the presence of multicollinearity, this problem may seriously affect the standard errors of the coefficients, invalidating confidence intervals [43]. Thus, models with low multicollinearity should be selected whenever possible [44]. To evaluate the presence of multicollinearity among variables in the models analyzed, the variance inflation factors (VIF) of all the independent variables were calculated, and values up to 5 were accepted [45]: i is the multiple determination coefficient obtained when the ith independent variable X i is regressed against all of the remaining independent variables in the model.

Model Evaluation and Validation
Four statistical criteria obtained from the residuals (e i ) were examined to test the accuracy and precision of ME estimates: bias (E), mean absolute error (|E|), root mean square error (RMSE) and the adjusted coefficient of determination (R 2 adj ). These expressions can by summarized as follows: where e i = y i − y i and y i ,ŷ i and y are the measured, predicted and average values of the dependent variables, respectively; n is the total number of observations used to fit the model and p is the number of model parameters.
Because the quality of the fit does not necessarily reflect the quality of the prediction, cross-validation of each model was carried out, and the residual for one observation was estimated by fitting the model without that observation. MEF, equivalent to the R 2 of the fitting phase, was then calculated from the residuals from the cross-validation.
Finally, the estimates obtained with the general model (fitted to all data) when applied to each of the five vegetation types and to the different stations are compared with the estimates obtained with the models fitted specifically for those situations in order to estimate the error of using the general model when it is not possible to discriminate in the database between vegetation types or seasons.

Nutritional Characteristics, Digestibility and ME (Based on IVDOM as Sole Predictor)
The chemical composition and digestibility of the vegetation types and plants studied here are summarized in Appendix A. A high variation in observed ME concentrations in continental and Atlantic heathlands was found (CV = 35.9 and 45.6%, respectively), suggesting a large range in the observed ME concentration for the main plants included in these vegetation types.
Observed ME values for continental heathlands ranged between 1.7-7.2 MJ/kg DM, whereas ME concentrations of plants in Atlantic heathlands extended to a wider 1.8-12.2 MJ/kg DM interval. The lowest ME concentrations in heathlands were found in heaths, followed by Halimium spp. and leguminous shrubs. On the other hand, the variation in ME values in understory plants of O-Bl and O-Ru oak woodlands was smaller than in heathlands (CV 28.7, 32.8%), and ranged from 2.2-12.6 and 3.6-11.2 MJ/kg DM, respectively. Only the degree of variation in ME values in O-My (CV = 38.3%; 2.8-13.6 MJ/kg DM) exceeded the one found for continental heathlands. ME concentrations of oak woodland understory is expected to be higher than heathlands because of the superior ME contribution from plants like Hedera helix, Ilex aquifolium, Lonicera periclymenum, leguminous brooms and Frangula alnus. Forbs and grasses, when present, would contribute considerably to increasing the ME concentrations in all vegetation types (see Tables A1-A4 in Appendix A). The variation in ADF, lignin and cellulose contents for each vegetation type was smaller than the ME oscillations, especially for heathlands (CV between 9.6-17%). The CP contents and CP variation found in continental heathlands (50.9 and 86.8 g/kg DM, CV = 14%) were more minor than in oak woodlands and Atlantic heathlands (CP minimum values around 50 g/kg DM and maximum around 160-205 g/kg DM, CV between 27-30%).

Relationships between Observed Metabolizable Energy and Nutritional Parameters
Pearson correlation coefficients were conducted to evaluate the relationships between nutritional composition and observed ME values ( Figure 2). The correlation matrices for each vegetation type revealed that there was a negative and highly significant relationship between ME and total fiber (ADF) or its fraction cellulose (r values ranging mostly from −0.50 to −0.67; p < 0.001). Conversely, no significant correlation was found between ME and lignin. The relationship of silica and ME was also significant, though it was much weaker in general (r = −0.24 to −0.41; p < 0.01; p < 0.05). Dry matter contents of plants were also highly significant and negatively correlated with ME (r values varying from −0.48 to −0.70; p < 0.001). CP, the only parameter positively correlated with ME, exhibited a weaker relationship in oak woodlands, especially in O-RU (r = 0.30; p < 0.05), although this relationship was highly significant in Atlantic heathlands (r = 0.77, p < 0.001).
Overall, the correlation coefficients analyses performed for each vegetation type showed similar relationships between ME and nutritional parameters, except for continental heathlands. For this vegetation type, silica and dry matter contents were not significantly related to ME. Pearson correlation coefficients between ME (MJ/kg DM) and nutritional parameters (g/kg DM) for pooled data (n = 297). *** indicates a statistical significance < 0.001 and * indicates a statistical significance < 0.05.

Prediction Equations of ME from Nutritional Parameters
Prediction equations for ME concentration in oak woodland understory and heathlands developed in the present study are shown in Table 1. The stepwise regression analysis indicated that ADF and log DM were the best single predictors of ME in oak woodland understory. The accuracy of ME prediction for this vegetation type improved by adding lignin and silica (as well as their natural logarithm transformations) as predictor variables, and final R 2 reached values between 0.61-0.71. In O-Ru and O-My, the final R 2 value with additional predictors rather than one single predictor rose by over 0.36 and 0.27 units, respectively. Conversely, however, the magnitude of increase for O-Bl was only 0.07 units higher.
On the other hand, equations from multi-step regression analysis for heathlands revealed CP as the best single predictor for ME ( Table 1). The accuracy of prediction for ME could be improved by adding ADF and silica in the latter equations (R 2 increased from 0.60-0.70 and from 0.30-0.41 in Atlantic and continental heathlands, respectively). For all vegetation types, each parameter included on the finally selected models had a highly significant effect on the relationship (p < 0.05), except the intercept of the model fitted for CH (Equation (5b)); therefore, the model was refitted without intercept (Equation (5c)).
Stepwise regression analysis of the data uniquely by season (not taking into account the particular vegetation type) indicated that the best single predictors of ME were DM for nutritional data obtained in Spring and Summer, and ADF for those from Fall and Winter ( Table 2). The accuracy of prediction increased by adding CP and fractions of fiber (lignin and silica) in the latter equations. Those equations showed the lowest R 2 values for Summer (0.57) and better adjustments for Spring and Winter (R 2 0.71 and 0.73, respectively). Again, for all seasons, each parameter had a highly significant effect on the relationship (p < 0.05).
The study of residuals allowed the comparison of the degree of change in the accuracy of ME estimations among the prediction equations obtained in this study. The use of the prediction equation obtained from total pooled data (i.e., general model or Equation (6d)) to estimate the ME in each vegetation type showed RMSE values higher than those found when operating with specific prediction equations for each vegetation type (Table 3). The highest increase of error was shown for Atlantic heathlands (29%) and the lowest for continental heathlands (2%). The rise in error when using the Equation (6d) to estimate ME for each type of oak woodland (i.e., instead of its specific equation) was around 11%. Similarly, the use of prediction equations obtained by season to estimate ME showed a lower RMSE than when the general model was applied. The expected loss of accuracy from using the general model instead of the prediction equations obtained for each season is expected to be especially high for estimates of ME concentrations in the winter (almost 25% RMSE increment). For ME estimates of vegetation in the Fall and Summer, the use of a general model versus a seasonal specific model would represent a loss of accuracy of 4%, which becomes only 1.8% in Spring.

Comparison of Predicted ME from Equations Developed in the Present Study and Using IVOMD as a Sole Predictor
The range of predicted ME concentrations in oak woodlands narrowed with respect to observed ME concentrations (3.4-12 MJ/kg DM for O-Bl, 4.3-11 MJ/kg DM for O-Ru), except for O-My (2.3-13.8 MJ/kg DM). In Atlantic heathlands the minimum and maximum values increased for predicted ME (2.7-12.9 MJ/kg DM), and for continental heathlands the amplitude of range was reduced to values between 2.6 and 6.5 MJ/kg DM.
The results of graphing the observed ME against predicted ME (equations developed in the present study) revealed the degree of disparity expected in the ME values obtained when using the prediction equations as opposed to one based on its IVOMD as a sole predictor (Figures 3-5).       A clear tendency towards overestimation was observed for the plant group of heaths with the use of the general model ( Figure 4A) that is partially corrected with the use of the specific models ( Figure 4B) for each vegetation type (i.e., Equations (1c), (2e), (3c), (4d) and (5b)). This same tendency can be followed by season in Figure 4C.

Establishment of the Prediction Equations
In the current study, we found that the content of ADF in combination with lignin or silica, as well as with DM, can be optimal predictors of ME for the understory of oak woodlands; in the case of heathlands, CP joined ADF and silica as a significant ME predictor. This is in agreement with a variety of research conducted on the effect of chemical composition on forage digestibility and metabolizable energy, which showed that fiber content and CP can be optimal predictors of ME [35,[46][47][48][49][50][51]. Because digestibility and CP content are important components of quality in diets, whereas fiber is usually detrimentally associated with digestibility [52][53][54], CP is often positively associated with ME while fiber is negatively related with ME in the equations [35,46,47,50]. Phuong et al., (2013) reported that fibrous parameters of feed, especially ADF (but also NDF), reduce energy-use efficiency because they contain limited energy to sustain microbial growth, decrease availability of nutrients in high fiber diets and tend to accelerate rumen filling and limit feed intake [35]. The relationships of ADF, silica and CP with ME in the prediction equations here obtained showed a consistency with those findings. Nevertheless, the positive relation of lignin with ME in our study was not coherent with the fact of it being an undigestible fraction in plants. According to Van Soest (2006), correlation of lignin contents and degradability manifest poorly across a large amount of research papers, and it is not uncommon that the association of lignin with digestibility can be difficult to demonstrate [55].
Studies on the prediction of the digestibility of forages (and therefore extended to the prediction of forage energy concentrations) underline stepwise regression analyses in which the digestibility of organic matter is most accurately predicted by equations that include cell-wall constituents and acid-detergent lignin and silica [35,46,47]. The relationship between silica and the digestibility of forages showed that silica can be a major factor in limiting nutritional value [46], and the compilation of studies revealed an evident negative association of silica with IVOMD [55]. An average decline of 3.0 units of digestibility per unit of silica in the dry matter was found in grasses [46], and Smith et al. (1971) reported that silica decreases the digestibility of organic matter by approximately 1% for each percentage increase in the silica content of the forage [47]. This negative relationship between silica and digestibility led to the consideration of acid-insoluble ash (mainly silica) as a good internal marker to estimate dietary digestibility in mammals and poultry (Van Keulen and Young, 1977, Vogtmann et al., 1975) [56]. The use of other internal markers, such as crude fiber and hydrolysis-resistant organic matter, have also been recommended for measurement of apparent digestibility coefficients (De Silva and Anderson, 1995) [56]. Our results showed a noticeable increase of 0.20-0.30 units in R 2 for oak woodland prediction equations (Equations (2b) and (3b)) when DM was entered as a predictor in the equations. Silica was incorporated later in O-My and in Atlantic heathlands, improving slightly the accuracy of the ME final equations. Although the contribution of silica in the raising of R 2 in our prediction equations was small, this parameter should not be overlooked, because estimates of dietary digestibility obtained using the acid-insoluble ash method were in tolerable agreement with those obtained using the standard total fecal collection method [57], which represents a laborious method of determination of in vivo diet digestibility.
The success of acid-insoluble ash as a suitable digestibility marker has depended on the amount of silica content in the forage, with recommendations of a minimum of 7.5 g/kg on a dry matter (DM) basis to get accurate measurements [58]. Based on our silica results, most of the plant species included in the prediction equations surpass that mean value, except some forbs such as Asphodelus albus (associated with forest stands previously burnt) and grasses like Agrostis curtisii. Therefore, if those ME equations with silica as a predictor were to be applied to vegetation types where these plant species are very dominant, some caution should be taken, and perhaps alternative equations with no silica as a predictor could be more suitable (i.e., Equations (4b) and (2c)). Finally, the inclusion of silica in the equations could have some interesting management implications as acid-insoluble ash has been also used to determine feed intake in grazing animals [58], a factor that is related to the continuance or permanence of pasture.

Accuracy and Applicability of the Prediction Equations
In predicting a numerical outcome with a statistical model, the residuals provide an idea of the differences between predicted and observed values (i.e., in our study values of ME using the AFRC formula with IVOMD as a sole predictor). The analysis of the residuals to test for the accuracy of the prediction equations obtained in this study revealed that the use of the general model (Equation (6d) from pooled data) as opposed to the use of the specific prediction equations for each vegetation type (i.e., Equations (1c), (2e), (3c), (4d) and (5c)) or for each season, would mean a loss of precision in the ME estimation (Table 3). Accuracy, expressed as a combination of RMSE and bias and MAE values (a measure of errors between paired observations), was higher when the general model was applied, and the magnitude of the loss of precision was superior for Atlantic heathlands (RMSE increase of 29%). ME estimation in continental heathlands and oak woodlands would lose 1 and 11% of accuracy, respectively, if the general model were chosen. On the other hand, the degree of total loss of accuracy in predicting seasonal ME with the general model as opposed to the use of the specific equation for each season would be much smaller (5.4%). Only Winter ME estimates are expected to be substantially influenced by the use of one or another (the expected loss of precision being around 25%). RMSE increased by only 1.8% for Spring prediction when applying the general model, whereas for Fall and Summer the precision loss with the general model is expected to be around 4%. Because the quality of the fit does not necessarily reflect the quality of the prediction, the cross-validation of each model was examined. The equivalent estimate to the R 2 of the fitting phase, namely MEF in the cross-validation, followed a similar tendency, and showed again for those models fitted specifically for each vegetation type higher values than those where the general model was applied. Therefore, the quality of the prediction is expected to improve when using the specific prediction equations as opposed to the general model. Conversely, however, cross-validation performed for each season revealed MFE values approaching closely those obtained from applying the general model (and only still differentiated for Winter), thus confirming once again that only Winter ME estimates are expected to be substantially improved by the use of the specific season model as opposed to the general model. Shrubs, and especially those in the Ericaceae family, show a nutritional composition characterized by a very small seasonal variation [17].
The range of predicted ME concentrations in oak woodlands and heathlands narrowed with respect to the observed ME concentrations because the minimum values of the range increased (+0.6 for O-Ru, +1.2 for O-Bl, +0.8 for both heathlands) and maximum values were still similar to observed ME (around 13 MJ/kg DM). The results of graphing the observed ME against the predicted ME revealed the degree of dispersion of the data when using a prediction equation based on nutritional parameters (ADF, lignin, silica, DM, and CP) as opposed to the one based solely on IVOMD. Our attempt to determine whether plant groups (forbs, grasses, heaths, other shrubs and trees) were differently affected by the usual discrepancy between observed ME and predicted ME showed a clear tendency towards the overestimation of heaths when using the general model. Such overestimation could be partially corrected with the use of the specific models for each vegetation type. This is important to take into account in the case of ME predictions performed with the general model for heathlands that are highly dominated by heaths, because the observed ME data obtained from AFRC (1993) could already be overestimating the actual ME concentration of heaths. Yan and Agnew (2004) reported an over-prediction of ME concentration of 1.2 MJ/kg DM when using AFRC ME estimations for poor quality silages, the under-prediction for very good quality silage being 1.3 MJ/kg DM. As an example of possible implications for management, efficient reducers of fuel load like goats have ME requirements for maintenance around 7.4 MJ/kg DM (0.424 MJ/kg BW 0.75 , taken from NRC 1981), for an average animal of 45 kg. In the case of beef cattle, similar values of 7.6 MJ/kg DM have been approached for Winter season, raising in general as an average of 8.6 MJ/kg DM (NRC 1996) for a mature animal of 1000 lb (453,592 kg). In both continental and Atlantic heathlands, the observed ME average values were under the minimum requirements for maintenance because of a high percentage of heaths and Halimium spp., and such overestimation could result in a substantial divergence facing management decisions (for example of the optimum animal stoking rates).
There is little literature discussing ME concentrations of browse considered in our study. Most data on pastures' ME concentrations are relative to fresh grass forages with ME values between 7.5-12.6 MJ/kg DM [28,35], or grass silages with ME around 13.6 MJ/kg DM [33]. Over 2/3 of the plants included in Atlantic and continental heathlands showed observed ME values below 5 MJ/kg DM, which are similar observed ME figures to those observed in the Ericaceae plants reported by Paton (2003) [50]. Other groups of plants studied here showed similar ME concentrations to those found in the same study. For example, gorse such as U. minor and Ulex europaeus included in Atlantic heathlands reached ME concentrations of 8.5 MJ/kg DM and between 6.9-7.2 MJ/kg DM, respectively. Gorse is highly selected by horses in heathlands, decreasing its biomass and dominance and leading to the reduction of fire risk as well as increasing plant diversity [4][5][6]31]. Unlike in heathlands, observed ME concentrations in the three types of oak woodlands revealed more records with values superior to 7 MJ/kg DM. ME concentrations of around 8 MJ/kg DM obtained for the leguminous shrubs Cytisus striatus and Genista florida in our study were similar to those found in plants of the same genera in Castro et al. (2021) [12]. The highest number of plants reaching minimum ME concentration for maintenance occurred in O-Bl, and were mainly associated with the shrubs Hedera helix, Lonicera periclymenum and Ilex aquifolium; to the forbs Asphodelus albus and Teucrium scorodonia and to grasses. These are relevant plant species highly preferent in roe deer diets [59,60], and observed ME concentrations in the present study are also in the range reported by Paton (2003) [50]. Some forbs in O-Bl and O-My showed higher ME concentrations because of superior nutrient content and better digestibility, and were close to reaching in Spring the minimum ME requirements of goats for late pregnancy (13 MJ/kg DM). This would balance the low ME values provided by other wooded plant species.
The stepwise regression analysis of the present study generated prediction equations of ME based on ordinary nutritional parameters that are easily obtained from laboratories (i.e., ADF, lignin, crude protein, dry matter), which could facilitate ME estimations in grazing management scenarios where absolute accuracy is not essential. This could greatly simplify grazing management decisions in an extensive forest grazing setting with a wide variety of plants, being an alternative to in vivo digestibility trials and in vitro digestibility analysis, which are not always available and are more difficult to conduct.
There has been a growing interest in predicting the energy value of diets by establishing relationships with chemical parameters that are easily analyzed in laboratories, but, to the best of our knowledge, this is the first time that regression equations are presented in order to predict the ME of the main Atlantic pedunculate oak woodlands and heathlands.
The wide-ranging information on the nutritional characteristics of the plant species from which the ME predictive equations have been constructed provides a reference source when these equations are to be applied in similar vegetation types. Prediction equations in the current study have been developed by using a great range of browse quality, also including some seasonal variation of nutritional characteristics, so that they can be broadly applicable in similar vegetation types depending on the available data on chemical composition. Our results showed that the prediction accuracy of ME is reduced when the general model is applied, as opposed to the specific prediction equations for each vegetation type, especially in the case of Atlantic heathlands. The general model tends to overestimate the ME concentrations in heaths when compared to the ME values obtained from IVOMD as a predictor (i.e., AFRC ME estimation), and this divergence could be corrected to some extent by applying specific prediction equations for each vegetation type. Although the use of prediction equations by season would improve accuracy in the case of a Winter scenario, using the general model as opposed to the prediction equations for Spring, Summer or Fall would represent a much smaller loss of accuracy.
Finally, some of the prediction equations here obtained included silica as a predictor, which has been reported as a reliable digestibility internal marker, especially when the content in plants is at least 7.0 g/kg DM. The range of silica values provided for different species in the present study will also serve as a reference to decide on the adequacy of applying equations that include silica for predicting ME in similar vegetation types. Data Availability Statement: Data sharing is not applicable to this article.

Acknowledgments:
The previous laboratory and field work, from which the current study was later based on, was funded by the Rural Development Agency of Galician Government (Xunta de Galicia, Spain).

Conflicts of Interest:
The authors declare no conflict of interest.
Appendix A   Table A1. Nutritional composition, in vitro digestibility and metabolizable energy (g/kg DM, g/kg, MJ/kg DM) of plant species in oak woodlands of NW Spain. For each species, first value corresponds to mean value, followed by standard deviation and coefficient of variation (in bold).