Next Article in Journal
Thermotherapy and Storage Temperature Manipulations Limit the Production of Reactive Oxygen Species in Stored Pedunculate Oak Acorns
Next Article in Special Issue
Long-Term Ecosystem Nutrient Pool Status for Aspen Forest Harvest Simulations on Glacial Till and Sandy Outwash Soils
Previous Article in Journal
Plant Regeneration via Somatic Embryogenesis in Larix principis-rupprechtii Mayr
Previous Article in Special Issue
How to Evaluate Downed Fine Woody Debris Including Logging Residues?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Biomass Yield and Economic, Energy and Carbon Balances of Ulmus pumila L., Robinia pseudoacacia L. and Populus × euroamericana (Dode) Guinier Short-Rotation Coppices on Degraded Lands under Mediterranean Climate

Department of Agroforestry Sciences, School of Engineering, University of Huelva, 21071 Huelva, Spain
*
Author to whom correspondence should be addressed.
Submission received: 5 September 2021 / Revised: 21 September 2021 / Accepted: 28 September 2021 / Published: 29 September 2021

Abstract

:
The steadily increasing demand for energy and concerns about climate change have prompted countries to promote the use of renewable energy sources, including lignocellulosic biomass. In this context, this work aims to assess the biomass production for energy purposes in crops with short rotation, as well as its effect on soil properties. Deciduous tree species were used, mainly Siberian elm (Ulmus pumila L.), black locust (Robinia pseudoacacia L.) and a hybrid poplar clone (Populus × euroamericana (Dode) Guinier, clone ‘AF2’). Four field trials were implemented, under two different types of Mediterranean climate, where highly productive taxa were tested, in addition to the mixed planting of a nitrogen-fixing species with a non-fixing one. Short-rotation coppices (SRCs) of these taxa yield about 12–14 t ha−1 year−1 of high-quality dry woody biomass, when fertilizers and irrigation water are supplied; generate 205–237 GJ ha−1 year−1 net and earnings of about EUR 1.5 per EUR 1 invested; and sequester into the soil 0.36–0.83 t ha−1 year−1 of C and 57 kg ha−1 year−1 of N. Therefore, these species raised as SRCs could improve degraded soils if the crop is properly managed, resulting in favorable economic, energy and CO2 emission balances. The use of mixed plantations can bring economic and environmental gains, and the biomass transformation into high-quality chips or pellets gives it added value.

1. Introduction

In the context of climate change, international institutions have promoted the use of clean and renewable energy sources, including biomass, to protect the global environment and reduce greenhouse gas emissions [1,2]. Currently, biomass accounts for more than half of all global renewable energy production and about 11% of all energy sources [3]. For instance, the Energy Roadmap 2050 of the European Commission has set as a goal to reduce CO2 emissions that renewable energy sources increase their contribution to gross final energy consumption by at least 55% [4]. As a consequence, the use of lignocellulosic biomass is steadily expanding and can play an essential role in achieving the renewable energy goal in 2050. It is regarded as a flexible primary energy due to the value-added products derived from it, such as chips, pellets, biogas, bio-based materials and chemicals [5,6,7], apart from traditional uses. Additionally, the increasing demand for woody biomass for industrial purposes such as wood pulp and wood-based panels [8] also contributes to expanding its use. Therefore, in order to meet the growing global demand, plantations of fast-growing species could supplement the additional woody biomass that traditional agroforestry systems cannot supply, and constitute a fundamental complement in the biomass market, since they provide an alternative use for abandoned land, contributing to rural development while reducing CO2 emissions [5,9,10,11,12]. Furthermore, these plantations could reduce some negative effects of biomass exploitation in traditional forest systems, such as the depletion of soil fertility caused by nutrient removal in whole-tree harvesting, as well as the alteration of habitat and biodiversity in sensitive sites [13,14,15,16]. Short-rotation coppice (SRC) plantations consist of a woody cultivation in short cycles (usually less than 7 years) of fast-growing species capable of sprouting from the stump after being cut, and provide a suitable means of obtaining this raw material in short periods of time [17].
Currently, different plant taxa are managed worldwide as SRC crops (e.g., the genera Populus, Salix, Eucalyptus, Paulownia, Gmelina, Robinia, Casuarina, Leucaena, etc.), and also some species such as Ulmus pumila L., Platanus × hispanica Mill. ex Münchh. and Ailanthus altissima (Mill) Swingle. have been studied for woody biomass production [18,19,20,21]. Most of the above-mentioned species are multiple-purpose species providing numerous goods and services (e.g., biofuel, carbon sequestration, soil restoration, wood pulp, forage, etc.). Under optimal environmental conditions, commercial plantations can reach up to 25 t ha−1 year−1 dry matter with poplars and willows, around 40 t ha−1 year−1 with eucalypts, Paulownia, Casuarina and Leucaena, and 15 t ha−1 year−1 with black locust (Robinia pseudoacacia L.) [11,12,18,22,23]. However, adverse environmental conditions, such as the winter frosts of inland and high-altitude continental areas as well as the typical summer drought of Mediterranean climates, compromise the growth and survival of many of these taxa to such an extent that if they survive, the production can drop well below 5 t ha−1 year−1. Moreover, most of these woody energy crops will be established on degraded soils, and they will be required to produce a large amount of biomass profitably [12,24]. Farmland degradation is a serious problem today. Currently, in Europe alone, 45% of farmland contains less than 2% organic matter, 15% suffer from problems derived from excess inorganic nitrogen fertilization, and significant signs of erosion can be observed in more than 147 million hectares [25,26]. At least 70% of all these degraded and low-fertility lands are under a climate where Robinia pseudoacacia, Ulmus pumila and/or Popupulus sp. can grow. However, special attention must also be paid to the quality of the biomass grown on contaminated soils (e.g., uncultivated land or post-industrial soils), since apart from nitrogen contamination, other elements, including heavy metals, may affect plant growth, biomass quality and industrial processes for transforming it into various types of fuels [27].
Recent studies have shown that Robinia pseudoacacia, Ulmus pumila and Populus sp. stand out in areas under a Mediterranean continental climate with cold winters, in Central Europe and North America, as they can tolerate winter frosts well below −20 °C [19,28,29,30,31]. The latter has a great capacity to adapt to different environmental conditions, as well as a high hybridization capacity and ease of vegetative propagation. These characteristics, together with its rapid growth, have contributed to the development of a very extensive clonal offer, which requires the prior selection of each site and specific use. However, this genus requires wet soils, while R. pseudoacacia and U. pumila are alternatives due to their great drought tolerance [18,19,32], but they must be managed with care since they can become invasive and alter the native plant habitats [33,34]. Ulmus pumila and Robinia pseudoacacia have not been subjected to advanced genetic improvement programs, nor do they offer selected clonal plants. For the latter species, studies in this regard were started a few years ago [35,36] in central Europe, but the clones are not adapted to the arid edaphoclimatic conditions of the Mediterranean basin, nor are they available on the market. Consequently, these two species require selection, improvement and breeding programs due to the high growth variability shown by R. pseudoacacia in trial plots established up to now, and the difficulties of vegetative propagation of U. pumila, mainly in vitro [37,38].
Thus, it is necessary to evaluate the biomass yield and the environmental effect of an intensive production system such as SRC on degraded soils [7,39,40], since the knowledge in Mediterranean environment is still scarce and depends on the cultivated species, cultivation practices and environmental factors [11,12,41]. In this respect, a N-fixing species such as Robinia pseudoacacia can be used in SRC energy crops [19,42], since it can contribute to decreasing inorganic nitrogen fertilization, and as a result to reducing economic and energy costs [43]. Moreover, many studies have shown that mixtures of species can be more productive than monocultures because species’ interactions can increase resource availability as well as resource use efficiency (i.e., less intense inter-specific than intraspecific competition). An example of this is the mixture of nitrogen-fixing and non-nitrogen-fixing species, mainly if soil nitrogen is limited [19,22,44,45].
This manuscript compiles four field trials carried out in the Mediterranean environment that assessed the biomass production and the evolution of the top layer of the soil after being subjected to intensive SRC energy crops. It aims to evaluate whether SRC plantations of fast-growing tree species can improve a degraded soil while producing high-quality biomass. Soil properties were assessed under hybrid poplar (Populus × euroamericana (Dode) Guinier), black locust and Siberian elm (Ulmus pumila L.) subjected to different cultivation treatments. The main objectives are: (a) to evaluate plant growth and biomass production of hybrid poplars, black locust and Siberian elm subjected to different SRC treatments in two contrasting Mediterranean environments; (b) to assess their appropriateness for use as bioenergy; (c) to determine the changes of soil properties after a cropping period of 5 to 8 years; and (d) to analyze carbon sequestration and the economic and environmental viability of this production system.

2. Materials and Methods

2.1. Experimental Design and Plant Material

This article shows the research carried out in four independent field trials, which have several aspects in common: the location was the same for three of them; two species, Robinia pseudoacacia L. and Populus × euroamericana (Dode) Guinier, were used in the four trials, and Ulmus pumila L. in three of them; the crops were developed in Spain under a Mediterranean climate; they were designed as short-rotation coppices (SRC) for energy purposes; and both the irrigation and fertilization also coincided to a great extent. Different intraspecific taxa were used for poplars and black locust because the tolerance to environmental stresses and the nutrient requirements of these species can vary significantly [19,31,36]. Two sites with marginal agricultural lands were selected, located in the provinces of Granada (Los Morales—Mo; UTM, zone 30S, X:535457, Y:4191514, 1084 m.a.s.l.) and Huelva (La Rábida—Ra; UTM, zone 29S, X: 685054, Y: 4119168, 7 m.a.s.l.). Three field trials were established at Los Morales (Mo1, Mo2 and Mo3) and one at La Rábida (Ra4). The main climate conditions and soil characteristics are shown in Table 1 and Table 2 and Table S1, respectively. Both sites have a typical Mediterranean climate with hot and dry summers, but with cold winters in Los Morales and mild winters in La Rábida, which makes the vegetative period more than two months longer in Ra than in Mo. The soils were permeable alkaline soils (pH = 8.2–8.7) with low salinity (EC ≤ 331 μS cm−1), appropriate C/N ratio, high Ca and Mg contents, deficient in P, Fe, Mn and Zn, and with medium (Los Morales) or low (La Rábida) levels of OM, CEC, N, K and active limestone. The stoniness was 19.2% and 23.4%, and the bulk density was 1.32 and 1.55 kg dm−3, for Los Morales and La Rábida, respectively. Previously, the land had been dedicated to the cultivation of rain-fed barley, and it did not present problems of contamination by heavy metals or other nutrients. The N, P and K contents in the irrigation water were 6.1, 0.15 and 0.22 mg L−1, respectively.
A randomized complete block design was established in each field trial, with four blocks per trial. Each taxon was present in each block within an experimental unit. Therefore, there were four experimental units in each field trial, which were called replicates. The experimental unit was of a different size in each of the four field trials. It consisted of 3–4 parallel lines containing 6–30 plants per line, with a spacing of 0.6–1.0 m between plants and 2.0–3.3 m between lines; the crop density ranged, therefore, from 3787 to 5556 plants per hectare (Table 3, Figures S1–S3). In the Mo1 and Mo2 trials, three adjacent uncropped areas with natural vegetation (predominantly weeds) were selected as controls for the purpose of taking soil samples during the study period (Figure S1).
Before planting, soil preparation and the eradication of pre-existing weeds consisted of a deep plowing at 50 cm, shallow tillage (harrowing), pre-emergent herbicide application (Oxifluorfen 24%, 2.5 L ha−1), and a standard fertilization with N (75 kg ha−1), K (62.3 kg ha−1), P (32.7 kg ha−1) and S (24 kg ha−1), corresponding to a 15–15–15 (12 S) ratio fertilizer (N–P2O5–K2O, plus 12 SO3). Subsequently, after planting, fertilizers were added by fertigation: N (50 kg ha−1 year−1), K (41.5 kg ha−1 year−1), P (21.8 kg ha−1 year−1), S (16.0 kg ha−1 year−1), 3 kg ha−1 year−1 of Fe–Zn–Mn (4.5%–1.5%–0.5%) chelated by EDDHA, 100 g ha−1 year−1 of Cu (14%) chelated by EDTA, 50 g ha−1 year−1 of H3BO3 and 5 g ha−1 year−1 of (NH4)2MoO4. Sulfur was supplied for facilitating the absorption of Fe, Zn and Mn, since they seem to be somewhat locked-away from plants, as symptoms of foliar interveinal chlorosis were observed. The only experiment that deviated from this general fertilization guideline was Mo2, since it was designed with the main purpose of evaluating treatments that include a mixture of a nitrogen-fixing species (Robinia pseudoacaia) with a non-nitrogen-fixing species (Populus × euroamericana). However, it was decided to supply nitrogen in the fertilization schedule of the Mo2 trial by reducing the amount applied by half (i.e., 37.5 kg ha−1 of N before planting and 25 kg ha−1 year−1 from the second year on), because the viability of the R. pseudoacacia-nitrifying bacteria’s symbiosis was not totally guaranteed, since it was the first time this species was planted on this site; the proportions of nitrogen in the plant derived from the atmosphere vary among seed sources [46]; and plant growth had to be facilitated in these nutrient-poor soils. The other nutrients were supplied at the same quantities as in the other trials. The experiments Mo1, Mo3 and Ra4, on the contrary, had as their main objective to compare the productivity of different taxa, belonging to the same or different species.
The Mo1 and Mo2 experiments were established in 2011, while Mo3 and Ra4 were in 2015 (Table 3). At the time of planting (from mid-March to mid-May), the vegetative material consisted of:
  • Bare root nursery seedlings, 50–70 cm in length, of two Robinia pseudoacacia L. taxa (the improved cultivar “Nyírségi” and commercial nursery plants from Germany), Populus nigra L., Ailanthus altissima (Mill) Swingle., and Ulmus pumila L;
  • Paulownia plants of 15–20 cm in height (Paulownia fortunei (Seem.) Hemsl., clone “UHU”) coming from root cuttings;
  • Bare root nursery plants, 50–70 cm in length, coming from rooted hardwood cuttings of Platanus × hispanica Mill. ex Münchh;
  • Poplar hardwood cuttings of 20–25 cm in length belonging to five hybrid clones—four clones of Populus × euroamericana (Dode) Guinier (clones “Adige”, “AF2”, “Oudenberg”, “I214”), and one clone of Populus × interamericana van Broekhuizen (clone “Raspalje”).
Weeds were mechanically controlled again three months after planting, but fortunately it was unnecessary to control them thereafter, due to the high planting density and the rapid growth of the new sprouts after harvesting. Every year, from early June to mid-September and according to the rainfall, 160–480 mm of water was supplied by drip irrigation in order to offset summer drought and achieve an annual water supply (rainfall + irrigation) of around 700–800 mm. The adjacent unmanaged areas did not receive any agronomic practices (soil preparation, weeds control, fertilization, irrigation, etc.), but they were left under their natural vegetation. The application of pesticides was unnecessary in the Ra4 trial; however, in the Mo1 and Mo2 trials, occasionally during the whole period studied (only in two years) a pest caused by Xanthogaleruca luteola that affected Ulmus pumila had to be controlled, by spraying 0.5 L ha−1 of Deltamethrin at 2.5% every time.

2.2. Shoot Growth and Biomass Assessment

Stem diameter (D, measured 10 cm above ground), height (H) and the number of stems from the stump (NSt) were annually measured in 4–20 plants per replicate depending on the trial. In order to avoid the edge effect, only plants belonging to the middle line were measured (Table 3, Figures S1–S3). Once D and NSt were measured, the quadratic mean diameter (QMD, i.e., the average diameter of the plant per plant basal area) and the plant basal area at the measurement height of D (BA, m2 ha−1) were obtained. The basal area of the plant (BA) was calculated by adding the cross-sectional area of all of the main stems of the plant and dividing by the area of land occupied by the plant (Aplant), where Aplant is the result of dividing 10000 (m2 ha−1) by plantation density (plants ha−1).
For the Mo1, Mo2 and Ra4 trials, two short rotations were analyzed by annually measuring plant growth for 3–4 years after planting (1st rotation), and also measuring growth every year in re-sprouted shoots for 3–4 years after cutting (2nd rotation), whereas only the 1st rotation was analyzed for Mo3 (Table 3). In addition, phenology was assessed on a monthly basis by recording the presence of reproductive organs (developing or ripe), and the status of the twigs and buds (developing or resting).
The plant age was identified according to the following convention: RijSk, where (i) indicates how many times a stump has been cut, (j) the age in years of the stump-roots, and (k) the age of the main shoots at the time of measuring. For instance, in the Mo1 experiment, R04S4 means a 4-year-old plant that has never been cut after planting, while R18S4 means a plant whose roots are 8 years old, but it was cut after the fourth year, so the aboveground part is only 4 years old.
Biomass production was assessed at the end of winter (February to March), before spring bud break, from the measured data of D and H. For this purpose, specific allometric equations relating D and H to the shoot dry weight were developed for these plants, similar to those performed in other studies [20,47,48,49], by cutting and weighing a sample of shoots with a stem diameter (D) ranging from 10 to 140 mm. The dry weights of aboveground woody (AGWB) and non-woody (AGNWB) parts were assessed separately, after being oven-dried at 65 °C. As leaves had fallen several months ago, non-woody biomass was reduced to dehydrated ripe legumes, and was only present in Robinia pseudoacacia. Total dry weight (AGTB) was determined as AGWB + AGNWB. Since all studied taxa belonged to deciduous species, none of them maintained leaves on the measurement dates (February to March), which coincide with the dates on which the harvest takes place. Different equations were tested with combinations of D, H and/or D2 as predictor variables, and most of them yielded an excellent fit. However, since the stem diameter is very easy to measure and the power equation using D as the only predictor (AGWB = a·Db) also showed a very good fit (Table 4), it was thus used for estimating AGWB. The strong correlation (r > 0.90, p < 0.001) obtained between D and H for the type of stems assessed (Table S2) allowed us to remove H from the biomass weight estimation equation. A single allometric relationship was used to evaluate the AGWB, on a dry weight basis, of the same taxon in any of the trials (Table 4).

2.3. Physical–Chemical Properties of Soil and Plant Biomass

Soil samples were taken at the beginning of the assays (original soil) and also at the end of every rotation, in March–April, under four taxa (Populus × euroamericana “AF2”, Ulmus pumila, Robinia pseudoacacia and Ailanthus altissima). Four soil subsamples per replicate and taxa (0–30 cm depth) were taken with a soil core between planted lines, 1 m apart from the central line. The litterfall and soil mineral layer were then separated, and the subsamples belonging to each taxon/treatment and replicate were pooled. For the Mo1 and Mo2 experiments, soil samples in adjacent uncultivated areas (control areas, Figure S1) were also taken. The dry weight of the litterfall layer (i.e., leaves, shoots and reproductive organs partially decomposed) was therefore also quantified.
Litterfall was washed with distilled water, oven-dried at 65 °C to a constant weight, ground and stored at room temperature (15–20 °C) in sealed containers until being analyzed. Soil samples were air-dried at room temperature for one month, milled and sieved (2 mm mesh). Concerning the plant biomass (roots and aboveground part), four shoots per replicate of Populus × euroamericana “AF2”, Ulmus pumila, and Robinia pseudoacacia were randomly sampled at the end of every rotation and pooled for analysis. In the same way as for litter, these samples were washed, oven-dried at 65 °C, ground and stored at room temperature, and their properties were analyzed. The physical–chemical characterization of litterfall, plant biomass and mineral soil was performed following standardized methods (Table S3).

2.4. Valuation of Biomass and Determination of the Economic, Energy and CO2 Balances of the Production Process

Aboveground woody biomass was ground (Woodstock 3PH, Smartec®, Italy) and sieved to a 0.2–5.0 mm particle size in order to create homogenous samples. Then, for the purpose of increasing the economic value of the biomass yielded, pellets were manufactured using a pelleting press (PLT-400, Smartec®, Italy). The moisture content was set to 130 g kg−1 with a bulk density of 183, 268 and 216 kg m−3 for “AF2”, R. pseudoacacia and U. pumila, respectively, and an operating temperature from 95 to 105 °C. Die channels six millimeters in diameter were used; the inlet part of the channel was a cone-shaped opening 3.5 mm deep with 70° angles; the active part was 26 mm long for “AF2” and 22 mm for R. pseudoacacia and U. pumila, so the compression ratios were, respectively, 4.33 and 3.67. After that, the ash content, bulk density, length and diameter, mechanical durability, chemical composition and moisture content of the pellets were assessed according to the standard ISO 17225-2:2014 [50].
Without attempting to delve deeper into the socio-economic and environmental analysis, this paper provides a brief environmental, energy and economic analysis of a biomass supply chain based on an experimental poplar/black locust/Siberian elm plantation grown as an SRC in an area under a Mediterranean climate. Calculations consider a plantation cycle of 15–17 years with 4–5 rotations (harvests) every 3–5 years. The rates and wood price, fuel consumption, and machinery performance (Table 5) have been obtained from the companies TRAGSA [51] and ENCE [52], two of the largest Spanish forestry companies, and from recent studies carried out by other authors and institutions [11,12,43,53,54,55,56,57]. We have also taken into account the calculation methodology used by previous studies [58,59,60,61,62,63,64]. Soil preparation, herbicide, plants and planting, harvesting, materials and the installation of irrigation system, grubbing and land clearing were considered as fixed costs, whereas irrigation (water and energy), fertilizer and transport were taken as variable costs. Financial costs, land, and extra costs have not been considered.

2.5. Data Analysis

The dependent variables (DV; i.e., D, H, NSt, etc.) were evaluated in the same plants for several years during the first two rotations after planting; therefore, our data structure resulted in repeated measurements for each plant, and repeated measures ANOVAs were applied [65]. Hence, as the within-plant observations were autocorrelated, a model in which the plant (within-replicate) was considered a random effect was used. The first and second rotations were analyzed separately. Taxa (TX, in the Mo1, Mo3 and Ra4 trials) or treatment (TR, in the Mo2 trial) were included as fixed effects. The annual growth was evaluated by introducing the evaluation period (year, Y) as a fixed effect. The replicate (Rep) growth differences were also assessed, which represents the block effect. The model also included the pairwise and triple interactions between main effects as fixed effects. It was tested in advance that the data met the assumptions of normality and sphericity. To evaluate the among taxa/treatments comparisons, the Tukey HSD test was used in order to differentiate the homogeneous groups. Significant differences were established at p < 0.05. Thus, the full models had the following structure:
DVijkl = Repi + Yj + TXk + (Rep × Y)ij + (Rep × TX)ik + (Y × TX)jk + (Rep × Y × TX)ijk + εijkl
or,
DVijkl = Repi + Yj + TRk + (Rep × Y)ij + (Rep × TR)ik + (Y × TR)jk + (Rep × Y × TR)ijk + εijkl
Soil properties were analyzed by a two-way ANOVA (Mo1, Mo2 and Mo3 trials) by considering the fixed effects of taxon/treatment (TX—each of the taxa evaluated and natural vegetation; TR—each of the mixture treatment assessed in Mo2 trial) and rotation (Year). A three-way ANOVA was conducted for the Ra4 assay by considering TX, soil layer (SLY: 0–15 cm and 15–30 cm) and rotation (Year) as fixed effects. Their pairwise and triple interactions were also included in the model as fixed effects. Since soil sampling did not take place at exactly the same point every year, a repeated measures ANOVA was not applied. A one-way ANOVA was applied to the nutrient content of litterfall and biomass sampled at the end of the studied periods, the taxon (TX) being the fixed effect. Significant differences were also established at p < 0.05. To evaluate the among-taxa/treatments and year comparisons the Bonferroni test was used. The free R software and the SPSS 19.0 software (IBM® SPSS Statistics®) were used for data analysis.

3. Results

3.1. Shoot Growth and Biomass Production

Mortality was very low in all trials (<5%) during the study period, except for Robinia pseudoacacia, which was 12% as a whole. It took place primarily during the first summer after planting, so the failed plants were replaced the following fall. No measurement took place on replanted plants. Except for Populus × euroamericana “AF2” and Ulmus pumila, all taxa showed slight symptoms of interveinal foliar chlorosis [66], which were more pronounced for R. pseudoacacia and Populus × interamericana “Raspalje”. This deficiency was corrected by applying chelates of Fe, Zn and Mn. As expected for deciduous species like these, the usual winter frosts in Los Morales (−7 to −9 °C) did not cause damage to the plants. However, an occasional frost of −13.6 °C in January 2012 severely damaged the aerial part of Paulownia fortunei “UHU”, although the plants did not die and new shoots sprouted the following spring. Ailanthus altissima and Ulmus pumila were the taxa that best withstood water stress (caused by a short heat wave event or a short delay in irrigation). The other taxa suffered partial wilting and defoliation.
The increase in height (H) and stem diameter (D), as well as biomass production, differed significantly both between taxa (p < 0.001) and between years (p < 0.001) for all the trials, while no significant differences were obtained between replicates within a trial or for the main effect interactions (p > 0.150). All experiments considered, the average heights achieved by the more productive taxa were 5–7 m, while average stem diameters exceeded 8 cm, reaching up to 12 cm in the Mo3 trial, and maximum diameters of 15 cm were obtained (Figure 1, Table S4). Populus × euroamericana “AF2”, Robinia pseudoacacia and Ulmus pumila stood out for their biomass production in the Mo1 trial.
Overall, the H and D increments were greater in the first two years after planting or after cutting, but this was not the case with the increase in biomass, which was greater during the second–fourth years after planting, or the second–third years after cutting (Figure 1, Figure 2, Figure 3 and Figure 4 and Figure S4). It should be noted that biomass production the first year after cutting (second rotation) was much more vigorous than that obtained the first year after planting (first rotation). In spite of the significantly smaller diameters in the second rotation than in the first one (Table S4 and Figure S5), the number of stems sprouted from the stump (NSt) increased in the second rotation, resulting in larger QMD and basal area and, consequently, increasing the biomass production. On average, NSt ranged from 1.02–1.88 at the end of the first rotation to 2.71–4.89 in the second (only diameters greater than 15 mm were counted). Therefore, QMD was quite similar to D in the first rotation, but more than 80% greater than the mean diameter in the second one. On average, for all species as a whole and the four field trials, BA increased by 55% from the first to the second rotation (from 35.3 to 55.7 m2 ha−1).
Considering Populus × euroamericana “AF2”, Robinia pseudoacacia and Ulmus pumila, the mean annual increment (MAI) was about 9–11 t ha−1 year−1 at the end of the first rotation, but increased to 12.5–14.5 t ha−1 year−1 at the end of the second one. However, current annual increments (CAI) up to 15-20 t ha−1 year−1 were obtained (Figure 4 and Figure S4). On average, the second rotation increased biomass production by 10–40% compared to the first (Figure 1, Figure 2 and Figure 3). Regarding the treatments with a mixture of taxa in the same planting line (Mo2 trial), it should be highlighted that the treatments with a 50–75% proportion of the nitrogen-fixing species improved the total biomass production the first rotation, although this positive effect was a bit diluted in the second rotation (Figure 2). From the second year (Robinia pseudoacacia) and the third year (Ulmus pumila) after planting, the production of reproductive organs began to be important, representing up to 8% of aboveground biomass at the end of the rotation. At harvest time, the litter layers averaged 942 (115), 1027 (103) and 1054 (96) kg ha−1 year−1 of dry matter for Populus × euroamericana “AF2”, Robinia pseudoacacia and Ulmus pumila, respectively, without significant differences among field trials.

3.2. Physical–Chemical Characterization of Soil and Plant Biomass

Many of the physical–chemical characteristics analyzed in the uppermost layer (0–30 cm) of the initial soils (Table 2 and Table S1) did not show significant differences when they were measured again after several years of cultivation (first and second rotation), neither for the main effects nor for their interactions (0.052 < p < 0.970). The only exception was the significant TX × Year interaction for organic matter in the Mo1 trial (p = 0.013). However, some of the properties did show significant differences (Table 6), mainly for Year, but also for TX or TR and soil depth (SLY).
Considering the evolution of the upper 30 cm of the soil during the first and second rotations, for all taxa as a whole, the following trends can be pointed out (Table 2 and Table 7):
  • A significant increase in exchangeable sodium percentage (ESP), organic carbon (OC), N content and available K content. The increase in ESP was evident at the end of the first rotation (3–4 years), but those of OC, N and K took 5–8 years to differentiate from the initial soil;
  • The opposite trend was obtained for C/N ratio and the contents of P and active limestone—that is, a significant decrease over time;
  • For the rest of the properties analyzed in the initial soils, no significant differences were obtained for the studied period (5–8 years, depending on the experiment);
  • In the Ra4 assay, those cases in which differences between the two soil layers (0–15 and 15–30 cm) were significant were due to a higher mineral nutrient content (+ 40–75%), greater OC (+ 130–160%) and lower ESP (−28%) in the shallowest horizon (0–15 cm) than in the deepest.
Finally, regarding the differences in soil properties due to the cultivation of different taxa, we can conclude that: (a) in the Mo1 and Mo2 trials, the few cases in which the differences were significant were due to differences between the cultivated taxa and the non-cultivated land, while the soil under the different cultivated taxa did not differ from each other for the studied period; (b) no significant differences were obtained between taxa in the Mo3 trial for nutrient and OC contents; and (c) in the Ra4 trial, higher N and P contents, greater ESP and lower C/N ratio for Robinia pseudoacacia (mean values 0.07%, 36.6 mg kg−1, 17.1% and 7.92, for N, P, ESP and C/N, respectively, at the end of the second rotation) than for Ulmus pumila and Populus × euroamericana “AF2” (0.06%, 27.3 mg kg−1, 10.1% and 9.35, respectively) were obtained.
Regarding the physical–chemical and energetic properties of the harvested biomass, as well as the roots and the litterfall (Table 8), it should be noted that:
  • No significant differences were found between different trials (p > 0.060), and therefore the mean value from the four trials as a whole has been shown for each taxon;
  • Significant differences between taxa were obtained for N, Ca, S, Fe, Cl and Ash contents (p < 0.040), although these differences were not very high in absolute terms. Robinia pseudoacacia stands out for its high N and low Ca contents, and high wood density;
  • According to the ISO 17225-4 standard [67], the chipped woody biomass is of high quality, since it is harvested without leaves and has a low bark percentage (12% for Robinia pseudoacacia, 15% for Ulmus pumila, and 14% for Populus × euroamericana “AF2”);
  • Due to their physical–mechanical and energetic properties (heating values; bulk density, mechanical durability and moisture of pellets), the highest-quality pellets (i.e., ENplus-A1 for domestic use) could be obtained for the three species studied according to standard ISO 17225-2 [50]. However, the chemical properties (N, S, Cl or ash contents) devalue their quality to commercial ENplus-B or Industrial grade pellets. ENplus-A2 quality pellets could be manufactured from debarked wood of these three species, since it contains ≤0.5% N, ≤0.05% S, ≤0.02% Cl and ≤1.2% ash;
  • Considering the Mo1, Mo3 and Ra4 trials as a whole, the average amounts of N, P and K removed with the aboveground woody biomass harvested during the studied period were 53.7, 10.8 and 35.0 kg ha−1 year−1, respectively. The nutrient contents in the litterfall (11.5, 1.3 and 2.2 kg ha−1 year−1, respectively) and roots (47.1, 4.6 and 15.6 kg ha−1 year−1, respectively), as well as the new amounts accumulated in the uppermost layer of the soil (0–30 cm) (57.9, −5.5, 36.7 kg ha−1 year−1, respectively), must be added to those contained in the AGWD, which amounted to 170.2, 11.2 and 89.5 kg ha−1 year−1 during the study period. Since in this study the root biomass has not been evaluated, in order to estimate it, averages value of the AGWB/root dry matter ratio of 2.03, 1.70 and 2.10 have been considered for Robinia pseudoacacia, Ulmus pumila and Populus × euroamericana “AF2”, respectively, according to the results reported by other authors for plants of similar size [68,69,70];
  • Taking into account that, on average for all the trials as a whole, the amounts of N, P and K supplied by fertilization (62.3, 27.2 and 51.7 kg ha−1 year−1, respectively) and by the irrigation water (24.4, 0.6 and 0.9 kg ha−1 year−1, respectively) were 86.7, 27.8 and 52.6 kg ha−1 year−1, respectively—for instance, in the case of N there were 83.5 kg ha−1 year−1 in the plants and in the soil layer, whose origin has not been determined.

3.3. Economic, Energy and CO2 Emission Balances of the Production Process

The technical characteristics of machinery, cultivation practices, frequencies and costs, as well as energy equivalent and materials used in the cultivation process, are shown in Table 5 and Table 9. In order to elaborate an approximation of the economic and energy costs and the CO2 emissions that the biomass production entails (Table 10), the following assumptions have been taken into account:
  • The average biomass production of Populus × euroamericana “AF2”, Robinia pseudoacacia and Ulmus pumila as a whole in the SRC field trials carried out was considered for the estimates (Table 10);
  • Although only one to two rotations have been studied in each field trial, the calculation has been extrapolated to four to five rotations, to complete a cycle of at least 15 years (Table 9). Biomass production and cultivation works from the third to the fourth–fifth rotations have been considered the same as for the second rotation;
  • The biomass produced will be harvested, chipped and transported by truck over a distance of 25 km;
  • The equivalent selling price of the chipped dry biomass, taken to the factory, has been estimated at EUR 82 t−1;
  • The valves, filters, dosing pump, irrigation pump, pressure gauges, plastic pipes (PE, PVC), installation, etc., of the irrigation system have been taken into account in the economic costs. However, only plastic pipes have been considered in the assessment of energy costs and CO2 emissions;
  • Diesel combustion releases 2.65 kg of CO2 per liter; the average C content of the harvested biomass is 50.5%; 1 tonne of C is equivalent to 3.67 tonnes of CO2; 1 kg of diesel is equivalent to 42.71 MJ, and its density is 0.85 kg L−1; the energy to CO2 conversion factor, corresponding to the Spanish electricity mix of 2019, is 241 g of CO2 released per kWh generated.
According to Table 10 and the assumptions aforementioned, the Ra4 trial improved the economic, energy and CO2 emissions ratios (that is the incomes/costs, energy generated/consumed, and CO2 fixed/released ratios) by around 2.9%, 8.1% and 3.5%, respectively, compared to Mo1; the Mo3 trial increased the economic ratio by 2.3% compared to Mo1, with similar environmental costs; in the case of the Mo2 trial (whose cultivation practices were similar to Mo1 in all aspects except in the nitrogen supplied, which was reduced by half), if it is considered that the Rp50–AF50 and Rp75–AF25 mixing treatments average the same amount of biomass as the Mo1 trial, then the Mo2 cultivation practices reduce economic costs by 4.5%, energy consumption by 12.4%, and CO2 emissions by 12.5% compared to Mo1.
In the Mo1 and Mo3 trials, the first rotation ended with negative cash flow (−570.0 and −356.4 EUR ha−1 year−1, respectively), but the investment showed profits after harvesting the second rotation (57.2 and 270.9 EUR ha−1 year−1), that is, at the 8th and 9th years after planting, respectively. However, in the Ra4 trial, economic losses persisted until the second rotation (–968.7 the 1st and −236.1 EUR ha−1 year−1 the second rotation), for which the profits were not achieved until the third rotation, in the 9th year.
Considering all the cultivation practices together, fertilizers accounted for 27–31% of the total cost, harvesting 17–21%, seedlings 17–18%, irrigation material 13–14%, transportation around 8.3%, stump removal and clearing at the end of the cycle over 4.4%, and the remaining 10% for the rest of mechanical works and water pumping. For energy consumption and CO2 emissions, the manufacture of fertilizers accounted for 35–38% of the total costs, the manufacture of irrigation plastic materials 24–26%, water pumping 21–22%, the harvest 7–9%, transportation 4–5%, and the rest of activities and materials around 4.5%.

4. Discussion

4.1. Plant Growth and Biomass Production

All the studied taxa are well adapted to the experimental growing conditions, as demonstrated by the growth rates and survival (>95%). The size reached by the plants (i.e., D, H) and the annual biomass yielded in these field trials (mean annual increments of about 12–14 t ha−1 year−1 for 3–5 years rotations, and CAI up to 15–20 t ha−1 year−1 according to Figure 1, Figure 2, Figure 3 and Figure 4) are lower than the highest mean biomass production (20–50 t ha−1 year−1) reported in the Iberian Peninsula with other woody crops, such as Eucalyptus, kenaf or Leucaena [11,12,30,59,71,72], but it is within the highest production range reported for the taxa studied [18,19,20,23,73,74]. However, under conditions of high water and nutrient availability, mean biomass yields of around 20 t ha−1 year−1 can be obtained with hybrid poplar clones [18,75,76]. It should be noted that in this study, the total aboveground dry biomass, which includes branches and bark and excludes leaves, is assessed. A larger H and D did not always imply more biomass production, since the latter is affected by wood density. For example, despite Paulownia fortunei “UHU” being one of the taxa that achieved the largest size, the biomass produced was overcome by many other taxa (Figure 1), because of its low wood density (0.282 kg dm−3) and the hole inside the stem. Bringing all this together, and the damage caused by frosts, this species is not recommended under a climate with cold winters; at the very least, it would be necessary to test the clones first before starting a commercial plantation for energy purposes. Similarly, due to their low productivity, the use of Ailanthus altissima and Platanus × hispanica is not recommended. On the contrary, the establishment of SRC commercial plantations of Populus × euroamericana, Robinia pseudoacaia and Ulmus pumila under these environmental conditions is recommended. However, in areas with mild winters, such as La Rábida, biomass production can be exceeded by evergreen species such as Eucalyptus sp. [11].
Considering the differences among taxa, on the one hand, Populus × euroamericana “AF2”, Robinia pseudoacaia and Ulmus pumila stood out in terms of biomass yield from the other deciduous species assessed (Figure 1). However, the differences between these three outstanding taxa were small in absolute value, and not significant in many cases, although they were in some others. For this reason, the discussion of results is shown together in most of the cases, differentiating between taxa only for those properties in which the differences were significant. Regarding Populus, the differences between clones make it advisable to carry out a previous analysis to select the clones best adapted to the climate and soil of the area to be planted, with special attention paid to drought tolerance [76,77]. It seems that it showed a tendency toward greater biomass production than R. pseudoacacia and U. pumila under conditions of good water availability (Figure 4 left), in addition to having a better aerial structure that facilitates the mechanization of the harvest [64,75]. Regarding Robinia pseudoacacia and Ulmus pumila, selection and improvement programs must still be developed based on survival, growth, delay in the change of phase from juvenile to mature, straightness of the main trunk, scarce branching, small thorns in the case of black locust, etc. [78,79], and, if possible, sterile hybrids that eliminate their invasive character. Ulmus pumila showed greater tolerance than Robinia pseudoacacia and Populus × euroamericana during the short water stress events suffered by the plants, showing physiological plasticity in adapting to drought by increasing leaf thickness and thus reducing water loss [80], hence its greater drought tolerance than the other two species. The climate change scenarios predict an increase in temperatures and aridity in the Mediterranean environment, which will undoubtedly affect plant growth and mortality, as well as sustainability and economic balance, since small differences in tree mortality have an important effect on land expectation value [15,81]. Therefore, the choice of the appropriate taxon and the development of genetic improvement programs are of vital importance.
On the other hand, with regard to the crop treatments assessed (monoculture or mixture of species), there was evidence that a mixed plantation that includes 50–75% of the N-fixing species may improve biomass production, as other authors have also reported for these species [19,80]. The positive effect of the mixtures 50Rp–50AF and 75Rp–25AF from our study in terms of biomass yield (7–32% higher than that of the Robinia pseudoacacia and Populus × euroamericana “AF2” monocultures) indicates the presence of competitive reduction or facilitation as mechanisms that increase biomass production under certain mixture ratios [44], and points to the importance of mixing the two species within the same row in this type of plantation [19]. However, in our study, the effect of the mixed plantation could have been slightly diminished, probably due to the N supplied with the fertilization, the N already present in the soil of the plantation site (5120 kg ha−1 in the 0–30 cm of the upper layer at the Mo2 trial), as well as the N provided by the irrigation water. This aspect should be studied in greater depth, testing different mixture ratios of R. pseudoacaia/P. × euroamericana or R. pseudoacacia/U. pumila depending on whether the growing conditions are of greater or lesser water availability, respectively.
The lower planting density in the Mo3 trial, compared to Mo1, did not increase biomass yield, but it did increase the diameters of the trunks of the plants, which would improve the quality of the biomass due to its lower proportion of bark. Under cold winter Mediterranean climates (Mo trials), rotations of 5 years (the first) followed by 4 years (the second and subsequent) could be recommended, in order to optimize biomass yield and save harvest costs. However, in areas with a climate with mild winters and a longer growing season (Ra trials), rotations can be shortened by up to 3 years without impairing biomass yield. The biomass produced by the one-year-old sprouts after cutting (R11S1) was always greater than that one year after planting (R01S1), with slight but significant differences among taxa, probably due to the planting shock, which slows growth for a few months until the root system expands and fully contacts the soil [82]. As a result, coppicing can invigorate plant growth after harvesting related to the breakdown of apical dominance and the following development of several new sprouts from dormant buds that had retained their juvenile stage. Not only the hormonal imbalance induced by cutting the shoots but also the higher proportion of earlywood vessels in the new stems might prove this response [12]. The proliferation of thinner stems after the st rotation, with higher stem number per stump, can also facilitate the mechanization of the harvest, since it largely depends on the D and H [83]. The fact that Robinia pseudoacaia and Ulmus pumila produce fruits from the 2nd or 3rd year (up to 8% of aboveground biomass) would reduce biomass production. Therefore, coppicing at the end of the rotation can improve vegetative growth for at least 1–2 years after harvesting by maintaining the juvenile stage of the plants, during which the production of reproductive organs is reduced, thus increasing the woody biomass yield in the following rotation.

4.2. Physical–Chemical Characteristics of Soil and Biomass and CO2 Sequestration

The mineral nutrient contents in plants (Table 8) are within the ranges reported by other authors for these plant species and cultivation practices, and the soils are typical farmlands with intermediate characteristics (Table 2 and Table 7), without any extreme value outside the normal range of a high-pH calcareous soil [41,80,84,85,86]. Regarding the mineral composition of the three species evaluated, it should be noted that Robinia pseudoacacia, the N-fixing species, showed higher N and lower Ca contents than the other two taxa (Ulmus pumila and Populus × euroamericana “AF2”).
Overall, for all trials and taxa as a whole, taking into account the tested experimental conditions, in which water and fertilizers were supplied, and despite the intensive production of lignocellulosic biomass, the cultivation practices made it possible to obtain a commercially demanded product, while improving soil properties. In just over 5–8 years, significant increases in parameters as important as the soil content of N, K and organic C were observed, compared with the original soil or with the uncultivated land. All this improves soil fertility, as well as the sequestration of C in stable forms in the soil, due to the cultivation system [11,12,41]. The average increase in C in the uppermost layer of the soil (0–30 cm) of this study (0.36–0.83 t ha−1 year−1 of C) is equivalent to sequestering 1.3 t ha−1 year−1 of atmospheric CO2 in the Ra4 trial, and about 3.0 t ha−1 year−1 in the Mo1, Mo2 and Mo3 trials, thus allowing the mitigation of the effects of climate change by acting as a true sink for C. The reason for the lower increase at Ra than at Mo could be due to the fact that the former has a warmer and more humid climate, which increases soil respiration, hence the C consumption and CO2 emission, in spite of presenting high biomass and litterfall production [87,88]. The capacity to sequester C obtained by SRC plantations of Robinia pseudoacacia, Ulmus pumila and Populus × euroamericana is within the average range reported for other crops—0.28–1.39 t ha−1 year−1 of C [89]. The positive effect of SRC on soil C and fertility contrasts with some negative effects of biomass production in traditional forest systems at sensitive sites [15,16,40,81]. These differences can be due to the lower level of organic matter in farmlands (<2%) with respect to forest soil and the fertilization designed to compensate for the extractions with the harvested biomass.
Additional atmospheric C has been fixed, but has not yet been mineralized into soil layers as it is in the litterfall and roots [48,86,90]. This accumulation of litterfall is close to that reported by other authors for SRC, and a good proportion of the litter C and other mineral nutrients represents a reservoir that will slowly be released into the soil [10,11,86,91], which, when added to those provided by the fine roots [92], will support the productivity of the sites. Only in litterfall, between 0.40 (Mo1) and 0.54 (Ra4) t ha−1 year−1 of C were accumulated, which is equivalent to 1.48 and 1.97 t ha−1 year−1 of CO2, respectively, although without significant differences among field trials. Likewise, the decrease in the C/N ratio to 8.6–10.4 is indicative of an adequate N release rate, and the decrease in active limestone will help to overcome the iron chlorosis observed in some taxa, as well as the decrease in the edaphic content of P [80,93,94]. The most negative aspect found in the field trials carried out has been the decrease in the soil P content (on average −5.5 kg ha−1 year−1), despite having supplied an average of 27.7 kg ha−1 year−1, and that the P contents in the harvested biomass (10.8 kg ha−1 year−1), the roots (4.6 kg ha−1 year−1) and the litter (1.3 kg ha−1 year−1) add up to 16.7 kg ha−1 year−1. Therefore, there are 16.5 kg ha−1 year−1 of P not absorbed by plants (i.e., 59.6% of the P supplied) that may have been locked by the active limestone or leached to deeper layers of the soil, out of reach of the roots. Special attention must be paid to phosphorous fertilization to ensure the sustainability of the production system. Another aspect to take into account is the slight increase in ESP, mainly at La Rábida, since this site presents values close to 15%, the threshold above which Na becomes a problem, affecting the mineral nutrition of other cations, soil dispersion and, consequently, soil crusting [95].
At Los Morales, the effects on the soil of the cultivated taxa and the mixture treatments did not differ significantly from each other in the period studied in any of the three trials (Mo1, Mo2 and Mo3). There was not even any distinction between the N-fixing species and the other species. This may be due to the fact that soil evolution processes take many years to show differences. However, the eight years for which the Mo1 and Mo2 trials lasted were enough to differentiate them from the uncultivated soil with natural herbaceous vegetation, demonstrating that a properly managed short-rotation coppice does not degrade the soil, but can improve it. However, at La Rábida, with a more humid and warmer climate than Los Morales, although the evolution of the soil under the four taxa was favorable for the 6 years studied, the effect of Robinia pseudoacacia differed from that of Ulmus pumila and Populus × euroamericana. In this case, R. pseudoacaia showed higher contents in N and P, higher ESP, and lower C/N ratio, than U. pumila and P. × euroamericana after the second rotation. The latter is consistent with other studies [80,96], and may be due to the fact that when black locust grows on nutrient-poor soils, it supplements soil nitrogen pools, increases nitrogen return in litterfall, and enhances soil nitrogen mineralization rates associated with an abundance of high nitrogen and low lignin leaf litter [96]. However, it is important to keep track of the soil’s evolution since the microorganism community could be modified [94,97].
It has been estimated that up to 76% of Robinia pseudoacaia N can come from atmospheric N2 fixation [34]. As a consequence, a large amount of inorganic N fertilizers can be substituted with the nitrogen fixed, thus reducing the high energy and economic costs related to their manufacturing [43]. Considering an average production of about 18 t ha−1 year−1 of dry biomass (AGWB + roots + leaves and fruits) and the biomass nitrogen content, around 93 kg ha−1 year−1 of N are needed to support plant growth, of which Robinia-Rhizobium root nodules can fix 71 kg ha−1 year−1. Consequently, the 71 kg of inorganic N fertilizer saved means saving about EUR 237 ha−1 year−1 (473 kg of fertilizer 15–15–15, EUR 0.5 kg−1) and 3195 MJ ha−1 year−1 during fertilizer manufacturing [43,53], and reducing the CO2 emissions into the atmosphere by about 214 kg ha−1 year−1. The role of Robinia pseudoacacia as a nitrogen-fixing species makes it an interesting species to be used in SRC energy crops, as a monoculture or mixed with other species, at least for degraded soils. However, its invasive nature suggests managing it with caution or looking for other alternative species that offer similar results. It was mentioned above that there are 83.5 kg ha−1 year−1 of N in the upper layer of the soil, the origin of which has not been verified. In the case of Robinia pseudoacacia, this could be explained by the fixation of atmospheric N. However, as other species are also evaluated in monoculture, other possible explanations are that it could come from soil layers deeper than 30 cm; N-fixing herbaceous species that grow from November to March, when the leaves of the trees have fallen; the experimental error itself, since a variation of only 0.01% in the soil N content, due to the sampling process or laboratory analysis, represents 54.5 kg ha−1 year−1 of N. As such, this will need to be assessed more precisely. In addition, the grasses that grow from November to March between the plantation lines (Figure S5) could be grazed by sheep, which would increase the production and the sustainability of the system.
The aboveground woody biomass from the three species (Robinia pseudoacacia, Ulmus pumila and Populus × euroamericana), which can be made into chips and pellets, has evidenced its appropriateness for energy use according to the international standards (Table 8), due to the energy and chemical characteristics, as well as the physical–mechanical properties, of the pellets and chips [98,99]. Its most immediate use is in thermal boilers or cogeneration plants after transformation into chips, offering a high-quality biomass since it is harvested without leaves and possesses a bark percentage of less than 15%, similar to the percentage reported for other poplar hybrid clones [100]. The resulting pellets possessed the characteristics of hardwood species [11,12,98,101,102,103], and were of a better quality than pellets from non-woody species, but did not reach the highest standardized quality of pellets derived from conifer debarked wood because of the chemical composition, which is one of the factors that affects the bonding quality of the biomass particles during pelletization [104,105]. Considering only the physical–mechanical properties of the pellets, they could provide the best quality for domestic use (ENplus-A1 commercial quality [99]). Nevertheless, the slightly high N, Cl, S and ash contents limit their suitability for domestic use to commercial ENplus-B quality, yet they are not suited for industrial purposes, according to international standards [99,106]. Therefore, this woody biomass presents high commercial quality thanks to the fact that it is harvested without leaves, the organs that possess the highest concentrations of harmful minerals, such as N, S, K and Cl, responsible for many undesirable reactions during combustion that damage furnaces and power boilers [107]. Consequently, transforming this biomass into pellets endows it with added commercial value and, taking into account the chemical properties, it would be possible to obtain pellets of the highest commercial quality (i.e., ENplus-A1) if 10–20% of this biomass is mixed with 90–80% of debarked coniferous wood (volume ratio), depending on the raw materials and the drying system. In addition, if possible, the aboveground woody biomass of Robinia pseudoacacia, Ulmus pumila and Populus × euroamericana should be debarked and used as raw material in the manufacture of pellets, in order to obtain ENplus-A2-quality pellets, or in combination with the debarked wood of species that present low enough nutrients and ash contents [102,103]. However, currently, the SRC production systems make debarking the plants difficult when they are harvested, due to their size and structure.
The determination of the rotation (harvest frequency) is of special relevance as it determines the distribution of the biomass (stems of different thickness, bark proportion, twigs) and the nutrient extraction. This last factor is relevant both in maintaining soil fertility [15,16,40,81] and in the proportion of some elements detrimental to its energy use (Cl, N, etc.) or its use in some industrial transformation processes [8,27].

4.3. Economic, Energy and CO2 Emission Balances

Overall, according to the results obtained and the proposed cultivation assumptions, with 15–17-year-long productive cycles and harvests every 3–5 years, the economic and environmental balances were shown to be positive (Table 10). The ratios greater than 1 (incomes/costs, energy generated/consumed, CO2 fixed/released) mean that the intensive cultivation of black locust, Siberian elm and hybrid poplar for energy purposes (SRC energy crop) under the Mediterranean climate can be considered sustainable over time, since it maintains or slightly improves the soil properties, while obtaining socio-economic and environmental benefits, such as the reuse of abandoned lands, rural development, the reduction of CO2 emissions and the mitigation of climate change [61,89,90,91,92,93,108]. The planning of multi-year cycles using species that sprout from the stump limits soil tillage by reducing the latter to the initial preparation of the land, and contributes to soil erosion control [59]. Despite the not too high yields, the investigated production system on degraded soils had a positive energy balance (205–237 GJ ha−1 year−1), producing 22–24 times more energy than it consumed, and in spite of the fact that the harvested biomass is burned, after a production cycle (15–17 years), 1.46–1.54 times more C remains sequestered in the production system than the production process releases into the atmosphere, in accordance with other studies [59,61,108]. With all this, the cost of the energy generated (EUR 3.1–3.2 GJ−1) is extremely competitive compared to other energy sources.
The resulting earnings (EUR 328–392 ha−1 year−1) are usually lower than those obtained from food crops when fertilizers and irrigation water are supplied, and are of the same order of magnitude as the current state subventions received by farmers. If we add to this that the initial investment does not begin to be recovered until the 8th–9th year after planting, national governments must take measures to incentivize the production of this renewable resource [58,62]. As the Mo3 trial increased the profit by 2.3% compared to Mo1, with similar environmental costs, it is possible to save some money on the initial investment by reducing the initial planting density from about 5500 to 3800 plants ha−1. Likewise, reducing the amount of inorganic N fertilizer, which can be saved in a culture of a N-fixing species mixed with another with high biomass production (Mo2 versus Mo1 trials), means reducing economic costs by at least 4.5%, and environmental costs by close to 12.5%. Therefore, in order to reduce the economic and environmental costs derived from the manufacture of N fertilizers, the inorganic fertilizer must be adjusted to the minimum amount necessary to support growth and replace the nutrients removed by the harvested biomass [11,63]. In the same way, any cultivation practice that promotes growth and reduces the supply of water and fertilizers as well as the cost of harvesting will have a major impact on the operation activities that have had the greatest weight in the economic, energy and CO2 emission balances (i.e., fertilization, irrigation, harvesting), as has also been reported by other authors [31,58,61,64]. So, to achieve higher annuities, several options could be proposed, such us the extension of the rotation cycle for 1–2 years when possible; the selection of highly efficient taxa in the use of water and nutrients; the optimization of water and nutrient supply, as well as the plantation density; the improvement of the stem structure (straightness, branching) that facilitates the harvest mechanization; the use of a mixed plantation that includes a N-fixing species; the use of a two-step harvesting system; transformation into quality chips or pellets that add value to biomass. In addition, it remains to be studied whether more than three or four cuts can be applied to the stumps, which would improve the expected earnings. As many countries have established guidelines for the sustainable harvesting of forest biomass in traditional forest systems in order to limit nutrient exports and the effect on diversity [7,16], it would be advisable to develop good practice guidelines for SRC in order to reduce the possible negative effects. Some of the aspects that should be addressed include the sites where SRC can be implanted, the distance to possible sensitive areas, the precautions to be taken with invasive species, the extension of the SRC harvest periods, the maximum area harvested per year, the correct management of fertigation, etc.

5. Conclusions

  • The biomass yield under a Mediterranean climate for short-rotation coppicing crops of Ulmus pumila, Robinia pseudoacacia and Populus × euroamericana “AF2” ranges from less than 5 to more than 20 t ha−1 year−1 (on average 12–14 t ha−1 year−1), depending on irrigation and soil quality;
  • At least 800 mm of water (rainfall + irrigation) and about 60 kg ha−1 of N (fertilization) are required annually;
  • This production system not only produces high-quality woody biomass, but, after a cycle of 15–17 years, on average it would also be able to:
    -
    Generate 205–237 GJ ha−1 year−1 net (equivalent to 57–66 MWh ha−1 year−1, or to the replacement of 5647–6528 L ha−1 year−1 of diesel);
    -
    Offer a profit of about EUR 1.5 per euro invested;
    -
    Sequester in the 0–30 cm layer of mineral soil 1.5 kg of C per kilogram released into the atmosphere.
  • After the first two rotations evaluated (5–8 years), the upper layer of the mineral soil (0–30 cm) compared to the original soil sequestered 0.36–0.83 t ha−1 year−1 of carbon (equivalent to 1.3–3.0 t ha−1 year−1 of CO2);
  • Nutrient inputs and outputs should be taken into consideration by the plantation managers in order to offset the outputs contained in the harvested biomass, and to prevent any loss in soil fertility and productivity;
  • Degraded soils could be improved if the crop is properly managed, and any cultivation practice that promotes growth, reduces irrigation, fertilization and the cost of harvesting, and increases the added value of the final marketed product will have a major impact on the economic, energy and CO2 emission balances;
  • The use of N-fixing species, such as Robinia pseudoacacia, and mixed plantations with the other two species are important aspects to be considered in future plantations.

Supplementary Materials

The following are available online at www.mdpi.com/article/10.3390/f12101337/s1, Figure S1: Schematic representation of the experimental design of the Mo1 (left) and Mo2 (right) field trials. Mo1 field trial: four hybrid clones of Populus × euroamericana (Dode) Guinier (clones “Adige” (Adg), “AF2” (AF2), “Oudenberg” (Odg), I214 (I214)); one hybrid clone of Populus × interamericana van Broekhuizen (clone “Raspalje” (Rj0); Pf, one clone of Paulownia fortunei (Seem.) Hemsl. (clone “UHU”); RpH (Robinia pseudoacacia, improved cultivar “Nyírségi”); commercial nursery plants of Populus nigra L. (Pn), Ailanthus altissima (Mill) Swingle. (Aa), Platanus × hispanica Mill. ex Münchh. (Pxh), and Ulmus pumila L. (Up). Mo2 treatments (combinations of the five different species compositions): 100% AF2 (Rp0-AF100); 25% Rp-G + 75% AF2 (Rp25-AF75); 50% Rp-G + 50% AF2 (Rp50-AF50); 75% Rp-G + 25% AF2 (Rp75-AF25); 100% Rp-G (Rp100-AF0). Rp-G (Robinia pseudoacacia, commercial nursery plants from Germany), Figure S2: Schematic representation of the experimental design of the Mo3 field trial. AF2 (Populus × euroamericana “AF2”), Rp-H and Rp-G (Robinia pseudoacacia cultivar “Nyírségi” and German provenance, respectively); Up (Ulmus pumila), Figure S3: Schematic representation of the experimental design of the Ra4 field trial. AF2 (Populus × euroamericana “AF2”), Rp-H and Rp-G (Robinia pseudoacacia cultivar “Nyírségi” and German provenance, respectively); Up (Ulmus pumila), Figure S4: Current annual increment (CAI) and mean annual increment (MAI) of the woody biomass dry weight of the aboveground part of the plant (AGWB) throughout the period studied for the 1st (left) and 2nd (right) rotations and the three taxa. Data from all field trials are shown as a whole. Taxa: Populus AF2 (Populus × euroamericana “AF2”), Robinia pseudoacacia (including cultivar “Nyírségi” and German provenance); Ulmus pumila, Figure S5: (a) 3.5-year-old Populus × euroamericana “AF2”, 1st rotation; (b) 3.5-year-old Ulmus pumila, 1st rotation—the litter layer can be seen on the soil surface; (c) 2.0-year-old Robinia pseudoacacia, 1st rotation—it can be seen the grass grown during the cold season, which can be grazed by sheep; (d) 2.0-year-old mixed plantation of Populus × euroamericana “AF2” and Robinia pseudoacacia, 1st rotation; (e, f, g) 3.0-year-old sprouts of Populus × euroamericana “AF2”, Ulmus pumila and Robinia pseudoacacia, respectively; (h) litterfall of Robinia pseudoacacia; (i, j) wood pellets from Ulmus pumila and Robinia pseudoacacia, respectively, Table S1: Physical–chemical characteristics of the upper soil layer at the beginning of the experiments (mean (SE)). Mo (Los Morales), Ra (La Rábida), Table S2: Relationships between stem diameter measured at 10 cm above ground (D, mm) and plant height (H, cm). H = a·D2 + b·D + c·D ranged from 6 to 170 mm. r = Pearson coefficient. Significance level, p < 0.001 for all cases. n = sample size, Table S3: Standard technical methods and instruments used to assess the physical–chemical properties of plant material, soil and litterfall, Table S4: Morphological characteristics (average (SE)) of the main stems emerging from the stumps at the end of the 1st (R04S4 plants) and the 2nd (R18S4 plants) rotation for the Mo1 assay. Dmax, maximum stem diameter among those emerging from a stump; NSt, number of stems; BA, basal area. The four variables showed significant differences between taxa and between rotations, p < 0.001.

Author Contributions

Conceptualization, M.F., J.A. and R.T.; methodology, S.P.A., M.F. and R.T.; software, S.P.A.; validation, M.F., R.T. and J.A.; formal analysis, S.P.A., R.T.; investigation, S.P.A.; resources, R.T., J.A.; data curation, S.P.A. and M.F.; writing—original draft preparation, S.P.A.; writing—review and editing, M.F., J.A. and R.T.; visualization, S.P.A., M.F. and J.A.; supervision, R.T. and M.F.; project administration, M.F.; funding acquisition, M.F., R.T. and J.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Science and Innovation Ministry (ref. AGL2010-16575) and the Economy and Competitiveness Ministry of Spain (ref. CTQ2013-46804-C2-1-R and CTQ2017-85251-C2-2-R), by FEDER funds of the EU, and by the company ENCE, energía y celulosa S.A. (8%, 6%, 6%, 70%, 10%, respectively).

Data Availability Statement

Not Applicable.

Acknowledgments

The authors want to thank the Diputación de Granada (Spain) for the donation of farmland for the establishment of the experimental plots in Los Morales, and for the help in the cultivation works; the company Tubocás S.L. for its contribution to the sampling and transport of biomass, as well as to the harvest at the end of each crop rotation; and Biopoplar Ibérica S.L. for the provision of some plant taxa.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. United Nations. In Proceedings of the Conference of the Parties on Its Twenty-First Session, Paris, France, 30 November–13 December 2015. Available online: http://unfccc.int/resource/docs/2015/cop21/eng/10a01.pdf (accessed on 29 July 2021).
  2. European Union. Directive (EU) 2018/2001 of the European Parliament and of the Council of 11 December 2018 on the Use of Energy from Renewable Sources. Available online: https://eur-lex.europa.eu/eli/dir/2018/2001/oj (accessed on 25 July 2021).
  3. REN21. Renewable. Global Status Report, Renewable Energy Policy Network for the 21st Century, France. Available online: https://www.ren21.net/wp-content/uploads/2019/05/GSR2017_Full-Report_English.pdf (accessed on 30 July 2021).
  4. European Commission. Communication from The Commission to The European Parliament, The Council, The European Economic and Social Committee and The Committee of The Regions Energy Roadmap 2050/* COM/2011/0885 final */. Available online: https://eur-lex.europa.eu/legal-content/EN/ALL/?uri=celex:52011DC0885 (accessed on 27 June 2021).
  5. International Energy Agency (IEA). World Energy Outlook. Available online: https://iea.blob.core.windows.net/assets/4a50d774-5e8c-457e-bcc9513357f9b2fb/World_Energy_Outlook_2017.pdf (accessed on 30 July 2021).
  6. Antczak, A.; Swierkosz, R.; Szeniawski, M.; Marchwicka, M.; Akus-Szylberg, F.; Przybysz, P.; Zawadzki, J. The comparison of acid and enzymatic hydrolysis of pulp obtained from poplar wood (Populus sp.) by the Kraft method. Drewno 2019, 63, 1–14. [Google Scholar] [CrossRef]
  7. Titus, B.D.; Brown, K.; Helmisaari, H.-S.; Vanguelova, E.; Stupak, I.; Evans, A.; Clarke, N.; Guidi, C.; Bruckman, V.J.; Varnagiryte-Kabasinskiene, I.; et al. Sustainable forest biomass: A review of current residue harvesting guidelines. Energy Sustain. Soc. 2021, 11, 1–32. [Google Scholar] [CrossRef]
  8. Małachowska, E.; Lipkiewicz, A.; Niemczyk, M.; Dubowik, M.; Boruszewski, P.; Przybysz, P. Influences of fiber and pulp properties on papermaking ability of cellulosic pulps produced from alternative fibrous raw materials. J. Nat. Fibers 2019, 18, 1–11. [Google Scholar] [CrossRef]
  9. BIC. Strategic innovation and research agenda (SIRA) Bio-based and renewable industries for development and growth in Europe. In A Public-Private Partnership on Bio-Based Industries; Biobased Industries Consortium, Ed.; Biobased Industries Consortium: Brussels, Belgium, 2013; Available online: https://ec.europa.eu/research/participants/data/ref/h2020/other/legal_basis/jtis/bbi/bbi-sira_en.pdf (accessed on 30 July 2021).
  10. Bogdanski, A.; Dubois, O.; Jamieson, C.; Krell, R. Making Integrated Food-Energy Systems Work for People and Climate; FAO: Rome, Italy, 2010; p. 116. Available online: http://www.fao.org/docrep/013/i2044e/i2044e.pdf (accessed on 27 June 2021).
  11. Fernández, M.; Alaejos, J.; Andivia, E.; Vázquez-Piqué, J.; Ruiz, F.; López, F.; Tapias, R. Eucalyptus x urograndis biomass production for energy purposes exposed to a Mediterranean climate under different irrigation and fertilisation regimes. Biomass Bioenergy 2018, 111, 22–30. [Google Scholar] [CrossRef]
  12. Fernández, M.; Alaejos, J.; Andivia, E.; Madejón, P.; Díaz, M.; Tapias, R. Short rotation coppice of leguminous tree Leucaena spp. improves soil fertility while producing high biomass yields in Mediterranean environment. Ind. Crop Prod. 2020, 157, 112911. [Google Scholar] [CrossRef]
  13. Riffell, S.; Verschuyl, J.; Miller, D.; Wigley, T.B. Biofuel harvests, coarse woody debris, and biodiversity—A meta-analysis. For. Ecol. Manag. 2011, 261, 878–887. [Google Scholar] [CrossRef]
  14. Bouget, C.; Lassauce, A.; Jonsell, M. Effects of fuelwood harvesting on biodiversity—A review focused on the situation in Europe. Can. J. For. Res. 2012, 42, 1421–1432. [Google Scholar] [CrossRef]
  15. Miettinen, J.; Ollikainen, M.; Nieminen, T.M.; Ukonmaanaho, L.; Laurén, A.; Hynynen, J.; Lehtonen, M.; Valsta, L. Whole-tree harvesting with stump removal versus stem-only harvesting in peatlands when water quality, biodiversity conservation and climate change mitigation matter. For. Policy Econ. 2013, 47, 25–35. [Google Scholar] [CrossRef]
  16. Bessaad, A.; Bilger, I.; Korboulewsky, N. Assessing Biomass Removal and Woody Debris in Whole-Tree Harvesting System: Are the Recommended Levels of Residues Ensured? Forests 2021, 12, 807. [Google Scholar] [CrossRef]
  17. Srirangan, K.; Akawi, L.; Moo-Young, M.; Chou, C.P. Towards sustainable production of clean energy carriers from biomass resources. Appl. Energy 2012, 100, 172–186. [Google Scholar] [CrossRef]
  18. Sixto, H.; Cañellas, I.; van Arendonk, J.; Ciria, P.; Camps, F.; Sánchez, M.; Sánchez-González, M. Growth potential of different species and genotypes for biomass production in short rotation in Mediterranean environments. For. Ecol. Manag. 2015, 354, 291–299. [Google Scholar] [CrossRef]
  19. Oliveira, N.; del Rio, M.; Forrester, D.; Rodríguez-Soalleiro, R.; Pérez-Cruzado, C.; Canellas, I.; Sixto, H. Mixed short rotation plantations of Populus alba and Robinia pseudoacacia for biomass yield. For. Ecol. Manag. 2018, 410, 48–55. [Google Scholar] [CrossRef]
  20. Seserman, D.-M.; Pohle, I.; Veste, M.; Freese, D. Simulating Climate Change Impacts on Hybrid-Poplar and Black Locust Short Rotation Coppices. Forests 2018, 9, 419. [Google Scholar] [CrossRef] [Green Version]
  21. Szabolcs, K.; Varga, D. Physical and mechanical properties of wood from invasive tree species. Maderas Cienc. Tecnol. 2020, 23. [Google Scholar] [CrossRef]
  22. Parrotta, J. Productivity, nutrient cycling, and succession in single-and mixed-species plantations of Casuarina equisetifolia, Eucalyptus robusta, and Leucaena leucocephala in Puerto Rico. For. Ecol. Manag. 1999, 124, 45–77. [Google Scholar] [CrossRef]
  23. Durán Zuazo, V.H.; Jiménez Bocanegra, J.A.; Perea Torres, F.; Rodríguez Pleguezuelo, C.R.; Francia Martínez, J.R. Biomass yield potential of Paulownia trees in a semi-arid Mediterranean environment (S Spain). Int. J. Renew. Energy Res. 2013, 3, 789–793. Available online: https://ijrer.org/ijrer/index.php/ijrer/article/view/844 (accessed on 30 August 2021).
  24. Lauri, P.; Havlík, P.; Kindermann, G.; Forsell, N.; Böttcher, H.; Obersteiner, M. Woody biomass energy potential in 2050. Energy Policy 2014, 66, 19–31. [Google Scholar] [CrossRef]
  25. Fischer, G.; Prieler, S.; van Velthuizen, H.; Berndes, G.; Faaij, A.; Londo, M.; de Wit, M. Biofuel production potentials in Europe: Sustainable use of cultivated land and pastures, Part II: Land use scenarios. Biomass Bioenergy 2010, 34, 173–187. [Google Scholar] [CrossRef] [Green Version]
  26. Food and Agriculture Organization of United Nations (FAO). Combating land degradation for food security and provision of soil ecosystem services in Europe and Central Asia. In Proceedings of the 39th Session of the European Commission on Agriculture, Budapest, Hungary, 22–23 September 2015; Available online: http://www.fao.org/europe/commissions/eca/eca-39/en/ (accessed on 27 June 2021).
  27. Zielenkiewicz, T.; Szadkowski, J.; Drożdżek, M.; Zielenkiewicz, A.; Kłosińska, T.; Antczak, A.; Zawadzki, J.; Gawron, J. Application of X-ray fluorescence technique for determination of heavy metals uptake by different species of poplar. Drewno 2016, 59, 113–126. [Google Scholar] [CrossRef]
  28. Geyer, W. Biomass production in the Central Great Plains USA under various coppice regimes. Biomass Bioenergy 2006, 30, 778–783. [Google Scholar] [CrossRef]
  29. O’Neill, M.K.; Shock, C.C.; Lombard, K.A.; Heyduck, R.F.; Feibert, E.B.G.; Smeal, D.; Arnold, R.N. Hybrid poplar (Populus ssp.) selections for arid and semi-arid intermountain regions of the western United States. Agrofor. Syst. 2010, 79, 409–418. [Google Scholar] [CrossRef]
  30. Pérez, I.; Pérez, J.; Carrasco, J.; Ciria, P. Siberian elm responses to different culture conditions under short rotation forestry in Mediterranean areas. Turk. J. Agric. For. 2014, 38, 652–662. [Google Scholar] [CrossRef]
  31. Vanbeveren, S.P.; De Francesco, F.; Ceulemans, R.; Spinelli, R. Productivity of mechanized whip harvesting with the Stemster MkIII in a short-rotation coppice established on farmland. Biomass Bioenergy 2018, 108, 323–329. [Google Scholar] [CrossRef]
  32. García, I.P. Evaluación de Ulmus pumila L. y Populus spp. Como Cultivos Energéticos en Corta Rotación. Ph.D. Thesis, Polytechnic University of Madrid, Madrid, Spain, 2016. [Google Scholar] [CrossRef]
  33. Sanz Elorza, M.; Dana Sánchez, E.D.; Sobrino Vesperinas, E. Atlas de las plantas alóctonas invasoras en España. In Dirección General Para la Biodiversidad; Ministerio de Medio Ambiente: Madrid, Spain, 2004; p. 378. Available online: https://www.miteco.gob.es/es/biodiversidad/temas/inventarios-nacionales/inventario-especies-terrestres/inventario-nacional-de-biodiversidad/ieet_flora_vasc_aloct_invas.aspx (accessed on 18 August 2021).
  34. Marron, N.; Gana, C.; Gerant, D.; Maillard, P.; Priault, P.; Epron, D. Estimating symbiotic N2 fixation in Robinia pseudoacacia. J. Plant Nutr. Soil Sci. 2018, 181, 296–304. [Google Scholar] [CrossRef] [Green Version]
  35. Rédei, K.; Osváth-Bujtás, Z.; Veperdi, I. Black locust (Robinia pseudoacacia L.). Improvement in Hungary: A review. Acta Silv. Lign. Hung. 2008, 4, 127–132. [Google Scholar]
  36. Rédei, K.; Keserű, Z.; Csiha, I.; Rásó, J.; Végh, K.; Antal, B. Juvenile Growth and Morphological Traits of Micropropagated Black Locust (Robinia pseudoacacia, L.) Clones under Arid Site Conditions. Acta Silv. Lignaria Hung. 2013, 9, 35–42. [Google Scholar] [CrossRef] [Green Version]
  37. Dini-Papanastasi, O. Effects of clonal selection on biomass production and quality in Robinia pseudoacacia var. monophylla Carr. For. Ecol. Manag. 2008, 256, 849–854. [Google Scholar] [CrossRef]
  38. Kreiser, M.; Giblin, C.; Murphy, R.; Fiesel, P.; Braun, L.; Johnson, G.; Wyse, D.; Cohen, J.D. Conversion of Indole-3-Butyric Acid to Indole-3-Acetic Acid in Shoot Tissue of Hazelnut (Corylus) and Elm (Ulmus). J. Plant Growth Regul. 2016, 35, 710–721. [Google Scholar] [CrossRef]
  39. Vanguelova, E.; Pitman, R. Impacts of short rotation forestry on soil sustainability. In Short Rotation Forestry: Review of Growth and Environmental Impacts, Forest Research Monograph, Forest Research; McKay, H., Ed.; Forest Research Alice Holt Lodge: Surrey, UK, 2011; p. 37. ISBN 978-0-85538-827-0. [Google Scholar]
  40. Augusto, L.; Achat, D.L.; Bakker, M.R.; Bernier, F.; Bert, D.; Danjon, F.; Khlifa, R.; Meredieu, C.; Trichet, P. Biomass and nutrients in tree root systems-sustainable harvesting of an intensively managed Pinus pinaster (Ait.) planted forest. GCB Bioenergy 2014, 7, 231–243. [Google Scholar] [CrossRef]
  41. Madejón, P.; Alaejos, J.; García-Álbala, J.; Fernández, M.; Madejón, E. Three-year study of fast-growing trees in degraded soils amended with composts: Effects on soil fertility and productivity. J. Environ. Manag. 2016, 169, 18–26. [Google Scholar] [CrossRef]
  42. Grünewald, H.; Böhm, C.; Quinkenstein, A.; Grundmann, P.; Eberts, J.; Von Wühlisch, G. Robinia pseudoacacia L.: A Lesser known tree species for biomass production. BioEnergy Res. 2009, 2, 123–133. [Google Scholar] [CrossRef]
  43. Giuntoli, J.; Agostini, A.; Edwards, R.; Marelli, L. Solid and gaseous bioenergy pathways: Input values and GHG emissions. In Calculated According to the Methodology Set in COM (2016) Version 2, EUR 27215 EN; Publications Office of the European Union: Luxembourg, 2017. [Google Scholar] [CrossRef]
  44. Forrester, D.I.; Bauhus, J.; Cowie, A.; Vanclay, J. Mixed-species plantations of Eucalyptus with nitrogen-fixing trees: A review. For. Ecol. Manag. 2006, 233, 211–230. [Google Scholar] [CrossRef] [Green Version]
  45. Bouillet, J.; Laclau, J.; Gonçalves, J.L.M.; Moreira, M.Z.; Trivelin, P.; Jourdan, C.; Silva, E.; Piccolo, M.D.C.; Tsai, S.; Galiana, A. Mixed-species plantations of Acacia mangium and Eucalyptus grandis in Brazil. For. Ecol. Manag. 2008, 255, 3918–3930. [Google Scholar] [CrossRef]
  46. Moshki, A.; Lamersdorf, N.P. Symbiotic nitrogen fixation in black locust (Robinia pseudoacacia L.) seedlings from four seed sources. J. For. Res. 2011, 22, 689–692. [Google Scholar] [CrossRef]
  47. Verlinden, M.; Broeckx, L.; Bulcke, J.V.D.; Van Acker, J.; Ceulemans, R. Comparative study of biomass determinants of 12 poplar (Populus) genotypes in a high-density short-rotation culture. For. Ecol. Manag. 2013, 307, 101–111. [Google Scholar] [CrossRef]
  48. Magar, L.B.; Khadka, S.; Joshi, J.R.R.; Pokharel, U.; Rana, N.; Thapa, P.; Sharma, K.R.S.R.; Khadka, U.; Marasini, B.P.; Parajuli, N. Total Biomass Carbon Sequestration Ability Under the Changing Climatic Condition by Paulownia tomentosa Steud. Int. J. Appl. Sci. Biotechnol. 2018, 6, 220–226. [Google Scholar] [CrossRef]
  49. Nilsen, E.T.; Huebner, C.D.; Carr, D.E.; Bao, Z. Interaction between Ailanthus altissima and Native Robinia pseudoacacia in Early Succession: Implications for Forest Management. Forests 2018, 9, 221. [Google Scholar] [CrossRef] [Green Version]
  50. UNE. Solid Biofuels-Fuel Specifications and Classes-Part 2: Graded Wood Pellets (ISO 17225-2:2014); UNE-EN ISO 17225-2; ISO: Geneva, Switzerland, 2014; Available online: https://www.une.org/encuentra-tu-norma/busca-tu-norma/norma/?c=N0053748 (accessed on 8 August 2021).
  51. TRAGSA (Empresa de Transformación Agraria, S.A.). Spain, Tarifas TRAGSA. Available online: https://www.tragsa.es/es/grupo-tragsa/regimen-juridico/tarifas/Paginas/default.aspx (accessed on 27 July 2021).
  52. ENCE (ENCE, Energía y Celulosa, S.A.). Spain, Tarifa ENCE. Available online: https://compramosmadera.es/images/tarifasPDF/20210301_TARIFAS_MADERA_WEB_Mar2021_v1.pdf (accessed on 27 July 2021).
  53. Audsley, E. Harmonisation of environmental life cycle assessment for agriculture, final report. In Concerted Action AIR3-CT94; European Commission, Directorate-General VI Agriculture: Brussels, Belgium, 2003; p. 129. Available online: https://cordis.europa.eu/project/id/AIR32028 (accessed on 27 June 2021).
  54. Boto Fidalgo, J.A.; Pastrana Santamarta, P.; Suárez de Cepeda Martínez, M. Consumos Energéticos en Las Operaciones Agrícolas en España; IDAE, Instituto Para la Diversificación y Ahorro de la Energía, Ministerio de Industria, Turismo y Comercio: Madrid, Spain, 2005; p. 74. Available online: https://www.idae.es/uploads/documentos/documentos_10255_Consumos_energeticos_operaciones_agricolas_Espana_05_d94c1676.pdf (accessed on 30 July 2021).
  55. López Martínez, J.M.; Sánchez Alejo, J. Energy Consumption and Emissions Associated with Transportation by Car and Truck. Monografías EnerTrans. 2009. n° 12. Available online: https://www.investigacion-ffe.es/documentos/enertrans/EnerTrans_Transporte_coche_y_cami%C3%B3n.pdf (accessed on 30 July 2021).
  56. López Martínez, J.M.; Sánchez Alejo, J.; Mora Sotomayor, A. Energy Consumption and Emissions Associated with Construction and Maintenance of Vehicles. Monografías EnerTrans. 2009. n° 11. Available online: https://www.investigacion-ffe.es/documentos/enertrans/EnerTrans_Construcci%C3%B3n_Vehiculos.pdf (accessed on 30 July 2021).
  57. Camacho, E.; Rodríguez, J.A.; Montesinos, P.; Carrillo, T. Ahorro de energía en el riego. Agricultura 2010, 930, 454–457. Available online: https://www.mapa.gob.es/ministerio/pags/biblioteca/revistas/pdf_Agri%5CAgri_2010_930_completa.pdf (accessed on 20 August 2021).
  58. Manzone, M.; Airoldi, G.; Balsari, P. Energetic and economic evaluation of a poplar cultivation for the biomass production in Italy. Biomass Bioenergy 2009, 33, 1258–1264. [Google Scholar] [CrossRef]
  59. Di Nasso, O.N.N.; Guidi, W.; Ragaglini, G.; Tozzini, C.; Bonari, E. Biomass production and energy balance of a 12-year-old short-rotation coppice poplar stand under different cutting cycles. GCB Bioenergy 2010, 2, 89–97. [Google Scholar] [CrossRef]
  60. De Wit, M.; Faaij, A. European biomass resource potential and costs. Biomass Bioenergy 2010, 34, 188–202. [Google Scholar] [CrossRef]
  61. Dillen, S.; Djomo, S.; Al Afas, N.; Vanbeveren, S.; Ceulemans, R.J. Biomass yield and energy balance of a short-rotation poplar coppice with multiple clones on degraded land during 16 years. Biomass Bioenergy 2013, 56, 157–165. [Google Scholar] [CrossRef] [Green Version]
  62. Schweier, J.; Becker, G. Economics of poplar short rotation coppice plantations on marginal land in Germany. Biomass Bioenergy 2013, 59, 494–502. [Google Scholar] [CrossRef]
  63. Miguel, G.S.; Corona, B.; Ruiz, D.; Landholm, D.; Laina, R.; Tolosana, E.; Sixto, H.; Canellas, I. Environmental, energy and economic analysis of a biomass supply chain based on a poplar short rotation coppice in Spain. J. Clean. Prod. 2015, 94, 93–101. [Google Scholar] [CrossRef]
  64. Vanbeveren, S.P.; Spinelli, R.; Eisenbies, M.; Schweier, J.; Mola-Yudego, B.; Magagnotti, N.; Acuna, M.; Dimitriou, I.; Ceulemans, R. Mechanised harvesting of short-rotation coppices. Renew. Sustain. Energy Rev. 2017, 76, 90–104. [Google Scholar] [CrossRef]
  65. McFarquhar, M. Modeling Group-level repeated measurements of neuroimaging data using the univariate general linear model. Front. Neurosci. 2019, 13, 352. [Google Scholar] [CrossRef]
  66. Römheld, V. Diagnosis of deficiency and toxicity of nutrients. In Mineral Nutrition of Higher Plants, 3rd ed.; Marschner, P., Ed.; Elsevier Ltd.: Amsterdam, The Netherlands, 2012; pp. 299–312. [Google Scholar] [CrossRef]
  67. UNE. UNE-EN ISO 17225-4: Solid Biofuels-Fuel Specifications and Classes-Part 4: Graded Wood Chips (ISO 17225-4:2014). Available online: https://www.une.org/encuentra-tu-norma/busca-tu-norma/norma/?c=N0053750 (accessed on 8 August 2021).
  68. Montero, G.; Ruiz-Peinado, R.; Muñoz, M. Producción de biomasa y fijación de CO2 por los bosques españoles. In Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria; Ministerio de Educación y Ciencia: Madrid, Spain, 2005; p. 270. Available online: http://libros.inia.es/libros/product_info.php?products_id=64 (accessed on 10 August 2021).
  69. Mokany, K.; Raison, R.J.; Prokushkin, A.S. Critical analysis of root: Shoot ratios in terrestrial biomes. Glob. Chang. Biol. 2005, 12, 84–96. [Google Scholar] [CrossRef]
  70. Luo, Y.; Wang, X.; Zhang, X.; Booth, T.H.; Lu, F. Root:shoot ratios across China’s forests: Forest type and climatic effects. For. Ecol. Manag. 2012, 269, 19–25. [Google Scholar] [CrossRef]
  71. Pérez, S.; Renedo, C.; Ortiz, A.; Mañana, M.; Delgado, F.; Tejedor, C. Energetic density of different forest species of energy crops in Cantabria (Spain). Biomass Bioenergy 2011, 35, 4657–4664. [Google Scholar] [CrossRef]
  72. Pérez-Cruzado, C.; Merino, A.; Rodríguez-Soalleiro, R. A management tool for estimating bioenergy production and carbon sequestration in Eucalyptus globulus and Eucalyptus nitens grown as short rotation woody crops in north-west Spain. Biomass Bioenergy 2011, 35, 2839–2851. [Google Scholar] [CrossRef]
  73. Stolarski, M.J.; Krzyżaniak, M.; Szczukowski, S.; Tworkowski, J.; Bieniek, A. Dendromass derived from agricultural land as energy feedstock. Pol. J. Environ. Stud. 2013, 22, 511–520. [Google Scholar]
  74. Berdón, J.B.; Calvo, A.J.M.; Barroso, L.R.; Alcobendas, A.I.P.; Cortés, J.G. Study of Paulownia’s Biomass Production in Mérida (Badajoz), Southwestern Spain. Environ. Ecol. Res. 2017, 5, 521–527. [Google Scholar] [CrossRef] [Green Version]
  75. Paris, P.; Mareschi, L.; Sabatti, M.; Pisanelli, A.; Ecosse, A.; Nardin, F.; Scarascia-Mugnozza, G. Comparing hybrid Populus clones for SRF across northern Italy after two biennial rotations: Survival, growth and yield. Biomass Bioenergy 2011, 35, 1524–1532. [Google Scholar] [CrossRef]
  76. Navarro, A.; Portillo-Estrada, M.; Arriga, N.; Vanbeveren, S.P.P.; Ceulemans, R. Genotypic variation in transpiration of coppiced poplar during the third rotation of a short-rotation bio-energy culture. GCB Bioenergy 2018, 10, 592–607. [Google Scholar] [CrossRef] [PubMed]
  77. Bloemen, J.; Fichot, R.; Horemans, J.A.; Broeckx, L.S.; Verlinden, M.; Zenone, T.; Ceulemans, R.J. Water use of a multigenotype poplar short-rotation coppice from tree to stand scale. GCB Bioenergy 2016, 9, 370–384. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Bautista, I.; Mauri, P.V.; Plaza, A.; Ruiz-Fernández, J.; Amorós, M.C.; Sanz, M.; Curt, M.D. Vegetative propagation of Ulmus pumila L. by stem cuttings with a view to the development of bred lines for woody biomass plantations. In Proceedings of the 25th European Biomass Conference and Exhibition, Session 1DV.2.6, Stockholm, Sweden, 12–15 June 2017; ETA-Florence Renewable Energies: Firenze, Italy, 2017; pp. 253–257. Available online: http://www.etaflorence.it/proceedings/?detail=13655 (accessed on 8 August 2021). [CrossRef]
  79. Amorós, M.D.L.C.; Mauri, P.V.; Curt, M.D. The influence of tree management practices on phenological growth stages of Ulmus pumila L. (Siberian elm). Ann. Appl. Biol. 2021, 179, 1–14. [Google Scholar] [CrossRef]
  80. Qin, J.; Xi, W.; Rahmlow, A.; Kong, H.; Zhang, Z.; Shangguan, Z. Effects of forest plantation types on leaf traits of Ulmus pumila and Robinia pseudoacacia on the Loess Plateau, China. Ecol. Eng. 2016, 97, 416–425. [Google Scholar] [CrossRef]
  81. Bessaad, A.; Terreaux, J.-P.; Korboulewsky, N. Assessing the land expectation value of even-aged vs coppice-with-standards stand management and long-term effects of whole-tree harvesting on forest productivity and profitability. Ann. For. Sci. 2021, 78, 1–15. [Google Scholar] [CrossRef]
  82. Parrotta, J.A. Leucaena leucocephala (Lam.) de Wit: Leucaena, tantan. Res. Note SO-ITFSM-USDA Forest Service; Southern Forest Experiment Estation: New Orleans, LA, USA, 1992; p. 8. Available online: https://www.fs.usda.gov/treesearch/pubs/50724 (accessed on 8 August 2021).
  83. Vanbeveren, S.; Schweier, J.; Berhongaray, G.; Ceulemans, R.J. Operational short rotation woody crop plantations: Manual or mechanised harvesting? Biomass Bioenergy 2015, 72, 8–18. [Google Scholar] [CrossRef]
  84. Jug, A.; Hofmann-Schielle, C.; Makeschin, F.; Rehfuess, K. Short-rotation plantations of balsam poplars, aspen and willows on former arable land in the Federal Republic of Germany. II. Nutritional status and bioelement export by harvested shoot axes. For. Ecol. Manag. 1999, 121, 67–83. [Google Scholar] [CrossRef]
  85. Romanyà, J.; Rovira, P.; Vallejo, R. Análisis del carbono en los suelos agrícolas de España. Aspectos relevantes en relación a la reconversión a la agricultura ecológica en el ámbito mediterráneo. Ecosistemas 2007, 16, 50–57. Available online: https://www.revistaecosistemas.net/index.php/ecosistemas/article/view/138 (accessed on 27 August 2021).
  86. González, I.; Sixto, H.; Rodríguez-Soalleiro, R.; Oliveira, N. Nutrient Contribution of Litterfall in a Short Rotation Plantation of Pure or Mixed Plots of Populus alba L. and Robinia pseudoacacia L. Forests 2020, 11, 1133. [Google Scholar] [CrossRef]
  87. Buerkert, A.; Joergensen, R.G.; Ludwig, B.; Schlecht, E. Nutrient and Carbon Fluxes in Terrestrial Agro-Ecosystems. In Mineral Nutrition of Higher Plants, 3rd ed.; Marschner, P., Ed.; Elsevier Ltd.: Amsterdam, The Netherlands, 2012; pp. 473–482. [Google Scholar] [CrossRef]
  88. Ferré, C.; Comolli, R. Comparison of soil CO2 emissions between short-rotation coppice poplar stands and arable lands. Iforest 2018, 11, 199–205. [Google Scholar] [CrossRef] [Green Version]
  89. Campos, R.; Pires, G.F.; Costa, M.H. Soil Carbon Sequestration in Rainfed and Irrigated Production Systems in a New Brazilian Agricultural Frontier. Agriculture 2020, 10, 156. [Google Scholar] [CrossRef]
  90. Berhongaray, G.; Verlinden, M.; Broeckx, L.S.; Janssens, I.; Ceulemans, R. Soil carbon and belowground carbon balance of a short-rotation coppice: Assessments from three different approaches. GCB Bioenergy 2016, 9, 299–313. [Google Scholar] [CrossRef] [Green Version]
  91. Andivia, E.; Rolo, V.; Jonard, M.; Formánek, P.; Ponette, Q. Tree species identity mediates mechanisms of top soil carbon sequestration in a Norway spruce and European beech mixed forest. Ann. For. Sci. 2016, 73, 437–447. [Google Scholar] [CrossRef] [Green Version]
  92. Jha, P.; Mohapatra, K.P. Leaf litterfall, fine root production and turnover in four major tree species of the semi-arid region of India. Plant Soil 2009, 326, 481–491. [Google Scholar] [CrossRef]
  93. George, E.; Horst, W.J.; Neumann, E. Adaptation of Plants to Adverse Chemical Soil Conditions. In Mineral Nutrition of Higher Plants, 3rd ed.; Marschner, P., Ed.; Elsevier Ltd.: Amsterdam, The Netherlands, 2012; pp. 409–472. [Google Scholar] [CrossRef]
  94. Li, K.; Han, X.; Ni, R.; Shi, G.; De-Miguel, S.; Li, C.; Shen, W.; Zhang, Y.; Zhang, X. Impact of Robinia pseudoacacia stand conversion on soil properties and bacterial community composition in Mount Tai, China. For. Ecosyst. 2021, 8, 1–12. [Google Scholar] [CrossRef]
  95. Laker, M.C.; Nortjé, G.P. Review of existing knowledge on soil crusting in South Africa. Adv. Agron. 2019, 115, 189–242. [Google Scholar] [CrossRef]
  96. Rice, S.K.; Westerman, B.; Federici, R. Impacts of the exotic, nitrogen-fixing black locust (Robinia pseudoacacia) on nitrogen cycling in a pine-oak ecosystem. Plant Ecol. 2004, 174, 97–107. [Google Scholar] [CrossRef]
  97. Castro-Díez, P.; Fierro-Brunnenmeister, N.; González-Muñoz, N.; Gallardo, A. Effects of exotic and native tree leaf litter on soil properties of two contrasting sites in the Iberian Peninsula. Plant Soil 2011, 350, 179–191. [Google Scholar] [CrossRef]
  98. Carmona, R.; Nuñez, T.; Alonso, M. Biomass yield and quality of an energy dedicated crop of poplar (Populus spp.) clones in the Mediterranean zone of Chile. Biomass Bioenergy 2015, 74, 96–102. [Google Scholar] [CrossRef]
  99. European Pellet Council (EPC). Available online: https://epc.bioenergyeurope.org/about-pellets/pellets-basics/pellet-standards (accessed on 27 August 2021).
  100. Morhart, C.; Sheppard, J.; Spiecker, H. Above Ground Leafless Woody Biomass and Nutrient Content within Different Compartments of a P. maximowicii × P. trichocarpa Poplar Clone. Forests 2013, 4, 471–487. [Google Scholar] [CrossRef] [Green Version]
  101. Fernández, M.; García-Albalá, J.; Andivia, E.; Alaejos, J.; Tapias, R.; Menéndez, J.; Martín, R.T. Sickle bush (Dichrostachys cinerea L.) field performance and physical–chemical property assessment for energy purposes. Biomass Bioenergy 2015, 81, 483–489. [Google Scholar] [CrossRef]
  102. Vasiliki, K.; Ioannis, B. Utilization perspectives of wood and bark of the invasive species of Ailanthus and Acacia in the production of pellets. In Proceedings of the International Forestry and Environment Symposium, Trabzon, Turkey, 7–10 November 2017; Available online: http://users.auth.gr/jbarb/Publications/IFES2017-209_Utilization.pdf (accessed on 8 August 2021).
  103. Nunes, L.; Rodrigues, A.; Loureiro, L.; Sá, L.; Matias, J. Energy Recovery from Invasive Species: Creation of Value Chains to Promote Control and Eradication. Recycling 2021, 6, 21. [Google Scholar] [CrossRef]
  104. Holm, J.K.; Henriksen, U.B.; Hustad, J.E.; Sørensen, L.H. Toward an Understanding of Controlling Parameters in Softwood and Hardwood Pellets Production. Energy Fuels 2006, 20, 2686–2694. [Google Scholar] [CrossRef]
  105. Stelte, W.; Holm, J.K.; Sanadi, A.R.; Barsberg, S.; Ahrenfeldt, J.; Henriksen, U.B. A study of bonding and failure mechanisms in fuel pellets from different biomass resources. Biomass Bioenergy 2011, 35, 910–918. [Google Scholar] [CrossRef] [Green Version]
  106. TRADERSbiomass. 2021. Available online: https://tradersbiomass.com/online-biomass-exchangetrading-with-tradersbiomass/trading-with-tradersbiomass/ (accessed on 27 August 2021).
  107. Jenkins, B.M.; Baxter, L.L.; Miles, T.R., Jr.; Miles, T.R. Combustion properties of biomass. Fuel Process. Technol. 1998, 54, 17–46. [Google Scholar] [CrossRef]
  108. Verlinden, M.; Broeckx, L.; Zona, D.; Berhongaray, G.; De Groote, T.; Camino-Serrano, M.; Janssens, I.; Ceulemans, R.J. Net ecosystem production and carbon balance of an SRC poplar plantation during its first rotation. Biomass Bioenergy 2013, 56, 412–422. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (Left) Mean annual values of plant height (H), quadratic mean diameter (QMD) and woody biomass dry weight of the aboveground part (AGWB) for the eleven taxa assessed in the Mo1 trial. The bars indicate the standard error (SE) at the end of the fourth year of each rotation (R04S4 the 1st rotation, and R18S4 the 2nd one); (Right) Boxplots of quadratic mean diameter (QMD) and plant height (H) in the Mo2 (a), Ra4 (b) and Mo3 (c) trials for the last year of each rotation. First rotation (grey colors) and second rotation (blue color). Different letters indicate significant differences between taxa or treatments (p < 0.001) at the end of each rotation (black letters for the 1st rotation; italic blue letters for the 2nd rotation). AF2 (Populus × euroamericana “AF2”), Rp-H and Rp-G (Robinia pseudoacacia cultivar “Nyírségi” and German provenance, respectively); Up (Ulmus pumila). Treatments: 100% AF2 (Rp0-AF100); 25% Rp-G + 75% AF2 (Rp25-AF75); 50% Rp-G + 50% AF2 (Rp50-AF50); 75% Rp-G + 25% AF2 (Rp75-AF25); 100% Rp-G (Rp100-AF0).
Figure 1. (Left) Mean annual values of plant height (H), quadratic mean diameter (QMD) and woody biomass dry weight of the aboveground part (AGWB) for the eleven taxa assessed in the Mo1 trial. The bars indicate the standard error (SE) at the end of the fourth year of each rotation (R04S4 the 1st rotation, and R18S4 the 2nd one); (Right) Boxplots of quadratic mean diameter (QMD) and plant height (H) in the Mo2 (a), Ra4 (b) and Mo3 (c) trials for the last year of each rotation. First rotation (grey colors) and second rotation (blue color). Different letters indicate significant differences between taxa or treatments (p < 0.001) at the end of each rotation (black letters for the 1st rotation; italic blue letters for the 2nd rotation). AF2 (Populus × euroamericana “AF2”), Rp-H and Rp-G (Robinia pseudoacacia cultivar “Nyírségi” and German provenance, respectively); Up (Ulmus pumila). Treatments: 100% AF2 (Rp0-AF100); 25% Rp-G + 75% AF2 (Rp25-AF75); 50% Rp-G + 50% AF2 (Rp50-AF50); 75% Rp-G + 25% AF2 (Rp75-AF25); 100% Rp-G (Rp100-AF0).
Forests 12 01337 g001
Figure 2. Mean annual values of woody biomass dry weight of the aboveground part of the plant (AGWB) for the four mixture treatments assessed in the Mo2 trial. First rotation (grey colors) and second rotation (blue color). The bars indicate the standard error (SE) at the end of the fourth year of each rotation (R04S4 for the 1st rotation, and R18S4 for the 2nd). Different letters indicate significant differences between treatments (p < 0.001) at the end of each rotation (black letters for the 1st rotation; italic blue letters for the 2nd rotation). Treatments: 100% AF2 (Rp0-AF100); 25% Rp-G + 75% AF2 (Rp25-AF75); 50% Rp-G + 50% AF2 (Rp50-AF50); 75% Rp-G + 25% AF2 (Rp75-AF25); 100% Rp-G (Rp100-AF0).
Figure 2. Mean annual values of woody biomass dry weight of the aboveground part of the plant (AGWB) for the four mixture treatments assessed in the Mo2 trial. First rotation (grey colors) and second rotation (blue color). The bars indicate the standard error (SE) at the end of the fourth year of each rotation (R04S4 for the 1st rotation, and R18S4 for the 2nd). Different letters indicate significant differences between treatments (p < 0.001) at the end of each rotation (black letters for the 1st rotation; italic blue letters for the 2nd rotation). Treatments: 100% AF2 (Rp0-AF100); 25% Rp-G + 75% AF2 (Rp25-AF75); 50% Rp-G + 50% AF2 (Rp50-AF50); 75% Rp-G + 25% AF2 (Rp75-AF25); 100% Rp-G (Rp100-AF0).
Forests 12 01337 g002
Figure 3. Mean annual values of woody biomass dry weight of the aboveground part of the plant (AGWB) for the four taxa assessed in the Ra4 (left) and Mo3 (right) trials. First rotation (grey colors) and second rotation (blue color). The bars indicate the standard error (SE) at the end of the third (Ra4 trial) or the fifth (Mo3 trial) year of each rotation. R03S3 for the 1st rotation and R16S3 for the 2nd one in the Ra4 trial; R05S5 for the Mo3 trial. Different letters indicate significant differences between taxa (p < 0.001) at the end of each rotation (black letters for the 1st rotation; italic blue letters for the 2nd rotation). AF2 (Populus × euroamericana “AF2”), Rp-H and Rp-G (Robinia pseudoacacia cultivar “Nyírségi” and German provenance, respectively); Up (Ulmus pumila).
Figure 3. Mean annual values of woody biomass dry weight of the aboveground part of the plant (AGWB) for the four taxa assessed in the Ra4 (left) and Mo3 (right) trials. First rotation (grey colors) and second rotation (blue color). The bars indicate the standard error (SE) at the end of the third (Ra4 trial) or the fifth (Mo3 trial) year of each rotation. R03S3 for the 1st rotation and R16S3 for the 2nd one in the Ra4 trial; R05S5 for the Mo3 trial. Different letters indicate significant differences between taxa (p < 0.001) at the end of each rotation (black letters for the 1st rotation; italic blue letters for the 2nd rotation). AF2 (Populus × euroamericana “AF2”), Rp-H and Rp-G (Robinia pseudoacacia cultivar “Nyírségi” and German provenance, respectively); Up (Ulmus pumila).
Forests 12 01337 g003
Figure 4. Current annual increment (CAI) and mean annual increment (MAI) of the woody biomass dry weight of the aboveground part of the plant (AGWB) throughout the period studied for the 1st and 2nd rotations. Data from all field trials and four taxa are shown as a whole. Taxa: AF2 (Populus × euroamericana “AF2”), Rp-H and Rp-G (Robinia pseudoacacia cultivar “Nyírségi” and German provenance, respectively); Up (Ulmus pumila). Equations of the fit lines: (i) 1st rotation, y = −2026x2 + 13,185x − 7043 (R2 = 0.69) for CAI, and y = −0.894x2 + 6902x − 2487 (R2 = 0.85) for MAI; (ii) 2nd rotation, y = -2427x2 + 13,381x − 2151 (R2 = 0.65) for CAI, and y = −0.824x2 + 5413x + 4280 (R2 = 0.62) for MAI.
Figure 4. Current annual increment (CAI) and mean annual increment (MAI) of the woody biomass dry weight of the aboveground part of the plant (AGWB) throughout the period studied for the 1st and 2nd rotations. Data from all field trials and four taxa are shown as a whole. Taxa: AF2 (Populus × euroamericana “AF2”), Rp-H and Rp-G (Robinia pseudoacacia cultivar “Nyírségi” and German provenance, respectively); Up (Ulmus pumila). Equations of the fit lines: (i) 1st rotation, y = −2026x2 + 13,185x − 7043 (R2 = 0.69) for CAI, and y = −0.894x2 + 6902x − 2487 (R2 = 0.85) for MAI; (ii) 2nd rotation, y = -2427x2 + 13,381x − 2151 (R2 = 0.65) for CAI, and y = −0.824x2 + 5413x + 4280 (R2 = 0.62) for MAI.
Forests 12 01337 g004
Table 1. Climatic conditions at the field plots during the three time periods studied. Mo (Los Morales), Ra (La Rábida).
Table 1. Climatic conditions at the field plots during the three time periods studied. Mo (Los Morales), Ra (La Rábida).
VariableField Trials
Mo1 and Mo2
Field Trial
Mo3
Field Trial
Ra4
Time period (years)2011–20182015–20192015–2020
Absolute maximum temperature (°C)44.844.840.6
Warmest month’s average maximum temperature (°C)34.134.731.9
Mean annual temperature (°C)13.914.116.8
Coldest month’s average minimum temperature (°C)−0.7−0.84.4
Absolute minimum temperature (°C)−13.7−8.4−4.0
Annual precipitation (mm)319.8332.8483.6
Range of annual precipitation214.2–476.3214.2–461.9236.6–646.2
Summer precipitation from June to August (mm)31.022.17.6
Vegetative period (>12.5 °C mean monthly temperature) (month)6.26.28.3
Annual reference crop evapotranspiration, ETo (mm)1396.81399.91074.4
Table 2. Physical–chemical characteristics of the upper soil layer (0–30 cm) at the beginning of the experiments (mean (SE)). Mo (Los Morales), Ra (La Rábida). See also Table S1.
Table 2. Physical–chemical characteristics of the upper soil layer (0–30 cm) at the beginning of the experiments (mean (SE)). Mo (Los Morales), Ra (La Rábida). See also Table S1.
VariableField Trials
Mo1 and Mo2
Field Trial
Mo3
Field Trial
Ra4
Soil sampling depth (cm)0–300–300–1515–30
Exchangeable Sodium Percentage (%)3.9 (0.2)5.4 (0.1)9.1 (0.3)12.7 (0.1)
Organic Matter (%)2.93 (0.09)2.67 (0.03)1.22 (0.02)0.53 (0.05)
C/N ratio10.64 (0.46)9.16 (0.04)11.27 (0.04)13.74 (0.76)
Organic Carbon (%)1.51 (0.04)1.46 (0.03)0.71 (0.02)0.31 (0.03)
Active limestone (%)7.57 (0.31)3.31 (0.03)1.27 (0.05)1.71 (0.08)
N (%)0.16 (0.01)0.15 (0.01)0.06 (0.01)0.04 (0.01)
Total P (mg kg−1) (1)88.8 (2.5)78.4 (5.0)59.6 (2.0)42.9 (1.5)
Available K (meq/100 g)1.37 (0.09)1.40 (0.07)0.27 (0.01)0.21 (0.01)
(1) The available P (Olsen Method) averaged: 5.21 (0.41) mg kg−1 (mean (SE)) in the Mo1 and Mo2 trials; 4.02 (0.36) in the Mo3; and 3.59 (0.19) and 2.35 (0.20) for the 0–15 and 15–30 cm layers, respectively, in the Ra4 trial.
Table 3. General characteristics of the four field trials at the two locations, Los Morales (Mo) and La Rábida (Ra). See also Figures S1–S3.
Table 3. General characteristics of the four field trials at the two locations, Los Morales (Mo) and La Rábida (Ra). See also Figures S1–S3.
VariableField Trial
Mo1
Field Trial
Mo2
Field Trial
Mo3
Field Trial
Ra4
Study period (years)2011–20182011–20182015–20192015–2020
Number of taxa/treatments used11 taxa5 treatments4 taxa4 taxa
Taxa/treatments codes *Adige; AF2; Oudenberg; I214; Raspalje; Pn; Pf; Aa; Pxh; Rp-H; Up.Rp0-AF100;
Rp25-AF75;
Rp50-AF50;
Rp75-AF25;
Rp100-AF0.
Rp-H; Rp-G; AF2; Up.Rp-H; Rp-G; AF2; Up.
Number of years per rotation4453
Rotations studied (number of years)2 (4 + 4)2 (4 + 4)1 (5)2 (3 + 3)
Experimental unit (ExpU):
Length (m)9.69.6246
Number of lines3343
Number of plants per line1616306
Space between lines (m)3.03.03.32.0
Space between plants into a line (m)0.60.60.81.0
Number of plants in the ExpU484812018
Area of the ExpU (m2)86.486.4316.836
Plantation density (plants ha−1)5556555637885000
Number of plants measured per ExpU and year1010204
Taxon/treatment replicates (number of ExpU)4444
Total number of plants in the plot21129601920288
* Four hybrid clones of Populus × euroamericana (Dode) Guinier (clones “Adige”, “AF2”, “Oudenberg”, I214); one hybrid clone of Populus × interamericana van Broekhuizen (clone “Raspalje”); Pf, one clone of Paulownia fortunei (Seem.) Hemsl. (clone “UHU”); Rp-H (Robinia pseudoacacia, improved cultivar “Nyírségi”); Rp-G (Robinia pseudoacacia, commercial nursery plants from Germany); commercial nursery plants of Populus nigra L. (Pn), Ailanthus altissima (Mill) Swingle. (Aa), Platanus × hispanica Mill. ex Münchh. (Pxh), and Ulmus pumila L. (Up). Treatments (combinations of the five different species compositions): 100% AF2 (Rp0-AF100); 25% Rp-G + 75% AF2 (Rp25-AF75); 50% Rp-G + 50% AF2 (Rp50-AF50); 75% Rp-G + 25% AF2 (Rp75-AF25); 100% Rp-G (Rp100-AF0).
Table 4. Allometric equations relating the stem diameter at 10 cm above ground (D, mm) to the aboveground dry weight (AGWB, g). AGWB = a·Db. Stems with D ranging from 10 to 140 mm were used. Significance level, p < 0.001 for all cases. n = sample size.
Table 4. Allometric equations relating the stem diameter at 10 cm above ground (D, mm) to the aboveground dry weight (AGWB, g). AGWB = a·Db. Stems with D ranging from 10 to 140 mm were used. Significance level, p < 0.001 for all cases. n = sample size.
TaxonabR2n
Populus × euroamericana “Adige”0.0432.70740.98130
Populus × euroamericana “AF2”0.0342.70130.98638
Populus × euroamericana “Oudenberg”0.0762.63210.97830
Populus × euroamericana “I214”0.0512.6900.89130
Populus × interamericana “Raspalje”0.052.7470.98130
Populus nigra L.0.0452.7000.99630
Paulownia fortunei “UHU” (D < 32.3 mm)0.0053.0330.97015
Paulownia fortunei “UHU” (D > 32.3 mm)(1)(1)0.97217
Robinia pseudoacacia0.1132.4820.98032
Ulmus pumila0.0652.6120.97832
Ailanthus altissima0.0093.0420.98030
Platanus × hispanica0.0732.6300.98930
(1)Paulownia fortunei “UHU”: when D ≥ 32.3 mm, the following equation was used: TDW = (0.003·D3.074 + (0.435·D2 + 29.981 D − 11.651))/2.
Table 5. Technical characteristics of machinery, operation inventory and costs, and energy equivalent.
Table 5. Technical characteristics of machinery, operation inventory and costs, and energy equivalent.
Cultivation Work/MaterialTractorCostFuelEnergy (9)CO2 (9)Efficiency
(kW)(EUR ha−1)(L h−1)(MJ h−1)(kg h−1)(h ha−1)
50 cm deep ploughing (1)11011515.078.54.042.00
Disk harrowing (1)74458.761.53.101.40
Rotary harrowing (1)74498.761.53.101.50
Mechanical fertilizing (broadcaster) (2)52553.532.41.490.33
Chemical weeding (sprayer) (2)52543.532.41.490.33
Mechanical planting (1)11025012.078.54.042.50
Mechanical weeding (2)52303.532.41.491.25
Harvesting (cutting and field stocking) (2)74908.061.53.101.70
Harvesting (chipping and trailer loading) (1)(8)13238520.578.54.042.30
Grubbing and land clearing (1)(4)13247520.578.54.042.50
(EUR t−1)(L km−1)(MJ km−1)(kg km−1)(t payload−1)
Chip transportation (truck with trailer, 25 km) (3)2574.50.4251.150.0622.0
(EUR kg−1) (MJ kg−1)(kg kg−1)(kg ha−1)
Fertilizer N/P/K (5) 0.5 45/30/53.0/2.0/0.3500/333.3
Chelated fertilizer 6.2 1006.7
(EUR L−1) (MJ L−1)(kg L−1)(L ha−1 y−1)
Herbicide 10.5 1067.12.50
(EUR ha−1) (MJ kg−1)(kg kg−1)(kg ha−1 y−1)
Irrigation (material) (PE/PVC) (6) 1490 78.6/575.3/3.8523.13
(EUR h−1) (MJ m−3)(kg m−3)(m3 ha−1 y−1)
Irrigation (water pumping) (7) 0.21 0.540.0364000
(EUR plant−1)
Plants 0.35
(1) Service life = 10,000 h. (2) Service life = 8000 h. (3) Service life = 450,000 km. (4) Stump removal and land clearing have been taken into account for the recovery of the soil at the end of the cultivation cycle (15–17 years). (5) 500 kg ha−1 at the beginning was applied by a broadcaster, plus 333.3 kg ha−1 as fertigation from the 1st to the last rotation. (6) To determine energy costs and CO2 emissions from the manufacture of irrigation equipment, only plastic materials (PE, PVC) have been taken into account, while the cost has taken into account all the irrigation materials. (7) A 10.3 kW irrigation pump has been taken into account. (8) Since the plants will have a stem diameter greater than 10 cm at the end of the rotation, cut-and-chip harvesters operate with difficulty (diameter too thick for them), and forest harvesters are very expensive (diameter too thin for them). That is why a two-phase harvest has been chosen, first cutting with a saw disc and field stocking until the moisture content drops below 35%, and secondly chipping and loading in the trailer simultaneously. (9) The energy consumed by the machinery and the CO2 emissions refer to the manufacturing and maintenance process during its service life, considering only the time (hours) worked per hectare, and this energy must be added to that contained in the fuel consumed; the CO2 emissions include emissions of other greenhouse gases, therefore they are CO2-equivalent emissions; in the case of fertilizers, plastic materials, herbicide and operation of the irrigation pump, the transformation of energy consumed to CO2 equivalent emissions has used the following equivalence: 241 g of CO2 per kWh used, which is the one corresponding to the Spanish electricity mix during 2019.
Table 6. Significance level (p values) of the fixed effects of the two- and three-way ANOVA used to evaluate soil characteristics for each assay carried out. Measurements took place at the end of the 1st and the 2nd rotations for Mo1, Mo2 and Ra4 trials, while for Mo3 trial it was measured in the middle (3rd year) and at the end (5th year) of its 1st rotation.
Table 6. Significance level (p values) of the fixed effects of the two- and three-way ANOVA used to evaluate soil characteristics for each assay carried out. Measurements took place at the end of the 1st and the 2nd rotations for Mo1, Mo2 and Ra4 trials, while for Mo3 trial it was measured in the middle (3rd year) and at the end (5th year) of its 1st rotation.
VariableFieldMain EffectEffect (1)
TrialTX or TRYearSLYInteractions
ESP (%)Mo1<0.001<0.001 0.133
Mo20.027<0.001 0.215
Mo30.010<0.001 0.052
Ra4<0.001<0.001<0.001≥0.200
C/N ratioMo10.051<0.001 0.390
Mo20.0100.001 0.076
Mo30.4140.002 0.638
Ra40.016<0.0010.910≥0.076
Organic carbon (%)Mo10.026<0.001 0.376
Mo20.132<0.001 0.748
Mo30.519<0.001 0.781
Ra40.9070.010<0.001≥0.289
Active limestone (%)Mo10.220<0.001 0.750
Mo20.459<0.001 0.593
Mo30.736<0.001 0.933
Ra40.0410.0020.302≥0.202
N (%)Mo10.0110.001 0.301
Mo20.0050.004 0.100
Mo30.487<0.001 0.752
Ra40.0030.009<0.001≥0.054
Total P (mg kg−1)Mo10.941<0.001 0.873
Mo20.829<0.001 0.785
Mo30.626<0.001 0.962
Ra40.004<0.001<0.001≥0.051
Available K (meq/100 g)Mo10.051<0.001 0.055
Mo20.0680.045 0.068
Mo30.7890.040 0.959
Ra40.292<0.001<0.001≥0.063
ESP, exchangeable sodium percentage; TX/TR, taxa (TX, for Mo1, Mo3 and Ra4) or cultivation treatments (TR, for Mo2); Year (at the end of the 1st and the 2nd rotations); SLY, soil layer (0–15 and 15–30 cm); the durations of the first two rotations for Mo1, Mo2 and Ra4 trials were 4 + 4, 4 + 4, and 3 + 3 years, respectively; (1) effect interactions include pairwise interactions (all field trials) and triple interactions (only Ra4 trial).
Table 7. Average values, for all the taxa as a whole, of some physical–chemical characteristics of the upper soil layer of the soils (0–30 cm) at the end of the 1st and 2nd rotations for Mo1, Mo2 and Ra4 trials, while for the Mo3 trial (1) these were measured in the middle (3rd year) and at the end (5th year) of the 1st rotation. The asterisks mean significant difference from the initial soil. In such cases, (+) and (−) mean that the measured parameter significantly increased or decreased, respectively, during the cultivation period. ns indicates non-significant differences with respect to the initial soil.
Table 7. Average values, for all the taxa as a whole, of some physical–chemical characteristics of the upper soil layer of the soils (0–30 cm) at the end of the 1st and 2nd rotations for Mo1, Mo2 and Ra4 trials, while for the Mo3 trial (1) these were measured in the middle (3rd year) and at the end (5th year) of the 1st rotation. The asterisks mean significant difference from the initial soil. In such cases, (+) and (−) mean that the measured parameter significantly increased or decreased, respectively, during the cultivation period. ns indicates non-significant differences with respect to the initial soil.
VariableRotation (1)Field Trial
Mo1Mo2Mo3Ra4
ESP (%)1st5.1 * (+)5.3 * (+)6.8 * (+)14.1 * (+)
2nd6.4 * (+)7.0 * (+)6.7 * (+)13.6 * (+)
C/N ratio1st10.56 ns10.62 ns9.16 ns8.31 * (−)
2nd9.49 * (−)9.48 * (−)10.35 * (+)8.59 * (−)
Organic carbon (%)1st1.68 * (+)1.70 * (+)1.36 ns0.38 ns
2nd1.71 * (+)1.70 * (+)1.59 * (+)0.57 (+)
Active limestone (%)1st6.35 * (−)6.29 * (−)2.67 * (−)1.28 * (−)
2nd4.91 * (−)4.85 * (−)2.88 * (−)1.30 * (−)
N (%)1st0.16 ns0.16 ns0.15 ns0.05 ns
2nd0.17 * (+)0.17 * (+)0.16 * (+)0.06 * (+)
Total P (mg kg−1)1st83.8 * (−)83.4 * (−)79.4 * (−)35.0 * (−)
2nd74.4 * (−)70.7 * (−)79.4 * (−)32.0 * (−)
Available K (meq/100 g)1st1.42 ns1.57 * (+)1.54 ns0.31 * (+)
2nd1.63 * (+)1.27 ns1.64 * (+)0.28 * (+)
ESP, exchangeable sodium percentage. The durations of the first two rotations for Mo1, Mo2 and Ra4 trials were 4 + 4, 4 + 4, and 3 + 3 years, respectively.
Table 8. Physical–chemical and energetic characteristics (mean (SE)) of Populus × euroamericana “AF2”, Robinia pseudoacacia and Ulmus pumila, expressed on a dry weight basis, for aboveground woody biomass (AGWB), roots and litterfall, and all the trials as a whole. HHV, LHV: high and low heating value, respectively. Bdp, pellets bulk density. MDup, pellets mechanical durability. Moisturep, pellets moisture. WdB, wood density including bark. na, not analyzed.
Table 8. Physical–chemical and energetic characteristics (mean (SE)) of Populus × euroamericana “AF2”, Robinia pseudoacacia and Ulmus pumila, expressed on a dry weight basis, for aboveground woody biomass (AGWB), roots and litterfall, and all the trials as a whole. HHV, LHV: high and low heating value, respectively. Bdp, pellets bulk density. MDup, pellets mechanical durability. Moisturep, pellets moisture. WdB, wood density including bark. na, not analyzed.
VariablePopulus × euroamericana “AF2”Robinia pseudoacaciaUlmus pumila
AGWBRootsLitterfallAGWBRootsLitterfallAGWBRootsLitterfall
C (%)50.9
(1.1)
48.8
(1.1)
54.1
(0.9)
50.7
(0.8)
49.4
(0.7)
53.0
(1.3)
50.1
(0.8)
48.5
(0.9)
52.6
(2.0)
N (%)0.38
(0.04) a
0.79
(0.05)
1.10
(0.06)
0.54
(0.03) b
0.74
(0.17)
1.22
(0.06)
0.47
(0.07) ab
0.80
(0.13)
1.10
(0.05)
P (%)0.09
(0.01)
0.09
(0.03)
0.18
(0.03)
0.06
(0.06)
0.10
(0.07)
0.11
(0.06)
0.08
(0.04)
0.08
(0.06)
0.11
(0.05)
K (%)0.30
(0.06)
0.28
(0.04)
0.23
(0.05)
0.31
(0.03)
0.24
(0.04)
0.22
(0.03)
0.29
(0.03)
0.24
(0.04)
0.20
(0.04)
Ca (%)0.57
(0.10) b
1.26
(0.11)
2.99
(0.35)
0.38
(0.02) a
0.48
(0.05)
2.41
(0.08)
0.57
(0.11) b
0.39
(0.08)
2.87
(0.16)
Mg (%)0.13
(0.04)
0.17
(0.05)
0.37
(0.04)
0.14
(0.05)
0.17
(0.04)
0.22
(0.04)
0.11
(0.04)
0.22
(0.04)
0.33
(0.06)
S (%)0.06
(0.01) b
0.08
(0.02)
0.02
(0.01)
0.04
(0.01) a
0.08
(0.02)
0.02
(0.01)
0.04
(0.01) a
0.10
(0.02)
0.02
(0.01)
Fe (mg kg−1)78 (12) b78 (21)303 (23)58 (10) ab82 (12)301 (15)53 (9) a66 (11)272 (23)
Mn (mg kg−1)7 (1)19 (2)34 (2)6 (1)12 (2)25 (3)5 (2)13 (2)26 (3)
Zn (mg kg−1)40 (10)29 (8)77 (27)16 (5)28 (4)46 (20)25 (5)28 (4)45 (15)
Cl (mg kg−1)890 (40) b817 (77)3360 (220)744 (23) a812 (64)2966 (123)869 (40) b843 (65)3132(160)
Ash (%)2.05 (0.16) anana2.04 (0.10) anana2.67 (0.13) bnana
HHV (MJ kg−1)19.2 (0.2)nana19.3 (0.2)nana19.4 (0.2)nana
LHV (MJ kg−1)17.8 (0.3)nana18.0 (0.2)nana18.1 (0.3)nana
Bdp (kg m−3) (1)671 (13)nana651 (12)nana636 (14)nana
MDup (%) (1)97.7 (1.5)nana98.5 (1.7)nana97.4 (1.8)nana
Moisturep (%) (1)6.6 (0.5)nana6.8 (0.4)nana6.8 (0.5)nana
LHVp (MJ kg−1) (2)16.5 (0.2)nana16.6 (0.2)nana16.7 (0.2)nana
WdB (kg m−3)0.42 (0.05) anana0.77 (0.07) bnana0.68 (0.06) bnana
Significance level (p-values) for AGWB of the taxa effect: C, 0.605; N, 0.022; P, 0.857; K, 0.278; Ca, 0.001; Mg, 0.288; S, 0.018; Fe, < 0.001; Mn, 0.516; Zn, 0.055; Cl, <0.001; Ash, 0.040; WdB, <0.001. (1) Assessed according to ISO 17225-2:2014; the bark was included in the pelletization process. (2) Low heating value of pellets refers to a wet basis (the moisture content of pellets as received, i.e., the moisture shown above). Different letters within a row denote significant differences between taxa for AGWB.
Table 9. Crop assumptions deduced from each SRC field trial carried out: duration of crop rotations and cycles, cultivation practices, used material and application frequencies.
Table 9. Crop assumptions deduced from each SRC field trial carried out: duration of crop rotations and cycles, cultivation practices, used material and application frequencies.
Crop AssumptionMo1 and Mo2Mo3Ra4
Year per cycle161715
1st rotation (years)453
2nd and following rotations (years per rotation)443
Rotations per cycle445
Frequency(depending on the rotation: the 1st or subsequent rotations)
1st2nd to 4th1st2nd to 4th1st2nd to 5th
50 cm-deep ploughing101010
Disk harrowing101010
Rotary harrowing101010
Mechanical fertilizing (broadcaster)101010
Chemical weeding (sprayer)101010
Mechanical planting101010
Mechanical weeding101010
Harvesting (cutting with saw disc, field stocking)131314
Harvesting (chipping and trailer loading)131314
Grubbing and land clearing010101
Chip transportation (truck with trailer, 25 km)131314
Fertilizer N/P/K (mechanical fertilizing)101010
Fertilizer N/P/K (fertigation)1 × 43 × 41 × 53 × 41 × 34 × 3
Chelated fertilizer1 × 43 × 41 × 53 × 41 × 34 × 3
Herbicide101010
Irrigation (material) (PE/PVC)101010
Irrigation (water pumping)1 × 43 × 41 × 53 × 41 × 34 × 3
Plants101010
Table 10. Economic, energy and CO2 emission balances according to the cultivation assumptions set out above and in Table 5 and Table 9. AGWB, aboveground woody biomass expressed on a dry weight basis.
Table 10. Economic, energy and CO2 emission balances according to the cultivation assumptions set out above and in Table 5 and Table 9. AGWB, aboveground woody biomass expressed on a dry weight basis.
Crop AssumptionMo1 (Cycle = 16 Years)Mo3 (Cycle = 17 Years)Ra4 (Cycle = 15 Years)
Biomass yielded according to rotation1st2nd to 4th1st2nd to 4th1st2nd to 5th
AGWB (t ha−1) (1)42.40151.9251.35151.9231.80174.36
AGWB (t ha−1 year−1) (1)10.6012.6610.2712.6610.6014.53
Economic balanceIncomesCostsIncomesCostsIncomesCosts
Gross incomes and costs (EUR ha−1)15,93410,68716,66810,92316,90511,020
Difference (Inc.–Cos.) (EUR ha−1)524757455885
Ratio (Incomes/Costs)1.4911.5261.534
Gross incomes and costs (EUR ha−1 year−1)9966689806431127735
Difference (Inc.–Cos.) (EUR ha−1 year−1)328337392
Energy balanceGeneratedConsumedGeneratedConsumedGeneratedConsumed
Energy (GJ ha−1) (2)349815936591663711156
Difference (Gen.–Con.) (GJ ha−1)333934933555
Ratio (Generated/Consumed)22.0022.0423.79
Energy (GJ ha−1 year−1) (1)218.610.0215.29.8247.410.4
Difference (Gen.–Con.) (GJ ha−1 year−1)208.6205.4237.0
CO2 emission balanceFixedReleasedFixedReleasedFixedReleased
CO2 fixed and released (t ha−1)564.010.5596.111.0572.610.3
Difference (Fix.–Rel.) (t ha−1) (3)553.5585.1562.3
Difference (Fix.–Rel.) (t ha−1) (3)193.4208.4180.2
Ratio (Fixed/Released) (3)53.7154.1955.59
Ratio (Fixed/Released) (3)1.5221.5381.459
CO2 emission balance (t ha−1 year−1)35.250.6635.070.6538.170.69
Difference (Fix.–Rel.) (t ha−1 year−1) (3)34.5934.4237.48
Difference (Fix.–Rel.) (t ha−1 year−1) (3)12.0812.2612.01
Cost of energy (EUR GJ−1) (4)3.203.133.10
Inc. = incomes; Cos. = costs; Gen. = generated; Con. = consumed; Fix. = fixed; Rel. = released. (1) According to Figure 1, Figure 2, Figure 3 and Figure 4. (2) An average low heating value (LHV) of 18 MJ kg−1, on a dry weight basis, has been considered in the estimates. (3) Two estimates have been shown because part of the fixed CO2 could be released back into the atmosphere after the energy transformation (i.e., AGWB combustion); the top row shows the fixed CO2 equivalent in the aboveground and underground woody biomass, the litter and the soil OM (that is, the equivalent of CO2 to the C content), while the bottom row shows the fixed CO2 equivalent from which the CO2 released after the combustion of aboveground woody biomass has been subtracted; the CO2 equivalents contained in AGWB were 360.1, 376.7 and 382.1 t ha−1, or 22.51, 22.16 and 25.47 t ha−1 year−1, in the Mo1, Mo3 and Ra4 trials, respectively. (4) The energy cost has been estimated as follows, costs/(energy generated–energy consumed).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alesso, S.P.; Tapias, R.; Alaejos, J.; Fernández, M. Biomass Yield and Economic, Energy and Carbon Balances of Ulmus pumila L., Robinia pseudoacacia L. and Populus × euroamericana (Dode) Guinier Short-Rotation Coppices on Degraded Lands under Mediterranean Climate. Forests 2021, 12, 1337. https://0-doi-org.brum.beds.ac.uk/10.3390/f12101337

AMA Style

Alesso SP, Tapias R, Alaejos J, Fernández M. Biomass Yield and Economic, Energy and Carbon Balances of Ulmus pumila L., Robinia pseudoacacia L. and Populus × euroamericana (Dode) Guinier Short-Rotation Coppices on Degraded Lands under Mediterranean Climate. Forests. 2021; 12(10):1337. https://0-doi-org.brum.beds.ac.uk/10.3390/f12101337

Chicago/Turabian Style

Alesso, Silvia Patricia, Raúl Tapias, Joaquín Alaejos, and Manuel Fernández. 2021. "Biomass Yield and Economic, Energy and Carbon Balances of Ulmus pumila L., Robinia pseudoacacia L. and Populus × euroamericana (Dode) Guinier Short-Rotation Coppices on Degraded Lands under Mediterranean Climate" Forests 12, no. 10: 1337. https://0-doi-org.brum.beds.ac.uk/10.3390/f12101337

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop