3.1. Framework to Assess the Impact of Mitigation Measures to Reduce Global Shipping Emissions
We estimate the future CO
2 emission from international shipping using the International Transport Forum’s International freight model (IFM) and the “ASIF” (Activity, Structure, Intensity, Emission Factor) method [
42].
The ITF’s International freight model is designed to project international freight transport activities (in tonne kilometers) for 19 commodities for all major transport modes and routes while taking into account different transport and economic policy measures (e.g., the development of new infrastructure networks, or the alleviation of trade barriers). The model is built on the four-steps freight transportation modelling approach and it takes the OECD trade projection as an input. The IFM is designed to be able to estimate the weight of commodities traded between countries, the choice between modes and transport routes used to transport these commodities based on transport networks characteristics, and relevant socio-economic variables such as transport costs and time. The model consists of the following components:
3.1.1. OECD International Trade Model
The OECD’s trade projection is produced using a Computable General Equilibrium (CGE) model called the ENV-Linkages model [
43]. The model is designed to estimate the dynamic evolution of international trade, in terms of both spatial patterns and commodity composition due to the changes in the global production and consumption of commodities. It is calibrated based on the macroeconomic trends of the OECD@100’s baseline scenario for the period 2013–2060 at sectorial and regional levels. As such, it projects international trade flows in values (US
$) for 26 regions and 25 commodities until 2060.
3.1.2. Trade Disaggregation Model
The underlying trade projections are disaggregated into 26 world regions. This level of resolution does not allow estimating transport flows with precision as it does not allow a proper discretization of the travel path used for different types of products. Therefore, we disaggregate the regional origin–destination (OD) trade flows into a larger number of production/consumption centroids. These centroids were calculated using an adapted
p-median procedure for all the cities around the world classified by United Nations in 2010 relative to their population (2539 cities). The objective function for this aggregation is based on the minimization of a distance function which includes two components: GDP density and geographical distance. The selection was also constrained by allowing one centroid within a 500 km radius in a country. This resulted in 333 centroids globally, with spatially balanced results also for all continents.
In Equation (1),
= trade values from centroid o to centroid d in year y for commodity k,
= trade values from origin region V to destination region L,
o, d = origin and destination centroids,
k = commodity k,
y = year of analysis,
k = centroid that belongs to the origin region V,
l = centroid that belongs to the destination region L.
3.1.3. Value-to-Weight Model
We used a Poisson regression model to estimate the rate of conversion of value units (dollars) into weight units of cargo (tonnes) by mode, calibrated using datasets from Eurostat and Economic Commission for Latin America and the Caribbean (ECLAC) data on value/weight ratios for different commodities.
We use the natural logarithm of the trade value in millions of dollars as the offset variable, with panel terms by commodity, a transport cost proxy variable (logsum calculation for maritime, road, rail, and air transport costs per ton between each pair of centroids), and geographical and cultural variables: binary variables for trade agreements and land borders used above and a binary variable identifying if two countries have the same official language. Moreover, economic profile variables were included to describe the trade relation between countries with different types of production sophistication and scale of trade intensity. We validate the output of the value-to-weight model using the UN Comtrade database that provides values and weights of all commodities traded between any countries worldwide.
Table A1 provides validation of the total values and weights of global trade produced by the model.
In Equations (2) and (3),
= weight of commodity k that is traded between origin o and destination d for year y (in tonnes),
= value of trade for commodity k between origin o and destination d for year y (in US$),
= value-to-weight conversion factor for commodity k, between origin o and destination d for year y (in tonnes/US$),
= GDP percentile of origin in year y,
= GDP percentile of destination in year y,
= GDP per capita percentile of origin in year y,
= GDP per capita percentile of destination in year y,
= natural logarithm of the ratio between GDP per capita of origin and GDP per capita of destination in year y,
= land contiguity between origin o and destination d, contig = (0, 1),
langod = shared language between origin o and destination d, lang = (0, 1),
rtaod = trade agreement between origin o and destination d, rta = (0,1),
logsum(costod) = logsum variable of transport costs using different modes between origin o and destination d, d
lgsk = logsum coefficient/panel term for commodity k.
3.1.4. Mode Choice Model
The mode share model (in weight) for international freight flows assigns the transport mode used for trade between any origin–destination pair of centroids. The mode attributed to each trade connection represents the longest transport section. All freight will require intermodal transport both at the origin and destination. This domestic component of international freight is usually not accounted for in the literature, but is included in our model. The model is estimated using a standard multinomial logit estimator including commodity type panel terms on travel times and cost. Both Eurostat and ECLAC datasets are used as sources of observation data for the volume of commodities and its mode of transport. Transport costs and travel times are estimated using the network model and observed data whenever available. Two geographical and economic context binary variables are added, one describing if the OD pair has a trade agreement and the other for the existence of a land border between trading partners. The mode choice model is validated by ensuring the mode share of the volume of goods transported is similar to the observed mode share for international transport in 2011 by weight. Additionally, the total tonne kilometers for all 4 major modes of transport (air, road, rail, sea) are also validated against the observed data. These observed data are obtained from reports of various organizations such as the International Maritime Organization (IMO), the International Civil Aviation Organization (ICAO), and the World Bank.
Table A1 and
Table A2 provide detailed descriptions of the validation result.
In Equations (4) and (5),
Pm = the choice probability of mode m,
= the choice utility of mode m for commodity k between origin o and destination d,
ascm = alternative specific constant for mode m,
CFk = transport cost coefficient for commodity k,
= transport cost for mode m between origin o and destination d,
TFk = travel time coefficient for commodity k,
= contiguity coefficient for mode m,
contigod = contiguity variable between origin o and destination d, contig = (0, 1),
Rt = trade agreement coefficient,
rtaod = trade agreement variable between origin o and destination d, rta = (0, 1).
3.1.5. Route Choice Model
We used a path size logit model in combination with a path generation method to assign the volume of freight transport across all possible international shipping routes between all origins and destinations. The model does this using a shortest path algorithm and choice set creation algorithm to identify the subsegments of the complete shortest route for each port-to-port segment of a shipping line. The model accounts both for maritime connections between two countries and for overland connections between the centroids. The route and port choice algorithms use a path size logit model which takes overlaps between the alternative routes into account and distinguishes the transport costs associated with these alternatives properly. The basis of this model can be found in [
44]. The model is calibrated by minimizing the difference between observed and modelled port throughputs for more than 400 major ports in the world. A detailed description on the model can be found in [
45] or in [
46]. The formal definition of the cost model is delineated below:
In Equation (6):
Cr = unit cost of route r from origin centroid to destination centroid (US$/Twenty-equivalent unit, TEU),
p = ports used by the route,
l = links used by the route,
Ap = unit cost of transhipment at port p (US$/TEU),
cl = unit cost of transportation over link l (US$/TEU),
Tp = time spent during transhipment at port p (days/TEU),
tl = time spent during transportation over link l (days/TEU),
α = value of transport time (US$/day).
The model accounts both for maritime connections between two countries and for overland connections between these countries. The route and port choice algorithms use a path size logit model which takes overlaps between the alternative routes into account and distinguishes the transport costs associated with these alternatives properly. The basis of this model can be found in [
44]. The following is the formal definition of the route choice model. The route probabilities are given by
while the path size overlap variable
S is defined as
In Equations (7) and (8):
Pr = the choice probability of route r,
Cr = generalized costs of route r,
Ch = generalized costs of route h within the choice set,
CS = the choice set with multiple routes,
h = path indicator/index, ,
μ = logit scale parameter,
a = link in route r,
Sr = degree of path overlap,
Lkr = set of links in route r,
Za = length of link a,
Zr = length of route r,
Nah = number of times link a is found in alternative routes.
3.1.6. CO2 Mitigation Impact Assessment Model
The ASIF framework is used to assess the impact of the maximum possible technical, operational, and alternative fuels measures on the total CO
2 emissions of international shipping (
Figure 1). The output of IFM provides the tonne kilometers for different commodities (“Activity”) and the projections for future transport demand scenarios which we assign to possible ship types to estimate the activities for each ship type. We consider 6 ship types in our model: dry bulk, container, oil tanker, gas, wet product and chemical, and general cargo. We estimate the vehicle kilometers for each ship type using the ship’s load factor data and the projection for the changes in ship size until 2035 (“Structure”). Furthermore, we compute the fuel consumptions of all ship types using engine efficiency improvement pathways together with the distribution of fuel types across different ship types (“Intensity”). The resulting fuel consumption for each ship type is then used to estimate the total CO
2 emissions using carbon factor and energy content data for different fuel types (“Emission factor”).
In Equations (9)–(11),
TCO2 = total CO2 emissions from international shipping,
Eqf = total CO2 emissions from ship type q using fuel type f (in tonnes CO2),
Aq = total annual activity for ship q (in tonne kilometers),
Sqf = total vehicle kilometers for ship type q which uses fuel type f (in vkm),
SHqf = share of ship type q which uses fuel type f (in %),
LFq = load factor of ship type q (in tonnes/vehicle),
Iq = engine intensity of ship type q (in MJ/vkm),
Ff = emission factor of fuel type f (in tonnes CO2/MJ).
3.2. Data
We use data from various datasets to estimate each component of the model and produce the baseline CO2 emissions for each ship type until 2035 using the ASIF framework.
The weight-to-value model and the mode choice model are estimated using observed international trade flows from Eurostat and ECLAC datasets. The Eurostat dataset registers trade flows between Europe and the rest of the world that are obtained from the customs of each EU country, and the ECLAC dataset records trade data between Latin American countries and countries worldwide. These datasets combined provide more than 17,427 observed trade flows in values (US$) and weights (tonnes) and their modes of transport. We use the UN Comtrade data that record trade flow data for all commodities, grouped in more than 97 product types (chapters) both in value and weight, to validate the output of the model. It is necessary to aggregate the trade flows both in values and weights at the country-to-country level for all commodities in the Harmonized Commodity Description and Coding Systems (HS) to enable a comparison with the output of our model. The GDP and population data, including their projection until 2050, were obtained from the environment directorate of the OECD.
Given the baseline transport demand projection for 6 ship types, the ships’ vehicle kilometers are estimated using load factor data for each ship type. The load factor data for each ship type are estimated by multiplying the average freight capacity and average utilization rate of each ship type. We obtained average freight capacity data from the UNCTAD Review of Maritime Transport 2013 [
47] and ships’ utilization rate data from the 3rd IMO GHG study [
1]. We estimated the evolution in ships’ load factor by taking into account the future evolution in ship size (
Figure A1). This projection of ship size is based on the observed historical pattern of ship sizes from 1996 to 2015. The pathways for ships’ engine efficiency for each ship type for the baseline scenario were obtained from a UMAS study [
7,
48]. Furthermore, the emission factor data for different fuel types were obtained from the International Energy Agency’s Mobility Model [
49] (
Table A3). By multiplying the engine efficiency (in MJ/vkm) with the emission factor for different fuel types (in CO
2/MJ), we obtain the carbon intensity of each ship type (
Figure 2).