Introduction

The expansion of renewable energy production offers an alternative to reduce net greenhouse gas emissions while contributing to socio-economic and environmental objectives1. The European Union’s Renewable Energy Directives—RED and RED II, henceforth collectively referred to as EU RED—set ambitious goals for the substantial growth in consumption of energy from renewable sources by its member countries and set a pathway to commitments under the 2015 Paris Agreement2,3,4. In 2018, renewable energy accounted for 18% of total final energy consumed in the EU-28, with its generation increasing by 64% over the previous 10 years. Renewable energy sources have central roles in the EU-28 electricity and heating-and-cooling sectors with shares in final consumption of 32% and 20%, respectively. Onshore wind (28%) was the leading source of renewable electricity, and solid biomass (83%) the dominant source of renewable heating-and-cooling in 20175.

Under the EU RED, individual countries have developed National Action Plans to meet EU-wide goals tailored to particular resource availability and unique energy markets6. Former and current EU members such as the United Kingdom, Belgium, and Denmark rely heavily on imported biomass materials to meet their national renewable energy goals7. Wood fibers compressed into pellets emerged as an important type of commercially traded biofuel partly due to multiple factors. These included: fiber malleability, high energy content, relative ease of combustion, and existing supply-chains that facilitate competitive transportation and procurement costs further enabled by pellet standardization and environmental certification of forests8,9,10. The most recent statistics show that total trade of wood pellets across the EU-28 exceeded 18.2 million tons in 2018; 10.3 million tons were imported from non-European nations with the US as the largest trade partner11,12 (Fig. 1).

Figure 1
figure 1

EU-28 wood pellet (commodity number: 440131) total trade and imports from top-three non-European partners11.

Implementation of EU RED National Action Plans triggered the establishment of new manufacturing facilities pelletizing wood fibers across the US coastal southeast13,14. It is important to note that EU RED National Action Plans do not target forests of any particular region, rather, the supply of US wood pellets developed organically partly due to cost competitive and comparative production advantages15. Prior to the adoption of EU RED in 2009, installed annual wood pellet manufacturing capacity along the coast of the southeastern US did not exceed 0.3 million tons but by 2017 it had expanded to 7.3 million tons (Fig. 2). Most recent estimates place total installed annual capacity in this region at 9.0 million tons16. Exponential growth in wood pellet manufacturing within the US coastal southeast has raised concerns over potential impacts on local forests. Insights into EU RED related policy effects on US forests have been inferred from a review of the extant literature17,18, case studies14, and land use projections19,20. Some argue that new demand for wood fibers is poised to cause major forestland losses or to degrade forests’ structure, composition and nutrient cycles21,22. Others posit that new energy markets could have a positive effect on forestland area by preventing deforestation and encouraging investments in multi-purpose tree plantations13,23.

Figure 2
figure 2

Location of wood pellet mills and corresponding procurement landscapes (concentric circles) in the eastern US distinguishing facilities of at least 100-thousand tons of annual installed capacity and region, 2005–2017. US coastal southeast includes the states: Alabama, Florida, Georgia, Louisiana, Mississippi, North Carolina, South Carolina, and Virginia. Maps generated in R using sf package93,94.

The emergence of a resource-based manufacturing sector in direct response to market and public policy-driven demand offers an opportunity to assess associated changes in local natural resource conditions. During the same period following adoption of the EU RED, although of a lower magnitude, the number and capacity of wood pellet mills in the rest of the eastern US that exclusively supplies the domestic market also expanded10,24 (Fig. 2). Growth was partly the result of increased price competitiveness of wood pellets over heating oil and domestic policies promoting renewable heating25. The rise of the wood pellet industry in different areas responding to market and policy signals meets the general conditions for a natural experiment26. We assess changes in selected forest conditions for the US coastal southeast as compared to the rest of the eastern US to discern overall industry- associated effects and compare changes between them to infer EU RED-associated effects. Differences may be attributable to this policy intervention but, as in any other natural experiment, causality cannot be conclusive.

We follow the development of the wood pellet industry in the eastern US along with contemporaneous attributes of timberland (forestland capable of producing in excess of 1.4 m3 per ha per year not legally withdrawn from timber production, with a minimum area classification of 0.41 ha) structure and carbon (C). We georeference wood pellet mills and infer timberland conditions within mill-specific concentric radii that defined fiber procurement landscapes. We include (Table 1) three fundamental attributes for structure (number of live trees, number of growing-stock trees, number of standing-dead trees) and estimates of three carbon stocks (live trees, standing-dead trees, and soils). These attributes were estimated every 3 years between 2005 and 2017 using plot-level forest inventory data. Variation in timberland attributes within wood pellet mill procurement areas is modeled as a function of mill descriptors (initial operation, years of operation, manufacturing capacity) in addition to regional differences (coastal southeast and rest of eastern US), identification of period post US recession and of intensive domestic bioenergy policy interventions (cross-sections after year 2009), population (total inhabitants), access to export markets (within close proximity to ports exporting wood products), severe weather (extreme drought), and indicators of competition for wood fibers (overlap with procurement landscapes of wood-using power plants, pulp mills, other wood pellet mills). Econometric spatial panel estimation controls for interaction effects between regional, manufacturing size, and extreme weather covariates; idiosyncratic errors within each procurement area; and the potential for spatial autocorrelation. Detailed information on spatially-explicit covariates and estimation is presented in “Methods” section.

Table 1 Selected fundamental descriptors of timberland structure and carbon stocks within wood pellet industry procurement areas.

Within the larger scope of industry-related impacts on forest resources, our ex post analysis of timberland structure and carbon stocks makes direct contributions to the extant literature by: (1) developing and implementing an analytical framework designed to assess changes in fundamental forest attributes of ecological, environmental, and commercial importance within industrial procurement landscapes; (2) identifying spatially-explicit renewable energy policy-induced effects on timberlands; (3) expanding the body of research assessing impacts of public policy interventions on forest conservation outcomes to actively managed forestlands; and (4) offering the first comprehensive and systematic assessment of EU RED-induced changes on US timberlands. Global expansion in national monitoring capacities27 offer the opportunity to comprehensively examine changes in fundamental forest attributes using a spatially-explicit analytical framework as we describe in this study. Our findings have direct implications for sustainable forest management, natural resource-based industries, and to the ongoing debate regarding the development of criteria and monitoring systems to deem biomass a sustainable and carbon neutral energy source28. Next, we describe our methods for the spatially-explicit analysis of timberland conditions within wood pellet mill procurement landscapes and industry-associated changes. We present and discuss our results first focusing on contemporaneous marginal changes in timberland conditions associated with wood pellet mill descriptors for annual operation, regional differences, and large manufacturing capacity. The latter helps discern significant associations linked to economies of scale needed to reach export markets13,14. We then discuss other influential factors related to anthropogenic, forest ecological region, and extreme weather effects.

Methods

Wood pellet industry and EU-induced effects on timberlands

The establishment and operation of wood pellet mills can impact timberland conditions as wood fibers are sourced directly from forests or indirectly as a by-product of other wood-using industries29. New fiber demand can support the commercial flow of roundwood, pulpwood and logging residues procured directly from timberlands, while utilization of other wood product manufacturing residues can cause indirect land effects when creating markets for previously merchantable and non-merchantable materials3. The ecological footprint30 associated with procurement of low-cost raw materials is geographically delimited by transportation costs and a facility’s manufacturing capacity in line with established industry location theory principles14,25,31. For instance, in the short-term, mills of larger manufacturing capacity source raw materials from a larger procurement area. After a review of the current literature, a procurement radius of 80 km (50 British imperial miles) was assigned to facilities with a nominal annual manufacturing capacity of at least 100 thousand tons and a procurement radius of 48 km (30 British imperial miles) to those of lower capacity based on prevalent transportation distances and regional road tortuosity32,33. Moreover, wood fibers used for commercial purposes cannot originate from forestland legally withdrawn from timber production. These two criteria, proximity and timberland status, were used to define procurement landscapes for each of the 123 locations in the eastern US that at any point during our study period had an operational wood pellet mill. A procurement landscape parallels the definition of ‘sourcing area’ in RED II’s Article 2 as a “geographically defined area from which the forest biomass feedstock is sourced, from which reliable and independent information is available and where conditions are sufficiently homogeneous…”34. Information on wood pellet mill location, operational status, and manufacturing capacity was obtained from the most current and comprehensive database of US wood pellet mills available at the time35. As of March 2017, there were 123 wood pellet manufacturers in the eastern US that were operational at least 1 year between 2005 and 2017 (Fig. 3).

Figure 3
figure 3

Forest ecological regions and centroids with corresponding procurement landscapes (concentric circles) of locations in the eastern US in which a wood pellet mill was operational between 2005 and 201735,81. Maps generated in R using sf package93,94.

Evidence gathered from regional manufacturing and trade information shows that the totality of wood pellets exported to the EU are confined to facilities and ports within the US southeast coastal region16,36. Regional delimitation of EU RED-induced effects lent itself to the identification of a non-treated policy scenario to discern related impacts. We identified other states in the eastern US, which experienced a major expansion in wood pellet manufacturing over 2005–2017 period although in a lesser magnitude18,37, as a contrasting non-treated policy scenario (Fig. 2). Recent assessments of policy interventions on forest conditions have often relied on the statistical-based matching of policy treated and non-treated areas38,39 but such a statistically-based approach to identify counterfactual observations was not necessary in our case due to clarity between US regions responding to EU RED.

Inferring timberland attributes within wood pellet mills’ procurement landscapes

We derived timberland conditions from plot-level information collected by the US Forest Service’s National Forest Inventory and Analysis (FIA) Program. Each plot covers an area of 674.48 m2 and is comprised of four equally-sized subplots each of 7.32 m radius. A plot is randomly located within each of the 24.28 km2 hexagons used to seamlessly sample the conterminous 48 US states. Full details on sampling design, data collection procedures, and attribute estimation are available in the FIA program documentation40,41. The FIA program applies remote and ground sampling to represent plot, condition, and tree measurements.

Timberland attributes were estimated for the 123 locations in the eastern US that at any point during our study period had an operational wood pellet mill. Timberland conditions assessed through selected attributes were estimated over 3-year windows between 2005 and 2017 from a total of 38,626 plots using the most recent available data from the FIA program42. Timberland attributes included number of live trees, number of growing-stock trees, number of standing-dead trees, above and belowground carbon in live trees (tons), above and belowground carbon in standing-dead trees (tons), and carbon in organic soil (tons). Chosen attributes reflected our interest in examining landscape level changes in structure (e.g. number of live trees), and addressing concerns over impacts on wildlife habitat (e.g. number of standing-dead trees) and major carbon pools (e.g., carbon above and belowground in standing live and dead trees, and soils). Area-based attributes were derived using data and SQL queries for area expansion factors from the FIA program40,42,43. Estimated average sampling errors for selected attributes were all less than 5%. Timberland attributes were accordingly estimated for the years 2005, 2008, 2011, 2014 and 2017 to construct our panel dataset. Supplementary Figure S1 shows box-plots for timberland attributes across all 123 observations by year and all 31 states in our study area. Supplementary Figure S2 shows average values of selected timberland attributes within procurement areas of wood pellet mills, and state-wide estimates, distinguishing between areas of the US coastal southeast and the rest of the eastern US. The study area encompasses a variety of forest conditions that allowed consideration of industrial biophysical factors influencing wood pellet mill procurement landscapes.

Ex post assessment of timberland attributes associated with EU RED and other explanatory factors

Variation in selected timberland attributes was modeled as a function of variables triggered by implementation of EU RED and other factors that can explain changes in conditions within procurement areas of wood pellet mills. Explanatory factors included wood pellet mill descriptors, geographic region, competing industries, forest ecological regions, post recession and intensive domestic bioenergy policy period, population, and extreme weather events. Descriptive statistics for all covariates are disclosed in Table 2. Wood pellet mill descriptors included time-variant information on whether a new pellet mill was established and operational in a given period, the number of years it had been operational and whether the procurement area centroid was within 121 km (75 miles) of a waterway port exporting wood pellets44. Time-invariant mill descriptors included a category denoting the procurement area for a wood pellet mill of large manufacturing capacity (at least 100 thousand tons per year) to capture a larger footprint over the landscape reflected by correspondingly longer fiber sourcing radii; forest ecological region of largest expansion within a mill’s procurement landscape; and whether the procurement area centroid was within 121 km of a waterway port trading forest products. The latter helped control for any effects associated with wood procurement by industries other than the ones included in our covariates that can also trade via waterways44.

Table 2 Descriptions and descriptive statistics for model explanatory variables (n = 615).

EU RED contemporaneous associations with selected timberland attributes were determined by statistically significant differences between wood pellet mills’ procurement areas in the coastal southeast over the rest of the eastern US, and between larger- and smaller-scale facilities. The rest of the US eastern region outside the coastal southeast lent itself for statistical ex post comparison due to the presence of mixed conifer and broadleaf dominated landscapes, similar wood pelletization processes, comparable growth rates in industrial installed capacity, being located within the US facing the same federal laws and regulations, and consistency in forest inventory data collection18,37. However, differences exist as it is impossible to find perfect alternative scenarios to discern public policy impacts on natural resources45. Differences include non-identical forest ecosystems, different levels of industrial growth in recent years and a degree of heterogeneity in state-level forest policies18. The EU RED has reportedly influenced the establishment of larger wood pellet mills to benefit from economies of scale and to meet demand volumes, with exporting mills located near commercial ports for access to export routes14. Accordingly, differences between these descriptors and smaller manufacturing capacity helped discern for net EU RED policy effects.

Procurement landscapes for industries competing for similar wood fibers were defined and their contemporaneous degree of overlap with those of wood pellet mills were estimated. Competing industries included power plants consuming wood to generate electricity, and pulp mills. Radii for wood-using power plants and pulp mills were set at 80 km and 121 km to define commercial procurement areas for biopower and pulp industries, respectively25,46. We identified a period following the 2006–2009 recession of the US wood products industry to capture larger macroeconomic trends potentially omitted in other explanatory variables47. This post 2009 period coincides with a time of intensive US bioenergy public policy interventions, hence, this variable is also capturing any role associated with the expansion of related federal- and state-level programs25. Other explanatory factors were human population as a well-known element associated with changes in forest conditions, and a 1-year lagged variable identifying procurement areas affected by extreme drought48,49. See Supplementary Notes 1 and 2, and Figure S3 for details on the estimation of values for population and procurement area intersection to gauge degree of industry competition.

The issue of confounding can challenge empirical assessments of public policy impacts. A potentially prominent confounding effect mentioned in the recent literature on policy interventions and forest conditions is that of anticipatory behavior38 which could involve tree harvesting in expectation of policy implementation50. However, anticipatory behavior in the form of preemptive timber harvesting is unlikely to be an issue in this assessment. Any type of harvesting associated with the supply of wood fibers for pellets only occurs concurrent with the beginning of manufacturing operations. While anticipatory timber harvesting would be a rare occurrence in the context of our study area, wood pellet mills can be deliberately located within areas that are already endowed with suitable standing timber (which explains the common use of the term ‘woodsheds’14). Hence, we attempted to control for initial conditions within wood pellet mill procurement landscapes in our econometric analysis by identifying whether in any of our time periods a pellet mill was operational. We then examine annualized contemporaneous changes by including the length of milling operations. Finally, the localized nature of wood fiber markets motivated the adoption of an econometric model that controls for spatially-explicit and idiosyncratic conditions within procurement areas.

Spatial econometric panel analysis

A model was specified to ex post examine changes in the \(Y\) th timberland attribute across \(N\) observations over \(t\) periods for landscapes where at any point between 2005 and 2017 a wood pellet mill was operational (Eq. 1). \({Y}_{N,t}\), the \(N\times 1\) vector of observations of the dependent variable in period year \(t\), was modeled as a function of an overall parameter intercept (α); \(N\times 2\) matrix of wood pellet mill time-variant variables \({X}_{N,t}\) inclusive of mills being in operation in year t, number of years of operation at year t; \(N\times 10\) matrix \({F}_{N}\) of time-invariant descriptors denoting access to forest products ports, mill annual manufacturing capacity (at least 100 thousand tons), coastal southeast areas, and forest ecological regions; \(N\times 3\) matrix \({C}_{N,t}\) of time-varying degrees of overlap in procurement areas of wood pellets mills, electricity-generating plants using wood as an energy feedstock, and pulp mills; vectors of 2011–2017 post recession period \({R}_{N,t}\), population \({P}_{N,t}\), 1-year lagged variable for extreme drought \({D}_{N,t-1}\) for \(N\) observations over period year \(t\); and finally \(N\times 11\) matrix of interaction terms \({INT}_{N,t}\) between wood pellet mills presence, coastal southeastern location and years of operation, large manufacturing capacity, and extreme weather, across observations and years. Extreme drought was included as a 1-year lagged variable to correctly infer associations due to recording of drought in the summer months and surveying of forest plots throughout the year. Moreover, this timeframe allowed for the expression of any response between extreme drought and changes in timberland attributes to occur. We also included the \(N\times 1\) vector of location-specific disturbance term \({u}_{N,t}\) including \({W}_{N}\), an \(N\times N\) time-invariant weight matrix of row-normalized inverse distances between procurement landscape centroids, \(\rho \) a scalar autoregressive parameter, and \({\varepsilon }_{N,t}\) a \(N\times 1\) vector of period \(t\) random effects25,51,52,53. It can be shown that the disturbance vector \({u}_{N,t}\) can be uniquely solved in terms of the random effects vector \({\varepsilon }_{N,t}\)53,54. As such, a spatial model was estimated as:

$${Y}_{N, t}=\alpha +{X}_{N,t}{\beta }_{1}+{F}_{N}{\beta }_{2}+{C}_{N,t}{\beta }_{3}+{R}_{N, t>2009}{\beta }_{4}+{P}_{N,t}{\beta }_{5}+{D}_{N,t-1}{\beta }_{6}+ {INT}_{N, t}{\beta }_{7}+{u}_{N,t},$$
(1)

and

$${u}_{N,t}={\left({I}_{N}-\rho {W}_{N}\right)}^{-1}{\varepsilon }_{N,t} ,$$
(2)

where \(N=123\), \(t periods=\{2005, 2008, 2011, 2014, 2017\}\), and βs correspond to \(K\times 1\) vectors of the respective explanatory variables. \({W}_{N}\) (Eq. 2) is a row-normalized weight matrix of the inverse distances between pellet mills \({\omega }_{N}\) truncated at a distance of 161 km (100 British imperial miles)8 as the limit for spatial autocorrelation such that:

$${\omega }_{N}=\left\{\begin{array}{l}\frac{1}{{d}_{ij}},\quad {d}_{ij}<161km\\ 0,\quad otherwise\end{array}\right., i,j=\mathrm{1,2},\dots ,123 \; and\; i\ne j$$
(3)

The random effects error vector \({\varepsilon }_{N,t} \sim IID(0,{\sigma }_{\varepsilon }^{2})\) including unobserved idiosyncratic procurement area effects, and cross-sectional fixed and time-variant components, were tested for correlation with explanatory variables for correct estimation as a fixed or random model. Hausman test55 results (see Supplementary Note 3 and Table S1) favored a random effects specification, hence, the spatial panel regression with random effects was estimated using maximum likelihood52. β coefficients were standardized to estimate relative association of the ith covariate on dependent variable y to ease comparison across models as follows51:

$$Standardized \; {\beta }_{i}= {\beta }_{i} \frac{{S}_{{x}_{i}}}{{S}_{y}}$$
(4)

where the Sx and Sy correspond to the sample standard deviations for the covariate and dependent variable, respectively. Coefficients were also adjusted to calculate average hectare-level effects within procurement areas to afford direct interpretation and comparability across models (see Supplementary Note 4 for details on conversion followed to derive associations from coefficients to area-adjusted effects).

Interpretation of coefficients and wood pellet industry-related associations

The linear nature of our panel regression allows for the direct interpretation of explanatory variables’ parameterized associations with particular timberland attributes, else constant. Here we offer a few examples to clarify their interpretation in relation to inference of wood pellet industry associations.

  • Overall timberland conditions: The intercept parameter α controls for overall baseline timberland conditions across all wood pellet mill procurement landscapes irrespective of regional location or if pelletization operations have started. The parameter for ‘Large manufacturing capacity’ captures an overall shift in conditions associated with larger procurement areas. As expected, α intercept and larger procurement area β parameters are strongly statistically significant from zero [p value < 0.001]. The parameter for the ‘Coastal southeast region’ captures overall differences between procurement areas in the coastal southeast and the rest of the eastern US. Coastal southeast regional differences specifically associated with wood pellet mill presence and operation, years of operation, large manufacturing capacity, export market access, extreme drought are captured in corresponding interaction terms.

  • Initial wood pellet industry associated timberland conditions: Any differences in conditions at the initiation of wood pellet mill operations is captured in the parameter for ‘Mill present and operational’. Initial differences in timberland conditions within procurement areas of the coastal southeast, as compared with those in the rest of the US east, are captured in the parameter for variable ‘Pellet mill present and operational × coastal southeast region’.

  • Wood pellet industry-related temporal changes in timberland conditions: The parameters and respective interactions for the variable ‘Years of operation’ capture differences in timberland conditions associated with an additional year of operation. Any significant deviations from this trend within procurement landscapes of wood pellet mills in the coastal southeast and of large capacity mills are captured in the coefficients for ‘Years of operation × coastal southeast region’ and ‘Years of operation × large manufacturing capacity’, respectively.

Results

Estimated regression coefficients, area-adjusted marginal effects, and p values from the random-effects spatial panel model are presented in Tables 3 and 4. Standardized coefficients, to denote the relative importance of selected covariates on changes in fundamental timberland attributes, are included in Fig. 4. Coefficients that tested the association between wood pellet manufacturing and descriptors linked to EU RED, US coastal southeast and large-scale manufacturing, exhibited several statistically-significant associations (Fig. 4a). Compared with procurement areas of wood pellet mills in the rest of the eastern US, the operation of wood pellet plants in the US coastal southeast over time was associated with fewer numbers of live (p = 0.04) and growing-stock trees (p = 0.02) and lower carbon pools in soils (p = 0.03) to the rate of 8.21 ± 3.97 million trees and 1.55 ± 0.69 million trees, and 234.95 ± 105.97 thousand tons of carbon per procurement area per year (Table 3). These are the equivalent of 11.22 and 2.12 trees/ha/year, and 0.32 tons C/ha/year, respectively. Among large-scale facilities, their annual operation was associated with a decreasing trend in the number of standing-dead trees (p = 0.00) and more C stocks in live trees (p = 0.01). Average associations were found at the rate of 0.23 fewer standing-dead trees/ha/year, and 0.20 more tons C/ha/year in live-tree pools.

Table 3 Panel regression coefficients of wood pellet mill descriptors and interaction terms, along area-adjusted marginal effects (n = 615).
Figure 4
figure 4

Standardized coefficients and p values for (a) variables capturing annual effects of wood pellet mill operation interacted with US coastal southeast region and large-scale manufacturing categories, (b) fiber competition from power, pulp operations and other wood product factors, and (c) population, forest product exporting access from US coastal southeast, and southern mixed forest ecological region. Standardized coefficients represent the expected change (standard deviation) of selected timberland descriptors and carbon stocks associated with one standard deviation change in a given covariate.

We found discernible trends on timberland conditions within procurement landscapes that overlapped with other pellet mills as well as with other industries competing for wood fibers (Fig. 4b). A one percentage-point increase in overlap of wood pellet mill procurement areas, denoting greater competition, was associated with larger C stocks in above and belowground live (2.07 tons C/ha; p = 0.00) and standing-dead (0.06 tons C/ha; p = 0.00) trees pools, and in soils (1.32 tons C/ha; p = 0.02). A direct association was also found between greater overlap in procurement areas and a larger number of growing-stock trees. All selected timberland conditions showed statistically significant (p < 0.05) positive associations with the greater overlap of procurement areas between wood pellet mills and power plants using biomass for electricity generation. The overlap between wood pellet mill procurement landscapes with those of pulp mills was associated with an increase in carbon within live tree pools (0.68 tons/ha; p = 0.00) possibly as a result of land management that incentivizes growing wood fibers within procurement distances (Table 4).

Table 4 Regression coefficients of competition, domestic energy policy, population, extreme weather, and forest ecological regions, along area-adjusted marginal effects (n = 615).

Across all models, coefficients capturing marginal effects of population and close proximity to a forest product exporting port in the US coastal southeast had some of the largest absolute and standardized values (Fig. 4c). Population levels were inversely associated with all timberland descriptors used in our analysis. Procurement areas that had 1000 more inhabitants than the average had 244.61 ± 82.01 thousand fewer live trees, 45.09 ± 11.42 thousand fewer growing-stock trees, 2.88 ± 1.41 thousand fewer standing-dead trees, as well as 5.07 ± 2.18 thousand fewer tons of C in live trees, 297 ± 122 fewer tons of C in standing-dead trees and 11.22 ± 2.76 thousand fewer tons of C in soils in a given year. Wood pellet mill procurement landscapes within close proximity (i.e. 121 km) to a forest product exporting port in the US coastal southeast had 405.49 ± 187.09 million fewer live trees than those in the rest of the US east at a rate of − 553.72 trees/ha (p = 0.03). Comparatively, the standardized average marginal association of population with lower carbon pools in soils was 53-fold greater than that of annual operation of a large-scale wood pellet mill (Fig. 4a,b). This highlights the greater extent to which population can influence changes in certain timberland attributes than does the wood products industry. We also mention the importance of controlling for overall ecological conditions where wood pellet mils are located. For instance, the Oak-Hickory forest ecoregion showed statistically significant mean differences in nearly all timberland attributes (Fig. 4c).

Coefficients capturing extreme drought effects on selected timberland conditions showed some statistically significant associations. Extreme drought (Fig. 5) was associated with more live and standing-dead trees, as well as an increase in C stocks in standing dead-trees within procurement areas of the US coastal southeast. The ecological relationships leading to this observation may not be intuitive. A larger number of live trees might be explained by resprouting, and more standing-dead trees and C found in dead and soils pools as a consequence to higher mortality rates56. It is worth mentioning that differences in localized drought effects might also be influenced by recent precipitation events, among other contemporary and legacy factors, such as species composition and soil characteristics (e.g. parental material, depth).

Figure 5
figure 5

Location of wood-using power plants, pulp mills, and areas of reported extreme weather conditions (lagged 1 year), 2005–2017. Maps generated in R using sf package93,94.

Discussion

EU RED and US coastal southeast timberland conditions

The EU recognizes environmental risks associated with greater demand for biomass for energy and has issued non-binding recommendations on biomass sustainability criteria to meet EU RED targets57. The European Commission requires monitoring so that biomass used for energy, regardless of origin, does not cause deforestation or degradation of habitats3,58. In order to ensure compliance with EU-wide sustainability recommendations and any EU member-level sustainability criteria, the Sustainable Biomass Program (SBP) emerged as a private certification system designed to evaluate the legal and sustainable sourcing of wood pellets and chips used in large-scale energy generation59. It relies heavily on existing sustainable forest management and chain-of-custody certificate programs such as the Forest Stewardship Council and the Programme for Endorsement of Forest Certification. In addition to forest management and product chain-of-custody compliance, SBP certification offers an assessment of risk associated with feedstock in regards to legality and sustainability compliance and methods to collect and communicate data along the supply chain inclusive of greenhouse gas budgets.

Results identify diverse trends within procurement landscapes of pellet mills of the US coastal southeast from which wood pellets are exported to the EU. Fewer numbers of live and growing-stock trees with no decline in their corresponding C pools likely reflect on land management practices that favor growth in merchantable trees denoting more intensive silvicultural treatments. This trend can be explained by how new demand for wood pellets creates market incentives to grow and regrow wood fibers13,19. Moreover, this trend for greater demand of fibers sourced directly from timberland is reflected on the composition of wood fibers being pelletized. National industry surveys show fiber types used in wood pellet production were largely comprised of sawmill residues (69%, by weight)10 prior to EU RED implementation. Sawmill residues accounted for about 18% of fibers used in wood pellet manufacturing by 201716. Most fibers (49%) came from residual biomass of little to no commercial value (e.g., bark, logging residues, wood chips, post-consumer wood, unmerchantable wood), and roundwood/pulpwood (20%). The remaining fibers (13%) were sourced from residues of wood product manufacturing.

A fewer number of standing-dead trees over time within procurement areas of large-scale wood pellet mills might also be an indication of more intensive management practices that remove previously non-commercial biomass. This trend had been previously reported within the procurement area for an exporting wood pellet mill in the state of Georgia14. This result is of direct relevance to policy makers as the expansion of international markets might make large-scale wood pellet manufacturing facilities more prevalent. Some have argued that management that removes dead biomass prone to ignition could reduce the risk and intensity of wildfires, hence, it might point to potentially desirable outcomes60. However, it is important to note that woody biomass extraction continuously exceeding natural disturbances has been linked with the disruption of nutrient cycles to the detriment of site productivity and saproxylic communities, and could also negatively impact ecosystem structure as it relates to wildlife habitat and forest biodiversity61,62,63. There were also negative discernible effects on carbon stocks in soil pools associated with the annual continued operation of wood pellet facilities in the coastal southeast (Table 3). Our results point to a flux in C stocks across different pools in procurement areas of large-scale wood pellet mills of the US coastal southeast. The former descriptor was associated with higher C levels in live trees (407.85 ± 166.05 thousand tons of C), while the latter exhibited less C in soil pools (234.95 ± 105.97 thousand tons of C) per procurement area. On the balance, there has been a net contemporaneous positive effect.

It is worth noting that our ex post analysis of systematic changes in timberland conditions is confined to the boundaries for current commercial procurement distances. Issues of leakage beyond procurement areas are not addressed in this research. For instance, it is plausible that growing wood pellet demand had spillover effects to lands well outside mills’ procurement radii by shifting competition for similar fibers to new supply areas. Further, our results are limited to the time period evaluated (2005–2017) with EU RED targets instituted in 2009 and sustainability recommendations first adopted in 2015. Within this timeline of public policy interventions and commercial timberland management rotations, and given the most recently available forest inventory data, our analysis of 12 years is still limited in its ability to identify longer-term shifts in timberland conditions and C pools. In addition to the core attributes presented in this manuscript to assess changes within procurement landscapes, we recommend to also study systematic changes in biodiversity indicators.

On 13 March 2019 the European Commission proposed to the European Parliament the formal assessment and continuous evaluation of C stocks for areas from which biomass is sourced from to qualify for EU RED targets64. Our results show an incremental trend in above and belowground C in live trees within procurement areas of wood pellet mills of large manufacturing capacity, and an inverse trend of C in soils in those of the US coastal southeast. Moreover, our empirical estimates can be used in greenhouse gas life-cycle assessments that to-date have relied on assumed forest management regimes within particular market scenarios to gauge the C neutrality of Trans-Atlantic wood pellet trading65,66.

Competing wood fiber sectors, domestic policy, population and extreme weather

We found that an increase in the overlap between procurement areas of wood pellet mills and wood-using power plants was associated with significantly greater carbon stocks across nearly all selected C pools. As noted previously, this can be indicative of more intensive land management practices. This finding highlights that assessment of any potential associations between EU RED with timberland conditions should take into consideration other explanatory factors such as local consumption of wood fibers and domestic policies that create incentives to procure similar fibers. For instance, over the 2009–2015 period the use of wood for electricity generation in the US grew by 35% reaching 257.3 PJ in 201567 with 83% generated by power plants in the eastern US.

The coefficient denoting a period of intense US bioenergy policy interventions following the Great Recession (Table 4) showed that there was a statistically significant increase in carbon stocks in live and standing-dead tree C pools. This signal might point to latent effects of other federal and state level interventions over this period that have encouraged the growth of biomass for energy, such as the US Department of Agriculture’s Biomass Crop Assistance Program68, and the widespread adoption of state-level best management practices for the removal of woody biomass beyond traditional timber harvesting largely implemented since 200969. Prospects for the expansion of US wood pellet manufacturing are positive and some suggest that if an upward export trajectory to the EU were to continue US pellets may contribute up to 60% of annual European import demand70. Nonetheless, the fruition of such expectations might be contingent on air quality regulations. PM2.5 and PM10 emissions across Europe between 2005 and 2017 have reportedly increased by 11% and 7% because of growth in biomass burning, respectively, which has resulted in calls for policy action to ameliorate potential health impacts71.

The individual and compounded effects of EU and US renewable energy policies deserve continued examination. The US Department of Energy’s most recent outlook suggest a 21% increase in domestic electricity generation using wood by 203072. The recent announcement by the US Environmental Protection Agency that in future federal regulatory actions biomass from managed forests will be treated as carbon neutral when used for energy production at stationary sources73 might support even greater growth in corresponding domestic wood fiber demand. Federal energy policies will continue to overlap with state-level goals such as renewable energy portfolio standards that classify wood as a renewable feedstock74. Continued spatially-explicit monitoring of changing timberland conditions seems to be warranted, particularly as localized demand intensifies.

Our analysis shows that population had strongly negative standardized associations with all descriptors of timberland structure and C stocks. For instance, procurement areas that hosted more inhabitants had fewer numbers of live trees and growing-stock trees likely reflecting on competing land uses. Past assessments and projections for the US eastern region consistently point to population being one of the greatest drivers of changes in timberland conditions75,76. Extreme drought effects showed systematic associations with several timberland conditions in the US coastal southeast. Extreme drought is likely to increase in frequency and intensity under a changing climate. Its sole and compounded effects with other anthropogenic processes should continue to be monitored explicitly.

Conclusions

The EU, among other governments around the world77, are increasingly adopting policy strategies that seek the sustainable use of natural resources to support their economic development. The EU RED is a core element of the EU Bioeconomy Strategy78 that seeks, among other objectives, the modernization of primary energy production systems using biological resources. Environmental and sustainability issues have long been at the heart of EU policy79,80, yet the most recent emphasis on biological resources places even greater pressure on their utilization which could possibly counter the same environmental goals it seeks to attain18. The role of EU RED in the emergence of a new wood pellet industry in the US coastal southeast region, and the expansion of the industry in other areas of the eastern US in response to domestic demand and policy signals, offered the opportunity to assess related impacts on local forest resources.

Analysis of timberland conditions within wood pellet industry procurement landscapes of the eastern US showed upward and downward associations. On one hand, operation of wood pellet manufacturing in the US coastal southeast showed no significant changes in C pools of live nor standing-dead trees, and a concentration of live C in fewer growing-stock trees. Procurement areas of large-scale manufacturing facilities had higher levels of C in live-trees. Moreover, a greater overlap of wood pellet mills’ procurement areas exhibited discernible increases of selected C pools. These trends suggest that wood pellet production in the US has promoted greater growth in trees and an expansion in carbon pools in live trees—which is compliant with current EU RED biofuel trade requirements for the preservation of C stocks in biomass sourcing areas3. On the other hand, procurement areas of continuously operating wood pellet mills of the US coastal southeast and of large-scale capacity showed, respectively, lower carbon stocks in soils and a downward trend in the number of standing-dead trees. Understanding potential impacts of these trends on wildlife habitats and biodiversity requires additional analyses. Despite the positive association of wood pellet manufacturing with carbon in live trees, its long-term sustainability could be challenged if downward trends in basic forest structure indicators and other carbon pools were to continue.

Our findings suggest largely positive trends in timberland conditions; however, the presence of some potentially negative trends suggests that continued monitoring of localized impacts of wood pellet mill operations is important. We emphasize that monitoring must be observant of other localized conditions that can have concurrent and even greater effects than wood pellet manufacturing. Specifically, monitoring of timberland conditions within industrial procurement areas must also take into consideration population changes, expansion in wood fiber demand from other competing sectors, and extreme weather, while also controlling for overall forest ecological region conditions. Future research should examine how changes in timberland attributes affect forest biodiversity, wildlife habitat, and human well-being.