1 Introduction

In the process of international negotiations that led to the agreement in Kyoto,Footnote 1 fears about possible carbon leakage effectsFootnote 2 from the Protocol were widespread (Barrett 1998). The United States, the major emitter of greenhouse gas at the time of negotiating, motivated its refusal to participate with leakage concerns (Hagel 1997). Models of international trade (Copeland and Taylor 2005) and computable general equilibrium models (Babiker 2005) were used to study the issue and predicted widely varying leakage rates. Predicted leakage was shown to be very sensitive to the choice of parameters (Burniaux and Oliveira Martins 2012).

Studies by Aichele and Felbermayr (2012, 2015) presented the first econometric evidence for carbon leakage via international trade. Their second paper examines leakage on a sectoral level. The authors argue that Kyoto ratification caused an increase of sectoral imports from non-member countries and a decrease of exports to non-members of 7.8%. Their estimation strategy derives from the building blocks of structural gravity models. Trade flows of embodied \({\hbox{CO}}_{2}\) are explained by a variable that captures the between-country difference in Kyoto ratification status. Regressions are conducted on pooled data of sectoral trade flows (40 countries, 12 sectors).

The aim of the present study is to estimate the partial equilibrium carbon leakage caused by the Kyoto Protocol while controlling for multilateral trade resistance.Footnote 3 More specifically, we are asking how Kyoto affected the global carbon trade balance of ratifiers versus non-ratifiers, i.e. the difference between exports and imports of embodied carbon dioxide. The concept of the carbon trade balance is also known in the literature as ‘\({\hbox{CO}}_{2}\) trade balance’ or ‘balance of embodied emissions in trade’ (Munksgaard and Pedersen 2001; Muradian et al. 2002) and was used to study emission transfers via international trade between Kyoto ratifiers and non-ratifiers (Peters et al. 2011). Examining carbon leakage using the concept of the carbon trade balance is a direct consequence of the ’strict’ definition of carbon leakage we employ [see the discussion in Peters and Hertwich (2008)]. Essentially the decision to ratify commitments under the Kyoto Protocol is a unilateral environmental policy even if the negotiating process was multilateral. The decision to ratify engages a country relative to its trading partners. At the same time, the decision to not ratify also amounts to a change of relative domestic environmental policy. The ratifying country will increase relative production costs and therefore substitute foreign for domestic products, while the non-ratifying country will decrease relative production costs and substitute domestic for foreign products. Thus we need to look at both exports and imports of Kyoto ratifiers (and non-ratifiers) to determine the amount of carbon leakage caused by the Protocol.

Fig. 1
figure 1

Carbon Trade Balance of Kyoto ratifiers with non-ratifiers. Notes The carbon trade balance consists of exports minus imports of embodied \(\hbox{CO}_{2}\). This graph is from the viewpoint of Kyoto ratifiers in the sample, where trading partners are all non-ratifiers in the sample

Looking at Fig. 1, we observe a noticeable kink in the development of the carbon trade balance between Kyoto ratifiers and non-ratifiers around the year 2002. Most member countries ratified their commitments in national parliaments that year. Committed countries increased their net imports of embodied carbon dioxide from around 100 to around 600 megatons. Several explanations for this development are possible, among them leakage effects, China joining the WTO,Footnote 4 an increasing number of regional trade agreements between members and non-members of Kyoto and so on. We follow the methodology of Aichele and Felbermayr (2015) of effectively treating the Protocol as influencing domestic production costs which will result in changes to domestic mill prices and consequently to multilateral trade resistance. Under the leakage hypothesis, the Kyoto Protocol leads to an increase in the effective price of carbon dioxide. This in turn will lead to a partial substitution of goods and services from non-ratifiers for those from ratifiers.

In order to meaningfully relate estimated leakage to the carbon trade balance we need to take into account sectoral differentiation of the effects and apply the estimated elasticities to sectoral emission trade flows. This is especially important if trade flows are very unequally distributed across sectors as is the case with services. While several studies have presented evidence on traded emissions including services sectors (Fernández-Amador et al. 2016, is one example), the use of the World Input–Output Database (WIOD) in combination with multi-regional input–output analysis allows us to present the first econometric evidence for leakage effects from the Kyoto Protocol through services trade. Carbon leakage from services turns out to be quite substantial in percentage terms, although less so in absolute terms.

Besides the estimation of Kyoto’s impact on the carbon trade balance, including the first estimates of services’ role in that, this study makes four further contributions. First, we estimate the treatment effect from Kyoto ratification via a two-stage procedure suggested by Head and Mayer (2014) in order to control for multilateral trade resistance while estimating the carbon trade impact of a unilateral policy. Second, we suggest an alternative specification of the Kyoto treatment variable that allows us to estimate asymmetric effects on Kyoto ratifers’ imports and exports. Third, we allow for time-varying effects of ratification. Fourth, we take into account the special role of Eastern Europe, a previously neglected issue. These countries are shown to significantly influence the results, however, we will raise some doubts about the estimates in their case.

The average annual estimated effect of the Kyoto Protocol on the carbon trade balance amounts to 88 megatons of \({\hbox{CO}}_{2}\), summing up to a total of 705 megatons in the period 2002–2009. Kyoto ratification is associated with an increase of 4.4% of total emissions embodied in imports, and 4.2% of total emissions embodied in exports from ratifiers to non-ratifiers. Depending on the years, between 11 and 54% of the negative carbon trade balance of Kyoto countries can thus be explained. We can quantify major sources of carbon leakage in transport sectors such as water or inland transport. The effect of services remains somewhat limited due to the low level of trade in these sectors. Energy-intensive manufacturing sectors such as metals, machinery and transport equipment exhibit the highest levels of leakage and are responsible for 76% of the total estimated effect on the carbon trade balance.

2 Data and construction of the carbon trade balance

The carbon trade balance consists of exports minus imports of carbon dioxide, i.e. flows of emissions embodied in the traded goods and services. To construct a data set of emission trade flows we use the World Input–Output Database (WIOD). WIOD offers annual multiregional input–output (MRIO) tables from 1995 to 2009.Footnote 5 These tables are needed in order to approximate average sectoral \({\hbox{CO}}_{2}\) emission intensities that include the upstream emissions from other sectors and countries.Footnote 6 The world in WIOD consists of 40 countries and one rest of the world category. Each economy is split into 34 sectors, 17 of which are services sectors. At the same harmonized country-sector resolution, this database also provides information on carbon dioxide emissions from fossil fuel combustion and non-energy usage.Footnote 7 This informs us about direct sectoral emissions. We use the MRIO tables to include indirect emissions (emissions from previous production stages) in the calculation of emission intensities of traded goods and services.Footnote 8 Bilateral sectoral trade flows are implicit in the MRIO tables and can be extracted thereof. For more details on the data and some descriptive statistics, consult “Appendix 1”.

The database is well-suited for our analysis due to the inclusion of the major Kyoto committed countries (with Switzerland and Norway being the most prominent missing ones) as well as the most important non-Kyoto countries.Footnote 9 The sample covers 85% of world GDP, slighly less of global carbon dioxide emissions (81–82%) and between 82–87% of total world trade in value terms. WIOD also carefully corrects problems in relation to re-exports and export-processing zones.Footnote 10

The world vector of emission intensities \(\Gamma _t\) with representative element \(\gamma _{j,s,t}\) in sector s of country j in year t is obtained by pre-multiplying the World Leontief inverse with the vector of direct emissions per unit of a sectors output \(R_t\).Footnote 11

$$\begin{aligned} \Gamma _t = R_t(I-A_t)^{-1} \end{aligned}$$
(1)

We combine the information in \(\Gamma _t\) with the value of imports \(M_{i,j,s,t}\) of importer i from exporter j to obtain the \({\hbox{CO}}_{2}\) emissions embodied in an import flow \(M^c_{i,j,s,t}\).Footnote 12

$$\begin{aligned} M^c_{i,j,s,t} = M_{i,j,s,t} \gamma _{j,s,t} \end{aligned}$$
(2)

We define the carbon trade balance \(CTB_{i,j,t}\) of country i with country j as the sum of sectoral exports \(E^c_{i,j,s,t} = M^c_{j,i,s,t}\) minus imports \(M^c_{i,j,s,t}\) of carbon dioxide.

$$\begin{aligned} CTB_{i,j,t} = \sum _{s \in S} M^c_{j,i,s,t} - \sum _{s \in S} M^c_{i,j,s,t} \end{aligned}$$
(3)
Fig. 2
figure 2

Country carbon trade balances: Largest \({\hbox{CO}}_{2}\) traders. Notes. Vertical axis: Megatons \({\hbox{CO}}_{2}\). Each panel shows the eight largest carbon dioxide traders among Kyoto ratifiers (upper panel) and non-ratifiers (lower panel). Each country is shown twice, its carbon trade balance with the group of ratifiers and non-ratifiers

In Fig. 1 we presented the carbon trade balance for all Kyoto countries with all non-Kyoto countries in WIOD. The picture would look mirror-inverted for non-Kyoto countries. Figure 2 depicts the carbon trade balances for individual countries. The core observation from the upper panel of these graphs is that the carbon trade balance with non-members decreased for all of the most important Kyoto ratifiers after 2002.Footnote 13 At the same time the carbon trade balance with other member countries does not show a clear pattern: decreasing for some (RUS, CAN), increasing for some (DEU, JPN, ITA), relatively constant for some (FRA, GBR, NLD). For all of these countries, except Russia, the carbon trade balance with other Kyoto members developed more positively compared to non-members. However, the difference between the development of member and non-member trade differs considerably across countries. The leakage hypothesis is consistent with developments like in Germany: The carbon trade balance with Kyoto non-members develops more negatively than with members.

For the most important non-ratifiers, we observe that on average the carbon trade balance increased by more for trade with Kyoto members. Interestingly, the largest emitter of carbon dioxide, China, exhibits no difference in the development of Kyoto and non-Kyoto trading patterns. The second largest emitter, the United States, displays the pattern we would expect under the leakage hypothesis: The carbon trade balance with Kyoto members develops more positively than with non-members.

3 Methodology

3.1 Estimating emission trade elasticities

In order to identify the effect of the Kyoto Protocol on the carbon trade balance we follow and extend the methodolgy proposed by Aichele and Felbermayr (2015). The logarithmic transformation of bilateral \({\hbox{CO}}_{2}\) trade flows \(M^c_{i,j,s,t}\) will act as the dependent variable in our regressions. We explain them by a measure of differential Kyoto treatment and control for time-varying pair-specific policy variables \(P_{i,j,t}\) (common membership in trade agreements, the European Union and the World Trade Organization). We control for country-pair idiosyncratic sectoral carbon trade levels \(\nu _{i,j,s}\), i.e. capturing path dependencies of particularly strong sectoral relations that are time-independent within our sample (Baltagi et al. 2003). Additionally, and in contrast to Aichele and Felbermayr, we control for multilateral resistance by including importer and exporter-specific time effects on carbon trade \(\theta _{i,s,t}\) and \(\zeta _{j,s,t}\). Thereby we capture if a country increases or decreases \({\hbox{CO}}_{2}\) imports or exports with all other countries in the sample over time.

Treatment (Kyoto participation) of importer i is modeled using an indicator variable \(K_{i,t}\) taking the value 1 if i has ratified binding emission constraints under Kyoto in year t. We propose two ways of including the difference in bilateral Kyoto status in our regression. Specification (i) uses the Aichele and Felbermayr formulation \(\Delta K_{i,j,t} = K_{i,t} - K_{j,t}\). Here, differential Kyoto treatment \(\Delta K_{i,j,t}\) can take the values \(\{-1,0,1\}\), with 1 if only the importer is treated and \(-1\) if only the exporter is treated. In specification (ii), we include two indicator variables: \(\overline{K^m}_{i,j,t}\) has value 1 whenever only the importer is treated, \(\overline{K^x}_{i,j,t}\) has value 1 whenever only the exporter is treated.

In a structural gravity framework one way to control for the multilateral resistance terms lies in the inclusion of importer- and exporter-specific time fixed effects \(\theta _{i,s,t}\) and \(\zeta _{j,s,t}\). As is well known in the gravity literature, the effects of unilateral policies cannot be directly estimated in most structural gravity models as importer and exporter-specific time fixed effects would wipe out all the variation with respect to the policy of interest (Yotov et al. 2016). From the viewpoint of each country the ratification of commitments under the Kyoto Protocol is a unilateral decision not varying over trading partners. This would lead to perfect collinearity between \(K_{i,t}\) and \(\theta _{i,t}\). Head and Mayer (2014) suggest a two-stage approach to estimate unilateral effects while controlling for multilateral resistance. In the first stage, a regression of logarithmic trade flows on all bilateral policy variables, standard trade cost variables and multilateral resistance terms is estimated. Standard trade costs (e.g. distance, contiguity, language, cross-border trade, colonial history) are captured by \(\nu _{i,j,s}\) and multilateral resistance is captured by \(\theta _{i,s,t}\) and \(\zeta _{j,s,t}\). In the second stage, the estimates \(\widehat{\theta _{i,s,t}}\) and \(\widehat{\zeta _{j,s,t}}\) from the first stage are regressed on all the it- and jt-specific variables including controls for importer carbon expenditure levels \(e_{i,s,t}\) and exporter carbon output levels \(y_{j,s,t}\).Footnote 14

$$\begin{aligned}&{\text{(1st stage) }} \nonumber \\&ln M^c_{i,j,s,t} = \alpha _0 + \delta P_{i,j,t} + \theta _{i,s,t} + \zeta _{j,s,t} + \nu _{i,j,s} + \epsilon _{i,j,s,t} \end{aligned}$$
(4)
$$\begin{aligned}&{\text{(2nd stage) }} \nonumber \\&\text {(i) }\widehat{\theta _{i,s,t}} + \widehat{\theta _{j,s,t}} = \alpha + \varvec{\sum _{g \in G} \sum _{s \in S} \sum _{t \in T} \left[ \beta _{g,\,s,\,t} \Delta K_{i,\,j,\,t} \right] } + \alpha _e ln e_{i,s,t} + \alpha _y ln y_{j,s,t} + \varepsilon _{i,j,s,t} \end{aligned}$$
(5)
$$\begin{aligned}&{\text {(ii) }}\widehat{\theta _{i,s,t}} + \widehat{\theta _{j,s,t}} = \dots + \varvec{\sum _{g \in G} \sum _{s \in S} \sum _{t \in T} \left[ \beta ^m_{g,\,s,\,t} \overline{K^m}_{i,\,j,\,t} + \beta ^x_{g,\,s,\,t} \overline{K^x}_{i,\,j,\,t} \right] } + \dots \end{aligned}$$
(6)

The motivation for specification (ii) is that we want to allow for asymmetric treatment effects on treated countries’ imports and exports. Specification (i) amounts to imposing the restriction \(\beta ^m_{g,s,t} = - \beta ^x_{g,s,t}\). This restriction might not be warranted a priori. As an example take a country that raises the environmental standard on some product group to a level that domestic producers would already fulfill but foreign producers would not. In this case effective trade costs for importing would increase while they would stay unaffected for exporting.

Generally, we allow treatment effects \(\beta _{g,s,t}\) to vary over time, sectors and countries/country groups. The set T contains all years between 2001 (first countries ratified Kyoto) and 2009 (last year in our sample). S contains all 34 sectors in the sample. G contains three country groups for which we expect differing effects (see discussion in Sect. 5). In addition we also report estimates with the restriction \(\beta _{g,s,t} = \bar{\beta }\)\(\forall g,s,t\) to get the average sectoral time-constant effect of Kyoto, as well as the estimates restricted by \(\beta _{g,s,t} = \tilde{\beta _t}\)\(\forall g,s\) to get the average sectoral treatment effect over time.

Following the tradition of the gravity literature, we estimate Eqs. 4 and 5 by pooled OLS. We perform several robustness checks in Sect. 7. With respect to omitted variable bias and consistency of the estimated treatment effect we trust that the included fixed effects take care of most of the unexplained variation in emission trade.Footnote 15 The cross-sectional importer–exporter-sector fixed effect \(\nu _{i,j,s}\) controls for time-invariant bilateral trade relations and thereby for the specialization of each economy which might affect selection into ratifying.Footnote 16

For a causal interpretation of the Kyoto effect, we need the common trend assumption: In the absence of Kyoto, treated country pairs would behave as their untreated counterparts. Some of the most important non-ratifiers are emerging economies (EECs) like China, India, Brazil, Mexico, Turkey and Indonesia, while most of the ratifiers are high-income economies (HECs). Sectoral trade growth of the EECs with all other countries is controlled for by \(\theta _{i,s,t}\) and \(\zeta _{j,s,t}\). Asymmetric trade growth of the EECs with the HECs might cause a bias to our estimate if it is timely correlated with Kyoto ratification. However, as Aichele and Felbermayr (2015) point out, the construction of their treatment variable \(\Delta K_{i,j,t}\) mitigates this potential bias. If asymmetric trade growth between EECs and HECs were solely the result of non-Kyoto related factors, this would bias \(\widehat{\beta }_{g,s,t}\) downwards in specification (i).Footnote 17 This is not necessarily true for \(\widehat{\beta ^m}_{g,s,t}\) and \(\widehat{\beta ^x}_{g,s,t}\) in specification (ii).Footnote 18 Thus specification (i) is preferable if we have reason to believe that trade growth between Kyoto ratifiers and non-ratifiers behaves significantly different from trade growth within each group (ratifiers, non-ratifiers) even in the absence of the Protocol.

3.2 Calculate effects on the carbon trade balance

Ultimately we are interested in the change of the carbon trade balance between Kyoto ratifiers (rat) and non-ratifiers (non). In the style of Eq. 3, we define

$$\begin{aligned} CTB_{rat,non,t}&= \sum _{i \in non} \sum _{j \in rat} \sum _{s \in S} M^c_{j,i,s,t} - \sum _{i \in rat} \sum _{j \in non} \sum _{s \in S} M^c_{i,j,s,t} \end{aligned}$$
(7)
$$\begin{aligned}&= \sum _{s \in S} M^c_{non,rat,s,t} - \sum _{s \in S} M^c_{rat,non,s,t} \end{aligned}$$
(8)

The estimated change in the carbon trade balance follows from applying the elasticities to actual emission trade flows.

$$\begin{aligned} \widehat{\Delta CTB_{rat,non,t}} = \sum _{g \in G} \sum _{s \in S} M^c_{non,g,s,t} \widehat{\beta ^x}_{g,s,t} - \sum _{g \in G} \sum _{s \in S} M^c_{g,non,s,t} \widehat{\beta ^m}_{g,s,t} \end{aligned}$$
(9)

When estimating specification (i) we would just use \(-\widehat{\beta ^x_{g,s,t}} = \widehat{\beta ^m_{g,s,t}} = \widehat{\beta _{g,s,t}}\).

4 Restricted estimates of Kyoto treatment

Estimating Eqs. 4 and 5 generates a large amount of output. The most important results are presented in Sect. 6. In the current section we discuss some of the results under relatively strict restrictions on the coefficient estimates. We start by looking at the coefficients under the constraint that one and the same estimate is applying to the whole sample.

Table 1 Kyoto treatment effects: time-, sector- and group-invariant
Table 2 Kyoto treatment effects: time-varying

Table 1 presents the estimates for time-, sector- and group-invariant Kyoto treatment effects. Columns 1–3 display results for specification (i), while columns 4–6 display results for specification (ii). For the full sample and specification (i), Kyoto ratification implies on average increasing sectoral imports and decreasing sectoral exports of 11.2%. Aichele and Felbermayr (2015) find a baseline estimate of 7.8%. What explains the difference? The latter study uses only twelve sectors, most of them manufacturing sectors. It also uses a slightly different sample of countries and a slightly shorter time span. If we exclude all tertiary sectors (services sectors 18–34) from the estimation the estimate is reduced to 6.4% (see column 2). The result from the pooled sample in column 1 is thus considerably driven by services sectors. Looking at column 4, specification (ii) reveals that on average the effect on exports is larger than that on imports. This is driven by primary (agriculture, mining) and secondary sectors (manufacturing), where imports from ratifiers increase by only 3.1% while exports to ratifiers decrease by 9.7%. For tertiary sectors, the difference between imports and exports is not statistically significant at the 5% level as the F-value of 0.85 confirms.

In Table 2 we examine the time-dimension of the effects. We can reject the hypothesis that coefficients are the same for all years at the 5% level. In specification (i), variation is relatively limited and there is no clear trend of growing or shrinking influence. Interestingly, the same is not true for specification (ii). The coefficients are presented graphically in Fig. 3. In the first years after ratification the negative effect on exports seems dominant. The F-tests in Table 2 inform us that for most of the years, effects on imports and exports are indeed significantly different. In the years after 2007 importer effects become dominant.

Fig. 3
figure 3

Kyoto treatment effects: Time-varying. Notes. Coefficients from Table 2. Coefficient for \(\beta ^x_t\) is sign-inverted

A further refinement of the estimates concerns the role of country groups, specifically the role of the Eastern European countries. In our main results we will distinguish estimates for three groups of Kyoto ratifying countries: old EU-15, new EU-10 from Eastern Europe and Non EU countries. In the following section we will discuss the motivation for doing so.

5 The role of Eastern Europe

The fact that Kyoto ratification meant widely varying constraints for individual countries is depicted in Fig. 4. It compares the emission target for each country with its actual emission growth rates before ratification. To the left of the graph we find countries facing relatively strict targets like Canada, Spain and Portugal. The more we move to the right end of the graph, the less binding are the emission targets. We should underline two main points here: First, the distribution of emission targets is much narrower than the distribution of the actual emission growth rates before ratification. Second, especially the Eastern European countries faced targets that meant almost no restrictions on their emission path. All countries with a particularly large negative difference of at least two percentage points are from that region. Even if we account for the fact that pre-ratification emissions are an imperfect measure of expected future business-as-usual emissions, the size of the differences is striking.

Fig. 4
figure 4

Country emission growth targets under Kyoto. Notes. For each country, the darker left bar shows the necessary average annual growth rate of emissions from 1990–2010 to stay below the binding emission target. The right bar shows the actual average annual growth rate of emissions from 1990 until the ratification of Kyoto. The difference (actual–target) indicates how strict a country’s target is, i.e. by how much would a country have to reduce its pre-ratification emission growth rate in order to achieve the target

Of the ten countries with a negative trend-versus-target difference of more than two percentage points, nine joined the European Union either in 2004 (Poland, Hungary, Czech Republic, Slovakia, Lithuania, Estonia, Latvia) or 2007 (Bulgaria, Romania). The only missing countries from this wave of EU enlargement are Cyprus and Malta—both of them without an emission target under Kyoto—, as well as Slovenia which faced a stricter Kyoto target being also the economically most developed Eastern European nation.

These observations lead us to allow for varying estimates of the Kyoto treatment effect over country groups. One group (NewEU) consists of the ten Eastern European EU accession countries with emission targets. We split the rest of the Kyoto committed countries into the group of the fifteen members of the European Union prior to 2004 (OldEU) and a third group that combines the Kyoto ratifiers outside of the European Union (NonEU), i.e. Russia, Canada, Australia and Japan. For each group we create an indicator variable that takes value 1 whenever a group member is either the importer or the exporter. We interact these dummies with our treatment variable to estimate the average sectoral time-constant effect of Kyoto for each country group. In terms of Eqs. 4 and 5 this amounts to imposing a restriction \(\beta _{g,s,t} = \check{\beta }_g\)\(\forall s,t\). A priori we would expect to find low leakage estimates for countries with non-binding targets. These countries would not be forced to undertake particularly strong measures in order to change their emission path.

Table 3 Kyoto treatment effects by country groups

Column 1 in Table 3 shows that, contrary to that expectation, the new Eastern European EU members exhibited the strongest effect on the baseline leakage estimate. The carbon trade balance of these countries is affected by 20.1%. Specification (ii) suggests a highly significant differentiation of this effect: We estimate that emission imports from non-Kyoto countries increased by 32.5% after the ratification while exports declined by only 8.7%. On the other hand, treatment effects for the Old- and Non-EU groups are estimated at 2.7 and 4.8% respectively and seem to be the result of decreases of emission exports with slightly negative leakage on the import side. Results are similar if we conduct the same exercise with the data used by Aichele and Felbermayr (2015). The treatment effect in the Old EU and Non-EU group is estimated at 5 respectively 3% while Eastern European countries’ carbon trade balance is affected by 16.5%. Again the major impact on Eastern European countries is an increase in \({\hbox{CO}}_{2}\) imports.

As Eastern Europe is found to be the main driver of leakage in our sample we should expect that a strong increase in domestic environmental regulation took place in Eastern Europe in the beginning of the 2000s. There is evidence for such a boost in regulation, however, the driving force behind this boost is more likely the accession to the European Union instead of the ratification of Kyoto. As we observed, Kyoto did not incentivize immediate environmental regulation for Eastern Europeans. On the other hand, environmental regulations as specified in the acquis communautaire were a prerequisite for joining the EU. As we proceed to estimate the disaggregated effects of Kyoto in the next section, we must control for the specific effects of Kyoto in Eastern Europe. Restricting \(\beta _{g,s,t}\) to be constant across groups would lead to an overestimation of the Kyoto effect for the Old EU and Non-EU countries, which would be particularly problematic given that these countries exhibit much higher trade volumes with the non-ratifiers. This would bias the magnitude of the leakage estimates upwards.

6 Kyoto’s impact on the carbon trade balance

We use Eq. 9 to calculate the effect of Kyoto ratification on the carbon trade balance between Kyoto ratifiers and non-ratifiers. The results presented here are obtained using specification (i) of the Kyoto treatment variable. This is our preferred specification as we might have reason to believe that trade between non-ratifiers and ratifiers would have grown by more than trade growth within each group due to long-run changes in specialization patterns.Footnote 19,Footnote 20 We use a significance level of 1% for the treatment effect. Figure 5 presents the overall development of the carbon trade balance between Kyoto ratifiers and non-ratifiers, indicating the estimated effect that the Kyoto Protocol exerted. Committed countries carbon trade deficit rises from 101.2 megatons in 2002 to 607.2 megatons in 2008. Kyoto is estimated to have affected the carbon trade balance negatively, with an absolute value of 76.9 megatons in 2002 and 152.9 megatons in 2008. Figure 5 also indicates that the dominant part of this effect originates from the Old EU and Non-EU countries. The Old EU countries experience the strongest leakage in the years after most members joined the Protocol between 2003 and 2005. Similarly the non-EU group experiences strong leakage after 2006, when Russia and Australia ratified the Protocol. The New EU countries display comparatively modest levels of net leakage across all years with around 14% of total leakage.

Fig. 5
figure 5

Effect of the Kyoto Protocol on the carbon trade balance, by country groups. Notes. The carbon trade balance of Kyoto ratifiers with non-ratifiers. Specification (i). Significance level 1%. The effects for three country groups of ratifiers (‘Non-EU’, ‘Old EU’ and ‘New EU’) are stacked

Fig. 6
figure 6

Effect of the Kyoto Protocol on the carbon trade balance, by sectors. Notes. The carbon trade balance of Kyoto ratifiers with non-ratifiers. Specification (i). Significance level 1%. The effects for different sectors (and groups of sectors) are stacked. Group ‘Other Manufacturing/Agriculture’: sectors 1–11, 14–17. Group ‘Services’: sectors 18–34. See sector list in Table 11, “Appendix 1

Figure 6 disaggregates these net leakage estimates according to sectors. The two sectors with the single largest impact on the overall effect are machinery (13) and basic and fabricated metals (12). Leakage from the machinery sector explains 30.5% of total net leakage between 2002 and 2009. Leakage from the basic and fabricated metals sector explains 26.4% of total net leakage between 2002 and 2009. A third substantial source of leakage is the transport equipment sector (15). The remaining primary and secondary sectors taken together exhibit a neutral effect on net leakage with -0.9% of the total. Service sectors (18–34) play a minor role compared to manufacturing, although exhibiting substantial leakage considering their overall lower volume of trade. Their share of net leakage lies at 24.7% of the total in the period 2002–2009, most of which is accounted for by the transport sectors (inland, water and air) as well as the renting of machinery and equipment sector.

Table 4 Kyoto treatment effects: summary table

Table 4 summarizes the main results from the viewpoint of Kyoto ratifiers. The Kyoto-induced increase in imports from non-committed countries lies at 4.4% of total carbon imports while the corresponding decrease of exports amounts to 4.2% of total carbon exports. Summing up exports and imports of carbon dioxide and comparing them with leakage yields an Kyoto leakage share in total carbon dioxide trade of 4.3% in the period 2002–2009, varying between 2.4% (2006/07) and 5.2% (2008/09). The corresponding Kyoto leakage share in total \({\hbox{CO}}_{2}\) emissions varies between 0.1% (2006) and 0.6% (2008), with a mean of 0.4%.

The sectoral structure of leakage yields some further interesting insights. Figure 7 depicts the country group-specific impact of the Kyoto Protocol on single sectors and sector groups. We present average annual leakage values for two time spans, the early period between 2002 and 2005 (upper panel) and the later period between 2006 and 2009 (lower panel). We observe a relatively constant leakage level on average. Average annual sectoral leakage from the Kyoto Protocol across all sectors decreases from 2.7 megatons for the early period to 2.5 megatons for the later period, a decrease of roughly 10%. At the same time the growth of overall \({\hbox{CO}}_{2}\) trade amounted to 20%.

Fig. 7
figure 7

Group-specific sectoral effects on the carbon trade balance. Notes. Kyoto-induced change in the carbon trade balance of Kyoto ratifiers with non-ratifiers. Specification (i). Group ‘Other Manufacturing’ (oth. manfct.) summarizes sectors 3–7, 9–11, 16–17. Group ‘Other Services’ (oth. serv.) summarizes sectors 18–22, 27–29, 31–34

Three sectors with a predominantly positive net impact on the carbon trade balance, especially in the period after 2005, are sectors 1, 2 and 8. On the one hand the coke, petroleum and petroleum products sector (8), which experiences decreasing imports and increasing exports to non-committed countries, with the latter effect dominating. These effects are especially pronounced in the Old EU countries and the Non-EU countries. On the other hand the primary sectors agriculture (1) and mining (2) also influence the carbon trade balance positively in the later period. This results from export increases in mining in the Non-EU group. Both of these observations would be consistent with a decrease in Kyoto-countries’ domestic demand for inputs of energy-intensive industries.

6.1 Manufacturing sources of carbon leakage

The four major sources of the Kyoto-induced decrease of the carbon trade balance are sectors 11, 12, 13 and 15. We discuss what is driving the results for each of them in the following listing, ordered by size of total leakage (in detail presented in Table 5).

Table 5 Kyoto treatment effects: manufacturing summary table
  • The machinery sector (13) explains 30.5% of total net leakage between 2002 and 2009. Both the Old EU and Non-EU countries are mostly affected. The machinery sector significantly affects the carbon trade balance on both the export and the import side equally in the early period while the import increase dominates in the later period. This is strongly driven by an average annual import increase of 12.5 megatons by the Non-EU group in the later period after 2.8 megatons in the early period. The effect for the New EU group is moderate with an average annual carbon trade balance decrease of 0.8 megatons (early period) respectively 2.3 megatons (later period).

  • The basic metals and fabricated metals sector (12) explains 26.4% of total net leakage between 2002 and 2009. It is the major source of leakage for the Old EU group throughout the 2000s. Emission imports from metals increase by 7 megatons annually in the early period, and 17 megatons in the later period. In contrast, emission exports from the Old EU group decrease by 4 megatons (early period) respectively 7.4 megatons (later period). In comparison moderate levels of leakage from this sector are observed for the New EU group, and even more moderate for the Non EU group.

  • The transport equipment sector (15) explains 19.4% of total net leakage between 2002 and 2009. The negative effect on the carbon trade balance in this sector stems almost entirely from the Non-EU group. In the early period, average annual emission imports increase by 2.7 megatons, while they increase by 12.7 megatons in the later period. Conversely, average annual emission exports decrease by 4.2 megatons in the early period and by 11.7 megatons in the later period.

  • The non-metallic mineral products sector (11) explains 6.3% of total net leakage between 2002 and 2009. All country groups experience negative impacts on their carbon trade balance. In the early period, average annual net emission imports increase by 2.7 megatons for the Old EU, 2.1 megatons for the Non-EU and 0.9 megatons for the New EU group. In the later period, the corresponding annual decrease in the carbon trade balance amounts to 2.6 megatons for the Old EU, 1.3 megatons for the Non-EU and 1.6 megatons for the New EU group.

With respect to the other manufacturing sectors sizable leakage effects in the early period are found in the electrical and optical equipment sector (14), the chemicals and chemical product sector (9), the manufacturing n.e.c. sector (16) as well as the textiles and textile products sector (4). In the later period many of these sectors experience a reversal, e.g. the electrical and optical equipment sector, with leakage coming to a halt or even reversing towards a positive effect on the carbon trade balance.

6.2 Services sources of carbon leakage

The two most important sources of leakage in services are sectors 23, 24, 25 and 30, which are presented in detail in Table 6.

Table 6 Kyoto treatment effects: services summary table
  • The renting of machinery and equipment sector (30) explains 9.1% of total net leakage between 2002 and 2009. The most important source region of leakage is the Old EU group, and here in particular the emission import side of the carbon trade balance. In the early period, average annual emission imports increase by 2.8 megatons while exports decrease by 1.6 megatons. In the later period, average annual emission imports increase by 7.4 megatons while exports decrease by 2.8 megatons. The remaining change in the carbon trade balance is accounted for by the New EU group.

  • The inland transport sector (23) explains 6.8% of total net leakage between 2002 and 2009. Average annual emission imports increase by 3.2 megatons in the early period and 7 megatons in the later period for the group of Old EU countries. Emission exports stay comparatively relatively unaffected. Again the groups of Non-EU and New EU countries play (almost) no role.

  • The water transport sector (24) explains 4.4% of total net leakage between 2002 and 2009. The picture is very similar to the inland transport sector. Almost the entire leakage effect is due to the Old EU countries increasing their average annual emission imports, in this case by 3.1 megatons (early period) respectively 3.3 megatons (later period).

  • The air transport sector (25) explains 2.1% of total net leakage between 2002 and 2009. Leakage is mostly restricted to the early period, when the Old EU group increased its average annual emission imports by 2 megatons and decreased its average annual emission exports by 1.5 megatons.

7 Robustness checks

To test the robustness of our methodological approach we perform several checks and compare the resulting restricted estimates on the Kyoto treatment variable. Table 7 summarizes the time-, sector- and group-invariant estimates on the Kyoto effect while Table 8 tests robustness for the country group-specific estimates. The baseline model proves to be robust to methodological changes. The Kyoto effect is found to be of very similar magnitude in the extended Aichele/Felbermayr model (R3) as well as in the Baier/Bergstrand (R4) model.

7.1 “Naive” gravity

We estimate a model without controlling for the unobserved fraction of bilateral trade costs and without controlling for multilateral resistance, i.e. excluding bilateral, importer- and exporter-specific fixed effects \(\nu _{i,j,s}\), \(\theta _{i,s,t}\) and \(\zeta _{j,s,t}\). Instead we control for bilateral trade costs by including the observable variables distance, contiguity, language, cross-border trade, colonial history, summarized as vector \(T_{i,j}\) in Eq. 10. To capture time trends of emission trade at the sectoral level we allow for a sectoral time trend \(\tau _{s,t}\). We call this model a “naive” gravity model.

$$\begin{aligned} ln M^c_{i,j,s,t} = \bar{\beta }^{na} \Delta K_{i,j,t} + \alpha _e ln e_{i,s,t} + \alpha _y ln y_{j,s,t} + \delta ^{na}_1 P_{i,j,t} + \delta ^{na}_2 T_{i,j} + \tau _{s,t} + \varepsilon ^{na}_{i,j,s,t} \end{aligned}$$
(10)

Row R1 of Table 7 depicts the resulting estimates for \(\bar{\beta }^{na}\). The “naive” model would lead us to conclude falsely that Kyoto resulted in inverse leakage.

Table 7 Robustness checks on \(\bar{\beta }\)

7.2 Aichele/Felbermayr (2015)

We test two specifications based on the work by Aichele and Felbermayr (2015). While the model in Eqs. 46 controls for importer- and exporter-specific sectoral time trends \(\theta _{i,s,t}\) and \(\zeta _{j,s,t}\), the model in Eq. 11 employs country-specific time trends \(\vartheta _{i,t}\) and \(\psi _{j,t}\). The latter capture the mean development of a country’s aggregate carbon trade as both importer and exporter.Footnote 21 While this allows for estimation of the model in one stage, the implication is that multilateral resistance is not fully controlled for. We estimate two versions of Aichele and Felbermayr’s model, the first being their original formulation, the second including additionally importer carbon expenditure levels \(e_{i,s,t}\) and exporter carbon output levels \(y_{j,s,t}\).

$$\begin{aligned} ln M^c_{i,j,s,t} = \bar{\beta }^{AF} \Delta K_{i,j,t} + \alpha _e ln e_{i,s,t} + \alpha _y ln y_{j,s,t} + \delta ^{AF}_1 P_{i,j,t} + \vartheta _{i,t} + \psi _{j,t} + \nu _{i,j,s} + \varepsilon ^{AF}_{i,j,s,t} \end{aligned}$$
(11)

Rows R2 and R3 of Table 7 depict the resulting estimates for \(\bar{\beta }^{AF}\). The estimates under the original Aichele and Felbermayr formulation (R2) are roughly 3 to 4 percentage points above those of the baseline model while including expenditure and output (R3) leads to similar results as in the baseline model. The same observations are true for the model allowing for country group variations of the Kyoto treatment effect, presented in Table 8. The differences between (R2) and the baseline model as well as (R3) suggest that taking into account sectoral differentiation of carbon output and expenditure trends is important.

Table 8 Robustness checks on \(\check{\beta }\)

7.3 Baier/Bergstrand (2009)

A further robustness check employs the methodology of Baier and Bergstrand (2009) in the context of carbon trade. These authors suggested using exogenous multilateral resistance terms \(MRT_{i,j}\) based on observable trade cost characteristics. We model trade costs as being determined by distance, contiguity, common language, colonial history and a dummy capturing within-country trade.Footnote 22

$$\begin{aligned} ln \frac{M^c_{i,j,s,t}}{e_{i,s,t} y_{j,s,t}} = \bar{\beta }^{BB} \Delta K_{i,j,t} + \delta ^{BB}_1 P_{i,j,t} + \delta ^{BB}_2 MRT_{i,j} + \tau _{s,t} + \varepsilon ^{BB}_{i,j,s,t} \end{aligned}$$
(12)

Row R4 of Table 7 depicts the resulting estimates for \(\bar{\beta }^{BB}\). The estimates are very close to the baseline model as well as the extended Aichele/Felbermayer model (R3), albeit with higher estimates for the manufacturing subsample and lower estimates for the service subsample. These same observations are true for the country group-specific estimates, see the estimates in Table 8.

8 Conclusion

In this study, we estimated carbon leakage as a consequence of the Kyoto Protocol by looking at the effects of ratification on the carbon trade balance. We argued that taking into account regional and sectoral disaggregation is essential in order to deal with very unequally distributed trade flows. The findings hint towards non-negligible leakage effects from unilateral environmental policy actions and underline the need to improve international coordination of climate change policies. In the absence of such coordination, border carbon adjustment might be a possible remedy to overcome the disincentive for unilateral climate policies resulting from evidence for carbon leakage.Footnote 23

An important shortcoming of the gravity methodology used in this and other studies is its partial equilibrium nature. Third-country effects are by construction ruled out. In general equilibrium, the change in the carbon trade balance between Germany and China after Germany ratified and China did not ratify would be expected to also affect their carbon trade balance with the United States. Future research will be directed towards incorporating these effects in the estimation of carbon leakage.