{'fraction_non_alphanumeric': 0.023242200950347774, 'fraction_numerical': 0.008952551477170993, 'mean_word_length': 4.874393203883495, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 1, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Impacts of the FLEGT Action Plan and the EU Timber Regulation on EU Trade in Timber Product

: The EU Timber Regulation (EUTR) is a key element in the efforts of the European Union to curb the trade in illegal timber products. This study helps to remedy the lack of systematic, statistical analysis of the EUTR’s potential impacts on international trade in timber products. Using cointegration intervention—or shock—models we quantify potential shifts in import prices and quantities of tropical hardwood lumber and oak lumber after the entry into effect of the EUTR. We further estimate import demand models to assess the relation between temperate and tropical hardwood products and whether there was a structural change in demand elasticities after the entry into force of the EUTR. The shock model analysis indicates, for most of the bilateral trade ﬂows where we observe cointegration and a signiﬁcant shock variable, increasing import prices and decreasing import quantities of tropical hardwood lumber following the EUTR start date, consistent with a contraction of the supply of tropical timber. The results of the import demand models do not give a clear indication as to whether oak lumber is a complementary or substitute product for tropical hardwood lumber, and there are no clear signs of structural changes in demand elasticities. Aside from the analysis, an important contribution of the paper is the procedure for building a long and homogeneous time series of tropical hardwood lumber.


Introduction
Illegal logging and the associated trade in timber products undermine legal timber trade and have detrimental environmental, economic, and social impacts [1]. The European Commission presented the Forest Law Enforcement, Governance, and Trade (FLEGT) Action Plan [2] as part of its endeavors to tackle illegal logging. Acknowledging the shared responsibility of exporters and importers, the objective is to eliminate illegal timber in international trade. A cornerstone of the FLEGT action plan is establishing Voluntary Partnership Agreements (VPAs) with timber producing and exporting countries, and the introduction of FLEGT-licensed timber [1]. A second key element of the plan is the EU Timber Regulation [3], which came into force on 3 March 2013. VPAs and the EUTR are meant to reinforce each other, addressing the supply (export) and demand (import) side of the timber product trade respectively (Figure 1).
A pertinent question then is what influences FLEGT and the EUTR, as well as other similar initiatives such as the US Lacey Act Amendment (LAA), have had on international trade in timber products. Estimates of the scale of illegal timber trade are surrounded by considerable uncertainty, since illegal trade is not recorded in trade databases [4]. Lawson and McFaul [5]-in assessing the effects of the LAA, FLEGT, and other efforts to reduce illegal logging through a combination of (i) surveys involving government officials, NGOs, and firms with (ii) the evaluation of trade and production data-found that illegal logging and exports of illegally sourced timber had decreased while timber product prices had risen since the beginning of the 2000s. Prestemon [6] performed econometric analysis on U.S. import data to assess the impact of the LAA. The results indicate general increases in prices and decreases in quantities of tropical hardwood imports to the U.S. from countries suspected of a high degree of illegal timber trade. Figure 1. The FLEGT Action Plan with its two main instruments and their main features. Source: [7].
Due to the rather short time the EUTR has been in force, there are yet no studies that have performed systematic, comprehensive statistical analysis regarding the impacts on international trade in timber products. Anecdotal evidence suggests that the main reason for declining tropical timber demand in the EU is substitution by temperate timber [8,9]. The European flooring industry in particular has reported a significant fall in the use of tropical timbers, to a large extent attributed to environmental concerns [10]-illegal timber trade is primarily associated with tropical hardwood [11]. It is claimed that EUTR may have been key in reinforcing this trend, through increased costs of conforming to the requirements of the EUTR [7][8][9]. Figure 2 gives some support to this claim, illustrating shrinking imports of tropical hardwoods from extra EU countries while temperate hardwood imports are rather stable through time. In addition, Figure 2 reports on intra EU trade to highlight three facts and a limitation of the Comext database. First, the intra EU trade of other sawnwood is an order of magnitude higher than that of oak and tropical hardwood lumber (note the differences in scales). Second, the intra EU trade of oak lumber has been increasing dramatically in recent years. Third there is a stable intra EU trade of tropical hardwood lumber throughout the period, and the limitation of our data source is that there is only data on the European partner country for such trade flows, but not the original tropical country of production. The extra-EU imports that are then re-exported within the EU are not accounted for in our analysis. This is why aggregating the EU as a whole leads to loss of information. It is not an issue for the demand model since in that case we are not investigating the source of tropical timber but the overall consumption of tropical timber.
The objective of this study is to assess trade effects of FLEGT and the EUTR: Have the patterns of EU imports of hardwood timber products changed as a result of FLEGT and the EUTR? We address this overriding question by analyzing the effects of FLEGT and the EUTR on imported quantities and prices of selected hardwood products. To accomplish this, we use two classes of statistical models. Hence, we use cointegration intervention models to quantify potential shifts-corresponding with the implementation of the EUTR-in the cointegrating relations of import prices and quantities of tropical hardwood lumber and oak lumber respectively. In addition, we estimate an import demand model to assess the relation between temperate and tropical hardwood products, trying to establish whether they are substitutes or complements and analyzing whether there was a structural change in demand elasticities after the entry into force of the EUTR. We hypothesize that the implementation of the EUTR resulted in a reduction of supply of tropical timber to EU member states, a result of the discontinuation of the supply of illegal timber and/or the disappearance of some of the legal supply due to the costs of conforming to the requirements of the EUTR. This backward shift of the supply curve would result in a new equilibrium market solution at the higher price, P2, and lower quantity, Q2, in the tropical timber market, ceteris paribus ( Figure 3). However, substitution by temperate timbers, such as oak, for tropical timber would in addition to increasing the equilibrium price and quantity in the temperate timber market through increased demand (Figure 3), result in a contraction of the demand (a backwards shift of the demand curve) for the tropical timber being substituted. While this would reinforce the decrease in the equilibrium quantity, it would counter the increase in the equilibrium price in the tropical timber market. The sign of the price movement will depend on the relative strength of the supply and demand shock respectively. The next section describes in detail the data and methods used in the analysis. Then results are reported, followed by the discussion and conclusions.

Shock Model
In an attempt to isolate the effect of the policy change from other economic variations, we analyze the price series of substitute or complementary products. Such products should be affected by economic variations in a similar way throughout the period, thus we look for a difference in how they relate to each other before and after the entry into force of the EUTR. Rather than estimating the effect of the shock directly with a dummy variable in a linear model, we follow the more robust approach of a cointegration intervention model [6]. The analysis starts by splitting the price series before and after the entry into force of the EUTR. We first estimate the relationship between the price of tropical products P rt and the price of temperate products P mt before the shock date (Equation (1)): ln P rt = γ 0 + γ 1 ln P mt + θ t . (1) The estimated coefficients are then used to generate pseudo errorsμ t after the shock date (Equation (2)). These pseudo errors contain information on how the relationship between the two products change after the shock. Third, theμ t vector is regressed onto a lagged version of itself to estimate λ, the size of the shock in Equation (3): The λ coefficient gives an indication as to how the shock impacted the relationship between tropical and temperate prices. Since we estimate the model with the log of prices, e λ − 1 can be interpreted as a percentage change of prices between a model without the shock and a model with the shock. The same shock model is used to estimate the relationship between the import weight series.

Demand Model
The theoretical background behind the import demand for forest products is explained in [12], based on the general formulation by [13]. The simplest formulation of a macroeconomic demand model relates Y rt , the import demand of tropical timber within one country at time t to the price of tropical timber P rt (proxied by the unit price of trade) and national revenue proxied by G t , the GDP at time t. It is estimated in its logarithm form (Equation (4)). This static model is sometimes estimated in a dynamic form by adding a lagged version of the demand Y r(t−1) as an additional regressor, but we only estimate the static version here: ln Y rt = β 0 + β 1 ln P rt + β 3 ln G t .
Some authors use simultaneous equations to estimate demand and supply in one model [14,15], but such estimations require additional exogenous variables as determinants of the supply side and are therefore generally done on a national level. International demand models on the other hand [12,16,17] are generally estimated with a single equation, independently of the supply. Thus making the assumption that the supply is infinitely elastic. This is a limitation of our model which we could not address at this stage.
To the simple demand model (Equation (4)), we add an additional explanatory variable representing the price of a substitute or complementary product (Equation (5)). The β 2 coefficient can be interpreted as a cross price elasticity of demand: We estimate the demand model (Equation (5)) from a balanced panel dataset. To capture the unobserved heterogeneity between countries, the estimation is performed with individual country fixed effects. Before performing the estimation, we test for the presence of autocorrelation issues with panel stationarity tests. In the case of a non stationary series, we also estimate the model with differenced variables to remove potential issues with autocorrelation. The problem with results based on differenced variables is that the coefficients cannot be interpreted as elasticities, still, the sign and strength of the coefficient give indications for model comparisons.

Data Preparation
It is preferable to use a narrow product group so that changing prices represent a market signal and not a change of the product mix among a large basket of products. CN codes of temperate sawnwood such as oak have remained the same throughout the whole period. Unfortunately, there has been several changes in tropical hardwood lumber product codes throughout the last 20 years and it is not possible to find a stable tropical timber product code at the CN8 level because product descriptions change through time. For example in 2007 product code 444072969 was split into three different codes (44072799, 44072899, and 44072968) described Table A1 in Appendix A. While in 2017, product code 44072960 was merged into 44072995. These changes are visible Figure A2 in the subplot showing the EU's imports from Africa. To build a long time series of tropical hardwood lumber, we selected products related to 44072969, which is the sawnwood product with the highest import value over the period. A correspondance table from EUROSTAT describes changes in the Combined Nomenclature of the EU's external trade statistics [18]. Based on that table, we built a procedure that looks forward in time for all products related to 44072969 and then backward for all products that have been related to these products in the past. From that list we excluded the products codes that are not specific to tropical timber. In particular we excluded product codes starting with 440799. These products are visible in a separate pane as shown in Figure A2. We exclude these product codes because they are not specific to tropical timber and represent mixes of timber of other origins. The procedure leads to an aggregated product category that represents the vast majority of tropical hardwood lumber export in most countries ( Figure A1). However there are exceptions such as Malaysia and Indonesia where the chosen product aggregation represents only a fourth of the trade value exported to the EU.
The shock model relates the bilateral import price (and weight) series of tropical timber to the import price (and weight) series of oak. The trade partner chosen for oak lumber imports is the USA, for most EU member states the most important extra-EU supplier. A monthly time series was used to increase the power of the estimation. The estimation techniques require a full time series, but many country pairs do not have a complete series. The detail number of missing data for all country pairs are visible in Table A2. Countries with low trade volumes (Table A3) are more likely to have missing data points. We selected country pairs which have zero missing values, in addition we performed linear interpolation for country pairs which had few missing values (less than 20 missing values out of 228 observations). We omitted the other country pairs because they did not have enough observations. Figure A3 illustrates the evolution of monthly prices, trade values, and weight for some of the largest EU tropical hardwood lumber importers from Cameroon. In Germany and Italy, the import weight and prices of tropical hardwood lumber have the same order of magnitude as the import weight and prices of oak lumber, for this reason the comparison makes most sense in those country pairs. France on the other hand imports very low volumes of oak from the USA, mainly because it is itself a large oak producer. Spain's imports of tropical and oak lumber have strongly decreased after the 2008 financial crisis and never recovered since then. The tropical hardwood lumber import series for Malaysia however show much lower import values compared to oak lumber imports in Germany, Italy, and the UK ( Figure A4). Note also the presence of outliers in the Netherlands, as well as the drastic decrease of lumber exports to Belgium. The variable pricew visible in those figures is the ratio of the trade value and weight. This trade unit value is used as a price proxy in the econometric models. It is also possible to compute a price based on the import quantity in cubic metres, but the quantity has more missing data or outliers. That is why we prefer to use the price per units of weight.
The demand model relates aggregate import demand to three variables: The import price of tropical timber, the import price of temperate timber, and GDP. GDP is expressed in constant prices of 2015 based on the chain linked volume method [19]. Monthly GDP data are not available, so we use yearly data. In addition, yearly data smooths seasonal variations which could affect the topical and temperate markets in different ways but the use of yearly data reduces the number of observations from 228 to 19 which negatively impacts the power of the estimation. A further restriction is imposed by the fact that the estimation methods require a balanced panel i.e., observations should be present in all years for all variables and all countries. We removed Cyprus and Malta as they have only 16 years of data and they have the lowest tropical timber import trade value of all countries. We started the series in 2002 to obtain a balanced panel with the largest cross section (N) dimension and the longest time dimension (T) possible. Descriptive statistics of the panel variables used in the demand model are visible in Table A6. Imported weights of tropical hardwood lumber are higher than those of oak in Belgium and France, while they are of a similar order of magnituted in large importing countries such as Germany and Italy ( Figure A6), most of the remaining countries import higher oak quantities than tropical hardwood. All import prices follow increasing trends ( Figure A7) with tropical prices generally higher than oak prices.

Shock Model Estimated on Time Series
Augmented Dickey Fuller (ADF) tests were performed on the period before the shock data to avoid the potential influence of the shock on the stationarity result. Tables A4 and A5 show that price series are non stationary (I(1)) for most country pairs. Therefore the use of cointegration methods is relevant to avoid spurious regression.
In the fifth column of Table 1 (and Table 2), the cointegration statistics are below the critical value of −3.37 for 28 country pairs (37 country pairs). For those country pairs we reject the null hypothesis that the residuals have unit roots, i.e., the price (weight) series of tropical and oak lumber are cointegrated. Results for the country pair, the Netherlands and Malaysia are not reliable because of the presence of outliers (see Figure A4).  For the country pairs Belgium-Ghana, France-Malaysia, France-Côte d'Ivoire, France-Ghana, Netherlands-Côte d'Ivoire, and Spain-Brazil, λ is not significantly different from zero, meaning that there is no significant change in the price relationship after the EUTR start date. In 22 of the cointegrated price series, the lambda coefficients are significant, and of these all but two, Belgium-Côte d'Ivoire and Netherlands-Malaysia, are positive (Table 1). This implies, with these two exceptions, increasing prices for tropical hardwood lumber following the introduction of the EUTR, ranging from 4% (Portugal-Cameroon) to 83% (Germany-Malaysia).
The results of a shock model based on a time series of import weights (Table 2)show more cases where λ is not significantly different from zero (for the country pairs Belgium-Congo DRC, Belgium-Malaysia, Denmark-Brazil, France-Gabon, Greec e-Côte d'Ivoire, Ireland-Cameroon, Italy-Congo, Italy-Gabon, Netherlands-Brazil, Netherlands-Cameroon, Spain-Brazil, Spain-Cameroon, and Utd. Kingdom-Cameroon). Eight country pairs-of which four Belgian import partners, all important in terms of trade volumes-have significant and positive lambdas, indicating increasing imports of tropical hardwood lumber relative to temperate lumber imports after the shock. The remaining 16 country pairs have significant and negative lambdas, implying decreasing imports of tropical lumber relative to oak lumber imports after the introduction of the EUTR. In percentage terms, German imports of tropical hardwood lumber from Malaysia and French imports from Côte d'Ivoire show notable reductions. In quantity terms, however, there are many other tropical hardwood lumber import flows that show stronger reductions after the shock, such as, e.g., German imports from Cameroon, French imports from Cameroon and Ghana, and Portuguese imports from Brazil and Cameroon. Figure A5 illustrates how a negative lambda captures some of the decrease in residuals after the EUTR start date, leading to a change of model from the dotted blue curve to the red curve (particularly visible in the quadrants of Germany for example). VPA member countries (third column in Tables 1 and 2) do not exhibit any clear pattern that would distinguish them from non-VPA countries.
Hence, our results from the shock model show considerable consistency with those of Prestemon [6], as there are general increases in the prices and decreases in the import quantities of tropical timber relative to temperate timber after the implementation of a policy instrument aimed at curbing the trade in illegal timber. This pattern is consistent with a reduction of the supply of tropical timber to these EU member states. There are however notable exceptions of coincedental decreasing prices and increases in import quantities of tropical sawnwod, which would be consistent with an increase in the supply of tropical timber in these cases.

Demand Model Estimated on Panel Data
Before performing the estimation of the demand models, we tested the stationarity of the individual time series and of the panel as a whole. Augmented Dickey-Fuller tests fail to reject the null hypothesis of stationarity for most time series in most countries (Table A7). Panel stationarity tests reported in Table A8 are less conclusive. The Maddala and Wu and Levin-Lin-Chu tests both assume a null hypothesis that some of the series in the panel have a unit root. Both tests reject the null hypothesis for the explained variable, i.e., import demand can be considered stationary. The Maddala and Wu test rejects the null hypothesis for both price variables but the Levin-Lin-Chu weakly fails to reject the null. The Hadri test rejects the null hypothesis of stationarity for all variables. It should be noted that stationarity tests can be biased in the presence of structural breaks. In the absence of a clear answer in favor of stationarity, we have estimated the model with variables in level and with differenced variables. Estimated coefficients (Table 3) are reported with clustered standard errors to adjust for the fact that in a panel fixed effect model, the composite error term is correlated over time [20]. Table 3 shows, first of all, that the own price elasticity of import demand is negative, as expected by the theory. It is however only significant, and elastic, which is considerably higher than the elasticity of −0.7 estimated by Turner and Buongiorno [12], over the whole period and the first time period, becoming insignificant for the period following the EUTR implementation. Second, the cross-price elasticity with respect to oak lumber is never significant, and the same applies to the income elasticity of demand. In general, the demand models fit poorly with the explanatory variables, a consequence of the important inter-country variability. We also estimated an oak import demand model ( Table 4). The own price elasticity is negative, again in accordance with theory, and stable through time. A price elasticity equal to −0.7 is similar to the value obtained in the fixed effect estimation of a static import demand model by Turner and Buongiorno [12]. Clearly the own price elasticity of oak import demand does not change over time. Furthermore, the cross-price elasticity with respect to tropical hardwood lumber is not significant over the whole time period and before the EUTR start date. It is significant and positive when estimated over the 2014-2019 period, indicating that tropical and oak lumber are substitutes, but it is not significant when estimated on differenced variables. The assessment as to whether tropical and oak lumber are substitute or complementary products is more plausible when their import volumes are of the same order of magnitude. This is true for many countries, but this is not the case in France which has a large domestic oak production and as a consequence imports small amounts of oak at higher prices ( Figure A3). An indication of the ratio of tropical to oak lumber import is given in the shock models results (Tables 1 and 2). In general, we observe that the import demand for oak lumber behaves similarly to a more general sawnwood demand import model [12].

Summary and Conclusions
The EU Timber Regulation (EUTR) is a key element of the Forest Law Enforcement, Governance and Trade (FLEGT) Action Plan aimed at curbing illegal logging and its associated trade. Due to the rather short time that the EUTR has been in force, there are yet no studies that have performed systematic, comprehensive statistical analysis as to impacts on international trade in timber products. This study aims to help fill this research gap by assessing whether the patterns of EU imports of hardwood timber products have changed as a result of the coming into force of the EUTR. In addition to the analysis itself, to the best of our knowledge the first attempt at a comprehensive statistical analysis of the effects of the EUTR, a major contribution of the paper is the procedure for tackling the issue of the frequent changes in tropical hardwood lumber product codes, in order to build as long and homogenous a time series as possible.
The cointegration intervention, or shock, models analysis, indicate-for most of the bilateral trade flows (EU member states imports) where we observe cointegration and a significant shock variable-that increasing import prices and decreasing import quantities of tropical hardwood lumber relative oak lumber. This is consistent with a predominant reduction of the supply of tropical timber in these cases. There are however notable exceptions from this pattern, in terms of coincedental decreases of prices and increases in import quantities of tropical sawnwod, which would instead be consistent with an increase in the supply of tropical timber in these cases. In addition, the model does not show any clear difference in results between VPA and non-VPA countries.
The results of the import demand models do not give a clear indication as to whether oak lumber is a complementary or substitute product to/for tropical hardwood sawnwood, and they do not give an unanimous picture as to any potential structural changes in demand elasticities after the entry into effect of the EUTR. It should be noted that price variables used in the demand model consider that the supply is infinitely elastic [12], but this assumption is a tad unrealistic, since the EUTR could be seen to introduce supply constraints, in the sense of costs of complying with the requirements of the legislation. An extension of this work could analyze the supply side of the market.
We estimate models before and after the EUTR start date respectively. However, changes in the market due to EUTR legislation are likely to happen gradually rather than discretely. Hence, market operators may have anticipated the introduction of the EUTR, as suggested by some studies [21]. Furthermore, the first legal cases of enforcement occurred several years after the EUTR entered into force. A growing body of legal proceedings on EUTR enforcement is likely to increase the dissuasive effect of the legislation, and to increase compliance of due diligence. These processes are ongoing, and their effects will be measurable over the next decade. It is therefore likely that the structural changes in demand and supply functions are not abrupt, but gradual, following a smooth transition from one model to a new post-EUTR model.
Even assuming that we have captured the most pertinent effects of the EUTR on EU imports of tropical hardwood sawnwood, we are not in a position to assess whether the regulation has been successful in curbing illegal logging and its associated trade. Thus, reduced supply of tropical timber to EU member states might indicate that exporters have seized exporting illegal timber to these destinations. However, this would not necessarily mean that illegal logging in the exporting nation has decreased, as the illegal timber could have been diverted to other less regulated markets, either abroad and/or domestic.
Further, it is important that the EUTR does not deter the export of legal timber products from developing countries to the EU. Addressing these issues also calls for a study of the supply of tropical timber products, in particular trade diversion, to complement the current analysis.

Acknowledgments:
The work described in this paper has (though not constituting its official output) been carried out in the context of the JRC Biomass assessment study https://ec.europa.eu/ knowledge4policy/projects-activities/jrc-biomass-study_en, accessed on 28 April 2021.

Conflicts of Interest:
The authors declare no conflict of interest. The opinions expressed herein are those of the authors and do not necessarily reflect the views of the European Commission. The scientific output does not imply a policy position of the European Commission.
Sawnwood product group tropical oak Figure A3. Comparison of tropical sawnwood imports from Cameroon and oak sawnwood imports the United States.
Tropical products related to the CN code 44072969 were used in the cointegration analysis.